B cell clonality in cancer

Carcinogenesis in the process of long-term co-evolution of tumor cells and immune environment essentially becomes possible due to incorrect decisions made, remembered, and reproduced by the immune system at the level of clonal populations of antigen-specific T- and B-lymphocytes. Tumor-immunity interaction determines the nature of such errors and, consequently, delineates the possible ways of successful immunotherapeutic intervention. It is generally recognized that tumor-infiltrating B cells (TIL-B) can play both pro-tumor and anti-tumor roles. However, the exact mechanisms that determine the contribution of clonal B cell lineages with different specificities and functions remain largely unclear. This is due to the variability of cancer types, the molecular heterogeneity of tumor cells, and, to a large extent, the individual pattern of each immune response. Further progress requires detailed investigation of the functional properties and phenotypes of clonally heterogeneous B cells in relation to their antigenic specificities, which determine the functionality of both effector B lymphocytes and immunoglobulins produced in the tumor environment. Based on a real understanding of the role of clonal antigen-specific populations of B lymphocytes in the tumor microenvironment, we need to learn how to develop new methods of targeted immunotherapy, as well as adapt existing treatment options to the specific needs of different patients and patient subgroups. In this review, we will cover B cells functional diversity and their multifaceted roles in the tumor environment.

Toolkit for mapping the clonal landscape of tumor-infiltrating B cells

Our current understanding of whether B cell involvement in the tumor microenvironment benefits the patient or the tumor - in distinct cancers, subcohorts and individual patients - is quite limited. Both statements are probably true in most cases: certain clonal B cell populations contribute to the antitumor response, while others steer the immune response away from the desired mechanics. To step up to a new level of understanding and managing B cell behaviors in the tumor microenvironment, we need to rationally discern these roles, which are cumulatively defined by B cell clonal functional programs, specificities of their B cell receptors, specificities and isotypes of the antibodies they produce, and their spatial interactions within the tumor environment. Comprehensive analysis of these characteristics of clonal B cell populations is now becoming feasible with the development of a whole arsenal of advanced technical approaches, which include (1) methods of single-cell and spatial transcriptomics, genomics, and proteomics; (2) methods of massive identification of B cell specificities; (3) methods of deep error-free profiling of B cell receptor repertoires. Here we overview existing techniques, summarize their current application for B cells studies and propose promising future directions in advancing B cells exploration.

Memory persistence and differentiation into antibody-secreting cells accompanied by positive selection in longitudinal BCR repertoires

The stability and plasticity of B cell-mediated immune memory ensures the ability to respond to the repeated challenges. We have analyzed the longitudinal dynamics of immunoglobulin heavy chain repertoires from memory B cells, plasmablasts, and plasma cells from the peripheral blood of generally healthy volunteers. We reveal a high degree of clonal persistence in individual memory B cell subsets, with inter-individual convergence in memory and antibody-secreting cells (ASCs). ASC clonotypes demonstrate clonal relatedness to memory B cells, and are transient in peripheral blood. We identify two clusters of expanded clonal lineages with differing prevalence of memory B cells, isotypes, and persistence. Phylogenetic analysis revealed signs of reactivation of persisting memory B cell-enriched clonal lineages, accompanied by new rounds of affinity maturation during proliferation and differentiation into ASCs. Negative selection contributes to both persisting and reactivated lineages, preserving the functionality and specificity of B cell receptors (BCRs) to protect against current and future pathogens.

Accounting for B-cell Behavior and Sampling Bias Predicts Anti-PD-L1 Response in Bladder Cancer

Cancer immunotherapy is predominantly based on T cell-centric approaches. At the same time, the adaptive immune response in the tumor environment also includes clonally produced immunoglobulins and clonal effector/memory B cells that participate in antigen-specific decisions through their interactions with T cells. Here, we investigated the role of infiltrating B cells in bladder cancer via patient dataset analysis of intratumoral immunoglobulin repertoires. We showed that the IgG1/IgA ratio is a prognostic indicator for several subtypes of bladder cancer and for the whole IMVigor210 anti-PD-L1 immunotherapy study cohort. A high IgG1/IgA ratio associated with the prominence of a cytotoxic gene signature, T-cell receptor signaling, and IL21-mediated signaling. Immunoglobulin repertoire analysis indicated that effector B-cell function, rather than clonally produced antibodies, was involved in antitumor responses. From the T-cell side, we normalized a cytotoxic signature against the extent of immune cell infiltration to neutralize the artificial sampling-based variability in immune gene expression. Resulting metrics reflected proportion of cytotoxic cells among tumor-infiltrating immune cells and improved prediction of anti-PD-L1 responses. At the same time, the IgG1/IgA ratio remained an independent prognostic factor. Integration of the B-cell, natural killer cell, and T-cell signatures allowed for the most accurate prediction of anti-PD-L1 therapy responses. On the basis of these findings, we developed a predictor called PRedIctive MolecUlar Signature (PRIMUS), which outperformed PD-L1 expression scores and known gene signatures. Overall, PRIMUS allows for reliable identification of responders among patients with muscle-invasive urothelial carcinoma, including the subcohort with the low-infiltrated "desert" tumor phenotype.

RNA-Seq-Based TCR Profiling Reveals Persistently Increased Intratumoral Clonality in Responders to Anti-PD-1 Therapy

Substantial effort is being invested in the search for peripheral or intratumoral T cell receptor (TCR) repertoire features that could predict the response to immunotherapy. Here we demonstrate the utility of MiXCR software for TCR and immunoglobulin repertoire extraction from RNA-Seq data obtained from sorted tumor-infiltrating T and B cells. We use this approach to extract TCR repertoires from RNA-Seq data obtained from sorted tumor-infiltrating CD4⁺ and CD8⁺ T cells in an HKP1 (Kras^G12Dp53^-/-) syngeneic mouse model of lung cancer after anti-PD-1 treatment. For both subsets, we demonstrate decreased TCR diversity in response to therapy. At a later time point, repertoire diversity is restored in progressing disease but remains decreased in responders to therapy in both CD4⁺ and CD8⁺ subsets. These observations complement previous studies and suggest that stably increased intratumoral CD4⁺ and CD8⁺ T cell clonality after anti-PD-1/PD-L1 therapy could serve as a predictor of long-term response.

Intratumoral immunoglobulin isotypes predict survival in lung adenocarcinoma subtypes

BACKGROUND

The role of tumor-infiltrating B-cells (TIBs) and intratumorally-produced antibodies in cancer-immunity interactions essentially remains terra incognita. In particular, it remains unexplored how driver mutations could be associated with distinct TIBs signatures and their role in tumor microenvironment.

METHODS

Here we analyzed associations of immunoglobulin isotypes and clonality with survival in TCGA RNA-Seq data for lung adenocarcinoma (LUAD), stratifying patients into 12 driver mutation and phenotypic tumor subgroups.

RESULTS

We revealed several unexpected associations between TIBs behavior and prognosis. Abundance and high proportion of IgG1 isotype, and low proportion of IgA among all intratumorally produced immunoglobulins were specifically associated with improved overall survival for KRAS^mut but not KRAS^wt LUAD, revealing the first link between a driver mutation and B-cell response. We found specific IgG1 signature associated with long survival, which suggests that particular specificities of IgG1+ TIBs could be beneficial in KRAS^mut LUAD. In contrast to our previous observations for melanoma, highly clonal IgG1 production by plasma cells had no meaningful effect on prognosis, suggesting that IgG1+ TIBs may exert a beneficial effect in KRAS^mut cases in an alternative way, such as efficient presentation of cognate antigens or direct B cell attack on tumor cells. Notably, a high proportion of the IgG1 isotype is positively correlated with the non-silent mutation burden both in the general LUAD cohort and in most patient subgroups, supporting a role for IgG1⁺ TIBs in antigen presentation. Complementing the recent finding that the presence of stromal IgG4-producing cells is associated with a favorable prognosis for patients with stage I squamous cell carcinoma, we show that the abundance of IgG4-producing TIBs likewise has a strong positive effect on overall survival in STK11^mut and proximal proliferative subgroups of LUAD patients. We hypothesize that the positive role of IgG4 antibodies in some of the lung cancer subtypes could be associated with reported inability of IgG4 isotype to form immune complexes, thus preventing immunosuppression via activation of the myeloid-derived suppressor cell (MDSC) phenotype.

CONCLUSIONS

We discover prominent and distinct associations between TIBs antibody isotypes and survival in lung adenocarcinoma carrying specific driver mutations. These findings indicate that particular types of tumor-immunity relations could be beneficial in particular driver mutation context, which should be taken into account in developing strategies of cancer immunotherapy and combination therapies. Specificity of protective B cell populations in specific cancer subgroups could become a clue to efficient targeted immunotherapies for appropriate cohorts of patients.

PiggyBac transposon tools for recessive screening identify B-cell lymphoma drivers in mice

B-cell lymphoma (BCL) is the most common hematologic malignancy. While sequencing studies gave insights into BCL genetics, identification of non-mutated cancer genes remains challenging. Here, we describe PiggyBac transposon tools and mouse models for recessive screening and show their application to study clonal B-cell lymphomagenesis. In a genome-wide screen, we discover BCL genes related to diverse molecular processes, including signaling, transcriptional regulation, chromatin regulation, or RNA metabolism. Cross-species analyses show the efficiency of the screen to pinpoint human cancer drivers altered by non-genetic mechanisms, including clinically relevant genes dysregulated epigenetically, transcriptionally, or post-transcriptionally in human BCL. We also describe a CRISPR/Cas9-based in vivo platform for BCL functional genomics, and validate discovered genes, such as Rfx7, a transcription factor, and Phip, a chromatin regulator, which suppress lymphomagenesis in mice. Our study gives comprehensive insights into the molecular landscapes of BCL and underlines the power of genome-scale screening to inform biology.

Comparative Analysis of B-Cell Receptor Repertoires Induced by Live Yellow Fever Vaccine in Young and Middle-Age Donors

Age-related changes can significantly alter the state of adaptive immune system and often lead to attenuated response to novel pathogens and vaccination. In present study we employed 5′RACE UMI-based full length and nearly error-free immunoglobulin profiling to compare plasma cell antibody repertoires in young (19–26 years) and middle-age (45–58 years) individuals vaccinated with a live yellow fever vaccine, modeling a newly encountered pathogen. Our analysis has revealed age-related differences in the responding antibody repertoire ranging from distinct IGH CDR3 repertoire properties to differences in somatic hypermutation intensity and efficiency and antibody lineage tree structure. Overall, our findings suggest that younger individuals respond with a more diverse antibody repertoire and employ a more efficient somatic hypermutation process than elder individuals in response to a newly encountered pathogen.

Reconstructing Antibody Repertoires from Error-Prone Immunosequencing Reads

Transforming error-prone immunosequencing datasets into Ab repertoires is a fundamental problem in immunogenomics, and a prerequisite for studies of immune responses. Although various repertoire reconstruction algorithms were released in the last 3 y, it remains unclear how to benchmark them and how to assess the accuracy of the reconstructed repertoires. We describe an accurate IgReC algorithm for constructing Ab repertoires from high-throughput immunosequencing datasets and a new framework for assessing the quality of reconstructed repertoires. Surprisingly, Ab repertoires constructed by IgReC from barcoded immunosequencing datasets in the blind mode (without using information about unique molecular identifiers) improved upon the repertoires constructed by the state-of-the-art tools that use barcoding. This finding suggests that IgReC may alleviate the need to generate repertoires using the barcoding technology (the workhorse of current immunogenomics efforts) because our computational approach to error correction of immunosequencing data is nearly as powerful as the experimental approach based on barcoding.

Application of nonsense-mediated primer exclusion (NOPE) for preparation of unique molecular barcoded libraries

Background

Recently we proposed efficient method to exclude undesirable primers at any stage of amplification reaction, here termed NOPE (NOnsense-mediated Primer Exclusion). According to this method, added oligonucleotide overlapping with the 3′-end of unwanted amplification primer (NOPE oligo) simultaneously provides a template for its elongation. This elongation disrupts specificity of unwanted primer, preventing its further participation in PCR. The suggested approach allows to rationally manage the course of PCR reactions in order to facilitate analysis of complex DNA mixtures as well as to perform multistage PCR bypassing intermediate purification steps.

Results

Here we apply NOPE method to DNA library preparation for the high-throughput sequencing (HTS) with the PCR-based introduction of unique molecular identifiers (UMI). We show that NOPE oligo efficiently neutralizes UMI-containing oligonucleotides after introduction of UMI into sample DNA molecules, thus allowing to proceed with further amplification steps without purification and associated loss of starting material. At the same time, NOPE oligo does not affect the efficiency of target PCR amplification.

Conclusion

We describe a simple, robust and cheap modification of UMI-labeled HTS libraries preparation procedure, that allows to bypass purification step and thus to preserve starting material which may be limited, e.g. circulating tumor DNA, circulating fetal DNA, or small amounts of isolated cells of interest. Furthermore, demonstrated simplicity and robustness of NOPE method should make it popular in various PCR protocols.

Dynamics of Individual T Cell Repertoires: From Cord Blood to Centenarians

The diversity, architecture, and dynamics of the TCR repertoire largely determine our ability to effectively withstand infections and malignancies with minimal mistargeting of immune responses. In this study, we have employed deep TCRβ repertoire sequencing with normalization based on unique molecular identifiers to explore the long-term dynamics of T cell immunity. We demonstrate remarkable stability of repertoire, where approximately half of all T cells in peripheral blood are represented by clones that persist and generally preserve their frequencies for 3 y. We further characterize the extremes of lifelong TCR repertoire evolution, analyzing samples ranging from umbilical cord blood to centenarian peripheral blood. We show that the fetal TCR repertoire, albeit structurally maintained within regulated borders due to the lower numbers of randomly added nucleotides, is not limited with respect to observed functional diversity. We reveal decreased efficiency of nonsense-mediated mRNA decay in umbilical cord blood, which may reflect specific regulatory mechanisms in development. Furthermore, we demonstrate that human TCR repertoires are functionally more similar at birth but diverge during life, and we track the lifelong behavior of CMV- and EBV-specific T cell clonotypes. Finally, we reveal gender differences in dynamics of TCR diversity constriction, which come to naught in the oldest age. Based on our data, we propose a more general explanation for the previous observations on the relationships between longevity and immunity.

VDJtools: Unifying Post-analysis of T Cell Receptor Repertoires

Despite the growing number of immune repertoire sequencing studies, the field still lacks software for analysis and comprehension of this high-dimensional data. Here we report VDJtools, a complementary software suite that solves a wide range of T cell receptor (TCR) repertoires post-analysis tasks, provides a detailed tabular output and publication-ready graphics, and is built on top of a flexible API. Using TCR datasets for a large cohort of unrelated healthy donors, twins, and multiple sclerosis patients we demonstrate that VDJtools greatly facilitates the analysis and leads to sound biological conclusions. VDJtools software and documentation are available at https://github.com/mikessh/vdjtools.

Age-related decrease in TCR repertoire diversity measured with deep and normalized sequence profiling

The decrease of TCR diversity with aging has never been studied by direct methods. In this study, we combined high-throughput Illumina sequencing with unique cDNA molecular identifier technology to achieve deep and precisely normalized profiling of TCR β repertoires in 39 healthy donors aged 6-90 y. We demonstrate that TCR β diversity per 10(6) T cells decreases roughly linearly with age, with significant reduction already apparent by age 40. The percentage of naive T cells showed a strong correlation with measured TCR diversity and decreased linearly up to age 70. Remarkably, the oldest group (average age 82 y) was characterized by a higher percentage of naive CD4(+) T cells, lower abundance of expanded clones, and increased TCR diversity compared with the previous age group (average age 62 y), suggesting the influence of age selection and association of these three related parameters with longevity. Interestingly, cross-analysis of individual TCR β repertoires revealed a set >10,000 of the most representative public TCR β clonotypes, whose abundance among the top 100,000 clones correlated with TCR diversity and decreased with aging.

Mother and child T cell receptor repertoires: deep profiling study

The relationship between maternal and child immunity has been actively studied in the context of complications during pregnancy, autoimmune diseases, and haploidentical transplantation of hematopoietic stem cells and solid organs. Here, we have for the first time used high-throughput Illumina HiSeq sequencing to perform deep quantitative profiling of T cell receptor (TCR) repertoires for peripheral blood samples of three mothers and their six children. Advanced technology allowed accurate identification of 5 × 10⁵ to 2 × 10⁶ TCR beta clonotypes per individual. We performed comparative analysis of these TCR repertoires with the aim of revealing characteristic features that distinguish related mother-child pairs, such as relative TCR beta variable segment usage frequency and relative overlap of TCR beta complementarity-determining region 3 (CDR3) repertoires. We show that thymic selection essentially and similarly shapes the initial output of the TCR recombination machinery in both related and unrelated pairs, with minor effect from inherited differences. The achieved depth of TCR profiling also allowed us to test the hypothesis that mature T cells transferred across the placenta during pregnancy can expand and persist as functional microchimeric clones in their new host, using characteristic TCR beta CDR3 variants as clonal identifiers.

Preparing unbiased T-cell receptor and antibody cDNA libraries for the deep next generation sequencing profiling

High-throughput sequencing has the power to reveal the nature of adaptive immunity as represented by the full complexity of T-cell receptor (TCR) and antibody (IG) repertoires, but is at present severely compromised by the quantitative bias, bottlenecks, and accumulated errors that inevitably occur in the course of library preparation and sequencing. Here we report an optimized protocol for the unbiased preparation of TCR and IG cDNA libraries for high-throughput sequencing, starting from thousands or millions of live cells in an investigated sample. Critical points to control are revealed, along with tips that allow researchers to minimize quantitative bias, accumulated errors, and cross-sample contamination at each stage, and to enhance the subsequent bioinformatic analysis. The protocol is simple, reliable, and can be performed in 1–2 days.

[Characterization of circulating RNA in plasma as potential tool for breast cancer diagnostics]

The representation patterns of 15 cytokines RNA in blood plasma and blood cells of patients with breast cancer and apparently healthy women were investigated. Relative levels of RNA IL-8 and IL-18 in plasma of breast cancer patients are significantly increased compared with control group. At the same time no obvious differences were found in relative concentrations of these transcripts in blood cells of patients and control groups. Relative concentration of IL-8 RNA was higher in blood plasma of locally advanced compared with early breast cancer patients.

Quantitative tracking of T cell clones after haematopoietic stem cell transplantation

Autologous haematopoietic stem cell transplantation is highly efficient for the treatment of systemic autoimmune diseases, but its consequences for the immune system remain poorly understood. Here, we describe an optimized RNA-based technology for unbiased amplification of T cell receptor beta-chain libraries and use it to perform the first detailed, quantitative tracking of T cell clones during 10 months after transplantation. We show that multiple clones survive the procedure, contribute to the immune response to activated infections, and form a new skewed and stable T cell receptor repertoire.

[Profile of RNA cytokines in blood plasma under conditions of normal physiological state of human body]

The level of representation of extracellular RNA 14 cytokines in blood plasma in a group of apparently healthy subjects was analyzed. The level of representation of the transcripts of these cytokines in extracellular medium is characterized by specific profile different from the profile of expression of the genes in blood cells.