High throughput sequencing reveals alterations in B cell receptor repertoires associated with the progression of hepatic cirrhosis to hepatocellular carcinoma

BACKGROUND

Hepatocellular carcinoma (HCC) is developed as a consequence of chronic liver cirrhosis, and both diseases are difficult to diagnose and differentiate. Accurate noninvasive biomarkers for HCC and liver cirrhosis are urgently needed.

METHODS

Here we used high-throughput sequencing to characterize the B cell receptor (BCR) repertoires from 36 HCC tumor samples and 10 liver cirrhosis (LC) tissue biopsies to understand the immune alterations during hepatic carcinogenesis.

RESULTS

The principal components analysis (PCA) showed that the pattern of BCR in HCC was distinct from that in LC. As measured by Clonality and Shannon indexes, the diversity of BCR repertoire was significantly lower in HCC than in LC (P < 0.01).

CONCLUSION

Our results corroborated that the BCR diversity and composition could be closely correlated with hepatic carcinogenesis. And BCR repertoire may be used to predict the progression of HCC and design targeting immunotherapy in the near future.

Comparative analysis of the immune repertoire between peripheral blood and bone marrow fluids in those infected by EBV and immunodeficiency: A retrospective case study

High-throughput immune repertoire (IR) sequencing provides direct insight into the diversity of B cell receptor (BCR) and T cell receptor (TCR), with great potential to revolutionize the diagnosis, monitoring, and prevention of immune system-related disorders. In this study, multiplex PCR was applied to amplify the complementarity-determining regions of BCR and TCR, followed by comprehensive analysis by high-throughput sequencing. We compare the TCR (BCR) of bone marrow fluid (BMF) and peripheral blood (PB) samples from 17 patients in the Epstein-Barr and immunodeficiency groups, respectively. Our study shows that the diversity of the IR of blood samples is very similar to that of bone marrow samples statistically. However, the distributions of the monoclonal genes are significantly different in these 2 samples of most patients. This suggests that the BMFs can be replaced by the PB samples in diversity detection of IR to monitor the immune status of the body, while the detection of the BMFs is unreplaceable when the monoclonal change occurs. We used high-throughput sequencing to assess the TCR and BCR of the patients and provide a basis for the clinical analysis of PB and bone marrow samples and selection of disease diagnosis and monitoring methods.

Gut Microbiota Defines Functional Direction of Colonic Regulatory T Cells with Unique TCR Repertoires

Intestinal microbiota and selected strains of commensal bacteria influence regulatory T (Treg) cell functionality in the colon. Nevertheless, whether and how microbiota changes the transcriptome profile and TCR specificities of colonic Tregs remain to be precisely defined. In this study, we have employed single-cell RNA sequencing and comparatively analyzed colonic Tregs from specific pathogen-free and germ-free (GF) mice. We found that microbiota shifts the activation trajectory of colonic Tregs toward a distinct phenotypic subset enriched in specific pathogen-free but not in GF mice. Moreover, microbiota induced the expansion of specific Treg clonotypes with shared transcriptional specificities. The microbiota-induced subset of colonic Tregs, identified as PD-1- CXCR3+ Tregs, displayed enhanced suppressive capabilities compared with colonic Tregs derived from GF mice, enhanced production of IL-10, and were the primary regulators of enteric inflammation in dextran sodium sulfate-induced colitis. These findings identify a hitherto unknown gut microbiota and immune cell interaction module that could contribute to the development of a therapeutic modality for intestinal inflammatory diseases.

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.

Benchmarking and integrating human B-cell receptor genomic and antibody proteomic profiling

Immunoglobulins (Ig), which exist either as B-cell receptors (BCR) on the surface of B cells or as antibodies when secreted, play a key role in the recognition and response to antigenic threats. The capability to jointly characterize the BCR and antibody repertoire is crucial for understanding human adaptive immunity. From peripheral blood, bulk BCR sequencing (bulkBCR-seq) currently provides the highest sampling depth, single-cell BCR sequencing (scBCR-seq) allows for paired chain characterization, and antibody peptide sequencing by tandem mass spectrometry (Ab-seq) provides information on the composition of secreted antibodies in the serum. Yet, it has not been benchmarked to what extent the datasets generated by these three technologies overlap and complement each other. To address this question, we isolated peripheral blood B cells from healthy human donors and sequenced BCRs at bulk and single-cell levels, in addition to utilizing publicly available sequencing data. Integrated analysis was performed on these datasets, resolved by replicates and across individuals. Simultaneously, serum antibodies were isolated, digested with multiple proteases, and analyzed with Ab-seq. Systems immunology analysis showed high concordance in repertoire features between bulk and scBCR-seq within individuals, especially when replicates were utilized. In addition, Ab-seq identified clonotype-specific peptides using both bulk and scBCR-seq library references, demonstrating the feasibility of combining scBCR-seq and Ab-seq for reconstructing paired-chain Ig sequences from the serum antibody repertoire. Collectively, our work serves as a proof-of-principle for combining bulk sequencing, single-cell sequencing, and mass spectrometry as complementary methods towards capturing humoral immunity in its entirety.

Composition and diversity analysis of the TCR CDR3 repertoire in patients with idiopathic orbital inflammation using high-throughput sequencing

BACKGROUND

Idiopathic orbital inflammation (IOI) is a nonspecific orbital inflammatory disease with the third highest prevalence among orbital diseases, and its pathogenesis is associated with T-cell-mediated immune responses. This study aimed to investigate the differences in T-cell receptor (TCR) expression between IOI patients and healthy subjects by high-throughput sequencing and to characterize TCR expression in patients with IOI and with respect to glucocorticoid response.

METHODS

A total of 19 subjects were enrolled in this study and were divided into the idiopathic orbital inflammation group (IOI group, n = 13) and the healthy control group (HC group, n = 6), and within the IOI group were further divided into the glucocorticoid therapy sensitive group (IOI(EF) group, n = 6) and the glucocorticoid therapy ineffective group (IOI(IN) group, n = 7) based on the degree of effectiveness to glucocorticoid therapy. High-throughput TCR sequencing was performed on peripheral blood mononuclear cells of IOI patients and healthy control individuals using 5' RACE technology combined with Unique Identifier (UID) digital tag correction technology. The TCR CDR3 region diversity, sharing patterns, and differential sequences between the IOI and HC groups, and between the IOI(EF) and IOI(IN) groups were analyzed.

RESULTS

It was found that the diversity of TCR CDR3 in the IOI group was significantly lower than that in the HC group, and the frequency of V gene use was significantly different between groups. The diversity of TCR CDR3 in patients in the IOI(EF) group was significantly lower than that in patients in the IOI(IN) group, and the frequency of V and J gene use was significantly different between the IOI(EF) group and the IOI(IN) group. Additionally, we found 133 nucleotide sequences shared in all IOI samples and screened two sequences with higher expression from them.

CONCLUSIONS

Our results suggested that abnormal clonal expansion of specific T-cells exists in IOI patients and that TCR diversity may had an impact on the prognosis of glucocorticoid-treated IOI. This study may contribute to a better understanding of the immune status of IOI and provide new insights for T-cell -associated IOI pathogenesis, diagnosis and treatment prediction.

The dynamic TRβ/IGH CDR3 repertoire features in patients with liver transplantation

OBJECTIVE

At present, little is known about the immune mechanism of liver transplantation caused by decompensated cirrhosis. Lymphocytes play an essential important role in the immune rejection of liver transplantation. In this study, we aimed to comprehensively analyze changes in complementary determinant 3 (CDR3) repertoire of T cell receptor β chain (TRβ) and immunoglobulin heavy chain (IGH) in liver transplantation patients and healthy controls (HC).

METHODS

High-throughput sequencing technology was used to study the characteristics of TRβ/IGH CDR3 repertoire, and identify the amino acid sequences of TRβ and IGH associated with liver transplantation patients and HC.

RESULTS

We found that some TRβ and IGH CDR3 repertoire characteristics differed between liver transplant patients and HC. The diversity of TRβ CDR3 increased in the liver transplantation group. First and seven days after live transplantation patients showed a lower degree of T cell clone amplification compared to the HC group. The CDR3 repertoire of the TRβ/IGH chain was certainly biased in the use of some V, D, and J gene segments, TRβ/IGH V-J combined frequency was also skewed and TRβ CDR3 clonotypes were shared at a higher degree in the liver transplantation patients. Importantly, one amino acid sequence in the decompensated cirrhosis group was significantly higher than that in the healthy group. It should be noted that the frequency of some CDR3 sequences is closely correlated with the different stages of liver transplantation, and these sequences may play a key role in liver transplantation.

CONCLUSION

Based on the above results, we can better understand the dynamic changes of TCβ/IGH CDR3 repertoire in patients during liver transplantation.

Application and challenges of TCR and BCR sequencing to investigate T- and B-cell clonality in elastase-induced experimental murine abdominal aortic aneurysm

Background

An abdominal aortic aneurysm (AAA) is a life-threatening cardiovascular disease. Although its pathogenesis is still poorly understood, recent evidence suggests that AAA displays autoimmune disease characteristics. Particularly, T cells responding to AAA-related antigens in the aortic wall may contribute to an initial immune response. Single-cell RNA (scRNA) T cell receptor (TCR) and B cell receptor (BCR) sequencing is a powerful tool for investigating clonality. However, difficulties such as limited numbers of isolated cells must be considered during implementation and data analysis, making biological interpretation challenging. Here, we perform a representative single-cell immune repertoire analysis in experimental murine AAA and show a reliable bioinformatic processing pipeline highlighting opportunities and limitations of this approach.

Methods

We performed scRNA TCR and BCR sequencing of isolated lymphocytes from the infrarenal aorta of male C57BL/6J mice 3, 7, 14, and 28 days after AAA induction via elastase perfusion of the aorta. Sham-operated mice at days 3 and 28 and non-operated mice served as controls.

Results

Comparison of complementarity-determining region (CDR3) length distribution of 179 B cells and 796 T cells revealed neither differences between AAA and control nor between the disease stages. We found no clonal expansion of B cells in AAA. For T cells, we identified several clones in 11 of 16 AAA samples and one of eight control samples. Immune receptor repertoire comparison indicated that only a few clones were shared between the individual AAA samples. The most frequently used V-genes in the TCR beta chain in AAA were TRBV3, TRBV19, and the splicing variant TRBV12-2 + TRBV13-2.

Conclusion

We found no clonal expansion of B cells but evidence for clonal expansion of T cells in elastase-induced AAA in mice. Our findings imply that a more precise characterization of TCR and BCR distribution requires a more extensive number of lymphocytes to prevent undersampling and potentially detect rare clones. Thus, further experiments are necessary to confirm our findings. In summary, this paper examines TCR and BCR sequencing results, identifies limitations and pitfalls, and offers guidance for future studies.

The technological landscape and applications of single-cell multi-omics

Single-cell multi-omics technologies and methods characterize cell states and activities by simultaneously integrating various single-modality omics methods that profile the transcriptome, genome, epigenome, epitranscriptome, proteome, metabolome and other (emerging) omics. Collectively, these methods are revolutionizing molecular cell biology research. In this comprehensive Review, we discuss established multi-omics technologies as well as cutting-edge and state-of-the-art methods in the field. We discuss how multi-omics technologies have been adapted and improved over the past decade using a framework characterized by optimization of throughput and resolution, modality integration, uniqueness and accuracy, and we also discuss multi-omics limitations. We highlight the impact that single-cell multi-omics technologies have had in cell lineage tracing, tissue-specific and cell-specific atlas production, tumour immunology and cancer genetics, and in mapping of cellular spatial information in fundamental and translational research. Finally, we discuss bioinformatics tools that have been developed to link different omics modalities and elucidate functionality through the use of better mathematical modelling and computational methods.

Transcriptional and Clonal Characterization of Cytotoxic T Cells in Crescentic Glomerulonephritis

SIGNIFICANCE STATEMENT

T-cell infiltration is a hallmark of crescentic GN (cGN), often caused by ANCA-associated vasculitis. Pathogenic T-cell subsets, their clonality, and downstream effector mechanisms leading to kidney injury remain to be fully elucidated. Single-cell RNA sequencing and T-cell receptor sequencing revealed activated, clonally expanded cytotoxic CD4 + and CD8 + T cells in kidneys from patients with ANCA-associated cGN. In experimental cGN, kidney-infiltrating CD8 + T cells expressed the cytotoxic molecule, granzyme B (GzmB), which induced apoptosis in renal tissue cells by activation of procaspase-3, and aggravated disease pathology. These findings describe a pathogenic function of (clonally expanded) cytotoxic T cells in cGN and identify GzmB as a mediator and potential therapeutic target in immune-mediated kidney disease.

BACKGROUND

Crescentic GN (cGN) is an aggressive form of immune-mediated kidney disease that is an important cause of end stage renal failure. Antineutrophilic cytoplasmic antibody (ANCA)-associated vasculitis is a common cause. T cells infiltrate the kidney in cGN, but their precise role in autoimmunity is not known.

METHODS

Combined single-cell RNA sequencing and single-cell T-cell receptor sequencing were conducted on CD3 + T cells isolated from renal biopsies and blood of patients with ANCA-associated cGN and from kidneys of mice with experimental cGN. Functional and histopathological analyses were performed with Cd8a-/- and GzmB-/- mice.

RESULTS

Single-cell analyses identified activated, clonally expanded CD8 + and CD4 + T cells with a cytotoxic gene expression profile in the kidneys of patients with ANCA-associated cGN. Clonally expanded CD8 + T cells expressed the cytotoxic molecule, granzyme B (GzmB), in the mouse model of cGN. Deficiency of CD8 + T cells or GzmB ameliorated the course of cGN. CD8 + T cells promoted macrophage infiltration and GzmB activated procaspase-3 in renal tissue cells, thereby increasing kidney injury.

CONCLUSIONS

Clonally expanded cytotoxic T cells have a pathogenic function in immune-mediated kidney disease.

Insufficient epitope-specific T cell clones are responsible for impaired cellular immunity to inactivated SARS-CoV-2 vaccine in older adults

Aging is a critical risk factor for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) vaccine efficacy. The immune responses to inactivated vaccine for older adults, and the underlying mechanisms of potential differences to young adults, are still unclear. Here we show that neutralizing antibody production by older adults took a longer time to reach similar levels in young adults after inactivated SARS-CoV-2 vaccination. We screened SARS-CoV-2 variant strains for epitopes that stimulate specific CD8 T cell response, and older adults exhibited weaker CD8 T-cell-mediated responses to these epitopes. Comparison of lymphocyte transcriptomes from pre-vaccinated and post-vaccinated donors suggested that the older adults had impaired antigen processing and presentation capability. Single-cell sequencing revealed that older adults had less T cell clone expansion specific to SARS-CoV-2, likely due to inadequate immune receptor repertoire size and diversity. Our study provides mechanistic insights for weaker response to inactivated vaccine by older adults and suggests the need for further vaccination optimization for the old population.

Linking antigen specific T-cell dynamics in a microfluidic chip to single cell transcription patterns

A new non-invasive screening profile has been realized that can aid in determining T-cell activation state at single-cell level. Production of activated T-cells with good specificity and stable proliferation is greatly beneficial for advancing adoptive immunotherapy as innate immunological cells are not effective in recognizing and eliminating cancer as expected. The screening method is realized by relating intracellular Ca²⁺ intensity and motility of T-cells interacting with APC (Antigen Presenting Cells) in a microfluidic chip. The system is tested using APC pulsed with OVA_257-264 peptide and its modified affinities (N4, Q4, T4 and V4), and the T-cells from OT-1 mice. In addition, single cell RNA sequencing reveals the activation states of the cells and the clusters from the derived profiles can be indicative of the T-cell activation state. The presented system here can be versatile for a comprehensive application to proceed with T-cell-based immunotherapy and screen the antigen-specific T-cells with excellent efficiency and high proliferation.

B cell-related tertiary lymphoid structure may exert inhibitory effects on lung adenocarcinoma and SARS-COV-2

Background

The prognosis of lung adenocarcinoma (LUAD) is poor. Infection with coronavirus disease 2019 (COVID-19) may further worsen the outcome of LUAD. This study utilized the immune model and the COVID-19 receptor signal to identify the potential immune structure affecting the prognosis of COVID-19 and LUAD.

Methods

A prognostic model was established and verified. The correlation between immune cells and risk score was examined through a variety of immune calculation methods. Gene set variation analysis (GSVA) was used to explore the correlation between the immune signaling pathway, risk model, COVID-19 binding receptor (CO19ORS) signal, and different clinicopathological factors.

Results

The analysis showed that the prognosis of patients was better in the low-risk group versus the high-risk group. The tertiary lymphoid structure dominated by T and B cells (TLS1) can improve the prognosis of patients in the low-risk group. Interestingly, the CO19ORS was enriched only in females and aged >65 years. The age group >65 years is closely related to the tertiary lymphatic structure of the newborn (TLS2), while the female sex is closely related to the TLS2 and TLS1 signature. The two groups exhibited a high level of inflammation-related signal distribution. In the near future, I will collect LUAD and COVID-19 related organizations to verify the changes of 8 risk protein.

Conclusion

TLS1 structure may improve the prognosis of patients with LUAD and SARS-COV-2 (Severe acute respiratory syndrome coronavirus 2). This unexpected discovery provides new insight into the comprehensive treatment of patients with LUAD and SARS-COV-2.

The Role of Pericytes in Regulation of Innate and Adaptive Immunity

Pericytes are perivascular multipotent cells wrapping microvascular capillaries, where they support vasculature functioning, participate in tissue regeneration, and regulate blood flow. However, recent evidence suggests that in addition to traditionally credited structural function, pericytes also manifest immune properties. In this review, we summarise recent data regarding pericytes' response to different pro-inflammatory stimuli and their involvement in innate immune responses through expression of pattern-recognition receptors. Moreover, pericytes express various adhesion molecules, thus regulating trafficking of immune cells across vessel walls. Additionally, the role of pericytes in modulation of adaptive immunity is discussed. Finally, recent reports have suggested that the interaction with cancer cells evokes immunosuppression function in pericytes, thus facilitating immune evasion and facilitating cancer proliferation and metastasis. However, such complex and multi-faceted cross-talks of pericytes with immune cells also suggest a number of potential pericyte-based therapeutic methods and techniques for cancer immunotherapy and treatment of autoimmune and auto-inflammatory disorders.

Comprehensive analysis of T-cell receptor repertoires reveals antigen-driven T-cell clusters in patients with Behçet's syndrome

T lymphocytes are the major components of adaptive immunity in Behçet's syndrome (BS) pathology. However, the precise mechanism of T-cell-induced inflammatory condition remains to be determined. We applied bulk sequencing of the T-cell receptor (TCR) β chain in peripheral blood samples from 45 patients with BS and 10 healthy donors as controls. TCR repertoires in BS patients displayed more clonality and less diversity than in healthy donors. Male patients exhibited lower diversity metrics of TCR and had a larger proportion in the top 10 clones than females (p = 0.016). There were no TCR clonality differences in other clinical features, such as age, disease duration, organ involvement, disease severity, and activity. By "Grouping of Lymphocyte Interactions by Paratope Hotspots" (GLIPH2) for antigen prediction, we found distinct 2477 clusters of TCR-β sequences that potentially recognize similar antigens shared between BS patients. We observed clonal T-cell expansion in BS patients. Sexual differences in TCR clonal expansion and public TCR groups deserve further study to reveal the underline T-cell-mediated immunity in BS.

Modeling and predicting the overlap of B- and T-cell receptor repertoires in healthy and SARS-CoV-2 infected individuals

Adaptive immunity's success relies on the extraordinary diversity of protein receptors on B and T cell membranes. Despite this diversity, the existence of public receptors shared by many individuals gives hope for developing population-wide vaccines and therapeutics. Using probabilistic modeling, we show many of these public receptors are shared by chance in healthy individuals. This predictable overlap is driven not only by biases in the random generation process of receptors, as previously reported, but also by their common functional selection. However, the model underestimates sharing between repertoires of individuals infected with SARS-CoV-2, suggesting strong specific antigen-driven convergent selection. We exploit this discrepancy to identify COVID-associated receptors, which we validate against datasets of receptors with known viral specificity. We study their properties in terms of sequence features and network organization, and use them to design an accurate diagnostic tool for predicting SARS-CoV-2 status from repertoire data.

The use of non-functional clonotypes as a natural calibrator for quantitative bias correction in adaptive immune receptor repertoire profiling

High-throughput sequencing of adaptive immune receptor repertoires is a valuable tool for receiving insights in adaptive immunity studies. Several powerful TCR/BCR repertoire reconstruction and analysis methods have been developed in the past decade. However, detecting and correcting the discrepancy between real and experimentally observed lymphocyte clone frequencies are still challenging. Here, we discovered a hallmark anomaly in the ratio between read count and clone count-based frequencies of non-functional clonotypes in multiplex PCR-based immune repertoires. Calculating this anomaly, we formulated a quantitative measure of V- and J-genes frequency bias driven by multiplex PCR during library preparation called Over Amplification Rate (OAR). Based on the OAR concept, we developed an original software for multiplex PCR-specific bias evaluation and correction named iROAR: immune Repertoire Over Amplification Removal (https://github.com/smiranast/iROAR). The iROAR algorithm was successfully tested on previously published TCR repertoires obtained using both 5' RACE (Rapid Amplification of cDNA Ends)-based and multiplex PCR-based approaches and compared with a biological spike-in-based method for PCR bias evaluation. The developed approach can increase the accuracy and consistency of repertoires reconstructed by different methods making them more applicable for comparative analysis.

Clinical application of immune repertoire sequencing in solid organ transplant

Background

Measurement of T cell receptor (TCR) or B cell receptor (BCR) gene utilization may be valuable in monitoring the dynamic changes in donor-reactive clonal populations following transplantation and enabling adjustment in therapy to avoid the consequences of excess immune suppression or to prevent rejection with contingent graft damage and to indicate the development of tolerance.

Objective

We performed a review of current literature to examine research in immune repertoire sequencing in organ transplantation and to assess the feasibility of this technology for clinical application in immune monitoring.

Methods

We searched MEDLINE and PubMed Central for English-language studies published between 2010 and 2021 that examined T cell/B cell repertoire dynamics upon immune activation. Manual filtering of the search results was performed based on relevancy and predefined inclusion criteria. Data were extracted based on study and methodology characteristics.

Results

Our initial search yielded 1933 articles of which 37 met the inclusion criteria; 16 of these were kidney transplant studies (43%) and 21 were other or general transplantation studies (57%). The predominant method for repertoire characterization was sequencing the CDR3 region of the TCR β chain. Repertoires of transplant recipients were found to have decreased diversity in both rejectors and non-rejectors when compared to healthy controls. Rejectors and those with opportunistic infections were more likely to have clonal expansion in T or B cell populations. Mixed lymphocyte culture followed by TCR sequencing was used in 6 studies to define an alloreactive repertoire and in specialized transplant settings to track tolerance.

Conclusion

Methodological approaches to immune repertoire sequencing are becoming established and offer considerable potential as a novel clinical tool for pre- and post-transplant immune monitoring.

From immune equilibrium to immunodynamics

Objective

The immunology field has long been short of a universally applicable theoretical model that can quantitatively describe the immune response, and the theory of immune equilibrium (balance) is usually limited to the interpretation of the philosophical significance of immune phenomena. Therefore, it is necessary to establish a new immunological theory, namely, immunodynamic theory, to reanalyze the immune response.

Methods

By quantifying the immune dynamic equilibrium as the ratio of positive and negative immune power, the immune dynamic equilibrium equation was established. Then, the area under the curve of the positive and negative immune power was assumed to be equal in the whole process of immune response (regardless of correct or not), and through thought experiments based on this key hypothesis, a series of new concepts and expressions were derived, to establish a series of immunodynamic equations.

Results

New concepts of immune force and immune braking force and their expression equations, namely, the theoretical equations of immunodynamics, were derived through thought experiments, and the theoretical curves of immunodynamics were obtained according to these equations. Via the equivalent transformation of the theoretical equations and practical calculation of functional data, and by the methods of curve comparison and fitting, some practical equations of immunodynamics were established, and these practical equations were used to solve theoretical and practical problems that are related to the immunotherapy of infectious diseases and cancers.

Conclusion

The traditional theory of immune equilibrium has been mathematized and transformed from a philosophical category into a new concrete scientific theory, namely the theory of immunodynamics, which solves the dilemma that the traditional theory cannot guide individualized medical practice for a long time. This new theory may develop into one of the core theories of immunology in the future.

NoisET: Noise Learning and Expansion Detection of T-Cell Receptors

High-throughput sequencing of T- and B-cell receptors makes it possible to track immune repertoires across time, in different tissues, in acute and chronic diseases and in healthy individuals. However, quantitative comparison between repertoires is confounded by variability in the read count of each receptor clonotype due to sampling, library preparation, and expression noise. We review methods for accounting for both biological and experimental noise and present an easy-to-use python package NoisET that implements and generalizes a previously developed Bayesian method. It can be used to learn experimental noise models for repertoire sequencing from replicates, and to detect responding clones following a stimulus. We test the package on different repertoire sequencing technologies and data sets. We review how such approaches have been used to identify responding clonotypes in vaccination and disease data. Availability: NoisET is freely available to use with source code at github.com/statbiophys/NoisET.

B-Cell Receptor Repertoire Sequencing: Deeper Digging into the Mechanisms and Clinical Aspects of Immune-mediated Diseases

B cells play an essential role in adaptive immunity and are intimately correlated with pleiotropic immune-mediated diseases. Each B cell occupies a unique B cell receptor (BCR), and all BCRs throughout our body form “BCR repertoire.” With the development of sequencing technology and coupled bioinformatics, accumulating evidence indicates that BCR repertoire largely varies under physiological and pathological conditions. Therefore, comprehensive grasp of BCR repertoire will provide new insights into the pathogenesis of immune-mediated diseases and help exploit efficient diagnostic and treatment strategies. In this review, we start with an overview of BCR repertoire and related sequencing technologies and summarize their current applications in immune-mediated diseases. We also underscore the challenges of this emerging field and propose promising future directions in advancing BCR repertoire exploration.

Analysis of Gene Expression and TCR/B Cell Receptor Profiling of Immune Cells in Primary Sjögren's Syndrome by Single-Cell Sequencing

Primary Sjögren's syndrome (pSS) is a chronic autoimmune disease that is estimated to affect 35 million people worldwide and is characterized by lymphocytic infiltration, elevated circulating autoantibodies, and proinflammatory cytokines. The key immune cell subset changes and the TCR/BCR repertoire alterations in pSS patients remain unclear. In this study, we sought to comprehensively characterize the transcriptional changes in PBMCs of pSS patients by single-cell RNA sequencing and single-cell V(D)J sequencing. Naive CD8⁺ T cells and mucosal-associated invariant T cells were markedly decreased but regulatory T cells were increased in pSS patients. There were a large number of differentially expressed genes shared by multiple subpopulations of T cells and B cells. Abnormal signaling pathways, including Ag processing and presentation, the BCR signaling pathway, the TCR signaling pathway, and Epstein-Barr virus infection, were highly enriched in pSS patients. Moreover, there were obvious differences in the CD30, FLT3, IFN-II, IL-1, IL-2, IL-6, IL-10, RESISTIN, TGF-β, TNF, and VEGF signaling networks between pSS patients and healthy controls. Single-cell TCR and BCR repertoire analysis showed that there was a lower diversity of T cells in pSS patients than in healthy controls; however, there was no significant difference in the degree of clonal expansion, CDR3 length distribution, or degree of sequence sharing. Notably, our results further emphasize the functional importance of αβ pairing in determining Ag specificity. In conclusion, our analysis provides a comprehensive single-cell map of gene expression and TCR/BCR profiles in pSS patients for a better understanding of the pathogenesis, diagnosis, and treatment of pSS.

Intensive single-cell analysis reveals immune-cell diversity among healthy individuals

Immune responses are different between individuals and personal health histories and unique environmental conditions should collectively determine the present state of immune cells. However, the molecular systems underlying such heterogeneity remain elusive. Here, we conducted a systematic time-lapse single-cell analysis, using 171 single-cell libraries and 30 mass cytometry datasets intensively for seven healthy individuals. We found substantial diversity in immune-cell profiles between different individuals. These patterns showed daily fluctuations even within the same individual. Similar diversities were also observed for the T-cell and B-cell receptor repertoires. Detailed immune-cell profiles at healthy statuses should give essential background information to understand their immune responses, when the individual is exposed to various environmental conditions. To demonstrate this idea, we conducted the similar analysis for the same individuals on the vaccination of influenza and SARS-CoV-2. In fact, we detected distinct responses to vaccines between individuals, although key responses are common. Single-cell immune-cell profile data should make fundamental data resource to understand variable immune responses, which are unique to each individual.

Translatability and transferability of in silico models: Context of use switching to predict the effects of environmental chemicals on the immune system

Immunotoxicity hazard identification of chemicals aims to evaluate the potential for unintended effects of chemical exposure on the immune system. Perfluorinated alkylate substances (PFAS), such as perfluorooctane sulfonate (PFOS) and perfluorooctanoic acid (PFOA), are persistent, globally disseminated environmental contaminants known to be immunotoxic. Elevated PFAS exposure is associated with lower antibody responses to vaccinations in children and in adults. In addition, some studies have reported a correlation between PFAS levels in the body and lower resistance to disease, in other words an increased risk of infections or cancers. In this context, modelling and simulation platforms could be used to simulate the human immune system with the aim to evaluate the adverse effects that immunotoxicants may have. Here, we show the conditions under which a mathematical model developed for one purpose and application (e.g., in the pharmaceutical domain) can be successfully translated and transferred to another (e.g., in the chemicals domain) without undergoing significant adaptation. In particular, we demonstrate that the Universal Immune System Simulator was able to simulate the effects of PFAS on the immune system, introducing entities and new interactions that are biologically involved in the phenomenon. This also revealed a potentially exploitable pathway for assessing immunotoxicity through a computational model.

Side-Directed Release of Differential Extracellular Vesicle-associated microRNA Profiles from Bronchial Epithelial Cells of Healthy and Asthmatic Subjects

Extracellular vesicles (EVs) are released by virtually all cells and may serve as intercellular communication structures by transmitting molecules such as proteins, lipids, and nucleic acids between cells. MicroRNAs (miRNAs) are an abundant class of vesicular RNA playing a pivotal role in regulating intracellular processes. In this work, we aimed to characterize vesicular miRNA profiles released in a side-directed manner by bronchial epithelial cells from healthy and asthmatic subjects using an air−liquid interface cell culture model. EVs were isolated from a culture medium collected from either the basolateral or apical cell side of the epithelial cell cultures and characterized by nano-flow cytometry (NanoFCM) and bead-based flow cytometry. EV-associated RNA profiles were assessed by small RNA sequencing and subsequent bioinformatic analyses. Furthermore, miRNA-associated functions and targets were predicted and miRNA network analyses were performed. EVs were released at higher numbers to the apical cell side of the epithelial cells and were considerably smaller in the apical compared to the basolateral compartment. EVs from both compartments showed a differential tetraspanins surface marker expression. Furthermore, 236 miRNAs were differentially expressed depending on the EV secretion side, regardless of the disease phenotype. On the apical cell side, 32 miRNAs were significantly altered in asthmatic versus healthy conditions, while on the basolateral cell side, 23 differentially expressed miRNAs could be detected. Downstream KEGG pathway analysis predicted mTOR and MAPK signaling pathways as potential downstream targets of apically secreted miRNAs. In contrast, miRNAs specifically detected at the basolateral side were associated with processes of T and B cell receptor signaling. The study proves a compartmentalized packaging of EVs by bronchial epithelial cells supposedly associated with site-specific functions of cargo miRNAs, which are considerably affected by disease conditions such as asthma.

Altered T-Cell Receptor β-Chain and Lactate Dehydrogenase Are Associated With the Immune Pathogenesis of Biliary Atresia

Background: Biliary atresia (BA) is considered to be an autoimmune-mediating inflammatory injury. The pathogenesis of BA has been proposed with the clonal transformation of T cells expressing analogous T-cell receptor β-chain variable regions (TRBVs).

Methods: The TRBV profile of the peripheral blood mononuclear cells (PBMCs) in infants with BA and control infants (healthy donors, HDs), respectively, were characterized by using high-throughput sequencing (HTS). The diversity of T cells was analyzed based on the frequency of complementarity-determining region 3 (CDR3) or V(CDR3)J. Moreover, the correlation between absolute lymphocyte count (ALC) and lactate dehydrogenase (LDH) or diversity (clonality) indices, respectively, were analyzed for subjects with BA and HD.

Results: The diversity indices of CDR3, V(CDR3)J in BA are lower than those in subjects with HD, in addition, there are significantly different levels of neutrophile, neutrophile/lymphocyte ratio (NLR), and LDH between groups of BA and HD. The correlation between ALC and diversity index is significant in subjects with HD but is not for subjects with BA. Conversely, the relationship between ALC and LDH is significant in subjects with BA but is not for subjects with HD. Moreover, 12 CDR3 motifs are deficient or lower expression in BA compared with that in the HD group.

Conclusion: Our results demonstrate that the profile of TRBV repertoire is significantly different between subjects with BA and HD, and suggest that the immune imbalance and elevated LDH level are associated with the pathogenesis of BA. Moreover, the values of neutrophile, NLR, and LDH could be used for the differential diagnosis of BA.

Analysis of TCR Repertoire by High-Throughput Sequencing Indicates the Feature of T Cell Immune Response after SARS-CoV-2 Infection

Coronavirus disease 2019 (COVID-19) is a global infectious disease caused by the SARS-CoV-2 coronavirus. T cells play an essential role in the body’s fighting against the virus invasion, and the T cell receptor (TCR) is crucial in T cell-mediated virus recognition and clearance. However, little has been known about the features of T cell response in convalescent COVID-19 patients. In this study, using 5′RACE technology and PacBio sequencing, we analyzed the TCR repertoire of COVID-19 patients after recovery for 2 weeks and 6 months compared with the healthy donors. The TCR clustering and CDR3 annotation were exploited to discover groups of patient-specific TCR clonotypes with potential SARS-CoV-2 antigen specificities. We first identified CD4⁺ and CD8⁺ T cell clones with certain clonal expansion after infection, and then observed the preferential recombination usage of V(D) J gene segments in CD4⁺ and CD8⁺ T cells of COVID-19 patients with different convalescent stages. More important, the TRBV6-5-TRBD2-TRBJ2-7 combination with high frequency was shared between CD4⁺ T and CD8⁺ T cells of different COVID-19 patients. Finally, we found the dominant characteristic motifs of the CDR3 sequence between recovered COVID-19 and healthy control. Our study provides novel insights on TCR in COVID-19 with different convalescent phases, contributing to our understanding of the immune response induced by SARS-CoV-2.

High-Throughput immunogenetics for precision medicine in cancer

Cancer is characterized by an extremely complex biological background, which hinders personalized therapeutic interventions. Precision medicine promises to overcome this obstacle through integrating information from different 'subsystems', including the host, the external environment, the tumor itself and the tumor micro-environment. Immunogenetics is an essential tool that allows dissecting both lymphoid cancer ontogeny at both a cell-intrinsic and a cell-extrinsic level, i.e. through characterizing micro-environmental interactions, with a view to precision medicine. This is particularly thanks to the introduction of powerful, high-throughput approaches i.e. next generation sequencing, which allow the comprehensive characterization of immune repertoires. Indeed, NGS immunogenetic analysis (Immune-seq) has emerged as key to both understanding cancer pathogenesis and improving the accuracy of clinical decision making in oncology. Immune-seq has applications in lymphoid malignancies, assisting in the diagnosis e.g. through differentiating from reactive conditions, as well as in disease monitoring through accurate assessment of minimal residual disease. Moreover, Immune-seq facilitates the study of T cell receptor clonal dynamics in critical clinical contexts, including transplantation as well as innovative immunotherapy for solid cancers. The clinical utility of Immune-seq represents the focus of the present contribution, where we highlight what can be achieved but also what must be addressed in order to maximally realize the promise of Immune-seq in precision medicine in cancer.

The immune niche of the liver

The liver is an essential organ that is critical for the removal of toxins, the production of proteins, and the maintenance of metabolic homeostasis. Behind each liver functional unit, termed lobules, hides a heterogeneous, complex, and well-orchestrated system. Despite parenchymal cells being most commonly associated with the liver's primary functionality, it has become clear that it is the immune niche of the liver that plays a central role in maintaining both local and systemic homeostasis by propagating hepatic inflammation and orchestrating its resolution. As such, the immunological processes that are at play in healthy and diseased livers are being investigated thoroughly in order to understand the underpinnings of inflammation and the potential avenues for restoring homeostasis. This review highlights recent advances in our understanding of the immune niche of the liver and provides perspectives for how the implementation of new transcriptomic, multimodal, and spatial technologies can uncover the heterogeneity, plasticity, and location of hepatic immune populations. Findings from these technologies will further our understanding of liver biology and create a new framework for the identification of therapeutic targets.

The Value of Flow Cytometry Clonality in Large Granular Lymphocyte Leukemia

Simple Summary

Large granular lymphocyte (LGL) leukemia, a lymphoproliferative disease, is characterized by an increased frequency of large-sized lymphocytes with typical expression of T-cell receptor (TCR) αβ, CD3, CD8, CD16, CD45RA, and CD57, and with the expansion of one to three subfamilies of the TCR variable β chain reflecting gene rearrangements. Molecular analysis remains the gold standard for confirmation of TCR clonality; however, flow cytometry is time and labor saving, and can be associated with simultaneous investigation of other surface markers. Moreover, Vβ usage by flow cytometry can be employed for monitoring clonal kinetics during treatment and follow-up of LGL leukemia patients.

Abstract

Large granular lymphocyte (LGL) leukemia is a lymphoproliferative disorder of mature T or NK cells frequently associated with autoimmune disorders and other hematological conditions, such as myelodysplastic syndromes. Immunophenotype of LGL cells is similar to that of effector memory CD8⁺ T cells with T-cell receptor (TCR) clonality defined by molecular and/or flow cytometric analysis. Vβ usage by flow cytometry can identify clonal TCR rearrangements at the protein level, and is fast, sensitive, and almost always available in every Hematology Center. Moreover, Vβ usage can be associated with immunophenotypic characterization of LGL clone in a multiparametric staining, and clonal kinetics can be easily monitored during treatment and follow-up. Finally, Vβ usage by flow cytometry might identify LGL clones silently underlying other hematological conditions, and routine characterization of Vβ skewing might identify recurrent TCR rearrangements that might trigger aberrant immune responses during hematological or autoimmune conditions.

Profiling T cell receptor β-chain in responders after immunization with recombinant hepatitis B vaccine

BACKGROUND

T cells with edited T cell receptor β-chain variable (TRBV) are involved in the immune response to recombinant hepatitis B surface antigen (rHBsAg) vaccine and the production of hepatitis B surface antibody (HBsAb). The immune repertoire (IR) profile and mechanism of vaccination positive responders (VPR) with rHBsAg are not fully understood.

METHODS

The IR of six VPRs (HBsAb+, HbsAg-) with rHBsAg vaccination was established by the high throughput sequencing technique and bioinformatics analysis and compared with those in five vaccination negative responders (VNRs) (HbsAb-, HbsAg-) who were also inoculated with rHBsAg. The repertoire features of the BV, BJ and V (CDR3) J genes and immune diversity in peripheral blood mononuclear cells, respectively, were analyzed for each subject.

RESULTS

There was no significant difference in sequencing amplification indices of each sample. However, TRBV15/BJ2-3 demonstrated significantly high expression levels in VPR compared to those in the VNR group (both p < 0.05). Further results showed that the BV15/BJ2-5 level was significantly increased for VPR compared to that of VNR group. Interestingly, the motif of CDR3 in TRBV15/BJ2-5 was mostly expressed as "GGETQ" or "GETQ". Additionally, there was no remarkable difference between the two groups of distribution with respect to the different clone expression levels of V (CDR3) J.

CONCLUSIONS

The features of IR in the VPR and VNR will contribute to the exploration of the mechanism of the positive response to rHBsAg, and also contribute to development of optimized hepatitis B vaccine, in addition to providing a partial interpretation of the VNR who has a relatively low infection with HBV.

Biological controls for standardization and interpretation of adaptive immune receptor repertoire profiling

Use of adaptive immune receptor repertoire sequencing (AIRR-seq) has become widespread, providing new insights into the immune system with potential broad clinical and diagnostic applications. However, like many high-throughput technologies, it comes with several problems, and the AIRR Community was established to understand and help solve them. We, the AIRR Community's Biological Resources Working Group, have surveyed scientists about the need for standards and controls in generating and annotating AIRR-seq data. Here, we review the current status of AIRR-seq, provide the results of our survey, and based on them, offer recommendations for developing AIRR-seq standards and controls, including future work.

Generative models of T-cell receptor sequences

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

Toll-Like Receptor 9 Agonists in Cancer

Toll-like receptor 9 (TLR9) is a pattern recognition receptor that is predominantly located intracellularly in immune cells, including dendritic cells, macrophages, natural killer cells, and other antigen-presenting cells (APC). The primary ligands for TLR9 receptors are unmethylated cytidine phosphate guanosine (CpG) oligodinucleotides (ODN). TLR9 agonists induce inflammatory processes that result in the enhanced uptake and killing of microorganisms and cancer cells as well as the generation of adaptive immune responses. Preclinical studies of TLR9 agonists suggested efficacy both as monotherapy and in combination with several agents, which led to clinical trials in patients with advanced cancer. In these studies, intravenous, intratumoral, and subcutaneous routes of administration have been tested; with anti-tumor responses in both treated and untreated metastatic sites. TLR9 agonist monotherapy is safe, although efficacy is minimal in advanced cancer patients; conversely, combinations appear to be more promising. Several ongoing phase I and II clinical trials are evaluating TLR9 agonists in combination with a variety of agents including chemotherapy, radiotherapy, targeted therapy, and immunotherapy agents. In this review article, we describe the distribution, structure and signaling of TLR9; discuss the results of preclinical studies of TLR9 agonists; and review ongoing clinical trials of TLR9 agonists singly and in combination in patients with advanced solid tumors.

Comprehensive investigation of T and B cell receptor repertoires in an MC38 tumor model following murine anti‑PD‑1 administration

The MC38 (derived from carcinogen-induced colon adenocarcinoma) tumor model is sensitive to anti-programmed cell death-1 (anti-PD-1) treatment. However, there is no comprehensive description of the T and B cell receptor (TCR, BCR) repertoires of the MC38 tumor model following anti-PD-1 treatment, an improved understanding of which is highly important in the development of anti-PD-1 immunotherapy. The present study analyzed the TCR and BCR repertoires of three types of tissue, including tumor, spleen and tumor draining lymph node (DLN) from 20 MC38 syngeneic mice receiving murine anti-PD-1 (mDX400) treatment or mouse immunoglobulin G1 (mIgG1) control treatment. To obtain enough tissues for high-throughput sequencing, samples were collected on day 8 after the start of initial treatment. The usage frequencies of seven TCR β chain (TRB) V genes and one TRBJ gene were significantly different between mDX400- and mIgG1-group tumors. TCR repertoire diversity was significantly lower in mDX400-group tumors compared with mIgG1-group tumors, with the top 10 most frequent TCR clonotypes notably expanded in mDX400-group tumors. In addition, the proportion of high-frequency TCR clonotypes from mDX400-group tumors that were also present both in the DLN and spleen was significantly higher than that in mIgG1-group tumors. Among the highly expanded TCR clonotypes, one TCR clonotype was consistently expanded in >50% of the mDX400-group tumors compared with mIgG1-group tumors. Similarly, one BCR clonal family was highly expanded in >50% of mDX400-group tumor samples. The consistently expanded TCR and BCR clones were co-expanded in 29% of mDX400-group tumors. Moreover, mutation rates of immunoglobulin heavy chain sequences in the spleen within complementarity determining region 2 and framework region 3 were significantly higher in the mDX400 group than in the mIgG1 group. The findings of this study may contribute to an improved understanding of the molecular mechanisms of anti-PD-1 treatment.