Pharmacologic modulation of RNA splicing enhances anti-tumor immunity.

Although mutations in DNA are the best-studied source of neoantigens that determine response to immune checkpoint blockade, alterations in RNA splicing within cancer cells could similarly result in neoepitope production. However, the endogenous antigenicity and clinical potential of such splicing-derived epitopes have not been tested. Here, we demonstrate that pharmacologic modulation of splicing via specific drug classes generates bona fide neoantigens and elicits anti-tumor immunity, augmenting checkpoint immunotherapy. Splicing modulation inhibited tumor growth and enhanced checkpoint blockade in a manner dependent on host T cells and peptides presented on tumor MHC class I. Splicing modulation induced stereotyped splicing changes across tumor types, altering the MHC I-bound immunopeptidome to yield splicing-derived neoepitopes that trigger an anti-tumor T cell response in vivo. These data definitively identify splicing modulation as an untapped source of immunogenic peptides and provide a means to enhance response to checkpoint blockade that is readily translatable to the clinic.

Neoantigen quality predicts immunoediting in survivors of pancreatic cancer.

Cancer immunoediting1 is a hallmark of cancer2 that predicts that lymphocytes kill more immunogenic cancer cells to cause less immunogenic clones to dominate a population. Although proven in mice1,3, whether immunoediting occurs naturally in human cancers remains unclear. Here, to address this, we investigate how 70 human pancreatic cancers evolved over 10 years. We find that, despite having more time to accumulate mutations, rare long-term survivors of pancreatic cancer who have stronger T cell activity in primary tumours develop genetically less heterogeneous recurrent tumours with fewer immunogenic mutations (neoantigens). To quantify whether immunoediting underlies these observations, we infer that a neoantigen is immunogenic (high-quality) by two features-'non-selfness'  based on neoantigen similarity to known antigens4,5, and 'selfness'  based on the antigenic distance required for a neoantigen to differentially bind to the MHC or activate a T cell compared with its wild-type peptide. Using these features, we estimate cancer clone fitness as the aggregate cost of T cells recognizing high-quality neoantigens offset by gains from oncogenic mutations. With this model, we predict the clonal evolution of tumours to reveal that long-term survivors of pancreatic cancer develop recurrent tumours with fewer high-quality neoantigens. Thus, we submit evidence that that the human immune system naturally edits neoantigens. Furthermore, we present a model to predict how immune pressure induces cancer cell populations to evolve over time. More broadly, our results argue that the immune system fundamentally surveils host genetic changes to suppress cancer.

Predicting the spectrum of TCR repertoire sharing with a data-driven model of recombination.

Despite the extreme diversity of T-cell repertoires, many identical T-cell receptor (TCR) sequences are found in a large number of individual mice and humans. These widely shared sequences, often referred to as "public," have been suggested to be over-represented due to their potential immune functionality or their ease of generation by V(D)J recombination. Here, we show that even for large cohorts, the observed degree of sharing of TCR sequences between individuals is well predicted by a model accounting for the known quantitative statistical biases in the generation process, together with a simple model of thymic selection. Whether a sequence is shared by many individuals is predicted to depend on the number of queried individuals and the sampling depth, as well as on the sequence itself, in agreement with the data. We introduce the degree of publicness conditional on the queried cohort size and the size of the sampled repertoires. Based on these observations, we propose a public/private sequence classifier, "PUBLIC" (Public Universal Binary Likelihood Inference Classifier), based on the generation probability, which performs very well even for small cohort sizes.

Population variability in the generation and selection of T-cell repertoires.

The diversity of T-cell receptor (TCR) repertoires is achieved by a combination of two intrinsically stochastic steps: random receptor generation by VDJ recombination, and selection based on the recognition of random self-peptides presented on the major histocompatibility complex. These processes lead to a large receptor variability within and between individuals. However, the characterization of the variability is hampered by the limited size of the sampled repertoires. We introduce a new software tool SONIA to facilitate inference of individual-specific computational models for the generation and selection of the TCR beta chain (TRB) from sequenced repertoires of 651 individuals, separating and quantifying the variability of the two processes of generation and selection in the population. We find not only that most of the variability is driven by the VDJ generation process, but there is a large degree of consistency between individuals with the inter-individual variance of repertoires being about ∼2% of the intra-individual variance. Known viral-specific TCRs follow the same generation and selection statistics as all TCRs.

OLGA: fast computation of generation probabilities of B- and T-cell receptor amino acid sequences and motifs.

MOTIVATION: High-throughput sequencing of large immune repertoires has enabled the development of methods to predict the probability of generation by V(D)J recombination of T- and B-cell receptors of any specific nucleotide sequence. These generation probabilities are very non-homogeneous, ranging over 20 orders of magnitude in real repertoires. Since the function of a receptor really depends on its protein sequence, it is important to be able to predict this probability of generation at the amino acid level. However, brute-force summation over all the nucleotide sequences with the correct amino acid translation is computationally intractable. The purpose of this paper is to present a solution to this problem. RESULTS: We use dynamic programming to construct an efficient and flexible algorithm, called OLGA (Optimized Likelihood estimate of immunoGlobulin Amino-acid sequences), for calculating the probability of generating a given CDR3 amino acid sequence or motif, with or without V/J restriction, as a result of V(D)J recombination in B or T cells. We apply it to databases of epitope-specific T-cell receptors to evaluate the probability that a typical human subject will possess T cells responsive to specific disease-associated epitopes. The model prediction shows an excellent agreement with published data. We suggest that OLGA may be a useful tool to guide vaccine design. AVAILABILITY AND IMPLEMENTATION: Source code is available at https://github.com/zsethna/OLGA. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Genetic and environmental determinants of human TCR repertoire diversity.

T cell discrimination of self and non-self is the foundation of the adaptive immune response, and is orchestrated by the interaction between T cell receptors (TCRs) and their cognate ligands presented by major histocompatibility (MHC) molecules. However, the impact of host immunogenetic variation on the diversity of the TCR repertoire remains unclear. Here, we analyzed a cohort of 666 individuals with TCR repertoire sequencing. We show that TCR repertoire diversity is positively associated with polymorphism at the human leukocyte antigen class I (HLA-I) loci, and diminishes with age and cytomegalovirus (CMV) infection. Moreover, our analysis revealed that HLA-I polymorphism and age independently shape the repertoire in healthy individuals. Our data elucidate key determinants of human TCR repertoire diversity, and suggest a mechanism underlying the evolutionary fitness advantage of HLA-I heterozygosity.

Insights into immune system development and function from mouse T-cell repertoires.

The ability of the adaptive immune system to respond to arbitrary pathogens stems from the broad diversity of immune cell surface receptors. This diversity originates in a stochastic DNA editing process (VDJ recombination) that acts on the surface receptor gene each time a new immune cell is created from a stem cell. By analyzing T-cell receptor (TCR) sequence repertoires taken from the blood and thymus of mice of different ages, we quantify the changes in the VDJ recombination process that occur from embryo to young adult. We find a rapid increase with age in the number of random insertions and a dramatic increase in diversity. Because the blood accumulates thymic output over time, blood repertoires are mixtures of different statistical recombination processes, and we unravel the mixture statistics to obtain a picture of the time evolution of the early immune system. Sequence repertoire analysis also allows us to detect the statistical impact of selection on the output of the VDJ recombination process. The effects we find are nearly identical between thymus and blood, suggesting that our analysis mainly detects selection for proper folding of the TCR receptor protein. We further find that selection is weaker in laboratory mice than in humans and it does not affect the diversity of the repertoire.

Persisting fetal clonotypes influence the structure and overlap of adult human T cell receptor repertoires.

The diversity of T-cell receptors recognizing foreign pathogens is generated through a highly stochastic recombination process, making the independent production of the same sequence rare. Yet unrelated individuals do share receptors, which together constitute a "public" repertoire of abundant clonotypes. The TCR repertoire is initially formed prenatally, when the enzyme inserting random nucleotides is downregulated, producing a limited diversity subset. By statistically analyzing deep sequencing T-cell repertoire data from twins, unrelated individuals of various ages, and cord blood, we show that T-cell clones generated before birth persist and maintain high abundances in adult organisms for decades, slowly decaying with age. Our results suggest that large, low-diversity public clones are created during pre-natal life, and survive over long periods, providing the basis of the public repertoire.

repgenHMM: a dynamic programming tool to infer the rules of immune receptor generation from sequence data.

MOTIVATION: The diversity of the immune repertoire is initially generated by random rearrangements of the receptor gene during early T and B cell development. Rearrangement scenarios are composed of random events-choices of gene templates, base pair deletions and insertions-described by probability distributions. Not all scenarios are equally likely, and the same receptor sequence may be obtained in several different ways. Quantifying the distribution of these rearrangements is an essential baseline for studying the immune system diversity. Inferring the properties of the distributions from receptor sequences is a computationally hard problem, requiring enumerating every possible scenario for every sampled receptor sequence. RESULTS: We present a Hidden Markov model, which accounts for all plausible scenarios that can generate the receptor sequences. We developed and implemented a method based on the Baum-Welch algorithm that can efficiently infer the parameters for the different events of the rearrangement process. We tested our software tool on sequence data for both the alpha and beta chains of the T cell receptor. To test the validity of our algorithm, we also generated synthetic sequences produced by a known model, and confirmed that its parameters could be accurately inferred back from the sequences. The inferred model can be used to generate synthetic sequences, to calculate the probability of generation of any receptor sequence, as well as the theoretical diversity of the repertoire. We estimate this diversity to be [Formula: see text] for human T cells. The model gives a baseline to investigate the selection and dynamics of immune repertoires. AVAILABILITY AND IMPLEMENTATION: Source code and sample sequence files are available at https://bitbucket.org/yuvalel/repgenhmm/downloads CONTACT: elhanati@lpt.ens.fr or tmora@lps.ens.fr or awalczak@lpt.ens.fr.

Quantifying selection in immune receptor repertoires.

The efficient recognition of pathogens by the adaptive immune system relies on the diversity of receptors displayed at the surface of immune cells. T-cell receptor diversity results from an initial random DNA editing process, called VDJ recombination, followed by functional selection of cells according to the interaction of their surface receptors with self and foreign antigenic peptides. Using high-throughput sequence data from the ß-chain of human T-cell receptors, we infer factors that quantify the overall effect of selection on the elements of receptor sequence composition: the V and J gene choice and the length and amino acid composition of the variable region. We find a significant correlation between biases induced by VDJ recombination and our inferred selection factors together with a reduction of diversity during selection. Both effects suggest that natural selection acting on the recombination process has anticipated the selection pressures experienced during somatic evolution. The inferred selection factors differ little between donors or between naive and memory repertoires. The number of sequences shared between donors is well-predicted by our model, indicating a stochastic origin of such public sequences. Our approach is based on a probabilistic maximum likelihood method, which is necessary to disentangle the effects of selection from biases inherent in the recombination process.

Preoperative immune checkpoint inhibition and cryoablation in early-stage breast cancer.

Local cryoablation can engender systemic immune activation/anticancer responses in tumors otherwise resistant to immune checkpoint blockade (ICB). We evaluated the safety/tolerability of preoperative cryoablation plus ipilimumab and nivolumab in 5 early-stage/resectable breast cancers. The primary endpoint was met when all 5 patients underwent standard-of-care primary breast surgery undelayedly. Three patients developed transient hyperthyroidism; one developed grade 4 liver toxicity (resolved with supportive management). We compared this strategy with cryoablation and/or ipilimumab. Dual ICB plus cryoablation induced higher expression of T cell activation markers and serum Th1 cytokines and reduced immunosuppressive serum CD4+PD-1hi T cells, improving effector-to-suppressor T cell ratio. After dual ICB and before cryoablation, T cell receptor sequencing of 4 patients showed increased T cell clonality. In this small subset of patients, we provide preliminary evidence that preoperative cryoablation plus ipilimumab and nivolumab is feasible, inducing systemic adaptive immune activation potentially more robust than cryoablation with/without ipilimumab.

Dynamics of the cell-cycle network under genome-rewiring perturbations.

The cell-cycle progression is regulated by a specific network enabling its ordered dynamics. Recent experiments supported by computational models have shown that a core of genes ensures this robust cycle dynamics. However, much less is known about the direct interaction of the cell-cycle regulators with genes outside of the cell-cycle network, in particular those of the metabolic system. Following our recent experimental work, we present here a model focusing on the dynamics of the cell-cycle core network under rewiring perturbations. Rewiring is achieved by placing an essential metabolic gene exclusively under the regulation of a cell-cycle's promoter, forcing the cell-cycle network to function under a multitasking challenging condition; operating in parallel the cell-cycle progression and a metabolic essential gene. Our model relies on simple rate equations that capture the dynamics of the relevant protein-DNA and protein-protein interactions, while making a clear distinction between these two different types of processes. In particular, we treat the cell-cycle transcription factors as limited 'resources' and focus on the redistribution of resources in the network during its dynamics. This elucidates the sensitivity of its various nodes to rewiring interactions. The basic model produces the correct cycle dynamics for a wide range of parameters. The simplicity of the model enables us to study the interface between the cell-cycle regulation and other cellular processes. Rewiring a promoter of the network to regulate a foreign gene, forces a multitasking regulatory load. The higher the load on the promoter, the longer is the cell-cycle period. Moreover, in agreement with our experimental results, the model shows that different nodes of the network exhibit variable susceptibilities to the rewiring perturbations. Our model suggests that the topology of the cell-cycle core network ensures its plasticity and flexible interface with other cellular processes, without a need for an optimal setting of the kinetic parameters.

High-throughput interrogation of immune responses using the Human Immune Profiling Pipeline.

The current abundance of immunotherapy clinical trials presents an opportunity to learn about the underlying mechanisms and pharmacodynamic effects of novel drugs on the human immune system. Here, we present a protocol to study how these immune responses impact clinical outcomes using large-scale high-throughput immune profiling of clinical cohorts. We describe the Human Immune Profiling Pipeline, which comprises an end-to-end solution from flow cytometry results to computational approaches and unsupervised patient clustering based on lymphocyte landscape. For complete details on the use and execution of this protocol, please refer to Lyudovyk et al. (2022).1.

DRAG in situ barcoding reveals an increased number of HSPCs contributing to myelopoiesis with age.

Ageing is associated with changes in the cellular composition of the immune system. During ageing, hematopoietic stem and progenitor cells (HSPCs) that produce immune cells are thought to decline in their regenerative capacity. However, HSPC function has been mostly assessed using transplantation assays, and it remains unclear how HSPCs age in the native bone marrow niche. To address this issue, we present an in situ single cell lineage tracing technology to quantify the clonal composition and cell production of single cells in their native niche. Our results demonstrate that a pool of HSPCs with unequal output maintains myelopoiesis through overlapping waves of cell production throughout adult life. During ageing, the increased frequency of myeloid cells is explained by greater numbers of HSPCs contributing to myelopoiesis rather than the increased myeloid output of individual HSPCs. Strikingly, the myeloid output of HSPCs remains constant over time despite accumulating significant transcriptomic changes throughout adulthood. Together, these results show that, unlike emergency myelopoiesis post-transplantation, aged HSPCs in their native microenvironment do not functionally decline in their regenerative capacity.

Tissue-specific features of the T cell repertoire after allogeneic hematopoietic cell transplantation in human and mouse.

T cells are the central drivers of many inflammatory diseases, but the repertoire of tissue-resident T cells at sites of pathology in human organs remains poorly understood. We examined the site-specificity of T cell receptor (TCR) repertoires across tissues (5 to 18 tissues per patient) in prospectively collected autopsies of patients with and without graft-versus-host disease (GVHD), a potentially lethal tissue-targeting complication of allogeneic hematopoietic cell transplantation, and in mouse models of GVHD. Anatomic similarity between tissues was a key determinant of TCR repertoire composition within patients, independent of disease or transplant status. The T cells recovered from peripheral blood and spleens in patients and mice captured a limited portion of the TCR repertoire detected in tissues. Whereas few T cell clones were shared across patients, motif-based clustering revealed shared repertoire signatures across patients in a tissue-specific fashion. T cells at disease sites had a tissue-resident phenotype and were of donor origin based on single-cell chimerism analysis. These data demonstrate the complex composition of T cell populations that persist in human tissues at the end stage of an inflammatory disorder after lymphocyte-directed therapy. These findings also underscore the importance of studying T cell in tissues rather than blood for tissue-based pathologies and suggest the tissue-specific nature of both the endogenous and posttransplant T cell landscape.

Investigative needle core biopsies for multi-omics in Glioblastoma.

Glioblastoma (GBM) is a primary brain cancer with an abysmal prognosis and few effective therapies. The ability to investigate the tumor microenvironment before and during treatment would greatly enhance both understanding of disease response and progression, as well as the delivery and impact of therapeutics. Stereotactic biopsies are a routine surgical procedure performed primarily for diagnostic histopathologic purposes. The role of investigative biopsies - tissue sampling for the purpose of understanding tumor microenvironmental responses to treatment using integrated multi-modal molecular analyses ('Multi-omics") has yet to be defined. Secondly, it is unknown whether comparatively small tissue samples from brain biopsies can yield sufficient information with such methods. Here we adapt stereotactic needle core biopsy tissue in two separate patients. In the first patient with recurrent GBM we performed highly resolved multi-omics analysis methods including single cell RNA sequencing, spatial-transcriptomics, metabolomics, proteomics, phosphoproteomics, T-cell clonotype analysis, and MHC Class I immunopeptidomics from biopsy tissue that was obtained from a single procedure. In a second patient we analyzed multi-regional core biopsies to decipher spatial and genomic variance. We also investigated the utility of stereotactic biopsies as a method for generating patient derived xenograft models in a separate patient cohort. Dataset integration across modalities showed good correspondence between spatial modalities, highlighted immune cell associated metabolic pathways and revealed poor correlation between RNA expression and the tumor MHC Class I immunopeptidome. In conclusion, stereotactic needle biopsy cores are of sufficient quality to generate multi-omics data, provide data rich insight into a patient's disease process and tumor immune microenvironment and can be of value in evaluating treatment responses. One sentence summary: Integrative multi-omics analysis of stereotactic needle core biopsies in glioblastoma.

Langerhans dendritic cell vaccine bearing mRNA-encoded tumor antigens induces antimyeloma immunity after autotransplant.

Posttransplant vaccination targeting residual disease is an immunotherapeutic strategy to improve antigen-specific immune responses and prolong disease-free survival after autologous stem cell transplantation (ASCT) for multiple myeloma (MM). We conducted a phase 1 vaccine trial to determine the safety, toxicity, and immunogenicity of autologous Langerhans-type dendritic cells (LCs) electroporated with CT7, MAGE-A3, and Wilms tumor 1 (WT1) messenger RNA (mRNA), after ASCT for MM. Ten patients received a priming immunization plus 2 boosters at 12, 30, and 90 days, respectively, after ASCT. Vaccines contained 9 × 106 mRNA-electroporated LCs. Ten additional patients did not receive LC vaccines but otherwise underwent identical ASCT and supportive care. At 3 months after ASCT, all patients started lenalidomide maintenance therapy. Vaccinated patients developed mild local delayed-type hypersensitivity reactions after booster vaccines, but no toxicities exceeded grade 1. At 1 and 3 months after vaccines, antigen-specific CD4 and CD8 T cells increased secretion of proinflammatory cytokines (interferon-γ, interleukin-2, and tumor necrosis factor-α) above prevaccine levels, and also upregulated the cytotoxicity marker CD107a. CD4 and CD8 T-cell repertoire analysis showed a trend for increased clonal expansion in the vaccine cohort, which was more pronounced in the CD4 compartment. Although not powered to assess clinical efficacy, treatment responses favored the vaccine arm. Triple antigen-bearing mRNA-electroporated autologous LC vaccination initiated at engraftment after ASCT, in conjunction with standard lenalidomide maintenance therapy for MM, is safe and induces antigen-specific immune reactivity. This trial was registered at www.clinicaltrials.gov as #NCT01995708.

Impaired humoral immunity is associated with prolonged COVID-19 despite robust CD8 T cell responses.

How immune dysregulation affects recovery from COVID-19 infection in patients with cancer remains unclear. We analyzed cellular and humoral immune responses in 103 patients with prior COVID-19 infection, more than 20% of whom had delayed viral clearance. Delayed clearance was associated with loss of antibodies to nucleocapsid and spike proteins with a compensatory increase in functional T cell responses. High-dimensional analysis of peripheral blood samples demonstrated increased CD8+ effector T cell differentiation and a broad but poorly converged COVID-specific T cell receptor (TCR) repertoire in patients with prolonged disease. Conversely, patients with a CD4+ dominant immunophenotype had a lower incidence of prolonged disease and exhibited a deep and highly select COVID-associated TCR repertoire, consistent with effective viral clearance and development of T cell memory. These results highlight the importance of B cells and CD4+ T cells in promoting durable SARS-CoV-2 clearance and the significance of coordinated cellular and humoral immunity for long-term disease control.

Dynamic modeling of cooperative protein secretion in microorganism populations.

Interactions between microorganisms can have a crucial effect on their population dynamics. Typically, interactions are mediated through the environment by molecules and proteins that are products of cell metabolism and physiology; they therefore reflect the internal dynamics of the single cell. In this work we aim to integrate single-cell properties of gene expression that affect indirect interactions between microorganisms under challenging conditions, into a quantitative model of population dynamics. Specifically we address the problem of a microbial population secreting a protein that can actively extract a growth-limiting resource, such as a simple sugar or iron, from the environment. The genes coding for the protein can undergo random epigenetic transitions between active and silenced states, and can be repressed by the product of their reaction. We model cooperative and competitive interactions between protein producing and non-producing phenotypes by nonlinear dynamical systems and analyze them both in terms of asymptotic states and of transient dynamics. Our model shows that phenotypic transitions allow a stable coexistence of the two phenotypes, and enables us to make predictions regarding the conditions required for such coexistence and the typical timescales of transient dynamics. It also shows how repression by the reaction product induces a feedback at the population-environment level that can result in limit cycle dynamics. The relation of these results to experiments are discussed.

Generative models of T-cell receptor sequences.

T-cell receptors (TCR) are key proteins of the adaptive immune system, generated randomly in each individual, whose diversity underlies our ability to recognize infections and malignancies. Modeling the distribution of TCR sequences is of key importance for immunology and medical applications. Here, we compare two inference methods trained on high-throughput sequencing data: a knowledge-guided approach, which accounts for the details of sequence generation, supplemented by a physics-inspired model of selection; and a knowledge-free variational autoencoder based on deep artificial neural networks. We show that the knowledge-guided model outperforms the deep network approach at predicting TCR probabilities, while being more interpretable, at a lower computational cost.