Thioredoxin is a metabolic rheostat controlling regulatory B cells.

Metabolic programming is important for B cell fate, but the bioenergetic requirement for regulatory B (Breg) cell differentiation and function is unknown. Here we show that Breg cell differentiation, unlike non-Breg cells, relies on mitochondrial electron transport and homeostatic levels of reactive oxygen species (ROS). Single-cell RNA sequencing analysis revealed that TXN, encoding the metabolic redox protein thioredoxin (Trx), is highly expressed by Breg cells, unlike Trx inhibitor TXNIP which was downregulated. Pharmacological inhibition or gene silencing of TXN resulted in mitochondrial membrane depolarization and increased ROS levels, selectively suppressing Breg cell differentiation and function while favoring pro-inflammatory B cell differentiation. Patients with systemic lupus erythematosus (SLE), characterized by Breg cell deficiencies, present with B cell mitochondrial membrane depolarization, elevated ROS and fewer Trx+ B cells. Exogenous Trx stimulation restored Breg cells and mitochondrial membrane polarization in SLE B cells to healthy B cell levels, indicating Trx insufficiency underlies Breg cell impairment in patients with SLE.

Epithelial damage and tissue γδ T cells promote a unique tumor-protective IgE response.

IgE is an ancient and conserved immunoglobulin isotype with potent immunological function. Nevertheless, the regulation of IgE responses remains an enigma, and evidence of a role for IgE in host defense is limited. Here we report that topical exposure to a common environmental DNA-damaging xenobiotic initiated stress surveillance by γδTCR+ intraepithelial lymphocytes that resulted in class switching to IgE in B cells and the accumulation of autoreactive IgE. High-throughput antibody sequencing revealed that γδ T cells shaped the IgE repertoire by supporting specific variable-diversity-joining (VDJ) rearrangements with unique characteristics of the complementarity-determining region CDRH3. This endogenous IgE response, via the IgE receptor FcεRI, provided protection against epithelial carcinogenesis, and expression of the gene encoding FcεRI in human squamous-cell carcinoma correlated with good disease prognosis. These data indicate a joint role for immunosurveillance by T cells and by B cells in epithelial tissues and suggest that IgE is part of the host defense against epithelial damage and tumor development.

sciCSR infers B cell state transition and predicts class-switch recombination dynamics using single-cell transcriptomic data.

Class-switch recombination (CSR) is an integral part of B cell maturation. Here we present sciCSR (pronounced 'scissor', single-cell inference of class-switch recombination), a computational pipeline that analyzes CSR events and dynamics of B cells from single-cell RNA sequencing (scRNA-seq) experiments. Validated on both simulated and real data, sciCSR re-analyzes scRNA-seq alignments to differentiate productive heavy-chain immunoglobulin transcripts from germline 'sterile' transcripts. From a snapshot of B cell scRNA-seq data, a Markov state model is built to infer the dynamics and direction of CSR. Applying sciCSR on severe acute respiratory syndrome coronavirus 2 vaccination time-course scRNA-seq data, we observe that sciCSR predicts, using data from an earlier time point in the collected time-course, the isotype distribution of B cell receptor repertoires of subsequent time points with high accuracy (cosine similarity ~0.9). Using processes specific to B cells, sciCSR identifies transitions that are often missed by conventional RNA velocity analyses and can reveal insights into the dynamics of B cell CSR during immune response.

CD1d-mediated lipid presentation by CD11c+ cells regulates intestinal homeostasis.

Intestinal homeostasis relies on a continuous dialogue between the commensal bacteria and the immune system. Natural killer T (NKT) cells, which recognize CD1d-restricted microbial lipids and self-lipids, contribute to the regulation of mucosal immunity, yet the mechanisms underlying their functions remain poorly understood. Here, we demonstrate that NKT cells respond to intestinal lipids and CD11c+ cells (including dendritic cells (DCs) and macrophages) are essential to mediate lipid presentation within the gut ultimately controlling intestinal NKT cell homeostasis and activation. Conversely, CD1d and NKT cells participate in the control of the intestinal bacteria composition and compartmentalization, in the regulation of the IgA repertoire and in the induction of regulatory T cells within the gut. These changes in intestinal homeostasis require CD1d expression on DC/macrophage populations as mice with conditional deletion of CD1d on CD11c+ cells exhibit dysbiosis and altered immune homeostasis. These results unveil the importance of CD11c+ cells in controlling lipid-dependent immunity in the intestinal compartment and reveal an NKT cell-DC crosstalk as a key mechanism for the regulation of gut homeostasis.

Multi-Omics Reveals Mechanisms of Partial Modulation of COVID-19 Dysregulation by Glucocorticoid Treatment.

Treatments for COVID-19 infections have improved dramatically since the beginning of the pandemic, and glucocorticoids have been a key tool in improving mortality rates. The UK's National Institute for Health and Care Excellence guidance is for treatment to be targeted only at those requiring oxygen supplementation, however, and the interactions between glucocorticoids and COVID-19 are not completely understood. In this work, a multi-omic analysis of 98 inpatient-recruited participants was performed by quantitative metabolomics (using targeted liquid chromatography-mass spectrometry) and data-independent acquisition proteomics. Both 'omics datasets were analysed for statistically significant features and pathways differentiating participants whose treatment regimens did or did not include glucocorticoids. Metabolomic differences in glucocorticoid-treated patients included the modulation of cortisol and bile acid concentrations in serum, but no alleviation of serum dyslipidemia or increased amino acid concentrations (including tyrosine and arginine) in the glucocorticoid-treated cohort relative to the untreated cohort. Proteomic pathway analysis indicated neutrophil and platelet degranulation as influenced by glucocorticoid treatment. These results are in keeping with the key role of platelet-associated pathways and neutrophils in COVID-19 pathogenesis and provide opportunity for further understanding of glucocorticoid action. The findings also, however, highlight that glucocorticoids are not fully effective across the wide range of 'omics dysregulation caused by COVID-19 infections.

Immunoglobulin gene analysis as a tool for investigating human immune responses.

The human immunoglobulin repertoire is a hugely diverse set of sequences that are formed by processes of gene rearrangement, heavy and light chain gene assortment, class switching and somatic hypermutation. Early B cell development produces diverse IgM and IgD B cell receptors on the B cell surface, resulting in a repertoire that can bind many foreign antigens but which has had self-reactive B cells removed. Later antigen-dependent development processes adjust the antigen affinity of the receptor by somatic hypermutation. The effector mechanism of the antibody is also adjusted, by switching the class of the antibody from IgM to one of seven other classes depending on the required function. There are many instances in human biology where positive and negative selection forces can act to shape the immunoglobulin repertoire and therefore repertoire analysis can provide useful information on infection control, vaccination efficacy, autoimmune diseases, and cancer. It can also be used to identify antigen-specific sequences that may be of use in therapeutics. The juxtaposition of lymphocyte development and numerical evaluation of immune repertoires has resulted in the growth of a new sub-speciality in immunology where immunologists and computer scientists/physicists collaborate to assess immune repertoires and develop models of immune action.

BRepertoire: a user-friendly web server for analysing antibody repertoire data.

Antibody repertoire analysis by high throughput sequencing is now widely used, but a persisting challenge is enabling immunologists to explore their data to discover discriminating repertoire features for their own particular investigations. Computational methods are necessary for large-scale evaluation of antibody properties. We have developed BRepertoire, a suite of user-friendly web-based software tools for large-scale statistical analyses of repertoire data. The software is able to use data preprocessed by IMGT, and performs statistical and comparative analyses with versatile plotting options. BRepertoire has been designed to operate in various modes, for example analysing sequence-specific V(D)J gene usage, discerning physico-chemical properties of the CDR regions and clustering of clonotypes. Those analyses are performed on the fly by a number of R packages and are deployed by a shiny web platform. The user can download the analysed data in different table formats and save the generated plots as image files ready for publication. We believe BRepertoire to be a versatile analytical tool that complements experimental studies of immune repertoires. To illustrate the server's functionality, we show use cases including differential gene usage in a vaccination dataset and analysis of CDR3H properties in old and young individuals. The server is accessible under http://mabra.biomed.kcl.ac.uk/BRepertoire.

Significant association between TNFAIP3 inactivation and biased immunoglobulin heavy chain variable region 4-34 usage in mucosa-associated lymphoid tissue lymphoma.

Both antigenic drive and genetic change play critical roles in the development of mucosa-associated lymphoid tissue (MALT) lymphoma, but neither alone is sufficient for malignant transformation, and lymphoma development critically depends on their cooperation. However, which of these different events concur and how they cooperate in MALT lymphomagenesis is totally unknown. To explore this, we investigated somatic mutations of 17 genes and immunoglobulin heavy chain variable region (IGHV) usage in 179 MALT lymphomas from various sites. We showed that: (1) there was a significant association between the biased usage of IGHV4-34 (binds to the carbohydrate I/i antigens) and inactivating mutation of TNFAIP3 [encoding a global negative regulator of the canonical nuclear factor-κB (NF-κB) pathway] in ocular adnexal MALT lymphoma; (2) IGHV1-69 was significantly overrepresented (54%) in MALT lymphoma of the salivary gland, but was not associated with mutation in any of the 17 genes investigated; and (3) MALT lymphoma lacked mutations that are frequently seen in other B-cell lymphomas characterized by constitutive NF-κB activities, including mutations in CD79B, CARD11, MYD88, TNFRSF11A, and TRAF3. Our findings show, for the first time, a significant association between biased usage of autoreactive IGHV and somatic mutation of NF-κB regulators in MALT lymphoma, arguing for their cooperation in sustaining chronic B-cell receptor signalling and driving oncogenesis in lymphoma development. Copyright © 2017 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

Aging affects B-cell antigen receptor repertoire diversity in primary and secondary lymphoid tissues.

The elderly immune system is characterized by reduced responses to infections and vaccines, and an increase in the incidence of autoimmune diseases and cancer. Age-related deficits in the immune system may be caused by peripheral homeostatic pressures that limit bone marrow B-cell production or migration to the peripheral lymphoid tissues. Studies of peripheral blood B-cell receptor spectratypes have shown that those of the elderly are characterized by reduced diversity, which is correlated with poor health status. In the present study, we performed for the first time high-throughput sequencing of immunoglobulin genes from archived biopsy samples of primary and secondary lymphoid tissues in old (74 ± 7 years old, range 61-89) versus young (24 ± 5 years old, range 18-45) individuals, analyzed repertoire diversities and compared these to results in peripheral blood. We found reduced repertoire diversity in peripheral blood and lymph node repertoires from old people, while in the old spleen samples the diversity was larger than in the young. There were no differences in somatic hypermutation characteristics between age groups. These results support the hypothesis that age-related immune frailty stems from altered B-cell homeostasis leading to narrower memory B-cell repertoires, rather than changes in somatic hypermutation mechanisms.

A Reassessment of IgM Memory Subsets in Humans.

From paired blood and spleen samples from three adult donors, we performed high-throughput VH sequencing of human B cell subsets defined by IgD and CD27 expression: IgD(+)CD27(+) ("marginal zone [MZ]"), IgD(-)CD27(+) ("memory," including IgM ["IgM-only"], IgG and IgA) and IgD(-)CD27(-) cells ("double-negative," including IgM, IgG, and IgA). A total of 91,294 unique sequences clustered in 42,670 clones, revealing major clonal expansions in each of these subsets. Among these clones, we further analyzed those shared sequences from different subsets or tissues for VH gene mutation, H-CDR3-length, and VH/JH usage, comparing these different characteristics with all sequences from their subset of origin for which these parameters constitute a distinct signature. The IgM-only repertoire profile differed notably from that of MZ B cells by a higher mutation frequency and lower VH4 and higher JH6 gene usage. Strikingly, IgM sequences from clones shared between the MZ and the memory IgG/IgA compartments showed a mutation and repertoire profile of IgM-only and not of MZ B cells. Similarly, all IgM clonal relationships (among MZ, IgM-only, and double-negative compartments) involved sequences with the characteristics of IgM-only B cells. Finally, clonal relationships between tissues suggested distinct recirculation characteristics between MZ and switched B cells. The "IgM-only" subset (including cells with its repertoire signature but higher IgD or lower CD27 expression levels) thus appear as the only subset showing precursor-product relationships with CD27(+) switched memory B cells, indicating that they represent germinal center-derived IgM memory B cells and that IgM memory and MZ B cells constitute two distinct entities.

Redemption of autoantibodies on anergic B cells by variable-region glycosylation and mutation away from self-reactivity.

The best-understood mechanisms for achieving antibody self/non-self discrimination discard self-reactive antibodies before they can be tested for binding microbial antigens, potentially creating holes in the repertoire. Here we provide evidence for a complementary mechanism: retaining autoantibodies in the repertoire displayed as low levels of IgM and high IgD on anergic B cells, masking a varying proportion of autoantibody-binding sites with carbohydrates, and removing their self-reactivity by somatic hypermutation and selection in germinal centers (GCs). Analysis of human antibody sequences by deep sequencing of isotype-switched memory B cells or in IgG antibodies elicited against allogeneic RhD+ erythrocytes, vaccinia virus, rotavirus, or tetanus toxoid provides evidence for reactivation of anergic IgM(low) IgD+ IGHV4-34+ B cells and removal of cold agglutinin self-reactivity by hypermutation, often accompanied by mutations that inactivated an N-linked glycosylation sequon in complementarity-determining region 2 (CDR2). In a Hy10 antibody transgenic model where anergic B cells respond to a biophysically defined lysozyme epitope displayed on both foreign and self-antigens, cell transfers revealed that anergic IgM(low) IgD+ B cells form twice as many GC progeny as naïve IgM(hi) IgD+ counterparts. Their GC progeny were rapidly selected for CDR2 mutations that blocked 72% of antigen-binding sites with N-linked glycan, decreased affinity 100-fold, and then cleared the binding sites of blocking glycan. These results provide evidence for a mechanism to acquire self/non-self discrimination by somatic mutation away from self-reactivity, and reveal how varying the efficiency of N-glycosylation provides a mechanism to modulate antibody avidity.

Immunoglobulin light chain allelic inclusion in systemic lupus erythematosus.

The principles of allelic exclusion state that each B cell expresses a single light and heavy chain pair. Here, we show that B cells with both kappa and lambda light chains (Igκ and Igλ) are enriched in some patients with the systemic autoimmune disease systemic lupus erythematosus (SLE), but not in the systemic autoimmune disease control granulomatosis with polyangiitis. Detection of dual Igκ and Igλ expression by flow cytometry could not be abolished by acid washing or by DNAse treatment to remove any bound polyclonal antibody or complexes, and was retained after two days in culture. Both surface and intracytoplasmic dual light chain expression was evident by flow cytometry and confocal microscopy. We observed reduced frequency of rearrangements of the kappa-deleting element (KDE) in SLE and an inverse correlation between the frequency of KDE rearrangement and the frequency of dual light chain expressing B cells. We propose that dual expression of Igκ and Igλ by a single B cell may occur in some patients with SLE when this may be a consequence of reduced activity of the KDE.

Yersinia enterocolitica provides the link between thyroid-stimulating antibodies and their germline counterparts in Graves' disease.

Graves' disease results from thyroid-stimulating Abs (TSAbs) activating the thyrotropin receptor (TSHR). How TSAbs arise from early precursor B cells has not been established. Genetic and environmental factors may contribute to pathogenesis, including the bacterium Yersinia enterocolitica. We developed two pathogenic monoclonal TSAbs from a single experimental mouse undergoing Graves' disease, which shared the same H and L chain germline gene rearrangements and then diversified by numerous somatic hypermutations. To address the Ag specificity of the shared germline precursor of the monoclonal TSAbs, we prepared rFab germline, which showed negligible binding to TSHR, indicating importance of somatic hypermutation in acquiring TSAb activity. Using rFab chimeras, we demonstrate the dominant role of the H chain V region in TSHR recognition. The role of microbial Ags was tested with Y. enterocolitica proteins. The monoclonal TSAbs recognize 37-kDa envelope proteins, also recognized by rFab germline. MALDI-TOF identified the proteins as outer membrane porin (Omp) A and OmpC. Using recombinant OmpA, OmpC, and related OmpF, we demonstrate cross-reactivity of monoclonal TSAbs with the heterogeneous porins. Importantly, rFab germline binds recombinant OmpA, OmpC, and OmpF confirming reactivity with Y. enterocolitica. A human monoclonal TSAb, M22 with similar properties to murine TSAbs, also binds recombinant porins, showing cross-reactivity of a spontaneously arising pathogenic Ab with Y. enterocolitica. The data provide a mechanistic framework for molecular mimicry in Graves' disease, where early precursor B cells are expanded by Y. enterocolitica porins to undergo somatic hypermutation to acquire a cross-reactive pathogenic response to TSHR.

Granulocyte colony stimulating factor exacerbates antineutrophil cytoplasmic antibody vasculitis.

OBJECTIVES: Granulocyte colony stimulating factor (GCSF) is important in mobilising neutrophils from the bone marrow but also has a range of proinflammatory effects. We therefore decided to investigate the role of GCSF in antineutrophil cytoplasmic antibody (ANCA) vasculitis. METHODS: We measured GCSF levels in the serum of 38 patients with active ANCA vasculitis compared with 31 age-matched controls, and assessed the effect of GCSF priming on the response of human neutrophils to ANCA. We also examined the effect of exogenous GCSF administration in a murine model of antimyeloperoxidase (anti-MPO) vasculitis, and the effect of GCSF on murine neutrophil activation. RESULTS: The serum levels of GCSF in patients with active ANCA vasculitis were significantly higher than those of age matched healthy controls (mean 38.04 vs 18.35 pg/ml, p<0.001). Furthermore, we demonstrated that GCSF primed human neutrophils in vitro for a respiratory burst in response to anti-MPO ANCA. In an anti-MPO antibody transfer model, mice given GCSF had more crescents (mean 29.1% vs 5.8% per glomerular cross section, p<0.05), more macrophages (mean 3.2 vs 1.2 per glomerular cross-section, p<0.01), higher serum creatines (mean 13.6 vs 8.3 µmol/l, p<0.05) and more haematuria (p<0.05) compared with controls. In vivo administration of GCSF with lipopolysaccharide (LPS), but not LPS alone, led to upregulation of CD11c on murine neutrophils. CONCLUSIONS: These data suggest that GCSF, which is raised in patient serum, may play an important role in exacerbating disease in ANCA vasculitis. In addition, GCSF therapy for neutropenia should be used with caution in these patients.

BCR CDR3 length distributions differ between blood and spleen and between old and young patients, and TCR distributions can be used to detect myelodysplastic syndrome.

Complementarity-determining region 3 (CDR3) is the most hyper-variable region in B cell receptor (BCR) and T cell receptor (TCR) genes, and the most critical structure in antigen recognition and thereby in determining the fates of developing and responding lymphocytes. There are millions of different TCR Vß chain or BCR heavy chain CDR3 sequences in human blood. Even now, when high-throughput sequencing becomes widely used, CDR3 length distributions (also called spectratypes) are still a much quicker and cheaper method of assessing repertoire diversity. However, distribution complexity and the large amount of information per sample (e.g. 32 distributions of the TCRα chain, and 24 of TCRß) calls for the use of machine learning tools for full exploration. We have examined the ability of supervised machine learning, which uses computational models to find hidden patterns in predefined biological groups, to analyze CDR3 length distributions from various sources, and distinguish between experimental groups. We found that (a) splenic BCR CDR3 length distributions are characterized by low standard deviations and few local maxima, compared to peripheral blood distributions; (b) healthy elderly people's BCR CDR3 length distributions can be distinguished from those of the young; and (c) a machine learning model based on TCR CDR3 distribution features can detect myelodysplastic syndrome with approximately 93% accuracy. Overall, we demonstrate that using supervised machine learning methods can contribute to our understanding of lymphocyte repertoire diversity.

PCR amplification and high throughput sequencing of immunoglobulin heavy chain genes from formalin-fixed paraffin-embedded human biopsies.

The use of high throughput sequencing (HTS) technologies in biomedicine is expanding in a variety of fields in recent years. The 454 system is an HTS platform that is ideally suited to characterize B cell receptor (BCR) repertoires by sequencing of immunoglobulin (Ig) genes, as it is able to sequence stretches of several hundred nucleotides. Most studies that used this platform for antibody repertoire analyses have started from fresh or frozen tissues or peripheral blood samples, and rely on starting with optimal quality DNA. In this paper we demonstrate that BCR repertoire analysis can be done using DNA from formalin-fixed paraffin-embedded (FFPE) human tissue samples. The heterogeneity of BCR repertoires we obtained confirms the plausibility of HTS of DNA from FFPE specimens. The establishment of experimental protocols and computational tools that enable sequence data analysis from the low quality DNA of FFPE tissues is important for enabling research, as it would enable the use of the rich source of preserved samples in clinical biobanks and biopsy archives.

Generation of Immunoglobulin diversity in human gut-associated lymphoid tissue.

The organised gut associated lymphoid tissue (GALT) exists adjacent to an extensive and diverse luminal flora. The follicle associated epithelium and associated dendritic cells and lymphocytes form a tightly fortified gateway between the flora and the host that permits connectivity between them and chronic activation of the lymphoid compartment. As a consequence, plasma cell precursors are generated continuously, and in abundance, in GALT by clonal proliferation. Clonal proliferation alone on this scale would reduce the spectrum of B cell specificity. To compensate, GALT also houses molecular machinery that diversifies the receptor repertoire by somatic hypermutation, class switch recombination and receptor revision. These three processes of enhancing the diversity of mature B cells ensure that although clonally related plasma cells may secrete immunoglobulin side by side in the mucosa they rarely have identical antigen binding sites.

Diversification of immunoglobulin genes by gene conversion in the domestic chicken (Gallus gallus domesticus).

Sustainable modern poultry production depends on effective protection against infectious diseases and a diverse range of antibodies is key for an effective immune response. In the domestic chicken, somatic gene conversion is the dominant process in which the antibody immunoglobulin genes are diversified. Affinity maturation by somatic hypermutation (SHM) also occurs, but the relative contribution of gene conversion versus somatic hypermutation to immunoglobulin (Ig) gene diversity is poorly understood. In this study, we use high throughput long-read sequencing to study immunoglobulin diversity in multiple immune-associated tissues in Rhode Island Red chickens. To better understand the impact of genetic diversification in the chicken, a novel gene conversion identification software was developed (BrepConvert). In this study, BrepConvert enabled the identification of over 1 million gene conversion events. Mapping the occurrence of putative somatic gene conversion (SGC) events throughout the variable gene region revealed repetitive and highly restricted patterns of genetic insertions in both the antibody heavy and light chains. These patterns coincided with the locations of genetic variability in available pseudogenes and align with antigen binding sites, predominately the complementary determining regions (CDRs). We found biased usage of pseudogenes during gene conversion, as well as immunoglobulin heavy chain diversity gene (IGHD) preferences during V(D)J gene rearrangement, suggesting that antibody diversification in chickens is more focused than the genetic potential for diversity would suggest.

B cell profiles, antibody repertoire and reactivity reveal dysregulated responses with autoimmune features in melanoma.

B cells are known to contribute to the anti-tumor immune response, especially in immunogenic tumors such as melanoma, yet humoral immunity has not been characterized in these cancers to detail. Here we show comprehensive phenotyping in samples of circulating and tumor-resident B cells as well as serum antibodies in melanoma patients. Memory B cells are enriched in tumors compared to blood in paired samples and feature distinct antibody repertoires, linked to specific isotypes. Tumor-associated B cells undergo clonal expansion, class switch recombination, somatic hypermutation and receptor revision. Compared with blood, tumor-associated B cells produce antibodies with proportionally higher levels of unproductive sequences and distinct complementarity determining region 3 properties. The observed features are signs of affinity maturation and polyreactivity and suggest an active and aberrant autoimmune-like reaction in the tumor microenvironment. Consistent with this, tumor-derived antibodies are polyreactive and characterized by autoantigen recognition. Serum antibodies show reactivity to antigens attributed to autoimmune diseases and cancer, and their levels are higher in patients with active disease compared to post-resection state. Our findings thus reveal B cell lineage dysregulation with distinct antibody repertoire and specificity, alongside clonally-expanded tumor-infiltrating B cells with autoimmune-like features, shaping the humoral immune response in melanoma.