Search

Age-associated B cells (ABCs) are formed under inflammatory conditions and are considered a type of memory B cell (MBC) expressing the transcription factor T-bet. In SLE, ABC frequency is correlated with disease, and they are thought to be the source of autoantibody-secreting cells. However, in inflammatory conditions, whether autoreactive B cells can become resting MBCs is uncertain. Further, the phenotypic identity of ABCs and their relationship to other B cell subsets, such as plasmablasts, is unclear. Whether ABCs directly promote disease is untested. Here we report, in the MRL/lpr SLE model, unexpected heterogeneity among ABC-like cells for expression of the integrins CD11b and CD11c, T-bet, and memory or plasmablast markers. Transfer and labeling studies demonstrated that ABCs are dynamic, rapidly turning over. scRNA-seq identified B cell clones present in multiple subsets, revealing that ABCs can be plasmablast precursors or undergo cycles of reactivation. Deletion of CD11c-expressing B cells revealed a direct role for ABC-like B cells in lupus pathogenesis.

Food allergy is caused by allergen-specific IgE antibodies but little is known about the B cell memory of persistent IgE responses. Here we describe in human pediatric peanut allergy CD23+IgG1+memory B cells arising in type 2 responses that contain peanut specific clones and generate IgE cells on activation. These ‘type2-marked’ IgG1+memory B cells differentially express IL-4/IL-13 regulated genesFCER2/CD23, IL4R, and germlineIGHEand carry highly mutated B cell receptors (BCRs). Further, high affinity memory B cells specific for the main peanut allergen Ara h 2 mapped to the population of ‘type2-marked’ IgG1+memory B cells and included convergent BCRs across different individuals. Our findings indicate that CD23+IgG1+memory B cells transcribing germlineIGHEare a unique memory population containing precursors of pathogenic IgE.One-Sentence SummaryWe describe a unique population of IgG+memory B cells poised to switch to IgE that contains high affinity allergen-specific clones in peanut allergy.

Multicellular life requires altruistic cooperation between cells. The adaptive immune system is a notable exception, wherein germinal center B cells compete vigorously for limiting positive selection signals. Studying primary human lymphomas and developing new mouse models, we found that mutations affecting BTG1 disrupt a critical immune gatekeeper mechanism that strictly limits B cell fitness during antibody affinity maturation. This mechanism converted germinal center B cells into supercompetitors that rapidly outstrip their normal counterparts. This effect was conferred by a small shift in MYC protein induction kinetics but resulted in aggressive invasive lymphomas, which in humans are linked to dire clinical outcomes. Our findings reveal a delicate evolutionary trade-off between natural selection of B cells to provide immunity and potentially dangerous features that recall the more competitive nature of unicellular organisms.

Most studies of adaptive immunity to SARS-CoV-2 infection focus on peripheral blood, which may not fully reflect immune responses at the site of infection. Using samples from 110 children undergoing tonsillectomy and adenoidectomy during the COVID-19 pandemic, we identified 24 samples with evidence of previous SARS-CoV-2 infection, including neutralizing antibodies in serum and SARS-CoV-2-specific germinal center and memory B cells in the tonsils and adenoids. Single-cell B cell receptor (BCR) sequencing indicated virus-specific BCRs were class-switched and somatically hypermutated, with overlapping clones in the two tissues. Expanded T cell clonotypes were found in tonsils, adenoids and blood post-COVID-19, some with CDR3 sequences identical to previously reported SARS-CoV-2-reactive T cell receptors (TCRs). Pharyngeal tissues from COVID-19-convalescent children showed persistent expansion of germinal center and antiviral lymphocyte populations associated with interferon (IFN)-γ-type responses, particularly in the adenoids, and viral RNA in both tissues. Our results provide evidence for persistent tissue-specific immunity to SARS-CoV-2 in the upper respiratory tract of children after infection.

T cell-B cell interaction is the key immune response to protect the host from severe viral infection. However, how T cells support B cells to exert protective humoral immunity in humans is not well understood. Here, we use COVID-19 as a model of acute viral infections and analyze CD4

The poor efficacy of seasonal influenza virus vaccines is often attributed to pre-existing immunity interfering with the persistence and maturation of vaccine-induced B cell responses. We previously showed that a subset of vaccine-induced B cell lineages are recruited into germinal centers (GCs) following vaccination, suggesting that affinity maturation of these lineages against vaccine antigens can occur. However, it remains to be determined whether seasonal influenza vaccination stimulates additional evolution of vaccine-specific lineages, and previous work has found no significant increase in somatic hypermutation among influenza-binding lineages sampled from the blood following seasonal vaccination in humans. Here, we investigate this issue using a phylogenetic test of measurable immunoglobulin sequence evolution. We first validate this test through simulations and survey measurable evolution across multiple conditions. We find significant heterogeneity in measurable B cell evolution across conditions, with enrichment in primary response conditions such as HIV infection and early childhood development. We then show that measurable evolution following influenza vaccination is highly compartmentalized: while lineages in the blood are rarely measurably evolving following influenza vaccination, lineages containing GC B cells are frequently measurably evolving. Many of these lineages appear to derive from memory B cells. We conclude from these findings that seasonal influenza virus vaccination can stimulate additional evolution of responding B cell lineages, and imply that the poor efficacy of seasonal influenza vaccination is not due to a complete inhibition of vaccine-specific B cell evolution., eLife digest When the immune system encounters a disease-causing pathogen, it releases antibodies that can bind to specific regions of the bacterium or virus and help to clear the infection. These proteins are generated by B cells which, upon detecting the pathogen, can begin to mutate and alter the structure of the antibody they produce: the better the antibody is at binding to the pathogen, the more likely the B cell is to survive. This process of evolution produces B cells that make more effective antibodies. After the infection, some of these cells become ‘memory B cells’ which can be stimulated in to action when the pathogen invades again. Many vaccines also depend on this process to trigger the production of memory B cells that can fight off a specific disease-causing agent. However, it is unclear to what extent memory B cells that already exist are able to continue to evolve and modify their antibodies. This is particularly important for the flu vaccine, as the virus that causes influenza rapidly mutates. To provide high levels of protection, the memory B cells formed following the vaccine may therefore need to evolve to make different antibodies that recognize mutated forms of the virus. It is thought that the low effectiveness of the flu vaccine is partially because the response it triggers does not stimulate additional evolution of memory B cells. To test this theory, Hoehn et al. developed a computational method that can detect the evolution of B cells over time. The tool was applied to samples collected from the blood and lymph nodes (organ where immune cells reside) of people who recently received the flu vaccine. The results were then compared to B cells taken from people after different infections, vaccinations, and other conditions. Hoehn et al. found the degree to which B cells evolve varies significantly between conditions. For example, B cells produced during chronic HIV infections frequently evolved over time, while such evolution was rarely observed during the autoimmune disease myasthenia gravis. The analysis also showed that memory B cells produced by the flu vaccine were able to evolve if recruited to the lymph nodes, but this was rarely detected in B cells in the blood. These findings suggest the low efficacy of the flu vaccine is not due to a complete lack of B cell evolution, but likely due to other factors. For instance, it is possible the evolutionary process it stimulates is not as robust as in other conditions, or is less likely to produce long-lived B cells that release antibodies. More research is needed to explore these ideas and could lead to the development of more effective flu vaccines.

Elevated N-linked glycosylation of immunoglobulin G variable regions (IgG-V(N-Glyc)) is an emerging molecular phenotype associated with autoimmune disorders. To test the broader specificity of elevated IgG-V(N-Glyc), we studied patients with distinct subtypes of myasthenia gravis (MG), a B cell-mediated autoimmune disease. Our experimental design focused on examining the B cell repertoire and total IgG. It specifically included adaptive immune receptor repertoire sequencing to quantify and characterize N-linked glycosylation sites in the circulating B cell receptor repertoire, proteomics to examine glycosylation patterns of the total circulating IgG, and an exploration of human-derived recombinant autoantibodies, which were studied with mass spectrometry and antigen binding assays to respectively confirm occupation of glycosylation sites and determine whether they alter binding. We found that the frequency of IgG-V(N-Glyc) motifs was increased in the total B cell receptor repertoire of MG patients when compared to healthy donors. The elevated frequency was attributed to both biased V gene segment usage and somatic hypermutation. IgG-V(N-Glyc) could be observed in the total circulating IgG in a subset of MG patients. Autoantigen binding, by four patient-derived MG autoantigen-specific monoclonal antibodies with experimentally confirmed presence of IgG-V(N-Glyc), was not altered by the glycosylation. Our findings extend prior work on patterns of immunoglobulin variable region N-linked glycosylation in autoimmunity to MG subtypes.

Diverse antibody repertoires spanning multiple lymphoid organs (e.g., 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 (e.g., 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.

BTG1 somatic mutations are exclusively found in mature B cell malignancies, ~12% diffuse large B cell lymphoma (DLBCL) and are particularly enriched in the MCD/cluster 5 subtype of ABC-DLBCL, characterized by extranodal dissemination and poor clinical outcome. However, the mechanism of action and clinical relevance of BTG1 mutations remain unknown. We find that BTG1 mutations score among the top mutations with DLBCL driver potential, using a rigorous genomic and epigenomic covariates analysis. Most notably, BTG1 mutant patients presented with inferior clinical outcome (p=0.0011) in ABC-DLBCL cases from 5 cohorts and BTG1 mutation was independently associated with lower overall survival in a multivariable Cox regression analysis (p=0.0190) DLBCL originates from mature B cells having experienced the germinal center (GC) reaction. We therefore generated a conditional knockin mouse model to express the most frequent Btg1 Q36H in B cells. Btg1 Q36H GC B cells showed a dramatic fitness advantage in in vivo competitive assays. This effect was specific to the GC compartment and was dependent on T cells. RNAseq showed that Btg1Q36H GC B cells were markedly enriched for genes normally transiently induced upon positive selection by T cells, including MYC targets. The same signatures were enriched in BTG1 mutant DLBCL patients and isogenic BTG1Q36H vs BTG1WT human DLBCL cell lines. We observed a higher proportion of MYC-expressing cells in Btg1Q36H GCs without an increase in maximal MYC levels per cell, also confirmed in human DLBCL lines and primary human tonsillar B cells, suggesting a lowered threshold for MYC induction in mutant cells. Mechanistically, ~800 transcripts associated with BTG1WT, but not BTG1Q36H. These belong to the same gene sets that characterize positively selected GC B cell, including MYC itself. BTG1Q36H therefore showed reduced association with the MYC mRNA. We found that this associated with facilitated MYC protein synthesis. Polysome profiling also showed higher translation capacity in BTG1Q36H DLBCL cells. Crossing our Btg1 Q36H mice to the VavP-Bcl2 model showed that VavP-Bcl2+Btg1Q36H mice displayed shorter survival, earlier onset of lymphoma and dysplastic B cell infiltration into non lymphoid organs. Immunoglobulin genes sequencing showed that VavP-Bcl2+Btg1Q36H lymphoma B cells were highly clonal, extensively mutated and selected. Importantly, the lack of BTG1 deletions suggested that BTG1 missense mutations do not cause a full loss-of-function of the protein. Indeed, we observed that shRNA-mediated knockdown of BTG1 resulted in apoptosis. Collectively, BTG1 mutations contribute to the formation of aggressive lymphomas through an entirely novel mechanism, by lowering the threshold to MYC induction in response to T cell selection signals, conferring dramatic fitness. This effect appears to correspond to a partial loss of function disrupting a novel GC context-specific check point, whereby BTG1 normally attenuates spurious MYC translation to tightly restrict fitness potential. Citation Format: Coraline Mlynarczyk, Matt Teater, Juhee Pae, Ling Wang, Jonatan Ersching, Antonin Papin, Darko Barisic, Ersilia Barin, Kenneth B. Hoehn, Zhengming Chen, Diu T. T. Nguyen, Chiara Evans, Ashley S. Doane, Michael G. Kharas, David W. Scott, Gabriel Victora, Ari Melnick. BTG1 mutations confer a fitness advantage and promote aggressive B cell lymphoma development by lowering the threshold for MYC induction [abstract]. In: Proceedings of the Third AACR International Meeting: Advances in Malignant Lymphoma: Maximizing the Basic-Translational Interface for Clinical Application; 2022 Jun 23-26; Boston, MA. Philadelphia (PA): AACR; Blood Cancer Discov 2022;3(5_Suppl):Abstract nr A24.

B cells undergo rapid mutation and selection for antibody binding affinity when producing antibodies capable of neutralizing pathogens. This evolutionary process can be intermixed with migration between tissues, differentiation between cellular subsets, and switching between functional isotypes. B cell receptor (BCR) sequence data has the potential to elucidate important information about these processes. However, there is currently no robust, generalizable framework for making such inferences from BCR sequence data. To address this, we develop three parsimony-based summary statistics to characterize migration, differentiation, and isotype switching along B cell phylogenetic trees. We use simulations to demonstrate the effectiveness of this approach. We then use this framework to infer patterns of cellular differentiation and isotype switching from high throughput BCR sequence datasets obtained from patients in a study of HIV infection and a study of food allergy. These methods are implemented in the R packagedowser, available athttps://bitbucket.org/kleinstein/dowser.Author summaryB cells produce high affinity antibodies through an evolutionary process of mutation and selection during adaptive immune responses. Migration between tissues, differentiation to cellular subtypes, and switching between different antibody isotypes can be important factors in shaping the role B cells play in response to infection, autoimmune disease, and allergies. B cell receptor (BCR) sequence data has the potential to elucidate important information about these processes. However, there is currently no robust, generalizable framework for making such inferences from BCR sequence data. Here, we develop three parsimony-based summary statistics to characterize migration, differentiation, and isotype switching along B cell phylogenetic trees. Using simulations, we confirm the effectiveness of our approach, as well as identify some caveats. We further use these summary statistics to investigate patterns of cellular differentiation in three HIV patients, and patterns of isotype switching in an individual with food allergies. Our methods are released in the R packagedowser:https://bitbucket.org/kleinstein/dowser.

Significance High-affinity antibodies that protect us from infection are produced by B cells through an evolutionary process of mutation and selection during adaptive immune responses. B cell repertoire sequencing combined with phylogenetic methods has provided unprecedented potential to study B cells as an evolving population. However, phylogenetic models operate on individual lineages rather than the thousands of lineages often found in B cell repertoires. Here, we develop an evolutionary framework that incorporates B cell-specific features and combines information across lineages to characterize mutation and selection dynamics of entire repertoires. We use this technique to demonstrate evidence of age-associated changes in somatic hypermutation targeting and uncover a general trend within our datasets toward negative selection over the course of affinity maturation., In order to produce effective antibodies, B cells undergo rapid somatic hypermutation (SHM) and selection for binding affinity to antigen via a process called affinity maturation. The similarities between this process and evolution by natural selection have led many groups to use phylogenetic methods to characterize the development of immunological memory, vaccination, and other processes that depend on affinity maturation. However, these applications are limited by the fact that most phylogenetic models are designed to be applied to individual lineages comprising genetically diverse sequences, while B cell repertoires often consist of hundreds to thousands of separate low-diversity lineages. Further, several features of affinity maturation violate important assumptions in standard phylogenetic models. Here, we introduce a hierarchical phylogenetic framework that integrates information from all lineages in a repertoire to more precisely estimate model parameters while simultaneously incorporating the unique features of SHM. We demonstrate the power of this repertoire-wide approach by characterizing previously undescribed phenomena in affinity maturation. First, we find evidence consistent with age-related changes in SHM hot-spot targeting. Second, we identify a consistent relationship between increased tree length and signs of increased negative selection, apparent in the repertoires of recently vaccinated subjects and those without any known recent infections or vaccinations. This suggests that B cell lineages shift toward negative selection over time as a general feature of affinity maturation. Our study provides a framework for undertaking repertoire-wide phylogenetic testing of SHM hypotheses and provides a means of characterizing dynamics of mutation and selection during affinity maturation.

Multisystem inflammatory syndrome in children (MIS-C) is a life-threatening post-infectious complication occurring unpredictably weeks after mild or asymptomatic SARS-CoV-2 infection. We profiled MIS-C, adult COVID-19, and healthy pediatric and adult individuals using single-cell RNA sequencing, flow cytometry, antigen receptor repertoire analysis, and unbiased serum proteomics, which collectively identified a signature in MIS-C patients that correlated with disease severity. Despite having no evidence of active infection, MIS-C patients had elevated S100A-family alarmins and decreased antigen presentation signatures, indicative of myeloid dysfunction. MIS-C patients showed elevated expression of cytotoxicity genes in NK and CD8+ T cells and expansion of specific IgG-expressing plasmablasts. Clinically severe MIS-C patients displayed skewed memory T cell TCR repertoires and autoimmunity characterized by endothelium-reactive IgG. The alarmin, cytotoxicity, TCR repertoire, and plasmablast signatures we defined have potential for application in the clinic to better diagnose and potentially predict disease severity early in the course of MIS-C., Graphical abstract, Multisystem inflammatory syndrome in children (MIS-C) is a life-threatening and unpredictable condition of unknown etiology. Ramaswamy et al. use peripheral blood single-cell transcriptomic profiling along with other techniques to define key innate and adaptive signatures that characterize MIS-C.

SummaryA dysregulated immune response against coronavirus-2 (SARS-CoV-2) plays a critical role in the outcome of patients with coronavirus disease 2019 (COVID-19). A significant increase in circulating plasmablasts is characteristic of COVID-19 though the underlying mechanisms and its prognostic implications are not known. Here, we demonstrate that in the acute phase of COVID-19, activated PD-1highCXCR5−CD4+ T cells, peripheral helper T cells, (Tph) are significantly increased and promote inflammatory tissue-homing plasmablasts in patients with stable but not severe COVID-19. Analysis of scRNA-seq data revealed that plasmablasts in stable patients express higher levels of tissue-homing receptors including CXCR3. The increased Tph cells exhibited “B cell help” signatures similar to that of circulating T follicular helper (cTfh) cells and promoted B cell differentiation in vitro. Compared with cTfh cells, Tph cells produced more IFNγ, inducing tissue-homing chemokine receptors on plasmablasts. Finally, expansion of activated Tph cells was correlated with the frequency of CXCR3+ plasmablasts in the acute phase of patients with stable disease. Our results demonstrate a novel role for Tph cells in acute viral immunity by inducing ectopic, antibody secreting plasmablasts.

Elevated N-linked glycosylation of immunoglobulin G variable regions (IgG-VN-Glyc) is an emerging molecular phenotype associated with autoimmune disorders. To test the broader specificity of elevated IgG-VN-Glyc, we studied patients with distinct subtypes of myasthenia gravis (MG), a B cell-mediated autoimmune disease. Our experimental design included adaptive immune receptor repertoire sequencing to quantify and characterize N-glycosylation sites in the global B cell receptor repertoire, proteomics to examine glycosylation patterns of the circulating IgG, and production of human-derived recombinant autoantibodies, which were studied with mass spectrometry and antigen binding assays to confirm occupation of glycosylation sites and determine whether they alter binding. We found that the frequency of IgG-VN-Glyc motifs was increased in the B cell repertoire of MG patients when compared to healthy donors. Motifs were introduced by both biased V gene segment usage and somatic hypermutation. IgG-VN-Glyc could be observed in the circulating IgG in a subset of MG patients. Autoantigen binding, by patient-derived MG autoantigen-specific monoclonal antibodies with experimentally confirmed presence of IgG-VN-Glyc, was not altered by the glycosylation. Our findings extend prior work on patterns of variable region N-linked glycosylation in autoimmunity to MG subtypes. Although occupied IgG-VN-Glyc motifs are found on MG autoantigen-specific monoclonal antibodies, they are not required for binding to the autoantigen in this disease.

The skin lesion erythema migrans (EM) is an initial sign of the Ixodes tick-transmitted Borreliella spirochetal infection known as Lyme disease. T cells and innate immune cells have previously been shown to predominate the EM lesion and promote the reaction. Despite the established importance of B cells and antibodies in preventing infection, the role of B cells in the skin immune response to Borreliella is unknown. Here, we used single-cell RNA-Seq in conjunction with B cell receptor (BCR) sequencing to immunophenotype EM lesions and their associated B cells and BCR repertoires. We found that B cells were more abundant in EM in comparison with autologous uninvolved skin; many were clonally expanded and had circulating relatives. EM-associated B cells upregulated the expression of MHC class II genes and exhibited preferential IgM isotype usage. A subset also exhibited low levels of somatic hypermutation despite a gene expression profile consistent with memory B cells. Our study demonstrates that single-cell gene expression with paired BCR sequencing can be used to interrogate the sparse B cell populations in human skin and reveals that B cells in the skin infection site in early Lyme disease expressed a phenotype consistent with local antigen presentation and antibody production.

Protective immunity against COVID-19 likely depends on the production of SARS-CoV-2–specific plasma cells and memory B cells postinfection or postvaccination. Previous work has found that germinal center reactions are disrupted in severe COVID-19. This may adversely affect long-term immunity against reinfection. Consistent with an extrafollicular B cell response, patients with severe COVID-19 have elevated frequencies of clonally expanded, class-switched, unmutated plasmablasts. However, it is unclear whether B cell populations in individuals with mild COVID-19 are similarly skewed. In this study, we use single-cell RNA sequencing of B cells to show that in contrast to patients with severe COVID-19, subjects with mildly symptomatic COVID-19 have B cell repertoires enriched for clonally diverse, somatically hypermutated memory B cells ∼30 d after the onset of symptoms. This provides evidence that B cell responses are less disrupted in mild COVID-19 and result in the production of memory B cells.

Poor efficacy of seasonal influenza virus vaccines is often attributed to pre-existing immunity interfering with the persistence and maturation of vaccine-induced B cell responses.1 Consistent with this notion, no significant increase in somatic hypermutation (SHM) among circulating influenza-binding lineages was detected following seasonal vaccination in humans.2 A more recent study showed that at least a subset of vaccine-induced B cell lineages are recruited into germinal centers (GCs) following vaccination, suggesting that affinity maturation of these lineages can occur.3 Crucially, however, it has not been demonstrated whether these GC-engaged lineages actually accumulate additional SHM. Here, we address this point using a phylogenetic test of measurable evolution. We first validate this test through simulations and demonstrate measurable B cell evolution in known examples of affinity maturation such as the response to HIV infection. We then show that lineages in the blood are rarely measurably evolving following influenza vaccination, but that GC-engaged lineages - likely derived from memory B cells - are frequently measurably evolving. These findings confirm that seasonal influenza virus vaccination can stimulate additional SHM among responding B cell lineages.

SUMMARYMultisystem inflammatory syndrome in children (MIS-C) is a life-threatening post-infectious complication occurring unpredictably weeks after mild or asymptomatic SARS-CoV2 infection in otherwise healthy children. Here, we define immune abnormalities in MIS-C compared to adult COVID-19 and pediatric/adult healthy controls using single-cell RNA sequencing, antigen receptor repertoire analysis, unbiased serum proteomics, andin vitroassays. Despite no evidence of active infection, we uncover elevated S100A-family alarmins in myeloid cells and marked enrichment of serum proteins that map to myeloid cells and pathways including cytokines, complement/coagulation, and fluid shear stress in MIS-C patients. Moreover, NK and CD8 T cell cytotoxicity genes are elevated, and plasmablasts harboring IgG1 and IgG3 are expanded. Consistently, we detect elevated binding of serum IgG from severe MIS-C patients to activated human cardiac microvascular endothelial cells in culture. Thus, we define immunopathology features of MIS-C with implications for predicting and managing this SARS-CoV2-induced critical illness in children.

A dysregulated immune response against the SARS-CoV-2 virus plays a critical role in severe COVID-19. However, the molecular and cellular mechanisms by which the virus causes lethal immunopathology are poorly understood. Here, we utilize multiomics single-cell analysis to probe dynamic immune responses in patients with stable or progressive manifestations of COVID-19, and assess the effects of tocilizumab, an anti-IL-6 receptor monoclonal antibody. Coordinated profiling of gene expression and cell lineage protein markers reveals a prominent type-1 interferon response across all immune cells, especially in progressive patients. An anti-inflammatory innate immune response and a pre-exhaustion phenotype in activated T cells are hallmarks of progressive disease. Skewed T cell receptor repertoires in CD8+ T cells and uniquely enriched V(D)J sequences are also identified in COVID-19 patients. B cell repertoire and somatic hypermutation analysis are consistent with a primary immune response, with possible contribution from memory B cells. Our in-depth immune profiling reveals dyssynchrony of the innate and adaptive immune interaction in progressive COVID-19, which may contribute to delayed virus clearance and has implications for therapeutic intervention.

Rituximab, a B cell-depleting therapy, is indicated for treating a growing number of autoantibody-mediated autoimmune disorders. However, relapses can occur after treatment, and autoantibody-producing B cell subsets may be found during relapses. It is not understood whether these autoantibody-producing B cell subsets emerge from the failed depletion of preexisting B cells or are generated de novo. To further define the mechanisms that cause postrituximab relapse, we studied patients with autoantibody-mediated muscle-specific kinase (MuSK) myasthenia gravis (MG) who relapsed after treatment. We carried out single-cell transcriptional and B cell receptor profiling on longitudinal B cell samples. We identified clones present before therapy that persisted during relapse. Persistent B cell clones included both antibody-secreting cells and memory B cells characterized by gene expression signatures associated with B cell survival. A subset of persistent antibody-secreting cells and memory B cells were specific for the MuSK autoantigen. These results demonstrate that rituximab is not fully effective at eliminating autoantibody-producing B cells and provide a mechanistic understanding of postrituximab relapse in MuSK MG.

Somatic missense mutations of BTG1 are exclusive to germinal center (GC)-derived B cell lymphomas (~12% of DLBCLs) and are most prevalent in ABC-DLBCL (p=0.0184 vs GCB-DLBCL), particularly in the MCD/cluster 5 subtype, which features extranodal dissemination and unfavorable outcome. However, the relevance, mechanism of action and biological contribution of BTG1 mutations have not been studied. Using a rigorous genomic covariate analysis, we identified BTG1 mutations as a top genetic driver in DLBCL. Furthermore, molecular dynamics simulations indicated that BTG1 recurrent mutations, including the most frequent Q36H, disrupted the protein structure, with likely deleterious functional consequences. To investigate the effect of BTG1 mutation in GC B cells, we generated a conditional Btg1Q36H knock in mouse crossed to the B cell specific Cd19Cre line. Surprisingly, there was no apparent phenotype in GC B cells or other B cell populations. However, placing Btg1 Q36H and WT GC B cells in competition within the same mouse through adoptive transfer revealed a dramatic competitive advantage of Btg1 Q36H cells, virtually taking over the GC reaction. To gain further insight into this striking fitness advantage, we performed RNAseq in Btg1 Q36H GCs, which showed marked enrichment for genes induced in positively selected GC B cells, including MYC targets and biosynthetic pathways. The same genes were also enriched in BTG1 mutant DLBCL patients in 2 independent cohorts. Furthermore, Btg1 Q36H GC B cells displayed greater RNA content and cell size, reflecting increased fitness. Positive selection normally triggers a brief Myc pulse in GC B cells. We therefore crossed our Btg1Q36H mice to the MycGFPprotein fusion reporter and observed higher proportion of Myc GFP+ cells in Btg1 Q36H GCs. For mechanistic studies, we generated isogenic BTG1 Q36H or BTG1 WT human DLBCL cell lines. BTG1 Q36H cells exhibited enrichment for the same positively selected GC B and MYC target genes, as well as greater RNA content and cell size. BTG1 family members were suggested to interact with RNA. Performing RNA-immunoprecipitation, we discovered that ~800 transcripts associated with BTG1 WT, but not BTG1 Q36H. Notably, these corresponded to the same positively selected GC B and MYC target genes, including MYC itself. BTG1 was shown to regulate mRNA stability in other cell types. However, BTG1 Q36H did not alter MYC mRNA stability and instead facilitated MYC protein synthesis, thus disrupting a novel GC context-specific checkpoint mechanism, whereby BTG1 normally attenuates spurious MYC translation to tightly restrict fitness potential. In GC B cells, Myc induction coincides with S phase entry, but G2/M progression requires re-entry into the proliferative dark zone. To characterize GC dynamics in vivo, we performed targeted single cell RNAseq in competing Btg1 Q36H and WT GC B cells and noted earlier and higher proportion of positively selected Btg1 Q36H GC B cells having committed to G2/M and the proliferative program. We confirmed faster S phase completion in competing Btg1 Q36H GC B cells by in vivo EdU/BrdU labelling and greater re-entry into the proliferative dark zone by in vivo antigen delivery to synchronize GC B cells at the time of positive selection. Given that MCD-DLBCLs express high levels of BCL2, we crossed our Btg1Q36H mice to the VavP-Bcl2 model. As compared to VavP-Bcl2, VavP-Bcl2+Btg1 Q36H mice displayed shorter survival (p=0.0005), earlier onset of lymphoma, dysplastic B cell infiltration into non lymphoid organs and they contained highly mutated, selected and clonal tumor B cells. Moribund VavP-Bcl2+Btg1 Q36H mice uniquely featured sheets of large, immunoblastic lymphoma cells, characteristic of ABC-DLBCLs. Most notably, examining ABC-DLBCLs from 5 independent cohorts showed inferior clinical outcome for BTG1 mutant patients (p=0.0011) and independent association of BTG1 mutation with inferior overall survival by multivariable Cox regression (p=0.0190). Collectively, we find that BTG1 mutations mediate lymphomagenesis through an entirely novel mechanism of action that recapitulates the embryonic MYC-dependent "super-competitive" phenotype originally described in Drosophila imaginal disc cells. In the GC, "super-competition" is provided by BTG1 mutation via a subtle acceleration of MYC induction and GC dynamics, conferring dramatic fitness and the potential to transform into aggressive lymphomas. Disclosures Hoehn: Prellis Biologics: Consultancy. Elemento: Janssen: Research Funding; Freenome: Consultancy, Other: Current equity holder in a privately-held company; Volastra Therapeutics: Consultancy, Other: Current equity holder, Research Funding; Owkin: Consultancy, Other: Current equity holder; Champions Oncology: Consultancy; One Three Biotech: Consultancy, Other: Current equity holder; Eli Lilly: Research Funding; AstraZeneca: Research Funding; Johnson and Johnson: Research Funding. Scott: NanoString Technologies: Patents & Royalties: Patent describing measuring the proliferation signature in MCL using gene expression profiling.; BC Cancer: Patents & Royalties: Patent describing assigning DLBCL COO by gene expression profiling--licensed to NanoString Technologies. Patent describing measuring the proliferation signature in MCL using gene expression profiling. ; AstraZeneca: Consultancy; Abbvie: Consultancy; Celgene: Consultancy; Incyte: Consultancy; Janssen: Consultancy, Research Funding; Rich/Genentech: Research Funding. Melnick: Constellation: Consultancy; Epizyme: Consultancy; Daiichi Sankyo: Research Funding; Sanofi: Research Funding; Janssen Pharmaceuticals: Research Funding; KDAC Pharma: Membership on an entity's Board of Directors or advisory committees.

Nearly all chronic human infections are associated with alterations in the memory B cell (MBC) compartment, including a large expansion of CD19(hi)T-bet(hi) MBC in the peripheral blood of HIV-infected individuals with chronic viremia. Despite their prevalence, it is unclear how these B cells arise and whether they contribute to the inefficiency of antibody-mediated immunity in chronic infectious diseases. We addressed these questions by characterizing T-bet-expressing B cells in lymph nodes (LN) and identifying a strong T-bet signature among HIV-specific MBC associated with poor immunologic outcome. Confocal microscopy and quantitative imaging revealed that T-bet(hi) B cells in LN of HIV-infected chronically viremic individuals distinctly accumulated outside germinal centers (GC), which are critical for optimal antibody responses. In single-cell analyses, LN T-bet(hi) B cells of HIV-infected individuals were almost exclusively found among CD19(hi) MBC and expressed reduced GC-homing receptors. Furthermore, HIV-specific B cells of infected individuals were enriched among LN CD19(hi)T-bet(hi) MBC and displayed a distinct transcriptome, with features similar to CD19(hi)T-bet(hi) MBC in blood and LN GC B cells (GCBC). LN CD19(hi)T-bet(hi) MBC were also related to GCBC by B cell receptor (BCR)–based phylogenetic linkage but had lower BCR mutation frequencies and reduced HIV-neutralizing capacity, consistent with diminished participation in GC-mediated affinity selection. Thus, in the setting of chronic immune activation associated with HIV viremia, failure of HIV-specific B cells to enter or remain in GC may help explain the rarity of high-affinity protective antibodies.

Rituximab, an anti-CD20 B cell-depleting therapy, is indicated to treat a growing number of autoimmune disorders. However, disease relapses can occur when the B cell compartment is re-established after treatment. To explore the mechanism of relapse, we studied patients with muscle-specific kinase (MuSK) myasthenia gravis (MG), a paradigm for B cell-mediated autoimmune diseases in which pathology is directly associated with autoantibody production. Whether the failed depletion of specific clones leads to the persistence of autoantibody producing B cell subsets has yet to be elucidated. Here we show that rituximab (RTX) therapy is not fully effective at eliminating disease-associated B cells and provide a mechanistic understanding of relapse. We carried out single-cell transcriptional and B cell receptor (BCR) profiling on longitudinal B cell samples from the peripheral blood of MuSK MG patients who relapsed after rituximab therapy. Computational analysis of these data identified individual B cells that were refractory to rituximab depletion by linking them to clones that were present prior to therapy and continued to persist after therapy. These persistent B cell clones were present at high frequency, and included both antibody-secreting cells and memory B cells. They were disproportionately isotype switched with elevated somatic hypermutation frequencies and a gene expression signature associated with antigen-experience and B cell clonal expansion. We further identified both antibody-secreting cells and memory B cells with specificity for MuSK autoantigen that were clonally expanded during relapse and persistent. Our results provide evidence that post-rituximab relapse in MuSK MG is associated with the failed depletion of autoantigen-specific B cell subsets. This work provides insight into the durability of human B cells during treatment with therapeutic CD20-specific monoclonal antibodies. The transcriptional signatures associated with B cell resistance may represent new therapeutic avenues for treatment and biomarkers of depletion efficacy.

International audience; Follicular lymphomas (FLs) are slow-growing, indolent tumors containing extensive follicular dendritic cell (FDC) networks and recurrent EZH2 gain-of-function mutations. Paradoxically, FLs originate from highly proliferative germinal center (GC) B cells with proliferation strictly dependent on interactions with T follicular helper cells. Herein, we show that EZH2 mutations initiate FL by attenuating GC B cell requirement for T cell help and driving slow expansion of GC centrocytes that become enmeshed with and dependent on FDCs. By impairing T cell help, mutant EZH2 prevents induction of proliferative MYC programs. Thus, EZH2 mutation fosters malignant transformation by epigenetically reprograming B cells to form an aberrant immunological niche that reflects characteristic features of human FLs, explaining how indolent tumors arise from GC B cells.

Summary To what extent can natural selection act on groupings above the species level? Despite extensive theoretical discussion and growing practical concerns over increased rates of global ecological turnover, the question has largely evaded empirical resolution. A flexible and robust hypothesis-testing framework for detecting the phenomenon could facilitate significant progress in resolving this issue. We introduce a permutation-based approach, implemented in the R package perspectev, which provides an explicit test of whether empirical patterns of correlation between upper level trait values and survivorship are reducible to correlations manifested at lower levels. The package is applicable to virtually any nested set of upper- and lower level groupings, a wide variety of upper level traits, and both historical and contemporary occurrence data. We apply this approach to five paleontological data sets that represent different magnitudes of extinction and differ in taxonomic breadth, geological timing and geographic extent. Using simulations, we demonstrate that this method is a robust means of detecting irreducibility in the relationship between upper level traits and survivorship, and outline circumstances in which the method is less effective. We also find evidence consistent with previous findings of selection above the species level for geographic range size in North American K-Pg molluscs and show that this phenomenon was evident for the same molluscan genera globally. Ultimately, we conclude that at certain points in history, some higher level taxonomic groups have survived differentially with respect to geographic range size in a manner that is not explained by the same trait at the species level, and we show that evidence for this phenomenon varies across taxa and extinction events. We release our method as a flexible and easy-to-use R package that will allow others to help determine the relative frequency of this macroevolutionary phenomenon, both in the fossil record and in estimates of contemporary extinction risk.

Phylogenetic methods have shown promise in understanding the development of broadly neutralizing antibody lineages (bNAbs). However, the mutational process that generates these lineages – somatic hypermutation (SHM) – is biased by hotspot motifs, which violates important assumptions in most phylogenetic substitution models. Here, we develop a modified GY94-type substitution model that partially accounts for this context-dependency while preserving independence of sites during calculation. This model shows a substantially better fit to three well-characterized bNAb lineages than the standard GY94 model. We show through simulations that accounting for hotspot motifs can lead to reduced bias of other substitution parameters, and more accurate ancestral state reconstructions. We also demonstrate how our model can be used to test hypotheses concerning the roles of different hotspot and coldspot motifs in the evolution of B-cell lineages. Further, we explore the consequences of the idea that the number of hotspot motifs – and perhaps the mutation rate in general – is expected to decay over time in individual bNAb lineages.

Phylogenetic methods have shown promise in understanding the development of broadly neutralizing antibody lineages (bNAbs). However, the mutational process that generates these lineages, somatic hypermutation, is biased by hotspot motifs which violates important assumptions in most phylogenetic substitution models. Here, we develop a modified GY94-type substitution model that partially accounts for this context dependency while preserving independence of sites during calculation. This model shows a substantially better fit to three well-characterized bNAb lineages than the standard GY94 model. We also demonstrate how our model can be used to test hypotheses concerning the roles of different hotspot and coldspot motifs in the evolution of B-cell lineages. Further, we explore the consequences of the idea that the number of hotspot motifs, and perhaps the mutation rate in general, is expected to decay over time in individual bNAb lineages.

Phylogenetic methods have shown great promise in understanding the development of broadly neutralizing antibody lineages (bNAbs). However, mutational process for generating these lineages - somatic hypermutation (SHM) - is biased by hotspot motifs, which violates important assumptions in most phylogenetic substitution models. Here, we develop a modified GY94-type substitution model which partially accounts for this context-dependency while preserving independence of sites in calculation. This model shows a substantially better fit to three well-characterized bNAb lineages than the standard GY94 model. We show through simulations that accounting for this can lead to reduced bias of other substitution parameters, and more accurate ancestral state reconstructions. We further explore other implications of this model; namely, that the number of hotspot motifs - and therefore likely the mutation rate in general - is expected to decay over time in individual bNAb lineages.

B-cell receptors (BCRs) are membrane-bound immunoglobulins that recognize and bind foreign proteins (antigens). BCRs are formed through random somatic changes of germline DNA, creating a vast repertoire of unique sequences that enable individuals to recognize a diverse range of antigens. After encountering antigen for the first time, BCRs undergo a process of affinity maturation, whereby cycles of rapid somatic mutation and selection lead to improved antigen binding. This constitutes an accelerated evolutionary process that takes place over days or weeks. Next-generation sequencing of the gene regions that determine BCR binding has begun to reveal the diversity and dynamics of BCR repertoires in unprecedented detail. Although this new type of sequence data has the potential to revolutionize our understanding of infection dynamics, quantitative analysis is complicated by the unique biology and high diversity of BCR sequences. Models and concepts from molecular evolution and phylogenetics that have been applied successfully to rapidly evolving pathogen populations are increasingly being adopted to study BCR diversity and divergence within individuals. However, BCR dynamics may violate key assumptions of many standard evolutionary methods, as they do not descend from a single ancestor, and experience biased mutation. Here, we review the application of evolutionary models to BCR repertoires and discuss the issues we believe need be addressed for this interdisciplinary field to flourish.

The scuttle fly, Megaselia scalaris, is often cited as a model in which to study early sex chromosome evolution because of its homomorphic sex chromosomes, low but measurable molecular differentiation between sex chromosomes, and occasional transposition of the male-determining element to different chromosomes in laboratory cultures. Counterintuitively, natural isolates consistently show sex linkage to the second chromosome. Frequent natural transposition of the male-determining element should lead to the loss of male specificity of any nontransposed material on the previous sex-linked chromosome pair. Using next-generation sequencing data from a newly obtained natural isolate of M. scalaris, we show that even highly conservative estimates for the size of the male-specific genome are likely too large to be contained within a transposable element. This result strongly suggests that transposition of the male-determining region either is extremely rare or has not persisted recently in natural populations, allowing for differentiation of the sex chromosomes of this species.

Premature termination codon (PTC) mutations can have dramatic effects--both adaptive and deleterious--on gene expression and function. Here, we examine the number and selective effects of PTC mutations within the Drosophila pseudoobscura subclade using 18 resequenced genomes aligned to the reference genome. We located and characterized 1,679 PTC mutations in 605 genes across each of these genomes relative to the D. pseudoobscura reference genome, and use RT-PCR to confirm transcription of a subset of these genes containing PTC mutations. We confirm previous findings that genes containing PTC mutations are less selectively constrained and less broadly expressed than non-PTC-containing genes, suggesting that the most of these mutations are at least mildly deleterious. Further, we find highly significant codon usage bias in regions downstream of the PTC in 38 of these PTC-containing genes, suggesting that some of these PTC mutations--if not alternatively spliced out of the transcript--have neutral effects. Ultimately, these analyzes support the view that the PTC mutations are mostly detrimental, but are nonetheless common enough in genomes that a subset could be effectively neutral.

Chromosomal inversions disrupt recombination in heterozygotes by both reducing crossing-over within inverted regions and increasing it elsewhere in the genome. The reduction of recombination in inverted regions facilitates the maintenance of hybridizing species, as outlined by various models of chromosomal speciation. We present a comprehensive comparison of the effects of inversions on recombination rates and on nucleotide divergence. Within an inversion differentiating Drosophila pseudoobscura and Drosophila persimilis, we detected one double recombinant among 9,739 progeny from F(1) hybrids screened, consistent with published double-crossover frequencies observed within species. Despite similar rates of exchange within and between species, we found no sequence-based evidence of ongoing gene exchange between species within this inversion, but significant exchange was inferred within species. We also observed greater differentiation at regions near inversion breakpoints between species versus within species. Moreover, we observed strong "interchromosomal effect" (higher recombination in inversion heterozygotes between species) with up to 9-fold higher recombination rates along collinear segments of chromosome two in hybrids. Further, we observed that regions most susceptible to changes in recombination rates corresponded to regions with lower recombination rates in homokaryotypes. Finally, we showed that interspecies nucleotide divergence is lower in regions with greater increases in recombination rate, potentially resulting from greater interspecies exchange. Overall, we have identified several similarities and differences between inversions segregating within versus between species in their effects on recombination and divergence. We conclude that these differences are most likely due to lower frequency of heterokaryotypes and to fitness consequences from the accumulation of various incompatibilities between species. Additionally, we have identified possible effects of inversions on interspecies gene exchange that had not been considered previously.

Duplications and deletions of large-scale genomic regions, including whole chromosomes, are a hallmark of many malignant tumors. Gross cytogenetic changes of this type were among the earliest genomic lesions to be observed in tumors and represent an important predictor of outcome for many tumor types. In multiple myeloma, a core set of structural variants are regularly used for disease prognosis (e.g., del17p, t(4;14), t(14,16)). The presence of one or more of these lesions is predictive of a more aggressive disease type and poor outcomes for patients. The application of next-generation sequencing (NGS) to cancer samples has enhanced our ability to detect copy number aberrations (CNAs) and identify driver mutations. The NGS approach can generate higher resolution maps of these mutations than traditional methods, which in turn enables faster identification of the functionally relevant genes they harbor. Most NGS-based methods for the identification of CNAs in tumors require whole genome sequencing (WGS) as input. Relatively fewer methods exist for the detection of CNAs from exome sequencing and, to our knowledge, these methods all rely on having matched normal data. However, WES is still far more common than WGS and matched normal samples are frequently not available. Therefore, we have developed a novel method for the identification of CNAs from single-sample WES data. Our method combines LOWESS smoothing & discrete-wavelet-transformation to normalize the exome coverage data with an HMM for coverage segmentation. We have applied our new method to WES data from 32 myeloma cell-lines for which aCGH copy number calls were also available for training and validation. Our method shows a high correlation with the aCGH calls (mean = 95%). We are currently applying the method to a set of 102 CD138+ tumor cell samples that were collected during patient screening on three Onyx-sponsored Phase 2 clinical trials of the proteasome inhibitor carfilzomib. We will report on the frequency of known and novel lesions in these samples, comparing to what we’ve observed in cell lines and other publically available CNA datasets. Disclosures: Degenhardt: Onyx Pharmaceuticals: Employment, Equity Ownership. Hoehn:Onyx Pharmaceuticals: Employment. Kwei:Onyx Pharmaceuticals: Employment, Equity Ownership. Kirk:Onyx Pharmaceuticals: Employment, Equity Ownership. Tuch:Onyx Pharmaceuticals: Employment, Equity Ownership.

Erratum to: J Mol Evol (2012) 75:141–150DOI 10.1007/s00239-012-9528-xOur original paper neglected to acknowledge that SteveSchaeffer and Andre Wallace genotyped and provided thePikes Peak lines of Drosophila pseudoobscura forsequencing. We thank them for this contribution andapologies for neglecting to mention it directly.