Human leukocyte antigen class-I variation is associated with atopic dermatitis: A case-control study.

Atopic dermatitis (AD) is a common immune-medicated skin disease. Previous studies have explored the relationship between Human Leukocyte Antigen (HLA) allelic variation and AD with conflicting results. The aim was to examine HLA Class I genetic variation, specifically peptide binding groove variation, and associations with AD. A case-control study was designed to evaluate HLA class I allelic variation and binding pocket polymorphisms, using next generation sequencing on 464 subjects with AD and 388 without AD. Logistic regression was used to evaluate associations with AD by estimating odds ratios (95% confidence intervals). Significant associations were noted with susceptibility to AD (B*53:01) and protection from AD (A*01:01, A*02:01, B*07:02 and C*07:02). Evaluation of polymorphic residues in Class I binding pockets revealed six amino acid residues conferring protection against AD: A9F (HLA-A, position 9, phenylalanine) [pocket B/C], A97I [pocket C/E], A152V [pocket E], A156R [pocket D/E], B163E [pocket A] and C116S [pocket F]. These findings demonstrate that specific HLA class I components are associated with susceptibility or protection from AD. Individual amino acid residues are relevant to protection from AD and set the foundation for evaluating potential HLA Class I molecules in complex with peptides/antigens that may initiate or interfere with T-cell responses.

Novel and Haplotype Specific MicroRNAs Encoded by the Major Histocompatibility Complex.

The MHC is recognized for its importance in human health and disease. However, many disease-associated variants throughout the region remain of unknown significance, residing predominantly within non-coding regions of the MHC. The characterization of non-coding RNA transcripts throughout the MHC is thus central to understanding the genetic contribution of these variants. Therefore, we characterize novel miRNA transcripts throughout the MHC by performing deep RNA sequencing of two B lymphoblastoid cell lines with completely characterized MHC haplotypes. Our analysis identifies 89 novel miRNA transcripts, 48 of which undergo Dicer-dependent biogenesis and are loaded onto the Argonaute silencing complex. Several of the identified mature miRNA and pre-miRNA transcripts are unique to specific MHC haplotypes and overlap common SNPs. Furthermore, 43 of the 89 identified novel miRNA transcripts lie within linkage disequilibrium blocks that contain a disease-associated SNP. These disease associated SNPs are associated with 65 unique disease phenotypes, suggesting that these transcripts may play a role in the etiology of numerous diseases associated with the MHC. Additional in silico analysis reveals the potential for thousands of putative pre-miRNA encoding loci within the MHC that may be expressed by different cell types and at different developmental stages.

Determining performance characteristics of an NGS-based HLA typing method for clinical applications.

This study presents performance specifications of an in-house developed human leukocyte antigen (HLA) typing assay using next-generation sequencing (NGS) on the Illumina MiSeq platform. A total of 253 samples, previously characterized for HLA-A, -B, -C, -DRB1 and -DQB1 were included in this study, which were typed at high-resolution using a combination of Sanger sequencing, sequence-specific primer (SSP) and sequence-specific oligonucleotide probe (SSOP) technologies and recorded at the two-field level. Samples were selected with alleles that cover a high percentage of HLA specificities in each of five different race/ethnic groups: European, African-American, Asian Pacific Islander, Hispanic and Native American. Sequencing data were analyzed by two software programs, Omixon's target and GenDx's NGSengine. A number of metrics including allele balance, sensitivity, specificity, precision, accuracy and remaining ambiguity were assessed. Data analyzed by the two software systems are shown independently. The majority of alleles were identical in the exonic sequences (third field) with both programs for HLA-A, -B, -C and -DQB1 in 97.7% of allele determinations. Among the remaining discrepant genotype calls at least one of the analysis programs agreed with the reference typing. Upon additional manual analysis 100% of the 2530 alleles were concordant with the reference HLA genotypes; the remaining ambiguities did not exceed 0.8%. The results demonstrate the feasibility and significant benefit of HLA typing by NGS as this technology is highly accurate, eliminates virtually all ambiguities, provides complete sequencing information for the length of the HLA gene and forms the basis for utilizing a single methodology for HLA typing in the immunogenetics labs.

The dichotomy between disease phenotype databases and the implications for understanding complex diseases involving the major histocompatibility complex.

Many genes related to innate and adaptive immunity reside within the major histocompatibility complex (MHC) and have been associated with a multitude of complex, immune-related disorders. Despite years of genetic study, this region has seen few causative determinants discovered for immune-mediated diseases. Reported associations have been curated in various databases including the Genetic Association Database, NCBI database of clinically relevant variants (ClinVar) and the Human Gene Mutation Database and together capture genetic associations and annotated pathogenic loci within the MHC and across the genome for a variety of complex, immune-mediated diseases. A review of these three distinct databases reveals disparate annotations between associated genes and pathogenic loci, alluding to the polygenic, multifactorial nature of immune-mediated diseases and the pleiotropic character of genes within the MHC. The technical limitations and inherent biases imposed by current approaches and technologies in studying the MHC create a strong case for the need to perform targeted deep sequencing of the MHC and other immunologically relevant loci in order to fully elucidate and study the causative elements of complex immune-mediated diseases.

Towards allele-level human leucocyte antigens genotyping - assessing two next-generation sequencing platforms: Ion Torrent Personal Genome Machine and Illumina MiSeq.

Human leucocyte antigens (HLA) typing has been a challenge due to extreme polymorphism of the HLA genes and limitations of the current technologies and protocols used for their characterization. Recently, next-generation sequencing techniques have been shown to be a well-suited technology for the complete characterization of the HLA genes. However, a comprehensive assessment of the different platforms for HLA typing, describing the limitations and advantages of each of them, has not been presented. We have compared the Ion Torrent Personal Genome Machine (PGM) and Illumina MiSeq, currently the two most frequently used platforms for diagnostic applications, for a number of metrics including total output, quality score per position across the reads and error rates after alignment which can all affect the accuracy of HLA genotyping. For this purpose, we have used one homozygous and three heterozygous well-characterized samples, at HLA-A, HLA-B, HLA-C, HLA-DRB1 and HLA-DQB1. The total output of bases produced by the MiSeq was higher, and they have higher quality scores and a lower overall error rate than the PGM. The MiSeq also has a higher fidelity when sequencing through homopolymer regions up to 9 bp in length. The need to set phase between distant polymorphic sites was more readily achieved with MiSeq using paired-end sequencing of fragments that are longer than those obtained with PGM. Additionally, we have assessed the workflows of the different platforms for complexity of sample preparation, sequencer operation and turnaround time. The effects of data quality and quantity can impact the genotyping results; having an adequate amount of good quality data to analyse will be imperative for confident HLA genotyping. The overall turnaround time can be very comparable between the two platforms; however, the complexity of sample preparation is higher with PGM, while the actual sequencing time is longer with MiSeq.

Identification and characterization of a novel HLA-B hybrid allele, B*08:132 with Next Generation Sequencing.

The new allele is a hybrid between B*08:01:01 and B*42:01:01.

Identification of a novel HLA-DQB1*02 variant allele, DQB1*02:01:06.

The new HLA-DQB1 allele most closely resembles DQB1*02:01:01, differing at a single position 141 (exon 2, codon 15.3).

16(th) IHIW: global distribution of extended HLA haplotypes.

This report describes the project to identify the global distribution of extended HLA haplotypes, a component of 16th International HLA and Immunogenetics Workshop (IHIW), and summarizes the initial analyses of data collected. The project aims to investigate extended HLA haplotypes, compare their distribution among different populations, assess their frequency in hematopoietic stem cell unrelated donor registries and initiate an international family studies database and DNA repository to be made publicly available. HLA haplotypes compiled in immunogenetics laboratories during the evaluation of transplant candidates and related potential donors were analysed. Haplotypes were determined using the pedigree analysis tool publicly available from the National Marrow Donor Program (NMDP) website. Nineteen laboratories from 10 countries (11 laboratories from North America, five from Asia, two from Latin America and one from Australia) contributed data on a total of 1719 families comprised of 7474 individuals. We identified 10393 HLA haplotypes, of which 1682 haplotypes included high-resolution typing at HLA-A, B, C, DRB1 and DQB1 loci. We also present haplotypes containing MICA and other HLA loci and haplotypes containing rare alleles seen in these families. The project will be extended through the 17th IHIW, and investigators interested in joining the project may communicate with the first author.

16(th) IHIW : review of HLA typing by NGS.

Human leucocyte antigen (HLA) genes play an important role in the success of organ transplantation and are associated with autoimmune and infectious diseases. Current DNA-based genotyping methods, including Sanger sequence-based typing (SSBT), have identified a high degree of polymorphism. This level of polymorphism makes high-resolution HLA genotyping challenging, resulting in ambiguous typing results due to an inability to resolve phase and/or defining polymorphisms lying outside the region amplified. Next-generation sequencing (NGS) may resolve the issue through the combination of clonal amplification, which provides phase information, and the ability to sequence larger regions of genes, including introns, without the additional effort or cost associated with current methods. The NGS HLA sequencing project of the 16IHIW aimed to discuss the different approaches to (i) template preparation including short- and long-range PCR amplicons, exome capture and whole genome; (ii) sequencing platforms, including GS 454 FLX, Ion Torrent PGM, Illumina MiSeq/HiSeq and Pacific Biosciences SMRT; (iii) data analysis, specifically allele-calling software. The pilot studies presented at the workshop demonstrated that although individual sequencers have very different performance characteristics, all produced sequence data suitable for the resolution of HLA genotyping ambiguities. The developments presented at this workshop clearly highlight the potential benefits of NGS in the HLA laboratory.

Characterization of a new HLA-C*04 allele, C*04:112.

HLA-C*04:112 differs from HLA-C*04:01:01:01 by one nucleotide at position 270 resulting in an amino acid change, Lysine to Asparagine, at codon 66 of exon 2.

Filling the gaps - the generation of full genomic sequences for 15 common and well-documented HLA class I alleles using next-generation sequencing technology.

Many common and well-documented (CWD) HLA alleles have only been partially characterized. The DNA sequence of these incomplete alleles, as published in the IMGT/HLA database, is most often limited to exons that code for the extracellular domains of the mature protein. Here we describe the application of next-generation sequencing technology to obtain full length genomic sequence from a single long-range PCR amplicon for 15 common and well-documented HLA Class I alleles. This technology is well suited to fill in the gaps of the current HLA allele sequence database which is largely incomplete. A more comprehensive catalog of HLA allele sequences would be beneficial in the evaluation of mismatches in transplantation, studies of population genetics, the evolution of HLAs, regulatory mechanisms and HLA expression, and issues related to the genomic organization of the MHC.

Identification of a novel HLA-A*23 variant allele, A*23:50.

HLA-A*23:50 differs from HLA-A*23:01:01 by one nucleotide at position 112 resulting in an amino acid change, Arginine to Tryptophan, at codon 14 of exon 2.

A multi-site study using high-resolution HLA genotyping by next generation sequencing.

The high degree of polymorphism at human leukocyte antigen (HLA) class I and class II loci makes high-resolution HLA typing challenging. Current typing methods, including Sanger sequencing, yield ambiguous typing results because of incomplete genomic coverage and inability to set phase for HLA allele determination. The 454 Life Sciences Genome Sequencer (GS FLX) next generation sequencing system coupled with conexio atf software can provide very high-resolution HLA genotyping. High-throughput genotyping can be achieved by use of primers with multiplex identifier (MID) tags to allow pooling of the amplicons generated from different individuals prior to sequencing. We have conducted a double-blind study in which eight laboratory sites performed amplicon sequencing using GS FLX standard chemistry and genotyped the same 20 samples for HLA-A, -B, -C, DPB1, DQA1, DQB1, DRB1, DRB3, DRB4, and DRB5 (DRB3/4/5) in a single sequencing run. The average sequence read length was 250 base pairs and the average number of sequence reads per amplicon was 672, providing confidence in the allele assignments. Of the 1280 genotypes considered, assignment was possible in 95% of the cases. Failure to assign genotypes was the result of researcher procedural error or the presence of a novel allele rather than a failure of sequencing technology. Concordance with known genotypes, in cases where assignment was possible, ranged from 95.3% to 99.4% for the eight sites, with overall concordance of 97.2%. We conclude that clonal pyrosequencing using the GS FLX platform and CONEXIO ATF software allows reliable identification of HLA genotypes at high resolution.

Identification of a novel HLA-B*35 variant allele, B*35:151.

HLA-B*35:151 differs from HLA-B*35:68:02 by one nucleotide at codon 116 (TCC>TTC) resulting in an amino acid substitution, serine to phenylalanine.

HLA class II DPB1 and DRB1 polymorphisms associated with genetic susceptibility to beryllium toxicity.

OBJECTIVES: Chronic beryllium disease (CBD) is a hypersensitivity granulomatous pulmonary disease caused by exposure to the metal beryllium (Be²âº). Our objective was to extend current knowledge of the genetics of beryllium disease by examining all HLA-DPB1 and HLA-DPR1 gene polymorphisms and the interactions between them. METHODS: DNA-based typing of HLA-DPB1 and HLA-DRB1 loci at the allele level was performed on 65 CBD, 44 beryllium sensitised (BeS) but without CBD and 288 non-affected, beryllium exposed controls. RESULTS: The DPßE69 residue regardless of zygosity, but particularly if present on non-*0201 alleles, was of primary importance for the development of CBD and BeS, while other negatively charged residues DPßDE55, 56 and DPßDE84, 85 incrementally increased, although not independently, the risk. The DPßE69 positive alleles with charge -7 or -9 were associated with both CBD and BeS. The polymorphic residues DPßE69, DPßDE55, 56 and DPßDE84, 85 were responsible for the -9 charge and the first two residues for the -7 charge. CONCLUSIONS: In the absence of DPßE69, DRßE71 is a risk factor for CBD and BeS. DPßE69 and DRßE71 are adjacent to other amino acids that are also negatively charged, suggesting that the positively charged Be²âº modifies the local environment of the epitopes in a way that promotes interactions between peptides and T cells and results in CBD. Finally, the protective effect of the DPB1*0201 positive haplotype may involve particular polymorphisms outside of the DPB1 gene.

Next-generation sequencing: the solution for high-resolution, unambiguous human leukocyte antigen typing.

Human leukocyte antigen (HLA) typing has been a challenge for more than 50 years. Current methods (Sanger sequencing, sequence-specific primers [SSP], sequence-specific oligonucleotide probes [SSOP]) continue to generate ambiguities that are time-consuming and expensive to resolve. However, next-generation sequencing (NGS) overcomes ambiguity through the combination of clonal amplification, which provides on-phase sequence and a high level of parallelism, whereby millions of sequencing reads are produced enabling an expansion of the HLA regions sequenced. We explored HLA typing using NGS through a three-step process. First, HLA-A, -B, -C, -DRB1, and -DQB1 were amplified with long-range PCR. Subsequently, amplicons were sequenced using the 454 GS-FLX platform. Finally, sequencing data were analyzed with Assign-NG software. In a single experiment, four individual samples and two mixtures were sequenced producing >75 Mb of sequence from >300,000 individual sequence reads (average length, 244 b). The reads were aligned and covered 100% of the regions amplified. Allele assignment was 100% concordant with the known HLA alleles of our samples. Our results suggest this method can be a useful tool for complete genomic characterization of new HLA alleles and for completion of sequence for existing, partially sequenced alleles. NGS can provide complete, unambiguous, high-resolution HLA typing; however, further evaluation is needed to explore the feasibility of its routine use.

HLA and environmental interactions in sarcoidosis.

Sarcoidosis is a systemic granulomatosis of unknown etiology despite being described over 100 years ago. While both genetic predisposition and environmental exposures have been proposed as playing a role in this disease, there have not been any systematic investigations of gene-environmental interaction in this disease. In the ACCESS dataset, detailed environmental histories and high resolution HLA class II typing were performed on 476 cases of newly diagnosed sarcoidosis and 476 matched controls from the patients' community. We evaluated gene-environmental interactions in exposures or HLA class II alleles that were present in > 5% of the population and had an odd ratio of > 1.0. Four exposures and four HLA Class II alleles met these criteria and were evaluated. Significant interaction was observed between HLA DRB1*1101 and insecticide exposure at work (p < 0.10) and suggestive interaction was observed between HLA DRB1*1101 and exposure to mold and musty odors and DRB1*1501 and insecticide exposure at work (P < 0.15). In addition, HLA DRB1*1101 and insecticide exposure at work was associated with extrapulmonary sarcoidosis, specifically cardiac sarcoidosis and hypercalcemia (p<0.05) and HLA DRB1*1101 and exposure to molds and musty odors was associated with pulmonary only sarcoidosis (P < 0.05). These studies suggest that sarcoidosis is due to an interaction of genetic predisposition and environmental exposure in at least some cases of sarcoidosis. Future studies in defined phenotypes of sarcoidosis may be necessary to define environmental and genetic associations with sarcoidosis.

Identification of a novel HLA-DPB1 allele (DPB1*1902) by haplotype-specific extraction and nucleotide sequencing.

The DPB1*1902 allele resulted from a single-nucleotide substitution at codon 35 (TTC --> CTC/F --> L) of exon 2 of DPB1*0402 or 0602.

Unrelated donor or partially matched related donor peripheral stem cell transplant with CD34+ selection and CD3+ addback for pediatric patients with leukemias.

Unmodified peripheral stem cell transplants are associated with an increased risk of extensive chronic GVHD. T depletion may reduce this risk, but the risk of graft failure or relapse may increase. To decrease the risks of both extensive chronic GVHD and graft failure, we added back a defined dose of CD3+ cells to CD34+ selected PSCs. Twenty-four patients were evaluable for outcome analysis. Donors were unrelated (23) or related (1). Conditioning was thiotepa, cyclophosphamide, and total body irradiation. Cyclosporine was used post transplant. Following CD34+ selection, a total of 5 x 10(5)/kg CD3+ cells were infused. Donors were matched for 12 patients. The median CD34+ dose infused was 7.1 x 10(6)/kg. Engraftment occurred in all patients at a median of 14 days (10-19). Twelve patients are alive in remission 15-34 months (median, 25) post PSCT. GVHD occurred in 17 patients, but was >grade II in only 2. Chronic GVHD occurred in 61.5% of evaluable patients, but was limited to skin and perioral cavity. Two patients relapsed, and 10 patients died of non-relapse causes. This study demonstrates that PSCT with CD34+ selection and a defined dose of CD3+ results in prompt engraftment and may limit development of extensive chronic GVHD.

Oligoclonal T cells are infiltrating the brains of children with AIDS: sequence analysis reveals high proportions of identical beta-chain T-cell receptor transcripts.

We have recently described the presence of perivascular CD3+ CD45RO+ T cells infiltrating the brains of children with AIDS. To determine whether these infiltrates contain oligoclonal populations of T cells, we amplified by PCR beta-chain T-cell receptor (TCR) transcripts from autopsy brains of four paediatric patients with AIDS. The amplified transcripts were cloned and sequenced. Sequence analysis of the beta-chain TCR transcripts from all four patients revealed multiple identical copies of TCR beta-chain transcripts, suggesting the presence of oligoclonal populations of T-cells. These TCR transcripts were novel. The presence of oligoclonal populations of T cells in the brains of these four paediatric patients with AIDS suggests that these T cells have undergone antigen-driven proliferation and clonal expansion very likely in situ, in the brains of these AIDS patients, in response to viral or self-antigens. Although the specificity of the clonally expanded beta-chain TCR transcripts remains to be elucidated, none of the beta-chain TCR transcripts identified in this study were identical to those specific for HIV-1 antigens that are currently reported in the GENBANK/EMBL databases. Certain common CDR3 motifs were observed in brain-infiltrating T cells within and between certain patients. Large proportions (24 of 61; 39%) of beta-chain TCR clones from one patient (NP95-73) and 2 of 27 (7%) of another patient (NP95-184-O) exhibited substantial CDR3 homology to myelin basic protein (MBP)-specific TCR derived from normal donors or TCR expressed in the brain of patients with multiple sclerosis (MS) or with viral encephalitis. These two patients (NP95-73 and NP95-184-O) also shared HLA class II with the normal donors and the MS patients who expressed these homologous TCR. Pathologic examination at autopsy of the brains revealed the presence of myelin pallor only in patient NP95-73. T-cell clones identified in the brain of patients NP95-73 and NP95-184-O may recognize MBP or another CNS self antigen and this recognition may be restricted by either DRB1*15 or DQB1*0602 specificities.