The mechanism and regulation of chromosomal V(D)J recombination

Differential tumor immune microenvironment coupled with tumor progression or tumor eradication in HPV-antigen expressing squamous cell carcinoma (SCC) models

Human papilloma virus (HPV) is an etiological factor of head and neck squamous cell carcinoma (HNSCC). To investigate the role of HPV antigen in anti-tumor immunity, we established mouse models by expressing HPV16 E6 and E7 in a SCC tumor cell line. We obtained two HPV antigen-expressing clones (C-225 and C-100) transplantable into C57BL/6 recipients. We found that C-225 elicited complete eradication in C57BL/6 mice (eradicated), whereas C-100 grew progressively (growing). We examined immune tumor microenvironment (TME) using flow cytometry and found that eradicated or growing tumors exhibited differential immune profiles that may influence the outcome of anti-tumor immunity. Surprisingly, the percentage of CD8 and CD4 tumor-infiltrating lymphocytes (TILs) was much higher in growing (C-100) than eradicated (C-225) tumor. However, the TILs upregulated PD-1 and LAG-3 more potently and exhibited impaired effector functions in growing tumor compared to their counterparts in eradicated tumor. C-225 TME is highly enriched with myeloid cells, especially polymorphonuclear (PMN) myeloid-derived suppressor cells (MDSC), whereas the percentage of M-MDSC and tumor-associated macrophages (TAMs) was much higher in C-100 TME, especially M2-TAMs (CD206+). The complete eradication of C-225 depended on CD8 T cells and elicited anti-tumor memory responses upon secondary tumor challenge. We employed DNA sequencing to identify differences in the T cell receptor of peripheral blood lymphocytes pre- and post-secondary tumor challenge. Lastly, C-225 and C-100 tumor lines harbored different somatic mutations. Overall, we uncovered differential immune TME that may underlie the divergent outcomes of anti-tumor immunity by establishing two SCC tumor lines, both of which express HPV16 E6 and E7 antigens. Our experimental models may provide a platform for pinpointing tumor-intrinsic versus host-intrinsic differences in orchestrating an immunosuppressive TME in HNSCCs and for identifying new targets that render tumor cells vulnerable to immune attack.

Protocol for fast clonal family inference and analysis from large-scale B cell receptor repertoire sequencing data

The expeditious identification and comprehensive analysis of clonal families from extensive B cell receptor (BCR) repertoire sequencing data are imperative for elucidating the intricacies of B cell immune responses. Here, we introduce a computational pipeline designed to swiftly deduce clonal families from bulk BCR heavy-chain sequencing data, accompanied by a suite of functional modules tailored to streamline post-clustering analysis. The outlined methodology encompasses guidelines for software installation, meticulous data preparation, and the systematic inference and analysis of clonal families. For complete details on the use and execution of this protocol, please refer to Wang et al.1.

Immunologic tumor microenvironment modulators for turning cold tumors hot

Tumors can be classified into distinct immunophenotypes based on the presence and arrangement of cytotoxic immune cells within the tumor microenvironment (TME). Hot tumors, characterized by heightened immune activity and responsiveness to immune checkpoint inhibitors (ICIs), stand in stark contrast to cold tumors, which lack immune infiltration and remain resistant to therapy. To overcome immune evasion mechanisms employed by tumor cells, novel immunologic modulators have emerged, particularly ICIs targeting cytotoxic T-lymphocyte-associated protein 4 (CTLA-4) and programmed cell death protein 1/programmed death-ligand 1(PD-1/PD-L1). These agents disrupt inhibitory signals and reactivate the immune system, transforming cold tumors into hot ones and promoting effective antitumor responses. However, challenges persist, including primary resistance to immunotherapy, autoimmune side effects, and tumor response heterogeneity. Addressing these challenges requires innovative strategies, deeper mechanistic insights, and a combination of immune interventions to enhance the effectiveness of immunotherapies. In the landscape of cancer medicine, where immune cold tumors represent a formidable hurdle, understanding the TME and harnessing its potential to reprogram the immune response is paramount. This review sheds light on current advancements and future directions in the quest for more effective and safer cancer treatment strategies, offering hope for patients with immune-resistant tumors.

Principles and therapeutic applications of adaptive immunity

Adaptive immunity provides protection against infectious and malignant diseases. These effects are mediated by lymphocytes that sense and respond with targeted precision to perturbations induced by pathogens and tissue damage. Here, we review key principles underlying adaptive immunity orchestrated by distinct T cell and B cell populations and their extensions to disease therapies. We discuss the intracellular and intercellular processes shaping antigen specificity and recognition in immune activation and lymphocyte functions in mediating effector and memory responses. We also describe how lymphocytes balance protective immunity against autoimmunity and immunopathology, including during immune tolerance, response to chronic antigen stimulation, and adaptation to non-lymphoid tissues in coordinating tissue immunity and homeostasis. Finally, we discuss extracellular signals and cell-intrinsic programs underpinning adaptive immunity and conclude by summarizing key advances in vaccination and engineering adaptive immune responses for therapeutic interventions. A deeper understanding of these principles holds promise for uncovering new means to improve human health.

Exploring the Potential of Structure-Based Deep Learning Approaches for T cell Receptor Design

Deep learning methods, trained on the increasing set of available protein 3D structures and sequences, have substantially impacted the protein modeling and design field. These advancements have facilitated the creation of novel proteins, or the optimization of existing ones designed for specific functions, such as binding a target protein. Despite the demonstrated potential of such approaches in designing general protein binders, their application in designing immunotherapeutics remains relatively unexplored. A relevant application is the design of T cell receptors (TCRs). Given the crucial role of T cells in mediating immune responses, redirecting these cells to tumor or infected target cells through the engineering of TCRs has shown promising results in treating diseases, especially cancer. However, the computational design of TCR interactions presents challenges for current physics-based methods, particularly due to the unique natural characteristics of these interfaces, such as low affinity and cross-reactivity. For this reason, in this study, we explored the potential of two structure-based deep learning protein design methods, ProteinMPNN and ESM-IF, in designing fixed-backbone TCRs for binding target antigenic peptides presented by the MHC through different design scenarios. To evaluate TCR designs, we employed a comprehensive set of sequence- and structure-based metrics, highlighting the benefits of these methods in comparison to classical physics-based design methods and identifying deficiencies for improvement.

When DNA-damage responses meet innate and adaptive immunity

When cells proliferate, stress on DNA replication or exposure to endogenous or external insults frequently results in DNA damage. DNA-Damage Response (DDR) networks are complex signaling pathways used by multicellular organisms to prevent DNA damage. Depending on the type of broken DNA, the various pathways, Base-Excision Repair (BER), Nucleotide Excision Repair (NER), Mismatch Repair (MMR), Homologous Recombination (HR), Non-Homologous End-Joining (NHEJ), Interstrand Crosslink (ICL) repair, and other direct repair pathways, can be activated separately or in combination to repair DNA damage. To preserve homeostasis, innate and adaptive immune responses are effective defenses against endogenous mutation or invasion by external pathogens. It is interesting to note that new research keeps showing how closely DDR components and the immune system are related. DDR and immunological response are linked by immune effectors such as the cyclic GMP-AMP synthase (cGAS)-Stimulator of Interferon Genes (STING) pathway. These effectors act as sensors of DNA damage-caused immune response. Furthermore, DDR components themselves function in immune responses to trigger the generation of inflammatory cytokines in a cascade or even trigger programmed cell death. Defective DDR components are known to disrupt genomic stability and compromise immunological responses, aggravating immune imbalance and leading to serious diseases such as cancer and autoimmune disorders. This study examines the most recent developments in the interaction between DDR elements and immunological responses. The DDR network's immune modulators' dual roles may offer new perspectives on treating infectious disorders linked to DNA damage, including cancer, and on the development of target immunotherapy.

Experimental Models to Study Immune Dysfunction in the Pathogenesis of Parkinson's Disease

Parkinson's disease (PD) is a chronic, age-related, progressive multisystem disease associated with neuroinflammation and immune dysfunction. This review discusses the methodological approaches used to study the changes in central and peripheral immunity in PD, the advantages and limitations of the techniques, and their applicability to humans. Although a single animal model cannot replicate all pathological features of the human disease, neuroinflammation is present in most animal models of PD and plays a critical role in understanding the involvement of the immune system (IS) in the pathogenesis of PD. The IS and its interactions with different cell types in the central nervous system (CNS) play an important role in the pathogenesis of PD. Even though culture models do not fully reflect the complexity of disease progression, they are limited in their ability to mimic long-term effects and need validation through in vivo studies. They are an indispensable tool for understanding the interplay between the IS and the pathogenesis of this disease. Understanding the immune-mediated mechanisms may lead to potential therapeutic targets for the treatment of PD. We believe that the development of methodological guidelines for experiments with animal models and PD patients is crucial to ensure the validity and consistency of the results.

Against all odds: The road to success in the development of human immune reconstitution mice

The mouse genome has a high degree of homology with the human genome, and its physiological, biochemical, and developmental regulation mechanisms are similar to those of humans; therefore, mice are widely used as experimental animals. However, it is undeniable that interspecies differences between humans and mice can lead to experimental errors. The differences in the immune system have become an important factor limiting current immunological research. The application of immunodeficient mice provides a possible solution to these problems. By transplanting human immune cells or tissues, such as peripheral blood mononuclear cells or hematopoietic stem cells, into immunodeficient mice, a human immune system can be reconstituted in the mouse body, and the engrafted immune cells can elicit human-specific immune responses. Researchers have been actively exploring the development and differentiation conditions of host recipient animals and grafts in order to achieve better immune reconstitution. Through genetic engineering methods, immunodeficient mice can be further modified to provide a favorable developmental and differentiation microenvironment for the grafts. From initially only being able to reconstruct single T lymphocyte lineages, it is now possible to reconstruct lymphoid and myeloid cells, providing important research tools for immunology-related studies. In this review, we compare the differences in immune systems of humans and mice, describe the development history of human immune reconstitution from the perspectives of immunodeficient mice and grafts, and discuss the latest advances in enhancing the efficiency of human immune cell reconstitution, aiming to provide important references for immunological related researches.

DNA-dependent protein kinase regulates cytosolic double-stranded DNA secretion from irradiated macrophages to increase radiosensitivity of tumors

BACKGROUND AND PURPOSE

To investigate the molecular mechanism by which irradiated macrophages secrete cytosolic double-stranded DNA (c-dsDNA) to increase radiosensitivity of tumors.

MATERIALS AND METHODS

Irradiated bone marrow-derived macrophages (BMDM) were co-incubated with irradiated EO771 or MC38 cancer cells to determine clonogenic survival. c-dsDNA were measured by agarose gel or enzyme-linked immunosorbent assay. BMDM or cancer cells were analyzed with immunostaining or western blot. Subcutaneously implanted MC38 cells in myeloid-specific Prkdc knockout (KO) mice or littermate control mice were irradiated with 8 Gy to determine radiosensitivity of tumors.

RESULTS

We observed that irradiated BMDM significantly increased radiosensitivity of cancer cells. By performing immunostaining, we found that there was a dose-dependent increase in the formation of c-dsDNA and phosphorylation in DNA-dependent protein kinase (DNA-PK) in irradiated BMDM. Importantly, c-dsDNA in irradiated BMDM could be secreted to the extracellular milieu and this process required DNA-PK, which phosphorylated myosin light chain to regulate the secretion. The secreted c-dsDNA from irradiated BMDM then activated toll-like receptor-9 and subsequent nuclear factor kappa-light-chain-enhancer of activated B cells signaling in the adjacent cancer cells inhibiting radiation-induced DNA double strand break repair. Lastly, we observed that irradiated tumors in vivo had a significantly increased number of tumor-associated macrophages (TAM) with phosphorylated DNA-PK expression in the cytosol. Furthermore, tumors grown in myeloid-specific Prkdc KO mice, in which TAM lacked phosphorylated DNA-PK expression were significantly more radioresistant than those of the wild-type control mice.

CONCLUSIONS

Irradiated macrophages can increase antitumor efficacy of radiotherapy through secretion of c-dsDNA under the regulation of DNA-PK.

The Cyclin D3 Protein Enforces Monogenic TCRβ Expression by Mediating TCRβ Protein-Signaled Feedback Inhibition of Vβ Recombination

In jawed vertebrates, adaptive immunity depends on the process of V(D)J recombination creating vast numbers of T and B lymphocytes that each expresses unique Ag receptors of uniform specificity. The asynchronous initiation of V-to-(D)J rearrangement between alleles and the resulting protein from one allele signaling feedback inhibition of V recombination on the other allele ensures homogeneous receptor specificity of individual cells. Upon productive Vβ-to-DβJβ rearrangements in noncycling double-negative thymocytes, TCRβ protein signals induction of the cyclin D3 protein to accelerate cell cycle entry, thereby driving proliferative expansion of developing αβ T cells. Through undetermined mechanisms, the inactivation of cyclin D3 in mice causes an increased frequency of αβ T cells that express TCRβ proteins from both alleles, producing lymphocytes of heterogeneous specificities. To determine how cyclin D3 enforces monogenic TCRβ expression, we used our mouse lines with enhanced rearrangement of specific Vβ segments due to replacement of their poor-quality recombination signal sequence (RSS) DNA elements with a better RSS. We show that cyclin D3 inactivation in these mice elevates the frequencies of αβ T cells that display proteins from RSS-augmented Vβ segments on both alleles. By assaying mature αβ T cells, we find that cyclin D3 deficiency increases the levels of Vβ rearrangements that occur within developing thymocytes. Our data demonstrate that a component of the cell cycle machinery mediates TCRβ protein-signaled feedback inhibition in thymocytes to achieve monogenic TCRβ expression and resulting uniform specificity of individual αβ T cells.

Inflammation in Development and Aging: Insights from the Zebrafish Model

Zebrafish are an emergent animal model to study human diseases due to their significant genetic similarity to humans, swift development, and genetic manipulability. Their utility extends to the exploration of the involvement of inflammation in host defense, immune responses, and tissue regeneration. Additionally, the zebrafish model system facilitates prompt screening of chemical compounds that affect inflammation. This study explored the diverse roles of inflammatory pathways in zebrafish development and aging. Serving as a crucial model, zebrafish provides insights into the intricate interplay of inflammation in both developmental and aging contexts. The evidence presented suggests that the same inflammatory signaling pathways often play instructive or beneficial roles during embryogenesis and are associated with malignancies in adults.

Locus folding mechanisms determine modes of antigen receptor gene assembly

The dynamic folding of genomes regulates numerous biological processes, including antigen receptor (AgR) gene assembly. We show that, unlike other AgR loci, homotypic chromatin interactions and bidirectional chromosome looping both contribute to structuring Tcrb for efficient long-range V(D)J recombination. Inactivation of the CTCF binding element (CBE) or promoter at the most 5'Vβ segment (Trbv1) impaired loop extrusion originating locally and extending to DβJβ CBEs at the opposite end of Tcrb. Promoter or CBE mutation nearly eliminated Trbv1 contacts and decreased RAG endonuclease-mediated Trbv1 recombination. Importantly, Trbv1 rearrangement can proceed independent of substrate orientation, ruling out scanning by DβJβ-bound RAG as the sole mechanism of Vβ recombination, distinguishing it from Igh. Our data indicate that CBE-dependent generation of loops cooperates with promoter-mediated activation of chromatin to juxtapose Vβ and DβJβ segments for recombination through diffusion-based synapsis. Thus, the mechanisms that fold a genomic region can influence molecular processes occurring in that space, which may include recombination, repair, and transcriptional programming.

Dclre1c-Mutation-Induced Immunocompromised Mice Are a Novel Model for Human Xenograft Research

Severe combined immunodeficient (SCID) mice serve as a critical model for human xenotransplantation studies, yet they often suffer from low engraftment rates and susceptibility to graft-versus-host disease (GVHD). Moreover, certain SCID strains demonstrate 'immune leakage', underscoring the need for novel model development. Here, we introduce an SCID mouse model with a targeted disruption of the dclre1c gene, encoding Artemis, which is essential for V(D)J recombination and DNA repair during T cell receptor (TCR) and B cell receptor (BCR) assembly. Artemis deficiency precipitates a profound immunodeficiency syndrome, marked by radiosensitivity and compromised T and B lymphocyte functionality. Utilizing CRISPR/Cas9-mediated gene editing, we generated dclre1c-deficient mice with an NOD genetic background. These mice exhibited a radiosensitive SCID phenotype, with pronounced DNA damage and defective thymic, splenic and lymph node development, culminating in reduced T and B lymphocyte populations. Notably, both cell lines and patient-derived tumor xenografts were successfully engrafted into these mice. Furthermore, the human immune system was effectively rebuilt following peripheral blood mononuclear cells (PBMCs) transplantation. The dclre1c-knockout NOD mice described herein represent a promising addition to the armamentarium of models for xenotransplantation, offering a valuable platform for advancing human immunobiological research.

The Structure Characteristics and Function of Non B Cell-Derived Immunoglobulin

According to classical immunology theory, immunoglobulin (Ig) is exclusively produced by differentiated B lymphocytes, which exhibit a typical tetrapeptide chain structure and are predominantly present on the surface of B cells and in bodily fluids. B-Ig is one of the critical effector molecules for humoral immune responses specifically recognising antigens and eliminating them. However, mounting evidence has demonstrated that Ig is widely expressed in non B lineage cells, especially malignant ones (referred to as non B-Ig). Interestingly, non B-Ig mainly resides in the cytoplasm and secretion, but to some extent on the cell surface. Furthermore non B-Ig not only displays a tetrapeptide chain structure but also shows free heavy chains and free light chains (FLCs). Additionally, Ig derived from non B cancer cell typically displays unique glycosylation modifications. Functionally, non B-Ig demonstrated diversity and versatility, showing antibody activity and cellular biological activity, such as promoting cell proliferation and survival, and it is implicated in cancer progression and some immune-related diseases, such as renal diseases.

Structural analysis of light chain-driven bispecific antibodies targeting CD47 and PD-L1

In contrast to natural antibodies that rely mainly on the heavy chain to establish contacts with their cognate antigen, we have developed a bispecific antibody format in which the light chain (LC) drives antigen binding and specificity. To better understand epitope-paratope interactions in this context, we determined the X-ray crystallographic structures of an antigen binding fragment (Fab) in complex with human CD47 and another Fab in complex with human PD-L1. These Fabs contain a κ-LC and a λ-LC, respectively, which are paired with an identical heavy chain (HC). The structural analysis of these complexes revealed the dominant contribution of the LCs to antigen binding, but also that the common HC provides some contacts in both CD47 and PD-L1 Fab complexes. The anti-CD47 Fab was affinity optimized by diversifying complementary-determining regions of the LC followed by phage display selections. Using homology modeling, the contributions of the amino acid modification to the affinity increase were analyzed. Our results demonstrate that, despite a less prominent role in natural antibodies, the LC can mediate high affinity binding to different antigens and neutralize their biological function. Importantly, Fabs containing a common variable heavy (VH) domain enable the generation of bispecific antibodies retaining a truly native structure, maximizing their therapeutic potential.

TEPCAM: Prediction of T-cell receptor-epitope binding specificity via interpretable deep learning

The recognition of T-cell receptor (TCR) on the surface of T cell to specific epitope presented by the major histocompatibility complex is the key to trigger the immune response. Identifying the binding rules of TCR-epitope pair is crucial for developing immunotherapies, including neoantigen vaccine and drugs. Accurate prediction of TCR-epitope binding specificity via deep learning remains challenging, especially in test cases which are unseen in the training set. Here, we propose TEPCAM (TCR-EPitope identification based on Cross-Attention and Multi-channel convolution), a deep learning model that incorporates self-attention, cross-attention mechanism, and multi-channel convolution to improve the generalizability and enhance the model interpretability. Experimental results demonstrate that our model outperformed several state-of-the-art models on two challenging tasks including a strictly split dataset and an external dataset. Furthermore, the model can learn some interaction patterns between TCR and epitope by extracting the interpretable matrix from cross-attention layer and mapping them to the three-dimensional structures. The source code and data are freely available at https://github.com/Chenjw99/TEPCAM.

Somatic hypermutation introduces bystander mutations that prepare SARS-CoV-2 antibodies for emerging variants

Somatic hypermutation (SHM) drives affinity maturation and continues over months in SARS-CoV-2-neutralizing antibodies (nAbs). However, several potent SARS-CoV-2 antibodies carry no or only a few mutations, leaving the question of how ongoing SHM affects neutralization unclear. Here, we reverted variable region mutations of 92 antibodies and tested their impact on SARS-CoV-2 binding and neutralization. Reverting higher numbers of mutations correlated with decreasing antibody functionality. However, for some antibodies, including antibodies of the public clonotype VH1-58, neutralization of Wu01 remained unaffected. Although mutations were dispensable for Wu01-induced VH1-58 antibodies to neutralize Alpha, Beta, and Delta variants, they were critical for Omicron BA.1/BA.2 neutralization. We exploited this knowledge to convert the clinical antibody tixagevimab into a BA.1/BA.2 neutralizer. These findings broaden our understanding of SHM as a mechanism that not only improves antibody responses during affinity maturation but also contributes to antibody diversification, thus increasing the chances of neutralizing viral escape variants.

Characterizing the dynamics of BCR repertoire from repeated influenza vaccination

Yearly epidemics of seasonal influenza cause an enormous disease burden around the globe. An understanding of the rules behind the immune response with repeated vaccination still presents a significant challenge, which would be helpful for optimizing the vaccination strategy. In this study, 34 healthy volunteers with 16 vaccinated were recruited, and the dynamics of the BCR repertoire for consecutive vaccinations in two seasons were tracked. In terms of diversity, length, network, V and J gene segments usage, somatic hypermutation (SHM) rate and isotype, it was found that the overall changes were stronger in the acute phase of the first vaccination than the second vaccination. However, the V gene segments of IGHV4-39, IGHV3-9, IGHV3-7 and IGHV1-69 were amplified in the acute phase of the first vaccination, with IGHV3-7 dominant. On the other hand, for the second vaccination, the changes were dominated by IGHV1-69, with potential for coding broad neutralizing antibody. Additional analysis indicates that the application of V gene segment for IGHV3-7 in the acute phase of the first vaccination was due to the elevated usage of isotypes IgM and IgG3. While for IGHV1-69 in the second vaccination, it was contributed by isotypes IgG1 and IgG2. Finally, 41 public BCR clusters were identified in the vaccine group, with both IGHV3-7 and IGHV1-69 were involved and representative complementarity determining region 3 (CDR3) motifs were characterized. This study provides insights into the immune response dynamics following repeated influenza vaccination in humans and can inform universal vaccine design and vaccine strategies in the future.

Specific T-cell receptor beta-rearrangements of gluten-triggered CD8+ T-cells are enriched in celiac disease patients' duodenal mucosa

Celiac disease (CeD) is an autoimmune disorder affecting the small intestine with gluten as disease trigger. Infections including Influenza A, increase the CeD risk. While gluten-specific CD4+ T-cells, recognizing HLA-DQ2/DQ8 presented gluten-peptides, initiate and sustain the celiac immune response, CD8+ α/β intraepithelial T-cells elicit mucosal damage. Here, we subjected TCRs from a cohort of 56 CeD patients and 22 controls to an analysis employing 749 published CeD-related TCRβ-rearrangements derived from gluten-specific CD4+ T-cells and gluten-triggered peripheral blood CD8+ T-cells. We show, that in addition to TCRs from gluten-specific CD4+ T-cells, TCRs of gluten-triggered CD8+ T-cells are significantly enriched in CeD duodenal tissue samples. TCRβ-rearrangements of gluten-triggered CD8+ T-cells were even more expanded in patients than TCRs from gluten-specific CD4+ T-cells (p < 0.0002) and highest in refractory CeD. Sequence alignments with TCR-antigen databases suggest that a subgroup of these most likely indirectly gluten-triggered TCRs recognize microbial, viral, and autoantigens.

Fast clonal family inference from large-scale B cell repertoire sequencing data

Advances in high-throughput sequencing technologies have facilitated the large-scale characterization of B cell receptor (BCR) repertoires. However, the vast amount and high diversity of the BCR sequences pose challenges for efficient and biologically meaningful analysis. Here, we introduce fastBCR, an efficient computational approach for inferring B cell clonal families from massive BCR heavy chain sequences. We demonstrate that fastBCR substantially reduces the running time while ensuring high accuracy on simulated datasets with diverse numbers of B cell lineages and varying mutation rates. We apply fastBCR to real BCR sequencing data from peripheral blood samples of COVID-19 patients, showing that the inferred clonal families display disease-associated features, as well as corresponding antigen-binding specificity and affinity. Overall, our results demonstrate the advantages of fastBCR for analyzing BCR repertoire data, which will facilitate the identification of disease-associated antibodies and improve our understanding of the B cell immune response.

Rag1 Deficiency Impairs Arteriogenesis in Mice

Increasing evidence suggests that lymphocytes play distinct roles in inflammation-induced tissue remodeling and tissue damage. Arteriogenesis describes the growth of natural bypasses from pre-existing collateral arteries. This process compensates for the loss of artery function in occlusive arterial diseases. The role of innate immune cells is widely understood in the process of arteriogenesis, whereas the role of lymphocytes remains unclear and is the subject of the present study. To analyze the role of lymphocytes, we induced arteriogenesis in recombination activating gene-1 (Rag1) knockout (KO) mice by unilateral ligation of the femoral artery. The lack of functional lymphocytes in Rag1 KO mice resulted in reduced perfusion recovery as shown by laser Doppler imaging. Additionally, immunofluorescence staining revealed a reduced vascular cell proliferation along with a smaller inner luminal diameter in Rag1 KO mice. The perivascular macrophage polarization around the growing collateral arteries was shifted to more pro-inflammatory M1-like polarized macrophages. Together, these data suggest that lymphocytes are crucial for arteriogenesis by modulating perivascular macrophage polarization.

Differentiation and traffic of IgM+ B cells between focal dark spots in skeletal muscle of Atlantic salmon, lymphoid and adipose tissues

Focal dark spots (DS) in farmed Atlantic salmon fillets contain a significant number of B cells as revealed by the high abundance of immunoglobulin (Ig) transcripts in transcriptome data. The immune response in DS remains unknown while they represent a major problem in commercial aquaculture. Here, we characterized the diversity and clonal composition of B cells in DS. Sixteen gene markers of immune cells and antigen presentation were analyzed with RT-qPCR. All genes expression showed a positive correlation with DS area and intensity. The flatter the DS, the higher the expression of cd28, csfr, ctla, igt, and sigm, the lower expression of cd83 and btla, and the larger the cumulative frequency within DS. The expression of most of the analyzed immune genes, including three Ig types and markers of B cells was lower in DS than in the lymphatic organs, head kidney and spleen, but significantly higher compared to skeletal muscle. High levels of ctla4 and cd28 in DS might indicate the recruitment of T cells. Sequencing of IgM repertoire (Ig-seq) assessed migration of B cells by co-occurrence of identical CDR3 sequences in different tissues. The combination of gene expression and Ig-seq revealed the presence of several stages of B cell differentiation in DS. B cells at the earliest stage, with high ratio of membrane to secretory IgM (migm and sigm), showed minor Ig repertoire overlap with other tissues. Further differentiation stage (increased sigm to migm ratio and high expression of pax5 and cd79) was associated with active movement of B cells from DS towards lymphatic organs and visceral fat. Traffic and expression of immune genes decreased at later stages. These B cells could be involved in a response directed against viruses, pathogenic or opportunistic bacteria in DS. Seven of eight fish were positive for salmon alphavirus, and levels were higher in DS than in unstained muscle. PCR with universal primers to the 16S rRNA gene did not detect bacteria in DS. Although the evolution of DS most likely implies local exposure to antigens, neither this nor previous studies have found a necessary association between DS and pathogens or self-antigens.

Immune response of DNA vaccinated-gilthead seabream (Sparus aurata) against LCDV-Sa infection: relevance of the inflammatory process

Lymphocystis disease is one of the main viral pathologies affecting cultured gilthead seabream (Sparus aurata) in the Mediterranean region. Recently, we have developed a DNA vaccine based on the major capsid protein (MCP) of the Lymphocystis disease virus 3 (LCDV-Sa). The immune response triggered by either LCDV-Sa infection or vaccination have been previously studied and seem to be highly related to the modulation of the inflammatory and the IFN response. However, a comprehensive evaluation of immune-related gene expression in vaccinated fish after viral infection to identify immunogenes involved in vaccine-induced protection have not been carried out to date. The present study aimed to fulfill this objective by analyzing samples of head-kidney, spleen, intestine, and caudal fin from fish using an OpenArray® platform containing targets related to the immune response of gilthead seabream. The results obtained showed an increase of deregulated genes in the hematopoietic organs between vaccinated and non-vaccinated fish. However, in the intestine and fin, the results showed the opposite trend. The global effect of fish vaccination was a significant decrease (p<0.05) of viral replication in groups of fish previously vaccinated, and the expression of the following immune genes related to viral recognition (tlr9), humoral and cellular response (rag1 and cd48), inflammation (csf1r, elam, il1β, and il6), antiviral response (isg15, mx1, mx2, mx3), cell-mediated cytotoxicity (nccrp1), and apoptosis (prf1). The exclusive modulation of the immune response provoked by the vaccination seems to control the progression of the infection in the experimentally challenged gilthead seabream.

Large-scale antibody immune response mapping of splenic B cells and bone marrow plasma cells in a transgenic mouse model

Molecular characterization of antibody immunity and human antibody discovery is mainly carried out using peripheral memory B cells, and occasionally plasmablasts, that express B cell receptors (BCRs) on their cell surface. Despite the importance of plasma cells (PCs) as the dominant source of circulating antibodies in serum, PCs are rarely utilized because they do not express surface BCRs and cannot be analyzed using antigen-based fluorescence-activated cell sorting. Here, we studied the antibodies encoded by the entire mature B cell populations, including PCs, and compared the antibody repertoires of bone marrow and spleen compartments elicited by immunization in a human immunoglobulin transgenic mouse strain. To circumvent prior technical limitations for analysis of plasma cells, we applied single-cell antibody heavy and light chain gene capture from the entire mature B cell repertoires followed by yeast display functional analysis using a cytokine as a model immunogen. We performed affinity-based sorting of antibody yeast display libraries and large-scale next-generation sequencing analyses to follow antibody lineage performance, with experimental validation of 76 monoclonal antibodies against the cytokine antigen that identified three antibodies with exquisite double-digit picomolar binding affinity. We observed that spleen B cell populations generated higher affinity antibodies compared to bone marrow PCs and that antigen-specific splenic B cells had higher average levels of somatic hypermutation. A degree of clonal overlap was also observed between bone marrow and spleen antibody repertoires, indicating common origins of certain clones across lymphoid compartments. These data demonstrate a new capacity to functionally analyze antigen-specific B cell populations of different lymphoid organs, including PCs, for high-affinity antibody discovery and detailed fundamental studies of antibody immunity.

TPBTE: A model based on convolutional Transformer for predicting the binding of TCR to epitope

T cell receptors (TCRs) selectively bind to antigens to fight pathogens with specific immunity. Current tools focus on the nature of amino acids within sequences and take less into account the nature of amino acids far apart and the relationship between sequences, leading to significant differences in the results from different datasets. We propose TPBTE, a model based on convolutional Transformer for Predicting the Binding of TCR to Epitope. It takes epitope sequences and the complementary decision region 3 (CDR3) sequences of TCRβ chain as inputs. And it uses a convolutional attention mechanism to learn amino acid representations between different positions of the sequences based on learning local features of the sequences. At the same time, it uses cross attention to learn the interaction information between TCR sequences and epitope sequences. A comprehensive evaluation of the TCR-epitope data shows that the average area under the curve of TPBTE outperforms the baseline model, and demonstrate an intentional performance. In addition, TPBTE can give the probability of binding TCR to epitopes, which can be used as the first step of epitope screening, narrowing the scope of epitope search and reducing the time of epitope search.

TCR sequencing and cloning methods for repertoire analysis and isolation of tumor-reactive TCRs

T cell receptor (TCR) technologies, including repertoire analyses and T cell engineering, are increasingly important in the clinical management of cellular immunity in cancer, transplantation, and other immune diseases. However, sensitive and reliable methods for repertoire analyses and TCR cloning are still lacking. Here, we report on SEQTR, a high-throughput approach to analyze human and mouse repertoires that is more sensitive, reproducible, and accurate as compared with commonly used assays, and thus more reliably captures the complexity of blood and tumor TCR repertoires. We also present a TCR cloning strategy to specifically amplify TCRs from T cell populations. Positioned downstream of single-cell or bulk TCR sequencing, it allows time- and cost-effective discovery, cloning, screening, and engineering of tumor-specific TCRs. Together, these methods will accelerate TCR repertoire analyses in discovery, translational, and clinical settings and permit fast TCR engineering for cellular therapies.

Identification and tracking of HTLV-1-infected T cell clones in virus-associated neurologic disease

Human T lymphotropic virus type 1-assoicated (HTLV-1-associated) myelopathy/tropical spastic paraparesis (HAM/TSP) is a neuroinflammatory disease caused by the persistent proliferation of HTLV-1-infected T cells. Here, we performed a T cell receptor (TCR) repertoire analysis focused on HTLV-1-infected cells to identify and track the infected T cell clones that are preserved in patients with HAM/TSP and migrate to the CNS. TCRβ repertoire analysis revealed higher clonal expansion in HTLV-1-infected cells compared with noninfected cells from patients with HAM/TSP and asymptomatic carriers (ACs). TCR clonality in HTLV-1-infected cells was similar in patients with HAM/TSP and ACs. Longitudinal analysis showed that the TCR repertoire signature in HTLV-1-infected cells remained stable, and highly expanded infected clones were preserved within each patient with HAM/TSP over years. Expanded HTLV-1-infected clones revealed different distributions between cerebrospinal fluid (CSF) and peripheral blood and were enriched in the CSF of patients with HAM/TSP. Cluster analysis showed similarity in TCRβ sequences in HTLV-1-infected cells, suggesting that they proliferate after common antigen stimulation. Our results indicate that exploring TCR repertoires of HTLV-1-infected cells can elucidate individual clonal dynamics and identify potential pathogenic clones expanded in the CNS.

Profiling the epigenetic landscape of the antigen receptor repertoire: the missing epi-immunogenomics data

High-resolution sequencing methods that capture the epigenetic landscape within the T cell receptor (TCR) gene loci are pivotal for a fundamental understanding of the epigenetic regulatory mechanisms of the TCR repertoire. In our opinion, filling the gaps in our understanding of the epigenetic mechanisms regulating the TCR repertoire will benefit the development of strategies that can modulate the TCR repertoire composition by leveraging the dynamic nature of epigenetic modifications.

Regulation of loop extrusion on the interphase genome

In the human cell nucleus, dynamically organized chromatin is the substrate for gene regulation, DNA replication, and repair. A central mechanism of DNA loop formation is an ATPase motor cohesin-mediated loop extrusion. The cohesin complexes load and unload onto the chromosome under the control of other regulators that physically interact and affect motor activity. Regulation of the dynamic loading cycle of cohesin influences not only the chromatin structure but also genome-associated human disorders and aging. This review focuses on the recently spotlighted genome organizing factors and the mechanism by which their dynamic interactions shape the genome architecture in interphase.

Correcting inborn errors of immunity: From viral mediated gene addition to gene editing

Allogeneic hematopoietic stem cell transplantation is an effective treatment to cure inborn errors of immunity. Remarkable progress has been achieved thanks to the development and optimization of effective combination of advanced conditioning regimens and use of immunoablative/suppressive agents preventing rejection as well as graft versus host disease. Despite these tremendous advances, autologous hematopoietic stem/progenitor cell therapy based on ex vivo gene addition exploiting integrating γ-retro- or lenti-viral vectors, has demonstrated to be an innovative and safe therapeutic strategy providing proof of correction without the complications of the allogeneic approach. The recent advent of targeted gene editing able to precisely correct genomic variants in an intended locus of the genome, by introducing deletions, insertions, nucleotide substitutions or introducing a corrective cassette, is emerging in the clinical setting, further extending the therapeutic armamentarium and offering a cure to inherited immune defects not approachable by conventional gene addition. In this review, we will analyze the current state-of-the art of conventional gene therapy and innovative protocols of genome editing in various primary immunodeficiencies, describing preclinical models and clinical data obtained from different trials, highlighting potential advantages and limits of gene correction.

Molecular and functional characterization of two IgM subclasses in large yellow croaker (Larimichthys crocea)

As the predominant immunoglobulin (Ig) isotype, IgM plays a crucial role in the acquired immunity of vertebrates. There is only one Igμ gene in mammals, except cattle, while the number of Igμ gene varies among teleost fish. In the current study, we found two functional Igμ genes (Igμ1 and Igμ2) and a pseudo Cμ gene (ψIgμ) in large yellow croaker (Larimichthys crocea). Both Igμ1 and Igμ2 genes possessed two transcript variants, which encoded the heavy chains of secreted (sIgM1 and sIgM2) and membrane-bound IgM1 and IgM2 (mIgM1 and mIgM2), respectively. Both the heavy chains of sIgM1 and sIgM2 consisted of a variable Ig domain, four constant Ig domains (CH1, CH2, CH3 and CH4) and a secretory tail, while those of mIgM1 and mIgM2 consisted of a variable Ig domain, three constant Ig domains (CH1, CH2 and CH3), a transmembrane domain and a short cytoplasmic tail. Cysteine residues that are necessary for the formation of intrachain and interchain disulfide bonds and tryptophan residues that are important for the folding of the Ig superfamily domain were well conserved in large yellow croaker IgM1 and IgM2. Interestingly, large yellow croaker IgM2 had an extra cysteine (C94) in the CH1 domain compared with IgM1, which may cause the structural difference between IgM1 and IgM2. A liquid chromatography-tandem mass spectrometry analysis revealed that both IgM1 and IgM2 were present at the protein level in large yellow croaker serum. Both the Igμ1 and Igμ2 genes were mainly expressed in systemic immune tissues, such as head kidney and spleen, but the expression level of Igμ2 was much lower than that of Igμ1. After Pseudomonas plecoglossicida infection, the expression levels of Igμ1 and Igμ2 in both the spleen and head kidney were significantly upregulated, with a higher upregulation of Igμ2 than that of Igμ1. These results suggested that Igμ1 and Igμ2 may play a differential role in the immune response of large yellow croaker against bacterial infection.

Reconstructing B cell lineage trees with minimum spanning tree and genotype abundances

B cell receptor (BCR) genes exposed to an antigen undergo somatic hypermutations and Darwinian antigen selection, generating a large BCR-antibody diversity. This process, known as B cell affinity maturation, increases antibody affinity, forming a specific B cell lineage that includes the unmutated ancestor and mutated variants. In a B cell lineage, cells with a higher antigen affinity will undergo clonal expansion, while those with a lower affinity will not proliferate and probably be eliminated. Therefore, cellular (genotype) abundance provides a valuable perspective on the ongoing evolutionary process. Phylogenetic tree inference is often used to reconstruct B cell lineage trees and represents the evolutionary dynamic of BCR affinity maturation. However, such methods should process B-cell population data derived from experimental sampling that might contain different cellular abundances. There are a few phylogenetic methods for tracing the evolutionary events occurring in B cell lineages; best-performing solutions are time-demanding and restricted to analysing a reduced number of sequences, while time-efficient methods do not consider cellular abundances. We propose ClonalTree, a low-complexity and accurate approach to construct B-cell lineage trees that incorporates genotype abundances into minimum spanning tree (MST) algorithms. Using both simulated and experimental data, we demonstrate that ClonalTree outperforms MST-based algorithms and achieves a comparable performance to a method that explores tree-generating space exhaustively. Furthermore, ClonalTree has a lower running time, being more convenient for building B-cell lineage trees from high-throughput BCR sequencing data, mainly in biomedical applications, where a lower computational time is appreciable. It is hundreds to thousands of times faster than exhaustive approaches, enabling the analysis of a large set of sequences within minutes or seconds and without loss of accuracy. The source code is freely available at github.com/julibinho/ClonalTree.

Ataxia Telangiectasia Mutated and MSH2 Control Blunt DNA End Joining in Ig Class Switch Recombination

Class-switch recombination (CSR) produces secondary Ig isotypes and requires activation-induced cytidine deaminase (AID)-dependent DNA deamination of intronic switch regions within the IgH (Igh) gene locus. Noncanonical repair of deaminated DNA by mismatch repair (MMR) or base excision repair (BER) creates DNA breaks that permit recombination between distal switch regions. Ataxia telangiectasia mutated (ATM)-dependent phosphorylation of AID at serine 38 (pS38-AID) promotes its interaction with apurinic/apyrimidinic endonuclease 1 (APE1), a BER protein, suggesting that ATM regulates CSR through BER. However, pS38-AID may also function in MMR during CSR, although the mechanism remains unknown. To examine whether ATM modulates BER- and/or MMR-dependent CSR, Atm-/- mice were bred to mice deficient for the MMR gene mutS homolog 2 (Msh2). Surprisingly, the predicted Mendelian frequencies of Atm-/-Msh2-/- adult mice were not obtained. To generate ATM and MSH2-deficient B cells, Atm was conditionally deleted on an Msh2-/- background using a floxed ATM allele (Atmf) and B cell-specific Cre recombinase expression (CD23-cre) to produce a deleted ATM allele (AtmD). As compared with AtmD/D and Msh2-/- mice and B cells, AtmD/DMsh2-/- mice and B cells display a reduced CSR phenotype. Interestingly, Sμ-Sγ1 junctions from AtmD/DMsh2-/- B cells that were induced to switch to IgG1 in vitro showed a significant loss of blunt end joins and an increase in insertions as compared with wild-type, AtmD/D, or Msh2-/- B cells. These data indicate that the absence of both ATM and MSH2 blocks nonhomologous end joining, leading to inefficient CSR. We propose a model whereby ATM and MSH2 function cooperatively to regulate end joining during CSR through pS38-AID.

Cetuximab Responses in Patients with HNSCC Correlate to Clonal Expansion Feature of Peripheral and Tumor-Infiltrating T Cells with Top T-Cell Receptor Clonotypes

PURPOSE

Cetuximab is a standard-of-care treatment for head and neck squamous cell carcinoma (HNSCC). Well-defined correlative markers of therapeutic responses are still lacking. Characterizing dynamic changes of T-cell receptor (TCR) repertoire in peripheral blood and tumor tissue may facilitate developing markers for cetuximab response in HNSCCs.

EXPERIMENTAL DESIGN

We analyzed high-throughput TCRβ sequencing data generated with ImmunoSEQ platform using peripheral blood mononuclear cells (PBMC) and tumor-infiltrating lymphocytes (TIL) from patients with HNSCC before and after cetuximab treatment (pre-/post-PBMC vs. pre-/post-TIL). Multiple analytic approaches were employed to normalize sequencing data.

RESULTS

Normalized TCR richness was significantly lower in post-TIL than pre-TIL, suggesting that cetuximab reduced TCR diversity and promoted TCR expansion in TIL samples, regardless of response status. The magnitude of clonal expansion (defined as expansion rate) in top 20 TCR clonotypes was significantly higher in responder PBMC with or without normalization, and in responder TIL upon normalization, than nonresponder ones. Notably, the expanded top 20 or top 50 TCR clonotypes overlapped between PBMC and TIL samples, which occurred significantly more frequently in responders than nonresponders.

CONCLUSIONS

Patients with cetuximab-treated HNSCC harbor dynamic changes of TCR repertoires correlative to therapeutic responses. The expansion rate of top TCR clonotypes in peripheral blood may serve as a minimally invasive, readily accessible, and feasible marker for predicting cetuximab responses in HNSCCs and beyond, and the expansion rate of top TCR clonotypes in TILs and their overlapping probability between PBMC and TIL may serve as additional predictive markers. Our study also highlights the importance of data normalization for TCR repertoire analysis.

Construction of Peptide Library in Mammalian Cells by dsDNA-Based Strategy

While different display technologies, represented by phage display, have been widely used in drug discovery, they still can hardly achieve function-based peptide screening, which in most cases is performed in mammalian cells. And most attempts to screen functional peptides with mammalian platforms utilized plasmids to store coding information. Our previous work established double-stranded DNAs (dsDNAs) as innovative biological parts to implement AND-gate genetic circuits in mammalian cells. In the current study, we employ dsDNAs with terminal NNK degenerate codons to implement AND-gate genetic circuits and generate peptide libraries in mammalian cells. This dsDNA-based AND-gate (DBAG) peptide library construction strategy is easy to perform, requiring only PCR reaction and cell transfection. High-throughput sequencing (HTS) and single-cell sequencing results revealed both peptide length and amino acid sequence diversity of DBAG peptide libraries. Moreover, as a feasibility test of this strategy, we identified an MDM2-interacting peptide by applying the DBAG peptide library to a mammalian cell-based two-hybrid system. Our work establishes dsDNAs with terminal degenerate codons as biological parts to build peptide libraries in mammalian cells, which may have great application potential in the future.

Methods for Study of Mouse T Cell Receptor α and β Gene Rearrangements

Quantitative real-time PCR and next-generation sequencing (NGS) are invaluable techniques to analyze T cell receptor (Tcr) gene rearrangements in mouse lymphocyte populations. Although these approaches are powerful, they also have limitations that must be accounted for in experimental design and data interpretation. Here, we provide relevant background required for understanding these limitations and then outline established quantitative real-time PCR and NGS methods that can be used for analysis of mouse Tcra and Tcrb gene rearrangements in mice.

Characterization of chromatin accessibility patterns in different mouse cell types using machine learning methods at single-cell resolution

Chromatin accessibility is a generic property of the eukaryotic genome, which refers to the degree of physical compaction of chromatin. Recent studies have shown that chromatin accessibility is cell type dependent, indicating chromatin heterogeneity across cell lines and tissues. The identification of markers used to distinguish cell types at the chromosome level is important to understand cell function and classify cell types. In the present study, we investigated transcriptionally active chromosome segments identified by sci-ATAC-seq at single-cell resolution, including 69,015 cells belonging to 77 different cell types. Each cell was represented by existence status on 20,783 genes that were obtained from 436,206 active chromosome segments. The gene features were deeply analyzed by Boruta, resulting in 3897 genes, which were ranked in a list by Monte Carlo feature selection. Such list was further analyzed by incremental feature selection (IFS) method, yielding essential genes, classification rules and an efficient random forest (RF) classifier. To improve the performance of the optimal RF classifier, its features were further processed by autoencoder, light gradient boosting machine and IFS method. The final RF classifier with MCC of 0.838 was constructed. Some marker genes such as H2-Dmb2, which are specifically expressed in antigen-presenting cells (e.g., dendritic cells or macrophages), and Tenm2, which are specifically expressed in T cells, were identified in this study. Our analysis revealed numerous potential epigenetic modification patterns that are unique to particular cell types, thereby advancing knowledge of the critical functions of chromatin accessibility in cell processes.

Longitudinal transcriptional analysis of peripheral blood leukocytes in COVID-19 convalescent donors

BACKGROUND

SARS-CoV2 can induce a strong host immune response. Many studies have evaluated antibody response following SARS-CoV2 infections. This study investigated the immune response and T cell receptor diversity in people who had recovered from SARS-CoV2 infection (COVID-19).

METHODS

Using the nCounter platform, we compared transcriptomic profiles of 162 COVID-19 convalescent donors (CCD) and 40 healthy donors (HD). 69 of the 162 CCDs had two or more time points sampled.

RESULTS

After eliminating the effects of demographic factors, we found extensive differential gene expression up to 241 days into the convalescent period. The differentially expressed genes were involved in several pathways, including virus-host interaction, interleukin and JAK-STAT signaling, T-cell co-stimulation, and immune exhaustion. A subset of 21 CCD samples was found to be highly "perturbed," characterized by overexpression of PLAU, IL1B, NFKB1, PLEK, LCP2, IRF3, MTOR, IL18BP, RACK1, TGFB1, and others. In addition, one of the clusters, P1 (n = 8) CCD samples, showed enhanced TCR diversity in 7 VJ pairs (TRAV9.1_TCRVA_014.1, TRBV6.8_TCRVB_016.1, TRAV7_TCRVA_008.1, TRGV9_ENST00000444775.1, TRAV18_TCRVA_026.1, TRGV4_ENST00000390345.1, TRAV11_TCRVA_017.1). Multiplexed cytokine analysis revealed anomalies in SCF, SCGF-b, and MCP-1 expression in this subset.

CONCLUSIONS

Persistent alterations in inflammatory pathways and T-cell activation/exhaustion markers for months after active infection may help shed light on the pathophysiology of a prolonged post-viral syndrome observed following recovery from COVID-19 infection. Future studies may inform the ability to identify druggable targets involving these pathways to mitigate the long-term effects of COVID-19 infection.

TRIAL REGISTRATION

https://clinicaltrials.gov/ct2/show/NCT04360278 Registered April 24, 2020.

Single-cell genomics identifies distinct B1 cell developmental pathways and reveals aging-related changes in the B-cell receptor repertoire

Background

B1 cells are self-renewing innate-like B lymphocytes that provide the first line of defense against pathogens. B1 cells primarily reside in the peritoneal cavity and are known to originate from various fetal tissues, yet their developmental pathways and the mechanisms underlying maintenance of B1 cells throughout adulthood remain unclear.

Results

We performed high-throughput single-cell analysis of the transcriptomes and B-cell receptor repertoires of peritoneal B cells of neonates, young adults, and elderly mice. Gene expression analysis of 31,718 peritoneal B cells showed that the neonate peritoneal cavity contained many B1 progenitors, and neonate B cell specific clustering revealed two trajectories of peritoneal B1 cell development, including pre-BCR dependent and pre-BCR independent pathways. We also detected profound age-related changes in B1 cell transcriptomes: clear difference in senescence genetic program was evident in differentially aged B1 cells, and we found an example that a B1 subset only present in the oldest mice was marked by expression of the fatty-acid receptor CD36. We also performed antibody gene sequencing of 15,967 peritoneal B cells from the three age groups and discovered that B1 cell aging was associated with clonal expansion and two B1 cell clones expanded in the aged mice had the same CDR-H3 sequence (AGDYDGYWYFDV) as a pathogenically linked cell type from a recent study of an atherosclerosis mouse model.

Conclusions

Beyond offering an unprecedent data resource to explore the cell-to-cell variation in B cells, our study has revealed that B1 precursor subsets are present in the neonate peritoneal cavity and dissected the developmental pathway of the precursor cells. Besides, this study has found the expression of CD36 on the B1 cells in the aged mice. And the single-cell B-cell receptor sequencing reveals B1 cell aging is associated with clonal expansion.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13578-022-00795-6.

Perfluoroalkyl substances (PFASs) decrease the expression of recombination-activating genes (RAG1 and RAG2) in human B lymphoma Namalwa cells

Per- and polyfluoroalkyl substances (PFASs) are omnipresent and have been shown to induce a wide range of adverse effects, including hepatotoxicity, developmental toxicity and immunotoxicity. So far, little information is available about the mechanisms underlying the toxicity of PFASs, including those related to their immunotoxicity. Reported immunotoxic effects of PFASs include decreased antibody responses in experimental animals and humans, indicating that PFASs may, among others, affect B cell function. In the present study, we first assessed the effects of PFOA on the transcriptome of the human Namalwa B cell line using RNA seq analysis. Gene expression changes, analyzed using Ingenuity Pathway Analysis, pointed to various cellular processes affected by PFOA, including 'B cell development' and 'Primary immunodeficiency signaling'. Interestingly, PFOA decreased the expression of RAG1 and RAG2, genes involved in immunoglobulin and T cell receptor V(D)J recombination. As a next step, time- and concentration-dependent changes in the expression of RAG1 and RAG2 upon exposure to PFOA, PFNA, PFHxS and PFOS were studied through RT-qPCR analysis. Analysis with the concentration-response modeling software PROAST resulted in the following potency ranking: PFNA > PFOA > PFOS > PFHxS. Altogether, the present in vitro study provides insights into the effects of selected PFASs on B cells, identifying RAG1 and RAG2 expression as possible relevant targets that may play a role in the immunotoxicity of PFASs.

An Attention Based Bidirectional LSTM Method to Predict the Binding of TCR and Epitope

The T-cell epitope prediction has always been a long-term challenge in immunoinformatics and bioinformatics. Studying the specific recognition between T-cell receptor (TCR) and peptide-major histocompatibility complex (p-MHC) complexes can help us better understand the immune mechanism, it's also make a signification contribution in developing vaccines and targeted drugs. Meanwhile, more advanced methods are needed for distinguishing TCRs binding from different epitopes. In this paper, we introduce a hybrid model composed of bidirectional long short-term memory networks (BiLSTM), attention and convolutional neural networks (CNN) that can identified the binding of TCRs to epitopes. The BiLSTM can more completely extract amino acid forward and backward information in the sequence, and attention mechanism can focus on amino acids at certain positions from complex sequences to capture the most important feature, then CNN was used to further extract salient features to predict the binding of TCR-epitope. In McPAS dataset, the AUC value (the area under ROC curve) of naive TCR-epitope binding is 0.974 and specific TCR-epitope binding is 0.887. The model has achieved better prediction results than other existing models (TCRGP, ERGO, NetTCR), and some experiments are used to analyze the advantages of our model. The algorithm is available at https://github.com/bijingshu/BiAttCNN.git.

Divergent outcomes of anti-PD-L1 treatment coupled with host-intrinsic differences in TCR repertoire and distinct T cell activation states in responding versus non-responding tumors

Differential responses to immune checkpoint inhibitors (ICI) may be attributed to tumor-intrinsic factors or environmental cues; however, these mechanisms cannot fully explain the variable ICI responses in different individuals. Here, we investigate the potential contribution of immunological heterogeneity with a focus on differences in T-cell receptor (TCR) repertoire to ICI responses, which has not been defined previously. To reveal additional factors underlying heterogeneous responses to ICI, we employed a squamous cell carcinoma (SCC) mouse model in which tumor-bearing recipients unambiguously diverged into responders (R) or non-responders (NR) upon anti-PD-L1 treatment. Treatment efficacy absolutely required CD8 T-cells and correlated positively with effector functions of CD8 tumor-infiltrating lymphocytes (TILs). We showed that TCR repertoires exhibited a similar magnitude of clonal expansion in R vs. NR CD8 TILs. However, the top expanded TCR clonotypes appeared to be mutually exclusive between R and NR CD8 TILs, which also occurred in a recipient-specific manner, demonstrating preferential expansion of distinct TCR clonotypes against the same SCC tumor. Unexpectedly, R vs. NR CD8 TILs reached all activation clusters and did not exhibit substantial global differences in transcriptomes. By linking single-cell transcriptomic data with unique TCR clonotypes, CD8 TILs harboring top TCR clonotypes were found to occupy distinct activation clusters and upregulate genes favoring anti-tumor immunity to different extents in R vs. NR. We conclude that stochastic differences in CD8 TIL TCR repertoire and distinct activation states of top TCR clonotypes may contribute to differential anti-PD-L1 responses. Our study suggests that host-intrinsic immunological heterogeneity may offer a new explanation for differential ICI responses in different individuals, which could impact on strategies for personalized cancer immunotherapy.

Why responses to immune checkpoint inhibitors are heterogeneous in head and neck cancers: Contributions from tumor-intrinsic and host-intrinsic factors

Immune checkpoint inhibitors (ICIs) have revolutionized cancer treatment including in head and neck squamous cell carcinomas (HNSCCs); however, only a fraction of HNSCC patients respond to ICI, whereas the majority fail to do so. The mechanisms underlying such variable responses remain incompletely understood. A better understanding of such mechanisms may broaden the spectrum of responding patients and enhance the rate of ICI response. HNSCCs exhibit a high level of genetic heterogeneity, manifested as mutations or amplifications of oncogenes (e.g., PIK3CA) and mutations of tumor suppressor genes (e.g., TP53). The immune tumor microenvironment (TME) of HNSCCs also varies significantly in composition and in relative abundance of distinct immune subsets such as CD8 tumor-infiltrating lymphocytes (TILs) or tumor-associated macrophages (TAMs), which represents a high degree of immunological heterogeneity. Here, we briefly discuss how heterogeneous ICI responses may be attributed to tumor-intrinsic factors, including genetic, transcriptional, and functional variations in tumor cells, and host-intrinsic factors, including cellular composition of the TME (e.g., CD8 TILs and TAMs), and host-intrinsic differences in the T cell receptor (TCR) repertoire of CD8 TILs. We also discuss the potential impact of these factors on designing strategies for personalized immunotherapy of HNSCCs.

T cell Repertoire Profiling and the Mechanism by which HLA-B27 Causes Ankylosing Spondylitis

Purpose of Review

Ankylosing spondylitis (AS) is strongly associated with the HLA-B27 gene. The canonical function of HLA-B27 is to present antigenic peptides to CD8 lymphocytes, leading to adaptive immune responses. The ‘arthritogenic peptide’ theory as to the mechanism by which HLA-B27 induces ankylosing spondylitis proposes that HLA-B27 presents peptides derived from exogenous sources such as bacteria to CD8 lymphocytes, which subsequently cross-react with antigens at the site of inflammation of the disease, causing inflammation. This review describes findings of studies in AS involving profiling of T cell expansions and discusses future research opportunities based on these findings.

Recent Findings

Consistent with this theory, there is an expanding body of data showing that expansion of a restricted pool of CD8 lymphocytes is found in most AS patients yet only in a small proportion of healthy HLA-B27 carriers.

Summary

These exciting findings strongly support the theory that AS is driven by presentation of antigenic peptides to the adaptive immune system by HLA-B27. They point to new potential approaches to identify the exogenous and endogenous antigens involved and to potential therapies for the disease.

Ezrin Promotes Antigen Receptor Diversity during B Cell Development by Supporting Ig H Chain Variable Gene Recombination

Genome-level rearrangements of Ig genes during B cell development are critical for generation of a diverse repertoire of BCRs that bind to a multitude of foreign Ags and some self Ags. Bone marrow B cell development involves a variety of cell–cell interactions, cell migration, and receptor signaling that likely benefit from the activity of membrane-cytoskeletal reorganizing proteins. However, the specific contribution of such proteins toward BCR repertoire diversification is poorly understood. Ezrin is a membrane-cytoskeletal linker protein that regulates mature B cell activation through spatial organization of the BCR. We employed next-generation sequencing to investigate whether Ezrin plays a role in IgH rearrangements and generation of BCR diversity in developing bone marrow B cells. BCR repertoire development occurred stochastically in B cell progenitors from both control and B cell conditional Ezrin-deficient mice. However, the loss of Ezrin resulted in fewer unique CDRs (CDR3s) in the BCRs and reduced Shannon entropy. Ezrin-deficient pre-B cells revealed similar utilization of joining (J) genes but significantly fewer variable (V) genes, thereby decreasing V-J combinatorial diversity. V-J junctional diversity, measured by CDR3 length and nucleotide additions and deletions, was not altered in Ezrin-deficient pre-B cells. Mechanistically, Ezrin-deficient cells showed a marked decrease in RAG1 gene expression, indicating a less efficient DNA recombination machinery. Overall, our results demonstrate that Ezrin shapes the BCR repertoire through combinatorial diversification.

Evaluation and comparison of adaptive immunity through analyzing the diversities and clonalities of T-cell receptor repertoires in the peripheral blood

The adaptive immune system plays an important role in defending against different kinds of diseases, including infection and cancer. There has been a longtime need for a simple method to quantitatively evaluate the potency of adaptive immunity in our bodies. The tremendously diversified T-cell receptor (TCR) repertoires are the foundation of the adaptive immune system. In this study, we analyzed the expressed TCRβ repertoires in the peripheral blood of 582 healthy donors and 60 cancer patients. The TCR repertoire in each individual is different, with different usages of TCR Vβ and Jβ genes. Importantly, the TCR diversity and clonality change along with age and disease situation. Most elder individuals and cancer patients have elevated numbers of large TCRβ clones and reduced numbers of shared common clones, and thus, they have very low TCR diversity index (D₅₀) values. These results reveal the alteration of the expressed TCRβ repertoire with aging and oncogenesis, and thus, we hypothesize that the TCR diversity and clonality in the peripheral blood might be used to evaluate and compare the adaptive immunities among different individuals in clinical practice.

The RAG1 Ubiquitin Ligase Domain Stimulates Recombination of TCRβ and TCRα Genes and Influences Development of αβ T Cell Lineages

RAG1/RAG2 (RAG) endonuclease-mediated assembly of diverse lymphocyte Ag receptor genes by V(D)J recombination is critical for the development and immune function of T and B cells. The RAG1 protein contains a ubiquitin ligase domain that stabilizes RAG1 and stimulates RAG endonuclease activity in vitro. We report in this study that mice with a mutation that inactivates the Rag1 ubiquitin ligase in vitro exhibit decreased rearrangements and altered repertoires of TCRβ and TCRα genes in thymocytes and impaired thymocyte developmental transitions that require the assembly and selection of functional TCRβ and/or TCRα genes. These Rag1 mutant mice present diminished positive selection and superantigen-mediated negative selection of conventional αβ T cells, decreased genesis of invariant NK T lineage αβ T cells, and mature CD4+ αβ T cells with elevated autoimmune potential. Our findings reveal that the Rag1 ubiquitin ligase domain functions in vivo to stimulate TCRβ and TCRα gene recombination and influence differentiation of αβ T lineage cells, thereby establishing replete diversity of αβ TCRs and populations of αβ T cells while restraining generation of potentially autoreactive conventional αβ T cells.

A high-quality severe combined immunodeficiency (SCID) rat bioresource

Immunodeficient animals are valuable models for the engraftment of exogenous tissues; they are widely used in many fields, including the creation of humanized animal models, as well as regenerative medicine and oncology. Compared with mice, laboratory rats have a larger body size and can more easily undergo transplantation of various tissues and organs. Considering the absence of high-quality resources of immunodeficient rats, we used the CRISPR/Cas9 genome editing system to knock out the interleukin-2 receptor gamma chain gene (Il2rg) in F344/Jcl rats—alone or together with recombination activating gene 2 (Rag2)—to create a high-quality bioresource that researchers can freely use: severe combined immunodeficiency (SCID) rats. We selected one founder rat with frame-shift mutations in both Il2rg (5-bp del) and Rag2 ([1-bp del+2-bp ins]/[7-bp del+2-bp ins]), then conducted mating to establish a line of immunodeficient rats. The immunodeficiency phenotype was preliminarily confirmed by the presence of severe thymic hypoplasia in Il2rg-single knockout (sKO) and Il2rg/Rag2-double knockout (dKO) rats. Assessment of blood cell counts in peripheral blood showed that the white blood cell count was significantly decreased in sKO and dKO rats, while the red blood cell count was unaffected. The decrease in white blood cell count was mainly caused by a decrease in lymphocytes. Furthermore, analyses of lymphocyte populations via flow cytometry showed that the numbers of B cells (CD3^- CD45⁺) and natural killer cells (CD3^- CD161⁺) were markedly reduced in both knockout rats. In contrast, T cells were markedly reduced but showed slightly different results between sKO and dKO rats. Notably, our immunodeficient rats do not exhibit growth retardation or gametogenesis defects. This high-quality SCID rat resource is now managed by the National BioResource Project in Japan. Our SCID rat model has been used in various research fields, demonstrating its importance as a bioresource.

Genomic crossroads between non-Hodgkin’s lymphoma and common variable immunodeficiency

Common variable immunodeficiency (CVID) represents the largest group of primary immunodeficiencies that may manifest with infections, inflammation, autoimmunity, and cancer, mainly B-cell non-Hodgkin’s lymphoma (NHL). Indeed, NHL may result from chronic or recurrent infections and has, therefore, been recognized as a clinical phenotype of CVID, although rare. The more one delves into the mechanisms involved in CVID and cancer, the stronger the idea that both pathologies can be a reflection of the same primer events observed from different angles. The potential effects of germline variants on specific somatic modifications in malignancies suggest that it might be possible to anticipate critical events during tumor development. In the same way, a somatic alteration in NHL could be conditioning a similar response at the transcriptional level in the shared signaling pathways with genetic germline alterations in CVID. We aimed to explore the genomic substrate shared between these entities to better characterize the CVID phenotype immunodeficiency in NHL. By means of an in-silico approach, we interrogated the large, publicly available datasets contained in cBioPortal for the presence of genes associated with genetic pathogenic variants in a panel of 50 genes recurrently altered in CVID and previously described as causative or disease-modifying. We found that 323 (25%) of the 1,309 NHL samples available for analysis harbored variants of the CVID spectrum, with the most recurrent alteration presented in NHL occurring in PIK3CD (6%) and STAT3 (4%). Pathway analysis of common gene alterations showed enrichment in inflammatory, immune surveillance, and defective DNA repair mechanisms similar to those affected in CVID, with PIK3R1 appearing as a central node in the protein interaction network. The co-occurrence of gene alterations was a frequent phenomenon. This study represents an attempt to identify common genomic grounds between CVID and NHL. Further prospective studies are required to better know the role of genetic variants associated with CVID and their reflection on the somatic pathogenic variants responsible for cancer, as well as to characterize the CVID-like phenotype in NHL, with the potential to influence early CVID detection and therapeutic management.