Cartilaginous fish possess one of the most ancient adaptive immune systems, and they uniquely produce the heavy chain-only antibody, immunoglobulin novel antigen receptor (IgNAR). In this study, we explored the mRNA transcription of genes related to antibody production and IgNAR diversity in various organs in banded houndsharks. IgNAR and antibody production-related genes exhibited similar relative transcription levels, with the highest expression detected in the spleen. Subsequently, we examined the diversity of IgNAR using next-generation sequencing. The most frequent clones were dominant (25 %-40 %) in the epigonal organ and liver but less common in the spleen. Large individual variation was noted in the kidney and pancreas. The length of complementarity-determining region 3 ranged 2-39 amino acids. The region tended to have a narrow length distribution of approximately 13 amino acids in the epigonal organ and liver, whereas wider length variation was noted in the kidney, pancreas, and spleen. Type II IgNAR variable regions (VNARs) were predominant (60 %-96 %) in all organs, whereas Type IV and "other" not conventionally defined VNARs were present at low frequencies and in different proportion between organs. Type I VNARs were present in multiple organs. The VNAR sequences were commonly shared among the epigonal organ, liver, and/or pancreas, but few were shared in the kidney or spleen.

B cells have important functions in gut homeostasis, and dysregulated B cell populations are frequently observed in patients with inflammatory bowel diseases, including both ulcerative colitis (UC) and Crohn's disease (CD). How these B cell perturbations contribute to disease remains largely unknown. Here, we perform deep sequencing of the B cell receptor (BCR) repertoire in four cohorts of patients with CD, together with healthy controls and patients with UC. We identify BCR clones that are shared between patients with CD but not found in healthy individuals nor in patients with UC, indicating CD-associated B cell immune responses. Shared clones are present in the inflamed gut mucosa, draining intestinal lymph nodes and blood, suggesting the presence of common CD-associated antigens that drive B cell responses in CD patients.

The COVID-19 pandemic has spurred a global race for a preventive vaccine, with a few becoming available just one year after describing this novel coronavirus disease. Among these are inactivated virus vaccines like CoronaVac (Sinovac Biotech), which are used in several countries to reduce the pandemic's effects. However, its use was associated with low protection, particularly against novel virus variants that quickly appeared in the following months. Vaccines play a crucial role in activating the immune system to combat infections, with Memory B-cells being a key part of this mechanism, eliciting protective neutralizing antibodies. This work focused on studying B-cell memory repertoire after two consecutive doses of CoronaVac.

Memory B-cells were isolated from five CoronaVac vaccinated and five pre-pandemic individuals and subsequently stimulated in vitro before high-throughput Illumina sequencing of the Heavy Chain Variable repertoire.

We observed a shift in the VH repertoire with increased HCDR3 length and enrichment of IGVH 3-23, 3-30, 3-7, 3-72, and 3-74 for IgA BCRs and IGHV 4-39 and 4-59 for IgG BCRs. A high expansion of IgA-specific clonal populations was observed in vaccinated individuals relative to pre-pandemic controls, accompanied by shared IgA variable heavy chain (VH) sequences among memory B cells across different vaccine recipients of IgA clones was also observed in vaccinated individuals compared to pre-pandemic controls, with several IgA VH sharing between memory B cells from different vaccines. Moreover, a high convergence was observed among vaccinees and SARS-CoV-2 neutralizing antibody sequences found in the CoV-abDab database.

These data show the ability of CoronaVac to elicit antibodies with characteristics similar to those previously identified as neutralizing antibodies, supporting its protective efficacy. Furthermore, this analysis of the immunological repertoire in the context of viral infections reinforces the importance of immunization in generating convergent antibodies for the antiviral response.

Rare naive B cells have special pathogen-recognition features that enable outsized contributions to protective immunity but infrequently participate in immune responses. We investigatee how germline-targeting vaccine delivery and adjuvant selection affect priming of exceptionally rare BG18-like HIV broadly neutralizing antibody-precursor B cells (<1-in-50 million) in non-human primates. Only escalating dose (ED) priming immunization using the saponin adjuvant SMNP elicited detectable BG18-like cells in germinal centers (GCs) compared with other conditions. All groups had strong GC responses, but only ED+SMNP and bolus+SMNP induced BG18-like memory B cells in >50% of animals. One group had vaccine-specific GC responses equivalent to ED+SMNP but scarce BG18-like B cells. Following homologous boosting, BG18-like memory B cells were present in a bolus priming group but with lower somatic hypermutation and affinities than ED+SMNP. This outcome inversely associated with post-prime antibody titers, suggesting antibody feedback significantly influences rare precursor B cell responses. Thus, antigen and inflammatory stimuli extensively impact priming and affinity maturation of rare B cells.

Sjögren's disease (SjD) is a chronic autoimmune disorder characterized by increased circulating self-reactive B cells. While many of these self-reactive B cells emerge from the bone marrow, it is not known whether they are excluded from or enriched in specific developmental stages in the periphery. The aim of this study was to determine the immunophenotype of circulating self-reactive B cells in SjD to inform more precise therapeutic targeting. Five major B cell populations: transitional, mature naïve, switched memory, double negative and plasmablasts were single-cell sorted and cultured to produce IgG. Self-reactive IgG was identified by ELISA, flow cytometry of permeabilized HEK293 cells and HEp-2 indirect immunofluorescence. Immunoglobulin heavy chains were sequenced by Sanger and next-generation sequencing. Compared with healthy donor controls (HCs), SjD patients had higher frequencies of naïve and CD21low atypical memory B cell subsets, while antigen-experienced B cells expressed more Ki67 and CD86. B cells recognizing intracellular self-antigens were identified in all stages of peripheral B cell development for SjD and HCs, but frequencies of autoreactive B cells were up to 10-fold higher in SjD. Self-reactive transitional B cells expressed higher surface CD38 and lower surface IgM. An increase in self-reactive B cells throughout peripheral development in SjD compared with HCs suggests that counterselection of autoantibody-bearing B cells during central and peripheral tolerance checkpoints are reduced in SjD. Therapeutic strategies focused on depleting B cells based on B cell receptor specificity rather than the developmental stage would be more efficient to target self-reactive B cells in SjD.

Early immune dynamics during the initiation of fatal tularemia caused by Francisella tularensis infection remain unknown. Unto that end, we generated a transcriptomic map at single-cell resolution of the innate-like lymphocyte responses to F. tularensis live vaccine strain (LVS) infection of mice. We found that both interferon-γ (IFN-γ)-producing type 1 and interleukin-17 (IL-17)-producing type 3 innate-like lymphocytes expanded in the infected lungs. Natural killer (NK) and NKT cells drove the type 1 response, whereas mucosal-associated invariant T (MAIT) and γδ T cells drove the type 3 response. Furthermore, tularemia-like disease resistant NKT cell-deficient, Cd1d -/- mice accumulated more MAIT1 cells, MAIT17 cells, and cells with a hybrid phenotype between MAIT1 and MAIT17 cells than wild-type mice. Critically, adoptive transfer of LVS-activated MAIT cells from Cd1d -/- mice, which were enriched in MAIT17 cells, was sufficient to protect LVS-susceptible, immunodeficient RAG2 -/- mice from severe LVS infection-inflicted pathology. Collectively, our findings position MAIT cells as potential mediators of IL-17-dependent protection from pulmonary tularemia-like disease.

Rhesus macaques (RMs) are a vital model for studying human disease and invaluable to pre-clinical vaccine research, particularly for the study of broadly neutralizing antibody responses. Such studies require robust genetic resources for antibody-encoding genes within the immunoglobulin (IG) loci. The complexity of the IG loci has historically made them challenging to characterize accurately. To address this, we developed novel experimental and computational methodologies to generate the largest collection to date of integrated antibody repertoire and long-read genomic sequencing data in 106 Indian origin RMs. We created a comprehensive resource of IG heavy and light chain variable (V), diversity (D), and joining (J) alleles, as well as leader, intronic, and recombination signal sequences (RSSs), including the curation of 1474 novel alleles, unveiling tremendous diversity, and expanding existing IG allele sets by 60%. This publicly available, continually updated resource (https://vdjbase.org/reference_book/Rhesus_Macaque) provides the foundation for advancing RM immunogenomics, vaccine discovery, and translational research.

Myalgic encephalomyelitis/chronic fatigue syndrome (ME/CFS) is a common and debilitating chronic illness of unknown aetiology. Chronic infection and autoimmune responses have been proposed as two mechanisms that potentially underlie the pathogenesis of ME/CFS. To explore these disease hypotheses, we characterised the antigen-specific receptors of B cells using adaptive immune receptor repertoire sequencing. We compared the B-cell receptor (BCR) repertoires of 25 patients with mild/moderate ME/CFS, 36 patients with severe ME/CFS, 21 healthy controls, and 28 patients with multiple sclerosis (MS) to identify signatures of infection or autoimmune responses. ME/CFS patients did not display increased clonality or differential somatic hypermutation compared to healthy controls and patients with MS. One of two immunoglobulin heavy variable (IGHV) genes, IGHV3-30, reported to be increased in ME/CFS patients in a previous study, was replicated in patients with mild/moderate disease in our cohort. However, there was no evidence of ongoing adaptive responses in IGHV3-30 repertoires from mild/moderate ME/CFS patients with increased IGHV3-30 usage. There were no detectable repertoire signatures associated with infection or autoimmunity in repertoires from ME/CFS patients, but we observed skewing of the ratio of IgM to IgG BCRs in patients with mild/moderate ME/CFS, a preliminary finding that presents an opportunity for follow-up work.

The major histocompatibility complex (MHC) class I-related molecule MHC-class-I-related protein 1 (MR1) presents metabolites to distinct MR1-restricted T cell subsets, including mucosal-associated invariant T (MAIT) and MR1T cells. However, self-reactive MR1T cells and the nature of recognized antigens remain underexplored. Here, we report a cell endogenous carbonyl adduct of adenine (8-(9H-purin-6-yl)-2-oxa-8-azabicyclo[3.3.1]nona-3,6-diene-4,6-dicarbaldehyde [M3Ade]) sequestered in the A' pocket of MR1. M3Ade induced in vitro MR1-mediated stimulation of MR1T cell clones that bound MR1-M3Ade tetramers. MR1-M3Ade tetramers identified heterogeneous MR1-reactive T cells ex vivo in healthy donors, individuals with acute myeloid leukemia, and tumor-infiltrating lymphocytes from non-small cell lung adenocarcinoma and hepatocarcinoma. These cells displayed phenotypic, transcriptional, and functional diversity at distinct differentiation stages, indicating their adaptive nature. They were also polyclonal, with some preferential T cell receptor (TCRαβ) pair usage. Thus, M3Ade is an MR1-presented self-metabolite that enables stimulation and tracking of human-MR1T cells from blood and tissue, aiding our understanding of their roles in health and disease.

B cells with low or absent expression of CD21 (CD21lo B cells) gained attention due to their expansion in the peripheral blood of patients with immune-mediated, rheumatic diseases. This is not only observed in typical autoimmune diseases like systemic lupus erythematosus and Sjögren's disease (SjD) but also in radiographic axial spondyloarthritis (r-axSpA), which is considered an autoinflammatory disease. To gain more insight into the origins of the heterogeneous CD21lo B-cell population, and its relation to the plasmablast (PB) compartment, we profiled the B-cell-receptor (BCR) repertoire in CD27- and CD27+ fractions of CD21lo B cells and early PBs using next-generation sequencing. Populations were sorted from peripheral blood of healthy individuals, SjD patients, and r-axSpA patients (n = 10 for each group). In healthy individuals and both patient groups, our findings indicate that CD27-CD21lo B cells, which exhibit few mutations in their BCR, may develop into CD27+CD21lo B cells and PBs, both marked by considerably more mutations. Given the known expansion of circulating CD27-CD21lo B cells in SjD and r-axSpA patients and clonal relationships with both CD27+CD21lo B cells and early PBs, these cells might actively contribute to (pathological) immune responses in rheumatic diseases with autoimmune and/or autoinflammatory characteristics.

Lymphatic filariasis is caused by infections of thread-like filarial worms, namely Wuchereria bancrofti, Brugia Malayi and Brugia timori. However, in-depth analysis of the antibody repertoire against Lymphatic filariasis is lacking. Using high-throughput sequencing of antibody repertoires, immunome analysis of IgG (LG) and IgM (LM) repertoires were studied. Despite significant differences between LG and LM in V(D)J gene usage, IGHV4-34, IGHV6-1, IGHD3-10 and IGHJ4 were preferred in both repertoires. The CDR3 in the LG repertoire showed a longer length than LM. Higher SHM level were observed in LG sequences and presence of oligoclonal sequences indicates the extent of clonal expansion. The prevalence of rare clonotypes in LM repertoire depicts the high clonal diversity when compared to LG repertoire. Monoclonal antibodies against closely related parasitic infections were present within the LG repertoire suggesting that immune repertoires may not be as exclusive and biased against the target infection as initially thought. The characterization of the immune repertoire can provide critical insight into the antibody response patterns in disease state, antibody generation process during infections and future antibody designs.

The distinctive characteristics of an individual's T cell receptor repertoire are crucial in recognizing and responding to a diverse array of antigens, contributing to immune specificity and adaptability. The repertoire, famously vast due to a series of cellular mechanisms, can be quantified using repertoire sequencing. In this study, we sampled the repertoire of 85 women: ovarian cancer patients (OC) and healthy donors (HD), generating a dataset of T cell clones and their abundance. For the alpha chain we obtained 6.4·106 reads, with an average of 75936 clones per sample, and an average of 30607 clonotypes per sample. For the beta chain we obtained 13.6·106 reads, with an average of 160400 clones per sample, and an average of 70071 clonotypes per sample. The changes in dynamics of the repertoire can be observed in response to disease, with specific clones undergoing clonal expansion and contraction. The data provided here offers a unique view of immune system behavior in health and disease and can be used to stratify OC and HD.

In the past decade, single-cell (sc) transcriptomics has overcome the limitations of bulk analysis by measuring gene expression in individual cells, not just a population average. This can identify diverse cell types and states within a sample with high resolution, even without prior purification. Various technologies exist, each with its own capture, barcoding, and library preparation methods. This chapter focuses on the analysis of normal and malignant mature B cells using the 10× Genomics 5' sc-gene expression in parallel with B cell immune repertoire profiling. By integrating the gene expression data from similar cells, the complete transcriptome for each population can be reconstructed, while the identification of the expressed immunoglobulin genes allows investigating clonotype evolution and the detection of tumor clones that share the same clonally rearranged B cell receptor sequence. Researchers are guided through both the experimental protocols and data analysis with a comprehensive, step-by-step walkthrough of how to use some of the more popular single-cell software tools.

Allelic variability in the adaptive immune receptor loci, which harbor the gene segments that encode B cell and T cell receptors (BCR/TCR), is of critical importance for immune responses to pathogens and vaccines. Adaptive immune receptor repertoire sequencing (AIRR-seq) has become widespread in immunology research making it the most readily available source of information about allelic diversity in immunoglobulin (IG) and T cell receptor (TR) loci. Here, we present a novel algorithm for extrasensitive and specific variable (V) and joining (J) gene allele inference, allowing the reconstruction of individual high-quality gene segment libraries. The approach can be applied for inferring allelic variants from peripheral blood lymphocyte BCR and TCR repertoire sequencing data, including hypermutated isotype-switched BCR sequences, thus allowing high-throughput novel allele discovery from a wide variety of existing data sets. The developed algorithm is a part of the MiXCR software. We demonstrate the accuracy of this approach using AIRR-seq paired with long-read genomic sequencing data, comparing it to a widely used algorithm, TIgGER. We applied the algorithm to a large set of IG heavy chain (IGH) AIRR-seq data from 450 donors of ancestrally diverse population groups, and to the largest reported full-length TCR alpha and beta chain (TRA and TRB) AIRR-seq data set, representing 134 individuals. This allowed us to assess the genetic diversity within the IGH, TRA, and TRB loci in different populations and to establish a database of alleles of V and J genes inferred from AIRR-seq data and their population frequencies with free public access through VDJ.online database.

Diverse antibody repertoires spanning multiple lymphoid organs (i.e., bone marrow, spleen, lymph nodes) form the foundation of protective humoral immunity. Changes in their composition across lymphoid organs are a consequence of B-cell selection and migration events leading to a highly dynamic and unique physiological landscape of antibody repertoires upon antigenic challenge (e.g., vaccination). However, to what extent B cells encoding identical or similar antibody sequences (clones) are distributed across multiple lymphoid organs and how this is shaped by the strength of a humoral response remains largely unexplored. Here, we performed an in-depth systems analysis of antibody repertoires across multiple distinct lymphoid organs of immunized mice and discovered that organ-specific antibody repertoire features (i.e., germline V-gene usage and clonal expansion profiles) equilibrated upon a strong humoral response (multiple immunizations and high serum titers). This resulted in a surprisingly high degree of repertoire consolidation, characterized by highly connected and overlapping B-cell clones across multiple lymphoid organs. Finally, we revealed distinct physiological axes indicating clonal migrations and showed that antibody repertoire consolidation directly correlated with antigen specificity. Our study uncovered how a strong humoral response resulted in a more uniform but redundant physiological landscape of antibody repertoires, indicating that increases in antibody serum titers were a result of synergistic contributions from antigen-specific B-cell clones distributed across multiple lymphoid organs. Our findings provide valuable insights for the assessment and design of vaccine strategies.

SARS-CoV-2 mRNA vaccines induce robust and persistent germinal centre (GC) B cell responses in humans. It remains unclear how the continuous evolution of the virus impacts the breadth of the induced GC B cell response. Using ultrasound-guided fine needle aspiration, we examined draining lymph nodes of nine healthy adults following bivalent booster immunization. We show that 77.8% of the B cell clones in the GC expressed as representative monoclonal antibodies recognized the spike protein, with a third (37.8%) of these targeting the receptor binding domain (RBD). Strikingly, only one RBD-targeting mAb, mAb-52, neutralized all tested SARS-CoV-2 strains, including the recent KP.2 variant. mAb-52 utilizes the IGHV3-66 public clonotype, protects hamsters challenged against the EG.5.1 variant and targets the class I/II RBD epitope, closely mimicking the binding footprint of ACE2. Finally, we show that the remarkable breadth of mAb-52 is due to the somatic hypermutations accumulated within vaccine-induced GC reaction.

The human gut includes plasma cells (PCs) expressing immunoglobulin A1 (IgA1) or IgA2, two structurally distinct IgA subclasses with elusive regulation, function, and reactivity. We show here that intestinal IgA1+ and IgA2+ PCs co-emerged early in life, comparably accumulated somatic mutations, and were enriched within short-lived CD19+ and long-lived CD19- PC subsets, respectively. IgA2+ PCs were extensively clonally related to IgA1+ PCs and a subset of them presumably emerged from IgA1+ precursors. Of note, secretory IgA1 (SIgA1) and SIgA2 dually coated a large fraction of mucus-embedded bacteria, including Akkermansia muciniphila. Disruption of homeostasis by inflammatory bowel disease (IBD) was associated with an increase in actively proliferating IgA1+ plasmablasts, a depletion in long-lived IgA2+ PCs, and increased SIgA1+SIgA2+ gut microbiota. Such increase featured enhanced IgA1 reactivity to pathobionts, including Escherichia coli, combined with depletion of beneficial A. muciniphila. Thus, gut IgA1 and IgA2 emerge from clonally related PCs and show unique changes in both frequency and reactivity in IBD.

Somatic hypermutation (SHM) is the diversity-generating process in antibody affinity maturation. Probabilistic models of SHM are needed for analyzing rare mutations, for understanding the selective forces guiding affinity maturation, and for understanding the underlying biochemical process. High throughput data offers the potential to develop and fit models of SHM on relevant data sets. In this paper we model SHM using modern frameworks. We are motivated by recent work suggesting the importance of a wider context for SHM, however, assigning an independent rate to each k-mer leads to an exponential proliferation of parameters. Thus, using convolutions on 3-mer embeddings, we develop "thrifty" models of SHM that have fewer free parameters than a 5-mer model and yet have a significantly wider context. These offer a slight performance improvement over a 5-mer model. We also find that a per-site effect is not necessary to explain SHM patterns given nucleotide context. Also, the two current methods for fitting an SHM model - on out-of-frame sequence data and on synonymous mutations - produce significantly different results, and augmenting out-of-frame data with synonymous mutations does not aid out-of-sample performance.

Whereas severe COVID-19 is often associated with elevated autoantibody titers, the underlying mechanism behind their generation has remained unclear. Here we report clonal composition and diversity of autoantibodies in humoral response to SARS-CoV-2. Immunoglobulin repertoire analysis and characterization of plasmablast-derived monoclonal antibodies uncovered clonal expansion of plasmablasts producing cardiolipin (CL)-reactive autoantibodies. Half of the expanded CL-reactive clones exhibited strong binding to SARS-CoV-2 antigens. One such clone, CoV1804, was reactive to both CL and viral nucleocapsid (N), and further showed anti-nucleolar activity in human cells. Notably, antibodies sharing genetic features with CoV1804 were identified in COVID-19 patient-derived immunoglobulins, thereby constituting a novel public antibody. These public autoantibodies had numerous mutations that unambiguously enhanced anti-N reactivity, when causing fluctuations in anti-CL reactivity along with the acquisition of additional self-reactivities, such as anti-nucleolar activity, in the progeny. Thus, potentially CL-reactive precursors may have developed multiple self-reactivities through clonal selection, expansion, and somatic hypermutation driven by viral antigens. Our results revealed the nature of autoantibody production during COVID-19 and provided novel insights into the origin of virus-induced autoantibodies.

Multiple alterations in the T cell repertoire were identified in many chronic inflammatory diseases such as inflammatory bowel disease, multiple sclerosis, and rheumatoid arthritis, suggesting that T cells might, directly or indirectly, be implicated in these pathologies. This has sparked a deep interest in characterizing disease-associated T cell clonotypes as well as in identifying and quantifying their contribution to the pathophysiology of different autoimmune and inflammation-mediated diseases. Bulk T cell repertoire sequencing (TCR-Seq) has emerged as a powerful method to profile the T cell repertoire of a sample in a high throughput fashion. Given the increasing utilization of TCR-Seq, we aimed here to provide a comprehensive, up-to-date review of the method, its extensions, and its ability to investigate chronic and autoimmune diseases. Specifically, we started by introducing the immunological basis of TCR repertoire generation and features, followed by discussing different experimental approach to perform TCR-Seq, then we describe different methods and frameworks for analyzing the generated datasets. Subsequently, different experimental techniques for investigating the antigenicity of T cell clonotypes are described. Lastly, we discuss recent studies that utilized TCR-Seq to understand different inflammation-mediated diseases, discuss fallbacks of the technology and potential future directions to overcome current limitations.

Autoimmune diseases such as lupus are characterized by polyclonal B cell activation, leading to the production of autoantibodies. The mechanism leading to B cell dysregulation is unclear; however, the defect may lie in selection within germinal centers (GCs). GC B cells cycle between proliferation and mutation in the dark zone and selection in the light zone (LZ). Temporal assessment of GCs from mice with either persistent infection or lupus showed an accumulation of LZ B cells. Yet, only in lupus, GC B cells exhibited reduced proliferation and progressive loss of MYC and FOXO1, which regulate zonal recycling and differentiation. As lupus progressed, decreased mutational frequency and repertoire diversity were associated with reduced responsiveness to CD40 signaling, despite accumulation of plasma cells. Collectively, these findings suggest that lupus disease progression coincides with an intrinsic defect in LZ B cell signaling, altering the zonal recycling, selection, and differentiation of autoreactive B cells.

Mycosis fungoides (MF) represents the most prevalent entity of cutaneous T cell lymphoma (CTCL). The MF aetiopathogenesis is incompletely understood, due to significant transcriptomic heterogeneity and conflicting views on whether oncologic transformation originates in early thymocytes or mature effector memory T cells. Recently, using clinical specimens, our group showed that the skin microbiome aggravates disease course, mainly driven by an outgrowing, pathogenic S. aureus strain carrying the virulence factor spa, which was shown by others to activate the T cell signalling pathway NF-κB.

To explore the role of the skin microbiome in MF aetiopathogenesis, we here performed RNA sequencing, multi-omic data integration of the skin microbiome and skin transcriptome using Multi-Omic Factor Analysis (MOFA), virome profiling, and T cell receptor (TCR) sequencing in 10 MF patients from our previous study group.

We observed that inter-patient transcriptional heterogeneity may be largely attributed to differential activation of T cell signalling pathways. Notably, the MOFA model resolved the heterogenous activation pattern of T cell signalling after denoising the transcriptome from microbial influence. The MOFA model suggested that the outgrowing S. aureus strain evoked signalling by non-canonical NF-κB and IL-1B, which in turn may have fuelled the aggravated disease course. Further, the MOFA model indicated aberrant pathways of early thymopoiesis alongside enrichment of antiviral innate immunity. In line with this, viral prevalence, particularly of Epstein-Barr virus (EBV), trended higher in both lesional skin and the blood compared to nonlesional skin. Additionally, TCRs in both MF skin lesions and the blood were significantly more likely to recognize EBV peptides involved in latent infection.

First, our findings suggest that S. aureus with its virulence factor spa fuels MF progression through non-canonical NF-κB and IL-1B signalling. Second, our data provide insights into the potential role of viruses in MF aetiology. Last, we propose a model of microbiome-driven MF aetiopathogenesis: Thymocytes undergo initial oncologic transformation, potentially caused by viruses. After maturation and skin infiltration, an outgrowing, pathogenic S. aureus strain evokes activation and maturation into effector memory T cells, resulting in aggressive disease. Further studies are warranted to verify and extend our data, which are based on computational analyses.

Rare B cells can have special pathogen-recognition features giving them the potential to make outsized contributions to protective immunity. However, rare naive B cells infrequently participate in immune responses. We investigated how germline-targeting vaccine antigen delivery and adjuvant selection affect priming of exceptionally rare BG18-like HIV broadly neutralizing antibody-precursor B cells (~1 in 50 million) in non-human primates. Only escalating dose (ED) priming immunization using the saponin adjuvant SMNP elicited detectable BG18-like cells in germinal centers (GCs). All groups had strong GC responses, but only ED+SMNP and bolus+SMNP induced BG18-like memory B cells in >50% of animals. One group had vaccine-specific GC responses equivalent to ED+SMNP, but BG18-like memory B cells were rarely detected. Following homologous boosting, BG18-like memory B cells were more frequent in a bolus priming group, but had lower somatic hypermutation and affinities. This outcome was inversely associated with post-prime antibody titers, suggesting antibody feedback can significantly influence rare precursor B cell responses.

Primary antiphospholipid syndrome (PAPS) is a life-threatening clotting disorder mediated by pathogenic autoantibodies. Here we dissect the origin of self-reactive B cells in human PAPS using peripheral blood and bone marrow of patients with triple-positive PAPS via combined single-cell RNA sequencing, B cell receptors (BCR) repertoire profiling, CITEseq analysis and single cell immortalization. We find that antiphospholipid (aPL)-specific B cells are present in the naive compartment, polyreactive, and derived from the natural repertoire. Furthermore, B cells with aPL specificities are not eliminated in patients with PAPS, persist until the memory and long-lived plasma cell stages, likely after defective germinal center selection, while becoming less polyreactive. Lastly, compared with the non-PAPS cells, PAPS B cells exhibit distinct IFN and APRIL signature as well as dysregulated mTORC1 and MYC pathways. Our findings may thus elucidate the survival mechanisms of these autoreactive B cells and suggest potential therapeutic targets for the treatment of PAPS.

Single-cell RNA sequencing (scRNA-seq) can resolve transcriptional features from individual cells, but scRNA-seq techniques capable of resolving the variable regions of B cell receptors (BCRs) remain limited, especially from widely-used 3'-barcoded libraries. Here, we report a method that can recover paired, full-length variable region sequences of BCRs from 3'-barcoded scRNA-seq libraries. We first verify this method (B3E-seq) can produce accurate, full-length BCR sequences. We then apply this method to profile B cell responses elicited against the capsular polysaccharide of Streptococcus pneumoniae serotype 3 (ST3) by glycoconjugate vaccines in five infant rhesus macaques. We identify BCR features associated with specificity for the ST3 antigen which are present in multiple vaccinated monkeys, indicating a convergent response to vaccination. These results demonstrate the utility of our method to resolve key features of the B cell repertoire and profile antigen-specific responses elicited by vaccination.

Immune checkpoint blockade (ICB) enhances T cell responses against cancer, leading to long-term survival in a fraction of patients. CD8+ T cell differentiation in response to chronic antigen stimulation is highly complex, and it remains unclear precisely which T cell differentiation states at which anatomic sites are critical for the response to ICB. We identified an intermediate-exhausted population in the white pulp of the spleen that underwent substantial expansion in response to ICB and gave rise to tumor-infiltrating clonotypes. Increased systemic antigen redirected differentiation of this population toward a more circulatory exhausted KLR state, whereas a lack of cross-presented tumor antigen reduced its differentiation in the spleen. An analogous population of exhausted KLR CD8+ T cells in human blood samples exhibited diminished tumor-trafficking ability. Collectively, our data demonstrate the critical role of antigen density within the spleen for the differentiation and expansion of T cell clonotypes in response to ICB.

Bats are reservoirs of many zoonotic viruses that are fatal in humans but do not cause disease in bats. Moreover, bats generate low neutralizing antibody titers in response to experimental viral infection, although more robust antibody responses have been observed in wild-caught bats during times of food stress. Here, we compared the antibody titers and B cell receptor (BCR) diversity of Jamaican fruit bats (Artibeus jamaicensis; JFBs) and BALB/c mice generated in response to T-dependent and T-independent antigens. We then manipulated the diet of JFBs and challenged them with H18N11 influenza A-like virus or a replication incompetent Nipah virus VSV (Nipah-riVSV). Under standard housing conditions, JFBs generated a lower avidity antibody response and possessed more BCR mRNA diversity compared to BALB/c mice. However, withholding protein from JFBs improved serum neutralization in response to Nipah-riVSV and improved serum antibody titers specific to H18 but reduced BCR mRNA diversity.

The sequencing of PCR amplicons is a core application of high-throughput sequencing technology. Using unique molecular identifiers (UMIs), individual amplified molecules can be sequenced to very high accuracy on an Illumina sequencer. However, Illumina sequencers have limited read length and are therefore restricted to sequencing amplicons shorter than 600 bp unless using inefficient synthetic long-read approaches. Native long-read sequencers from Pacific Biosciences and Oxford Nanopore Technologies can, using consensus read approaches, match or exceed Illumina quality while achieving much longer read lengths. Using a circularization-based concatemeric consensus sequencing approach (R2C2) paired with UMIs (R2C2 + UMI), we show that we can sequence an ∼550-nt antibody heavy chain (Immunoglobulin heavy chain - IGH) and an ∼1,500-nt 16S amplicons at accuracies up to and exceeding Q50 (<1 error in 100,000 sequenced bases), which exceeds accuracies of UMI-supported Illumina-paired sequencing as well as synthetic long-read approaches.

Mycosis fungoides (MF) is the most common entity of Cutaneous T cell lymphomas (CTCL) and is characterized by the presence of clonal malignant T cells in the skin. The role of the skin microbiome for MF development and progression are currently poorly understood. Using shotgun metagenomic profiling, real-time qPCR, and T cell receptor sequencing, we compared lesional and nonlesional skin of 20 MF patients with early and advanced MF. Additionally, we isolated Staphylococcus aureus and other bacteria from MF skin for functional profiling and to study the S. aureus virulence factor spa. We identified a subgroup of MF patients with substantial dysbiosis on MF lesions and concomitant outgrowth of S. aureus on plaque-staged lesions, while the other MF patients had a balanced microbiome on lesional skin. Dysbiosis and S. aureus outgrowth were accompanied by ectopic levels of cutaneous antimicrobial peptides (AMPs), including adaptation of the plaque-derived S. aureus strain. Furthermore, the plaque-derived S. aureus strain showed a reduced susceptibility towards antibiotics and an upregulation of the virulence factor spa, which may activate the NF-κB pathway. Remarkably, patients with dysbiosis on MF lesions had a restricted T cell receptor repertoire and significantly lower event-free survival. Our study highlights the potential for microbiome-modulating treatments targeting S. aureus to prevent MF progression.

B-cell repertoires are characterized by a diverse set of receptors of distinct specificities generated through two processes of somatic diversification: V(D)J recombination and somatic hypermutations. B-cell clonal families stem from the same V(D)J recombination event, but differ in their hypermutations. Clonal families identification is key to understanding B-cell repertoire function, evolution, and dynamics. We present HILARy (high-precision inference of lineages in antibody repertoires), an efficient, fast, and precise method to identify clonal families from single- or paired-chain repertoire sequencing datasets. HILARy combines probabilistic models that capture the receptor generation and selection statistics with adapted clustering methods to achieve consistently high inference accuracy. It automatically leverages the phylogenetic signal of shared mutations in difficult repertoire subsets. Exploiting the high sensitivity of the method, we find the statistics of evolutionary properties such as the site frequency spectrum and dN/dS ratio do not depend on the junction length. We also identify a broad range of selection pressures spanning two orders of magnitude.

Germinal centers (GC) are microanatomical lymphoid structures where affinity-matured memory B cells and long-lived bone marrow plasma cells are primarily generated. It is unclear how the maturation of B cells within the GC impacts the breadth and durability of B cell responses to influenza vaccination in humans. We used fine needle aspiration of draining lymph nodes to longitudinally track antigen-specific GC B cell responses to seasonal influenza vaccination. Antigen-specific GC B cells persisted for at least 13 wk after vaccination in two out of seven individuals. Monoclonal antibodies (mAbs) derived from persisting GC B cell clones exhibit enhanced binding affinity and breadth to influenza hemagglutinin (HA) antigens compared with related GC clonotypes isolated earlier in the response. Structural studies of early and late GC-derived mAbs from one clonal lineage in complex with H1 and H5 HAs revealed an altered binding footprint. Our study shows that inducing sustained GC reactions after influenza vaccination in humans supports the maturation of responding B cells.