A modular concept of HLA for comprehensive peptide binding prediction

A variety of algorithms have been successful in predicting human leukocyte antigen (HLA)-peptide binding for HLA variants for which plentiful experimental binding data exist. Although predicting binding for only the most common HLA variants may provide sufficient population coverage for vaccine design, successful prediction for as many HLA variants as possible is necessary to understand the immune response in transplantation and immunotherapy. However, the high cost of obtaining peptide binding data limits the acquisition of binding data. Therefore, a prediction algorithm, which applies the binding information from well-studied HLA variants to HLA variants, for which no peptide data exist, is necessary. To this end, a modular concept of class I HLA-peptide binding prediction was developed. Accurate predictions were made for several alleles without using experimental peptide binding data specific to those alleles. We include a comparison of module-based prediction and supertype-based prediction. The modular concept increased the number of predictable alleles from 15 to 75 of HLA-A and 12 to 36 of HLA-B proteins. Under the modular concept, binding data of certain HLA alleles can make prediction possible for numerous additional alleles. We report here a ranking of HLA alleles, which have been identified to be the most informative. Modular peptide binding prediction is freely available to researchers on the web at http://www.peptidecheck.org .     

HLA Typing for the Next Generation

Allele-level resolution data at primary HLA typing is the ideal for most histocompatibility testing laboratories. Many high-throughput molecular HLA typing approaches are unable to determine the phase of observed DNA sequence polymorphisms, leading to ambiguous results. The use of higher resolution methods is often restricted due to cost and time limitations. Here we report on the feasibility of using Pacific Biosciences’ Single Molecule Real-Time (SMRT) DNA sequencing technology for high-resolution and high-throughput HLA typing. Seven DNA samples were typed for HLA-A, -B and -C. The results showed that SMRT DNA sequencing technology was able to generate sequences that spanned entire HLA Class I genes that allowed for accurate allele calling. Eight novel genomic HLA class I sequences were identified, four were novel alleles, three were confirmed as genomic sequence extensions and one corrected an existing genomic reference sequence. This method has the potential to revolutionize the field of HLA typing. The clinical impact of achieving this level of resolution HLA typing data is likely to considerable, particularly in applications such as organ and blood stem cell transplantation where matching donors and recipients for their HLA is of utmost importance.     

Evidence for recessive and against dominant inheritance at the HLA-"linked" locus in coeliac disease.

It has been proposed that gluten sensitive enteropathy (GSE) results from the interaction of two loci: one locus linked to HLA and associated with dominant inheritance, and the other, a non-HLA-linked GSE-associated B-cell alloantigen, exhibiting recessive inheritance. We have shown in previous analyses that a two-locus, dominant-recessive model is less compatible with the existing population prevalence and observed familial segregation data than is a recessive-recessive two-locus model. Here we present additional analyses of reported population and familial HLA data that support the recessive mode of inheritance for the HLA-linked disease locus. Reported data from HLA typing of affected sib pairs, the association of GSE with DR3 and DR7 in different populations, and the proportions of different HLA phenotypes and genotypes were compared with expected data derived by three different methods. The HLA data analyses consistently reject a dominant mode of inheritance for the presumed HLA-linked disease allele but do not reject a recessive model. The affected sib-pair data also support a recessive model. These analyses are consistent with our previous prediction that the HLA-"linked" disease allele in GSE is recessive inherited.     

Signatures of demographic history and natural selection in the human major histocompatibility complex Loci

Many lines of evidence show that several HLA loci have experienced balancing selection. However, distinguishing among demographic and selective explanations for patterns of variation observed with HLA genes remains a challenge. In this study we address this issue using data from a diverse set of human populations at six classical HLA loci and, employing a comparative genomics approach, contrast results for HLA loci to those for non-HLA markers. Using a variety of analytic methods, we confirm and extend evidence for selection acting on several HLA loci. We find that allele frequency distributions for four of the six HLA loci deviate from neutral expectations and show that this is unlikely to be explained solely by demographic factors. Other features of HLA variation are explained in part by demographic history, including decreased heterozygosity and increased LD for populations at greater distances from Africa and a similar apportionment of genetic variation for HLA loci compared to putatively neutral non-HLA loci. On the basis of contrasts among different HLA loci and between HLA and non-HLA loci, we conclude that HLA loci bear detectable signatures of both natural selection and demographic history.     

Pathogen-driven selection and worldwide HLA class I diversity

The human leukocyte antigen (HLA; known as MHC in other vertebrates) plays a central role in the recognition and presentation of antigens to the immune system and represents the most polymorphic gene cluster in the human genome [1]. Pathogen-driven balancing selection (PDBS) has been previously hypothesized to explain the remarkable polymorphism in the HLA complex, but there is, as yet, no direct support for this hypothesis [2 and 3]. A straightforward prediction coming out of the PDBS hypothesis is that populations from areas with high pathogen diversity should have increased HLA diversity in relation to their average genomic diversity. We tested this prediction by using HLA class I genetic diversity from 61 human populations. Our results show that human colonization history explains a substantial proportion of HLA genetic diversity worldwide. However, between-population variation at the HLA class I genes is also positively correlated with local pathogen richness (notably for the HLA B gene), thus providing support for the PDBS hypothesis. The proportion of variations explained by pathogen richness is higher for the HLA B gene than for the HLA A and HLA C genes. This is in good agreement with both previous immunological and genetic data suggesting that HLA B could be under a higher selective pressure from pathogens.     

Efficient incorporation of HLA class II onto human immunodeficiency virus type 1 requires envelope glycoprotein packaging

HLA class II DR is one of the most abundant cell surface proteins incorporated onto human immunodeficiency virus type 1 (HIV-1) during budding. The mechanism for HLA class II protein incorporation is not known and may involve a viral protein. To determine whether Env affects HLA class II protein incorporation, HIV-1 virions, either with or without Env on their surface, were produced from HLA class II-expressing cells and analyzed by whole-virus immunoprecipitation with antisera against HLA class II proteins. HLA class II proteins were detected on virions only when wild-type Env was incorporated, while similar experiments showed that HLA class I proteins were incorporated independent of Env packaging. Therefore, the packaging of HIV-1 Env protein is required for the efficient incorporation of HLA class II but not class I proteins into the virion. Analysis of two Env mutants revealed that the presence of a 43-amino-acid sequence between amino acids 708 and 750 in the gp41(TM) cytoplasmic tail was required for efficient incorporation of HLA class II proteins. These data show that HIV-1 actively incorporates HLA class II proteins in a process that, either directly or indirectly, requires Env.     

Statistical resolution of ambiguous HLA typing data

High-resolution HLA typing plays a central role in many areas of immunology, such as in identifying immunogenetic risk factors for disease, in studying how the genomes of pathogens evolve in response to immune selection pressures, and also in vaccine design, where identification of HLA-restricted epitopes may be used to guide the selection of vaccine immunogens. Perhaps one of the most immediate applications is in direct medical decisions concerning the matching of stem cell transplant donors to unrelated recipients. However, high-resolution HLA typing is frequently unavailable due to its high cost or the inability to re-type historical data. In this paper, we introduce and evaluate a method for statistical, in silico refinement of ambiguous and/or low-resolution HLA data. Our method, which requires an independent, high-resolution training data set drawn from the same population as the data to be refined, uses linkage disequilibrium in HLA haplotypes as well as four-digit allele frequency data to probabilistically refine HLA typings. Central to our approach is the use of haplotype inference. We introduce new methodology to this area, improving upon the Expectation-Maximization (EM)-based approaches currently used within the HLA community. Our improvements are achieved by using a parsimonious parameterization for haplotype distributions and by smoothing the maximum likelihood (ML) solution. These improvements make it possible to scale the refinement to a larger number of alleles and loci in a more computationally efficient and stable manner. We also show how to augment our method in order to incorporate ethnicity information (as HLA allele distributions vary widely according to race/ethnicity as well as geographic area), and demonstrate the potential utility of this experimentally. A tool based on our approach is freely available for research purposes at http://microsoft.com/science.     

MiDAS—Meaningful Immunogenetic Data at Scale

Human immunogenetic variation in the form of HLA and KIR types has been shown to be strongly associated with a multitude of immune-related phenotypes. However, association studies involving immunogenetic loci most commonly involve simple analyses of classical HLA allelic diversity, resulting in limitations regarding the interpretability and reproducibility of results. We here present MiDAS, a comprehensive R package for immunogenetic data transformation and statistical analysis. MiDAS recodes input data in the form of HLA alleles and KIR types into biologically meaningful variables, allowing HLA amino acid fine mapping, analyses of HLA evolutionary divergence as well as experimentally validated HLA-KIR interactions. Further, MiDAS enables comprehensive statistical association analysis workflows with phenotypes of diverse measurement scales. MiDAS thus closes the gap between the inference of immunogenetic variation and its efficient utilization to make relevant discoveries related to immune and disease biology. It is freely available under a MIT license.     

Sequence analysis of three novel DRw14-DRB1 alleles

The nucleotide sequence data reported in this paper have been submitted to the GenBank nucleotide sequence database and have been assigned the accession number M74030 for HLA DRB1*14.7, M74031 for HLA DRB1*14.8, and M74032 for HLA DRB1*14.6.     

以HLAⅡ类转基因鼠研究类风湿关节炎的发病机制

类风湿关节炎（RA）是一种常见的慢性自身免疫性疾病,至今发病原因不明。许多研究表明遗传因素是导致RA发病的最主要原因,而在遗传因素中约35%来自于人类白细胞抗原（HLA）Ⅱ类复合体。因此,对于HLA Ⅱ类复合体参与RA发病的分子机制研究一直是人们关注的热点,而表达HLA Ⅱ类分子的转基因鼠是研究HLA Ⅱ类复合体与RA发病关系最佳的平台。目前,国外已建立了几个HLA Ⅱ类转基因鼠品系,为RA发病机制的研究奠定了很好的基础。本文对HLA Ⅱ类转基因鼠及以此为基础的RA相关研究进行综述。     

Evaluation of photochemical properties of compost humic-like materials

Humic-like acids (HLA) were extracted from compost at the beginning and after 70, 130 and 730 days of maturation in order to be investigated for their ability to induce the transformation of 2,4,6-trimethylphenol under irradiation at 365 nm. The rate of 2,4,6-trimethylphenol phototransformation in the presence of HLA (25 mg l(-1)) varied within HLA 0相似文献   

On estimating HLA/disease association with application to a study of aplastic anemia

A multiplicative model is described relating HLA typing information to disease incidence. A likelihood-based method for estimating parameters in this model is proposed for use with data sets in which HLA haplotype information is available on a series of cases and their parents. This approach is extended to incorporate information from a matched control series for the purpose of estimating HLA and environmental risk factor effects simultaneously. The method is applied to data from aplastic anemia patients treated by bone marrow transplantation and the results are compared to unmatched case-control analyses using the same case series and several different control series.     

Distribution of HLA class I altered phenotypes in colorectal carcinomas: high frequency of HLA haplotype loss associated with loss of heterozygosity in chromosome region 6p21

HLA class I loss or down-regulation is a widespread mechanism used by tumor cells to avoid tumor recognition by cytotoxic T lymphocytes, and thus favor tumor immune escape. Multiple mechanisms are responsible for these HLA class I alterations. In different epithelial tumors, loss of heterozygosity (LOH) at chromosome region 6p21.3, leading to HLA haplotype loss, occurs in 6–50% of all cases depending on the tumor entity. In this paper we report the frequency of LOH at 6p21 in 95 colorectal carcinomas (CRC) previously analyzed for altered HLA class I expression with immunohistological techniques. We used PCR microsatellite amplification of selected STR markers located on Chromosome 6 to identify LOH with DNA from microdissected tumor tissues and the surrounding stroma. Sequence-specific oligonucleotide analysis was performed in microdissected stroma and tumor cells for HLA typing, and to detect HLA haplotype loss. A high frequency (40%) of HLA haplotype loss was found in CRC. Eight tumors showed microsatellite instability. We sometimes observed two or more mechanisms responsible for HLA alteration within the same HLA-altered phenotype, such as LOH and HLA class I total loss. In 25 tumors (26%) no HLA class I alteration could be identified. These data are potentially relevant for CRC patients undergoing T-cell-based immunotherapy.     

A statistical framework for modeling HLA-dependent T cell response data

The identification of T cell epitopes and their HLA (human leukocyte antigen) restrictions is important for applications such as the design of cellular vaccines for HIV. Traditional methods for such identification are costly and time-consuming. Recently, a more expeditious laboratory technique using ELISpot assays has been developed that allows for rapid screening of specific responses. However, this assay does not directly provide information concerning the HLA restriction of a response, a critical piece of information for vaccine design. Thus, we introduce, apply, and validate a statistical model for identifying HLA-restricted epitopes from ELISpot data. By looking at patterns across a broad range of donors, in conjunction with our statistical model, we can determine (probabilistically) which of the HLA alleles are likely to be responsible for the observed reactivities. Additionally, we can provide a good estimate of the number of false positives generated by our analysis (i.e., the false discovery rate). This model allows us to learn about new HLA-restricted epitopes from ELISpot data in an efficient, cost-effective, and high-throughput manner. We applied our approach to data from donors infected with HIV and identified many potential new HLA restrictions. Among 134 such predictions, six were confirmed in the lab and the remainder could not be ruled as invalid. These results shed light on the extent of HLA class I promiscuity, which has significant implications for the understanding of HLA class I antigen presentation and vaccine development.     

HLA Typing from 1000 Genomes Whole Genome and Whole Exome Illumina Data

Specific HLA genotypes are known to be linked to either resistance or susceptibility to certain diseases or sensitivity to certain drugs. In addition, high accuracy HLA typing is crucial for organ and bone marrow transplantation. The most widespread high resolution HLA typing method used to date is Sanger sequencing based typing (SBT), and next generation sequencing (NGS) based HLA typing is just starting to be adopted as a higher throughput, lower cost alternative. By HLA typing the HapMap subset of the public 1000 Genomes paired Illumina data, we demonstrate that HLA-A, B and C typing is possible from exome sequencing samples with higher than 90% accuracy. The older 1000 Genomes whole genome sequencing read sets are less reliable and generally unsuitable for the purpose of HLA typing. We also propose using coverage % (the extent of exons covered) as a quality check (QC) measure to increase reliability.     

HLA class I antigen and HLA-A, -B, and -C haplotype frequencies in Uruguayans

HLA class I antigens were determined for 959 unrelated Uruguayans. The predominant HLA alleles were A2, Cw4, and B35, and the most frequently observed two-loci haplotypes were A2-B44 and B35-Cw4. The most frequent three-loci HLA haplotype was A2-Cw5-B44. We compared the Uruguayan sample with similar data from other populations.     

HLA homozygous stem cell lines derived from human parthenogenetic blastocysts

Individual HLA homozygous parthenogenetic human stem cell (hpSC-Hhom) lines have the potential for cell-based therapy in a significant number of individuals, provided the HLA haplotype is prevalent. We report the successful derivation of four stable hpSC-Hhom lines from both HLA homozygous and HLA heterozygous donors. Of these, the hpSC-Hhom-4 line carries the HLA haplotype found most commonly within the U.S. population, and is shared by different racial groups. These hpSC-Hhom lines demonstrate typical human embryonic stem cell morphology, expressing appropriate stem cell markers and possessing high levels of alkaline phosphatase and telomerase activity. Additionally, injection of these cell lines into immunodeficient animals leads to teratoma formation. G-banded karyotyping demonstrates a normal 46,XX karyotype in lines hpSC-Hhom-1 and hpSC-Hhom-4, and chromosomal anomalies in lines hpSC-Hhom-2 and hpSC-Hhom-3, both derived from the same donor. HLA genotyping of all four hpSC-Hhom lines demonstrates that they are HLA homozygous. Furthermore, in the case of HLA heterozygous donors, the hpSC-Hhom lines inherit the haplotype from only one of the donor's parents. Single-nucleotide polymorphism (SNP) data analysis suggests that hpSC-Hhom lines derived from HLA heterozygous oocyte donors are homozygous throughout the genome as assessed by SNP analysis. The protocol used for deriving these HLA homozygous stem cell lines minimizes the use of animal-derived components, which makes them more appealing for potential clinical application.     

Cutting edge: resistance to HIV-1 infection among African female sex workers is associated with inhibitory KIR in the absence of their HLA ligands

NK cells are regulated in part by killer Ig-like receptors (KIR) that interact with HLA molecules on potential target cells. KIR and HLA loci are highly polymorphic and certain KIR/HLA combinations were found to protect against HIV disease progression. We show in this study that KIR/HLA interactions also influence resistance to HIV transmission. HIV-exposed but seronegative female sex workers in Abidjan, C?te d'Ivoire, frequently possessed inhibitory KIR genes in the absence of their cognate HLA genes: KIR2DL2/KIR2DL3 heterozygosity in the absence of HLA-C1 and KIR3DL1 homozygosity in the absence of HLA-Bw4. HIV-seropositive female sex workers were characterized by corresponding inhibitory KIR/HLA pairings: KIR2DL3 homozygosity together with HLA-C1 and a trend toward KIR3DL1/HLA-Bw4 homozygosity. Absence of ligands for inhibitory KIR could lower the threshold for NK cell activation. In addition, exposed seronegatives more frequently possessed AB KIR genotypes, which contain more activating KIR. The data support an important role for NK cells and KIR/HLA interactions in antiviral immunity.     

Data‐Independent Acquisition of HLA Class I Peptidomes on the Q Exactive Mass Spectrometer Platform

The characterization of peptides presented by human leukocyte antigen (HLA) class I molecules is crucial for understanding immune processes, biomarker discovery, and the development of novel immunotherapies or vaccines. Mass spectrometry allows the direct identification of thousands of HLA‐bound peptides from cell lines, blood, or tissue. In recent years, data‐independent acquisition (DIA) mass spectrometry methods have evolved, promising to increase reproducibility and sensitivity over classical data‐dependent acquisition (DDA) workflows. Here, we describe a DIA setup on the Q Exactive mass spectrometer, optimized regarding the unique properties of HLA class I peptides. The methodology enables sensitive and highly reproducible characterization of HLA peptidomes from individual cell lines. From up to 16 DDA analyses of 100 million human cells, more than 10 000 peptides could be confidently identified, serving as basis for the generation of spectral libraries. This knowledge enabled the subsequent interrogation of DIA data, leading to the identification of peptide sets with >90% overlap between replicate samples, a prerequisite for the comparative study of closely related specimens. Furthermore, >3000 peptides could be identified from just one million cells after DIA analysis using a library generated from 300 million cells. The reduction in sample quantity and the high reproducibility of DIA‐based HLA peptidome analysis should facilitate personalized medicine applications.     

中国壮族和裕固族群体HLAⅠ类区域Alu插入多态性及其与HLA-A位点的相关性

近年来研究发现: 位于HLAⅠ类基因区域的Alu插入是研究不同群体HLAⅠ类基因区域祖先单倍型和HLAⅠ类基因多样性产生、进化和重组的理想工具。文章对中国壮族和裕固族群体HLAⅠ类基因区域5个Alu插入多态性(AluMICB、AluTF、AluHJ、AluHG和AluHF)进行研究, 结合HLA基因分型数据, 分析壮族、裕固族、哈尼族、布朗族和傣族5个民族群体中Alu插入与HLA-A等位基因的关系。研究结果显示: (1)壮族和裕固族人群中5个Alu插入频率范围分别为1.5%~35.8%和9.2~34.8%, AluMICB、AluTF和AluHF插入频率在这两个群体中有统计学差异(P<0.05); (2)在5个研究的群体中, AluHG插入与HLA-A*02的不同亚型关联; AluHJ插入与HLA-A*2402在5个群体中都关联, 但AluHJ与HLA-A*1101和HLA-A*2407只在布朗族中关联。表明不同群体HLAⅠ类基因区域内Alu插入具有各自的特征, 且Alu插入与不同的HLA-A等位基因相关联。这种Alu插入及其与HLA-A的关联特征可作为研究群体中HLAⅠ类基因和单倍型系谱变化的重要遗传标记。