Exploring the impact of clonal definition on B-cell diversity: implications for the analysis of immune repertoires

The adaptive immune system has the extraordinary ability to produce a broad range of immunoglobulins that can bind a wide variety of antigens. During adaptive immune responses, activated B cells duplicate and undergo somatic hypermutation in their B-cell receptor (BCR) genes, resulting in clonal families of diversified B cells that can be related back to a common ancestor. Advances in high-throughput sequencing technologies have enabled the high-throughput characterization of B-cell repertoires, however, the accurate identification of clonally related BCR sequences remains a major challenge. In this study, we compare three different clone identification methods on both simulated and experimental data, and investigate their impact on the characterization of B-cell diversity. We observe that different methods lead to different clonal definitions, which affects the quantification of clonal diversity in repertoire data. Our analyses show that direct comparisons between clonal clusterings and clonal diversity of different repertoires should be avoided if different clone identification methods were used to define the clones. Despite this variability, the diversity indices inferred from the repertoires’ clonal characterization across samples show similar patterns of variation regardless of the clonal identification method used. We find the Shannon entropy to be the most robust in terms of the variability of diversity rank across samples. Our analysis also suggests that the traditional germline gene alignment-based method for clonal identification remains the most accurate when the complete information about the sequence is known, but that alignment-free methods may be preferred for shorter sequencing read lengths. We make our implementation freely available as a Python library cdiversity

EPIdemiology of Surgery-Associated Acute Kidney Injury (EPIS-AKI) : Study protocol for a multicentre, observational trial

More than 300 million surgical procedures are performed each year. Acute kidney injury (AKI) is a common complication after major surgery and is associated with adverse short-term and long-term outcomes. However, there is a large variation in the incidence of reported AKI rates. The establishment of an accurate epidemiology of surgery-associated AKI is important for healthcare policy, quality initiatives, clinical trials, as well as for improving guidelines. The objective of the Epidemiology of Surgery-associated Acute Kidney Injury (EPIS-AKI) trial is to prospectively evaluate the epidemiology of AKI after major surgery using the latest Kidney Disease: Improving Global Outcomes (KDIGO) consensus definition of AKI. EPIS-AKI is an international prospective, observational, multicentre cohort study including 10 000 patients undergoing major surgery who are subsequently admitted to the ICU or a similar high dependency unit. The primary endpoint is the incidence of AKI within 72 hours after surgery according to the KDIGO criteria. Secondary endpoints include use of renal replacement therapy (RRT), mortality during ICU and hospital stay, length of ICU and hospital stay and major adverse kidney events (combined endpoint consisting of persistent renal dysfunction, RRT and mortality) at day 90. Further, we will evaluate preoperative and intraoperative risk factors affecting the incidence of postoperative AKI. In an add-on analysis, we will assess urinary biomarkers for early detection of AKI. EPIS-AKI has been approved by the leading Ethics Committee of the Medical Council North Rhine-Westphalia, of the Westphalian Wilhelms-University Münster and the corresponding Ethics Committee at each participating site. Results will be disseminated widely and published in peer-reviewed journals, presented at conferences and used to design further AKI-related trials. Trial registration number NCT04165369

Base Excision Repair in the Immune System: Small DNA Lesions With Big Consequences

The integrity of the genome is under constant threat of environmental and endogenous agents that cause DNA damage. Endogenous damage is particularly pervasive, occurring at an estimated rate of 10,000–30,000 per cell/per day, and mostly involves chemical DNA base lesions caused by oxidation, depurination, alkylation, and deamination. The base excision repair (BER) pathway is primary responsible for removing and repairing these small base lesions that would otherwise lead to mutations or DNA breaks during replication. Next to preventing DNA mutations and damage, the BER pathway is also involved in mutagenic processes in B cells during immunoglobulin (Ig) class switch recombination (CSR) and somatic hypermutation (SHM), which are instigated by uracil (U) lesions derived from activation-induced cytidine deaminase (AID) activity. BER is required for the processing of AID-induced lesions into DNA double strand breaks (DSB) that are required for CSR, and is of pivotal importance for determining the mutagenic outcome of uracil lesions during SHM. Although uracils are generally efficiently repaired by error-free BER, this process is surprisingly error-prone at the Ig loci in proliferating B cells. Breakdown of this high-fidelity process outside of the Ig loci has been linked to mutations observed in B-cell tumors and DNA breaks and chromosomal translocations in activated B cells. Next to its role in preventing cancer, BER has also been implicated in immune tolerance. Several defects in BER components have been associated with autoimmune diseases, and animal models have shown that BER defects can cause autoimmunity in a B-cell intrinsic and extrinsic fashion. In this review we discuss the contribution of BER to genomic integrity in the context of immune receptor diversification, cancer and autoimmune diseases

Exploring the impact of clonal definition on B-cell diversity: implications for the analysis of immune repertoires

The adaptive immune system has the extraordinary ability to produce a broad range of immunoglobulins that can bind a wide variety of antigens. During adaptive immune responses, activated B cells duplicate and undergo somatic hypermutation in their B-cell receptor (BCR) genes, resulting in clonal families of diversified B cells that can be related back to a common ancestor. Advances in high-throughput sequencing technologies have enabled the high-throughput characterization of B-cell repertoires, however, the accurate identification of clonally related BCR sequences remains a major challenge. In this study, we compare three different clone identification methods on both simulated and experimental data, and investigate their impact on the characterization of B-cell diversity. We observe that different methods lead to different clonal definitions, which affects the quantification of clonal diversity in repertoire data. Our analyses show that direct comparisons between clonal clusterings and clonal diversity of different repertoires should be avoided if different clone identification methods were used to define the clones. Despite this variability, the diversity indices inferred from the repertoires’ clonal characterization across samples show similar patterns of variation regardless of the clonal identification method used. We find the Shannon entropy to be the most robust in terms of the variability of diversity rank across samples. Our analysis also suggests that the traditional germline gene alignment-based method for clonal identification remains the most accurate when the complete information about the sequence is known, but that alignment-free methods may be preferred for shorter sequencing read lengths. We make our implementation freely available as a Python library cdiversity

Convergent evolution and B-cell recirculation in germinal centers in a human lymph node

Germinal centers (GCs) play a central role in generating an effective immune response against infectious pathogens, and failures in their regulating mechanisms can lead to the development of autoimmune diseases and cancer. Although previous works study experimental systems of the immune response with mouse models that are immunized with specific antigens, our study focused on a real-life situation, with an ongoing GC response in a human lymph node (LN) involving multiple asynchronized GCs reacting simultaneously to unknown antigens. We combined laser capture microdissection of individual GCs from human LN with next-generation repertoire sequencing to characterize individual GCs as distinct evolutionary spaces. In line with well-characterized GC responses in mice, elicited by immunization with model antigens, we observe a heterogeneous clonal diversity across individual GCs from the same human LN. Still, we identify shared clones in several individual GCs, and phylogenetic tree analysis combined with paratope modeling suggest the re-engagement and rediversification of B-cell clones across GCs and expanded clones exhibiting shared antigen responses across distinct GCs, indicating convergent evolution of the GCs.ISSN:2575-107

Presentation_1_Exploring the impact of clonal definition on B-cell diversity: implications for the analysis of immune repertoires.pdf

The adaptive immune system has the extraordinary ability to produce a broad range of immunoglobulins that can bind a wide variety of antigens. During adaptive immune responses, activated B cells duplicate and undergo somatic hypermutation in their B-cell receptor (BCR) genes, resulting in clonal families of diversified B cells that can be related back to a common ancestor. Advances in high-throughput sequencing technologies have enabled the high-throughput characterization of B-cell repertoires, however, the accurate identification of clonally related BCR sequences remains a major challenge. In this study, we compare three different clone identification methods on both simulated and experimental data, and investigate their impact on the characterization of B-cell diversity. We observe that different methods lead to different clonal definitions, which affects the quantification of clonal diversity in repertoire data. Our analyses show that direct comparisons between clonal clusterings and clonal diversity of different repertoires should be avoided if different clone identification methods were used to define the clones. Despite this variability, the diversity indices inferred from the repertoires’ clonal characterization across samples show similar patterns of variation regardless of the clonal identification method used. We find the Shannon entropy to be the most robust in terms of the variability of diversity rank across samples. Our analysis also suggests that the traditional germline gene alignment-based method for clonal identification remains the most accurate when the complete information about the sequence is known, but that alignment-free methods may be preferred for shorter sequencing read lengths. We make our implementation freely available as a Python library cdiversity.</p

Systematic evaluation of B-cell clonal family inference approaches

Abstract The reconstruction of clonal families (CFs) in B-cell receptor (BCR) repertoire analysis is a crucial step to understand the adaptive immune system and how it responds to antigens. The BCR repertoire of an individual is formed throughout life and is diverse due to several factors such as gene recombination and somatic hypermutation. The use of Adaptive Immune Receptor Repertoire sequencing (AIRR-seq) using next generation sequencing enabled the generation of full BCR repertoires that also include rare CFs. The reconstruction of CFs from AIRR-seq data is challenging and several approaches have been developed to solve this problem. Currently, most methods use the heavy chain (HC) only, as it is more variable than the light chain (LC). CF reconstruction options include the definition of appropriate sequence similarity measures, the use of shared mutations among sequences, and the possibility of reconstruction without preliminary clustering based on V- and J-gene annotation. In this study, we aimed to systematically evaluate different approaches for CF reconstruction and to determine their impact on various outcome measures such as the number of CFs derived, the size of the CFs, and the accuracy of the reconstruction. The methods were compared to each other and to a method that groups sequences based on identical junction sequences and another method that only determines subclones. We found that after accounting for data set variability, in particular sequencing depth and mutation load, the reconstruction approach has an impact on part of the outcome measures, including the number of CFs. Simulations indicate that unique junctions and subclones should not be used as substitutes for CF and that more complex methods do not outperform simpler methods. Also, we conclude that different approaches differ in their ability to correctly reconstruct CFs when not considering the LC and to identify shared CFs. The results showed the effect of different approaches on the reconstruction of CFs and highlighted the importance of choosing an appropriate method

Additional file 1 of Systematic evaluation of B-cell clonal family inference approaches

Additional file 1: Supplementary Figure 1. Data simulation pipeline. Simulation approach is an integration of ImmuneSim, Alakazam and SHazaM tools and equally use the data of CF groupings obtained from each of the 10 CF inference approaches. Supplementary Figure 2. Determination of the number of TP, TN, FP, and FN. Three simulated CFs (2 singletons) and two inferred CFs are shown. Supplementary Figure 3. Overall correlation between the log10(number of CFs) and the standardized sequence depth for all combinations of approach (except SCOPer; A7, A8) and dataset. Supplementary Figure 4. Overall trend between the log10(number of CFs) and the standardized mutation load for all combinations of approach (except SCOPer; A7, A8) and dataset. Supplementary Figure 5. Summary of significant pairwise comparisons between Approaches. Supplementary Figure 6. Number of TP, TN, FP, and FN cases produced by the ten approaches when applied to six samples from three simulated datasets (D10, D11, D12). Supplementary Figure 7. Normalized number of TP, TN, FP, and FN cases produced by the ten approaches when applied to six samples from three simulated datasets (D10, D11, D12)

Understanding repertoire sequencing data through a multiscale computational model of the germinal center

Sequencing of B-cell and T-cell immune receptor repertoires helps us to understand the adaptive immune response, although it only provides information about the clonotypes (lineages) and their frequencies and not about, for example, their affinity or antigen (Ag) specificity. To further characterize the identified clones, usually with special attention to the particularly abundant ones (dominant), additional time-consuming or expensive experiments are generally required. Here, we present an extension of a multiscale model of the germinal center (GC) that we previously developed to gain more insight in B-cell repertoires. We compare the extent that these simulated repertoires deviate from experimental repertoires established from single GCs, blood, or tissue. Our simulations show that there is a limited correlation between clonal abundance and affinity and that there is large affinity variability among same-ancestor (same-clone) subclones. Our simulations suggest that low-abundance clones and subclones, might also be of interest since they may have high affinity for the Ag. We show that the fraction of plasma cells (PCs) with high B-cell receptor (BcR) mRNA content in the GC does not significantly affect the number of dominant clones derived from single GCs by sequencing BcR mRNAs. Results from these simulations guide data interpretation and the design of follow-up experiments