Resolving incomplete single nucleotide polymorphism tagging of HLA-DQ2.2 for coeliac disease genotyping using digital droplet PCR.

A hallmark of coeliac disease (CD) is the exceptionally strong genetic association with HLA-DQ2.5, DQ8, and DQ2.2. HLA typing provides information on CD risk important to both clinicians and researchers. A method that enables simple and fast detection of all CD risk genotypes is particularly desirable for the study of large populations. Single nucleotide polymorphism (SNP)-based HLA typing can detect the CD risk genotypes by detecting a combination of six SNPs but this approach can struggle to resolve HLA-DQ2.2, seen in 4% of European CD patients, because of the low resolution of one negatively predicting SNP. We sought to optimise SNP-based HLA typing by harnessing the additional resolution of digital droplet PCR to resolve HLA-DQ2.2. Here we test this two-step approach in an unselected sample of Mexican DNA and compare its accuracy to DNA typed using traditional exon detection. The addition of digital droplet PCR for samples requiring negative prediction of HLA-DQ2.2 enabled HLA-DQ2.2 to be accurately typed. This technique is a simple addition to a SNP-based typing strategy and enables comprehensive definition of all at-risk HLA genotypes in CD in a timely and cost-effective manner.

Appropriate clinical use of human leukocyte antigen typing for coeliac disease: an Australasian perspective.

The past decade has seen human leukocyte antigen (HLA) typing emerge as a remarkably popular test for the diagnostic work-up of coeliac disease with high patient acceptance. Although limited in its positive predictive value for coeliac disease, the strong disease association with specific HLA genes imparts exceptional negative predictive value to HLA typing, enabling a negative result to exclude coeliac disease confidently. In response to mounting evidence that the clinical use and interpretation of HLA typing often deviates from best practice, this article outlines an evidence-based approach to guide clinically appropriate use of HLA typing, and establishes a reporting template for pathology providers to improve communication of results.

Higher body mass index in adults at diagnosis of the slowly progressive form of type 1 diabetes mellitus is associated with lower risk HLA genes.

We hypothesised that higher body weight, a proposed risk factor for type 1 diabetes mellitus, would be associated with increased penetrance of lower risk genes. In adults at diagnosis of the slowly progressive form of type 1 diabetes mellitus we found that higher body mass index was associated with the absence of the highest risk HLA genes.

Common and well-documented HLA alleles: 2012 update to the CWD catalogue.

We have updated the catalogue of common and well-documented (CWD) human leukocyte antigen (HLA) alleles to reflect current understanding of the prevalence of specific allele sequences. The original CWD catalogue designated 721 alleles at the HLA-A, -B, -C, -DRB1, -DRB3/4/5, -DQA1, -DQB1, and -DPB1 loci in IMGT (IMmunoGeneTics)/HLA Database release 2.15.0 as being CWD. The updated CWD catalogue designates 1122 alleles at the HLA-A, -B, -C, -DRB1, -DRB3/4/5, -DQA1, -DQB1, -DPA1 and -DPB1 loci as being CWD, and represents 14.3% of the HLA alleles in IMGT/HLA Database release 3.9.0. In particular, we identified 415 of these alleles as being 'common' (having known frequencies) and 707 as being 'well-documented' on the basis of ~140,000 sequence-based typing observations and available HLA haplotype data. Using these allele prevalence data, we have also assigned CWD status to specific G and P designations. We identified 147/151 G groups and 290/415 P groups as being CWD. The CWD catalogue will be updated on a regular basis moving forward, and will incorporate changes to the IMGT/HLA Database as well as empirical data from the histocompatibility and immunogenetics community. This version 2.0.0 of the CWD catalogue is available online at cwd.immunogenomics.org, and will be integrated into the Allele Frequencies Net Database, the IMGT/HLA Database and National Marrow Donor Program's bioinformatics web pages.

A community standard for immunogenomic data reporting and analysis: proposal for a STrengthening the REporting of Immunogenomic Studies statement.

Modern high-throughput HLA and KIR typing technologies are generating a wealth of immunogenomic data with the potential to revolutionize the fields of histocompatibility and immune-related disease association and population genetic research, much as SNP-based approaches have revolutionized association research. The STrengthening the REporting of Genetic Association studies (STREGA) statement provides community-based data reporting and analysis standards for genomic disease-association studies, identifying specific areas in which adoption of reporting guidelines can improve the consistent interpretation of genetic studies. While aspects of STREGA can be applied to immunogenomic studies, HLA and KIR research requires additional consideration, as the high levels of polymorphism associated with immunogenomic data pose unique methodological and computational challenges to the synthesis of information across datasets. Here, we outline the principle challenges to consistency in immunogenomic studies, and propose that an immunogenomic-specific analog to the STREGA statement, a STrengthening the REporting of Immunogenomic Studies (STREIS) statement, be developed as part of the 16th International HLA and Immunogenetics Workshop. We propose that STREIS extends at least four of the 22 elements of the STREGA statement to specifically address issues pertinent to immunogenomic data: HLA and KIR nomenclature, data-validation, ambiguity resolution, and the analysis of highly polymorphic genetic systems. As with the STREGA guidelines, the intent behind STREIS is not to dictate the design of immunogenomic studies, but to ensure consistent and transparent reporting of research, facilitating the synthesis of HLA and KIR data across studies.

HLA-DRB1 associations with disease susceptibility and clinical course in Australians with multiple sclerosis.

Human leucocyte antigen (HLA)-DRB1*1501 and other class II alleles influence susceptibility to multiple sclerosis (MS), but their contribution if any to the clinical course of MS remains uncertain. Here, we have investigated DRB1 alleles in a large sample of 1230 Australian MS cases, with some enrichment for subjects with primary progressive (PPMS) disease (n = 246) and 1210 healthy controls. Using logistic regression, we found that DRB1*1501 was strongly associated with risk (P = 7 x 10(-45)), as expected, and after adjusting for DRB1*1501, a predisposing effect was also observed for DRB1*03 (P = 5 x 10(-7)). Individuals homozygous for either DRB1*15 or DRB1*03 were considerably more at risk of MS than heterozygotes and non-carriers. Both the DRB1*04 and the DRB1*01/DRB1*15 genotype combination, respectively, protected against PPMS in comparison to subjects with relapsing disease. Together, these data provide further evidence of heterogeneity at the DRB1 locus and confirm the importance of HLA variants in the phenotypic expression of MS.

Replication of KIAA0350, IL2RA, RPL5 and CD58 as multiple sclerosis susceptibility genes in Australians.

A recent genome-wide association study (GWAS) conducted by the International Multiple Sclerosis Genetics Consortium (IMSGC) identified a number of putative MS susceptibility genes. Here we have performed a replication study in 1134 Australian MS cases and 1265 controls for 17 risk-associated single nucleotide polymorphisms (SNPs) reported by the IMSGC. Of 16 SNPs that passed quality control filters, four, each corresponding to a different non-human leukocyte antigen (HLA) gene, were associated with disease susceptibility: KIAA0350 (rs6498169) P=0.001, IL2RA (rs2104286) P=0.033, RPL5 (rs6604026) P=0.041 and CD58 (rs12044852) P=0.042. There was no association (P=0.58) between rs6897932 in the IL7R gene and the risk of MS. No interactions were detected between the replicated IMSGC SNPs and HLA-DRB1*15, gender, disease course, disease progression or age-at-onset. We used a novel Bayesian approach to estimate the extent to which our data increased or decreased evidence for association with the six most-associated IMSGC loci. These analyses indicated that even modest P-values, such as those reported here, can contribute markedly to the posterior probability of 'true' association in replication studies. In conclusion, these data provide support for the involvement of four non-HLA genes in the pathogenesis of MS, and combined with previous data, increase to genome-wide significance (P=3 x 10(-8)) evidence of an association between KIAA0350 and risk of disease.

SNP mapping and candidate gene sequencing in the class I region of the HLA complex: searching for multiple sclerosis susceptibility genes in Tasmanians.

This study is an extension to previously published work that has linked variation in the human leukocyte antigen (HLA) class I region with susceptibility to multiple sclerosis (MS) in Australians from the Island State of Tasmania. Single nucleotide polymorphism (SNP) mapping was performed on an 865-kb candidate region (D6S1683-D6S265) in 166 Tasmanian MS families, and seven candidate genes [ubiquitin D (UBD), olfactory receptor 2H3 (OR2H3), gamma-aminobutyric acid B receptor 1 (GABBR1), myelin oligodendrocyte glycoprotein (MOG), HLA-F, HLA complex group 4 (HCG4) and HLA-G] were resequenced. SNPs tagging the extended MS susceptibility haplotype were genotyped in an independent sample of 356 Australian MS trios and SNPs in the MOG gene were significantly over-transmitted to MS cases. We identified significant effects on MS susceptibility of HLA-A*2 (OR: 0.51; P = 0.05) and A*3 (OR: 2.85; P = 0.005), and two coding polymorphisms in the MOG gene (V145I: P = 0.01, OR: 2.2; V142L: P = 0.04, OR: 0.45) after full conditioning on HLA-DRB1. We have therefore identified plausible candidates for the causal MS susceptibility allele, and although not conclusive at this stage, our data provide suggestive evidence for multiple class I MS susceptibility genes.

Clinical, laboratory and genetic markers associated with erosions and remission in patients with early inflammatory arthritis: a prospective cohort study.

We investigated the relationship between clinical, laboratory and genetic markers and outcome measures in 159 patients with recent onset of inflammatory arthritis (IA). The majority of patients were managed in community-based rheumatology practice. Median duration of arthritis at baseline was 3 months with median follow-up of 4.0 years (range 0-10). Markers of disease activity and 1987 ACR criteria for rheumatoid arthritis (RA) were estimated every 6 months for the first 2 years and annually thereafter. Presence of shared epitopes (SE) was established by PCR-based method. Main outcome variables were attainment of remission and presence of erosions on X-rays of hands and feet at 3 years. Remission was seen in 34.3% of patients and was independently related to age 60 and older (odds ratio (OR) 3.2; 95% confidence interval (CI), 1.2-8.7) and inversely to the presence of rheumatoid factor (RF) (OR 8.3; 95% CI, 3.2-21.3 for persistent arthritis). Patients with two SE were likely to have persistent arthritis (P=0.006), but this was not significant when corrected for RF. Independent predictors for erosions at 3 years were RF (OR 7.5; 95% CI, 1.9-29.5) and area under the curve for number of swollen joints (OR 1.08; 95% CI, 1.02-1.16). SE status was not predictive of erosions at 3 years (OR 1.6; 95% CI, 0.7-3.7). In univariate analysis, patients possessing DERAA motif on DRB1 were less likely to have erosive disease than without this motif at 4 years (OR 0.21; 95% CI, 0.0-0.9, P=0.037) but this finding was partly explained by adjusting for RF (adjusted OR 0.24; 95% CI 0.04-1.37). In this study of recent onset IA, active disease and RF were associated with poor outcome. Whilst SE did not predict erosive disease, patients with DERAA motif may be protected against erosions whilst the presence of two SE alleles suggests persistence of arthritis.

Identification of a novel DPB1* allele (DPB1*6901) and the occurrence of HLA haplotypes that extend to DPB1.

The sequence of a novel DPB1 allele, DPB1*6901, observed in a Caucasian bone marrow donor phenotype HLA A2; Cw*0501,*1601; B*4402, *4403; DRB1*0401; *07; DQB1*02; *0301; DPB1*0401; *6901, is described. The sequence is consistent with that previously described for DPB1*0601 with the exception of codon 69. The sequence at this codon is consistent with that previously observed only in the DPB1*1101 and *1501 alleles. It is suggested that DPB1*6901 may have arisen as a result of a recombination event occurring between codons 58 and 64 between DPB1*0601 and DPB1*1101. The sequence of DPB1*1501 from codon 64 is not consistent with DPB1*6901. A linkage disequilibria analysis that examined 212 potential bone marrow recipients in which HLA-A to DQ haplotypes had been established by family studies showed linkage disequilibrium between HLA-B, DRB1 and DPB1 in some haplotypes and not others.

Polymorphism in the regulatory regions of the HLA-DPB1 gene.

Nucleotide sequencing of approximately 400 basepairs upstream from exon 1 of the DPB1 gene and sequence specific oligonucleotide hybridisation identified eight nucleotide positions to be polymorphic which were in linkage disequilibrium (LD) with DPA1 and DPB1 alleles. Substitutions at two sites (-230 and -224) formed three genotypes (DP-PRO1-3). DP-PRO 1 was the most common genotype and was in LD with DPA1*0103, *0202 and DPB1*0401, *0501. DP-PRO 2 was observed in LD with DPB1*02012, *1601, *1701 and DPA1*0104. DP-PRO3 was in LD with DPB1*0901, *1001 and DPA1*0201. Electrophoretic Mobility Shift Assays (EMSA) performed with restriction enzyme fragments showed substitutions at -230 and -224 not to be involved in binding nuclear proteins. Six substitutions were found on a single genotype (DP-PRO4) which was observed in seven samples; 67% of DP-PRO4 inferred haplotypes were HLA-A2-B46, DRB1*0901, DQB1*03032, DPA1*0401, DPB1*1301. Three of the substitutions occurred in conserved regulatory region boxes, W', X and Y, and three in the signal and leader sequences. EMSA competitive binding assays performed with oligonucleotide probes for the substitutions showed no difference in binding affinity for W' and X probes. The DP-PRO4 Y box had a decreased nuclear protein binding affinity compared to DP-PRO1-3. Whether the sum of the differences in DP-PRO4 relate to a change in the cell surface expression of HLA-DP is yet to be determined.

HLA-DPB1 alleles in house dust mite allergic patients.

Associations were sought between the presence of allergic sensitivity to Der p 1, a major allergen from the house dust mite, and the HLA-DPB1 genotype. Whilst allergic patients did not differ from controls in DPB1 allelic distribution, there was a correlation of DPB1*11011 with strongly reactive T-cell proliferative responses to Der p 1 and high titre specific IgE to Dermatophagoides pteronyssinus.

Matching for HLA DPA1 and DPB1 alleles in unrelated bone marrow transplantation.

The impact of donor-recipient DPA1 and DPB1 matching was examined in 122 unrelated bone marrow transplant pairs. All pairs were serologically matched at the time of transplantation for HLA class I and II and a majority also DRB1 allele matched. Retrospective A, B, C, DRB1, DQA1, DQB1 in addition to DPA1 and DPB1 allele matching was performed by molecular techniques. The percentage of pairs that were allele matched was as follows; HLA-A = 91% (n = 80), HLA-B = 94% (n = 80), HLA-C = 78% (n = 80), HLA-DRB1 = 96% (n = 122), HLA-DQA1 = 99% (n = 80), HLA-DQB1 = 92% (n = 122). 92 recipient/donor pairs with informative clinical data were available for analysis. DPA1 identity (no incompatibility in either direction) was observed in 57% and DPA1 compatibility in 76% of pairs with no apparent beneficial effect of matching on patient survival or Graft Versus Host Disease (GVHD). DPB1 identity was observed in 11% and compatibility in 27% of pairs. A significant improvement in patient survival was observed in DPB1 matched compared to one DPB1 mismatch (p < 0.01) and combined one and two DPB1 mismatched transplants (p = 0.03). This beneficial effect remained when allele mismatches at HLA-A, B, C, DRB1, DQA1, DQB1 were excluded (p = 0.05, p = 0.03, respectively). There was a significant association of increased frequency of severe GVHD (grades III-IV) compared to mild GVHD (grades I-II) with DPB1 mismatched transplants compared to DPB1 matched transplants (p = 0.04). In DPB1 mismatched transplants an association between patient survival and matching for individual DPB1 polymorphic regions was not observed; however in the HLA-A, B, DRB1, DQA1, DQB1 allele matched transplants a non significant increase in the frequency of Grade IV GVHD was observed in recipients who were negative compared to those who were positive for DPB1 alleles coding for glutamic acid at position 69.

A novel strategy for unambiguous DPB1 typing.

A DPB1 typing method is described that assigns DPB1 alleles into six groups based on polymorphism at amino acid positions 8-9 and 84-87 using sequence-specific priming (SSP). The results obtained allow the selection of primers for a subsequent sequence-specific oligonucleotide (SSO) hybridization procedure which permits DPB1 alleles to be analysed separately in a heterozygote individual. This has greatly reduced the occurrence of typing reaction patterns consistent with multiple combinations of DPB1 alleles seen in other DPB1 typing methods.

Protein structure-based design, synthesis, and biological evaluation of 5-thia-2,6-diamino-4(3H)-oxopyrimidines: potent inhibitors of glycinamide ribonucleotide transformylase with potent cell growth inhibition.

The design, synthesis, biochemical, and biological evaluation of a novel series of 5-thia-2,6-diamino-4(3H)-oxopyrimidine inhibitors of glycinamide ribonucleotide transformylase (GART) are described. The compounds were designed using the X-ray crystal structure of human GART. The monocyclic 5-thiapyrimidinones were synthesized by coupling an alkyl thiol with 5-bromo-2, 6-diamino-4(3H)-pyrimidinone, 20. The bicyclic compounds were prepared in both racemic and diastereomerically pure forms using two distinct synthetic routes. The compounds were found to have human GART KiS ranging from 30 microM to 2 nM. The compounds inhibited the growth of both L1210 and CCRF-CEM cells in culture with potencies down to the low nanomolar range and were found to be selective for the de novo purine biosynthesis pathway. The most potent inhibitors had 2,5-disubstituted thiophene rings attached to the glutamate moiety. Placement of a methyl substituent at the 4-position of the thiophene ring to give compounds 10, 18, and 19 resulted in inhibitors with significantly decreased mFBP affinity.

Structure-based design of substituted diphenyl sulfones and sulfoxides as lipophilic inhibitors of thymidylate synthase.

Six new diphenyl sulfoxide and five new diphenyl sulfones were designed, synthesized, and tested for their inhibition of human and Escherichia coli thymidylate synthase (TS) and of the growth of cells in tissue culture. The best sulfoxide inhibitor of human TS was 3-chloro-N-((3,4-dihydro-2-methyl-4-oxo-6-quinazolinyl)methyl)-4- (phenylsulfinyl)-N-(prop-2-ynyl)-aniline (7c) that had a Ki of 27 nM. No sulfone improved on TS inhibition by the previously reported 4-(N-((3,4-dihydro-2-methyl-6-quinazolinyl)methyl)-N-prop-2- ynylamino)phenyl phenyl sulfone (Ki = 12 nM). Nevertheless, one sulfone, 4-((2-chlorophenyl)sulfonyl)-N-((3,4-dihydro-2-methyl-4-oxo-6- quinazolinyl)methyl)-N-(prop-2-ynyl)aniline, was selected, on the basis of its inhibition of both TS and cell growth, for antitumor testing; it gave a 61% increase in life span to mice bearing the thymidino kinase-deficient L5178Y (TK-) lymphoma. A crystal structure of N-((3,4-dihydro-2-methyl-4-oxo-6-quinazolinyl)methyl)-4-((2- methylphenyl)sulfinyl)-N-(prop-2-ynyl)aniline complexed with E. coli TS was solved and revealed selective binding of one sulfoxide enantiomer. AMBER calculations showed that the enantioselection was due to asymmetric electrostatic effects at the mouth of the active site. In contrast, a similar crystal structure of the sulfoxide 7c, along with AMBER calculations, indicated that both enantiomers bound, but with different affinities. The side chain of Phe176 shifted in order to structurally accommodate the chlorine of the more weakly bound enantiomer.

An approach to rapid estimation of relative binding affinities of enzyme inhibitors: application to peptidomimetic inhibitors of the human immunodeficiency virus type 1 protease.

This report describes a method for rapid assessment of the binding affinities of a series of analogous ligands to an enzyme. This approach is based on two variables (scores), representing (i) the enthalpy of binding and (ii) the strength of hydrophobic interaction. The method is then used to evaluate the binding of 11 different peptidomimetic inhibitors to the HIV-1 protease. Three-dimensional structures of these enzyme-inhibitor complexes are modeled based on the crystal structures of HIV-1 protease complexes with the known inhibitors. These structures are minimized using the AMBER force field, and the scores of binding enthalpy for each of the ligands are calculated. A second score to represent the hydrophobic interaction between a pair of atoms uses an exponential function of distance between the atoms and the product of their atomic hydrophobicity constants. This exponential function is used to assess the hydrophobic interaction energy between an enzyme and its inhibitor and also to compute and display a 'molecular hydrophobicity map' as a 3D visualization tool. These methods are then applied to obtain trends in relative binding affinities of pairs of analogous inhibitors. Calculated scores agree well with corresponding results from thermodynamic cycle perturbation (TCP) simulations as well as experimental binding data. Since the proposed calculations are computationally cheaper and faster than TCP calculations, it is suggested that these scores can form the basis for rapid, preliminary theoretical screening of proposed derivatives of an inhibitor prior to TCP analysis, synthesis, and testing.

Structure-based design of lipophilic quinazoline inhibitors of thymidylate synthase.

To develop novel lipophilic thymidylate synthase (TS) inhibitors, the X-ray structure of Escherichia coli TS in ternary complex with FdUMP and the inhibitor 10-propargyl-5,8-dideazafolic acid (CB3717) was used as a basis for structure-based design. A total of 31 novel lipophilic TS inhibitors, lacking a glutamate residue, were synthesized; 26 of them had in common a N-((3,4-dihydro-2-methyl-6-quinazolinyl)methyl)-N-prop-2-ynylaniline+ ++ structure in which the aniline was appropriately substituted with simple lipophilic substituents either in position 3 or 4, or in both. Compounds were tested for their inhibition of E. coli TS and human TS and also for their inhibition of the growth in tissue culture of a murine leukemia, a human leukemia, and a thymidine kinase-deficient human adenocarcinoma. The crystal structures of five inhibitors complexed with E. coli TS were determined. Five main conclusions are drawn from this study. (i) A 3-substituent such as CF(3), iodo, or ethynyl enhances binding by up to 1 order of magnitude and in the case of CF(3) was proven to fill a nearby pocket in the enzyme. (ii) A simple strongly electron-withdrawing substituent such as NO(2) or CF(3)SO(2) in the 4-position enhances binding by 2 orders of magnitude; it is hypothesized that the transannular dipole so induced interacts favorably with the protein. (iii) Attempts to combine the enhancements of i and ii in the same molecule were generally unsuccessful (iv) A 4-C(6)H(5)SO(2) substituent provided both electron withdrawal and a van der Waal's interaction of the phenyl group with a hydrophobic surface at the mouth of the active site. The inhibition (K(is) = 12 nM) of human TS by this compound, 7n, showed that C(6)H(5)SO(2) provided virtually as much binding affinity as the CO-glutamate which it had replaced. (v) The series of compounds were poorly water soluble, and also the potent TS inhibition shown by several of them did not translate into good cytotoxicity. Compounds with large cyclic groups linked to position 4 by an SO or SO(2) group did, however, have IC(50)'s in the range 1-5 microM. Of these, 4-(N-((3,4-dihydro-2-methyl-6-quinazolinyl)methyl)-N-prop-2-ynylamino )phenyl phenyl sulfone, 7n, had IC(50)'s of about 1 microM and was chosen for further elaboration.

HLA-DRB1*0401 is associated with susceptibility to insulin-dependent diabetes mellitus independently of the DQB1 locus.

The distribution of DRB1*04 alleles and DRB1/DQB1 haplotypes was analysed in 57 DR4+ caucasoid subjects with insulin-dependent diabetes mellitus (IDDM) and 96 DR4+ healthy controls selected on the basis of DR serology, and the findings were analysed in relation to age at diagnosis of IDDM. DNA samples were amplified using specific DR and DQ primers and hybridized with sequence-specific oligonucleotide probes. A significantly increased combined frequency of DRB1*0401 and 0402 was observed in IDDM subjects aged < or = 12 years at diagnosis (allele frequency 88.4% compared with 62.0% in controls, P < 0.025). There was a non-significant increase in DRB1*0401 and 0402 in IDDM subjects < or = 12 years when compared with IDDM subjects > 12 years (P < 0.1). DRB1*0404 was decreased in the total IDDM subject group compared with controls (4.8% vs. 19.0%, P < 0.025) but did not reach statistical significance in the individual age at diagnosis groups. In contrast, the frequency of DQB1*0302 was increased uniformly across both ages at diagnosis groups. In controls DRB1*0401 occurred in haplotype association with DQB1*0301 in a significantly greater frequency than with DQB1*0302. However, 95.0% of DRB1*0401 IDDM subjects were DQB1*0302. DRB1*0404, which was decreased in frequency in IDDM subjects, occurred in association significantly more frequently with DQB1*0302 in controls. These results imply that DRB1 and DQB1 have independent roles as HLA susceptibility genes in IDDM. DQB1 may have a permissive role whereas DRB1 could influence the rate at which underlying disease progresses to clinical IDDM.