Antigen Receptor Galaxy: A User-Friendly, Web-Based Tool for Analysis and Visualization of T and B Cell Receptor Repertoire Data.

Antigen Receptor Galaxy (ARGalaxy) is a Web-based tool for analyses and visualization of TCR and BCR sequencing data of 13 species. ARGalaxy consists of four parts: the demultiplex tool, the international ImMunoGeneTics information system (IMGT) concatenate tool, the immune repertoire pipeline, and the somatic hypermutation (SHM) and class switch recombination (CSR) pipeline. Together they allow the analysis of all different aspects of the immune repertoire. All pipelines can be run independently or combined, depending on the available data and the question of interest. The demultiplex tool allows data trimming and demultiplexing, whereas with the concatenate tool multiple IMGT/HighV-QUEST output files can be merged into a single file. The immune repertoire pipeline is an extended version of our previously published ImmunoGlobulin Galaxy (IGGalaxy) virtual machine that was developed to visualize V(D)J gene usage. It allows analysis of both BCR and TCR rearrangements, visualizes CDR3 characteristics (length and amino acid usage) and junction characteristics, and calculates the diversity of the immune repertoire. Finally, ARGalaxy includes the newly developed SHM and CSR pipeline to analyze SHM and/or CSR in BCR rearrangements. It analyzes the frequency and patterns of SHM, Ag selection (including BASELINe), clonality (Change-O), and CSR. The functionality of the ARGalaxy tool is illustrated in several clinical examples of patients with primary immunodeficiencies. In conclusion, ARGalaxy is a novel tool for the analysis of the complete immune repertoire, which is applicable to many patient groups with disturbances in the immune repertoire such as autoimmune diseases, allergy, and leukemia, but it can also be used to address basic research questions in repertoire formation and selection.

XLF deficiency results in reduced N-nucleotide addition during V(D)J recombination.

Repair of DNA double-strand breaks (DSBs) by the nonhomologous end-joining pathway (NHEJ) is important not only for repair of spontaneous breaks but also for breaks induced in developing lymphocytes during V(D)J (variable [V], diversity [D], and joining [J] genes) recombination of their antigen receptor loci to create a diverse repertoire. Mutations in the NHEJ factor XLF result in extreme sensitivity for ionizing radiation, microcephaly, and growth retardation comparable to mutations in LIG4 and XRCC4, which together form the NHEJ ligation complex. However, the effect on the immune system is variable (mild to severe immunodeficiency) and less prominent than that seen in deficiencies of NHEJ factors ARTEMIS and DNA-dependent protein kinase catalytic subunit, with defects in the hairpin opening step, which is crucial and unique for V(D)J recombination. Therefore, we aimed to study the role of XLF during V(D)J recombination. We obtained clinical data from 9 XLF-deficient patients and performed immune phenotyping and antigen receptor repertoire analysis of immunoglobulin (Ig) and T-cell receptor (TR) rearrangements, using next-generation sequencing in 6 patients. The results were compared with XRCC4 and LIG4 deficiency. Both Ig and TR rearrangements showed a significant decrease in the number of nontemplated (N) nucleotides inserted by terminal deoxynucleotidyl transferase, which resulted in a decrease of 2 to 3 amino acids in the CDR3. Such a reduction in the number of N-nucleotides has a great effect on the junctional diversity, and thereby on the total diversity of the Ig and TR repertoire. This shows that XLF has an important role during V(D)J recombination in creating diversity of the repertoire by stimulating N-nucleotide insertion.

Rapid Evolution of the CD8+ TCR Repertoire in Neonatal Mice.

Currently, there is little consensus regarding the most appropriate animal model to study acute infection and the virus-specific CD8(+) T cell (CTL) responses in neonates. TCRß high-throughput sequencing in naive CTL of differently aged neonatal mice was performed, which demonstrated differential Vß family gene usage. Using an acute influenza infection model, we examined the TCR repertoire of the CTL response in neonatal and adult mice infected with influenza type A virus. Three-day-old mice mounted a greatly reduced primary NP(366-374)-specific CTL response when compared with 7-d-old and adult mice, whereas secondary CTL responses were normal. Analysis of NP(366-374)-specific CTL TCR repertoire revealed different Vß gene usage and greatly reduced public clonotypes in 3-d-old neonates. This could underlie the impaired CTL response in these neonates. To directly test this, we examined whether controlling the TCR would restore neonatal CTL responses. We performed adoptive transfers of both nontransgenic and TCR-transgenic OVA(257-264)-specific (OT-I) CD8(+) T cells into influenza-infected hosts, which revealed that naive neonatal and adult OT-I cells expand equally well in neonatal and adult hosts. In contrast, nontransgenic neonatal CD8(+) T cells when transferred into adults failed to expand. We further demonstrate that differences in TCR avidity may contribute to decreased expansion of the endogenous neonatal CTL. These studies highlight the rapid evolution of the neonatal TCR repertoire during the first week of life and show that impaired neonatal CTL immunity results from an immature TCR repertoire, rather than intrinsic signaling defects or a suppressive environment.

Antigen receptor sequencing of paired bone marrow samples shows homogeneous distribution of acute lymphoblastic leukemia subclones.

In B-cell precursor acute lymphoblastic leukemia, the initial leukemic cells share the same antigen receptor gene rearrangements. However, due to ongoing rearrangement processes, leukemic cells with different gene rearrangement patterns can develop, resulting in subclone formation. We studied leukemic subclones and their distribution in the bone marrow and peripheral blood at diagnosis. Antigen receptor gene rearrangements (IGH, IGK, TRG, TRD, TRB) were analyzed by next-generation sequencing in seven paired bone marrow samples and five paired bone marrow-peripheral blood samples. Background-thresholds were defined, which enabled identification of leukemic gene rearrangements down to very low levels. Paired bone marrow analysis showed oligoclonality in all 7 patients and up to 34 leukemic clones per patient. Additional analysis of evolutionary-related IGH gene rearrangements revealed up to 171 leukemic clones per patient. Interestingly, overall 86% of all leukemic gene rearrangements, including small subclones, were present in both bone marrow samples (range per patient: 72-100%). Paired bone marrow-peripheral blood analysis showed that 83% of all leukemic gene rearrangements in bone marrow were also found in peripheral blood (range per patient: 81-100%). Remarkably, in the paired bone marrow samples and paired bone marrow-peripheral blood samples the vast majority of leukemic gene rearrangements had a similar frequency (<5-fold frequency difference) (96% and 96% of leukemic rearrangements, respectively). Together, these results indicate that B-cell precursor acute lymphoblastic leukemia is generally highly oligoclonal. Nevertheless, the vast majority of leukemic clones, even the minor antigen receptor-defined subclones, are homogeneously distributed throughout the bone marrow and peripheral blood compartment.

Identification of Associated Genes and Diseases in Patients With Congenital Upper-Limb Anomalies: A Novel Application of the OMT Classification.

PURPOSE: Congenital upper-limb anomalies (CULA) can present as a part of a syndrome or association. There is a wide spectrum of CULA, each of which might be related to different diseases. The structure provided by the Oberg, Manske, and Tonkin (OMT) classification could aid in differential diagnosis formulation in patients with CULA. The aims of this study were to review the Human Phenotype Ontology (HPO) project database for diseases and causative genes related to the CULA described in the OMT classification and to develop a methodology for differential diagnosis formulation based on the observed congenital anomalies, CulaPhen. METHODS: We reviewed the HPO database for all diseases, including causative genes related to CULA. All CULA were classified according to the OMT classification; associated non-hand phenotypes were classified into 12 anatomical groups. We analyzed the contribution of each anatomical group to a given disease and developed a tool for differential diagnosis formulation based on these contributions. We compared our results with cases from the literature and with a current HPO tool, Phenomizer. RESULTS: In total, 514 hand phenotypes were obtained, 384 of which could be classified in the OMT classification. A total of 1,403 diseases could be related to those CULA. A comparison with 10 recently published cases with CULA revealed that the presented phenotype matched the descriptions in our dataset. The differential diagnosis produced using our methodology was more accurate than Phenomizer in 4 of 5 examples. CONCLUSIONS: The OMT classification can be used to describe hand anomalies that may present in over 1,400 diseases. CulaPhen was developed to provide a (hand) phenotype-based differential diagnosis. Differential diagnosis formulation based on the proposed system outperforms the system in current use. CLINICAL RELEVANCE: This study illustrates that the OMT diagnoses, either individually or combined, can be cross-referenced with different diseases and syndromes. Therefore, use of the OMT classification can aid differential diagnosis formulation for CULA patients.

ImmunoGlobulin galaxy (IGGalaxy) for simple determination and quantitation of immunoglobulin heavy chain rearrangements from NGS.

BACKGROUND: Sequence analysis of immunoglobulin heavy chain (IGH) gene rearrangements and frequency analysis is a powerful tool for studying the immune repertoire, immune responses and immune dysregulation in health and disease. The challenge is to provide user friendly, secure and reproducible analytical services that are available for both small and large laboratories which are determining VDJ repertoire using NGS technology. RESULTS: In this study we describe ImmunoGlobulin Galaxy (IGGalaxy)- a convenient web based application for analyzing next-generation sequencing results and reporting IGH gene rearrangements for both repertoire and clonality studies. IGGalaxy has two analysis options one using the built in igBLAST algorithm and the second using output from IMGT; in either case repertoire summaries for the B-cell populations tested are available. IGGalaxy supports multi-sample and multi-replicate input analysis for both igBLAST and IMGT/HIGHV-QUEST. We demonstrate the technical validity of this platform using a standard dataset, S22, used for benchmarking the performance of antibody alignment utilities with a 99.9 % concordance with previous results. Re-analysis of NGS data from our samples of RAG-deficient patients demonstrated the validity and user friendliness of this tool. CONCLUSIONS: IGGalaxy provides clinical researchers with detailed insight into the repertoire of the B-cell population per individual sequenced and between control and pathogenic genomes. IGGalaxy was developed for 454 NGS results but is capable of analyzing alternative NGS data (e.g. Illumina, Ion Torrent). We demonstrate the use of a Galaxy virtual machine to determine the VDJ repertoire for reference data and from B-cells taken from immune deficient patients. IGGalaxy is available as a VM for download and use on a desktop PC or on a server.

Deciphering imprints of impaired memory B-cell maturation in germinal centers of three patients with common variable immunodeficiency.

Common variable immunodeficiency (CVID), characterized by recurrent infections, low serum class-switched immunoglobulin isotypes, and poor antigen-specific antibody responses, comprises a heterogeneous patient population in terms of clinical presentation and underlying etiology. The diagnosis is regularly associated with a severe decrease of germinal center (GC)-derived B-cell populations in peripheral blood. However, data from B-cell differentiation within GC is limited. We present a multiplex approach combining histology, flow cytometry, and B-cell receptor repertoire analysis of sorted GC B-cell populations allowing the modeling of distinct disturbances in GCs of three CVID patients. Our results reflect pathophysiological heterogeneity underlying the reduced circulating pool of post-GC memory B cells and plasmablasts in the three patients. In patient 1, quantitative and qualitative B-cell development in GCs is relatively normal. In patient 2, irregularly shaped GCs are associated with reduced somatic hypermutation (SHM), antigen selection, and class-switching, while in patient 3, high SHM, impaired antigen selection, and class-switching with large single clones imply increased re-cycling of cells within the irregularly shaped GCs. In the lymph nodes of patients 2 and 3, only limited numbers of memory B cells and plasma cells are formed. While reduced numbers of circulating post GC B cells are a general phenomenon in CVID, the integrated approach exemplified distinct defects during GC maturation ranging from near normal morphology and function to severe disturbances with different facets of impaired maturation of memory B cells and/or plasma cells. Integrated dissection of disturbed GC B-cell maturation by histology, flow cytometry, and BCR repertoire analysis contributes to unraveling defects in the essential steps during memory formation.

Next-generation antigen receptor sequencing of paired diagnosis and relapse samples of B-cell acute lymphoblastic leukemia: Clonal evolution and implications for minimal residual disease target selection.

Antigen receptor gene rearrangements are frequently applied as molecular targets for detection of minimal residual disease (MRD) in B-cell precursor acute lymphoblastic leukemia patients. Since such targets may be lost at relapse, appropriate selection of antigen receptor genes as MRD-PCR target is critical. Recently, next-generation sequencing (NGS) - much more sensitive and quantitative than classical PCR-heteroduplex approaches - has been introduced for identification of MRD-PCR targets. We evaluated 42 paired diagnosis-relapse samples by NGS (IGH, IGK, TRG, TRD, and TRB) to evaluate clonal evolution patterns and to design an algorithm for selection of antigen receptor gene rearrangements most likely to remain stable at relapse. Overall, only 393 out of 1446 (27%) clonal rearrangements were stable between diagnosis and relapse. If only index clones with a frequency >5% at diagnosis were taken into account, this number increased to 65%; including only index clones with an absolute read count >10,000, indicating truly major clones, further increased the stability to 84%. Over 90% of index clones at relapse were also present as index clone at diagnosis. Our data provide detailed information about the stability of antigen receptor gene rearrangements, based on which we propose an algorithm for selecting stable MRD-PCR targets, successful in >97% of patients.

KREAP: an automated Galaxy platform to quantify in vitro re-epithelialization kinetics.

Background: In vitro scratch assays have been widely used to study the influence of bioactive substances on the processes of cell migration and proliferation that are involved in re-epithelialization. The development of high-throughput microscopy and image analysis has enabled scratch assays to become compatible with high-throughput research. However, effective processing and in-depth analysis of such high-throughput image datasets are far from trivial and require integration of multiple image processing and data extraction software tools. Findings: We developed and implemented a kinetic re-epithelialization analysis pipeline (KREAP) in Galaxy. The KREAP toolbox incorporates freely available image analysis tools and automatically performs image segmentation and feature extraction of each image series, followed by automatic quantification of cells inside and outside the scratched area over time. The enumeration of infiltrating cells over time is modeled to extract three biologically relevant parameters that describe re-epithelialization kinetics. The output of the tools is organized, displayed, and saved in the Galaxy environment for future reference. Conclusions: The KREAP toolbox in Galaxy provides an open-source, easy-to-use, web-based platform for reproducible image processing and data analysis of high-throughput scratch assays. The KREAP toolbox could assist a broad scientific community in the discovery of compounds that are able to modulate re-epithelialization kinetics.

Identification of CVID Patients With Defects in Immune Repertoire Formation or Specification.

Common variable immune deficiency disorder (CVID) is the most clinically relevant cause of antibody failure. It is a highly heterogeneous disease with different underlying etiologies. CVID has been associated with a quantitative B cell defect, however, little is known about the quality of B cells present. Here, we studied the naïve and antigen selected B-cell receptor (BCR) repertoire in 33 CVID patients using next generation sequencing, to investigate B cells quality. Analysis for each individual patient revealed whether they have a defect in immune repertoire formation [V(D)J recombination] or specification (somatic hypermutation, subclass distribution, or selection). The naïve BCR repertoire was normal in most of the patients, although alterations in repertoire diversity and the junctions were found in a limited number of patients indicating possible defects in early B-cell development or V(D)J recombination in these patients. In contrast, major differences were found in the antigen selected BCR repertoire. Here, most patients (15/17) showed a reduced frequency of somatic hypermutation (SHM), changes in subclass distribution and/or minor alterations in antigen selection. Together these data show that in our CVID cohort only a small number of patients have a defect in formation of the naïve BCR repertoire, whereas the clear majority of patients have disturbances in their antigen selected repertoire, suggesting a defect in repertoire specification in the germinal centers of these patients. This highlights that CVID patients not only have a quantitative B cell defect, but that also the quality of, especially post germinal center B cells, is impaired.

Public Clonotypes and Convergent Recombination Characterize the Naïve CD8+ T-Cell Receptor Repertoire of Extremely Preterm Neonates.

Respiratory support improvements have aided survival of premature neonates, but infection susceptibility remains a predominant problem. We previously reported that neonatal mice have a rapidly evolving T-cell receptor (TCR) repertoire that impairs CD8+ T cell immunity. To understand the impact of prematurity on the human CD8+ TCR repertoire, we performed next-generation sequencing of the complementarity-determining region 3 (CDR3) from the rearranged TCR variable beta (Vß) in sorted, naïve CD8+ T cells from extremely preterm neonates (23-27 weeks gestation), term neonates (37-41 weeks gestation), children (16-56 months), and adults (25-50 years old). Strikingly, preterm neonates had an increased frequency of public clonotypes shared between unrelated individuals. Public clonotypes identified in preterm infants were encoded by germline gene sequences, and some of these clonotypes persisted into adulthood. The preterm neonatal naïve CD8+ TCR repertoire exhibited convergent recombination, characterized by different nucleotide sequences encoding the same amino acid CDR3 sequence. As determined by Pielou's evenness and iChao1 metrics, extremely preterm neonates have less clonality, and a much lower bound for the number of unique TCR within an individual preterm neonate, which indicates a less rich and diverse repertoire, as compared to term neonates, children, and adults. This suggests that T cell selection in the preterm neonate may be less stringent or different. Our analysis is the first to compare the TCR repertoire of naïve CD8+ T cells between viable preterm neonates and term neonates. We find preterm neonates have a repertoire immaturity which potentially contributes to their increased infection susceptibility. A developmentally regulated, evenly distributed repertoire in preterm neonates may lead to the inclusion of public TCR CDR3ß sequences that overlap between unrelated individuals in the preterm repertoire.

Evaluation of the Antigen-Experienced B-Cell Receptor Repertoire in Healthy Children and Adults.

Upon antigen recognition via their B cell receptor (BR), B cells migrate to the germinal center where they undergo somatic hypermutation (SHM) to increase their affinity for the antigen, and class switch recombination (CSR) to change the effector function of the secreted antibodies. These steps are essential to create an antigen-experienced BR repertoire that efficiently protects the body against pathogens. At the same time, the BR repertoire should be selected to protect against responses to self-antigen or harmless antigens. Insights into the processes of SHM, selection, and CSR can be obtained by studying the antigen-experienced BR repertoire. Currently, a large reference data set of healthy children and adults, which ranges from neonates to the elderly, is not available. In this study, we analyzed the antigen-experienced repertoire of 38 healthy donors (HD), ranging from cord blood to 74 years old, by sequencing IGA and IGG transcripts using next generation sequencing. This resulted in a large, freely available reference data set containing 412,890 IGA and IGG transcripts. We used this data set to study mutation levels, SHM patterns, antigenic selection, and CSR from birth to elderly HD. Only small differences were observed in SHM patterns, while the mutation levels increase in early childhood and stabilize at 6 years of age at around 7%. Furthermore, comparison of the antigen-experienced repertoire with sequences from the naive immune repertoire showed that features associated with autoimmunity such as long CDR3 length and IGHV4-34 usage are reduced in the antigen-experienced repertoire. Moreover, IGA2 and IGG2 usage was increased in HD in higher age categories, while IGG1 usage was decreased. In addition, we studied clonal relationship in the different samples. Clonally related sequences were found with different subclasses. Interestingly, we found transcripts with the same CDR1-CDR3 sequence, but different subclasses. Together, these data suggest that a single antigen can provoke a B-cell response with BR of different subclasses and that, during the course of an immune response, some B cells change their isotype without acquiring additional SHM or can directly switch to different isotypes.

Strategies for B-cell receptor repertoire analysis in primary immunodeficiencies: from severe combined immunodeficiency to common variable immunodeficiency.

The antigen receptor repertoires of B- and T-cells form the basis of the adaptive immune response. The repertoires should be sufficiently diverse to recognize all possible pathogens. However, careful selection is needed to prevent responses to self or harmless antigens. Limited antigen receptor repertoire diversity leads to immunodeficiency, whereas unselected or misdirected repertoires can result in autoimmunity. The antigen receptor repertoire harbors information about abnormalities in many immunological disorders. Recent developments in next generation sequencing allow the analysis of the antigen receptor repertoire in much greater detail than ever before. Analyzing the antigen receptor repertoire in patients with mutations in genes responsible for the generation of the antigen receptor repertoire will give new insights into repertoire formation and selection. In this perspective, we describe strategies and considerations for analysis of the naive and antigen-selected B-cell repertoires in primary immunodeficiency patients with a focus on severe combined immunodeficiency and common variable immunodeficiency.

Automated Selection of Hotspots (ASH): enhanced automated segmentation and adaptive step finding for Ki67 hotspot detection in adrenal cortical cancer.

BACKGROUND: In prognosis and therapeutics of adrenal cortical carcinoma (ACC), the selection of the most active areas in proliferative rate (hotspots) within a slide and objective quantification of immunohistochemical Ki67 Labelling Index (LI) are of critical importance. In addition to intratumoral heterogeneity in proliferative rate i.e. levels of Ki67 expression within a given ACC, lack of uniformity and reproducibility in the method of quantification of Ki67 LI may confound an accurate assessment of Ki67 LI. RESULTS: We have implemented an open source toolset, Automated Selection of Hotspots (ASH), for automated hotspot detection and quantification of Ki67 LI. ASH utilizes NanoZoomer Digital Pathology Image (NDPI) splitter to convert the specific NDPI format digital slide scanned from the Hamamatsu instrument into a conventional tiff or jpeg format image for automated segmentation and adaptive step finding hotspots detection algorithm. Quantitative hotspot ranking is provided by the functionality from the open source application ImmunoRatio as part of the ASH protocol. The output is a ranked set of hotspots with concomitant quantitative values based on whole slide ranking. CONCLUSION: We have implemented an open source automated detection quantitative ranking of hotspots to support histopathologists in selecting the 'hottest' hotspot areas in adrenocortical carcinoma. To provide wider community easy access to ASH we implemented a Galaxy virtual machine (VM) of ASH which is available from http://bioinformatics.erasmusmc.nl/wiki/Automated_Selection_of_Hotspots . VIRTUAL SLIDES: The virtual slide(s) for this article can be found here: http://www.diagnosticpathology.diagnomx.eu/vs/13000_2014_216.