Dysfunctional tumor-infiltrating Vδ1 + T lymphocytes in microsatellite-stable colorectal cancer.

Although γδ T cells are known to participate in immune dysregulation in solid tumors, their relevance to human microsatellite-stable (MSS) colorectal cancer (CRC) is still undefined. Here, using integrated gene expression analysis and T cell receptor sequencing, we characterized γδ T cells in MSS CRC, with a focus on Vδ1 + T cells. We identified Vδ1+ T cells with shared motifs in the third complementarity-determining region of the δ-chain, reflective of antigen recognition. Changes in gene and protein expression levels suggested a dysfunctional effector state of Vδ1+ T cells in MSS CRC, distinct from Vδ1+ T cells in microsatellite-instable (MSI). Interaction analysis highlighted an immunosuppressive role of fibroblasts in the dysregulation of Vδ1+ T cells in MSS CRC via the TIGIT-NECTIN2 axis. Blocking this pathway with a TIGIT antibody partially restored cytotoxicity of the dysfunctional Vδ1 phenotype. These results define an operative pathway in γδ T cells in MSS CRC.

De novo generation of SARS-CoV-2 antibody CDRH3 with a pre-trained generative large language model.

Artificial Intelligence (AI) techniques have made great advances in assisting antibody design. However, antibody design still heavily relies on isolating antigen-specific antibodies from serum, which is a resource-intensive and time-consuming process. To address this issue, we propose a Pre-trained Antibody generative large Language Model (PALM-H3) for the de novo generation of artificial antibodies heavy chain complementarity-determining region 3 (CDRH3) with desired antigen-binding specificity, reducing the reliance on natural antibodies. We also build a high-precision model antigen-antibody binder (A2binder) that pairs antigen epitope sequences with antibody sequences to predict binding specificity and affinity. PALM-H3-generated antibodies exhibit binding ability to SARS-CoV-2 antigens, including the emerging XBB variant, as confirmed through in-silico analysis and in-vitro assays. The in-vitro assays validate that PALM-H3-generated antibodies achieve high binding affinity and potent neutralization capability against spike proteins of SARS-CoV-2 wild-type, Alpha, Delta, and the emerging XBB variant. Meanwhile, A2binder demonstrates exceptional predictive performance on binding specificity for various epitopes and variants. Furthermore, by incorporating the attention mechanism inherent in the Roformer architecture into the PALM-H3 model, we improve its interpretability, providing crucial insights into the fundamental principles of antibody design.

Chemical complementarity of tumor resident, T-cell receptor CDR3s and renalase-1 correlates with increased melanoma survival.

Overexpression of the secretory protein renalase-1 negatively impacts the survival of melanoma and pancreatic cancer patients, while inhibition of renalase-1 signaling drives tumor rejection by promoting T-cell activation. Thus, we investigated the chemical complementarity between melanoma-resident, T-cell receptor (TCR) complementarity-determining region 3 (CDR3) amino acid sequences (AAs) and the renalase-1 protein. Increasing complementarity of TCR CDR3s to renalase-1 AAs, as assessed by a chemical complementarity scoring algorithm, was associated with improved overall survival (OS) in melanoma patients. The expression levels of several immune signature genes were significantly, positively correlated with increasing TCR CDR3-renalase-1 complementarity scores. Additionally, the survival association observed with high complementarity of TCR CDR3s to renalase-1 AAs was more robust in cases with low renalase-1 gene expression levels. Mapping of TCR CDR3-renalase-1 in silico interaction sites identified major epitope candidates including RP220, the signaling module of the renalase-1 protein, consistent with the fact that a monoclonal antibody to RP220 is a potent inhibitor of melanoma growth. These findings indicate that renalase-1 is a potential antigen for TCR recognition in melanoma and could be considered as a target for immunotherapy.

Comparison of "framework Shuffling" and "CDR Grafting" in humanization of a PD-1 murine antibody.

Introduction: Humanization is typically adopted to reduce the immunogenicity of murine antibodies generated by hybridoma technology when used in humans. Methods: Two different strategies of antibody humanization are popularly employed, including "complementarity determining region (CDR) grafting" and "framework (FR) shuffling" to humanize a murine antibody against human programmed death-1 (PD-1), XM PD1. In CDR-grafting humanization, the CDRs of XM PD-1, were grafted into the human FR regions with high homology to the murine FR counterparts, and back mutations of key residues were performed to retain the antigen-binding affinities. While in FR-shuffling humanization, a combinatorial library of the six murine CDRs in-frame of XM PD-1 was constructed to a pool of human germline FRs for high-throughput screening for the most favorable variants. We evaluated many aspects which were important during antibody development of the molecules obtained by the two methods, including antibody purity, thermal stability, binding efficacy, predicted humanness, and immunogenicity, along with T cell epitope prediction for the humanized antibodies. Results: While the ideal molecule was not achieved through CDR grafting in this particular instance, FR-shuffling proved successful in identifying a suitable candidate. The study highlights FR-shuffling as an effective complementary approach that potentially increases the success rate of antibody humanization. It is particularly noted for its accessibility to those with a biological rather than a computational background. Discussion: The insights from this comparison are intended to assist other researchers in selecting appropriate humanization strategies for drug development, contributing to broader application and understanding in the field.

A single-domain antibody library based on a stability-engineered human VH3 scaffold.

Using conventional immunoglobulin G (IgG) molecules as therapeutic agents presents several well-known disadvantages owing to their large size and structural complexity, negatively impacting development and production efficiency. Single-domain antibodies (sdAbs) are the smallest functional antibody format (~ 15 kDa) and represent a viable alternative to IgG in many applications. However, unlike natural single-domain antibodies, such as camelid VHH, the variable domains of conventional antibodies show poor physicochemical properties when expressed as sdAbs. This report identified stable sdAb variants of human VH3-23 from a framework region 2-randomized human VH library by phage display selection under thermal challenge. Synthetic complementarity determining region diversity was introduced to one of the selected variants with high thermal stability, expression level, and monomeric content to construct a human VH sdAb library. The library was validated by panning against a panel of antigens, and target-specific binders were identified and characterized for their affinity and biophysical properties. The results of this study suggest that a synthetic sdAb library based on a stability-engineered human VH scaffold could be a facile source of high-quality sdAb for many practical applications.

B cells expressing mutated IGHV1-69-encoded antigen receptors related to virus neutralization show lymphoma-like transcriptomes in patients with chronic HCV infection.

BACKGROUND: Chronic HCV infection leads to a complex interplay with adaptive immune cells that may result in B cell dyscrasias like cryoglobulinemia or lymphoma. While direct-acting antiviral therapy has decreased the incidence of severe liver damage, its effect on extrahepatic HCV manifestations such as B cell dyscrasias is still unclear. METHODS: We sequenced B cell receptor (BCR) repertoires in patients with chronic HCV mono-infection and patients with HCV with a sustained virological response (SVR) after direct-acting antiviral therapy. This data set was mined for highly neutralizing HCV antibodies and compared to a diffuse large B cell lymphoma data set. The TKO model was used to test the signaling strength of selected B-BCRs in vitro. Single-cell RNA sequencing of chronic HCV and HCV SVR samples was performed to analyze the transcriptome of B cells with HCV-neutralizing antigen receptors. RESULTS: We identified a B cell fingerprint with high richness and somatic hypermutation in patients with chronic HCV and SVR. Convergence to specific immunoglobulin genes produced high-connectivity complementarity-determining region 3 networks. In addition, we observed that IGHV1-69 CDR1 and FR3 mutations characterizing highly neutralizing HCV antibodies corresponded to recurrent point mutations found in clonotypic BCRs of high-grade lymphomas. These BCRs did not show autonomous signaling but a lower activation threshold in an in vitro cell model for the assessment of BCR signaling strength. Single-cell RNA sequencing revealed that B cells carrying these point mutations showed a persisting oncogenic transcriptome signature with dysregulation in signaling nodes such as CARD11, MALT1, RelB, MAPK, and NFAT. CONCLUSIONS: We provide evidence that lymphoma-like cells derive from the anti-HCV immune response. In many patients, these cells persist for years after SVR and can be interpreted as a mechanistic basis for HCV-related B cell dyscrasias and increased lymphoma risk even beyond viral elimination.

Characteristics of T cell premature senescence in maintenance hemodialysis patients.

BACKGROUND: Uremia-associated immunodeficiency, mainly characterized by T cell dysfunction, exists in patients on maintenance hemodialysis (MHD) and promotes systemic inflammation. However, T cell senescence, one of the causes of T cell dysfunction, has not been clearly revealed yet. In this cross-sectional research, we aimed to study the manifestation of T cell premature senescence in MHD patients and further investigate the associated clinical factors. METHODS: 76 MHD patients including 33 patients with cardiovascular diseases (CVD) and 28 patients with arteriovenous fistula (AVF) event history were enrolled in this study. Complementarity determining region 3 (CDR3) of T cell receptor (TCR) was analyzed by immune repertoire sequencing (IR-Seq). CD28- T cell subsets and expression of senescence marker p16 and p21 genes were detected by multicolor flow cytometry and RT-qPCR, respectively. RESULTS: MHD patients had significantly decreased TCR diversity (P < 0.001), increased CDR3 clone proliferation (P = 0.001) and a left-skewed CDR3 length distribution. The proportion of CD4 + CD28- T cells increased in MHD patients (P = 0.014) and showed a negative correlation with TCR diversity (P = 0.001). p16 but not p21 expression in T cells was up-regulated in MHD patients (P = 0.039). Patients with CVD exhibited increased expression of p16 and p21 genes (P = 0.010 and 0.004, respectively), and patients with AVF events showed further TCR diversity and evenness reduction (P = 0.002 and 0.017, respectively) compared to patients without the comorbidities. Moreover, age, average convection volume, total cholesterol, high-density lipoprotein cholesterol and transferrin saturation were associated with TCR diversity or CD4 + CD28- T cell proportion (P < 0.05). CONCLUSIONS: MHD patients undergo T cell premature senescence characterized by significant TCR diversity reduction and repertoire skew, as well as accumulation of the CD4 + CD28- subset and up-regulation of p16 gene. Patients with CVD or AVF events show higher level of immunosenescence. Furthermore, T cell senescence in MHD patients is associated with blood cholesterol and uremic toxin retention, suggesting potential intervention strategies in the future.

IGH Complementarity Determining Region-3-Cytomegalovirus Protein Chemical Complementarity Linked to Better Overall Survival Probabilities for Glioblastoma.

Cytomegalovirus (CMV) has long been thought to have an association with glioblastoma multiforme (GBM), although the exact role of CMV and any subsequent implications for treatment have yet to be fully understood. This study addressed whether IGH complementarity determining region-3 (CDR3)-CMV protein chemical complementarity, with IGH CDR3s representing both tumor resident and blood-sourced IGH recombinations, was associated with overall survival (OS) distinctions. IGH recombination sequencing reads were obtained from (a) the Clinical Proteomic Tumor Analysis Consortium, tumor RNAseq files; and (b) the cancer genome atlas, blood exome-derived files. The Adaptive Match web tool was used to calculate chemical complementarity scores (CSs) based on hydrophobic interactions, and those scores were used to group GBM cases and assess survival probabilities. We found a higher OS probability for cases whose hydrophobic IGH CDR3-CMV protein chemical complementarity scores (Hydro CSs) were in the upper 50th percentile for several CMV proteins, including UL99 and UL123, as well as for CSs based on known B cell epitopes representing these proteins. We also identified multiple immune signature genes, including CD79A and TNFRSF17, for which higher RNA expression was associated with higher Hydro CSs. Results were consistent with the idea that stronger immunoglobulin responses to CMV are associated with better OS probabilities for GBM.

Large-scale data mining of four billion human antibody variable regions reveals convergence between therapeutic and natural antibodies that constrains search space for biologics drug discovery.

The naïve human antibody repertoire has theoretical access to an estimated > 1015 antibodies. Identifying subsets of this prohibitively large space where therapeutically relevant antibodies may be found is useful for development of these agents. It was previously demonstrated that, despite the immense sequence space, different individuals can produce the same antibodies. It was also shown that therapeutic antibodies, which typically follow seemingly unnatural development processes, can arise independently naturally. To check for biases in how the sequence space is explored, we data mined public repositories to identify 220 bioprojects with a combined seven billion reads. Of these, we created a subset of human bioprojects that we make available as the AbNGS database (https://naturalantibody.com/ngs/). AbNGS contains 135 bioprojects with four billion productive human heavy variable region sequences and 385 million unique complementarity-determining region (CDR)-H3s. We find that 270,000 (0.07% of 385 million) unique CDR-H3s are highly public in that they occur in at least five of 135 bioprojects. Of 700 unique therapeutic CDR-H3, a total of 6% has direct matches in the small set of 270,000. This observation extends to a match between CDR-H3 and V-gene call as well. Thus, the subspace of shared ('public') CDR-H3s shows utility for serving as a starting point for therapeutic antibody design.

CDR grafting and site-directed mutagenesis approach for the generation and affinity maturation of Anti-CD20 nanobody.

BACKGROUND: Recently, new and advanced techniques have been adopted to design and produce nanobodies, which are used in diagnostic and immunotherapy treatments. Traditionally, nanobodies are prepared from camelid immune libraries that require animal treatments. However, such approaches require large library sizes and complicated selection procedures. The current study has employed CDR grafting and site-directed mutagenesis techniques to create genetically engineered nanobodies against the tumor marker CD20 (anti-CD20 nanobodies) used in leukemia treatment. METHODS AND RESULTS: In this study, we utilized the swapping method to graft CDRs from the VH Rituximab antibody to VHH CDRs. We aimed to enhance the binding affinity of the nanobodies by substituting the amino acids (Y101R-Y102R-Y107R) in the VHH-CDR3. To assess the binding capacity of the mutated nanobodies, we conducted an ELISA test. Moreover, through flow cytometry analysis, we compared the fluorescence intensity of the grafted CD20 and mutant nanobodies with that of the commercially available human anti-CD20 in Raji cells. The results showed a significant difference in the fluorescence intensity of the grafted nanobodies and mutant nanobodies when compared to the commercially available human anti-CD20. CONCLUSION: The approach we followed in this study makes it possible to create multiple anti-CD20 nanobodies with varying affinities without the need for extensive selection efforts. Additionally, our research has demonstrated that computational tools are highly reliable in designing functional nanobodies.

Identification of the immune repertoire of γδ T-cell receptors in psoriasis.

The γδ T cells are a subpopulation of T cells that are abundantly found in the skin and mucous membranes. Their reactivity to self-antigens and ability to secrete various cytokines make them a key component in psoriasis development. Although the correlation between the immune repertoire (IR) of γδ T-cell receptors and the occurrence and severity of psoriasis remains incompletely explored, high-throughput sequencing of γδ T cells has led to a deeper understanding of IR in psoriasis. This study investigated the differences between γδ T cells in patients with psoriasis and healthy controls. The γδ T cells were identified via immunofluorescence staining and a correlation analysis was performed according to the psoriasis area and severity index (PASI) scores. The IR sequencing method was used to detect IR in the γδ T-cell receptors. The findings demonstrated more skin γδ T cells in patients with psoriasis, which were positively correlated with the PASI score. There were subtle differences in most variable (V), diversity (D) and joining (J) gene segments and VJ/VDJ combination segments between patients with psoriasis and healthy controls. However, a higher diversity of complementarity-determining region 3 (CDR3) was observed in patients with psoriasis. In summary, the IR of skin γδ T cells was significantly altered in patients with psoriasis, and the diversity in the cell's CDR3 population is a promising biomarker for assessment of psoriasis severity.

Local chicken breeds exhibit abundant TCR-V segments but similar repertoire diversity.

The thymus-derived lymphocytes of jawed vertebrates have four T-cell receptor (TCR) chains that play a significant role in immunity. As chickens have commercial value, their immune systems require a great deal of attention. Local chicken breeds are an essential part of poultry genetic resources in China. Here, we used high-throughput sequencing to analyze the TCRα and TCRß repertoires and their relative expression levels in the native chicken breeds Baier Buff, Longyou Partridge, Xiaoshan, and Xianju. We found that TCR Vα and TCR Vß were expressed and included 17, 19, 17, and six segments of the Vα2, Vα3, Vß1, and Vß2 subgroups, respectively. V-J pairing was biased; Jα11 was utilized by nearly all Vα segments and was the most commonly used. Breed-specific V segments and V-J pairings were detected as well. The results of the principal coordinate analysis (PCoA) as well as the V-J pairing and CDR3 diversity analyses suggested that the four local chicken breeds did not significantly differ in terms of TCR diversity. Hence, they expressed not significant differentiation, and they are rich genetic resources for the development and utilization of immune-related poultry breeding.

T-cell receptor diversity and allergen sensitivity in childhood asthma and atopic dermatitis.

BACKGROUND: Childhood allergies of asthma and atopic dermatitis (AD) involve an overactive T-cell immune response triggered by allergens. However, the impact of T-cell receptor (TCR) repertoires on allergen sensitization and their role in mediating different phenotypes of asthma and AD in early childhood remains unclear. METHODS: A total of 78 children, comprising 26 with asthma alone, 26 with AD alone, and 26 healthy controls (HC), were enrolled. TCR repertoire profiles were determined using a unique molecular identifier system for next-generation sequencing. Integrative analyses of their associations with allergen-specific IgE levels and allergies were performed. RESULTS: The diversity in TCR alpha variable region (TRAV) genes of TCR repertoires and complementarity determining region 3 (CDR3) clonality in TRAV/TRBV (beta) genes were significantly higher in children with AD compared with those with asthma and HC (p < .05). Compared with HC, the expression of TRAV13-1 and TRAV4 genes was significantly higher in both asthma and AD (p < .05), with a significant positive correlation with mite-specific IgE levels (p < .01). In contrast, TRBV7-9 gene expression was significantly lower in both asthma and AD (p < .01), with this gene showing a significant negative correlation with mite-specific IgE levels (p < .01). Furthermore, significantly higher TRAV8-3 gene expression, positively correlated with food-specific IgE levels, was found in children with AD compared with those with asthma (p < .05). CONCLUSION: Integrated TCR repertoires analysis provides clinical insights into the diverse TCR genes linked to antigen specificity, offering potential for precision immunotherapy in childhood allergies.

An integrated technology for quantitative wide mutational scanning of human antibody Fab libraries.

Antibodies are engineerable quantities in medicine. Learning antibody molecular recognition would enable the in silico design of high affinity binders against nearly any proteinaceous surface. Yet, publicly available experiment antibody sequence-binding datasets may not contain the mutagenic, antigenic, or antibody sequence diversity necessary for deep learning approaches to capture molecular recognition. In part, this is because limited experimental platforms exist for assessing quantitative and simultaneous sequence-function relationships for multiple antibodies. Here we present MAGMA-seq, an integrated technology that combines multiple antigens and multiple antibodies and determines quantitative biophysical parameters using deep sequencing. We demonstrate MAGMA-seq on two pooled libraries comprising mutants of nine different human antibodies spanning light chain gene usage, CDR H3 length, and antigenic targets. We demonstrate the comprehensive mapping of potential antibody development pathways, sequence-binding relationships for multiple antibodies simultaneously, and identification of paratope sequence determinants for binding recognition for broadly neutralizing antibodies (bnAbs). MAGMA-seq enables rapid and scalable antibody engineering of multiple lead candidates because it can measure binding for mutants of many given parental antibodies in a single experiment.

Invariant γδTCR natural killer-like effector T cells in the naked mole-rat.

The naked mole-rat (Heterocephalus glaber) is a long-lived rodent species showing resistance to the development of cancer. Although naked mole-rats have been reported to lack natural killer (NK) cells, γδ T cell-based immunity has been suggested in this species, which could represent an important arm of the immune system for antitumor responses. Here, we investigate the biology of these unconventional T cells in peripheral tissues (blood, spleen) and thymus of the naked mole-rat at different ages by TCR repertoire profiling and single-cell gene expression analysis. Using our own TCR annotation in the naked mole-rat genome, we report that the γδ TCR repertoire is dominated by a public invariant Vγ4-2/Vδ1-4 TCR, containing the complementary-determining-region-3 (CDR3)γ CTYWDSNYAKKLF / CDR3δ CALWELRTGGITAQLVF that are likely generated by short-homology-repeat-driven DNA rearrangements. This invariant TCR is specifically found in γδ T cells expressing genes associated with NK cytotoxicity and is generated in both the thoracic and cervical thymus of the naked mole-rat until adult life. Our results indicate that invariant Vγ4-2/Vδ1-4 NK-like effector T cells in the naked mole-rat can contribute to tumor immunosurveillance by γδ TCR-mediated recognition of a common molecular signal.

Understanding the Specific Implications of Amino Acids in the Antibody Development.

As the demand for immunotherapy to treat and manage cancers, infectious diseases and other disorders grows, a comprehensive understanding of amino acids and their intricate role in antibody engineering has become a prime requirement. Naturally produced antibodies may not have the most suitable amino acids at the complementarity determining regions (CDR) and framework regions, for therapeutic purposes. Therefore, to enhance the binding affinity and therapeutic properties of an antibody, the specific impact of certain amino acids on the antibody's architecture must be thoroughly studied. In antibody engineering, it is crucial to identify the key amino acid residues that significantly contribute to improving antibody properties. Therapeutic antibodies with higher binding affinity and improved functionality can be achieved through modifications or substitutions with highly suitable amino acid residues. Here, we have indicated the frequency of amino acids and their association with the binding free energy in CDRs. The review also analyzes the experimental outcome of two studies that reveal the frequency of amino acids in CDRs and provides their significant correlation between the outcomes. Additionally, it discusses the various bond interactions within the antibody structure and antigen binding. A detailed understanding of these amino acid properties should assist in the analysis of antibody sequences and structures needed for designing and enhancing the overall performance of therapeutic antibodies.

Adaptive immune receptor germline gene variation.

Recognition of antigens by T cell receptors (TCRs) and B cell receptors (BCRs) is a key step in lymphocyte activation. T and B cells mediate adaptive immune responses, which protect us against infections and provide immunological memory, and also, in some instances, drive pathogenic responses in autoimmune diseases. TCRs and BCRs are encoded within loci that are known to be genetically diverse. However, the extent and functional impact of this variation, both in humans and model animals used in immunological research, remain largely unknown. Experimental and genetic evidence has demonstrated that the complementarity determining regions 1 and 2 (HCDR1 and HCDR2), encoded by the variable (V) region of TCRs and BCRs, also often make critical contacts with the targeted antigen. Thus, knowledge about allelic variation in the genes encoding TCRs and BCRs is critically important for understanding adaptive immune responses in outbred populations and to define responder and non-responder phenotypes.

High-throughput quantitative assessments of the chemical complementarity of celiac disease-related IGH CDR3s and a gliadin epitope.

The long-term value of efficient antigen discovery includes gaining insights into the variety of potential cancer neoantigens, effective vaccines lacking adverse effects, and adaptive immune receptor (IR) targets for blocking adaptive IR-antigen interactions in autoimmunity. While the preceding goals have been partially addressed via big data approaches to HLA (human leukocyte antigen)-epitope binding, there has been little such progress in the big data setting for adaptive IR-epitope binding. This delay in progress for the latter is likely due to, among other things, the much more complicated adaptive IR repertoire in an individual compared to individual HLA alleles. Thus, results described here represent the application of an algorithm for efficient assessment of immunoglobulin heavy chain complementarity determining region-3 (IGH CDR3)-gliadin epitope interactions, with a focus on epitopes known to be associated with an immune response in celiac disease. The hydrophobic, chemical complementarity between celiac case IGH CDR3s and known celiac epitopes was found to be greater in comparison to the hydrophobic, chemical complementarity between the same celiac case IGH CDR3s and a series of control epitopes. Thus, the approaches indicated here likely offer guidance for the development of conveniently applied algorithms for antigen verification and discovery.

[Biological characteristics and clinical significance of stereotyped B-cell receptor in chronic lymphocytic leukemia].

Chronic lymphocytic leukemia (CLL) is the most common leukemia in the Western adults, although the incidence of CLL is relatively low in Asian populations. However, with the aging population, the incidence of CLL is increasing in China. The interaction between CLL cells and the microenvironment plays a crucial role in the recognition of antigens by the B-cell receptor immunoglobulin (BCR IG). The mutational status of the immunoglobulin heavy variable region (IGHV) is a classical prognostic marker for CLL. Over 40% of CLL patients exhibit biased usage of IGHV and highly similar amino acid sequences in the heavy complementarity-determining region 3 (HCDR3), known as the BCR stereotypy. Different subgroups of stereotyped BCR exhibit distinct biological and clinical features. Among them, subset #2 with mutated IGHV and poor prognosis, as well as the subset #8 with a high risk of Richter transformation, have been recommended by the European Research Initiative on CLL to be included in clinical reports on IGHV mutational status. This review summarizes the definition, distribution, biological characteristics, and clinical significance of clonality patterns of the BCR in CLL.

Identification of apolipoprotein B-reactive CDR3 motifs allows tracking of atherosclerosis-related memory CD4+T cells in multiple donors.

Introduction: Atherosclerosis is a major pathological condition that underlies many cardiovascular diseases (CVDs). Its etiology involves breach of tolerance to self, leading to clonal expansion of autoreactive apolipoprotein B (APOB)-reactive CD4+T cells that correlates with clinical CVD. The T-cell receptor (TCR) sequences that mediate activation of APOB-specific CD4+T cells are unknown. Methods: In a previous study, we had profiled the hypervariable complementarity determining region 3 (CDR3) of CD4+T cells that respond to six immunodominant APOB epitopes in most donors. Here, we comprehensively analyze this dataset of 149,065 APOB-reactive and 199,211 non-reactive control CDR3s from six human leukocyte antigen-typed donors. Results: We identified 672 highly expanded (frequency threshold > 1.39E-03) clones that were significantly enriched in the APOB-reactive group as compared to the controls (log10 odds ratio ≥1, Fisher's test p < 0.01). Analysis of 114,755 naïve, 91,001 central memory (TCM) and 29,839 effector memory (TEM) CDR3 sequences from the same donors revealed that APOB+ clones can be traced to the complex repertoire of unenriched blood T cells. The fraction of APOB+ clones that overlapped with memory CDR3s ranged from 2.2% to 46% (average 16.4%). This was significantly higher than their overlap with the naïve pool, which ranged from 0.7% to 2% (average 1.36%). CDR3 motif analysis with the machine learning-based in-silico tool, GLIPHs (grouping of lymphocyte interactions by paratope hotspots), identified 532 APOB+ motifs. Analysis of naïve and memory CDR3 sequences with GLIPH revealed that ~40% (209 of 532) of these APOB+ motifs were enriched in the memory pool. Network analysis with Cytoscape revealed extensive sharing of the memory-affiliated APOB+ motifs across multiple donors. We identified six motifs that were present in TCM and TEM CDR3 sequences from >80% of the donors and were highly enriched in the APOB-reactive TCR repertoire. Discussion: The identified APOB-reactive expanded CD4+T cell clones and conserved motifs can be used to annotate and track human atherosclerosis-related autoreactive CD4+T cells and measure their clonal expansion.