Cost-efficient multiplex PCR for routine genotyping of up to nine classical HLA loci in a single analytical run of multiple samples by next generation sequencing.

BACKGROUND: HLA genotyping by next generation sequencing (NGS) requires three basic steps, PCR, NGS, and allele assignment. Compared to the conventional methods, such as PCR-sequence specific oligonucleotide primers (SSOP) and -sequence based typing (SBT), PCR-NGS is extremely labor intensive and time consuming. In order to simplify and accelerate the NGS-based HLA genotyping method for multiple DNA samples, we developed and evaluated four multiplex PCR methods for genotyping up to nine classical HLA loci including HLA-A, HLA-B, HLA-C, HLA-DRB1/3/4/5, HLA-DQB1, and HLA-DPB1. RESULTS: We developed multiplex PCR methods using newly and previously designed middle ranged PCR primer sets for genotyping different combinations of HLA loci, (1) HLA-DRB1/3/4/5, (2) HLA-DQB1 (3.8 kb to 5.3 kb), (3) HLA-A, HLA-B, HLA-C, and (4) HLA-DPB1 (4.6 kb to 7.2 kb). The primer sets were designed to genotype polymorphic exons to the field 3 level or 6-digit typing. When we evaluated the PCR method for genotyping all nine HLA loci (9LOCI) using 46 Japanese reference subjects who represented a distribution of more than 99.5% of the HLA alleles at each of the nine HLA loci, all of the 276 alleles genotyped, except for HLA-DRB3/4/5 alleles, were consistent with known alleles assigned by the conventional methods together with relevant locus balance and no excessive allelic imbalance. One multiplex PCR method (9LOCI) was able to provide precise genotyping data even when only 1 ng of genomic DNA was used for the PCR as a sample template. CONCLUSIONS: In this study, we have demonstrated that the multiplex PCR approach for NGS-based HLA genotyping could serve as an alternative routine HLA genotyping method, possibly replacing the conventional methods by providing an accelerated yet robust amplification step. The method also could provide significant merits for clinical applications with its ability to amplify lower quantity of samples and the cost-saving factors.

Evolution and single-nucleotide polymorphisms in methicillin-resistant Staphylococcus aureus strains with reduced susceptibility to vancomycin and daptomycin, based on determination of the complete genome.

We obtained a series of methicillin-resistant Staphylococcus aureus isolates, including both daptomycin-susceptible strain TD1 and daptomycin-resistant strain TD4, from a patient. We determined the complete genome sequences of TD1 and TD4 using next-generation sequencing, and only four single-nucleotide polymorphisms (SNPs) were identified, one each in capB (E58K), rpoB (H481Y), lytN (I16V), and mprF (V351E). We determined that these four SNPs were sufficient to cause the strains to develop daptomycin, vancomycin, and rifampin resistance.

Discovery of novel MHC-class I alleles and haplotypes in Filipino cynomolgus macaques (Macaca fascicularis) by pyrosequencing and Sanger sequencing: Mafa-class I polymorphism.

Although the low polymorphism of the major histocompatibility complex (MHC) transplantation genes in the Filipino cynomolgus macaque (Macaca fascicularis) is expected to have important implications in the selection and breeding of animals for medical research, detailed polymorphism information is still lacking for many of the duplicated class I genes. To better elucidate the degree and types of MHC polymorphisms and haplotypes in the Filipino macaque population, we genotyped 127 unrelated animals by the Sanger sequencing method and high-resolution pyrosequencing and identified 112 different alleles, 28 at cynomolgus macaque MHC (Mafa)-A, 54 at Mafa-B, 12 at Mafa-I, 11 at Mafa-E, and seven at Mafa-F alleles, of which 56 were newly described. Of them, the newly discovered Mafa-A8*01:01 lineage allele had low nucleotide similarities (<86%) with primate MHC class I genes, and it was also conserved in the Vietnamese and Indonesian populations. In addition, haplotype estimations revealed 17 Mafa-A, 23 Mafa-B, and 12 Mafa-E haplotypes integrated with 84 Mafa-class I haplotypes and Mafa-F alleles. Of these, the two Mafa-class I haplotypes, F/A/E/B-Hp1 and F/A/E/B-Hp2, had the highest haplotype frequencies at 10.6 and 10.2%, respectively. This suggests that large scale genetic screening of the Filipino macaque population would identify these and other high-frequency Mafa-class I haplotypes that could be used as MHC control animals for the benefit of biomedical research.

Association of HLA-DRB1*15:02:01, DQB1*05:01:24 and DPB1*13:01:01 in Thai patients with systemic sclerosis.

HLA studies in patients with systemic sclerosis (SSc) have shown variable results. This study aimed to examine the association of HLA class I and II risk alleles in Thai SSc patients, and clarify the contribution of risk HLA alleles to the pathogenesis and clinical manifestations. Blood samples from 92 SSc patients and 135 healthy controls (HCs) were collected. Eleven loci of the HLA class I (HLA-A, B, and C) and class II (HLA-DR, DP, and DQ) genes were determined by a 3-field (6-digit) analysis using the Next Generation DNA Sequencing (NGS) method. Anti-topoisomerase-I antibodies (ATA) and anti-centromere antibodies (ACA) were identified by ELISA methods. Allele frequencies (AFs) of HLA-DRB1*15:02:01, DRB5*01:02:01, DQB1*05:01:24, DPB1*13:01:01, and DQA1*01:01:01 were increased significantly in the whole SSc and SSc patients with positive ATA, but with negative ACA (SSc/ATA+/ACA-). Of these, DPB1*13:01:01 was the most susceptible allele. The DRB1*15:02:01, DQB1:05:01:24, and DPB1*13:01:01 alleles were estimated to locate on the unique haplotype, and haplotype frequency was estimated to be significantly higher than those in the HCs (p = 0.002). The linkage analysis of DRB1*15/16 revealed that most of the DRB1*15:02:01 alleles were linked to DRB5*01:02:01 or DRB5*01:08:01N. The linkage of DRB1*16:02:01 to DRB5*01:01:01 was observed frequently. The associations of risk alleles with several SSc clinical features were observed. HLA-DRB1*15:02:01, DRB5*01:02:01, DQB1*05:01:24, and DPB1*13:01:01 on the unique haplotype were associated with the pathogenesis and clinical features of SSc in Thai patients. The linkage of DRB1*15:02:01 to DRB5*01:08:01N was observed commonly in northern Thai patients.

A new chromosome 14-based human artificial chromosome (HAC) vector system for efficient transgene expression in human primary cells.

The use of non-integrating human artificial chromosomes (HACs) in gene therapy possibly allows for safe and reliable genetic modification of human cells without insertional mutagenesis and/or unexpected oncogene activations. Although we previously demonstrated that the HAC provides long-term therapeutic erythropoietin (EPO) production in normal human primary fibroblasts (hPFs), the expression level of EPO was too low to provide medical benefits for human therapy. Thus, the next challenge for the application of this system in therapeutic purposes is to improve the transgene expression on HACs. Here, we newly constructed chromosome 14-based HACs and examined the effects of the telomere and promoter regions on the expression level of the tansgene in hPFs. We showed that the use of natural telomere/sub-telomere and enhancers within the 5' untranslated region of the human ubiquitin C gene greatly increased (over 1000-fold) the EPO production in hPFs. Furthermore, we demonstrated the reprogramming of mouse embryonic fibroblasts by HAC-mediated introduction of four transcription factors, and established induced pluripotent stem cells with no trace of the HACs carrying multiple expression cassettes with large genome fragments. These results indicate that this HAC system could allow us to manipulate multiple transgenes efficiently in human primary cells, providing a promising tool not only for gene therapy but also for investigating genome functions in drug discoveries.

Identification of HLA-A*02:06:01 as the primary disease susceptibility HLA allele in cold medicine-related Stevens-Johnson syndrome with severe ocular complications by high-resolution NGS-based HLA typing.

Stevens-Johnson syndrome (SJS) and toxic epidermal necrolysis (TEN) are life-threatening acute inflammatory vesiculobullous reactions of the skin and mucous membranes. These severe cutaneous drug reactions are known to be caused by inciting drugs and infectious agents. Previously, we have reported the association of HLA-A*02:06 and HLA-B*44:03 with cold medicine (CM)-related SJS/TEN with severe ocular complications (SOCs) in the Japanese population. However, the conventional HLA typing method (PCR-SSOP) sometimes has ambiguity in the final HLA allele determination. In this study, we performed HLA-disease association studies in CM-SJS/TEN with SOCs at 3- or 4-field level. 120 CM-SJS/TEN patients with SOCs and 817 Japanese healthy controls are HLA genotyped using the high-resolution next-generation sequencing (NGS)-based HLA typing of HLA class I genes, including HLA-A, HLA-B, and HLA-C. Among the alleles of HLA class I genes, HLA-A*02:06:01 was strongly associated with susceptibility to CM-SJS/TEN (p = 1.15 × 10-18, odds ratio = 5.46). Four other alleles (HLA-A*24:02:01, HLA-B*52:01:01, HLA-B*46:01:01, and HLA-C*12:02:02) also demonstrated significant associations. HLA haplotype analyses indicated that HLA-A*02:06:01 is primarily associated with susceptibility to CM-SJS/TEN with SOCs. Notably, there were no specific disease-causing rare variants among the high-risk HLA alleles. This study highlights the importance of higher resolution HLA typing in the study of disease susceptibility, which may help to elucidate the pathogenesis of CM-SJS/TEN with SOCs.

Reference Grade Characterization of Polymorphisms in Full-Length HLA Class I and II Genes With Short-Read Sequencing on the ION PGM System and Long-Reads Generated by Single Molecule, Real-Time Sequencing on the PacBio Platform.

Although NGS technologies fuel advances in high-throughput HLA genotyping methods for identification and classification of HLA genes to assist with precision medicine efforts in disease and transplantation, the efficiency of these methods are impeded by the absence of adequately-characterized high-frequency HLA allele reference sequence databases for the highly polymorphic HLA gene system. Here, we report on producing a comprehensive collection of full-length HLA allele sequences for eight classical HLA loci found in the Japanese population. We augmented the second-generation short read data generated by the Ion Torrent technology with long amplicon spanning consensus reads delivered by the third-generation SMRT sequencing method to create reference grade high-quality sequences of HLA class I and II gene alleles resolved at the genomic coding and non-coding level. Forty-six DNAs were obtained from a reference set used previously to establish the HLA allele frequency data in Japanese subjects. The samples included alleles with a collective allele frequency in the Japanese population of more than 99.2%. The HLA loci were independently amplified by long-range PCR using previously designed HLA-locus specific primers and subsequently sequenced using SMRT and Ion PGM sequencers. The mapped long and short-reads were used to produce a reference library of consensus HLA allelic sequences with the help of the reference-aware software tool LAA for SMRT Sequencing. A total of 253 distinct alleles were determined for 46 healthy subjects. Of them, 137 were novel alleles: 101 SNVs and/or indels and 36 extended alleles at a partial or full-length level. Comparing the HLA sequences from the perspective of nucleotide diversity revealed that HLA-DRB1 was the most divergent among the eight HLA genes, and that the HLA-DPB1 gene sequences diverged into two distinct groups, DP2 and DP5, with evidence of independent polymorphisms generated in exon 2. We also identified two specific intronic variations in HLA-DRB1 that might be involved in rheumatoid arthritis. In conclusion, full-length HLA allele sequencing by third-generation and second-generation technologies has provided polymorphic gene reference sequences at a genomic allelic resolution including allelic variations assigned up to the field-4 level for a stronger foundation in precision medicine and HLA-related disease and transplantation studies.

Phenotypic difference between Δ(srl-recA)306 and ΔrecA::Km elucidated by next-generation sequencing combined with a long-PCR system.

Many significant gene mutations in E. coli have contributed to the development of genetics. Among these, a commonly used recA mutation, Δ(srl-recA)306 has been sequenced by a next-generation sequencer combined with a long PCR. An original report described that Δ(srl-recA)306 cells were deleted from srlR to recA genes in their genome. The next-generation sequencer enables more accurate details to be determined. We ask whether both surrounding genes from hypF to norV for srlR and alaS for recA is there first. The long PCR was carried out with primers, norR and alaS, and amplified DNA fragments differed in length from wild to Δ(srl-recA)306 cells, suggesting that an entire Δ(srl-recA)306 mutation was included. Sequences of those DNA fragments indicated that 9147 bp, from srlR to recA including 10 genes, were replaced by a Tn10 DNA sequence. Junction points at both srlR-Tn10 and Tn10-recA were determined precisely. The results indicate that the first 97% of recA gene sequences were lost with a downstream recX gene remaining intact. The phenotypic difference between Δ(srl-recA)306 and ΔrecA::Km is discussed.