Global characterization of B cell receptor repertoire in COVID-19 patients by single-cell V(D)J sequencing.

The world is facing a pandemic of Corona Virus Disease 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Adaptive immune responses are essential for SARS-CoV-2 virus clearance. Although a large body of studies have been conducted to investigate the immune mechanism in COVID-19 patients, we still lack a comprehensive understanding of the BCR repertoire in patients. In this study, we used the single-cell V(D)J sequencing to characterize the BCR repertoire across convalescent COVID-19 patients. We observed that the BCR diversity was significantly reduced in disease compared with healthy controls. And BCRs tend to skew toward different V gene segments in COVID-19 and healthy controls. The CDR3 sequences of heavy chain in clonal BCRs in patients were more convergent than that in healthy controls. In addition, we discovered increased IgG and IgA isotypes in the disease, including IgG1, IgG3 and IgA1. In all clonal BCRs, IgG isotypes had the most frequent class switch recombination events and the highest somatic hypermutation rate, especially IgG3. Moreover, we found that an IgG3 cluster from different clonal groups had the same IGHV, IGHJ and CDR3 sequences (IGHV4-4-CARLANTNQFYDSSSYLNAMDVW-IGHJ6). Overall, our study provides a comprehensive characterization of the BCR repertoire in COVID-19 patients, which contributes to the understanding of the mechanism for the immune response to SARS-CoV-2 infection.

Non-coding Class Switch Recombination-Related Transcription in Human Normal and Pathological Immune Responses.

Antibody class switch recombination (CSR) to IgG, IgA, or IgE is a hallmark of adaptive immunity, allowing antibody function diversification beyond IgM. CSR involves a deletion of the IgM/IgD constant region genes placing a new acceptor Constant gene, downstream of the VDJH exon. CSR depends on non-coding (CSRnc) transcription of donor Iµ and acceptor IH exons, located 5' upstream of each CH coding gene. Although, our knowledge of the role of CSRnc transcription has advanced greatly, its extension and importance in healthy and diseased humans is scarce. We analyzed CSRnc transcription in 70,603 publicly available RNA-seq samples, including GTEx, TCGA, and the Sequence Read Archive using recount2, an online resource consisting of normalized RNA-seq gene and exon counts, as well as, coverage BigWig files that can be programmatically accessed through R. CSRnc transcription was validated with a qRT-PCR assay for Iµ, Iγ3, and Iγ1 in humans in response to vaccination. We mapped IH transcription for the human IGH locus, including the less understood IGHD gene. CSRnc transcription was restricted to B cells and is widely distributed in normal adult tissues, but predominant in blood, spleen, MALT-containing tissues, visceral adipose tissue and some so-called "immune privileged" tissues. However, significant Iγ4 expression was found even in non-lymphoid fetal tissues. CSRnc expression in cancer tissues mimicked the expression of their normal counterparts, with notable pattern changes in some common cancer subsets. CSRnc transcription in tumors appears to result from tumor infiltration by B cells, since CSRnc transcription was not detected in corresponding tumor-derived immortal cell lines. Additionally, significantly increased Iδ transcription in ileal mucosa in Crohn's disease with ulceration was found. In conclusion, CSRnc transcription occurs in multiple anatomical locations beyond classical secondary lymphoid organs, representing a potentially useful marker of effector B cell responses in normal and pathological immune responses. The pattern of IH exon expression may reveal clues of the local immune response (i.e., cytokine milieu) in health and disease. This is a great example of how the public recount2 data can be used to further our understanding of transcription, including regions outside the known transcriptome.

VDJServer: A Cloud-Based Analysis Portal and Data Commons for Immune Repertoire Sequences and Rearrangements.

Background: Recent technological advances in immune repertoire sequencing have created tremendous potential for advancing our understanding of adaptive immune response dynamics in various states of health and disease. Immune repertoire sequencing produces large, highly complex data sets, however, which require specialized methods and software tools for their effective analysis and interpretation. Results: VDJServer is a cloud-based analysis portal for immune repertoire sequence data that provide access to a suite of tools for a complete analysis workflow, including modules for preprocessing and quality control of sequence reads, V(D)J gene segment assignment, repertoire characterization, and repertoire comparison. VDJServer also provides sophisticated visualizations for exploratory analysis. It is accessible through a standard web browser via a graphical user interface designed for use by immunologists, clinicians, and bioinformatics researchers. VDJServer provides a data commons for public sharing of repertoire sequencing data, as well as private sharing of data between users. We describe the main functionality and architecture of VDJServer and demonstrate its capabilities with use cases from cancer immunology and autoimmunity. Conclusion: VDJServer provides a complete analysis suite for human and mouse T-cell and B-cell receptor repertoire sequencing data. The combination of its user-friendly interface and high-performance computing allows large immune repertoire sequencing projects to be analyzed with no programming or software installation required. VDJServer is a web-accessible cloud platform that provides access through a graphical user interface to a data management infrastructure, a collection of analysis tools covering all steps in an analysis, and an infrastructure for sharing data along with workflows, results, and computational provenance. VDJServer is a free, publicly available, and open-source licensed resource.

Insights into immune system development and function from mouse T-cell repertoires.

The ability of the adaptive immune system to respond to arbitrary pathogens stems from the broad diversity of immune cell surface receptors. This diversity originates in a stochastic DNA editing process (VDJ recombination) that acts on the surface receptor gene each time a new immune cell is created from a stem cell. By analyzing T-cell receptor (TCR) sequence repertoires taken from the blood and thymus of mice of different ages, we quantify the changes in the VDJ recombination process that occur from embryo to young adult. We find a rapid increase with age in the number of random insertions and a dramatic increase in diversity. Because the blood accumulates thymic output over time, blood repertoires are mixtures of different statistical recombination processes, and we unravel the mixture statistics to obtain a picture of the time evolution of the early immune system. Sequence repertoire analysis also allows us to detect the statistical impact of selection on the output of the VDJ recombination process. The effects we find are nearly identical between thymus and blood, suggesting that our analysis mainly detects selection for proper folding of the TCR receptor protein. We further find that selection is weaker in laboratory mice than in humans and it does not affect the diversity of the repertoire.

IMonitor: A Robust Pipeline for TCR and BCR Repertoire Analysis.

The advance of next generation sequencing (NGS) techniques provides an unprecedented opportunity to probe the enormous diversity of the immune repertoire by deep sequencing T-cell receptors (TCRs) and B-cell receptors (BCRs). However, an efficient and accurate analytical tool is still on demand to process the huge amount of data. We have developed a high-resolution analytical pipeline, Immune Monitor ("IMonitor") to tackle this task. This method utilizes realignment to identify V(D)J genes and alleles after common local alignment. We compare IMonitor with other published tools by simulated and public rearranged sequences, and it demonstrates its superior performance in most aspects. Together with this, a methodology is developed to correct the PCR and sequencing errors and to minimize the PCR bias among various rearranged sequences with different V and J gene families. IMonitor provides general adaptation for sequences from all receptor chains of different species and outputs useful statistics and visualizations. In the final part of this article, we demonstrate its application on minimal residual disease detection in patients with B-cell acute lymphoblastic leukemia. In summary, this package would be of widespread usage for immune repertoire analysis.

Reversing Tolerance in Isotype Switch-Competent Anti-Insulin B Lymphocytes.

Autoreactive B lymphocytes that escape central tolerance and mature in the periphery are a liability for developing autoimmunity. IgG insulin autoantibodies that predict type 1 diabetes and complicate insulin therapies indicate that mechanisms for tolerance to insulin are flawed. To examine peripheral tolerance in anti-insulin B cells, we generated C57BL/6 mice that harbor anti-insulin VDJH-125 site directed to the native IgH locus (VH125(SD)). Class switch-competent anti-insulin B cells fail to produce IgG Abs following T cell-dependent immunization of VH125(SD) mice with heterologous insulin, and they exhibit markedly impaired proliferation to anti-CD40 plus insulin in vitro. In contrast, costimulation with LPS plus insulin drives robust anti-insulin B cell proliferation. Furthermore, VH125(SD) mice produce both IgM and IgG2a anti-insulin Abs following immunization with insulin conjugated to type 1 T cell-independent Brucella abortus ring test Ag (BRT). Anti-insulin B cells undergo clonal expansion in vivo and emerge as IgM(+) and IgM(-) GL7(+)Fas(+) germinal center (GC) B cells following immunization with insulin-BRT, but not BRT alone. Analysis of Igκ genes in VH125(SD) mice immunized with insulin-BRT reveals that anti-insulin Vκ from the preimmune repertoire is selected into GCs. These data demonstrate that class switch-competent anti-insulin B cells remain functionally silent in T cell-dependent immune responses, yet these B cells are vulnerable to reversal of anergy following combined BCR/TLR engagement that promotes Ag-specific GC responses and Ab production. Environmental factors that lead to infection and inflammation could play a critical yet underappreciated role in driving loss of tolerance and promoting autoimmune disease.

Somatic hypermutation leads to diversification of the heavy chain immunoglobulin repertoire in cattle.

The availability of unique variable (VH), diversity (D), and joining (JH) gene segments in the vertebrate germline determines the extent to which a primary immunoglobulin (Ig) repertoire can be generated through combinatorial rearrangement. Although bovine D segments possess unusual properties, the diversity of the primary Ig heavy chain (IgH) repertoire in cattle is restricted by the dominance of a single family of germline VH genes of limited number and diversity. Cattle therefore must employ other diversification strategies in order to generate a functional IgH repertoire, the main candidates being gene conversion and somatic hypermutation. In considering these possibilities, we predicted that if somatic hypermutation was active during B lymphocyte development, the process would introduce nucleotide substitutions to the VDJ exon and also non-coding region lying downstream of the rearranged JH segment. In contrast, our expectation was that gene conversion would show a greater tendency to confine modification to the IgH coding sequence, leaving intron regions substantially unmodified. An analysis of rearranged IgH sequences from cattle of different ages revealed that the diversification of germline sequences could be observed in very young calves and that substitution frequency increased with age. The age-dependent accumulation of mutations was particularly apparent in the second IgH complementarity-determining region (CDR2). Single base substitutions were found to predominate, with purines targeted more frequently than pyrimidines and transitions favoured over transversions. In non-coding regions, mutations were detected at a normalised frequency that was indistinguishable from that observed in CDR2. These data are consistent with a process of IgH diversification driven predominantly by somatic hypermutation.

Selection of T-cell receptors with a recurrent CDR3ß peptide-contact motif within the repertoire of alloreactive CD8(+) T cells.

Peptide/MHC complexes recognized by alloreactive T lymphocytes (TLs) have been identified, but their contribution to in vivo allo-rejection is not known. We previously characterized the peptide pBM1, highly represented among endogenous H-2K(b) (K(b) )-associated peptides and critically required to induce full activation of H-2(k) monoclonal CD8(+) TLs expressing the cognate TCR-BM3.3. Here, we asked whether a pBM1/K(b) -specific TL subset could be detected within a polyclonal TL population rejecting allogeneic cells in vivo. We show that the proportion of pBM1/K(b) -binding CD8(+) TLs increased from <0.04% in naïve mice to 3% of activated CD44(+) CD8(+) TLs in H-2(k) mice rejecting K(b) -expressing cells. Among these, TCR-Vß2 usage was greatly enriched, and 75% of them shared a TCR-Vß2 CDR3ß motif with the prototype TCR-BM3.3. Fewer than 5% of K(b) -reactive CD44(+) CD8(+) TLs not binding pBM1/K(b) displayed this CDR3ß motif. We found that the recurrent CDR3ß motif of pBM1/K(b) -binding TLs was assembled from distinct V/D/J recombination events, suggesting that it is recruited upon immunization for its optimal TCR-peptide/MHC fit. Thus, a CDR3ß motif generated by a process akin to "convergent recombination" accounts for a sizable fraction of the alloreactive anti-K(b) TCR repertoire.

Biased dA/dT somatic hypermutation as regulated by the heavy chain intronic iEmu enhancer and 3'Ealpha enhancers in human lymphoblastoid B cells.

Somatic hypermutation (SHM) in immunoglobulin gene (Ig) variable (V) regions is critical for the maturation of the antibody response. It is dependent on the expression of activation-induced cytidine deaminase (AID) and translesion DNA polymerases in germinal center B cells as well as Ig V transcription, as regulated by the Ig heavy chain (H) intronic enhancer (iEmu) and the 3' enhancer (3'Ealpha) region. We analyzed the role of these cis elements in SHM by stably transfecting Ramos human lymphoblastoid B cells with a rearranged human IgH chain VD (diversity) J (joining) DNA construct containing a V(H) promoter at the 5' end and C(H)1 and C(H)2 exons of Cgamma1 at the 3' end. In this construct, mutations preferentially targeted dA/dT basepairs in the RGYW/WRCY hotspot. Most of the dA/dT mutations and accompanying dC/dG mutations were transitions. Deletion of iEmu resulted in decreased SHM which could be partially restored by insertion of the IgH hs1,2 enhancer. Other two 3'Ealpha enhancers, hs3-hs4, did not significantly increase the mutation frequency, but further strengthened the dA/dT bias. The frequency and spectrum of the mutations were independent of the genomic integration of the transgene or V gene transcription level. Thus, we have established a novel in vitro system to analyze SHM and identify the role of multiple cis-regulatory elements in regulating dA/dT biased SHM. This model system will be useful to further address the role of other cis-regulating elements and recruited trans-acting factors in expressing the modalities of SHM.