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

INTRODUCTION

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.

METHODS

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.

RESULTS

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.

CONCLUSION

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. This article is protected by copyright. All rights reserved.

Decoding the diversity of killer immunoglobulin-like receptors by deep sequencing and a high-resolution imputation method

The killer cell immunoglobulin-like receptor (KIR) recognizes human leukocyte antigen (HLA) class I molecules and modulates the function of natural killer cells. Despite its role in immunity, the complex genomic structure has limited a deep understanding of the KIR genomic landscape. Here we conduct deep sequencing of 16 KIR genes in 1,173 individuals. We devise a bioinformatics pipeline incorporating copy number estimation and insertion or deletion (indel) calling for high-resolution KIR genotyping. We define 118 alleles in 13 genes and demonstrate a linkage disequilibrium structure within and across KIR centromeric and telomeric regions. We construct a KIR imputation reference panel (n_reference = 689, imputation accuracy = 99.7%), apply it to biobank genotype (n_total = 169,907), and perform phenome-wide association studies of 85 traits. We observe a dearth of genome-wide significant associations, even in immune traits implicated previously to be associated with KIR (the smallest p = 1.5 × 10^-4). Our pipeline presents a broadly applicable framework to evaluate innate immunity in large-scale datasets.

Proteogenomic identification of an immunogenic HLA class I neoantigen in mismatch repair-deficient colorectal cancer tissue

Although CD8⁺ T cells recognize neoantigens that arise from somatic mutations in cancer, only a small fraction of nonsynonymous mutations give rise to clinically relevant neoantigens. In this study, HLA class I ligandomes of a panel of human colorectal cancer (CRC) and matched normal tissues were analyzed using mass spectrometry–based proteogenomic analysis. Neoantigen presentation was rare; however, the analysis detected a single neoantigen in a mismatch repair–deficient CRC (dMMR-CRC) tissue sample carrying 3967 nonsynonymous mutations, where abundant tumor-infiltrating lymphocytes (TILs) and inflamed gene expression status were observed in the tumor microenvironment (TME). Using the HLA class I ligandome data and gene expression profiles, a set of nonmutated tumor-associated antigen (TAA) candidates was concomitantly identified. Interestingly, CD8⁺ TILs predominantly recognized the detected neoantigen over the array of TAA candidates. Neoantigen-reactive CD8⁺ TILs showed PD-1 positivity and exhibited functional and specific responses. Moreover, T cell receptor (TCR) profiling identified the sequence of the neoantigen-reactive TCR clonotype and showed its expansion in the TME. Transduction of the sequenced TCR conferred neoantigen specificity and cytotoxicity to peripheral blood lymphocytes. The proteogenomic approach revealed the antigenic and reactive T cell landscape in dMMR-CRC, demonstrating the presence of an immunogenic neoantigen and its potential therapeutic applications.

Association of Bovine Leukemia Virus-Induced Lymphoma with BoLA-DRB3 Polymorphisms at DNA, Amino Acid, and Binding Pocket Property Levels

Bovine leukemia virus (BLV) causes enzootic bovine leucosis, a malignant B-cell lymphoma in cattle. The DNA sequence polymorphisms of bovine leukocyte antigen (BoLA)-DRB3 have exhibited a correlation with BLV-induced lymphoma in Holstein cows. However, the association may vary between different cattle breeds. Furthermore, little is known about the relationship between BLV-induced lymphoma and DRB3 at the amino acid and structural diversity levels. Here, we comprehensively analyzed the correlation between BLV-induced lymphoma and DRB3 at DNA, amino acid, and binding pocket property levels, using 106 BLV-infected asymptomatic and 227 BLV-induced lymphoma Japanese black cattle samples. DRB3*011:01 was identified as a resistance allele, whereas DRB3*005:02 and DRB3*016:01 were susceptibility alleles. Amino acid association studies showed that positions 9, 11, 13, 26, 30, 47, 57, 70, 71, 74, 78, and 86 were associated with lymphoma susceptibility. Structure and electrostatic charge modeling further indicated that binding pocket 9 of resistance DRB3 was positively charged. In contrast, alleles susceptible to lymphoma were neutrally charged. Altogether, this is the first association study of BoLA-DRB3 polymorphisms with BLV-induced lymphoma in Japanese black cattle. In addition, our results further contribute to understanding the mechanisms regarding how BoLA-DRB3 polymorphisms mediate susceptibility to BLV-induced lymphoma.

Genetic control of CCL24, POR, and IL23R contributes to the pathogenesis of sarcoidosis

Sarcoidosis is a genetically complex systemic inflammatory disease that affects multiple organs. We present a GWAS of a Japanese cohort (700 sarcoidosis cases and 886 controls) with replication in independent samples from Japan (931 cases and 1,042 controls) and the Czech Republic (265 cases and 264 controls). We identified three loci outside the HLA complex, CCL24, STYXL1-SRRM3, and C1orf141-IL23R, which showed genome-wide significant associations (P < 5.0 × 10⁻⁸) with sarcoidosis; CCL24 and STYXL1-SRRM3 were novel. The disease-risk alleles in CCL24 and IL23R were associated with reduced CCL24 and IL23R expression, respectively. The disease-risk allele in STYXL1-SRRM3 was associated with elevated POR expression. These results suggest that genetic control of CCL24, POR, and IL23R expression contribute to the pathogenesis of sarcoidosis. We speculate that the CCL24 risk allele might be involved in a polarized Th1 response in sarcoidosis, and that POR and IL23R risk alleles may lead to diminished host defense against sarcoidosis pathogens.

Akira Meguro et al. report a genome-wide association study for sarcoidosis—a systemic inflammatory disease—in the Japanese population. They identify 3 non-HLA loci with genome-wide significance, 2 of which have not been previously associated with sarcoidosis in any population.

Genome-Wide Association Study in Asians Identifies Novel Loci for High Myopia and Highlights a Nervous System Role in Its Pathogenesis

PURPOSE

To identify novel susceptibility loci for high myopia.

DESIGN

Genome-wide association study (GWAS) followed by replication and meta-analysis.

PARTICIPANTS

A total of 14 096 samples from East and Southeast Asian populations (2549 patients with high myopia and 11 547 healthy controls).

METHODS

We performed a GWAS in 3269 Japanese individuals (1668 with high myopia and 1601 control participants), followed by replication analysis in a total of 10 827 additional samples (881 with high myopia and 9946 control participants) from Japan, Singapore, and Taiwan. To confirm the biological role of the identified loci in the pathogenesis of high myopia, we performed functional annotation and Gene Ontology (GO) analyses.

MAIN OUTCOME MEASURES

We evaluated the association of single nucleotide polymorphisms with high myopia and GO terms enriched among genes identified in the current study.

RESULTS

We identified 9 loci with genome-wide significance (P < 5.0 × 10^-8). Three loci were previously reported myopia-related loci (ZC3H11B on 1q41, GJD2 on 15q14, and RASGRF1 on 15q25.1), and the other 6 were novel (HIVEP3 on 1p34.2, NFASC/CNTN2 on 1q32.1, CNTN4/CNTN6 on 3p26.3, FRMD4B on 3p14.1, LINC02418 on 12q24.33, and AKAP13 on 15q25.3). The GO analysis revealed a significant role of the nervous system related to synaptic signaling, neuronal development, and Ras/Rho signaling in the pathogenesis of high myopia.

CONCLUSIONS

The current study identified 6 novel loci associated with high myopia and demonstrated an important role of the nervous system in the disease pathogenesis. Our findings give new insight into the genetic factors underlying myopia, including high myopia, by connecting previous findings and allowing for a clarified interpretation of the cause and pathophysiologic features of myopia at the molecular level.

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⁻¹⁸, 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.

Identification of HLA-DRB1*04:10 allele as risk allele for Japanese moyamoya disease and its association with autoimmune thyroid disease: A case-control study

Background and purpose

Moyamoya disease (MMD) is a progressive cerebrovascular disease with unknown etiology. Growing evidence suggest its involvement of autoimmune and genetic mechanisms in the pathogenesis of MMD. This study aims to clarify the association between HLA allele and MMD.

Methods

Case-control study: the DNA of 136 MMD patients in Japan was extracted and the genotype of human leukocyte antigen (HLA) from this DNA was determined by super-high-resolution single-molecule sequence-based typing using next-generation sequencing. Next, the frequency of each HLA allele (HLA-A, HLA-B, HLA-C, HLA-DRB1, HLA-DQB1, and HLA-DPB1) was compared with those in the Japanese control database. In addition, haplotype estimation was performed using the expectation maximization algorithm.

Results

The frequencies of the HLA-DRB1*04:10 allele (4.77% vs. 1.47% in the control group; P = 1.7 × 10⁻³; odds ratio [OR] = 3.35) and of the HLA-DRB1*04:10–HLA-DQB1*04:02 haplotype (haplotype frequency 4.41% vs. 1.35% in the control group; P = 2.0 × 10⁻³; OR = 3.37) significantly increased. The frequency of thyroid diseases, such as Graves’ disease and Hashimoto thyroiditis, increased in HLA-DRB1*04:10-positive MMD patients compared with that in HLA-DRB1*04:10-negative MMD patients.

Conclusions

HLA-DRB1*04:10 is a risk allele and HLA-DRB1*04:10–HLA-DQB1*04:02 a risk haplotype for MMD. In addition, HLA-DRB1*04:10 is associated with thyroid disease in MMD patients.

Analysis of STAT1, STAT2 and STAT3 mRNA expression levels in the blood of patients with multiple sclerosis

BACKGROUND

Multiple sclerosis (MS) is the most common chronic, inflammatory, autoimmune disease of the central nervous system (CNS) maintained by the secretion of a large number of cytokines [1]. The signal transducer and activator of transcription (STAT) family has an essential role in transmitting many of the cytokine-mediated signals and failure in the signaling process contributes to the etiopathogenesis of MS.

METHODS

This study aimed to assess STAT1, STAT2 and STAT3 gene expression in the blood of 50 relapsing-remitting MS (RR-MS) patients and 50 healthy controls by TaqMan Quantitative Real-Time PCR.

RESULTS

The results showed that STAT1 gene expression was significantly up-regulated (p= 0.023), whereas STAT2 gene expression was significantly down-regulated (p< 0.0001) in MS patients compared to controls. On the other hand, there was no significant difference between MS patients and controls for STAT3 gene expression (p= 0.837). In addition, there was no significant correlation between the expression of STAT1, STAT2, STAT3 genes and clinical findings, such as the level of physical disability in MS patients (according to the Kurtzke Expanded Disability Status Scale (EDSS) criterion) and disease duration.

CONCLUSION

A significant positive correlation was demonstrated between STAT1 and STAT2 and also between STAT1 and STAT3. This study shows for the first time that a comparison of the relative quantitative expression of three different STAT genes in the blood cells of MS patients compared to controls revealed marked differences in the expression of the STAT family genes that might reflect their different roles in the pathogenesis of MS. These transcripts might be useful biomarkers for evaluating the efficacy of IFN treatment of the MS patients.

Genome-wide association study identifies eight risk loci and implicates metabo-psychiatric origins for anorexia nervosa

Characterized primarily by a low body-mass index, anorexia nervosa is a complex and serious illness1, affecting 0.9-4% of women and 0.3% of men2-4, with twin-based heritability estimates of 50-60%5. Mortality rates are higher than those in other psychiatric disorders6, and outcomes are unacceptably poor7. Here we combine data from the Anorexia Nervosa Genetics Initiative (ANGI)8,9 and the Eating Disorders Working Group of the Psychiatric Genomics Consortium (PGC-ED) and conduct a genome-wide association study of 16,992 cases of anorexia nervosa and 55,525 controls, identifying eight significant loci. The genetic architecture of anorexia nervosa mirrors its clinical presentation, showing significant genetic correlations with psychiatric disorders, physical activity, and metabolic (including glycemic), lipid and anthropometric traits, independent of the effects of common variants associated with body-mass index. These results further encourage a reconceptualization of anorexia nervosa as a metabo-psychiatric disorder. Elucidating the metabolic component is a critical direction for future research, and paying attention to both psychiatric and metabolic components may be key to improving outcomes.

Identification and Structure of an MHC Class I-Encoded Protein with the Potential to Present N-Myristoylated 4-mer Peptides to T Cells

Similar to host proteins, N-myristoylation occurs for viral proteins to dictate their pathological function. However, this lipid-modifying reaction creates a novel class of "lipopeptide" Ags targeted by host CTLs. The primate MHC class I-encoded protein, Mamu-B*098, was previously shown to bind N-myristoylated 5-mer peptides. Nevertheless, T cells exist that recognize even shorter lipopeptides, and much remains to be elucidated concerning the molecular mechanisms of lipopeptide presentation. We, in this study, demonstrate that the MHC class I allele, Mamu-B*05104, binds the N-myristoylated 4-mer peptide (C14-Gly-Gly-Ala-Ile) derived from the viral Nef protein for its presentation to CTLs. A phylogenetic tree analysis indicates that these classical MHC class I alleles are not closely associated; however, the high-resolution x-ray crystallographic analyses indicate that both molecules share lipid-binding structures defined by the exceptionally large, hydrophobic B pocket to accommodate the acylated glycine (G1) as an anchor. The C-terminal isoleucine (I4) of C14-Gly-Gly-Ala-Ile anchors at the F pocket, which is distinct from that of Mamu-B*098 and is virtually identical to that of the peptide-presenting MHC class I molecule, HLA-B51. The two central amino acid residues (G2 and A3) are only exposed externally for recognition by T cells, and the methyl side chain on A3 constitutes a major T cell epitope, underscoring that the epitopic diversity is highly limited for lipopeptides as compared with that for MHC class I-presented long peptides. These structural features suggest that lipopeptide-presenting MHC class I alleles comprise a distinct MHC class I subset that mediates an alternative pathway for CTL activation.

Association between human leucocyte antigen alleles and risk of stroke in Iranian population

Previous studies have demonstrated associations between human leucocyte antigen (HLA) and some types of ischaemic stroke. In the present study, we genotyped HLA-A,-B and -DRB1 alleles in 140 Iranian patients with history of ischaemic stroke and 140 age-/sex-matched healthy subjects. No significant difference has been found in the distribution of HLA-A and B alleles between cases and controls. The DRB1*16 allele was significantly over-represented in patient group compared with control group (Adjusted p value = 0.048). Other HLA-DRB1 alleles were not associated with stroke risk. The HLA-B*35,B*52 genotype was significantly more prevalent among patients compared with controls (Adjusted p value = 0.03, OR [95% CI] = 9.3 [1.3, 407.2]). Several HLA haplotypes were associated with risk of stroke in the assessed population. The current study provides further evidences for participation of HLA in conferring risk of ischaemic stroke.

Super High Resolution for Single Molecule-Sequence-Based Typing of Classical HLA Loci Using Ion Torrent PGM

Super high resolution-single molecule-sequence-based typing (SS-SBT) is an HLA DNA typing method to the field 4 level of allelic resolution (formerly known as 8-digit typing) to efficiently detect novel and null alleles without phase ambiguity by combination of long ranged PCR amplification and next-generation sequencing (NGS) technologies. In this chapter, we describe three basic steps, long ranged PCR, NGS, and allele assignment.

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.

Association of HLA alleles with autism

BACKGROUND

Autism spectrum disorders (ASD) are a group of heterogeneous neurodevelopmental disorders known by impaired social interaction and activities and abnormal repetitive behavior. As a multifactorial disorder, several genetic and immunological factors have been shown to be implicated in its pathogenesis.

METHODS

Among them are certain human leukocyte antigen (HLA) alleles. In the current study, we genotyped HLA-A, -B & DRB alleles in 103 Iranian ASD patients and 180 age, gender, and ethnic-matched healthy controls.

RESULTS

After Boferroni correction no allele or haplotype was associated with genetic susceptibility to ASD in Iranian population.

CONCLUSION

Future studies are needed to assess contribution of immunological factors such as HLA alleles in ASD pathogenesis.

A new MHC-linked susceptibility locus for primary Sjögren's syndrome: MICA

The association of primary Sjögren's syndrome (pSS) with Major Histocompatibility Complex (MHC) alleles is quintessential of MHC-disease associations. Indeed, although disease associations with classical HLA class I and II alleles/haplotypes are amply documented, further dissection is often prevented by the strong linkage disequilibrium across the entire MHC complex. Here we study the association of pSS, not with HLA genes, but with the non-conventional MHC encoded class I gene, MICA (MHC class I chain-related gene A). MICA is selectively expressed within epithelia, and is the major ligand for the activatory receptor, NKG2D, both attributes relevant to pSS' etiology. MICA-pSS association was studied in two independent (French and UK) cohorts representing a total of 959 cases and 1,043 controls. MICA*008 allele was shown to be significantly associated with pSS (pcor=2.61 × 10-35). A multivariate logistic regression showed that this association was independent of all major known MHC-linked risk loci/alleles, as well as other relevant candidate loci that are in linkage disequilibrium with MICA*008 i.e. HLA-B*08:01, rs3131619 (T), MICB*008, TNF308A, HLA-DRB1*03:01 and HLA-DRB1*15:01 (P = 1.84 × 10-04). Furthermore, independently of the MICA*008 allele, higher levels of soluble MICA proteins were detected in sera of pSS patients compared to healthy controls. This study hence defines MICA as a new, MHC-linked, yet HLA-independent, pSS risk locus and opens a new front in our understanding of the still enigmatic pathophysiology of this disease. The fact that the soluble MICA protein is further amplified in MICA*008 carrying individuals, might also be relevant in other auto-immune diseases and cancer.

Significant Locus and Metabolic Genetic Correlations Revealed in Genome-Wide Association Study of Anorexia Nervosa

OBJECTIVE

The authors conducted a genome-wide association study of anorexia nervosa and calculated genetic correlations with a series of psychiatric, educational, and metabolic phenotypes.

METHOD

Following uniform quality control and imputation procedures using the 1000 Genomes Project (phase 3) in 12 case-control cohorts comprising 3,495 anorexia nervosa cases and 10,982 controls, the authors performed standard association analysis followed by a meta-analysis across cohorts. Linkage disequilibrium score regression was used to calculate genome-wide common variant heritability (single-nucleotide polymorphism [SNP]-based heritability [h²_SNP]), partitioned heritability, and genetic correlations (r_g) between anorexia nervosa and 159 other phenotypes.

RESULTS

Results were obtained for 10,641,224 SNPs and insertion-deletion variants with minor allele frequencies >1% and imputation quality scores >0.6. The h²_SNP of anorexia nervosa was 0.20 (SE=0.02), suggesting that a substantial fraction of the twin-based heritability arises from common genetic variation. The authors identified one genome-wide significant locus on chromosome 12 (rs4622308) in a region harboring a previously reported type 1 diabetes and autoimmune disorder locus. Significant positive genetic correlations were observed between anorexia nervosa and schizophrenia, neuroticism, educational attainment, and high-density lipoprotein cholesterol, and significant negative genetic correlations were observed between anorexia nervosa and body mass index, insulin, glucose, and lipid phenotypes.

CONCLUSIONS

Anorexia nervosa is a complex heritable phenotype for which this study has uncovered the first genome-wide significant locus. Anorexia nervosa also has large and significant genetic correlations with both psychiatric phenotypes and metabolic traits. The study results encourage a reconceptualization of this frequently lethal disorder as one with both psychiatric and metabolic etiology.

Interleukin 7 Receptor Alpha Gene Variants Are Correlated with Gene Expression in Patients with Relapsing-remitting Multiple Sclerosis

The association of single nucleotide polymorphisms (SNPs) of the IL-7Rα gene with multiple sclerosis (MS) have been documented in various populations. This study aimed to evaluate the genotype distributions of two SNPs, rs6897932 and rs201084372, and the functional association of rs6897932 in relation to IL-7Rα gene expression in a group of Iranian relapsing-remitting MS (RRMS) patients.  Genotyping for both SNPs in the IL7Rα gene and relative quantification of mRNA expression for both isoforms of IL-7Rα were performed in 100 RRMS patients and 100 ethnic-matched healthy controls. Higher significant frequencies of the T allele and TT genotype for rs6897932 (C/T) were observed in patients comparing to controls (p=0.006). Higher frequencies of the T allele and the TT and TG genotypes and lower frequencies of the G allele and GG genotypes for rs201084372 (G/A) were found in patients comparing to controls (p<0.0001).  A decreased level of mRNA expression for the membrane-bound IL-7Rα (mbIL-7Rα) and an increased level of mRNA for the soluble IL-7Rα (sIL-7Rα) were observed in patients versus controls (p=0.005 and p=0.002 respectively). A significant decreased level of mRNA expression for mbIL-7Rα (p=0.01) and an increased level of mRNA for sIL-7Rα (p=0.008) were observed in RRMS patients compared to healthy controls carrying the TT+CT genotypes. The higher levels of mRNA expression for the sIL-7Rα isoform in MS patients carrying the IL7R*TT genotype is a new finding not previously reported in studies on the genotype-induced effects of IL-7Rα expression in multiple sclerosis.

Single nucleotide polymorphisms in the FOXP3 gene are associated with increased risk of relapsing-remitting multiple sclerosis

BACKGROUND

Although Multiple Sclerosis (MS) is an autoimmune multifactorial disease with unknown etiology, various genetic and environmental factors are known to contribute to the pathogenesis of the disease.

OBJECTIVE

Recent studies have confirmed that the suppressive function of regulatory T cells (T (reg)) is impaired in MS patients and that the FOXP3 gene is a crucial transcription factor in the regulation of CD4+CD25+FOXP3+ Treg cells. Polymorphisms in the promoter region of the FOXP3 gene may alter the gene expression level and, therefore, contribute to the disease susceptibility.

METHODS

The present study aimed to investigate the possible association between single nucleotide polymorphisms (SNPs) rs3761548 and rs2232365 in the FOXP3 gene and predisposition to MS. We conducted a case-control study on 410 patients with sporadic MS and 446 healthy controls. Genotyping was performed using polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP).

RESULTS

Significant differences in distribution of both rs3761548 and rs2232365 A allele were found in MS patients in comparison to controls. Haplotype frequencies were also different among the studied groups. The A-A and C-G haplotype blocks showed a significant difference between case and controls.

CONCLUSION

we have provided further evidence for the association between genetic variations and haplotypes in FOXP3 and MS in Iranian population.

Inverse Relation between MxA Gene Expression and Age in Multiple Sclerosis Patients Reveals a Gender Difference in Response to Interferon Therapy

Multiple sclerosis (MS) is an inflammatory, multifocal, immune-mediated disease of the central nervous system that women are at a higher risk to acquire than men. Myxovirus resistance protein A (MxA) is used as a predictive marker of bioactivity of interferon-beta (IFN-β) therapy in MS patients. This study was undertaken in west of Iran to investigate gender differences in the expression level of MxA in relapsing-remitting MS (RRMS) patients receiving IFN-β therapy, compared with untreated normal individuals. The expression level of the MxA gene in RRMS samples were compared to untreated normal individuals using the extracted RNA from whole blood of 50 RRMS patients (31 females and 19 males) and 50 normal controls (29 females and 21 males). All patients were HLA-DRB1*15 negative and responded to IFN-β with a normal vitamin D level. The level of MxA gene expression was measured by quantitative RT-PCR. The levels of gene expression were decreased in RRMS patients compared with normal counterparts (p=0.025). This decrease was significant in females (p=0.009) compared to males (p>0.05). The level of expression varied across different female age-groups with no significant difference in women younger than 30 years, but a significant decrease in expression in women between 30 to 40 years or above 40 years of age was seen. There was neither linear correlation between the MxA expression level and risk of expanded disability status scale of Kurtzke (EDSS); nor were there any significant correlation between expression status of MxA and duration of the disease. In conclusion, the decrease in the level of MxA expression in MS patients treated with IFN-β when compared to normal individuals was significantly lower in females than males.  This demonstrated a gender bias in the response to IFN-β therapy that will need to be confirmed and further investigated in more detail.

Production of a Locus- and Allele-Specific Monoclonal Antibody for the Characterization of SLA-1*0401 mRNA and Protein Expression Levels in MHC-Defined Microminipigs

The class I major histocompatibility complex (MHC) presents self-developed peptides to specific T cells to induce cytotoxity against infection. The MHC proteins are encoded by multiple loci that express numerous alleles to preserve the variability of the antigen-presenting ability in each species. The mechanism regulating MHC mRNA and protein expression at each locus is difficult to analyze because of the structural and sequence similarities between alleles. In this study, we examined the correlation between the mRNA and surface protein expression of swine leukocyte antigen (SLA)-1*0401 after the stimulation of peripheral blood mononuclear cells (PBMCs) by Staphylococcus aureus superantigen toxic shock syndrome toxin-1 (TSST-1). We prepared a monoclonal antibody (mAb) against a domain composed of Y102, L103 and L109 in the α2 domain. The Hp-16.0 haplotype swine possess only SLA-1*0401, which has the mAb epitope, while other haplotypes possess 0 to 3 SLA classical class I loci with the mAb epitopes. When PBMCs from SLA-1*0401 homozygous pigs were stimulated, the SLA-1*0401 mRNA expression level increased until 24 hrs and decreased at 48 hrs. The kinetics of the interferon regulatory transcription factor-1 (IRF-1) mRNA level were similar to those of the SLA-1*0401 mRNA. However, the surface protein expression level continued to increase until 72 hrs. Similar results were observed in the Hp-10.0 pigs with three mAb epitopes. These results suggest that TSST-1 stimulation induced both mRNA and surface protein expression of class I SLA in the swine PBMCs differentially and that the surface protein level was sustained independently of mRNA regulation.

Carboxypeptidase E, Identified As a Direct Interactor of Growth Hormone, Is Important for Efficient Secretion of the Hormone

We have identified 88 interactor candidates for human growth hormone (GH) by the yeast two-hybrid assay. Among those, we focused our efforts on carboxypeptidase E (CPE), which has been thought to play a key role in sorting prohormones, such as pro-opiomelanocortin (POMC), to regulated secretory vesicles. We found that CPE co-localizes with and interacts with GH in AtT20 pituitary cells. Downregulation of CPE led to decreased levels of GH secretion, consistent with involvement of CPE in GH sorting/secretion. Our binding assay in vitro with bacterially expressed proteins suggested that GH directly interacts with CPE but in a manner different from POMC.

RAR-related orphan receptor A (RORA): A new susceptibility gene for multiple sclerosis

Retinoic acid receptor-related orphan receptor alpha (RORA) is proposed to promote Th17 cells differentiation that play a crucial role in many inflammatory diseases, including multiple sclerosis (MS). The gene is also involved in regulation of inflammatory responses and neuronal cell development. The aim of the present study is to determine if any relation exists between RORA rs11639084 and rs4774388 gene polymorphisms on the individual susceptibility of multiple sclerosis. 410 patients with clinically definite MS and 500 ethnically-matched healthy controls participated in this study. Genotyping was performed using tetra primer-amplification refractory mutation system-PCR (4P-ARMS-PCR) method for the mentioned polymorphisms in the RORA gene. Both variants showed significant differences in allele and genotype distributions between the studied groups. Genotypes were risk associated in additive (P-value of 0.0003 and odds ratio equal to 1.7 (95% CI: 1.27-2.26)), dominant (P-value of <0.0001 and odds ratio equal to 0.55 (95% CI: 0.41-0.73)) and recessive (P-value of 0.04 and odds ratio equal to 0.33 (95% CI: (0.12-0.96)) models for rs11639084. However, the rs4774388 genotypes were risk associated in recessive model with a P-value of 0.036 and an odds ratio of 0.62 (95% CI: (0.4-0.97)). To the best of our knowledge this is the first report concerning the association between RORΑ gene polymorphisms and MS. The further study of RORΑ related pathways and gene networks might result in the better understanding of the pathophysiology of MS and related symptoms.

Comparative genomics of the human, macaque and mouse major histocompatibility complex

The MHC is a highly polymorphic genomic region that encodes the transplantation and immune regulatory molecules. It receives special attention for genetic investigation because of its important role in the regulation of innate and adaptive immune responses and its strong association with numerous infectious and/or autoimmune diseases. The MHC locus was first discovered in the mouse and for the past 50 years it has been studied most intensively in both mice and humans. However, in recent years the macaque species have emerged as some of the more important and advanced experimental animal models for biomedical research into MHC with important human immunodeficiency virus/simian immunodeficiency virus and transplantation studies undertaken in association with precise MHC genotyping and haplotyping methods using Sanger sequencing and next-generation sequencing. Here, in this special issue on 'Macaque Immunology' we provide a short review of the genomic similarities and differences among the human, macaque and mouse MHC class I and class II regions, with an emphasis on the association of the macaque class I region with MHC polymorphism, haplotype structure and function.

HLA alleles and haplotypes in Burmese (Myanmarese) and Karen in Thailand

This is the first report on human leukocyte antigen (HLA) allele and haplotype frequencies at three class I loci and two class II loci in unrelated healthy individuals from two ethnic groups, 170 Burmese and 200 Karen, originally from Burma (Myanmar), but sampled while residing in Thailand. Overall, the HLA allele and haplotype frequencies detected by polymerase chain reaction sequence-specific primer (PCR-SSP) at five loci (A, B, C, DRB1 and DRQB1) at low resolution showed distinct differences between the Burmese and Karen. In Burmese, five HLA-B*15 haplotypes with different HLA-A and HLA-DR/DQ combinations were detected with three of these not previously reported in other Asian populations. The data are important in the fields of anthropology, transplantation and disease-association studies.

HLA genes as modifiers of response to IFN-β-1a therapy in relapsing-remitting multiple sclerosis

AIMS

This study investigated the influence of HLA class-I and -II genes in the response to IFN-β in relapsing-remitting multiple sclerosis (MS) patients.

PATIENTS & METHODS

In this cohort, 231 relapsing-remitting MS patients who are classified into IFN-β responders (n = 146) and nonresponders (n = 85) and 180 ethnic-matched healthy controls were analyzed. Clinical outcome of IFN-β therapy particularly Expanded Disability Status Scale scores were evaluated in relation to HLA-A, -B and -DRB1 alleles and haplotypes.

RESULTS

Increased frequencies of HLA-DRB1*04 allele and HLA-A*03-B*44-DRB1*04 haplotype, and decreased frequency of HLA-B*15 were associated with better response to IFN-β treatment.

CONCLUSION

The possibility of genetic screening particularly HLA typing prior to starting IFN-β therapy for MS may permit the identification of likely responders or nonresponders.

SLC1A1 Gene Variants and Normal Tension Glaucoma: An Association Study

BACKGROUND

It has been hypothesized that dysfunction of the solute carrier family 1, member1 gene (SLC1A1), which encodes the glutamate aspartate transporter, may play a role in normal tension glaucoma. In this study we investigate whether SLC1A1 is associated with normal tension glaucoma in Japanese patients.

METHODS

A total of 292 Japanese patients with normal tension glaucoma and 500 healthy control subjects were recruited. We genotyped 12 single-nucleotide polymorphisms in SLC1A1. We also performed an imputation analysis to evaluate the potential association of un-genotyped SLC1A1 single-nucleotide polymorphisms, and 165 single-nucleotide polymorphisms were imputed.

RESULTS

We observed an increased frequency of the G allele of rs10739062 in patients compared to controls (p = 0.043, OR = 1.25). The rs10739062 polymorphism exhibited a dominant effect: individuals with genotype GG and GC showed a 1.91-fold increase in risk compared to genotype CC (p = 0.0082). However, the statistical significance disappeared after Bonferroni correction for multiple testing (pc > 0.05). We did not find any significant association between any of the remaining 176 single-nucleotide polymorphisms and disease risk.

CONCLUSIONS

Our study showed a lack of association between SLC1A1 variants and normal tension glaucoma in Japanese patients, suggesting that the SLC1A1 gene does not play a critical role in the development of the disorder in this patient population. However, further genetic studies with larger sample sizes are needed to clarify whether SLC1A1 may make some contribution that affects the risk of developing normal tension glaucoma.

Crystal structure of the N-myristoylated lipopeptide-bound MHC class I complex

The covalent conjugation of a 14-carbon saturated fatty acid (myristic acid) to the amino-terminal glycine residue is critical for some viral proteins to function. This protein lipidation modification, termed N-myristoylation, is targeted by host cytotoxic T lymphocytes (CTLs) that specifically recognize N-myristoylated short peptides; however, the molecular mechanisms underlying lipopeptide antigen (Ag) presentation remain elusive. Here we show that a primate major histocompatibility complex (MHC) class I-encoded protein is capable of binding N-myristoylated 5-mer peptides and presenting them to specific CTLs. A high-resolution X-ray crystallographic analysis of the MHC class I:lipopeptide complex reveals an Ag-binding groove that is elaborately constructed to bind N-myristoylated short peptides rather than prototypic 9-mer peptides. The identification of lipopeptide-specific, MHC class I-restricted CTLs indicates that the widely accepted concept of MHC class I-mediated presentation of long peptides to CTLs may need some modifications to incorporate a novel MHC class I function of lipopeptide Ag presentation.

Lipid antigens have been added to the antigenic repertoire in recent years. Here, the authors have identified Mamu-B*098 as a lipopeptide antigen presenting molecule and using structural and biochemical analysis have shown that it has a different antigen binding pocket to previously analysed proteins.

Increased expression ratio of matrix metalloproteinase-9 (MMP9) and tissue inhibitor of matrix metalloproteinase (TIMP-1) RNA levels in Iranian multiple sclerosis patients

OBJECTIVES

Multiple sclerosis (MS) is an autoimmune disease involving the central nervous system (CNS) with unknown immunopathogenic mechanisms. Matrix metalloproteinase-9 (MMP-9) facilitates T-cell migration into the CNS while the tissue inhibitor matrix metalloproteinase-1 (TIMP-1) inhibits MMP-9 actions. The aim of this study was to evaluate the expression of TIMP-1 RNA and MMP-9/TIMP-1 RNA ratio in blood cells of Iranian patients with relapsing-remitting multiple sclerosis (RRMS) treated with IFNb.

MATERIAL AND METHODS

The study compared the expression level of TIMP-1 gene in RRMS samples with normal individuals in Iran and the results were compared using a ratio of MMP-9 to TIMP-1. All patients were HLA-DRB1*15 negative and were responders to interferon-beta with a normal vitamin D level.

RESULTS

The RRMS patients manifested a lower expression level of TIMP-1 RNA than their normal counterparts although the result was not significant (P= 0.06). Also, the ratio of MMP-9 to TIMP-1 RNA increased significantly (P= 0.009). There was no linear correlation between TIMP-1 expression level and risk of Expanded Disability Status Scale of Kurtzke (EDSS); nor was there any significant correlation between expression status of TIMP-1 and duration of the disease. Although there was no significant decrease in TIMP-1 expression level, the MMP-9/TIMP-1 RNA ratio in RRMS was significantly higher than normal subjects.

CONCLUSION

Further studies are recommended to compare MMP-9/TIMP-1 RNA ratio in patients before and after taking IFN-beta in order to find out if MMP-9/TIMP-1 RNA ratio can function as a proper marker of the bio efficacy of IFN-beta treatment of MS.

Differentiation ability of multipotent hematopoietic stem/progenitor cells detected by a porcine specific anti-CD117 monoclonal antibody

CD117 is a cytokine receptor expressed on the surface of hematopoietic stem cells with a likely role in cell survival, proliferation and differentiation. In order to study the differentiation activity of porcine CD117 hematopoietic cells in vitro and in vivo we prepared an anti-swine CD117 Mab (2A1) with high specificity for flow-cytometrical analysis. The 2A1 Mab did not recognize mouse or human mast cells suggesting that 2A1 is species-specific. Swine bone marrow (BM) CD117+ cells differentiated in vitro mainly into erythroid and monocyte lineages in the methylcellulose-based colony assay. When the swine BM CD117+ cells were transplanted in vivo into immunodeficient NOG (NOD/SCID/IL-2gc-null) mice, a significant amount of swine CD45+ leukocytes, including CD3 positive T cells, were developed in the mice. These results revealed that the swine BM CD117+ cells possess hematopoietic stem/progenitor activity and when monitored in immunodeficient mice or in vitro they can develop into lymphoid, erythroid, and myeloid cells efficiently with the new monoclonal antibody.

HLA-C*12:02 is a susceptibility factor in late-onset type of psoriasis in Japanese

Psoriasis is thought to be a multifactorial disease triggered by both genetic and environmental factors. The HLA-C locus on chromosome 6p21.33 remains the strongest susceptibility candidate locus in psoriasis. The strong association between psoriasis and the HLA-Cw6 allele has been well documented in various races. It is known that psoriatic patients with early onset are more likely to be familial and associated with HLA-Cw6. Familial occurrence of Japanese psoriasis is smaller than other populations. Furthermore, males are predominant over females in Japanese psoriasis. We investigated the relation between HLA-C alleles and age of onset, and in each gender for Japanese psoriasis, and discuss male predominance in the incidence of psoriasis in Japan. Four hundred forty six unrelated Japanese patients with psoriasis vulgaris and 557 sex- and age-matched unrelated Japanese healthy controls were investigated by genotyping. We confirmed the association between early-onset type of psoriasis with HLA-C*06:02 allele in Japanese. In addition, we detected the association between the late-onset type of psoriasis and the HLA-C*12:02 allele in Japanese. No significant differences in allele frequency were observed between females and males. Our results suggest that there is no genetic factor effect on male predominance in Japanese. In contract, the effect of environmental risk factors on the onset of Japanese psoriatic patients is stronger in males than in females. As a result, male predominant in psoriasis may occur in Japan.

Establishment of immortalized mesenchymal stem cells derived from the submandibular glands of tdTomato transgenic mice

Transgenic mice that overexpress the red fluorescent protein tdTomato (tdTomato mice) are well suited for use in regenerative medicine studies. Cultured cells from this murine model exhibit strong red fluorescence, enabling real-time in vivo imaging through the body surface of grafted animals. Mesenchymal stem cells (MSCs) have marked potential for use in cell therapy and regenerative medicine; however, the mechanisms that regulate their dynamics in vivo are poorly understood. In the present study, an MSC line was derived from the submandibular gland fibroblasts of tdTomato mice. The fluorescent signal from this cell line was observed in organs throughout the body, as well as in salivary glands. Primary culture cells derived from the submandibular gland were immortalized with SV40 large T antigen (GManSV cells); these cells exhibited increased migratory ability, as compared with those isolated from the sublingual gland. GManSV cells were tdTomato-positive and exhibited spindle-shaped fibroblastic morphology; they also robustly expressed mouse MSC markers: Stem cell antigen-1 (Sca-1), CD44, and CD90. This cell line retained multipotent stem cell characteristics, as evidenced by its ability to differentiate into both osteogenic and adipogenic lineages. These results indicate that Sca-1⁺/CD44⁺/CD90⁺-GManSV cells may be useful for kinetic studies of submandibular gland-derived MSCs in the context of in vitro co-culture with other types of salivary gland-derived cells. These cells may also be used for in vivo imaging studies, in order to identify novel cell therapy and regenerative medicine for the treatment of salivary gland diseases.

Assessment of Artificial MiRNA Architectures for Higher Knockdown Efficiencies without the Undesired Effects in Mice

RNAi-based strategies have been used for hypomorphic analyses. However, there are technical challenges to achieve robust, reproducible knockdown effect. Here we examined the artificial microRNA (amiRNA) architectures that could provide higher knockdown efficiencies. Using transient and stable transfection assays in cells, we found that simple amiRNA-expression cassettes, that did not contain a marker gene (−MG), displayed higher amiRNA expression and more efficient knockdown than those that contained a marker gene (+MG). Further, we tested this phenomenon in vivo, by analyzing amiRNA-expressing mice that were produced by the pronuclear injection-based targeted transgenesis (PITT) method. While we observed significant silencing of the target gene (eGFP) in +MG hemizygous mice, obtaining −MG amiRNA expression mice, even hemizygotes, was difficult and the animals died perinatally. We obtained only mosaic mice having both “−MG amiRNA” cells and “amiRNA low-expression” cells but they exhibited growth retardation and cataracts, and they could not transmit the –MG amiRNA allele to the next generation. Furthermore, +MG amiRNA homozygotes could not be obtained. These results suggested that excessive amiRNAs transcribed by −MG expression cassettes cause deleterious effects in mice, and the amiRNA expression level in hemizygous +MG amiRNA mice is near the upper limit, where mice can develop normally. In conclusion, the PITT-(+MG amiRNA) system demonstrated here can generate knockdown mouse models that reliably express highest and tolerable levels of amiRNAs.

Genome-wide association study of IgA nephropathy using 23 465 microsatellite markers in a Japanese population

Immunoglobulin A nephropathy (IgAN) is the most common form of primary glomerulonephritis in many parts of the world. Although previous genome-wide association studies (GWAS) identified the major susceptibility loci for IgAN, the causal genes currently remain unknown. We performed a GWAS using 23 465 microsatellite (MS) markers to identify genes related to IgAN in a Japanese population. A pooled sample analysis was conducted in three-stage screenings of three independent case-control populations, and after the final step of individual typing, 11 markers survived. Of these, we focused on two regions on 6p21 and 12q21 because they (i) showed the strongest relationship with IgAN, and (ii) appeared to be highly relevant to IgAN in view of several previous studies. These regions contained the HLA, TSPAN8 and PTPRR genes. This study on GWAS, using >20 000 MS markers, provides a new approach regarding susceptible genes for IgAN for investigators seeking new tools for the prevention and treatment of IgAN.

On the genetics of the Silk Route: association analysis of HLA, IL10, and IL23R-IL12RB2 regions with Behçet's disease in an Iranian population

Despite that the association of Behçet's disease (BD) with the HLA-B5 was first established in the 1970s, a number of recent genome-wide association studies have both confirmed and furthered this association--in various populations--to individual SNPs both inside and outside the HLA. The former include HLA-B, MICA, and HLA-A, while the latter encompass IL10 and IL23R-IL12RB2 regions. The present study examined whether some of these SNPs could be replicated in an Iranian population, where the prevalence of disease is amply documented. Eight SNPs were selected and tested in 552 patients and 417 controls. These were rs7539328, rs12119179, rs1495965, rs1518111, and rs1800871 in IL10 and IL23R-IL12RB2 regions and rs114854070, rs12525170, and rs76546355 (formerly rs116799036) in the HLA locus. The well-known BD-associated genes HLA-B and MICA were independently genotyped. Although we were not able to formally replicate the association with IL10 and IL23R-IL12RB2, we do report that BD in Iran is strongly associated with HLA-B*51, MICA-A6, and the three HLA-linked SNPs (odds ratio (OR) = 3.38, P = 6.21 × 10(-14); OR = 2.08, P = 1.58 × 10(-13); and OR = 1.67-4.05, P = 1.45 × 10(-04) to 4.79 × 10(-34), respectively). Our data further indicate that the robust HLA-B/MICA association may be explained by a single variant (rs76546355) between the two genes. Overall, these data contribute to a better appraisal of the intriguing linkage between BD and the ancient Silk Route, spanning from the Mediterranean shores to Japan.

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.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1514-4) contains supplementary material, which is available to authorized users.

No contribution of the ABCB11 p.444A polymorphism in Japanese patients with drug-induced cholestasis

European studies have revealed that the ABCB11 c.1331T>C (V444A) polymorphism (rs2287622) C-allele frequency is higher among patients with drug-induced cholestasis. Given the low incidence of this disease, however, this association has not been sufficiently elucidated. We aimed to investigate the significance of this polymorphism in Japanese patients. We determined ABCB11 V444A polymorphism frequencies and HLA genotypes in two independent drug-induced cholestasis cohorts. Expression and taurocholate transport activity of proteins from 444A variants were analyzed using Madin-Darby canine kidney II cells. In cohort 1 (n = 40), the V444A polymorphism C-allele frequency (66%) was lower than that in controls (n = 190, 78%), but this difference was not significant (P = 0.09). In cohort 2 (n = 119), comprising patients with cholestatic (n = 19), hepatocellular (n = 74), and mixed (n = 26) liver injuries, the C-allele frequency was lower among patients with cholestatic liver injury (68%) than among those with hepatocellular (75%) or mixed liver injury (83%), although this difference was not significant. In cohort 1, HLA-A*0201 was observed more frequently in patients (22%) than in controls [11%; P = 0.003; odds ratio, 2.4 (95% confidence interval, 1.4-4.0)]. Taurocholate transport activity of 444A-encoded protein was significantly lower than that of 444V-encoded protein (81% of 444V, P < 0.05) because of the reduced protein stability. In conclusion, ABCB11 444A had slightly reduced transport activity, but it did not contribute to the occurrence of drug-induced cholestasis in Japanese patients. Therefore, genetic susceptibility to acquired cholestasis may differ considerably by ethnicity.

Association of sick building syndrome with neuropathy target esterase (NTE) activity in Japanese

Sick building syndrome (SBS) is a set of several clinically recognizable symptoms reported by occupants of a building without a clear cause. Neuropathy target esterase (NTE) is a membrane bound serine esterase and its reaction with organophosphates (OPs) can lead to OP-induced delayed neuropathy (OPIDN) and nerve axon degeneration. The aim of our study was to determine whether there was a difference in NTE activity in the peripheral blood mononuclear cells (PBMCs) of Japanese patients with SBS and healthy controls and whether PNPLA6 (alias NTE) gene polymorphisms were associated with SBS. We found that the enzymatic activity of NTE was significantly higher (P < 0.0005) in SBS patients compared with controls. Moreover, population with an AA genotype of a single nucleotide polymorphism (SNP), rs480208, in intron 21 of the PNPLA6 gene strongly reduced the activity of NTE. Fifty-eight SNP markers within the PNPLA6 gene were tested for association in a case-control study of 188 affected individuals and 401 age-matched controls. Only one SNP, rs480208, was statistically different in genotype distribution (P = 0.005) and allele frequency (P = 0.006) between the cases and controls (uncorrected for testing multiple SNP sites), but these were not significant by multiple corrections. The findings of the association between the enzymatic activity of NTE and SBS in Japanese show for the first time that NTE activity might be involved with SBS.

A genome-wide association study of anorexia nervosa

Anorexia nervosa (AN) is a complex and heritable eating disorder characterized by dangerously low body weight. Neither candidate gene studies nor an initial genome-wide association study (GWAS) have yielded significant and replicated results. We performed a GWAS in 2907 cases with AN from 14 countries (15 sites) and 14 860 ancestrally matched controls as part of the Genetic Consortium for AN (GCAN) and the Wellcome Trust Case Control Consortium 3 (WTCCC3). Individual association analyses were conducted in each stratum and meta-analyzed across all 15 discovery data sets. Seventy-six (72 independent) single nucleotide polymorphisms were taken forward for in silico (two data sets) or de novo (13 data sets) replication genotyping in 2677 independent AN cases and 8629 European ancestry controls along with 458 AN cases and 421 controls from Japan. The final global meta-analysis across discovery and replication data sets comprised 5551 AN cases and 21 080 controls. AN subtype analyses (1606 AN restricting; 1445 AN binge-purge) were performed. No findings reached genome-wide significance. Two intronic variants were suggestively associated: rs9839776 (P=3.01 × 10(-7)) in SOX2OT and rs17030795 (P=5.84 × 10(-6)) in PPP3CA. Two additional signals were specific to Europeans: rs1523921 (P=5.76 × 10(-)(6)) between CUL3 and FAM124B and rs1886797 (P=8.05 × 10(-)(6)) near SPATA13. Comparing discovery with replication results, 76% of the effects were in the same direction, an observation highly unlikely to be due to chance (P=4 × 10(-6)), strongly suggesting that true findings exist but our sample, the largest yet reported, was underpowered for their detection. The accrual of large genotyped AN case-control samples should be an immediate priority for the field.