Pan-cancer γδ TCR analysis uncovers clonotype diversity and prognostic potential.

Gamma-delta T cells (γδ T cells) play a crucial role in both innate and adaptive immunity within tumors, yet their presence and prognostic value in cancer remain underexplored. This study presents a large-scale analysis of γδ T cell receptor (γδ TCR) reads from 11,000 tumor samples spanning 33 cancer types, utilizing the TRUST4 algorithm. Our findings reveal extensive diversity in γδ TCR clonality and gene expression, underscoring the potential of γδ T cells as prognostic biomarkers in various cancers. We further demonstrate the utility of TCR gamma (TRG) and delta (TRD) gene expression from standard RNA-sequencing (RNA-seq) data. This comprehensive dataset offers a valuable resource for advancing γδ T cell research, with implications for enhanced immunotherapy approaches or alternative therapeutic strategies. Additionally, our centralized database supports translational research into the therapeutic significance of γδ T cells.

New generalized metric based on branch length distance to compare B cell lineage trees.

The B cell lineage tree encapsulates the successive phases of B cell differentiation and maturation, transitioning from hematopoietic stem cells to mature, antibody-secreting cells within the immune system. Mathematically, this lineage can be conceptualized as an evolutionary tree, where each node represents a distinct stage in B cell development, and the edges reflect the differentiation pathways. To compare these lineage trees, a rigorous mathematical metric is essential. Analyzing B cell lineage trees mathematically and quantifying changes in lineage attributes over time necessitates a comparison methodology capable of accurately assessing and measuring these changes. Addressing the intricacies of multiple B cell lineage tree comparisons, this study introduces a novel metric that enhances the precision of comparative analysis. This metric is formulated on principles of metric theory and evolutionary biology, quantifying the dissimilarities between lineage trees by measuring branch length distance and weight. By providing a framework for systematically classifying lineage trees, this metric facilitates the development of predictive models that are crucial for the creation of targeted immunotherapy and vaccines. To validate the effectiveness of this new metric, synthetic datasets that mimic the complexity and variability of real B cell lineage structures are employed. We demonstrated the ability of the new metric method to accurately capture the evolutionary nuances of B cell lineages.

Immune repertoire profiling in myasthenia gravis.

Myasthenia gravis (MG) is the most frequent immune-mediated neurological disorder, characterized by fluctuating muscle weakness. Specific recognition of self-antigens by T-cell receptors (TCRs) and B-cell receptors (BCRs), coupled with T-B cell interactions, activates B cells to produce autoantibodies, which are critical for the initiation and perpetuation of MG. The immune repertoire comprises all functionally diverse T and B cells at a specific time point in an individual, reflecting the essence of immune selectivity. By sequencing the nucleotide sequences of TCRs and BCRs, it is possible to track individual T- and B-cell clones. This review delves into the generation of autoreactive TCRs and BCRs in MG and comprehensively examines the applications of immune repertoire sequencing in understanding disease pathogenesis, developing diagnostic and prognostic markers and informing targeted therapies. We also discuss the current limitations and future potential of this approach.

Limits on inferring T cell specificity from partial information.

A key challenge in molecular biology is to decipher the mapping of protein sequence to function. To perform this mapping requires the identification of sequence features most informative about function. Here, we quantify the amount of information (in bits) that T cell receptor (TCR) sequence features provide about antigen specificity. We identify informative features by their degree of conservation among antigen-specific receptors relative to null expectations. We find that TCR specificity synergistically depends on the hypervariable regions of both receptor chains, with a degree of synergy that strongly depends on the ligand. Using a coincidence-based approach to measuring information enables us to directly bound the accuracy with which TCR specificity can be predicted from partial matches to reference sequences. We anticipate that our statistical framework will be of use for developing machine learning models for TCR specificity prediction and for optimizing TCRs for cell therapies. The proposed coincidence-based information measures might find further applications in bounding the performance of pairwise classifiers in other fields.

Nucleos(T)ide Analogue Treatment Has a More Pronounced Impact on Immune Repertoires of CHB Patients Compared to HCC Patients.

Background: Nucleos(t)ide analogues (NAs) as the first-line treatment for chronic hepatitis B (CHB) have been shown to partially restore the antiviral immunity of the patients. However, hepatitis B virus (HBV) related hepatocellular carcinoma (HCC) patients have a relatively longer duration of HBV infection and lower level of HBV DNA. Whether NAs treatments have a different effect on their immune repertoires between CHB and HCC patients remains to be determined. Patients and Methods: In this study, 126 CHB patients and 85 HBV-related HCC patients who received or did not receive NAs treatment, as well as 361 healthy individuals were enrolled to analyze the effect of NAs treatment on T cell receptor ß chain (TCRß) and B cell receptor heavy chain (BCRh) repertoires in peripheral blood of the patients. Results: We found that after NAs therapy, the richness and evenness of TCRß and BCRh repertoires in CHB patients were significantly lower than those in untreated patients and healthy controls, while the diversity of TCRß and BCRh repertoires was stable in HCC patients. The alanine aminotransferase and HBV DNA levels were not correlated with the TCR or BCR diversity in CHB and HCC patients. Conclusion: The results suggest that NAs therapy could influence the overall T cell and B cell repertoires diversity in CHB patients but has minimal impact on HCC patients, indicating a significant difference in the potential to restore antiviral immunity between CHB and HCC patients by NAs treatment.

Similar usage of T-cell receptor ß-chain between tumor and adjacent normal tissue in hepatocellular carcinoma.

BACKGROUND: In this study, we comprehensively profiled the T-cell receptor (TCR) repertoire of the tumor and adjacent normal tissue in patients with HBV-associated hepatocellular carcinoma (HCC) and determined the baseline characteristics and clinical significance of TCR. METHODS: High-throughput sequencing was used to determine the profile of complementarity-determining region 3 (CDR3) of the TCR-ß chain variable (TRBV) in the tumor and normal tissue samples of 14 HCC patients. At the same time, TRBV diversity and differences in expression between tumor and normal tissues were investigated. The cumulative frequency of top 100 CDR3 (CF100), clonality, and Shannon entropy as indices to evaluate diversity, RESULTS: The diversity of TRBV CDR3 showed no significant difference between tumor and normal tissues. Of the 58 V gene segments in TRBV, TRBV16 and TRBV7-6 had a significantly higher frequency in the tumor group than in the normal group (p < 0.05). The frequency of 14 J gene segments showed no significant difference between tumor and normal tissues. In contrast, the frequency of 22 TRBVx/BJx combinations was significantly higher in the tumor than in the normal tissue. In addition, the length and type of TRBV CDR3 were similar in tumor and normal tissues, and a Gaussian distribution was observed in both groups. CONCLUSION: This study provided a large amount of information about the TCR lineage in HBV-associated HCC, laying the foundation for further research. In addition, the fact that the immune repertoire (TRBV CDR3) hardly differs between tumor and adjacent normal tissue provides a new clue for exploring the mechanism of the liver as an organ with immune privileges.

High throughput long-read sequencing of circulating lymphocytes of the evolutionarily distant sea lamprey reveals diversity and common elements of the variable lymphocyte receptor B (VLRB) repertoire.

The leucine-rich repeat-based variable lymphocyte receptor B (VLRB) antibody system of jawless vertebrates is capable of generating an antibody repertoire equal to or exceeding the diversity of antibody repertoires of jawed vertebrates. Unlike immunoglobulin-based immune repertoires, the VLRB repertoire diversity is characterized by variable lengths of VLRB encoding transcripts, rendering conventional immunoreceptor repertoire sequencing approaches unsuitable for VLRB repertoire sequencing. Here we demonstrate that long-read single-molecule real-time (SMRT) sequencing (PacBio) approaches permit the efficient large-scale assessment of the VLRB repertoire. We present a computational pipeline for sequence data processing and provide the first repertoire-based analysis of VLRB protein characteristics including properties of its subunits and regions of diversity within each structural leucine-rich repeat subunit. Our study provides a template to explore changes in the VLRB repertoire during immune responses and to establish large scale VLRB repertoire databases for computational approaches aimed at isolating monoclonal VLRB reagents for biomedical research and clinical applications.

Electroacupuncture Ameliorates Knee Osteoarthritis By Rebalancing T Cell Homeostasis as Revealed By Immune Repertoire (IR) Sequencing.

BACKGROUND: In this study, we used immune repertoire (IR) sequencing technology to profile the diversity of peripheral blood T cell receptors and used transcriptomics to profile the gene expression of peripheral blood neutrophil mRNA in patients with mild-moderate knee osteoarthritis (KOA) before and after electroacupuncture (EA) treatment. METHODS: An 8-week intervention with EA was performed on 3 subjects with KOA. IR sequencing of complementarity determining region 3 (CDR3) was performed using RNA extracted from peripheral blood T cells of KOA subjects prior to and at the end of the intervention, as well as healthy volunteers (controls) who matched the subjects in sex and age. Neutrophils were extracted from the plasma of healthy individuals, pretreatment patients, and posttreatment patients for further transcriptome sequencing. RESULTS: The D50, diversity index (DI), and Shannon entropy values of circulatory T-cells were significantly lower in pretreatment KOA patients compared to healthy controls. Posttreatment KOA samples displayed significant decreases in serum proinflammatory factors, IL-8 and IL-18 (P < 0.01), as well as a substantial reduction in serum matrix MMP-3 and MMP-13 (P < 0.01, P < 0.05). Transcriptome analysis revealed that the expression of CXCL2, IRF8, and PEAR1 (P < 0.05) was significantly higher in patients before the treatment than in the healthy population and was significantly down-regulated after the treatment. In contrast, the expression of SMPD3 (P < 0.05) showed the opposite trend. CONCLUSION: EA may alleviate KOA by rebalancing T-cell homeostasis and improving systemic inflammation. At the same time, EA treatment can significantly enhance TCR diversity, reduce levels of proinflammatory factors, and increase levels of anti-inflammatory factors, thereby achieving therapeutic effects.

Distinct adaptive immune receptor feature of adipose-derived mesenchymal stem cells (AD-MSCs) treatment of psoriasis.

Psoriasis (Ps) is one of the most common chronic inflammatory skin disorders with its pathogenesis correlated with dysregulated innate and adaptive system. Even though biological agents have advanced the treatment of psoriasis, however, there are huge limitations, like high adverse reactions and relapse rate. Therefore, it is of great interest in searching clinical resolutions with better safety and efficacy. In the current study, we utilized the adipose-derived mesenchymal stem cell (AD-MSCs) to treat moderate/severe cases of psoriasis in a single-arm clinical study. This AD-MSC treatment has proven to be clinically safe and effective. Interestingly, a trend of adaptome improvement, including increased diversity, elevated uCDR3s and decreased large clone after AD-MSC treatment in a short (2 weeks) and long (12 weeks) terms. In conclusion, allogenic AD-MSC treatment has shown a good safety and efficacy in treating Ps and can effectively improve the compromised adaptive immune system of Ps patients.

Computational estimation of clonal diversity in autoimmunity.

Diversity is the cornerstone of the adaptive immune system, crucial for its effectiveness against constantly evolving pathogens that pose threats to higher vertebrates. Accurately measuring and interpreting this diversity presents challenges for immunologists, as changes in diversity and clonotype composition can tip the balance between protective immunity and autoimmunity. In this review, we present the current methods commonly used to measure diversity from single-cell T-cell receptor and B-cell receptor sequencing. We also discuss two case studies where single-cell sequencing and diversity estimations have led to breakthroughs in autoimmune disease discovery and therapeutic innovation, and reflect upon the necessity and importance of accurately defining and measuring lymphocyte diversity in these contexts.

Reconstitution of peripheral blood T cell receptor ß immune repertoire in immune checkpoint inhibitors associated myocarditis.

PURPOSE: Immune checkpoint inhibitors (ICIs)-associated myocarditis was a rare yet severe complication observed in individuals undergoing immunotherapy. This study investigated the immune status and characteristics of patients diagnosed with ICIs- associated myocarditis. METHODS: A total of seven patients diagnosed with ICIs-associated myocarditis were included in the study, while five tumor patients without myocarditis were recruited as reference controls. Additionally, 30 healthy individuals were recruited as blank controls. Biochemical indices, electrocardiogram, and echocardiography measurements were obtained both prior to and following the occurrence of myocarditis. High-throughput sequencing of T cell receptor (TCR) was employed to assess the diversity and distribution characteristics of TCR CDR3 length, as well as the diversity of variable (V) and joining (J) genes of T lymphocytes in peripheral blood. RESULTS: In the seven patients with ICIs-associated myocarditis, Troponin T (TNT) levels exhibited a significant increase following myocarditis, while other parameters such as brain natriuretic peptide (BNP), QTc interval, and left ventricular ejection fraction (LVEF) did not show any significant differences. Through sequencing, it was observed that the diversity and uniformity of CDR3 in the ICIs-associated myocarditis patients were significantly diminished. Additionally, the distribution of CDR3 nucleotides deviated from normality, and variations in the utilization of V and J gene segments. CONCLUSION: The reconstitution of the TCR immune repertoire may play a pivotal role in the recognition of antigens in patients with ICIs-associated myocarditis.

The TCR ß chain CDR3 insights of bovine liver immune repertoire under heat stress.

Objective: The liver plays a dual role in regulating temperature and immune responses. Examining the influence of Heat stress (HS) on liver T cells contributes significantly to understanding the intricate interplay between the immune system and hepatic tissues under thermal stress. This study focused on investigating the characteristics of the T-cell receptor (TCR) ß chain CDR3 repertoire in bovine liver samples under both HS and pair-fed (PF) environmental conditions. Methods: Sequencing data from six samples sourced from the GEO database underwent annotation. Utilizing immunarch and VDJtool software, the study conducted comprehensive analyses encompassing basic evaluation, clonality assessment, immune repertoire comparison, diversity estimation, gene usage profiling, VJ gene segment pairing scrutiny, clonal tracking, and Kmers analysis. Results: All four TCR chains, namely α, ß, γ, and δ, were detected, with the α chains exhibiting the highest detection frequency, followed closely by the ß chains. The prevalence of αß TCRs in bovine liver samples underscored their crucial role in governing hepatic tissue's physiological functions. The TCR ß CDR3 repertoire showcased substantial inter-individual variability, featuring diverse clonotypes exhibiting distinct amino acid lengths. Intriguingly, HS cattle displayed heightened diversity and clonality, suggesting potential peripheral T cell migration into the liver under environmental conditions. Notably, differential VJ gene pairings were observed in HS cattle compared to the PF, despite individual variations in V and J gene utilization. Additionally, while most high-frequency amino acid 5-mers remained consistent between the HS and PF, GELHF and YDYHF were notably prevalent in the HS group. Across all samples, a prevalent trend of high-frequency 5mers skewed towards polar and hydrophobic amino acids was evident. Conclusion: This study elucidates the characteristics of liver TCR ß chain CDR3 repertoire under HS conditions, enhancing our understanding of HS implications.

Improved detection of tryptic immunoglobulin variable region peptides by chromatographic and gas-phase fractionation techniques.

The polyclonal repertoire of circulating antibodies potentially holds valuable information about an individual's humoral immune state. While bottom-up proteomics is well suited for serum proteomics, the vast number of antibodies and dynamic range of serum challenge this analysis. To acquire the serum proteome more comprehensively, we incorporated high-field asymmetric waveform ion-mobility spectrometry (FAIMS) or two-dimensional chromatography into standard trypsin-based bottom-up proteomics. Thereby, the number of variable region (VR)-related spectra increased 1.7-fold with FAIMS and 10-fold with chromatography fractionation. To match antibody VRs to spectra, we combined de novo searching and BLAST alignment. Validation of this approach showed that, as peptide length increased, the de novo accuracy decreased and BLAST performance increased. Through in silico calculations on antibody repository sequences, we determined the uniqueness of tryptic VR peptides and their suitability as antibody surrogate. Approximately one-third of these peptides were unique, and about one-third of all antibodies contained at least one unique peptide.

RAIN: a Machine Learning-based identification for HIV-1 bNAbs.

Broadly neutralizing antibodies (bNAbs) are promising candidates for the treatment and prevention of HIV-1 infection. Despite their critical importance, automatic detection of HIV-1 bNAbs from immune repertoire is still lacking. Here, we developed a straightforward computational method for Rapid Automatic Identification of bNAbs (RAIN) based on Machine Learning methods. In contrast to other approaches using one-hot encoding amino acid sequences or structural alignment for prediction, RAIN uses a combination of selected sequence-based features for accurate prediction of HIV-1 bNAbs. We demonstrate the performance of our approach on non-biased, experimentally obtained sequenced BCR repertoires from HIV-1 immune donors. RAIN processing leads to the successful identification of novel HIV-1 bNAbs targeting the CD4-binding site of the envelope glycoprotein. In addition, we validate the identified bNAbs using in vitro neutralization assay and we solve the structure of one of them in complex with the soluble native-like heterotrimeric envelope glycoprotein by single-particle cryo-electron microscopy (cryo-EM). Overall, we propose a method to facilitate and accelerate HIV-1 bNAbs discovery from non-selected immune repertoires.

Deciphering Autoimmune Diseases: Unveiling the Diagnostic, Therapeutic, and Prognostic Potential of Immune Repertoire Sequencing.

Autoimmune diseases (AIDs) are immune system disorders where the body exhibits an immune response to its own antigens, causing damage to its own tissues and organs. The pathogenesis of AIDs is incompletely understood. However, recent advances in immune repertoire sequencing (IR-seq) technology have opened-up a new avenue to study the IR. These studies have revealed the prevalence in IR alterations, potentially inducing AIDs by disrupting immune tolerance and thereby contributing to our comprehension of AIDs. IR-seq harbors significant potential for the clinical diagnosis, personalized treatment, and prognosis of AIDs. This article reviews the application and progress of IR-seq in diseases, such as multiple sclerosis, systemic lupus erythematosus, rheumatoid arthritis, and type 1 diabetes, to enhance our understanding of the pathogenesis of AIDs and offer valuable references for the diagnosis and treatment of AIDs.

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.

Elucidating immunological characteristics of the adenoma-carcinoma sequence in colorectal cancer patients in South Korea using a bioinformatics approach.

Colorectal cancer (CRC) is one of the top five most common and life-threatening malignancies worldwide. Most CRC develops from advanced colorectal adenoma (ACA), a precancerous stage, through the adenoma-carcinoma sequence. However, its underlying mechanisms, including how the tumor microenvironment changes, remain elusive. Therefore, we conducted an integrative analysis comparing RNA-seq data collected from 40 ACA patients who visited Dongguk University Ilsan Hospital with normal adjacent colons and tumor samples from 18 CRC patients collected from a public database. Differential expression analysis identified 21 and 79 sequentially up- or down-regulated genes across the continuum, respectively. The functional centrality of the continuum genes was assessed through network analysis, identifying 11 up- and 13 down-regulated hub-genes. Subsequently, we validated the prognostic effects of hub-genes using the Kaplan-Meier survival analysis. To estimate the immunological transition of the adenoma-carcinoma sequence, single-cell deconvolution and immune repertoire analyses were conducted. Significant composition changes for innate immunity cells and decreased plasma B-cells with immunoglobulin diversity were observed, along with distinctive immunoglobulin recombination patterns. Taken together, we believe our findings suggest underlying transcriptional and immunological changes during the adenoma-carcinoma sequence, contributing to the further development of pre-diagnostic markers for CRC.

Comprehensive analysis of juvenile idiopathic arthritis patients' immune characteristics based on bulk and single-cell sequencing data.

Background: The pathogenesis of juvenile idiopathic arthritis (JIA) is strongly influenced by an impaired immune system. However, the molecular mechanisms underlying its development and progression have not been elucidated. In this study, the computational methods TRUST4 were used to construct a T-cell receptor (TCR) and B-cell receptor (BCR) repertoire from the peripheral blood of JIA patients via bulk RNA-seq data, after which the clonality and diversity of the immune repertoire were analyzed. Results: Our findings revealed significant differences in the frequency of clonotypes between the JIA and healthy control groups in terms of the TCR and BCR repertoires. This work identified specific V genes and J genes in TCRs and BCRs that could be used to expand our understanding of JIA. After single-cell RNA analysis, the relative percentages of CD14 monocytes were significantly greater in the JIA group. Cell-cell communication analysis revealed the significant role of the MIF signaling pathway in JIA. Conclusion: In conclusion, this work describes the immune features of both the TCR and BCR repertoires under JIA conditions and provides novel insight into immunotherapy for JIA.

Immune checkpoint therapy responders display early clonal expansion of tumor infiltrating lymphocytes.

Immune checkpoint therapy (ICT) causes durable tumour responses in a subgroup of patients, but it is not well known how T cell receptor beta (TCRß) repertoire dynamics contribute to the therapeutic response. Using murine models that exclude variation in host genetics, environmental factors and tumour mutation burden, limiting variation between animals to naturally diverse TCRß repertoires, we applied TCRseq, single cell RNAseq and flow cytometry to study TCRß repertoire dynamics in ICT responders and non-responders. Increased oligoclonal expansion of TCRß clonotypes was observed in responding tumours. Machine learning identified TCRß CDR3 signatures unique to each tumour model, and signatures associated with ICT response at various timepoints before or during ICT. Clonally expanded CD8+ T cells in responding tumours post ICT displayed effector T cell gene signatures and phenotype. An early burst of clonal expansion during ICT is associated with response, and we report unique dynamics in TCRß signatures associated with ICT response.

A comparative analysis of TCR immune repertoire in COVID-19 patients.

The coronavirus disease 2019 (COVID-19) has merged as a global health threat since its outbreak in December 2019. Despite widespread recognition, there has been a paucity of studies focusing on the T cell receptor (TCR) bias in adaptive immunity induced by SARS-CoV-2. This research conducted a comparative analysis of the TCR immune repertoire to identify notable αß TCR bias sequences associated with the SARS-CoV-2 virus antigen. The present study encompassed 73 symptomatic COVID-19 patients, categorized as moderate/mild or severe/critical, along with 9 healthy controls. Our findings revealed specific TCR chains prominently utilized by moderate and severe patients, identified as TRAV30-J34-TRBV3-1-J2-7 and TRAV12-3-J6-TRBV28-J1-1, respectively. Additionally, our research explored critical TCR preferences in the bronchoalveolar lavage fluid (BALF) of COVID-19 patients at various disease stages. Indeed, monitoring the dynamics of immune repertoire changes in COVID-19 patients could serve as a crucial biomarker for predicting disease progression and recovery. Furthermore, the study explored TCR bias in both peripheral blood mononuclear cells (PBMCs) and BALF. The most common αß VJ pair observed in BALF was TRAV12-3-J18-TRBV7-6-J2-7. In addition, a comparative analysis with the VDJdb database indicated that the HLA-A*02:01 allele exhibited the widest distribution and highest frequency in COVID-19 patients across different periods. This comprehensive examination provided a global characterization of the TCR immune repertoire in COVID-19 patients, contributing significantly to our understanding of TCR bias induced by SARS-CoV-2.