1
|
Collins AM, Ohlin M, Corcoran M, Heather JM, Ralph D, Law M, Martínez-Barnetche J, Ye J, Richardson E, Gibson WS, Rodriguez OL, Peres A, Yaari G, Watson CT, Lees WD. AIRR-C IG Reference Sets: curated sets of immunoglobulin heavy and light chain germline genes. Front Immunol 2024; 14:1330153. [PMID: 38406579 PMCID: PMC10884231 DOI: 10.3389/fimmu.2023.1330153] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/30/2023] [Accepted: 12/27/2023] [Indexed: 02/27/2024] Open
Abstract
Introduction Analysis of an individual's immunoglobulin (IG) gene repertoire requires the use of high-quality germline gene reference sets. When sets only contain alleles supported by strong evidence, AIRR sequencing (AIRR-seq) data analysis is more accurate and studies of the evolution of IG genes, their allelic variants and the expressed immune repertoire is therefore facilitated. Methods The Adaptive Immune Receptor Repertoire Community (AIRR-C) IG Reference Sets have been developed by including only human IG heavy and light chain alleles that have been confirmed by evidence from multiple high-quality sources. To further improve AIRR-seq analysis, some alleles have been extended to deal with short 3' or 5' truncations that can lead them to be overlooked by alignment utilities. To avoid other challenges for analysis programs, exact paralogs (e.g. IGHV1-69*01 and IGHV1-69D*01) are only represented once in each set, though alternative sequence names are noted in accompanying metadata. Results and discussion The Reference Sets include less than half the previously recognised IG alleles (e.g. just 198 IGHV sequences), and also include a number of novel alleles: 8 IGHV alleles, 2 IGKV alleles and 5 IGLV alleles. Despite their smaller sizes, erroneous calls were eliminated, and excellent coverage was achieved when a set of repertoires comprising over 4 million V(D)J rearrangements from 99 individuals were analyzed using the Sets. The version-tracked AIRR-C IG Reference Sets are freely available at the OGRDB website (https://ogrdb.airr-community.org/germline_sets/Human) and will be regularly updated to include newly observed and previously reported sequences that can be confirmed by new high-quality data.
Collapse
Affiliation(s)
- Andrew M. Collins
- School of Biotechnology and Biomolecular Sciences, University of New South Wales, Sydney, NSW, Australia
| | - Mats Ohlin
- Department of Immunotechnology, and SciLifeLab, Lund University, Lund, Sweden
| | - Martin Corcoran
- Department of Microbiology, Tumor and Cell Biology, Karolinska Institute, Stockholm, Sweden
| | - James M. Heather
- Mass General Cancer Center, Massachusetts General Hospital, Charlestown, MA, United States
- Department of Medicine, Harvard Medical School, Boston, MA, United States
| | - Duncan Ralph
- Fred Hutchinson Cancer Research Center, Seattle, WA, United States
| | - Mansun Law
- Department of Immunology and Microbiology, The Scripps Research Institute, La Jolla, CA, United States
| | - Jesus Martínez-Barnetche
- Centro de Investigación Sobre Enfermedades Infecciosas, Instituto Nacional de Salud Pública, Cuernavaca, Morelos, Mexico
| | - Jian Ye
- National Center for Biotechnology Information, National Library of Medicine, National Institutes of Health, Bethesda, MD, United States
| | - Eve Richardson
- La Jolla Institute for Immunology, San Diego, CA, United States
| | - William S. Gibson
- Department of Biochemistry and Molecular Genetics, School of Medicine, University of Louisville, Louisville, KY, United States
| | - Oscar L. Rodriguez
- Department of Biochemistry and Molecular Genetics, School of Medicine, University of Louisville, Louisville, KY, United States
| | - Ayelet Peres
- Bioengineering Program, Faculty of Engineering, Bar-Ilan University, Ramat Gan, Israel
| | - Gur Yaari
- Bioengineering Program, Faculty of Engineering, Bar-Ilan University, Ramat Gan, Israel
| | - Corey T. Watson
- Department of Biochemistry and Molecular Genetics, School of Medicine, University of Louisville, Louisville, KY, United States
| | - William D. Lees
- Institute of Structural and Molecular Biology, Birkbeck College, London, United Kingdom
- Human-Centered Computing and Information Science, Institute for Systems and Computer Engineering, Technology and Science, Porto, Portugal
| |
Collapse
|
2
|
Chen Q, Guo Q, Wang D, Zhu S, Wu D, Wang Z, Lu Y. Diagnosis and prognosis of pancreatic cancer with immunoglobulin heavy constant delta blood marker. J Cancer Res Clin Oncol 2023; 149:12977-12992. [PMID: 37466798 DOI: 10.1007/s00432-023-05161-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/30/2023] [Accepted: 07/09/2023] [Indexed: 07/20/2023]
Abstract
BACKGROUND Pancreatic cancer (PC) is highly malignant and difficult to detect, while few blood markers are currently available for diagnosing PC. METHODS We obtained differential expression genes (DEGs) from GEO (gene expression omnibus) database and assessed by quantitative real-time polymerase chain reaction (qRT-PCR), receiver operating characteristic (ROC), univariate and multifactorial regression analysis, and survival analysis in our clinic center. Through the TCGA (the cancer genome atlas) database, we analyzed functional enrichment, different risk groups with survival analysis, immunological features, and the risk score established by the Cox regression model and constructed a nomogram. RESULT Immunoglobulin heavy constant delta (IGHD) was remarkably upregulated in peripheral blood from PC patients, and IGHD was a potential independent biomarker for PC diagnosis (ROC sensitivity, 76.0%; specificity, 74.2%; area under the curve (AUC) = 0.817; univariate logistic regression analysis: odds ratio (OR) 1.488; 95% confidence interval (CI) 1.182-1.872; P < 0.001; multiple logistic: OR 2.097; 95% CI 1.276-3.389, P = 0.003). In addition, the IGHD expression was remarkably reduced after resectioning the primary tumor. High IGHD expression indicated higher lymphocyte infiltration and increased activities of immunological pathways in PC patients. KRAS and SMAD were observed with a prominent difference among top mutated genes between the two groups. The risk score predicted reliable clinical prognosis and drug responses. Furthermore, a nomogram with the risk score and clinical characteristics was constructed, showing a better predictive performance. CONCLUSION IGHD is a valuable PC diagnosis, prognosis, and therapeutic response marker.
Collapse
Affiliation(s)
- Qiyang Chen
- Department of Hepatobiliary and Pancreatic Surgery, Affiliated Hospital of Nantong University, Medical School of Nantong University, Nantong, 226001, Jiangsu, China
| | - Qingsong Guo
- Department of Hepatobiliary and Pancreatic Surgery, Affiliated Hospital of Nantong University, Medical School of Nantong University, Nantong, 226001, Jiangsu, China
| | - Dongzhi Wang
- Department of Hepatobiliary and Pancreatic Surgery, Affiliated Hospital of Nantong University, Medical School of Nantong University, Nantong, 226001, Jiangsu, China
| | - Shajun Zhu
- Department of Hepatobiliary and Pancreatic Surgery, Affiliated Hospital of Nantong University, Medical School of Nantong University, Nantong, 226001, Jiangsu, China
| | - Di Wu
- Department of Hepatobiliary and Pancreatic Surgery, Affiliated Hospital of Nantong University, Medical School of Nantong University, Nantong, 226001, Jiangsu, China
| | - Zhiwei Wang
- Department of Hepatobiliary and Pancreatic Surgery, Affiliated Hospital of Nantong University, Medical School of Nantong University, Nantong, 226001, Jiangsu, China
| | - Yuhua Lu
- Department of Hepatobiliary and Pancreatic Surgery, Affiliated Hospital of Nantong University, Medical School of Nantong University, Nantong, 226001, Jiangsu, China.
| |
Collapse
|
3
|
Qiu Q, Zhang P, Zhang N, Shen Y, Lou S, Deng J. Development of a Prognostic Nomogram for Acute Myeloid Leukemia on IGHD Gene Family. Int J Gen Med 2021; 14:4303-4316. [PMID: 34408473 PMCID: PMC8364394 DOI: 10.2147/ijgm.s317528] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/29/2021] [Accepted: 07/15/2021] [Indexed: 11/29/2022] Open
Abstract
Purpose Acute myeloid leukaemia (AML) is a common haematological disease in adults. The overall survival (OS) remains unsatisfactory. It is critical to identify potential prognostic biomarkers and develop a nomogram that predicts overall survival in patients with AML. Patients and Methods We used gene expression dataset and clinical data from The Cancer Genome Atlas (TCGA) and Genotype-Tissue Expression (GTEx) to identify differential expression analysis, survival analysis, and prognostic value of IGHD gene family (IGHDs) in AML patients. A risk score model was built through Lasso analysis and multivariate Cox regression. We also developed a nomogram and evaluated its accuracy with Harrell’s Harmony Index (C-index) and calibration curve. Last, the Therapeutically Applicable Research to Generate Effective Treatments (TARGET) database was used for external validation. Results IGHD1-20 mRNA expression level was an independent prognostic factor for patients with AML by multivariate analysis. After Lasso analysis and multivariate Cox regression, we constructed a 3-gene model (IGHD1-1, IGHD1-20, IGHD3-16) associated with OS in AML. Risk score and age were validated as independent risk factors for prognosis and were used to build a nomogram. The C index and calibration curve results show that its ability to predict 1-year, 3-year and 5-year overall survival is accurate. Conclusion The mRNA level of IGHDs was increased in AML patients. IGHD1-20 was an independent risk factor for OS in AML patients. The IGHDs risk model (IGHD1-1, IGHD1-20, IGHD3-16) relates to the OS of AML patients. The nomogram, including risk score and age, can conveniently and effectively predict the overall survival rate of patients.
Collapse
Affiliation(s)
- Qunxiang Qiu
- Department of Hematology, The Second Affiliated Hospital of Chongqing Medical University, Chongqing, 400010, People's Republic of China
| | - Ping Zhang
- Hematology Laboratory, The Second Affiliated Hospital of Chongqing Medical University, Chongqing, 400010, People's Republic of China
| | - Nan Zhang
- Department of Hematology, The Second Affiliated Hospital of Chongqing Medical University, Chongqing, 400010, People's Republic of China
| | - Yan Shen
- Department of Hematology, The Second Affiliated Hospital of Chongqing Medical University, Chongqing, 400010, People's Republic of China
| | - Shifeng Lou
- Department of Hematology, The Second Affiliated Hospital of Chongqing Medical University, Chongqing, 400010, People's Republic of China
| | - Jianchuan Deng
- Department of Hematology, The Second Affiliated Hospital of Chongqing Medical University, Chongqing, 400010, People's Republic of China
| |
Collapse
|
4
|
Ralph DK, Matsen FA. Per-sample immunoglobulin germline inference from B cell receptor deep sequencing data. PLoS Comput Biol 2019; 15:e1007133. [PMID: 31329576 PMCID: PMC6675132 DOI: 10.1371/journal.pcbi.1007133] [Citation(s) in RCA: 39] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2017] [Revised: 08/01/2019] [Accepted: 05/28/2019] [Indexed: 11/26/2022] Open
Abstract
The collection of immunoglobulin genes in an individual's germline, which gives rise to B cell receptors via recombination, is known to vary significantly across individuals. In humans, for example, each individual has only a fraction of the several hundred known V alleles. Furthermore, the currently-accepted set of known V alleles is both incomplete (particularly for non-European samples), and contains a significant number of spurious alleles. The resulting uncertainty as to which immunoglobulin alleles are present in any given sample results in inaccurate B cell receptor sequence annotations, and in particular inaccurate inferred naive ancestors. In this paper we first show that the currently widespread practice of aligning each sequence to its closest match in the full set of IMGT alleles results in a very large number of spurious alleles that are not in the sample's true set of germline V alleles. We then describe a new method for inferring each individual's germline gene set from deep sequencing data, and show that it improves upon existing methods by making a detailed comparison on a variety of simulated and real data samples. This new method has been integrated into the partis annotation and clonal family inference package, available at https://github.com/psathyrella/partis, and is run by default without affecting overall run time.
Collapse
Affiliation(s)
- Duncan K. Ralph
- Fred Hutchinson Cancer Research Center, Seattle, Washington, United States of America
| | - Frederick A. Matsen
- Fred Hutchinson Cancer Research Center, Seattle, Washington, United States of America
| |
Collapse
|
5
|
Safonova Y, Pevzner PA. De novo Inference of Diversity Genes and Analysis of Non-canonical V(DD)J Recombination in Immunoglobulins. Front Immunol 2019; 10:987. [PMID: 31134072 PMCID: PMC6516046 DOI: 10.3389/fimmu.2019.00987] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2019] [Accepted: 04/16/2019] [Indexed: 12/03/2022] Open
Abstract
The V(D)J recombination forms the immunoglobulin genes by joining the variable (V), diversity (D), and joining (J) germline genes. Since variations in germline genes have been linked to various diseases, personalized immunogenomics aims at finding alleles of germline genes across various patients. Although recent studies described algorithms for de novo inference of V and J genes from immunosequencing data, they stopped short of solving a more difficult problem of reconstructing D genes that form the highly divergent CDR3 regions and provide the most important contribution to the antigen binding. We present the IgScout algorithm for de novo D gene reconstruction and apply it to reveal new alleles of human D genes and previously unknown D genes in camel, an important model organism in immunology. We further analyze non-canonical V(DD)J recombination that results in unusually long CDR3s with tandem fused IGHD genes and thus expands the diversity of the antibody repertoires. We demonstrate that tandem CDR3s represent a consistent and functional feature of all analyzed immunosequencing datasets, reveal ultra-long CDR3s, and shed light on the mechanism responsible for their formation.
Collapse
Affiliation(s)
- Yana Safonova
- Center for Information Theory and Applications, University of California, San Diego, San Diego, CA, United States
| | - Pavel A Pevzner
- Department of Computer Science and Engineering, University of California, San Diego, San Diego, CA, United States
| |
Collapse
|
6
|
On being the right size: antibody repertoire formation in the mouse and human. Immunogenetics 2017; 70:143-158. [DOI: 10.1007/s00251-017-1049-8] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/26/2017] [Accepted: 12/04/2017] [Indexed: 01/01/2023]
|
7
|
Hou D, Chen C, Seely EJ, Chen S, Song Y. High-Throughput Sequencing-Based Immune Repertoire Study during Infectious Disease. Front Immunol 2016; 7:336. [PMID: 27630639 PMCID: PMC5005336 DOI: 10.3389/fimmu.2016.00336] [Citation(s) in RCA: 35] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/06/2016] [Accepted: 08/19/2016] [Indexed: 11/13/2022] Open
Abstract
The selectivity of the adaptive immune response is based on the enormous diversity of T and B cell antigen-specific receptors. The immune repertoire, the collection of T and B cells with functional diversity in the circulatory system at any given time, is dynamic and reflects the essence of immune selectivity. In this article, we review the recent advances in immune repertoire study of infectious diseases, which were achieved by traditional techniques and high-throughput sequencing (HTS) techniques. HTS techniques enable the determination of complementary regions of lymphocyte receptors with unprecedented efficiency and scale. This progress in methodology enhances the understanding of immunologic changes during pathogen challenge and also provides a basis for further development of novel diagnostic markers, immunotherapies, and vaccines.
Collapse
Affiliation(s)
- Dongni Hou
- Department of Pulmonary Medicine, Zhongshan Hospital, Fudan University , Shanghai , China
| | - Cuicui Chen
- Department of Pulmonary Medicine, Zhongshan Hospital, Fudan University , Shanghai , China
| | - Eric John Seely
- Department of Medicine, Division of Pulmonary and Critical Care Medicine, University of California San Francisco , San Francisco, CA , USA
| | - Shujing Chen
- Department of Pulmonary Medicine, Zhongshan Hospital, Fudan University , Shanghai , China
| | - Yuanlin Song
- Department of Pulmonary Medicine, Zhongshan Hospital, Fudan University , Shanghai , China
| |
Collapse
|
8
|
Ralph DK, Matsen FA. Consistency of VDJ Rearrangement and Substitution Parameters Enables Accurate B Cell Receptor Sequence Annotation. PLoS Comput Biol 2016; 12:e1004409. [PMID: 26751373 PMCID: PMC4709141 DOI: 10.1371/journal.pcbi.1004409] [Citation(s) in RCA: 82] [Impact Index Per Article: 10.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/04/2015] [Accepted: 06/20/2015] [Indexed: 11/18/2022] Open
Abstract
VDJ rearrangement and somatic hypermutation work together to produce antibody-coding B cell receptor (BCR) sequences for a remarkable diversity of antigens. It is now possible to sequence these BCRs in high throughput; analysis of these sequences is bringing new insight into how antibodies develop, in particular for broadly-neutralizing antibodies against HIV and influenza. A fundamental step in such sequence analysis is to annotate each base as coming from a specific one of the V, D, or J genes, or from an N-addition (a.k.a. non-templated insertion). Previous work has used simple parametric distributions to model transitions from state to state in a hidden Markov model (HMM) of VDJ recombination, and assumed that mutations occur via the same process across sites. However, codon frame and other effects have been observed to violate these parametric assumptions for such coding sequences, suggesting that a non-parametric approach to modeling the recombination process could be useful. In our paper, we find that indeed large modern data sets suggest a model using parameter-rich per-allele categorical distributions for HMM transition probabilities and per-allele-per-position mutation probabilities, and that using such a model for inference leads to significantly improved results. We present an accurate and efficient BCR sequence annotation software package using a novel HMM "factorization" strategy. This package, called partis (https://github.com/psathyrella/partis/), is built on a new general-purpose HMM compiler that can perform efficient inference given a simple text description of an HMM.
Collapse
Affiliation(s)
- Duncan K. Ralph
- Fred Hutchinson Cancer Research Center, Seattle, Washington, United States of America
| | - Frederick A. Matsen
- Fred Hutchinson Cancer Research Center, Seattle, Washington, United States of America
- * E-mail:
| |
Collapse
|
9
|
Kidd MJ, Jackson KJL, Boyd SD, Collins AM. DJ Pairing during VDJ Recombination Shows Positional Biases That Vary among Individuals with Differing IGHD Locus Immunogenotypes. THE JOURNAL OF IMMUNOLOGY 2015; 196:1158-64. [PMID: 26700767 DOI: 10.4049/jimmunol.1501401] [Citation(s) in RCA: 28] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/19/2015] [Accepted: 10/20/2015] [Indexed: 11/19/2022]
Abstract
Human IgH diversity is influenced by biases in the pairing of IGHD and IGHJ genes, but these biases have not been described in detail. We used high-throughput sequencing of VDJ rearrangements to explore DJ pairing biases in 29 individuals. It was possible to infer three contrasting IGHD-IGHJ haplotypes in nine of these individuals, and two of these haplotypes include deletion polymorphisms involving multiple contiguous IGHD genes. Therefore, we were able to explore how the underlying genetic makeup of the H chain locus influences the formation of particular DJ pairs. Analysis of nonproductive rearrangements demonstrates that 3' IGHD genes tend to pair preferentially with 5' IGHJ genes, whereas 5' IGHD genes pair preferentially with 3' IGHJ genes; the relationship between IGHD gene pairing frequencies and IGHD gene position is a near linear one for each IGHJ gene. However, striking differences are seen in individuals who carry deletion polymorphisms in the D locus. The absence of different blocks of IGHD genes leads to increases in the utilization frequencies of just a handful of genes, and these genes have no clear positional relationships to the deleted genes. This suggests that pairing frequencies may be influenced by additional complex positional relationships that perhaps arise from chromatin structure. In contrast to IGHD gene usage, IGHJ gene usage is unaffected by the IGHD gene-deletion polymorphisms. Such an outcome would be expected if the recombinase complex associates with an IGHJ gene before associating with an IGHD gene partner.
Collapse
Affiliation(s)
- Marie J Kidd
- School of Biotechnology and Biomolecular Sciences, University of New South Wales, Kensington, Sydney, New South Wales 2052, Australia; and
| | - Katherine J L Jackson
- School of Biotechnology and Biomolecular Sciences, University of New South Wales, Kensington, Sydney, New South Wales 2052, Australia; and Department of Pathology, Stanford University, Stanford, CA 94305
| | - Scott D Boyd
- Department of Pathology, Stanford University, Stanford, CA 94305
| | - Andrew M Collins
- School of Biotechnology and Biomolecular Sciences, University of New South Wales, Kensington, Sydney, New South Wales 2052, Australia; and
| |
Collapse
|
10
|
Old questions, new tools: does next-generation sequencing hold the key to unraveling intestinal B-cell responses? Mucosal Immunol 2015; 8:29-37. [PMID: 25389035 DOI: 10.1038/mi.2014.103] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/28/2014] [Accepted: 09/18/2014] [Indexed: 02/04/2023]
Abstract
Analysis of the intestinal B-cell system and properties of immunoglobulin A, the main antibody isotype produced in the gut, has dominated the rise of mucosal immunology as a discipline. Seminal work established concepts describing the induction, transport, and function of mucosal antibodies. Still, open questions remain and we lack a comprehensive view of how the various sites and pathways of immunoglobulin A induction are integrated to respond to gut antigens. Next-generation sequencing (NGS) offers a novel approach to study B-cell responses, which might substantially enhance our tool box to answer key questions in the field and to take the next steps toward therapeutic exploitation of the mucosal B-cell system. In this review we discuss the potential, challenges, and emerging solutions for gut B-cell repertoire analysis by NGS.
Collapse
|
11
|
Mehr R, Sternberg-Simon M, Michaeli M, Pickman Y. Models and methods for analysis of lymphocyte repertoire generation, development, selection and evolution. Immunol Lett 2012; 148:11-22. [PMID: 22902400 DOI: 10.1016/j.imlet.2012.08.002] [Citation(s) in RCA: 27] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/26/2011] [Revised: 08/01/2012] [Accepted: 08/03/2012] [Indexed: 01/10/2023]
Abstract
T and B cell receptor repertoires are diversified by variable region gene rearrangement and selected based on functionality and lack of self-reactivity. Repertoires can also be defined based on phenotype and function rather than receptor specificity - such as the diversity of T helper cell subsets. Natural killer (NK) cell repertoires, in which each cell expresses a randomly chosen subset of its inhibitory receptor genes, and is educated based on self-MHC recognition by yet unknown mechanisms, are also phenotypic repertoires. Studying the generation, development and selection of lymphocyte repertoires, and their functions during immune responses, is essential for understanding the function of the immune system in healthy individuals and in immune deficient, autoimmune or cancer patients. The study of lymphocyte repertoires will enable clinical immunologists to develop better therapeutic monoclonal antibodies, vaccines, transplantation donor-recipient matching protocols, and other immune intervention strategies. The recent development of high-throughput methods for repertoire data collection - from multicolor flow cytometry through single-cell imaging to deep sequencing - presents us now, for the first time, with the ability to analyze and compare large samples of lymphocyte repertoires in health, aging and disease. The exponential growth of these datasets, however, challenges the theoretical immunology community to develop methods for data organization and analysis. Furthermore, the need to test hypotheses regarding immune function, and generate predictions regarding the outcomes of medical interventions, necessitates the development of complex mathematical and computational models, covering processes on multiple scales, from the genetic and molecular to the cellular and system scales.
Collapse
Affiliation(s)
- Ramit Mehr
- The Mina and Everard Goodman Faculty of Life Sciences, Bar-Ilan University, Ramat-Gan 52900, Israel.
| | | | | | | |
Collapse
|
12
|
Kidd MJ, Chen Z, Wang Y, Jackson KJ, Zhang L, Boyd SD, Fire AZ, Tanaka MM, Gaëta BA, Collins AM. The inference of phased haplotypes for the immunoglobulin H chain V region gene loci by analysis of VDJ gene rearrangements. THE JOURNAL OF IMMUNOLOGY 2011; 188:1333-40. [PMID: 22205028 DOI: 10.4049/jimmunol.1102097] [Citation(s) in RCA: 80] [Impact Index Per Article: 6.2] [Reference Citation Analysis] [Abstract] [Subscribe] [Scholar Register] [Indexed: 01/30/2023]
Abstract
The existence of many highly similar genes in the lymphocyte receptor gene loci makes them difficult to investigate, and the determination of phased "haplotypes" has been particularly problematic. However, V(D)J gene rearrangements provide an opportunity to infer the association of Ig genes along the chromosomes. The chromosomal distribution of H chain genes in an Ig genotype can be inferred through analysis of VDJ rearrangements in individuals who are heterozygous at points within the IGH locus. We analyzed VDJ rearrangements from 44 individuals for whom sufficient unique rearrangements were available to allow comprehensive genotyping. Nine individuals were identified who were heterozygous at the IGHJ6 locus and for whom sufficient suitable VDJ rearrangements were available to allow comprehensive haplotyping. Each of the 18 resulting IGHV│IGHD│IGHJ haplotypes was unique. Apparent deletion polymorphisms were seen that involved as many as four contiguous, functional IGHV genes. Two deletion polymorphisms involving multiple contiguous IGHD genes were also inferred. Three previously unidentified gene duplications were detected, where two sequences recognized as allelic variants of a single gene were both inferred to be on a single chromosome. Phased genomic data brings clarity to the study of the contribution of each gene to the available repertoire of rearranged VDJ genes. Analysis of rearrangement frequencies suggests that particular genes may have substantially different yet predictable propensities for rearrangement within different haplotypes. Together with data highlighting the extent of haplotypic variation within the population, this suggests that there may be substantial variability in the available Ab repertoires of different individuals.
Collapse
Affiliation(s)
- Marie J Kidd
- School of Biotechnology and Biomolecular Sciences, University of New South Wales, Kensington, Sydney, New South Wales 2052, Australia
| | | | | | | | | | | | | | | | | | | |
Collapse
|
13
|
Genomic screening by 454 pyrosequencing identifies a new human IGHV gene and sixteen other new IGHV allelic variants. Immunogenetics 2011; 63:259-65. [PMID: 21249354 DOI: 10.1007/s00251-010-0510-8] [Citation(s) in RCA: 41] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/04/2010] [Accepted: 12/31/2010] [Indexed: 02/02/2023]
Abstract
Complete and accurate knowledge of the genes and allelic variants of the human immunoglobulin gene loci is critical for studies of B cell repertoire development and somatic point mutation, but evidence from studies of VDJ rearrangements suggests that our knowledge of the available immunoglobulin gene repertoire is far from complete. The reported repertoire has changed little over the last 15 years. This is, in part, a consequence of the inefficiencies involved in searching for new members of large, multigenic gene families by cloning and sequencing. The advent of high-throughput sequencing provides a new avenue by which the germline repertoire can be explored. In this report, we describe pyrosequencing studies of the heavy chain IGHV1, IGHV3 and IGHV4 gene subgroups in ten Papua New Guineans. Thousands of 454 reads aligned with complete identity to 51 previously reported functional IGHV genes and allelic variants. A new gene, IGHV3-NL1*01, was identified, which differs from the nearest previously reported gene by 15 nucleotides. Sixteen new IGHV alleles were also identified, 15 of which varied from previously reported functional IGHV genes by between one and four nucleotides, while one sequence appears to be a functional variant of the pseudogene IGHV3-25. BLAST searches suggest that at least six of these new genes are carried within the relatively well-studied populations of North America, Europe or Asia. This study substantially expands the known immunoglobulin gene repertoire and demonstrates that genetic variation of immunoglobulin genes can now be efficiently explored in different human populations using high-throughput pyrosequencing.
Collapse
|
14
|
Jackson KJL, Boyd S, Gaëta BA, Collins AM. Benchmarking the performance of human antibody gene alignment utilities using a 454 sequence dataset. ACTA ACUST UNITED AC 2010; 26:3129-30. [PMID: 21036814 DOI: 10.1093/bioinformatics/btq604] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Abstract
MOTIVATION Immunoglobulin heavy chain genes are formed by recombination of genes randomly selected from sets of IGHV, IGHD and IGHJ genes. Utilities have been developed to identify genes that contribute to observed VDJ rearrangements, but in the absence of datasets of known rearrangements, the evaluation of these utilities is problematic. We have analyzed thousands of VDJ rearrangements from an individual (S22) whose IGHV, IGHD and IGHJ genotype can be inferred from the dataset. Knowledge of this genotype means that the Stanford_S22 dataset can serve to benchmark the performance of IGH alignment utilities. RESULTS We evaluated the performance of seven utilities. Failure to partition a sequence into genes present in the S22 genome was considered an error, and error rates for different utilities ranged from 7.1% to 13.7%. AVAILABILITY Supplementary data includes the S22 genotypes and alignments. The Stanford_S22 dataset and an evaluation tool is available at http://www.emi.unsw.edu.au/~ihmmune/IGHUtilityEval/.
Collapse
Affiliation(s)
- Katherine J L Jackson
- School of Computer Science and Engineering, University of New South Wales, Sydney, NSW, Australia.
| | | | | | | |
Collapse
|
15
|
Boyd SD, Gaëta BA, Jackson KJ, Fire AZ, Marshall EL, Merker JD, Maniar JM, Zhang LN, Sahaf B, Jones CD, Simen BB, Hanczaruk B, Nguyen KD, Nadeau KC, Egholm M, Miklos DB, Zehnder JL, Collins AM. Individual variation in the germline Ig gene repertoire inferred from variable region gene rearrangements. THE JOURNAL OF IMMUNOLOGY 2010; 184:6986-92. [PMID: 20495067 DOI: 10.4049/jimmunol.1000445] [Citation(s) in RCA: 195] [Impact Index Per Article: 13.9] [Reference Citation Analysis] [Abstract] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
Abstract
Individual variation in the Ig germline gene repertoire leads to individual differences in the combinatorial diversity of the Ab repertoire, but the study of such variation has been problematic. The application of high-throughput DNA sequencing to the study of rearranged Ig genes now makes this possible. The sequencing of thousands of VDJ rearrangements from an individual, either from genomic DNA or expressed mRNA, should allow their germline IGHV, IGHD, and IGHJ repertoires to be inferred. In addition, where previously mere glimpses of diversity could be gained from sequencing studies, new large data sets should allow the rearrangement frequency of different genes and alleles to be seen with clarity. We analyzed the DNA of 108,210 human IgH chain rearrangements from 12 individuals and determined their individual IGH genotypes. The number of reportedly functional IGHV genes and allelic variants ranged from 45 to 60, principally because of variable levels of gene heterozygosity, and included 14 previously unreported IGHV polymorphisms. New polymorphisms of the IGHD3-16 and IGHJ6 genes were also seen. At heterozygous loci, remarkably different rearrangement frequencies were seen for the various IGHV alleles, and these frequencies were consistent between individuals. The specific alleles that make up an individual's Ig genotype may therefore be critical in shaping the combinatorial repertoire. The extent of genotypic variation between individuals is highlighted by an individual with aplastic anemia who appears to lack six contiguous IGHD genes on both chromosomes. These deletions significantly alter the potential expressed IGH repertoire, and possibly immune function, in this individual.
Collapse
Affiliation(s)
- Scott D Boyd
- Department of Pathology, Stanford University, Stanford, CA 94305, USA
| | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | |
Collapse
|
16
|
Collins AM, Wang Y, Singh V, Yu P, Jackson KJ, Sewell WA. The reported germline repertoire of human immunoglobulin kappa chain genes is relatively complete and accurate. Immunogenetics 2008; 60:669-76. [PMID: 18712520 DOI: 10.1007/s00251-008-0325-z] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/19/2008] [Accepted: 07/16/2008] [Indexed: 01/31/2023]
Abstract
We describe a bioinformatic analysis of germline and rearranged immunoglobulin kappa chain (IGK) gene sequences, performed in order to assess the completeness and reliability of the reported IGK repertoire. In contrast to the reported heavy-chain gene repertoire, which includes many dubious sequences, only five IGK variable gene (IGKV) alleles appear to have been reported in error. There was, however, insufficient evidence to justify removing these IGKV genes from the germline repertoire. Bioinformatic analysis of apparent mismatches between reported germline genes and 1,863 expressed IGK sequences suggested the existence of two unreported IGKV polymorphisms. Genomic screening of 12 individuals led to the confirmation of both of these polymorphisms, IGKV1-16*02 and IGKV2-30*02. We also show that in contrast to the heavy chain, the IGK repertoire is dominated by sequences that use just a handful of kappa variable (IGKV) and junction (IGKJ) gene pairs. There is also little modification of IGKV and IGKJ genes by the processes of exonuclease removal and N nucleotide addition. The expressed IGK repertoire therefore lacks diversity and the junction region is particularly constrained. Remarkably, the analysis of a dataset of 435 relatively unmutated rearranged kappa genes showed that ten amino acid sequences account for almost 10% of the rearrangements, with identical sequences being derived from as many as seven independent sources. Such dominant sequences are likely to have important roles in the operation of the humoral immune response.
Collapse
Affiliation(s)
- Andrew M Collins
- School of Biotechnology and Biomolecular Sciences, University of New South Wales, Kensington, Australia.
| | | | | | | | | | | |
Collapse
|
17
|
Many human immunoglobulin heavy-chain IGHV gene polymorphisms have been reported in error. Immunol Cell Biol 2007; 86:111-5. [PMID: 18040280 DOI: 10.1038/sj.icb.7100144] [Citation(s) in RCA: 51] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Abstract
The identification of the genes that make up rearranged immunoglobulin genes is critical to many studies. For example, the enumeration of mutations in immunoglobulin genes is important for the prognosis of chronic lymphocytic leukemia, and this requires the accurate identification of the germline genes from which a particular sequence is derived. The immunoglobulin heavy-chain variable (IGHV) gene repertoire is generally considered to be highly polymorphic. In this report, we describe a bioinformatic analysis of germline and rearranged immunoglobulin gene sequences which casts doubt on the existence of a substantial proportion of reported germline polymorphisms. We report a five-level classification system for IGHV genes, which indicates the likelihood that the genes have been reported accurately. The classification scheme also reflects the likelihood that germline genes could be incorrectly identified in mutated VDJ rearrangements, because of similarities to other alleles. Of the 226 IGHV alleles that have previously been reported, our analysis suggests that 104 of these alleles almost certainly include sequence errors, and should be removed from the available repertoire. The analysis also highlights the presence of common mismatches, with respect to the germline, in many rearranged heavy-chain sequences, suggesting the existence of twelve previously unreported alleles. Sequencing of IGHV genes from six individuals in this study confirmed the existence of three of these alleles, which we designate IGHV3-49*04, IGHV3-49*05 and IGHV4-39*07. We therefore present a revised repertoire of expressed IGHV genes, which should substantially improve the accuracy of immunoglobulin gene analysis.
Collapse
|
18
|
Jackson KJL, Gaëta BA, Collins AM. Identifying highly mutated IGHD genes in the junctions of rearranged human immunoglobulin heavy chain genes. J Immunol Methods 2007; 324:26-37. [PMID: 17553518 DOI: 10.1016/j.jim.2007.04.011] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/06/2006] [Revised: 04/19/2007] [Accepted: 04/19/2007] [Indexed: 10/23/2022]
Abstract
The reliable identification of IGHD genes within human immunoglobulin heavy chains is challenging with up to one third of rearrangements having no identifiable IGHD gene. The short, mutated IGHD genes are generally assumed to be indistinguishable from the N-REGIONS of non-template encoded nucleotides that surround them. In this study we have characterised N-REGIONS, demonstrating the importance of nucleotide composition biases in the addition process, including the formation of homopolymer tracts. We then use a simulation approach to determine the likelihood of misidentification of highly mutated IGHD genes among the JUNCTION nucleotides. These likelihoods provide general rules for the identification of mutated D-REGIONs, and suggest that longer D-REGIONs (>25 nucleotides) with as many as ten mutations can be identified with a low risk of error. Shorter D-REGIONs (>16 nucleotides) with as many as four mutations are also identifiable. The reliability of different alignments is dependent upon the junction length (combined N-REGIONs and D-REGION). Data is presented that can guide the alignment of sequences with junction lengths from 5 to 50 nucleotides, including explicit selection between two D-REGION possibilities. The use of such a statistically-based approach to the alignment of IGHD genes will improve the reliability of the partitioning of immunoglobulin sequences, and this in turn will facilitate the study of the many processes that contribute to the diversity of the immunoglobulin repertoire.
Collapse
Affiliation(s)
- Katherine J L Jackson
- School of Biotechnology and Biomolecular Sciences, University of New South Wales, Sydney, Australia
| | | | | |
Collapse
|
19
|
Gaëta BA, Malming HR, Jackson KJL, Bain ME, Wilson P, Collins AM. iHMMune-align: hidden Markov model-based alignment and identification of germline genes in rearranged immunoglobulin gene sequences. Bioinformatics 2007; 23:1580-7. [PMID: 17463026 DOI: 10.1093/bioinformatics/btm147] [Citation(s) in RCA: 105] [Impact Index Per Article: 6.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022] Open
Abstract
MOTIVATION Immunoglobulin heavy chain (IGH) genes in mature B lymphocytes are the result of recombination of IGHV, IGHD and IGHJ germline genes, followed by somatic mutation. The correct identification of the germline genes that make up a variable VH domain is essential to our understanding of the process of antibody diversity generation as well as to clinical investigations of some leukaemias and lymphomas. RESULTS We have developed iHMMune-align, an alignment program that uses a hidden Markov model (HMM) to model the processes involved in human IGH gene rearrangement and maturation. The performance of iHMMune-align was compared to that of other immunoglobulin gene alignment utilities using both clonally related and randomly selected IGH sequences. This evaluation suggests that iHMMune-align provides a more accurate identification of component germline genes than other currently available IGH gene characterization programs. AVAILABILITY iHMMune-align cross-platform Java executable and web interface are freely available to academic users and can be accessed at http://www.emi.unsw.edu.au/~ihmmune/.
Collapse
Affiliation(s)
- Bruno A Gaëta
- School of Biotechnology and Biomolecular Sciences, The University of New South Wales, Sydney, NSW 2052, Australia.
| | | | | | | | | | | |
Collapse
|
20
|
Lee CEH, Jackson KJL, Sewell WA, Collins AM. Use of IGHJ and IGHD gene mutations in analysis of immunoglobulin sequences for the prognosis of chronic lymphocytic leukemia. Leuk Res 2006; 31:1247-52. [PMID: 17169423 DOI: 10.1016/j.leukres.2006.10.013] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/21/2006] [Revised: 09/21/2006] [Accepted: 10/19/2006] [Indexed: 11/28/2022]
Abstract
The level of somatic point mutation in immunoglobulin genes is an important prognostic indicator for patients with chronic lymphocytic leukemia (CLL). Mutation analysis presently focuses solely upon the heavy chain IGHV gene, however mutation is a stochastic process that also targets IGHD and IGHJ genes. Here, we evaluate the completeness and reliability of the reported IGHJ gene repertoire, and demonstrate the likely consequences of the inclusion of IGHD and IGHJ mutations in CLL analysis, using a dataset of 607 sequences. Inclusion of these mutations would lead to the re-classification of many sequences, which should significantly improve the prognostic value of mutation analysis.
Collapse
Affiliation(s)
- Cindy E H Lee
- School of Biotechnology and Biomolecular Sciences, University of New South Wales, Kensington, Sydney, NSW 2052, Australia
| | | | | | | |
Collapse
|