Pan-cancer analysis of post-translational modifications reveals shared patterns of protein regulation.

Post-translational modifications (PTMs) play key roles in regulating cell signaling and physiology in both normal and cancer cells. Advances in mass spectrometry enable high-throughput, accurate, and sensitive measurement of PTM levels to better understand their role, prevalence, and crosstalk. Here, we analyze the largest collection of proteogenomics data from 1,110 patients with PTM profiles across 11 cancer types (10 from the National Cancer Institute's Clinical Proteomic Tumor Analysis Consortium [CPTAC]). Our study reveals pan-cancer patterns of changes in protein acetylation and phosphorylation involved in hallmark cancer processes. These patterns revealed subsets of tumors, from different cancer types, including those with dysregulated DNA repair driven by phosphorylation, altered metabolic regulation associated with immune response driven by acetylation, affected kinase specificity by crosstalk between acetylation and phosphorylation, and modified histone regulation. Overall, this resource highlights the rich biology governed by PTMs and exposes potential new therapeutic avenues.

Pan-cancer proteogenomics connects oncogenic drivers to functional states.

Cancer driver events refer to key genetic aberrations that drive oncogenesis; however, their exact molecular mechanisms remain insufficiently understood. Here, our multi-omics pan-cancer analysis uncovers insights into the impacts of cancer drivers by identifying their significant cis-effects and distal trans-effects quantified at the RNA, protein, and phosphoprotein levels. Salient observations include the association of point mutations and copy-number alterations with the rewiring of protein interaction networks, and notably, most cancer genes converge toward similar molecular states denoted by sequence-based kinase activity profiles. A correlation between predicted neoantigen burden and measured T cell infiltration suggests potential vulnerabilities for immunotherapies. Patterns of cancer hallmarks vary by polygenic protein abundance ranging from uniform to heterogeneous. Overall, our work demonstrates the value of comprehensive proteogenomics in understanding the functional states of oncogenic drivers and their links to cancer development, surpassing the limitations of studying individual cancer types.

Potent cross-reactive antibodies following Omicron breakthrough in vaccinees.

Highly transmissible Omicron variants of SARS-CoV-2 currently dominate globally. Here, we compare neutralization of Omicron BA.1, BA.1.1, and BA.2. BA.2 RBD has slightly higher ACE2 affinity than BA.1 and slightly reduced neutralization by vaccine serum, possibly associated with its increased transmissibility. Neutralization differences between sub-lineages for mAbs (including therapeutics) mostly arise from variation in residues bordering the ACE2 binding site; however, more distant mutations S371F (BA.2) and R346K (BA.1.1) markedly reduce neutralization by therapeutic antibody Vir-S309. In-depth structure-and-function analyses of 27 potent RBD-binding mAbs isolated from vaccinated volunteers following breakthrough Omicron-BA.1 infection reveals that they are focused in two main clusters within the RBD, with potent right-shoulder antibodies showing increased prevalence. Selection and somatic maturation have optimized antibody potency in less-mutated epitopes and recovered potency in highly mutated epitopes. All 27 mAbs potently neutralize early pandemic strains, and many show broad reactivity with variants of concern.

Antibody escape of SARS-CoV-2 Omicron BA.4 and BA.5 from vaccine and BA.1 serum.

The Omicron lineage of SARS-CoV-2, which was first described in November 2021, spread rapidly to become globally dominant and has split into a number of sublineages. BA.1 dominated the initial wave but has been replaced by BA.2 in many countries. Recent sequencing from South Africa's Gauteng region uncovered two new sublineages, BA.4 and BA.5, which are taking over locally, driving a new wave. BA.4 and BA.5 contain identical spike sequences, and although closely related to BA.2, they contain further mutations in the receptor-binding domain of their spikes. Here, we study the neutralization of BA.4/5 using a range of vaccine and naturally immune serum and panels of monoclonal antibodies. BA.4/5 shows reduced neutralization by the serum from individuals vaccinated with triple doses of AstraZeneca or Pfizer vaccine compared with BA.1 and BA.2. Furthermore, using the serum from BA.1 vaccine breakthrough infections, there are, likewise, significant reductions in the neutralization of BA.4/5, raising the possibility of repeat Omicron infections.

SARS-CoV-2 Omicron-B.1.1.529 leads to widespread escape from neutralizing antibody responses.

On 24th November 2021, the sequence of a new SARS-CoV-2 viral isolate Omicron-B.1.1.529 was announced, containing far more mutations in Spike (S) than previously reported variants. Neutralization titers of Omicron by sera from vaccinees and convalescent subjects infected with early pandemic Alpha, Beta, Gamma, or Delta are substantially reduced, or the sera failed to neutralize. Titers against Omicron are boosted by third vaccine doses and are high in both vaccinated individuals and those infected by Delta. Mutations in Omicron knock out or substantially reduce neutralization by most of the large panel of potent monoclonal antibodies and antibodies under commercial development. Omicron S has structural changes from earlier viruses and uses mutations that confer tight binding to ACE2 to unleash evolution driven by immune escape. This leads to a large number of mutations in the ACE2 binding site and rebalances receptor affinity to that of earlier pandemic viruses.

Evidence of escape of SARS-CoV-2 variant B.1.351 from natural and vaccine-induced sera.

The race to produce vaccines against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) began when the first sequence was published, and this forms the basis for vaccines currently deployed globally. Independent lineages of SARS-CoV-2 have recently been reported: UK, B.1.1.7; South Africa, B.1.351; and Brazil, P.1. These variants have multiple changes in the immunodominant spike protein that facilitates viral cell entry via the angiotensin-converting enzyme-2 (ACE2) receptor. Mutations in the receptor recognition site on the spike are of great concern for their potential for immune escape. Here, we describe a structure-function analysis of B.1.351 using a large cohort of convalescent and vaccinee serum samples. The receptor-binding domain mutations provide tighter ACE2 binding and widespread escape from monoclonal antibody neutralization largely driven by E484K, although K417N and N501Y act together against some important antibody classes. In a number of cases, it would appear that convalescent and some vaccine serum offers limited protection against this variant.

The antigenic anatomy of SARS-CoV-2 receptor binding domain.

Antibodies are crucial to immune protection against SARS-CoV-2, with some in emergency use as therapeutics. Here, we identify 377 human monoclonal antibodies (mAbs) recognizing the virus spike and focus mainly on 80 that bind the receptor binding domain (RBD). We devise a competition data-driven method to map RBD binding sites. We find that although antibody binding sites are widely dispersed, neutralizing antibody binding is focused, with nearly all highly inhibitory mAbs (IC50 < 0.1 µg/mL) blocking receptor interaction, except for one that binds a unique epitope in the N-terminal domain. Many of these neutralizing mAbs use public V-genes and are close to germline. We dissect the structural basis of recognition for this large panel of antibodies through X-ray crystallography and cryoelectron microscopy of 19 Fab-antigen structures. We find novel binding modes for some potently inhibitory antibodies and demonstrate that strongly neutralizing mAbs protect, prophylactically or therapeutically, in animal models.

Reduced neutralization of SARS-CoV-2 B.1.1.7 variant by convalescent and vaccine sera.

SARS-CoV-2 has caused over 2 million deaths in little over a year. Vaccines are being deployed at scale, aiming to generate responses against the virus spike. The scale of the pandemic and error-prone virus replication is leading to the appearance of mutant viruses and potentially escape from antibody responses. Variant B.1.1.7, now dominant in the UK, with increased transmission, harbors 9 amino acid changes in the spike, including N501Y in the ACE2 interacting surface. We examine the ability of B.1.1.7 to evade antibody responses elicited by natural SARS-CoV-2 infection or vaccination. We map the impact of N501Y by structure/function analysis of a large panel of well-characterized monoclonal antibodies. B.1.1.7 is harder to neutralize than parental virus, compromising neutralization by some members of a major class of public antibodies through light-chain contacts with residue 501. However, widespread escape from monoclonal antibodies or antibody responses generated by natural infection or vaccination was not observed.

Antibody evasion by the P.1 strain of SARS-CoV-2.

Terminating the SARS-CoV-2 pandemic relies upon pan-global vaccination. Current vaccines elicit neutralizing antibody responses to the virus spike derived from early isolates. However, new strains have emerged with multiple mutations, including P.1 from Brazil, B.1.351 from South Africa, and B.1.1.7 from the UK (12, 10, and 9 changes in the spike, respectively). All have mutations in the ACE2 binding site, with P.1 and B.1.351 having a virtually identical triplet (E484K, K417N/T, and N501Y), which we show confer similar increased affinity for ACE2. We show that, surprisingly, P.1 is significantly less resistant to naturally acquired or vaccine-induced antibody responses than B.1.351, suggesting that changes outside the receptor-binding domain (RBD) impact neutralization. Monoclonal antibody (mAb) 222 neutralizes all three variants despite interacting with two of the ACE2-binding site mutations. We explain this through structural analysis and use the 222 light chain to largely restore neutralization potency to a major class of public antibodies.

Reduced neutralization of SARS-CoV-2 B.1.617 by vaccine and convalescent serum.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has undergone progressive change, with variants conferring advantage rapidly becoming dominant lineages, e.g., B.1.617. With apparent increased transmissibility, variant B.1.617.2 has contributed to the current wave of infection ravaging the Indian subcontinent and has been designated a variant of concern in the United Kingdom. Here we study the ability of monoclonal antibodies and convalescent and vaccine sera to neutralize B.1.617.1 and B.1.617.2, complement this with structural analyses of Fab/receptor binding domain (RBD) complexes, and map the antigenic space of current variants. Neutralization of both viruses is reduced compared with ancestral Wuhan-related strains, but there is no evidence of widespread antibody escape as seen with B.1.351. However, B.1.351 and P.1 sera showed markedly more reduction in neutralization of B.1.617.2, suggesting that individuals infected previously by these variants may be more susceptible to reinfection by B.1.617.2. This observation provides important new insights for immunization policy with future variant vaccines in non-immune populations.

Dynamic mitochondrial transcription and translation in B cells control germinal center entry and lymphomagenesis.

Germinal center (GC) B cells undergo proliferation at very high rates in a hypoxic microenvironment but the cellular processes driving this are incompletely understood. Here we show that the mitochondria of GC B cells are highly dynamic, with significantly upregulated transcription and translation rates associated with the activity of transcription factor A, mitochondrial (TFAM). TFAM, while also necessary for normal B cell development, is required for entry of activated GC precursor B cells into the germinal center reaction; deletion of Tfam significantly impairs GC formation, function and output. Loss of TFAM in B cells compromises the actin cytoskeleton and impairs cellular motility of GC B cells in response to chemokine signaling, leading to their spatial disorganization. We show that B cell lymphoma substantially increases mitochondrial translation and that deletion of Tfam in B cells is protective against the development of lymphoma in a c-Myc transgenic mouse model. Finally, we show that pharmacological inhibition of mitochondrial transcription and translation inhibits growth of GC-derived human lymphoma cells and induces similar defects in the actin cytoskeleton.

Neutrophils and emergency granulopoiesis drive immune suppression and an extreme response endotype during sepsis.

Sepsis arises from diverse and incompletely understood dysregulated host response processes following infection that leads to life-threatening organ dysfunction. Here we showed that neutrophils and emergency granulopoiesis drove a maladaptive response during sepsis. We generated a whole-blood single-cell multiomic atlas (272,993 cells, n = 39 individuals) of the sepsis immune response that identified populations of immunosuppressive mature and immature neutrophils. In co-culture, CD66b+ sepsis neutrophils inhibited proliferation and activation of CD4+ T cells. Single-cell multiomic mapping of circulating hematopoietic stem and progenitor cells (HSPCs) (29,366 cells, n = 27) indicated altered granulopoiesis in patients with sepsis. These features were enriched in a patient subset with poor outcome and a specific sepsis response signature that displayed higher frequencies of IL1R2+ immature neutrophils, epigenetic and transcriptomic signatures of emergency granulopoiesis in HSPCs and STAT3-mediated gene regulation across different infectious etiologies and syndromes. Our findings offer potential therapeutic targets and opportunities for stratified medicine in severe infection.

Self-Reporting Transposons Enable Simultaneous Readout of Gene Expression and Transcription Factor Binding in Single Cells.

Cellular heterogeneity confounds in situ assays of transcription factor (TF) binding. Single-cell RNA sequencing (scRNA-seq) deconvolves cell types from gene expression, but no technology links cell identity to TF binding sites (TFBS) in those cell types. We present self-reporting transposons (SRTs) and use them in single-cell calling cards (scCC), a novel assay for simultaneously measuring gene expression and mapping TFBS in single cells. The genomic locations of SRTs are recovered from mRNA, and SRTs deposited by exogenous, TF-transposase fusions can be used to map TFBS. We then present scCC, which map SRTs from scRNA-seq libraries, simultaneously identifying cell types and TFBS in those same cells. We benchmark multiple TFs with this technique. Next, we use scCC to discover BRD4-mediated cell-state transitions in K562 cells. Finally, we map BRD4 binding sites in the mouse cortex at single-cell resolution, establishing a new method for studying TF biology in situ.

An immunodominant NP105-113-B*07:02 cytotoxic T cell response controls viral replication and is associated with less severe COVID-19 disease.

NP105-113-B*07:02-specific CD8+ T cell responses are considered among the most dominant in SARS-CoV-2-infected individuals. We found strong association of this response with mild disease. Analysis of NP105-113-B*07:02-specific T cell clones and single-cell sequencing were performed concurrently, with functional avidity and antiviral efficacy assessed using an in vitro SARS-CoV-2 infection system, and were correlated with T cell receptor usage, transcriptome signature and disease severity (acute n = 77, convalescent n = 52). We demonstrated a beneficial association of NP105-113-B*07:02-specific T cells in COVID-19 disease progression, linked with expansion of T cell precursors, high functional avidity and antiviral effector function. Broad immune memory pools were narrowed postinfection but NP105-113-B*07:02-specific T cells were maintained 6 months after infection with preserved antiviral efficacy to the SARS-CoV-2 Victoria strain, as well as Alpha, Beta, Gamma and Delta variants. Our data show that NP105-113-B*07:02-specific T cell responses associate with mild disease and high antiviral efficacy, pointing to inclusion for future vaccine design.

Single-Cell Analysis of Crohn's Disease Lesions Identifies a Pathogenic Cellular Module Associated with Resistance to Anti-TNF Therapy.

Clinical benefits of cytokine blockade in ileal Crohn's disease (iCD) are limited to a subset of patients. Here, we applied single-cell technologies to iCD lesions to address whether cellular heterogeneity contributes to treatment resistance. We found that a subset of patients expressed a unique cellular module in inflamed tissues that consisted of IgG plasma cells, inflammatory mononuclear phagocytes, activated T cells, and stromal cells, which we named the GIMATS module. Analysis of ligand-receptor interaction pairs identified a distinct network connectivity that likely drives the GIMATS module. Strikingly, the GIMATS module was also present in a subset of patients in four independent iCD cohorts (n = 441), and its presence at diagnosis correlated with failure to achieve durable corticosteroid-free remission upon anti-TNF therapy. These results emphasize the limitations of current diagnostic assays and the potential for single-cell mapping tools to identify novel biomarkers of treatment response and tailored therapeutic opportunities.

Natural Killer Cells Degenerate Intact Sensory Afferents following Nerve Injury.

Sensory axons degenerate following separation from their cell body, but partial injury to peripheral nerves may leave the integrity of damaged axons preserved. We show that an endogenous ligand for the natural killer (NK) cell receptor NKG2D, Retinoic Acid Early 1 (RAE1), is re-expressed in adult dorsal root ganglion neurons following peripheral nerve injury, triggering selective degeneration of injured axons. Infiltration of cytotoxic NK cells into the sciatic nerve by extravasation occurs within 3 days following crush injury. Using a combination of genetic cell ablation and cytokine-antibody complex stimulation, we show that NK cell function correlates with loss of sensation due to degeneration of injured afferents and reduced incidence of post-injury hypersensitivity. This neuro-immune mechanism of selective NK cell-mediated degeneration of damaged but intact sensory axons complements Wallerian degeneration and suggests the therapeutic potential of modulating NK cell function to resolve painful neuropathy through the clearance of partially damaged nerves.

Antibody Lineages with Vaccine-Induced Antigen-Binding Hotspots Develop Broad HIV Neutralization.

The vaccine-mediated elicitation of antibodies (Abs) capable of neutralizing diverse HIV-1 strains has been a long-standing goal. To understand how broadly neutralizing antibodies (bNAbs) can be elicited, we identified, characterized, and tracked five neutralizing Ab lineages targeting the HIV-1-fusion peptide (FP) in vaccinated macaques over time. Genetic and structural analyses revealed two of these lineages to belong to a reproducible class capable of neutralizing up to 59% of 208 diverse viral strains. B cell analysis indicated each of the five lineages to have been initiated and expanded by FP-carrier priming, with envelope (Env)-trimer boosts inducing cross-reactive neutralization. These Abs had binding-energy hotspots focused on FP, whereas several FP-directed Abs induced by immunization with Env trimer-only were less FP-focused and less broadly neutralizing. Priming with a conserved subregion, such as FP, can thus induce Abs with binding-energy hotspots coincident with the target subregion and capable of broad neutralization.

An Integrated TCGA Pan-Cancer Clinical Data Resource to Drive High-Quality Survival Outcome Analytics.

For a decade, The Cancer Genome Atlas (TCGA) program collected clinicopathologic annotation data along with multi-platform molecular profiles of more than 11,000 human tumors across 33 different cancer types. TCGA clinical data contain key features representing the democratized nature of the data collection process. To ensure proper use of this large clinical dataset associated with genomic features, we developed a standardized dataset named the TCGA Pan-Cancer Clinical Data Resource (TCGA-CDR), which includes four major clinical outcome endpoints. In addition to detailing major challenges and statistical limitations encountered during the effort of integrating the acquired clinical data, we present a summary that includes endpoint usage recommendations for each cancer type. These TCGA-CDR findings appear to be consistent with cancer genomics studies independent of the TCGA effort and provide opportunities for investigating cancer biology using clinical correlates at an unprecedented scale.

Cell-of-Origin Patterns Dominate the Molecular Classification of 10,000 Tumors from 33 Types of Cancer.

We conducted comprehensive integrative molecular analyses of the complete set of tumors in The Cancer Genome Atlas (TCGA), consisting of approximately 10,000 specimens and representing 33 types of cancer. We performed molecular clustering using data on chromosome-arm-level aneuploidy, DNA hypermethylation, mRNA, and miRNA expression levels and reverse-phase protein arrays, of which all, except for aneuploidy, revealed clustering primarily organized by histology, tissue type, or anatomic origin. The influence of cell type was evident in DNA-methylation-based clustering, even after excluding sites with known preexisting tissue-type-specific methylation. Integrative clustering further emphasized the dominant role of cell-of-origin patterns. Molecular similarities among histologically or anatomically related cancer types provide a basis for focused pan-cancer analyses, such as pan-gastrointestinal, pan-gynecological, pan-kidney, and pan-squamous cancers, and those related by stemness features, which in turn may inform strategies for future therapeutic development.

Machine Learning Identifies Stemness Features Associated with Oncogenic Dedifferentiation.

Cancer progression involves the gradual loss of a differentiated phenotype and acquisition of progenitor and stem-cell-like features. Here, we provide novel stemness indices for assessing the degree of oncogenic dedifferentiation. We used an innovative one-class logistic regression (OCLR) machine-learning algorithm to extract transcriptomic and epigenetic feature sets derived from non-transformed pluripotent stem cells and their differentiated progeny. Using OCLR, we were able to identify previously undiscovered biological mechanisms associated with the dedifferentiated oncogenic state. Analyses of the tumor microenvironment revealed unanticipated correlation of cancer stemness with immune checkpoint expression and infiltrating immune cells. We found that the dedifferentiated oncogenic phenotype was generally most prominent in metastatic tumors. Application of our stemness indices to single-cell data revealed patterns of intra-tumor molecular heterogeneity. Finally, the indices allowed for the identification of novel targets and possible targeted therapies aimed at tumor differentiation.