ePlatypus: an ecosystem for computational analysis of immunogenomics data.

MOTIVATION: The maturation of systems immunology methodologies requires novel and transparent computational frameworks capable of integrating diverse data modalities in a reproducible manner. RESULTS: Here, we present the ePlatypus computational immunology ecosystem for immunogenomics data analysis, with a focus on adaptive immune repertoires and single-cell sequencing. ePlatypus is an open-source web-based platform and provides programming tutorials and an integrative database that helps elucidate signatures of B and T cell clonal selection. Furthermore, the ecosystem links novel and established bioinformatics pipelines relevant for single-cell immune repertoires and other aspects of computational immunology such as predicting ligand-receptor interactions, structural modeling, simulations, machine learning, graph theory, pseudotime, spatial transcriptomics, and phylogenetics. The ePlatypus ecosystem helps extract deeper insight in computational immunology and immunogenomics and promote open science. AVAILABILITY AND IMPLEMENTATION: Platypus code used in this manuscript can be found at github.com/alexyermanos/Platypus.

Deep mutational learning predicts ACE2 binding and antibody escape to combinatorial mutations in the SARS-CoV-2 receptor-binding domain.

The continual evolution of SARS-CoV-2 and the emergence of variants that show resistance to vaccines and neutralizing antibodies threaten to prolong the COVID-19 pandemic. Selection and emergence of SARS-CoV-2 variants are driven in part by mutations within the viral spike protein and in particular the ACE2 receptor-binding domain (RBD), a primary target site for neutralizing antibodies. Here, we develop deep mutational learning (DML), a machine-learning-guided protein engineering technology, which is used to investigate a massive sequence space of combinatorial mutations, representing billions of RBD variants, by accurately predicting their impact on ACE2 binding and antibody escape. A highly diverse landscape of possible SARS-CoV-2 variants is identified that could emerge from a multitude of evolutionary trajectories. DML may be used for predictive profiling on current and prospective variants, including highly mutated variants such as Omicron, thus guiding the development of therapeutic antibody treatments and vaccines for COVID-19.

Learning what not to select for in antibody drug discovery.

Identifying antibodies with high affinity and target specificity is crucial for drug discovery and development; however, filtering out antibody candidates with nonspecific or polyspecific binding profiles is also important. In this issue of Cell Reports Methods, Saksena et al. report a computational counterselection method combining deep sequencing and machine learning for identifying nonspecific antibody candidates and demonstrate that it has advantages over more established molecular counterselection methods.

DeepSARS: simultaneous diagnostic detection and genomic surveillance of SARS-CoV-2.

BACKGROUND: The continued spread of SARS-CoV-2 and emergence of new variants with higher transmission rates and/or partial resistance to vaccines has further highlighted the need for large-scale testing and genomic surveillance. However, current diagnostic testing (e.g., PCR) and genomic surveillance methods (e.g., whole genome sequencing) are performed separately, thus limiting the detection and tracing of SARS-CoV-2 and emerging variants. RESULTS: Here, we developed DeepSARS, a high-throughput platform for simultaneous diagnostic detection and genomic surveillance of SARS-CoV-2 by the integration of molecular barcoding, targeted deep sequencing, and computational phylogenetics. DeepSARS enables highly sensitive viral detection, while also capturing genomic diversity and viral evolution. We show that DeepSARS can be rapidly adapted for identification of emerging variants, such as alpha, beta, gamma, and delta strains, and profile mutational changes at the population level. CONCLUSIONS: DeepSARS sets the foundation for quantitative diagnostics that capture viral evolution and diversity. DeepSARS uses molecular barcodes (BCs) and multiplexed targeted deep sequencing (NGS) to enable simultaneous diagnostic detection and genomic surveillance of SARS-CoV-2. Image was created using Biorender.com .

Profiling the specificity of clonally expanded plasma cells during chronic viral infection by single-cell analysis.

Plasma cells and their secreted antibodies play a central role in the long-term protection against chronic viral infection. However, due to experimental limitations, a comprehensive description of linked genotypic, phenotypic, and antibody repertoire features of plasma cells (gene expression, clonal frequency, virus specificity, and affinity) has been challenging to obtain. To address this, we performed single-cell transcriptome and antibody repertoire sequencing of the murine BM plasma cell population following chronic lymphocytic choriomeningitis virus infection. Our single-cell sequencing approach recovered full-length and paired heavy- and light-chain sequence information for thousands of plasma cells and enabled us to perform recombinant antibody expression and specificity screening. Antibody repertoire analysis revealed that, relative to protein immunization, chronic infection led to increased levels of clonal expansion, class-switching, and somatic variants. Furthermore, antibodies from the highly expanded and class-switched (IgG) plasma cells were found to be specific for multiple viral antigens and a subset of clones exhibited cross-reactivity to nonviral and autoantigens. Integrating single-cell transcriptome data with antibody specificity suggested that plasma cell transcriptional phenotype was correlated to viral antigen specificity. Our findings demonstrate that chronic viral infection can induce and sustain plasma cell clonal expansion, combined with significant somatic hypermutation, and can generate cross-reactive antibodies.

Echidna: integrated simulations of single-cell immune receptor repertoires and transcriptomes.

Motivation: Single-cell sequencing now enables the recovery of full-length immune receptor repertoires [B cell receptor (BCR) and T cell receptor (TCR) repertoires], in addition to gene expression information. The feature-rich datasets produced from such experiments require extensive and diverse computational analyses, each of which can significantly influence the downstream immunological interpretations, such as clonal selection and expansion. Simulations produce validated standard datasets, where the underlying generative model can be precisely defined and furthermore perturbed to investigate specific questions of interest. Currently, there is no tool that can be used to simulate single-cell datasets incorporating immune receptor repertoires and gene expression. Results: We developed Echidna, an R package that simulates immune receptors and transcriptomes at single-cell resolution with user-tunable parameters controlling a wide range of features such as clonal expansion, germline gene usage, somatic hypermutation, transcriptional phenotypes and spatial location. Echidna can additionally simulate time-resolved B cell evolution, producing mutational networks with complex selection histories incorporating class-switching and B cell subtype information. We demonstrated the benchmarking potential of Echidna by simulating clonal lineages and comparing the known simulated networks with those inferred from only the BCR sequences as input. Finally, we simulated immune repertoire information onto existing spatial transcriptomic experiments, thereby generating novel datasets that could be used to develop and integrate methods to profile clonal selection in a spatially resolved manner. Together, Echidna provides a framework that can incorporate experimental data to simulate single-cell immune repertoires to aid software development and bioinformatic benchmarking of clonotyping, phylogenetics, transcriptomics and machine learning strategies. Availability and implementation: The R package and code used in this manuscript can be found at github.com/alexyermanos/echidna and also in the R package Platypus (Yermanos et al., 2021). Installation instructions and the vignette for Echidna is described in the Platypus Computational Ecosystem (https://alexyermanos.github.io/Platypus/index.html). Publicly available data and corresponding sample accession numbers can be found in Supplementary Tables S2 and S3. Supplementary information: Supplementary data are available at Bioinformatics Advances online.

Platypus: an open-access software for integrating lymphocyte single-cell immune repertoires with transcriptomes.

High-throughput single-cell sequencing (scSeq) technologies are revolutionizing the ability to molecularly profile B and T lymphocytes by offering the opportunity to simultaneously obtain information on adaptive immune receptor repertoires (VDJ repertoires) and transcriptomes. An integrated quantification of immune repertoire parameters, such as germline gene usage, clonal expansion, somatic hypermutation and transcriptional states opens up new possibilities for the high-resolution analysis of lymphocytes and the inference of antigen-specificity. While multiple tools now exist to investigate gene expression profiles from scSeq of transcriptomes, there is a lack of software dedicated to single-cell immune repertoires. Here, we present Platypus, an open-source software platform providing a user-friendly interface to investigate B-cell receptor and T-cell receptor repertoires from scSeq experiments. Platypus provides a framework to automate and ease the analysis of single-cell immune repertoires while also incorporating transcriptional information involving unsupervised clustering, gene expression and gene ontology. To showcase the capabilities of Platypus, we use it to analyze and visualize single-cell immune repertoires and transcriptomes from B and T cells from convalescent COVID-19 patients, revealing unique insight into the repertoire features and transcriptional profiles of clonally expanded lymphocytes. Platypus will expedite progress by facilitating the analysis of single-cell immune repertoire and transcriptome sequencing.

TLR7 Signaling Shapes and Maintains Antibody Diversity Upon Virus-Like Particle Immunization.

Virus-like particles (VLPs) are used in different marketed vaccines and are able to induce potent antibody responses. The innate pattern recognition receptors TLR7/8 recognize single stranded (ss) RNA naturally packaged into some VLPs and have been shown to enhance the production of IgG antibodies upon immunization. Here we demonstrate that, upon immunization with RNA-loaded bacteriophage-derived VLP Qß, TLR7 signaling accelerates germinal center formation, promotes affinity/avidity maturation of VLP-specific IgG and isotype switching to IgG2b/2c. These findings extrapolated to antigens displayed on Qß; as Fel d 1, the major cat allergen, chemically attached to Qß also induced higher affinity/avidity IgG2b/2c antibodies in a TLR7-dependent fashion. Chimeric mice lacking TLR7-expression exclusively in B cells demonstrated that the enhanced IgG responses were driven by a B cell intrinsic mechanism. Importantly, deep sequencing of the BCR repertoire of antigen-specific B cells demonstrated higher diversity in mice with TLR7 signaling in B cells, suggesting that TLR7-signaling drives BCR repertoire development and diversity. Furthermore, the current data demonstrate that high levels of clonal diversity are reached early in the response and maintained by TLR7 signaling. In conclusion, TLR7 signaling enhances levels and quality of IgG antibodies, and this finding has major implications for vaccine design.

Simvastatin functions as a heat shock protein 90 inhibitor against triple-negative breast cancer.

Acetylation plays an important role in regulating the chaperone activity of heat shock protein 90 (Hsp90) during malignant transformation through the stabilization and conformational maturation of oncogenic proteins. However, the functional acetylation sites, potential anticancer drug targets, are still emerging. We found that acetylation at K292 in Hsp90α is critical for the development and treatment of breast cancer. Acetylation at K292 not only augments the affinity of Hsp90 to ATP, cochaperones, and client proteins but it also promotes cancer cell colony formation, migration, and invasion in vitro as well as tumor growth in vivo. Importantly, K292-acetylated Hsp90 has been validated as an exciting anticancer drug target by interfering with the complex formation between K292-acetylated Hsp90 and cochaperone Cdc37, leading to diminishment of kinase client maturation and proteasome-dependent degradation of kinase substrates. Furthermore, we showed that simvastatin prevented, whereas LBH589 promoted, the progression of Hsp90 chaperone cycling and client maturation, resulting in an increment of cell apoptosis by the combination of simvastatin and LBH589 in a mouse xenograft model. These data suggest that simvastatin is a novel Hsp90 inhibitor to disrupt the formation of the K292-acetylated Hsp90/Cdc37 complex in triple-negative breast cancer cells. The combination of simvastatin with LBH589 could be used as a novel therapeutic strategy for triple-negative breast cancer.