Epigenetic memory of coronavirus infection in innate immune cells and their progenitors.

Inflammation can trigger lasting phenotypes in immune and non-immune cells. Whether and how human infections and associated inflammation can form innate immune memory in hematopoietic stem and progenitor cells (HSPC) has remained unclear. We found that circulating HSPC, enriched from peripheral blood, captured the diversity of bone marrow HSPC, enabling investigation of their epigenomic reprogramming following coronavirus disease 2019 (COVID-19). Alterations in innate immune phenotypes and epigenetic programs of HSPC persisted for months to 1 year following severe COVID-19 and were associated with distinct transcription factor (TF) activities, altered regulation of inflammatory programs, and durable increases in myelopoiesis. HSPC epigenomic alterations were conveyed, through differentiation, to progeny innate immune cells. Early activity of IL-6 contributed to these persistent phenotypes in human COVID-19 and a mouse coronavirus infection model. Epigenetic reprogramming of HSPC may underlie altered immune function following infection and be broadly relevant, especially for millions of COVID-19 survivors.

Deficiency of metabolic regulator PKM2 activates the pentose phosphate pathway and generates TCF1+ progenitor CD8+ T cells to improve immunotherapy.

TCF1high progenitor CD8+ T cells mediate the efficacy of immunotherapy; however, the mechanisms that govern their generation and maintenance are poorly understood. Here, we show that targeting glycolysis through deletion of pyruvate kinase muscle 2 (PKM2) results in elevated pentose phosphate pathway (PPP) activity, leading to enrichment of a TCF1high progenitor-exhausted-like phenotype and increased responsiveness to PD-1 blockade in vivo. PKM2KO CD8+ T cells showed reduced glycolytic flux, accumulation of glycolytic intermediates and PPP metabolites and increased PPP cycling as determined by 1,2-13C glucose carbon tracing. Small molecule agonism of the PPP without acute glycolytic impairment skewed CD8+ T cells toward a TCF1high population, generated a unique transcriptional landscape and adoptive transfer of agonist-treated CD8+ T cells enhanced tumor control in mice in combination with PD-1 blockade and promoted tumor killing in patient-derived tumor organoids. Our study demonstrates a new metabolic reprogramming that contributes to a progenitor-like T cell state promoting immunotherapy efficacy.

TBL1XR1 Mutations Drive Extranodal Lymphoma by Inducing a Pro-tumorigenic Memory Fate.

The most aggressive B cell lymphomas frequently manifest extranodal distribution and carry somatic mutations in the poorly characterized gene TBL1XR1. Here, we show that TBL1XR1 mutations skew the humoral immune response toward generating abnormal immature memory B cells (MB), while impairing plasma cell differentiation. At the molecular level, TBL1XR1 mutants co-opt SMRT/HDAC3 repressor complexes toward binding the MB cell transcription factor (TF) BACH2 at the expense of the germinal center (GC) TF BCL6, leading to pre-memory transcriptional reprogramming and cell-fate bias. Upon antigen recall, TBL1XR1 mutant MB cells fail to differentiate into plasma cells and instead preferentially reenter new GC reactions, providing evidence for a cyclic reentry lymphomagenesis mechanism. Ultimately, TBL1XR1 alterations lead to a striking extranodal immunoblastic lymphoma phenotype that mimics the human disease. Both human and murine lymphomas feature expanded MB-like cell populations, consistent with a MB-cell origin and delineating an unforeseen pathway for malignant transformation of the immune system.

OCT2 pre-positioning facilitates cell fate transition and chromatin architecture changes in humoral immunity.

During the germinal center (GC) reaction, B cells undergo profound transcriptional, epigenetic and genomic architectural changes. How such changes are established remains unknown. Mapping chromatin accessibility during the humoral immune response, we show that OCT2 was the dominant transcription factor linked to differential accessibility of GC regulatory elements. Silent chromatin regions destined to become GC-specific super-enhancers (SEs) contained pre-positioned OCT2-binding sites in naive B cells (NBs). These preloaded SE 'seeds' featured spatial clustering of regulatory elements enriched in OCT2 DNA-binding motifs that became heavily loaded with OCT2 and its GC-specific coactivator OCAB in GC B cells (GCBs). SEs with high abundance of pre-positioned OCT2 binding preferentially formed long-range chromatin contacts in GCs, to support expression of GC-specifying factors. Gain in accessibility and architectural interactivity of these regions were dependent on recruitment of OCAB. Pre-positioning key regulators at SEs may represent a broadly used strategy for facilitating rapid cell fate transitions.

Smc3 dosage regulates B cell transit through germinal centers and restricts their malignant transformation.

During the germinal center (GC) reaction, B cells undergo extensive redistribution of cohesin complex and three-dimensional reorganization of their genomes. Yet, the significance of cohesin and architectural programming in the humoral immune response is unknown. Herein we report that homozygous deletion of Smc3, encoding the cohesin ATPase subunit, abrogated GC formation, while, in marked contrast, Smc3 haploinsufficiency resulted in GC hyperplasia, skewing of GC polarity and impaired plasma cell (PC) differentiation. Genome-wide chromosomal conformation and transcriptional profiling revealed defects in GC B cell terminal differentiation programs controlled by the lymphoma epigenetic tumor suppressors Tet2 and Kmt2d and failure of Smc3-haploinsufficient GC B cells to switch from B cell- to PC-defining transcription factors. Smc3 haploinsufficiency preferentially impaired the connectivity of enhancer elements controlling various lymphoma tumor suppressor genes, and, accordingly, Smc3 haploinsufficiency accelerated lymphomagenesis in mice with constitutive Bcl6 expression. Collectively, our data indicate a dose-dependent function for cohesin in humoral immunity to facilitate the B cell to PC phenotypic switch while restricting malignant transformation.

Stem Cell Lineage Infidelity Drives Wound Repair and Cancer.

Tissue stem cells contribute to tissue regeneration and wound repair through cellular programs that can be hijacked by cancer cells. Here, we investigate such a phenomenon in skin, where during homeostasis, stem cells of the epidermis and hair follicle fuel their respective tissues. We find that breakdown of stem cell lineage confinement-granting privileges associated with both fates-is not only hallmark but also functional in cancer development. We show that lineage plasticity is critical in wound repair, where it operates transiently to redirect fates. Investigating mechanism, we discover that irrespective of cellular origin, lineage infidelity occurs in wounding when stress-responsive enhancers become activated and override homeostatic enhancers that govern lineage specificity. In cancer, stress-responsive transcription factor levels rise, causing lineage commanders to reach excess. When lineage and stress factors collaborate, they activate oncogenic enhancers that distinguish cancers from wounds.

5' UTR m(6)A Promotes Cap-Independent Translation.

Protein translation typically begins with the recruitment of the 43S ribosomal complex to the 5' cap of mRNAs by a cap-binding complex. However, some transcripts are translated in a cap-independent manner through poorly understood mechanisms. Here, we show that mRNAs containing N(6)-methyladenosine (m(6)A) in their 5' UTR can be translated in a cap-independent manner. A single 5' UTR m(6)A directly binds eukaryotic initiation factor 3 (eIF3), which is sufficient to recruit the 43S complex to initiate translation in the absence of the cap-binding factor eIF4E. Inhibition of adenosine methylation selectively reduces translation of mRNAs containing 5'UTR m(6)A. Additionally, increased m(6)A levels in the Hsp70 mRNA regulate its cap-independent translation following heat shock. Notably, we find that diverse cellular stresses induce a transcriptome-wide redistribution of m(6)A, resulting in increased numbers of mRNAs with 5' UTR m(6)A. These data show that 5' UTR m(6)A bypasses 5' cap-binding proteins to promote translation under stresses.

The BRAF pseudogene functions as a competitive endogenous RNA and induces lymphoma in vivo.

Research over the past decade has suggested important roles for pseudogenes in physiology and disease. In vitro experiments demonstrated that pseudogenes contribute to cell transformation through several mechanisms. However, in vivo evidence for a causal role of pseudogenes in cancer development is lacking. Here, we report that mice engineered to overexpress either the full-length murine B-Raf pseudogene Braf-rs1 or its pseudo "CDS" or "3' UTR" develop an aggressive malignancy resembling human diffuse large B cell lymphoma. We show that Braf-rs1 and its human ortholog, BRAFP1, elicit their oncogenic activity, at least in part, as competitive endogenous RNAs (ceRNAs) that elevate BRAF expression and MAPK activation in vitro and in vivo. Notably, we find that transcriptional or genomic aberrations of BRAFP1 occur frequently in multiple human cancers, including B cell lymphomas. Our engineered mouse models demonstrate the oncogenic potential of pseudogenes and indicate that ceRNA-mediated microRNA sequestration may contribute to the development of cancer.

The nuclear factor ID3 endows macrophages with a potent anti-tumour activity.

Macrophage activation is controlled by a balance between activating and inhibitory receptors1-7, which protect normal tissues from excessive damage during infection8,9 but promote tumour growth and metastasis in cancer7,10. Here we report that the Kupffer cell lineage-determining factor ID3 controls this balance and selectively endows Kupffer cells with the ability to phagocytose live tumour cells and orchestrate the recruitment, proliferation and activation of natural killer and CD8 T lymphoid effector cells in the liver to restrict the growth of a variety of tumours. ID3 shifts the macrophage inhibitory/activating receptor balance to promote the phagocytic and lymphoid response, at least in part by buffering the binding of the transcription factors ELK1 and E2A at the SIRPA locus. Furthermore, loss- and gain-of-function experiments demonstrate that ID3 is sufficient to confer this potent anti-tumour activity to mouse bone-marrow-derived macrophages and human induced pluripotent stem-cell-derived macrophages. Expression of ID3 is therefore necessary and sufficient to endow macrophages with the ability to form an efficient anti-tumour niche, which could be harnessed for cell therapy in cancer.

The interplay of mutagenesis and ecDNA shapes urothelial cancer evolution.

Advanced urothelial cancer is a frequently lethal disease characterized by marked genetic heterogeneity1. In this study, we investigated the evolution of genomic signatures caused by endogenous and external mutagenic processes and their interplay with complex structural variants (SVs). We superimposed mutational signatures and phylogenetic analyses of matched serial tumours from patients with urothelial cancer to define the evolutionary dynamics of these processes. We show that APOBEC3-induced mutations are clonal and early, whereas chemotherapy induces mutational bursts of hundreds of late subclonal mutations. Using a genome graph computational tool2, we observed frequent high copy-number circular amplicons characteristic of extrachromosomal DNA (ecDNA)-forming SVs. We characterized the distinct temporal patterns of APOBEC3-induced and chemotherapy-induced mutations within ecDNA-forming SVs, gaining new insights into the timing of these mutagenic processes relative to ecDNA biogenesis. We discovered that most CCND1 amplifications in urothelial cancer arise within circular ecDNA-forming SVs. ecDNA-forming SVs persisted and increased in complexity, incorporating additional DNA segments and contributing to the evolution of treatment resistance. Oxford Nanopore Technologies long-read whole-genome sequencing followed by de novo assembly mapped out CCND1 ecDNA structure. Experimental modelling of CCND1 ecDNA confirmed its role as a driver of treatment resistance. Our findings define fundamental mechanisms that drive urothelial cancer evolution and have important therapeutic implications.

The intrinsic substrate specificity of the human tyrosine kinome.

Phosphorylation of proteins on tyrosine (Tyr) residues evolved in metazoan organisms as a mechanism of coordinating tissue growth1. Multicellular eukaryotes typically have more than 50 distinct protein Tyr kinases that catalyse the phosphorylation of thousands of Tyr residues throughout the proteome1-3. How a given Tyr kinase can phosphorylate a specific subset of proteins at unique Tyr sites is only partially understood4-7. Here we used combinatorial peptide arrays to profile the substrate sequence specificity of all human Tyr kinases. Globally, the Tyr kinases demonstrate considerable diversity in optimal patterns of residues surrounding the site of phosphorylation, revealing the functional organization of the human Tyr kinome by substrate motif preference. Using this information, Tyr kinases that are most compatible with phosphorylating any Tyr site can be identified. Analysis of mass spectrometry phosphoproteomic datasets using this compendium of kinase specificities accurately identifies specific Tyr kinases that are dysregulated in cells after stimulation with growth factors, treatment with anti-cancer drugs or expression of oncogenic variants. Furthermore, the topology of known Tyr signalling networks naturally emerged from a comparison of the sequence specificities of the Tyr kinases and the SH2 phosphotyrosine (pTyr)-binding domains. Finally we show that the intrinsic substrate specificity of Tyr kinases has remained fundamentally unchanged from worms to humans, suggesting that the fidelity between Tyr kinases and their protein substrate sequences has been maintained across hundreds of millions of years of evolution.

Long-molecule scars of backup DNA repair in BRCA1- and BRCA2-deficient cancers.

Homologous recombination (HR) deficiency is associated with DNA rearrangements and cytogenetic aberrations1. Paradoxically, the types of DNA rearrangements that are specifically associated with HR-deficient cancers only minimally affect chromosomal structure2. Here, to address this apparent contradiction, we combined genome-graph analysis of short-read whole-genome sequencing (WGS) profiles across thousands of tumours with deep linked-read WGS of 46 BRCA1- or BRCA2-mutant breast cancers. These data revealed a distinct class of HR-deficiency-enriched rearrangements called reciprocal pairs. Linked-read WGS showed that reciprocal pairs with identical rearrangement orientations gave rise to one of two distinct chromosomal outcomes, distinguishable only with long-molecule data. Whereas one (cis) outcome corresponded to the copying and pasting of a small segment to a distant site, a second (trans) outcome was a quasi-balanced translocation or multi-megabase inversion with substantial (10 kb) duplications at each junction. We propose an HR-independent replication-restart repair mechanism to explain the full spectrum of reciprocal pair outcomes. Linked-read WGS also identified single-strand annealing as a repair pathway that is specific to BRCA2 deficiency in human cancers. Integrating these features in a classifier improved discrimination between BRCA1- and BRCA2-deficient genomes. In conclusion, our data reveal classes of rearrangements that are specific to BRCA1 or BRCA2 deficiency as a source of cytogenetic aberrations in HR-deficient cells.

An atlas of substrate specificities for the human serine/threonine kinome.

Protein phosphorylation is one of the most widespread post-translational modifications in biology1,2. With advances in mass-spectrometry-based phosphoproteomics, 90,000 sites of serine and threonine phosphorylation have so far been identified, and several thousand have been associated with human diseases and biological processes3,4. For the vast majority of phosphorylation events, it is not yet known which of the more than 300 protein serine/threonine (Ser/Thr) kinases encoded in the human genome are responsible3. Here we used synthetic peptide libraries to profile the substrate sequence specificity of 303 Ser/Thr kinases, comprising more than 84% of those predicted to be active in humans. Viewed in its entirety, the substrate specificity of the kinome was substantially more diverse than expected and was driven extensively by negative selectivity. We used our kinome-wide dataset to computationally annotate and identify the kinases capable of phosphorylating every reported phosphorylation site in the human Ser/Thr phosphoproteome. For the small minority of phosphosites for which the putative protein kinases involved have been previously reported, our predictions were in excellent agreement. When this approach was applied to examine the signalling response of tissues and cell lines to hormones, growth factors, targeted inhibitors and environmental or genetic perturbations, it revealed unexpected insights into pathway complexity and compensation. Overall, these studies reveal the intrinsic substrate specificity of the human Ser/Thr kinome, illuminate cellular signalling responses and provide a resource to link phosphorylation events to biological pathways.

Punctuated evolution of prostate cancer genomes.

The analysis of exonic DNA from prostate cancers has identified recurrently mutated genes, but the spectrum of genome-wide alterations has not been profiled extensively in this disease. We sequenced the genomes of 57 prostate tumors and matched normal tissues to characterize somatic alterations and to study how they accumulate during oncogenesis and progression. By modeling the genesis of genomic rearrangements, we identified abundant DNA translocations and deletions that arise in a highly interdependent manner. This phenomenon, which we term "chromoplexy," frequently accounts for the dysregulation of prostate cancer genes and appears to disrupt multiple cancer genes coordinately. Our modeling suggests that chromoplexy may induce considerable genomic derangement over relatively few events in prostate cancer and other neoplasms, supporting a model of punctuated cancer evolution. By characterizing the clonal hierarchy of genomic lesions in prostate tumors, we charted a path of oncogenic events along which chromoplexy may drive prostate carcinogenesis.

Unique Immune Cell Coactivators Specify Locus Control Region Function and Cell Stage.

Locus control region (LCR) functions define cellular identity and have critical roles in diseases such as cancer, although the hierarchy of structural components and associated factors that drive functionality are incompletely understood. Here we show that OCA-B, a B cell-specific coactivator essential for germinal center (GC) formation, forms a ternary complex with the lymphoid-enriched OCT2 and GC-specific MEF2B transcription factors and that this complex occupies and activates an LCR that regulates the BCL6 proto-oncogene and is uniquely required by normal and malignant GC B cells. Mechanistically, through OCA-B-MED1 interactions, this complex is required for Mediator association with the BCL6 promoter. Densely tiled CRISPRi screening indicates that only LCR segments heavily bound by this ternary complex are essential for its function. Our results demonstrate how an intimately linked complex of lineage- and stage-specific factors converges on specific and highly essential enhancer elements to drive the function of a cell-type-defining LCR.

Comprehensive analysis of mRNA methylation reveals enrichment in 3' UTRs and near stop codons.

Methylation of the N(6) position of adenosine (m(6)A) is a posttranscriptional modification of RNA with poorly understood prevalence and physiological relevance. The recent discovery that FTO, an obesity risk gene, encodes an m(6)A demethylase implicates m(6)A as an important regulator of physiological processes. Here, we present a method for transcriptome-wide m(6)A localization, which combines m(6)A-specific methylated RNA immunoprecipitation with next-generation sequencing (MeRIP-Seq). We use this method to identify mRNAs of 7,676 mammalian genes that contain m(6)A, indicating that m(6)A is a common base modification of mRNA. The m(6)A modification exhibits tissue-specific regulation and is markedly increased throughout brain development. We find that m(6)A sites are enriched near stop codons and in 3' UTRs, and we uncover an association between m(6)A residues and microRNA-binding sites within 3' UTRs. These findings provide a resource for identifying transcripts that are substrates for adenosine methylation and reveal insights into the epigenetic regulation of the mammalian transcriptome.

Efficient direct reprogramming of mature amniotic cells into endothelial cells by ETS factors and TGFß suppression.

ETS transcription factors ETV2, FLI1, and ERG1 specify pluripotent stem cells into induced vascular endothelial cells (iVECs). However, iVECs are unstable and drift toward nonvascular cells. We show that human midgestation c-Kit(-) lineage-committed amniotic cells (ACs) can be reprogrammed into vascular endothelial cells (rAC-VECs) without transitioning through a pluripotent state. Transient ETV2 expression in ACs generates immature rAC-VECs, whereas coexpression with FLI1/ERG1 endows rAC-VECs with a vascular repertoire and morphology matching mature endothelial cells (ECs). Brief TGFß-inhibition functionalizes VEGFR2 signaling, augmenting specification of ACs into rAC-VECs. Genome-wide transcriptional analyses showed that rAC-VECs are similar to adult ECs in which vascular-specific genes are expressed and nonvascular genes are silenced. Functionally, rAC-VECs form stable vasculature in Matrigel plugs and regenerating livers. Therefore, short-term ETV2 expression and TGFß inhibition with constitutive ERG1/FLI1 coexpression reprogram mature ACs into durable rAC-VECs with clinical-scale expansion potential. Banking of HLA-typed rAC-VECs establishes a vascular inventory for treatment of diverse disorders.

The spatial landscape of lung pathology during COVID-19 progression.

Recent studies have provided insights into the pathology of and immune response to COVID-191-8. However, a thorough investigation of the interplay between infected cells and the immune system at sites of infection has been lacking. Here we use high-parameter imaging mass cytometry9 that targets the expression of 36 proteins to investigate the cellular composition and spatial architecture of acute lung injury in humans (including injuries derived from SARS-CoV-2 infection) at single-cell resolution. These spatially resolved single-cell data unravel the disordered structure of the infected and injured lung, alongside the distribution of extensive immune infiltration. Neutrophil and macrophage infiltration are hallmarks of bacterial pneumonia and COVID-19, respectively. We provide evidence that SARS-CoV-2 infects predominantly alveolar epithelial cells and induces a localized hyperinflammatory cell state that is associated with lung damage. We leverage the temporal range of fatal outcomes of COVID-19 in relation to the onset of symptoms, which reveals increased macrophage extravasation and increased numbers of mesenchymal cells and fibroblasts concomitant with increased proximity between these cell types as the disease progresses-possibly as a result of attempts to repair the damaged lung tissue. Our data enable us to develop a biologically interpretable landscape of lung pathology from a structural, immunological and clinical standpoint. We use this landscape to characterize the pathophysiology of the human lung from its macroscopic presentation to the single-cell level, which provides an important basis for understanding COVID-19 and lung pathology in general.

Histone H1 loss drives lymphoma by disrupting 3D chromatin architecture.

Linker histone H1 proteins bind to nucleosomes and facilitate chromatin compaction1, although their biological functions are poorly understood. Mutations in the genes that encode H1 isoforms B-E (H1B, H1C, H1D and H1E; also known as H1-5, H1-2, H1-3 and H1-4, respectively) are highly recurrent in B cell lymphomas, but the pathogenic relevance of these mutations to cancer and the mechanisms that are involved are unknown. Here we show that lymphoma-associated H1 alleles are genetic driver mutations in lymphomas. Disruption of H1 function results in a profound architectural remodelling of the genome, which is characterized by large-scale yet focal shifts of chromatin from a compacted to a relaxed state. This decompaction drives distinct changes in epigenetic states, primarily owing to a gain of histone H3 dimethylation at lysine 36 (H3K36me2) and/or loss of repressive H3 trimethylation at lysine 27 (H3K27me3). These changes unlock the expression of stem cell genes that are normally silenced during early development. In mice, loss of H1c and H1e (also known as H1f2 and H1f4, respectively) conferred germinal centre B cells with enhanced fitness and self-renewal properties, ultimately leading to aggressive lymphomas with an increased repopulating potential. Collectively, our data indicate that H1 proteins are normally required to sequester early developmental genes into architecturally inaccessible genomic compartments. We also establish H1 as a bona fide tumour suppressor and show that mutations in H1 drive malignant transformation primarily through three-dimensional genome reorganization, which leads to epigenetic reprogramming and derepression of developmentally silenced genes.

A molecular single-cell lung atlas of lethal COVID-19.

Respiratory failure is the leading cause of death in patients with severe SARS-CoV-2 infection1,2, but the host response at the lung tissue level is poorly understood. Here we performed single-nucleus RNA sequencing of about 116,000 nuclei from the lungs of nineteen individuals who died of COVID-19 and underwent rapid autopsy and seven control individuals. Integrated analyses identified substantial alterations in cellular composition, transcriptional cell states, and cell-to-cell interactions, thereby providing insight into the biology of lethal COVID-19. The lungs from individuals with COVID-19 were highly inflamed, with dense infiltration of aberrantly activated monocyte-derived macrophages and alveolar macrophages, but had impaired T cell responses. Monocyte/macrophage-derived interleukin-1ß and epithelial cell-derived interleukin-6 were unique features of SARS-CoV-2 infection compared to other viral and bacterial causes of pneumonia. Alveolar type 2 cells adopted an inflammation-associated transient progenitor cell state and failed to undergo full transition into alveolar type 1 cells, resulting in impaired lung regeneration. Furthermore, we identified expansion of recently described CTHRC1+ pathological fibroblasts3 contributing to rapidly ensuing pulmonary fibrosis in COVID-19. Inference of protein activity and ligand-receptor interactions identified putative drug targets to disrupt deleterious circuits. This atlas enables the dissection of lethal COVID-19, may inform our understanding of long-term complications of COVID-19 survivors, and provides an important resource for therapeutic development.