High-throughput Interpretation of Killer-cell Immunoglobulin-like Receptor Short-read Sequencing Data with PING.

The killer-cell immunoglobulin-like receptor (KIR) complex on chromosome 19 encodes receptors that modulate the activity of natural killer cells, and variation in these genes has been linked to infectious and autoimmune disease, as well as having bearing on pregnancy and transplant outcomes. The medical relevance and high variability of KIR genes makes short-read sequencing an attractive technology for interrogating the region, providing a high-throughput, high-fidelity sequencing method that is cost-effective. However, because this gene complex is characterized by extensive nucleotide polymorphism, structural variation including gene fusions and deletions, and a high level of homology between genes, its interrogation at high resolution has been thwarted by bioinformatic challenges, with most studies limited to examining presence or absence of specific genes. Here, we present the PING (Pushing Immunogenetics to the Next Generation) pipeline, which incorporates empirical data, novel alignment strategies and a custom alignment processing workflow to enable high-throughput KIR sequence analysis from short-read data. PING provides KIR gene copy number classification functionality for all KIR genes through use of a comprehensive alignment reference. The gene copy number determined per individual enables an innovative genotype determination workflow using genotype-matched references. Together, these methods address the challenges imposed by the structural complexity and overall homology of the KIR complex. To determine copy number and genotype determination accuracy, we applied PING to European and African validation cohorts and a synthetic dataset. PING demonstrated exceptional copy number determination performance across all datasets and robust genotype determination performance. Finally, an investigation into discordant genotypes for the synthetic dataset provides insight into misaligned reads, advancing our understanding in interpretation of short-read sequencing data in complex genomic regions. PING promises to support a new era of studies of KIR polymorphism, delivering high-resolution KIR genotypes that are highly accurate, enabling high-quality, high-throughput KIR genotyping for disease and population studies.

Common and well-documented (CWD) alleles of human leukocyte antigen-A, -B, -C, -DRB1, and -DQB1 loci for the Chinese Han population do not quite correlate with the ASHI CWD alleles.

Human leukocyte antigen (HLA), which is extremely polymorphic, plays an important role in stem cell transplantation. The Chinese Han comprise a large population of approximately 1.3 billion with diverse HLA alleles that need to be characterized. Data from 3,296 independent, unrelated Chinese Han individuals (1,457 recipients and 1,839 donors) were provided by the China Marrow Donor Program (CMDP) for donor-recipient confirmatory typing. Sequence-based typing, sequence-specific oligonucleotide probe (SSOP)/High Definition-SSOP, and sequence-specific primer methods were used to obtain 4-digit alleles. A total of 49, 86, 50, 63, and 24 HLA-A, -B, -C, -DRB1, and -DQB1 alleles were observed. Following American Society for Histocompatibility and Immunogenetics (ASHI) common and well-documented (CWD) criteria, CWD alleles for Chinese Han in our laboratory test and other laboratory reports do not quite correlate with the ASHI CWD alleles: A*11:53, A*02:34, A*02:53N, B*27:24, B*46:02, B*55:12, C*01:06, C*03:17, C*06:06, C*07:66, C*07:67, C*08:22, DRB1*12:10, DQB1*03:13, and DQB1*06:05 are CWD, but are not included in the ASHI CWD list. A series of alleles are well-documented alleles and are listed in the ASHI CWD list. Conversely, A*26:03, B*51:03, C*12:05, C*15:09, C*15:11, C*17:03, DRB1*11:07, DRB1*11:11, DRB1*13:05, DRB1*13:13, DRB1*14:06, DRB1*14:12, DRB1*14:22, DRB1*14:25, and DQB1*06:11 are rare alleles, but are included in the ASHI CWD list. HLA ethnic diversity is the main reason for the differences in HLA alleles worldwide. The ASHI HLA CWD alleles help reduce the workload and expenses in high-resolution donor registries and the HLA allele frequencies provide a basis from which to predict the chances of finding HLA matching donors. Our data will be meaningful for the CMDP, for other worldwide donor registries, and for an updated ASHI CWD allele list.

IMGT, the international ImMunoGeneTics information system for Immunoinformatics. Methods for querying IMGT databases, tools, and Web resources in the context of immunoinformatics.

IMGT, the international ImMunoGeneTics information system (http://imgt.cines.fr), was created in 1989 by the Laboratoire d'lmmunoGénétique Moléculaire (LIGM) (Université Montpellier II and CNRS) at Montpellier, France, in order to standardize and manage the complexity of immunogenetics data. IMGT is recognized as the international reference in immunogenetics and immunoinformatics. IMGT is a high quality integrated knowledge resource, specialized in (i) the immunoglobulin (IG), T cell receptors (TR), major histocompatibility complex (MHC) of human and other vertebrates; (ii) proteins that belong to the immunoglobulin superfamily (IgSF) and to the MHC superfamily (MhcSF); and (iii) related proteins of the immune systems (RPI) of any species. IMGT provides a common access to standardized data from genome, proteome, genetics, and three-dimensional (3D) structures for the IG, TR, MHC, IgSF, MhcSF, and RPI. IMGT interactive on-line tools are provided for genome, sequence, and 3D structure analysis. IMGT Web resources comprise 8,000 HTML pages of synthesis and knowledge (IMGT Scientific chart, IMGT Repertoire, IMGT Education, etc.) and external links (IMGT Bloc-notes and IMGT other accesses).

The IMGT/HLA database.

The human leukocyte antigen (HLA) complex is located within the 6p21.3 region on the short arm of human chromosome 6 and contains more than 220 genes of diverse function. Many of the genes encode proteins of the immune system and include many highly polymorphic HLA genes. The naming of new HLA genes and allele sequences and their quality control is the responsibility of the WHO Nomenclature Committee for Factors of the HLA System. The IMGT/HLA Database acts as the repository for these sequences and is recognized as the primary source of up-to-date and accurate HLA sequences. The IMGT/HLA website provides a number of tools for accessing the database: these include allele reports, sequence alignments, and sequence similarity searches. The website is updated every 3 months with all the new and confirmatory sequences submitted to the WHO Nomenclature Committee. Submission of HLA sequences to the committee is possible through the tools provided by the IMGT/HLA Database.

IPD: the Immuno Polymorphism Database.

The Immuno Polymorphism Database (IPD) (http://www.ebi.ac.uk/ipd/) is a set of specialist databases related to the study of polymorphic genes in the immune system. IPD currently consists of four databases: IPD-KIR, contains the allelic sequences of killer cell immunoglobulin-like receptors (KIRs); IPD-MHC, a database of sequences of the major histocompatibility complex (MHC) of different species; IPD-HPA, alloantigens expressed only on platelets; and IPD-ESTAB, which provides access to the European Searchable Tumour Cell Line Database, a cell bank of immunologically characterized melanoma cell lines. The IPD project works with specialist groups or nomenclature committees who provide and curate individual sections before they are submitted to IPD for online publication. The IPD project stores all the data in a set of related databases. Those sections with similar data, such as IPD-KIR and IPD-MHC, share the same database structure.

SYFPEITHI: database for searching and T-cell epitope prediction.

Reverse immunology has been used for about 12 years in order to identify T-cell epitopes from pathogens or tumor-associated antigens. In this chapter, we discuss the advantages and pitfalls of T-cell epitope prediction compared to classical experimental procedures such as epitope mapping and cloning experiments. We introduce our three established programs, SYFPEITHI, PAProc, and SNEP, which are freely accessible at no cost in the World Wide Web for the prediction of either HLA-peptide binding or proteasomal processing of antigens. We demonstrate the performance of our epitope prediction programs with several examples and in comparison to other epitope prediction programs available. We also reflect the actual possibilities and limitations of such computer-aided work.

Searching and mapping of T-cell epitopes, MHC binders, and TAP binders.

This chapter describes searching and mapping tools of MHCBN database, which is a curated database. It comprises over 23,000 peptide sequences, whose binding affinity with major histocompatibility complex (MHC) or transporter associated with antigen processing (TAP) molecules has been assayed experimentally. Each entry of the database provides full information (such as sequence, its MHC- or TAP-binding specificity, and source protein) about peptide whose binding affinity (IC50) and T-cell activity is experimentally determined. MHCBN has number of web-based tools for analyzing and retrieving information. In this chapter, we describe how to use web tools integrated in MHCBN that include (i) mapping of experimentally determined antigenic regions on the query sequence, (ii) creation of allele-specific peptide data set, and (iii) BLAST search against MHC or antigen databases.

Searching and mapping of B-cell epitopes in Bcipep database.

One of the major challenges in the field of subunit vaccine design is to identify the antigenic regions in an antigen, which can activate B cell. These antigenic regions are called B-cell epitopes. In this chapter, we describe how to use Bcipep, which is a database of experimentally determined linear B-cell epitopes of varying immunogenicity collected from literature and other publicly available databases. The current version of Bcipep database contains 3,031 entries that include 763 immunodominant, 1,797 immunogenic, and 471 null-immunogenic epitopes. The database provides a set of tools for analysis and extraction of data that includes keyword search, peptide mapping, and BLAST search. The database is available at http://www.imtech.res.in/raghava/bcipep/.

Searching haptens, carrier proteins, and anti-hapten antibodies.

Haptens are small molecules that are usually nonimmunogenic unless coupled to some carrier proteins. The generation of anti-hapten antibodies is important for the development of immunodiagnostics and therapeutics. Recently, our group has developed a database called HaptenDB, which provides comprehensive information about 1,087 haptens. In this chapter, we describe following web tools integrated in HaptenDB: (i) keyword search facility allows search on major fields, (ii) browsing service, to display all haptens, carrier proteins and antibodies, and (iii) structure similarity search, which allows the users to search their structure against hapten structures.

The classification of HLA supertypes by GRID/CPCA and hierarchical clustering methods.

Biological experiments often produce enormous amount of data, which are usually analyzed by data clustering. Cluster analysis refers to statistical methods that are used to assign data with similar properties into several smaller, more meaningful groups. Two commonly used clustering techniques are introduced in the following section: principal component analysis (PCA) and hierarchical clustering. PCA calculates the variance between variables and groups them into a few uncorrelated groups or principal components (PCs) that are orthogonal to each other. Hierarchical clustering is carried out by separating data into many clusters and merging similar clusters together. Here, we use an example of human leukocyte antigen (HLA) supertype classification to demonstrate the usage of the two methods. Two programs, Generating Optimal Linear Partial Least Square Estimations (GOLPE) and Sybyl, are used for PCA and hierarchical clustering, respectively. However, the reader should bear in mind that the methods have been incorporated into other software as well, such as SIMCA, statistiXL, and R.

Analysis of VH gene sequences using two web-based immunogenetics resources gives different results, but the affinity maturation status of chronic lymphocytic leukaemia clones as assessed from either of the resulting data sets has no prognostic significance.

Some cellular and molecular features of chronic lymphocytic leukaemia (CLL) cells that are associated with prognosis may reflect the context within which their progenitors encountered antigen. It follows that the nature of antigen drive in CLL could influence the clinical course and we were prompted to assess the impact, if any, of affinity maturation (an antigen-driven process) on prognosis. Statistical models for assessing affinity maturation status are typically applied to V(H) gene sequence data analysed using a web-based resource like IMGT or VBASE. Since these resources differ with respect to some key relevant features, we evaluated a cohort of CLL cases by applying statistical models to V(H) data derived from both IMGT and VBASE. Important differences between the resulting data sets became apparent. These resulted from database variance and because IMGT and VBASE define complementarity-determining and framework regions (CDRs, FRs) in different ways. Thus, the numbers of mutations identified and their distribution between CDRs/FRs varied between the data sets for the majority of clones. Consequently, two different but overlapping sets of cases with evidence of affinity maturation were defined. Notwithstanding their differences, no significant associations of affinity maturation status with CD38 expression, p53 functional status or survival were identifiable in either data set.

IMGT-ONTOLOGY and IMGT databases, tools and Web resources for immunogenetics and immunoinformatics.

The international ImMunoGeneTics information system (IMGT; http://imgt.cines.fr), is a high quality integrated information system specialized in immunoglobulins (IG), T cell receptors (TR), major histocompatibility complex (MHC), and related proteins of the immune system (RPI) of human and other vertebrates, created in 1989, by the Laboratoire d'ImmunoGénétique Moléculaire (LIGM; Université Montpellier II and CNRS) at Montpellier, France. IMGT provides a common access to standardized data which include nucleotide and protein sequences, oligonucleotide primers, gene maps, genetic polymorphisms, specificities, 2D and 3D structures. IMGT consists of several sequence databases (IMGT/LIGM-DB, IMGT/MHC-DB, IMGT/PRIMER-DB), one genome database (IMGT/GENE-DB) and one 3D structure database (IMGT/3Dstructure-DB), interactive tools for sequence analysis (IMGT/V-QUEST, IMGT/JunctionAnalysis, IMGT/PhyloGene, IMGT/Allele-Align), for genome analysis (IMGT/GeneSearch, IMGT/GeneView, IMGT/LocusView) and for 3D structure analysis (IMGT/StructuralQuery), and Web resources ("IMGT Marie-Paule page") comprising 8000 HTML pages. IMGT other accesses include SRS, FTP, search by BLAST, etc. By its high quality and its easy data distribution, IMGT has important implications in medical research (repertoire in autoimmune diseases, AIDS, leukemias, lymphomas, myelomas), veterinary research, genome diversity and genome evolution studies of the adaptive immune responses, biotechnology related to antibody engineering (single chain Fragment variable (scFv), phage displays, combinatorial libraries) and therapeutical approaches (grafts, immunotherapy). IMGT is freely available at http://imgt.cines.fr.

IMGT, the international ImMunoGeneTics information system, http://imgt.cines.fr.

IMGT, the international ImMunoGeneTics information system (http://imgt.cines.fr), is a high quality integrated knowledge resource specializing in immunoglobulins (IG), T cell receptors (TR) and major histocompatibility complexes (MHC) and related proteins of the immune system (RPI) of human and other vertebrates, created in 1989 by LIGM at the Université Montpellier II, CNRS, Montpellier, France. IMGT provides a common access to standardized data which include nucleotide and protein sequences, oligonucleotide primers, gene maps, genetic polymorphisms, specificities, and 2D and 3D structures. IMGT includes five databases (IMGT/LIGM-DB, IMGT/3Dstructure-DB, IMGT/MHC-DB, IMGT/PRIMER-DB, IMGT/GENE-DB) Web resources ('IMGT Marie-Paule page') and interactive tools (IMGT/V-QUEST, IMGT/JunctionAnalysis, IMGT/PhyloGene, IMGT/LocusView, IMGT/GeneView, IMGT/GeneSearch, IMGT/StructureQuery). IMGT data are expertly annotated according to the rules of the IMGT Scientific chart based on IMGTONTOLOGY. IMGT tools are particularly useful for the analysis of the IG and TR repertoires in physiological normal and pathological situations. IMGT has important applications in medical research (autoimmune diseases, AIDS, leukaemias, lymphomas, myelomas), biotechnology related to antibody engineering (phage displays, combinatorial libraries) and therapeutic approaches (graft, immunotherapy). IMGT is freely available at http://imgt.cines.fr.

HLA nomenclature and the IMGT/HLA sequence database.

Early in their study it was recognized that the genes encoding the HLA molecules were highly polymorphic and that there was a need for a systematic nomenclature. The result was the WHO Nomenclature Committee for Factors of the HLA System, which first met in 1968, and laid down the criteria for successive meetings. This committee meets regularly to discuss issues of nomenclature and has published 16 major reports documenting firstly the HLA antigens and more recently the genes and alleles. The standardization of HLA antigenic specificities has been controlled by the exchange of typing reagents and cells in the International Histocompatibility Workshops. Since 1989 when a large number of HLA allele sequences were first analysed and named, the job of curating and maintaining a database of sequences has been of prime importance. In 1998 the IMGT/HLA database became the official repository for HLA sequences. In addition to the nucleotide and protein sequences the database contains information of the cell from which the sequence was obtained. The database which provides tools for sequence analysis and the submission of new data, is updated quarterly and now contains over 1500 HLA allele sequences.

Immunogenomics: towards a digital immune system.

One of the major differences that set apart vertebrates from non-vertebrates is the presence of a complex immune system. Over the past 400-500 million years, many novel immune genes and gene families have emerged and their products form sophisticated pathways providing protection against most pathogens. The Human Genome Project has laid the foundation to study these genes and pathways in unprecedented detail. Members of the immunoglobulin (Ig) superfamily alone were found to make up over 2% of human genes possibly constituting the largest gene family in the human genome. A subgroup of these human immune genes, those (among others) involved in antigen processing and presentation, are encoded in a single region, the major histocompatibility complex (MHC) on the short arm of chromosome 6. My laboratory has a long-standing interest in understanding the molecular organization and evolution of the MHC. To this end, we have been generating a range of MHC genomic resources that we make available in the form of maps and databases. Much of the complex data of the immune system can be reduced to binary (on/off) information that can easily be made available and analysed by bioinformatics approaches, thus contributing to better understand immune function via a 'digital immune system'.

Immunogenetic heterogeneity of multiple sclerosis in Sardinia.

Multiple sclerosis (MS) predisposition is thought to be influenced by a complex, yet unclear interaction of genetic and environmental factors. Studying ethnically selected populations may reduce genetic and environmental heterogeneities and help clarify the underlying mechanisms of MS susceptibility. Sardinians kept a homogeneous genetic structure and have among the highest MS frequency rates worldwide. Interestingly, MS in Sardinia is linked to otherwise rare HLA alleles. In this light, recent findings from epidemiological and immunogenetic studies of Sardinian MS are presented. Results confirm that, likely due to significant genetic differences at a microgeographic level, even in this homogeneous population MS is immunogenetically heterogeneous and tends to preferentially cluster in some more archaic areas of the island.

Diabetes mellitus insulino-dependiente en Santiago de Chile: rol patogénico de factores inmunogenéticos y ambientales / Insulin-dependent diabetes mellitus in Santiago, Chile: pathogenic role of immunogenetic and environmental factors

The role of HLA class II alleles in the genetic susceptibility to develop insulin-dependent diabetes mellitus (IDDM) was examined by means of PCR and oligospecific probes in 63 IDDM children and 74 controls subjects. In diabetic patients we found a significant increase in the alleles frequency DR3, DR4, DQB1*0302 and DQA1*0301 compared to the control group, where the most prevalent alleles were DR2, DR14 (DRB1*1402), DQA1*0101 and DQA1*0201. All the risk genotypes in the diabetic group were similar than in other caucasian groups: DR3/DR4-DQB1*0201/0302-DQA1*0301/0501 and DR4/DR4-DQB1*0302/0302-DQA1*0301/0301. The homozygote character no asp57 conferred an absolute risk (AR) of 3.87 and the marker Arg52 an AR of 5.78/100.000 hab year. The homozygosis for both markers (no Asp57+Arg52) had an AR of 7.56/100.000 hab year. Regarding environmental factors associated with IDDM, our population under study showed a low prevalence of infectious agents (mainly mumps and rubella, specifically associated with IDDM) and a high prevalence of effective breast-feeding (over 3 months). These factors could be exercising a protector role in the development of IDDM. The factors that appear to be important in the low incidence of IDDM, the high percentage of breast-feeding children in the population, the reduced frequency of susceptible molecules as DR3, DQB1*0201 (compared to other caucasian groups) and the presence of protective genotypes related to DR13 and DR14 observed in the non diabetic children