Dataset Information

HLA-VBSeq: accurate HLA typing at full resolution from whole-genome sequencing data.

ABSTRACT: Human leucocyte antigen (HLA) genes play an important role in determining the outcome of organ transplantation and are linked to many human diseases. Because of the diversity and polymorphisms of HLA loci, HLA typing at high resolution is challenging even with whole-genome sequencing data.We have developed a computational tool, HLA-VBSeq, to estimate the most probable HLA alleles at full (8-digit) resolution from whole-genome sequence data. HLA-VBSeq simultaneously optimizes read alignments to HLA allele sequences and abundance of reads on HLA alleles by variational Bayesian inference. We show the effectiveness of the proposed method over other methods through the analysis of predicting HLA types for HLA class I (HLA-A, -B and -C) and class II (HLA-DQA1,-DQB1 and -DRB1) loci from the simulation data of various depth of coverage, and real sequencing data of human trio samples.HLA-VBSeq is an efficient and accurate HLA typing method using high-throughput sequencing data without the need of primer design for HLA loci. Moreover, it does not assume any prior knowledge about HLA allele frequencies, and hence HLA-VBSeq is broadly applicable to human samples obtained from a genetically diverse population.

SUBMITTER: Nariai N

PROVIDER: S-EPMC4331721 | biostudies-literature | 2015

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

HLA-VBSeq: accurate HLA typing at full resolution from whole-genome sequencing data.

Nariai Naoki N Kojima Kaname K Saito Sakae S Mimori Takahiro T Sato Yukuto Y Kawai Yosuke Y Yamaguchi-Kabata Yumi Y Yasuda Jun J Nagasaki Masao M

BMC genomics 20150121

<h4>Background</h4>Human leucocyte antigen (HLA) genes play an important role in determining the outcome of organ transplantation and are linked to many human diseases. Because of the diversity and polymorphisms of HLA loci, HLA typing at high resolution is challenging even with whole-genome sequencing data.<h4>Results</h4>We have developed a computational tool, HLA-VBSeq, to estimate the most probable HLA alleles at full (8-digit) resolution from whole-genome sequence data. HLA-VBSeq simultaneo ...[more]

PMID: 25708870

Similar Datasets

Project description:The incidence of narcolepsy type 1 (NT1) increased in Sweden following the 2009-2010 mass-vaccination with the influenza Pandemrix-vaccine. NT1 has been associated with Human leukocyte antigen (HLA) DQB1*06:02 but full high-resolution HLA-typing of all loci in vaccine-induced NT1 remains to be done. Therefore, here we performed HLA typing by sequencing HLA-DRB3, DRB4, DRB5, DRB1, DQA1, DQB1, DPA1 and DPB1 in 31 vaccine-associated NT1 patients and 66 of their first-degree relatives (FDR), and compared these data to 636 Swedish general population controls (GP). Previously reported disease-related alleles in the HLA-DRB5*01:01:01-DRB1*15:01:01-DQA1*01:02:01-DQB1*06:02:01 extended haplotype were increased in NT1 patients (34/62 haplotypes, 54.8%) compared to GP (194/1272 haplotypes, 15.3%, p = 6.17E-16). Indeed, this extended haplotype was found in 30/31 patients (96.8%) and 178/636 GP (28.0%). In total, 15 alleles, four extended haplotypes, and six genotypes were found to be increased or decreased in frequency among NT1 patients compared to GP. Among subjects with the HLA-DRB5*01:01:01-DRB1*15:01:01-DQA1*01:02-DQB1*06:02 haplotype, a second DRB4*01:03:01-DRB1*04:01:01-DQA1*03:02//*03:03:01-DQB1*03:01:01 haplotype (p = 2.02E-2), but not homozygosity for DRB1*15:01:01-DQB1*06:02:01 (p = 7.49E-1) conferred association to NT1. Alleles with increased frequency in DQA1*01:02:01 (p = 1.07E-2) and DQA1*03:02//*03:03:01 (p = 3.26E-2), as well as with decreased frequency in DRB3*01:01:02 (p = 8.09E-3), DRB1*03:01:01 (p = 1.40E-2), and DQB1*02:01:01 (p = 1.40E-2) were found among patients compared to their FDR. High-resolution HLA sequencing in Pandemrix-associated NT1 confirmed the strong association with the DQB1*06:02:01-containing haplotype but also revealed an increased association to the not previously reported extended HLA-DRB4*01:03:01-DRB1*04:01:01-DQA1*03:02//*03:03:01-DQB1*03:01:01 haplotype. High-resolution HLA typing should prove useful in dissecting the immunological mechanisms of vaccination-associated NT1.

Project description:The MLST scheme currently used for Enterococcus faecium typing was designed in 2002 and is based on putative gene functions and Enterococcus faecalis gene sequences available at that time. As a result, the original MLST scheme does not correspond to the real genetic relatedness of E. faecium strains and often clusters genetically distant strains to the same sequence types (ST). Nevertheless, typing has a significant impact on the subsequent epidemiological conclusions and introduction of appropriate epidemiological measures, thus it is crucial to use a more accurate MLST scheme. Based on the genome analysis of 1,843 E. faecium isolates, a new scheme, consisting of 8 highly discriminative loci, was created in this study. These strains were divided into 421 STs using the new MLST scheme, as opposed to 223 STs assigned by the original MLST scheme. The proposed MLST has a discriminatory power of D = 0.983 (CI95% 0.981 to 0.984), compared to the original scheme's D = 0.919 (CI95% 0.911 to 0.927). Moreover, we identified new clonal complexes with our newly designed MLST scheme. The scheme proposed here is available within the PubMLST database. Although whole genome sequencing availability has increased rapidly, MLST remains an integral part of clinical epidemiology, mainly due to its high standardization and excellent robustness. In this study, we proposed and validated a new MLST scheme for E. faecium, which is based on genome-wide data and thus reflects the tested isolates' more accurate genetic similarity. IMPORTANCE Enterococcus faecium is one of the most important pathogens causing health care associated infections. One of the main reasons for its clinical importance is a rapidly spreading resistance to vancomycin and linezolid, which significantly complicates antibiotic treatment of infections caused by such resistant strains. Monitoring the spread and relationships between resistant strains causing severe conditions represents an important tool for implementing appropriate preventive measures. Therefore, there is an urgent need to establish a robust method enabling strain monitoring and comparison at the local, national, and global level. Unfortunately, the current, extensively used MLST scheme does not reflect the real genetic relatedness between individual strains and thus does not provide sufficient discriminatory power. This can lead directly to incorrect epidemiological measures due to insufficient accuracy and biased results.

Dataset Information

HLA-VBSeq: accurate HLA typing at full resolution from whole-genome sequencing data.

Publications

HLA-VBSeq: accurate HLA typing at full resolution from whole-genome sequencing data.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets