Intratumoral immunoglobulin isotypes predict survival in lung adenocarcinoma subtypes
Journal for ImmunoTherapy of Cancer volume 7, Article number: 279 (2019)
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.
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.
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 KRASmut but not KRASwt 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 KRASmut 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 KRASmut 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 STK11mut 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.
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.
Recent work has revealed the importance of the antigenic specificity [1,2,3], clonality [4, 5], hypermutation [5, 6], and isotype [5,6,7,8] of TIBs, circulating plasmablasts , and serum self- and tumor-specific antibodies [10, 11] in tumor-immunity interactions. In particular, TIBs and tumor-infiltrating plasma cells—or more accurately, intra-tumoral B cells, since these may also be generated in local tertiary lymphoid structures (TLS) [1, 12,13,14,15]—are believed to play important roles in the tumor microenvironment. Their major modes of action include presentation of B-cell receptor (BCR)-cognate antigens to CD4+ and potentially CD8+ T-cells [15,16,17,18], production of cytokines that can stimulate or suppress anti-tumor response [19,20,21], and production of tumor-specific antibodies [1,2,3] that may enhance killing of tumor cells via ADCC [8, 22], enhance antigen capture and presentation by dendritic cells , or form immune complexes that promote activation of MDSCs [23, 24].
The involvement of specific BCRs and antibodies in both antitumor—and potentially, pro-tumor—reactions is becoming increasingly clear. Serum auto-antibodies against tumor-associated and self-antigens have been validated as biomarkers for the early detection of cancer [11, 25, 26] and could also serve as useful prognostic markers at later stages of disease [27, 28]. Blood plasmablast counts and BCR clonality were recently shown to be hallmarks of non-progressing cancer in patients treated with anti-CTLA4 immunotherapy, and the therapy itself led to an increase in these parameters . The amount of B-cells and plasma cells in a tumor is also associated with a good prognosis for diverse solid tumors [29,30,31,32]. On the other hand, several clinical studies of hepatocellular carcinoma , prostate , renal , and breast  cancer have indicated that high B-cell or plasma cell content may be associated with a negative prognosis. This highlights the possibility that B cells can also help maintain an immunosuppressive microenvironment, a function that is sometimes attributed to a certain population of B-cells referred as regulatory B-cells [36, 37].
Antibody functionality is strongly influenced by isotype , and this feature is likely to be an important piece of the puzzle in terms of understanding B-cell-tumor interactions. Particular isotypes could also be associated with specific B-cell functionalities. For example, IgA+ plasma cells have been described as preferentially producing the immunosuppressive IL-10 and PD-L1 in some cancers [21, 33].
Based on our analysis of RNA-Seq data from TCGA, we recently demonstrated that high expression levels, proportion, and clonality (i.e., focused expression of particular clonal variants) of cytotoxic IgG1 antibodies is associated with a markedly better prognosis in melanoma, while a high proportion of the IgA isotype is associated with shorter survival .
Here, we have investigated the role of different antibody isotypes and clonality in lung adenocarcinoma (LUAD), splitting the 442 patients available in the TCGA database into relevant subgroups based on the presence of key driver mutations or transcriptional subtypes. We reveal previously unexpected associations between dominating isotypes of TIBs and survival in specific subgroups of LUAD patients. Our results for the first time link driver mutations and B cell response in tumor microenvironment, and suggest that patient stratification for immunotherapies and design of combination therapies should take both these parameters in consideration.
Initial data filtering
Patient data from the TCGA LUAD project was obtained from the GDC portal repository (https://portal.gdc.cancer.gov/). We downloaded the HTSeq-FPKM files and transformed the transcript-level data into gene-level data by summarizing alternative transcripts. FPKM were then transformed to TPM. Samples from formalin-fixed paraffin-embedded (FFPE) tissue, normal tissue or metastatic lesions as well as entities with warnings were removed. One sample for each patient was selected in accordance with GDC recommendations (https://confluence.broadinstitute.org/display/GDAC/FAQ#FAQ-replicateFilteringQ%C2%A0Whatdoyoudowhenmultiplealiquotbarcodesexistforagivensampleportionanalytecombination) After such filtration, a cohort of 442 patients was formed (“general cohort” in our study).
Mutation information was also obtained from GDC portal. Mutations with low Variant Effect Predictor (VEP) impact that were not annotated by SIFT or PolyPhen as having impact were excluded from analysis. We identified the following numbers of patients with relevant genotypes: KRASmut, 122 patients; KRASwt, 320 patients; STK11mut, 73 patients; STK11wt, 369 patients; TP53mut, 220 patients; TP53wt, 222 patients; EGFRmut, 57 patients.
PD-L1 levels were characterized based on CD274 gene expression data. Samples that had CD274 expression higher than twice the mean value in the general cohort were assigned to the PD-L1high group (N = 51); the remaining samples formed the PD-L1low group (N = 391).
IGH, IGA and IGHG
Total BCR/antibody expression (IGH) was calculated as a sum of the expressions of IGHA1, IGHA2, IGHG1, IGHG2, IGHG3, IGHG4, IGHM, IGHD and IGHE genes. IgA expression was calculated as a sum of expression values for the IGHA1 and IGHA2 genes.
Principal component analysis (prcomp function in R) was used to determine IGHV-IGHJ profiles that explain the most variance across the LUAD cohort. This analysis was limited to IgG1 CDR3 clonotypes, samples with less than 500 IgG1 CDR3-covering reads were removed. High and low survival was determined by comparing to median survival in corresponding cohort. Note that we limited to IGHV families, so e.g. IGHV3–11 and IGHV3–13 were treated as the same family IGHV3.
Information about the expression-based classification of the samples in the general cohort was obtained from Ref. . Only 184 patients in the general cohort had their expression-based subtype annotated, with 65 of them belonging to the proximal inflammatory (PI) subtype, 51 - to the proximal proliferative (PP) subtype, and 68 - to the terminal respiratory unit (TRU) subtype.
IgG1 and IgA clonality
To obtain clonality data, we have downloaded the BAM files with reads aligned by STAR from the GDC portal, using the Genomic Data Commons Bioconductor R package (https://bioconductor.org/packages/release/bioc/html/GenomicDataCommons.html). BAM files were then sorted with samtools  and converted to Fastq files using the SamToFastq Picard tool (http://broadinstitute.github.io/picard/). MiXCR software  was used to extract CDR3 repertoires from Fastq files, and VDJtools  was used for the repertoire statistical analysis. Only samples that had more than 500 IgG1 or IgA CDR3-covering sequencing reads were included in the analysis. IgG1 and IgA CDR3 repertoires were downsampled to 500 randomly chosen reads for normalization purposes. Clonality was calculated as: 1 - normalized Shannon-Wiener index .
Survival plots were created using the Kaplan-Meier estimator. Plots were created using matplotlib  based on modified functions from the lifelines package (https://zenodo.org/record/2638135#.XMCtiegzaUl). We used a statistical significance threshold of p = 0.05. Data analysis was performed with Python2 and R. Multivariable analysis was performed with Cox proportional hazard regression.
Non-silent mutation burden
The non-silent mutation burden per megabase for each sample was obtained from Ref. . The correlation between IGHG1/IGH proportion and non-silent mutation burden was calculated using Spearman rank correlation coefficient and visualized with Seaborn.
Extrapolating from our previous results obtained with TCGA data for human melanoma , we expected to observe an association between a high proportion of IGHG1/IGH and long survival, where IGH is a sum of the expression of the IGHA1, IGHA2, IGHG1, IGHG2, IGHG3, IGHG4, IGHM, IGHD and IGHE genes. However, this was not the case for the TCGA LUAD cohort as a whole (Fig. 1a, hereinafter patient cohorts are split by median).
A broader analysis of each of the cancer types available in TCGA revealed that a high IGHG1 proportion is only associated with a significantly better prognosis for the full patient cohorts in non-papillary bladder cancer (Bonferroni adjusted p < 0.002, HR = 0.5) and melanoma (adjusted p < 0.02, HR = 0.6).
However, each cancer type is heterogeneous. It is expected that donors having the same cancer type can show huge variance in survival and immune response owing to the heterogeneity of mutation and gene expression profiles defining cancer phenotype. In an effort to reveal the distinct B-cell characteristics associated with specific types of LUAD and control for potential phenotypic differences, we split patients into 12 partially-overlapping genetic and phenotypic subgroups. These included 3 oncogenic driver status subgroups (KRASmut, KRASwt, EGFRmut), 4 tumor suppressor genes subgroups (p53mut, p53wt, STK11mut, STK11wt), PD-L1high and PD-L1low subgroups, and proximal inflammatory, proximal proliferative, and terminal respiratory unit transcriptional subtypes .
Remarkably, out of the 12 subgroups we investigated here, only KRASmut subgroup demonstrated a significant association of a high IGHG1/IGH proportion with overall survival (Fig. 1a, adjusted p = 0.01, HR = 0.36). A multivariable analysis using Cox proportional-hazards regression with adjustment for stage, gender, smoking, age, and infiltration (CD45 expression) confirmed that a high IGHG1/IGH proportion is associated with overall survival for the KRASmut subgroup (p = 0.01, HR = 0.38). In contrast, we observed an inverse - albeit not statistically significant - association in patients with KRASwt tumors (Fig. 1a).
In the IGHG1/IGHlow subgroup, patients with KRASmut tumor status had a worse prognosis compared to KRASwt (adjusted p = 0.006), while in the IGHG1/IGHhigh subgroup, there was no significant difference (Fig. 1b). These results suggest that a high IGHG1/IGH proportion plays a protective role in the KRASmut but not in the KRASwt context.
At the same time, in contrast to melanoma, high IGHG1/MS4A1 ratio, reflecting the relative abundance of IgG1-producing plasma cells compared to CD20+ (i.e., non-plasma) B-cells, is not associated with longer survival in KRASmut tumors, and tends to have a negative association for KRASwt tumors (data not shown). This observation argues against IgG1-mediated ADCC playing a major protective role in KRASmut LUAD.
Remarkably, the IGHG1/IGH proportion is positively correlated with the non-silent mutation burden both for the general LUAD cohort (R = 0.48, p = 6 × 10− 12, Fig. 1c) and for most LUAD subgroups (Additional file 1: Figure S1). This finding may support the point that IgG1 B cells are involved in the antigen presentation process.
B-cells versus plasma cells
Abundance of TIBs measured based on CD19 expression level was associated with a positive prognosis in general LUAD cohort (adjusted p = 0.03) and in most subgroups, in agreement with previous works based on immunohistochemical analysis [1, 46], tissue microarrays [47, 48], RNA-Seq [31, 32], and RNA expression microarrays [30, 49].
However, all of these prior studies have considered LUAD as a general cohort, while here our goal is to find distinct dependencies in LUAD subgroups that can potentially be characterized by different types of the balance in tumor-immunity interactions. Analysis of the 12 subgroups described above has revealed that B-cell infiltration as measured by CD19 expression level has an especially beneficial impact on survival for the proximal proliferative LUAD transcriptional subtype (Fig. 2a).
Although high IGH expression did not have any meaningful association with survival for most subgroups, it was significantly beneficial for proximal proliferative LUAD (adjusted p = 0.02, Fig. 2b). A multivariable analysis using Cox proportional-hazards regression with adjustment for stage, gender, smoking, age, and infiltration confirmed that high IGH expression is associated with prolonged overall survival in proximal proliferative LUAD (p = 0.006, HR = 0.08).
High SDC1 (encoding CD138, indicator of plasma cells) expression was associated with a neutral or non-significant negative effect on overall survival in all cohorts with the exception of proximal proliferative LUAD, where non-significant association with longer survival was observed (data not shown).
We also assessed the ratio of IGH to MS4A1 (encoding CD20) expression, which reflects the relative abundance of CD20-negative antibody-producing plasma cells compared to CD20-positive non-plasma B-cells. This ratio generally had a slightly negative or neutral effect in all groups, but we observed a non-significant association with positive prognosis in the proximal proliferative LUAD subgroup (Fig. 2c).
We concluded that, although tumor infiltration with CD19-positive B cells is generally a positive signature for most types of LUAD, the presence of antibody-producing plasma cells is specifically associated with better tumor immunosurveillance in the proximal proliferative LUAD subgroup.
IgG1 antibody production and clonality
We  and others [4, 50] have earlier identified an association of high IGH (mainly IgG1 ) “clonality” with better survival in melanoma patients, where this metric is calculated as [1 – the normalized Shannon-Wiener index] . In the T-cell world, this metric reflects the relative presence of large clonal expansions. For B-cells, this also reflects RNA expression levels that may differ dramatically between B-cells of differing functional status—with average expression varying by factors of as much as 2:5:500 for naïve, memory and plasma cells, respectively, according to our recent estimations . For highly-infiltrated tumors, antibody CDR3 regions are covered by a relatively large proportion of RNA-Seq reads, which makes it possible to efficiently extract intratumorally-produced immunoglobulin repertoires with MiXCR  and thereby analyze clonality metrics. We extracted IgG1 CDR3 repertoires from all patient samples from general LUAD cohort, but only used data from 283 out of 442 patients with > 500 sequencing reads covering IgG1 CDR3, as this is the minimum coverage that allows us to accurately assess clonality . Therefore, this analysis was performed only for tumor samples with relatively high IgG1 expression. For normalization, each dataset was down-sampled to 500 randomly-chosen CDR3-covering sequencing reads.
Notably, high IgG1 clonality, which reflects the presence of a focused IgG1 plasma cell response, did not influence the prognosis for KRASmut LUAD patients (Additional file 1: Figure S2). The neutral effect of both IGHG1/MS4A1 ratio and IgG1 clonality in KRASmut tumors suggests that IgG1-producing plasma B-cells do not play a prominent role as key drivers of anti-tumor response via ADCC in this subtype of LUAD. In contrast, there is evidence for such a model in melanoma, based on the correlation of large hypermutating clonal IgG1 expansions  and high IGHG1/MS4A1 ratio with survival (p = 0.006, HR = 0.7), and cytotoxic activity of tumor-specific IgG1 antibodies .
One possible interpretation is that in KRASmut tumors, the abundance and high proportion of IgG1 TIBs may play an active positive role via presentation of cognate antigens. Among other TAA and neoantigens, the IgG1 shift of TIBs could lead to more efficient presentation of the mutated KRAS peptide itself. Recent work by the Rosenberg group revealed KRASmut-specific CD4+ T-cells , and KRAS-specific tumor-infiltrating IgG B-cells were identified in patients with pancreatic cancer . Ability of lung tumor-infiltrating B-cells to present antigens and activate CD4+ T-cells was also reported .
Given the fact that a high IGHG1//IGH proportion is associated with longer survival in KRASmut, but not KRASwt LUAD cases, we have explicitly tested the repertoires for the presence of specific IgG1 sequence motifs that can be associated with survival. Analysis of IGHV-IGHJ profiles that are most variable across LUAD samples (see Methods section) has revealed the presence of a specific signature that is up-regulated in KRASmut cases with high survival (Additional file 1: Figure S3), which is characterized by high abundance of IGHV6-IGHJ4 and IGHV4-IGHJ3 clonotypes and low abundance of IGHV3-IGHJ1 and IGHV3-IGHJ2 clonotypes. These results may indicate that the response to particular tumor antigens is associated with tumor immunosurveillance in KRASmut LUAD, a hypothesis that will require further investigation to confirm.
IgA and patient survival
High IgA expression levels (the sum of IGHA1 and IGHA2 genes) was a neutral parameter in all LUAD subgroups, including KRASmut (data not shown). However, high IgA/IGH (reflecting the proportion of IgA among all intratumorally-produced antibodies) and IgA/MS4A1 (reflecting IgA production relative to non-plasma B cell abundance) proportions were associated with shorter survival in KRASmut but not in KRASwt patients (Additional file 1: Figure S4a). This dominant presence of IgA could be an indirect consequence of a deficiency in the relative proportion of IgG1 (Fig. 1) and IgG4 (see below), thus leading to a negative prognosis in the KRASmut subgroup. At the same time, the effect of the IGHG1/MS4A1 parameter was neutral, while increased IgA/MS4A1 had a negative effect on survival in the KRASmut subgroup (Additional file 1: Figure S4b). This observation supports the negative role of IgA-producing plasma cells in KRASmut LUAD, as was previously reported for hepatocellular carcinoma  and bladder cancer . Notably, the IgA clonality parameter remained neutral in all 12 analyzed subgroups (data not shown), suggesting that antigenic specificity of the antibodies produced by the IgA plasma cells does not play a critical role in survival. At the same time, IgA-positive B cells could be involved in antigen presentation, skewing CD4+ T cells towards functional phenotypes that are negative or suboptimal for anti-tumor response .
IgG4 expression and patient survival
In 2013, Fujimoto and coauthors used immunohistochemistry analysis to show that the presence of stromal B-cells producing high levels of IgG4 is associated with prominently better prognosis in patients with stage I squamous cell carcinoma . Our analysis extends this observation to LUAD, and delineates a group of LUAD patients that benefit in particular from the presence of IgG4-producing B cells.
High intra-tumoral IGHG4 expression levels were associated with a better prognosis for the general LUAD cohort (adjusted p = 0.06, HR = 0.64). This positive association was based on the positive effect of elevated IGHG4 levels in TP53wt (adjusted p = 0.04, HR = 0.49), PD-L1low (HR = 0.64), STK11mut (HR = 0.4), and proximal proliferative (HR = 0.26) LUAD (Fig. 3a), whereas no such beneficial effect was observed in other patient subgroups (not shown).
For proximal proliferative LUAD patients, we also found that the overall abundance of several other antibody isotypes, including IgM, IgG1, IgG2, IgG3 - and unexpectedly, IgE - was also associated with better prognosis (Additional file 1: Fig. S5a). This is in agreement with the beneficial role of CD19+ B-cells and high IGH expression in this subgroup (Fig. 2a,b). But in terms of relative proportion among all IGH, only high representation of IgG4 – and to some extent IgM, but not other isotypes – tended to be associated with favorable prognosis (Fig. 3b, Additional file 1: Figure S5b).
A high IGHG4/IGH proportion was also beneficial for the STK11mut subgroup (Fig. 3a,b). A multivariable analysis using Cox proportional-hazards regression with adjustment for stage, gender, smoking, age, and infiltration showed that a high IGHG4/IGH ratio was significantly associated with overall survival in STK11mut LUAD (p = 0.04, HR = 0.4).
Notably, proximal proliferative transcriptional portrait of LUAD tumors is associated with alteration of STK11 tumor suppressor gene , while STK11 alteration and absence of TP53 mutations are correlated with PD-L1low LUAD subtype [56, 57]. The patterns that we have observed here thus could represent components of the same picture, describing the cumulative tumor portrait of lung adenocarcinoma patients that mostly benefit from abundant intratumoral IgG4 expression as TP53wt/STK11mut/PD-L1low/proximal proliferative LUAD.
Indeed, for the small subgroup of 19 patients with both STK11mut and proximal proliferative LUAD, the positive association of survival with high IGHG4 expression (adjusted p = 0.05) and proportion (adjusted p = 0.04) was even more prominent (Fig. 3).
When we normalized to CD20—measuring relative IgG4 production compared to the abundance of non-plasma B-cells—we determined that a high IGHG4/MS4A1 ratio only tended to be associated with a positive prognosis in proximal proliferative LUAD (Fig. 3c).
Clearly, the impact of B-cells in cancer immunology is not black and white, and we cannot simply distinguish “positive” IgG1/IgG3 isotypes that initiate tumor-specific ADCC and immune responses from “negative” IgA/IgG4 isotypes as a signature of or precursor to immunosuppression. The particular antigenic specificities of intratumoral BCRs/antibodies - which include the surface versus intracellular localization of cognate tumor antigens - and associated phenotypes of antigen-presenting and cytokine-producing B-cells all contribute to the complex picture of tumor-immunity interactions.
Nevertheless, here we have delineated LUAD patient subgroups that can be characterized by striking dependencies between the abundance and proportion of particular intratumorally-produced BCR/antibody isotypes and survival. Association of IgG1 and IgG4 isotype dominance with favorable prognosis in KRASmut and STK11mut/proximal proliferative LUAD patients, respectively, hints at the existence of specific types of established tumor-immunity interaction profiles. The latter could either involve driver mutations themselves in antigen-specific response [52, 53, 58] or result from downstream pathways characteristic for specific driver mutation.
This discovery of existing links between driver mutations and TIB-mediated immunity complements recently described interconnections between driver mutations, T-cell behavior, and PD-L1 expression [56, 59]. In particular, the KRASG12D mutation and MEK/ERK pathway activation was shown to up-regulate production of IL-10 and TGF-β, thus promoting CD4 T cells conversion in Tregs in pancreatic cancer . In LUAD, KRAS mutation was associated with more intense immune cells infiltration . Several other driver mutations correlated with lower or higher leukocyte infiltration across all cancer types .
The positive influence of IgG1 TIBs in KRASmut tumors could be explained by presentation of BCR-cognate tumor antigens to CD4+ T-cells. Considered alongside recent reports revealing the importance of antigen-specific B-cells as cognate antigen presenters [6, 7, 18, 62, 63], these results support the concept that therapeutic vaccination using whole proteins or their encoding genes (including KRASmut)  could more efficiently exploit the antigen-presentation machinery of cognate B-cells. The hypothesis that mutant KRAS peptide itself is among the involved antigenic targets is especially attractive since, in contrast to other neoantigens, the driver mutation is a sensitive component of tumor survival. However, exploring this hypothesis in depth will require further investigation.
The reason for the observed association of high IgG4 production with a favorable prognosis in STK11mut and proximal proliferative LUAD remains unclear, and will require further progress in our fundamental understanding of the functionality of the IgG4 isotype.
We hypothesize the following explanation. In IgG4, inter-heavy chain disulfides are in equilibrium with intra-heavy chain disulfides , which enables heavy chain monomer exchange in vivo . As a result, IgG4 functions as a monovalent antibody, which is unable to cross-link antigen and form immune complexes . Notably, persistent immune complexes formed by tumor-specific antibodies may be associated with an unfavorable clinical outcome  due to their immunosuppressive action through the modulation of FcR-bearing myeloid cell activity, leading to a MDSC phenotype [23, 24]. Thus, a positive role of IgG4 in lung cancer could be connected with diminished formation of immune complexes and subsequent MDSC-related immunosuppressive reactions. At the same time, it should be noted that antibodies of IgG4 isotype may also have a negative impact on prognosis for some cancer types, as has been reported for human melanoma .
Our discovery of direct links existing between antibody isotypes and survival in lung adenocarcinoma carrying specific driver mutations strengthens the importance of TIBs as immune system players with multi-parametric roles in the battle with cancer. This may suggest prospective strategies for more rational design of combination approaches incorporating targeted therapies, immune checkpoint inhibitors, and vaccines. In particular, these results indicate that immunotherapy efforts must take into account B-cell component of the tumor microenvironment, which role, most importantly, may depend on particular context of driver mutations.
Availability of data and materials
Extracted IGH, IgA, and IgG1 CDR3 repertoires, metadata, expression and clonality metrics are deposited on Figshare: https://figshare.com/projects/BCR_profiling_in_lung_adenocarcinoma_TCGA_cohort/64106.
Antibody-Dependent Cellular Cytotoxicity
Formalin-Fixed Paraffin-Embedded tissue
Fragments Per Kilobase Million
Myeloid-Derived Suppressor Cell
The Cancer Genome Atlas
Germain C, Gnjatic S, Tamzalit F, Knockaert S, Remark R, Goc J, Lepelley A, Becht E, Katsahian S, Bizouard G, et al. Presence of B cells in tertiary lymphoid structures is associated with a protective immunity in patients with lung cancer. Am J Respir Crit Care Med. 2014;189(7):832–44.
Carmi Y, Spitzer MH, Linde IL, Burt BM, Prestwood TR, Perlman N, Davidson MG, Kenkel JA, Segal E, Pusapati GV, et al. Allogeneic IgG combined with dendritic cell stimuli induce antitumour T-cell immunity. Nature. 2015;521(7550):99–104.
Yasuda M, Mizukami M, Hanagiri T, Shigematsu Y, Fukuyama T, Nagata Y, So T, Ichiki Y, Sugaya M, Takenoyama M, et al. Antigens recognized by IgG derived from tumor-infiltrating B lymphocytes in human lung cancer. Anticancer Res. 2006;26(5A):3607–11.
Mose LE, Selitsky SR, Bixby LM, Marron DL, Iglesia MD, Serody JS, Perou CM, Vincent BG, Parker JS. Assembly-based inference of B-cell receptor repertoires from short read RNA sequencing data with V'DJer. Bioinformatics. 2016;32(24):3729–34.
Bolotin DAea: Antigen receptor repertoire profiling from RNA-seq data. Nat Biotechnol 2017.
Nielsen JS, Sahota RA, Milne K, Kost SE, Nesslinger NJ, Watson PH, Nelson BH. CD20+ tumor-infiltrating lymphocytes have an atypical CD27- memory phenotype and together with CD8+ T cells promote favorable prognosis in ovarian cancer. Clin Cancer Res. 2012;18(12):3281–92.
Shi JY, Gao Q, Wang ZC, Zhou J, Wang XY, Min ZH, Shi YH, Shi GM, Ding ZB, Ke AW, et al. Margin-infiltrating CD20(+) B cells display an atypical memory phenotype and correlate with favorable prognosis in hepatocellular carcinoma. Clin Cancer Res. 2013;19(21):5994–6005.
Gilbert AE, Karagiannis P, Dodev T, Koers A, Lacy K, Josephs DH, Takhar P, Geh JL, Healy C, Harries M, et al. Monitoring the systemic human memory B cell compartment of melanoma patients for anti-tumor IgG antibodies. PLoS One. 2011;6(4):e19330.
DeFalco J, Harbell M, Manning-Bog A, Baia G, Scholz A, Millare B, Sumi M, Zhang D, Chu F, Dowd C, et al. Non-progressing cancer patients have persistent B cell responses expressing shared antibody paratopes that target public tumor antigens. Clin Immunol. 2018;187:37–45.
Kobold S, Lutkens T, Cao Y, Bokemeyer C, Atanackovic D. Autoantibodies against tumor-related antigens: incidence and biologic significance. Hum Immunol. 2010;71(7):643–51.
Macdonald IK, Parsy-Kowalska CB, Chapman CJ. Autoantibodies: opportunities for early Cancer detection. Trends in cancer. 2017;3(3):198–213.
Pitzalis C, Jones GW, Bombardieri M, Jones SA. Ectopic lymphoid-like structures in infection, cancer and autoimmunity. Nat Rev Immunol. 2014;14(7):447–62.
Germain C, Gnjatic S, Dieu-Nosjean MC. Tertiary lymphoid structure-associated B cells are key players in anti-tumor immunity. Front Immunol. 2015;6:67.
Colbeck EJ, Ager A, Gallimore A, Jones GW. Tertiary lymphoid structures in Cancer: drivers of antitumor immunity, immunosuppression, or bystander sentinels in disease? Front Immunol. 2017;8:1830.
Zhu W, Germain C, Liu Z, Sebastian Y, Devi P, Knockaert S, Brohawn P, Lehmann K, Damotte D, Validire P, et al. A high density of tertiary lymphoid structure B cells in lung tumors is associated with increased CD4(+) T cell receptor repertoire clonality. Oncoimmunology. 2015;4(12):e1051922.
Rivera A, Chen CC, Ron N, Dougherty JP, Ron Y. Role of B cells as antigen-presenting cells in vivo revisited: antigen-specific B cells are essential for T cell expansion in lymph nodes and for systemic T cell responses to low antigen concentrations. Int Immunol. 2001;13(12):1583–93.
Bouaziz JD, Yanaba K, Venturi GM, Wang Y, Tisch RM, Poe JC, Tedder TF. Therapeutic B cell depletion impairs adaptive and autoreactive CD4+ T cell activation in mice. Proc Natl Acad Sci U S A. 2007;104(52):20878–83.
Bruno TC, Ebner PJ, Moore BL, Squalls OG, Waugh KA, Eruslanov EB, Singhal S, Mitchell JD, Franklin WA, Merrick DT, et al. Antigen-presenting Intratumoral B cells affect CD4(+) TIL phenotypes in non-small cell lung Cancer patients. Cancer immunology research. 2017;5(10):898–907.
Dang VD, Hilgenberg E, Ries S, Shen P, Fillatreau S. From the regulatory functions of B cells to the identification of cytokine-producing plasma cell subsets. Curr Opin Immunol. 2014;28:77–83.
Zhou J, Min Z, Zhang D, Wang W, Marincola F, Wang X. Enhanced frequency and potential mechanism of B regulatory cells in patients with lung cancer. J Transl Med. 2014;12:304.
Shalapour S, Font-Burgada J, Di Caro G, Zhong Z, Sanchez-Lopez E, Dhar D, Willimsky G, Ammirante M, Strasner A, Hansel DE, et al. Immunosuppressive plasma cells impede T-cell-dependent immunogenic chemotherapy. Nature. 2015;521(7550):94–8.
Kurai J, Chikumi H, Hashimoto K, Yamaguchi K, Yamasaki A, Sako T, Touge H, Makino H, Takata M, Miyata M, et al. Antibody-dependent cellular cytotoxicity mediated by cetuximab against lung cancer cell lines. Clin Cancer Res. 2007;13(5):1552–61.
Barbera-Guillem E, May KF Jr, Nyhus JK, Nelson MB. Promotion of tumor invasion by cooperation of granulocytes and macrophages activated by anti-tumor antibodies. Neoplasia. 1999;1(5):453–60.
Andreu P, Johansson M, Affara NI, Pucci F, Tan T, Junankar S, Korets L, Lam J, Tawfik D, DeNardo DG, et al. FcRgamma activation regulates inflammation-associated squamous carcinogenesis. Cancer Cell. 2010;17(2):121–34.
Dai L, Tsay JC, Li J, Yie TA, Munger JS, Pass H, Rom WN, Zhang Y, Tan EM, Zhang JY. Autoantibodies against tumor-associated antigens in the early detection of lung cancer. Lung Cancer. 2016;99:172–9.
Tang ZM, Ling ZG, Wang CM, Wu YB, Kong JL. Serum tumor-associated autoantibodies as diagnostic biomarkers for lung cancer: a systematic review and meta-analysis. PLoS One. 2017;12(7):e0182117.
Gnjatic S, Ritter E, Buchler MW, Giese NA, Brors B, Frei C, Murray A, Halama N, Zornig I, Chen YT, et al. Seromic profiling of ovarian and pancreatic cancer. Proc Natl Acad Sci U S A. 2010;107(11):5088–93.
Hirasawa Y, Kohno N, Yokoyama A, Kondo K, Hiwada K, Miyake M. Natural autoantibody to MUC1 is a prognostic indicator for non-small cell lung cancer. Am J Respir Crit Care Med. 2000;161(2 Pt 1):589–94.
Wouters MCA, Nelson BH. Prognostic significance of tumor-infiltrating B cells and plasma cells in human Cancer. Clin Cancer Res. 2018;24(24):6125–35.
Gentles AJ, Newman AM, Liu CL, Bratman SV, Feng W, Kim D, Nair VS, Xu Y, Khuong A, Hoang CD, et al. The prognostic landscape of genes and infiltrating immune cells across human cancers. Nat Med. 2015;21(8):938–45.
Iglesia MD, Parker JS, Hoadley KA, Serody JS, Perou CM, Vincent BG. Genomic Analysis of Immune Cell Infiltrates Across 11 Tumor Types. J Natl Cancer Inst. 2016:108(11).
Charoentong P, Finotello F, Angelova M, Mayer C, Efremova M, Rieder D, Hackl H, Trajanoski Z. Pan-cancer Immunogenomic analyses reveal genotype-Immunophenotype relationships and predictors of response to checkpoint blockade. Cell Rep. 2017;18(1):248–62.
Shalapour S, Lin XJ, Bastian IN, Brain J, Burt AD, Aksenov AA, Vrbanac AF, Li W, Perkins A, Matsutani T, et al. Inflammation-induced IgA+ cells dismantle anti-liver cancer immunity. Nature. 2017;551(7680):340–5.
Sjoberg E, Frodin M, Lovrot J, Mezheyeuski A, Johansson M, Harmenberg U, Egevad L, Sandstrom P, Ostman A. A minority-group of renal cell cancer patients with high infiltration of CD20+B-cells is associated with poor prognosis. Br J Cancer. 2018;119(7):840–6.
Mohammed ZM, Going JJ, Edwards J, Elsberger B, McMillan DC. The relationship between lymphocyte subsets and clinico-pathological determinants of survival in patients with primary operable invasive ductal breast cancer. Br J Cancer. 2013;109(6):1676–84.
Balkwill F, Montfort A, Capasso M. B regulatory cells in cancer. Trends Immunol. 2013;34(4):169–73.
Mauri C, Menon M. Human regulatory B cells in health and disease: therapeutic potential. J Clin Invest. 2017;127(3):772–9.
Schroeder HW Jr, Cavacini L. Structure and function of immunoglobulins. J Allergy Clin Immunol. 2010;125(2 Suppl 2):S41–52.
Bolotin DA, Poslavsky S, Davydov AN, Frenkel FE, Fanchi L, Zolotareva OI, Hemmers S, Putintseva EV, Obraztsova AS, Shugay M, et al. Antigen receptor repertoire profiling from RNA-seq data. Nat Biotechnol. 2017;35(10):908–11.
Cancer Genome Atlas Research N. Comprehensive molecular profiling of lung adenocarcinoma. Nature. 2014;511(7511):543–50.
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R. Genome project data processing S: the sequence alignment/map format and SAMtools. Bioinformatics. 2009;25(16):2078–9.
Bolotin DA, Poslavsky S, Mitrophanov I, Shugay M, Mamedov IZ, Putintseva EV, Chudakov DM. MiXCR: software for comprehensive adaptive immunity profiling. Nat Methods. 2015;12(5):380–1.
Shugay M, Bagaev DV, Turchaninova MA, Bolotin DA, Britanova OV, Putintseva EV, Pogorelyy MV, Nazarov VI, Zvyagin IV, Kirgizova VI, et al. VDJtools: unifying post-analysis of T cell receptor repertoires. PLoS Comput Biol. 2015;11(11):e1004503.
Tumeh PC, Harview CL, Yearley JH, Shintaku IP, Taylor EJ, Robert L, Chmielowski B, Spasic M, Henry G, Ciobanu V, et al. PD-1 blockade induces responses by inhibiting adaptive immune resistance. Nature. 2014;515(7528):568–71.
Hunter JD. Matplotlib: a 2D graphics environment. Computing In Science & Engineering. 2007;9(3):90–5.
Kinoshita T, Muramatsu R, Fujita T, Nagumo H, Sakurai T, Noji S, Takahata E, Yaguchi T, Tsukamoto N, Kudo-Saito C, et al. Prognostic value of tumor-infiltrating lymphocytes differs depending on histological type and smoking habit in completely resected non-small-cell lung cancer. Ann Oncol. 2016;27(11):2117–23.
Lohr M, Edlund K, Botling J, Hammad S, Hellwig B, Othman A, Berglund A, Lambe M, Holmberg L, Ekman S, et al. The prognostic relevance of tumour-infiltrating plasma cells and immunoglobulin kappa C indicates an important role of the humoral immune response in non-small cell lung cancer. Cancer Lett. 2013;333(2):222–8.
Al-Shibli KI, Donnem T, Al-Saad S, Persson M, Bremnes RM, Busund LT. Prognostic effect of epithelial and stromal lymphocyte infiltration in non-small cell lung cancer. Clin Cancer Res. 2008;14(16):5220–7.
Hernandez-Prieto S, Romera A, Ferrer M, Subiza JL, Lopez-Asenjo JA, Jarabo JR, Gomez AM, Molina EM, Puente J, Gonzalez-Larriba JL, et al. A 50-gene signature is a novel scoring system for tumor-infiltrating immune cells with strong correlation with clinical outcome of stage I/II non-small cell lung cancer. Clin Transl Oncol. 2015;17(4):330–8.
Selitsky SR, Mose LE, Smith CC, Chai S, Hoadley KA, Dittmer DP, Moschos SJ, Parker JS, Vincent BG. Prognostic value of B cells in cutaneous melanoma. Genome medicine. 2019;11(1):36.
Turchaninova MA, Davydov A, Britanova OV, Shugay M, Bikos V, Egorov ES, Kirgizova VI, Merzlyak EM, Staroverov DB, Bolotin DA, et al. High-quality full-length immunoglobulin profiling with unique molecular barcoding. Nat Protoc. 2016;11(9):1599–616.
Cafri G, Yossef R, Pasetto A, Deniger DC, Lu YC, Parkhurst M, Gartner JJ, Jia L, Ray S, Ngo LT, et al. Memory T cells targeting oncogenic mutations detected in peripheral blood of epithelial cancer patients. Nat Commun. 2019;10(1):449.
Meng Q, Valentini D, Rao M, Maeurer M. KRAS RENAISSANCE(S) in tumor infiltrating B cells in pancreatic Cancer. Front Oncol. 2018;8:384.
Welinder C, Jirstrom K, Lehn S, Nodin B, Marko-Varga G, Blixt O, Danielsson L, Jansson B. Intra-tumour IgA1 is common in cancer and is correlated with poor prognosis in bladder cancer. Heliyon. 2016;2(8):e00143.
Fujimoto M, Yoshizawa A, Sumiyoshi S, Sonobe M, Kobayashi M, Koyanagi I, Aini W, Tsuruyama T, Date H, Haga H. Stromal plasma cells expressing immunoglobulin G4 subclass in non-small cell lung cancer. Hum Pathol. 2013;44(8):1569–76.
Biton J, Mansuet-Lupo A, Pecuchet N, Alifano M, Ouakrim H, Arrondeau J, Boudou-Rouquette P, Goldwasser F, Leroy K, Goc J, et al. TP53, STK11, and EGFR mutations predict tumor immune profile and the response to anti-PD-1 in lung adenocarcinoma. Clin Cancer Res. 2018;24(22):5710–23.
Zhang XC, Wang J, Shao GG, Wang Q, Qu X, Wang B, Moy C, Fan Y, Albertyn Z, Huang X, et al. Comprehensive genomic and immunological characterization of Chinese non-small cell lung cancer patients. Nat Commun. 2019;10(1):1772.
Tran E, Robbins PF, Lu YC, Prickett TD, Gartner JJ, Jia L, Pasetto A, Zheng Z, Ray S, Groh EM, et al. T-cell transfer therapy targeting mutant KRAS in Cancer. N Engl J Med. 2016;375(23):2255–62.
Schabath MB, Welsh EA, Fulp WJ, Chen L, Teer JK, Thompson ZJ, Engel BE, Xie M, Berglund AE, Creelan BC, et al. Differential association of STK11 and TP53 with KRAS mutation-associated gene expression, proliferation and immune surveillance in lung adenocarcinoma. Oncogene. 2016;35(24):3209–16.
Cheng H, Fan K, Luo G, Fan Z, Yang C, Huang Q, Jin K, Xu J, Yu X, Liu C. Kras(G12D) mutation contributes to regulatory T cell conversion through activation of the MEK/ERK pathway in pancreatic cancer. Cancer Lett. 2019;446:103–11.
Thorsson V, Gibbs DL, Brown SD, Wolf D, Bortone DS, Ou Yang TH, Porta-Pardo E, Gao GF, Plaisier CL, Eddy JA, et al. The immune landscape of Cancer. Immunity. 2018;48(4):812–30 e814.
Rubtsov AV, Rubtsova K, Kappler JW, Jacobelli J, Friedman RS, Marrack P. CD11c-expressing B cells are located at the T cell/B cell border in spleen and are potent APCs. J Immunol. 2015;195(1):71–9.
Rossetti RAM, Lorenzi NPC, Yokochi K, Rosa M, Benevides L, Margarido PFR, Baracat EC, Carvalho JP, Villa LL, Lepique AP. B lymphocytes can be activated to act as antigen presenting cells to promote anti-tumor responses. PLoS One. 2018;13(7):e0199034.
Chaft JE, Litvak A, Arcila ME, Patel P, D'Angelo SP, Krug LM, Rusch V, Mattson A, Coeshott C, Park B, et al. Phase II study of the GI-4000 KRAS vaccine after curative therapy in patients with stage I-III lung adenocarcinoma harboring a KRAS G12C, G12D, or G12V mutation. Clinical lung cancer. 2014;15(6):405–10.
Schuurman J, Perdok GJ, Gorter AD, Aalberse RC. The inter-heavy chain disulfide bonds of IgG4 are in equilibrium with intra-chain disulfide bonds. Mol Immunol. 2001;38(1):1–8.
van der Neut KM, Schuurman J, Losen M, Bleeker WK, Martinez-Martinez P, Vermeulen E, den Bleker TH, Wiegman L, Vink T, Aarden LA, et al. Anti-inflammatory activity of human IgG4 antibodies by dynamic fab arm exchange. Science. 2007;317(5844):1554–7.
Gunderson AJ, Coussens LM. B cells and their mediators as targets for therapy in solid tumors. Exp Cell Res. 2013;319(11):1644–9.
Karagiannis P, Gilbert AE, Josephs DH, Ali N, Dodev T, Saul L, Correa I, Roberts L, Beddowes E, Koers A, et al. IgG4 subclass antibodies impair antitumor immunity in melanoma. J Clin Invest. 2013;123(4):1457–74.
We are grateful to Alexey N Davydov (Central European Institute of Technology, Brno), Felix E. Frenkel (BostonGene LLC), and Maxim A. Chelushkin (BostonGene LLC) for their valuable comments and help with the data analysis, and to Michael Eisenstein for his great help with manuscript edition.
The work was supported by grant of the Ministry of Education and Science of the Russian Federation (14.W03.31.0005).
Ethics approval and consent to participate
Consent for publication
Authors declare no competing financial interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Non-silent mutation burden correlates with IGHG1/IGH proportion in LUAD subgroups. Figure S2. IgG1 clonality. Figure S3. Exploring IGH motifs linked to IgG1-mediated survival in KRASmut LUAD. Figure S4. Role of IgA expression in LUAD. Figure S5. Immunoglobulin isotypes and proportions in proximal proliferative LUAD. (DOCX 1198 kb)
About this article
Cite this article
Isaeva, O.I., Sharonov, G.V., Serebrovskaya, E.O. et al. Intratumoral immunoglobulin isotypes predict survival in lung adenocarcinoma subtypes. j. immunotherapy cancer 7, 279 (2019). https://doi.org/10.1186/s40425-019-0747-1