Fine-mapping the MHC locus in juvenile idiopathic arthritis (JIA) reveals genetic heterogeneity corresponding to distinct adult inflammatory arthritic diseases.

OBJECTIVES: Juvenile idiopathic arthritis (JIA) is a heterogeneous group of diseases, comprising seven categories. Genetic data could potentially be used to help redefine JIA categories and improve the current classification system. The human leucocyte antigen (HLA) region is strongly associated with JIA. Fine-mapping of the region was performed to look for similarities and differences in HLA associations between the JIA categories and define correspondences with adult inflammatory arthritides. METHODS: Dense genotype data from the HLA region, from the Immunochip array for 5043 JIA cases and 14 390 controls, were used to impute single-nucleotide polymorphisms, HLA classical alleles and amino acids. Bivariate analysis was performed to investigate genetic correlation between the JIA categories. Conditional analysis was used to identify additional effects within the region. Comparison of the findings with those in adult inflammatory arthritic diseases was performed. RESULTS: We identified category-specific associations and have demonstrated for the first time that rheumatoid factor (RF)-negative polyarticular JIA and oligoarticular JIA are genetically similar in their HLA associations. We also observe that each JIA category potentially has an adult counterpart. The RF-positive polyarthritis association at HLA-DRB1 amino acid at position 13 mirrors the association in adult seropositive rheumatoid arthritis (RA). Interestingly, the combined oligoarthritis and RF-negative polyarthritis dataset shares the same association with adult seronegative RA. CONCLUSIONS: The findings suggest the value of using genetic data in helping to classify the categories of this heterogeneous disease. Mapping JIA categories to adult counterparts could enable shared knowledge of disease pathogenesis and aetiology and facilitate transition from paediatric to adult services.

Genetic variation in uncontrolled childhood asthma despite ICS treatment.

Genetic variation may partly explain asthma treatment response heterogeneity. We aimed to identify common and rare genetic variants associated with asthma that was not well controlled despite inhaled corticosteroid (ICS) treatment. Data of 110 children was collected in the Children Asthma Therapy Optimal trial. Associations of genetic variation with measures of lung function (FEV1%pred), airway hyperresponsiveness (AHR) to methacholine (Mch PD20) and treatment response outcomes were analyzed using the exome chip. The 17q12-21 locus (containing ORMDL3 and GSMDB) previously associated with childhood asthma was investigated separately. Single-nucleotide polymorphisms (SNPs) in the 17q12-21 locus were found nominally associated with the outcomes. The strongest association in this region was found for rs72821893 in KRT25 with FEV1%pred (P=3.75*10(-5)), Mch PD20 (P=0.00095) and Mch PD20-based treatment outcome (P=0.006). No novel single SNPs or burden tests were significantly associated with the outcomes. The 17q12-21 region was associated with FEV1%pred and AHR, and additionally with ICS treatment response.

Common variants associated with blood lipid levels do not affect carotid plaque composition.

INTRODUCTION: Although plasma lipid levels are known to influence the risk of cardiovascular disease (CVD), little is known about their effect on atherosclerotic plaque composition. To date, large-scale genome-wide association studies have identified 157 common single-nucleotide polymorphisms (SNPs) that influence plasma lipid levels, providing a powerful tool to investigate the effect of plasma lipid levels on atherosclerotic plaque composition. METHODS: In this study, we included 1443 carotid endarterectomy patients from the Athero-Express Biobank Study with genotype data. Plasma concentrations of high-density lipoprotein cholesterol (HDL-C), low-density lipoprotein cholesterol (LDL-C), total cholesterol (TC) and triglycerides (TG) were determined at the time of endarterectomy. Atherosclerotic plaques, obtained during surgery, were histologically examined. For all patients, we calculated weighted genetic burden scores (GBS) for all lipid traits on the basis of the available genotype data. Plasma lipid levels and GBS were tested for association with 7 histological features using linear and logistic regression models. RESULTS: All GBS were associated with their respective plasma lipid concentrations (pHDL-C = 2.4 × 10(-14), pLDL-C = 0.003, pTC = 2.1 × 10(-6), pTG = 3.4 × 10(-8)). Neither the measured plasma lipids, nor the GBS, were associated with histological features of atherosclerotic plaque composition. In addition, neither the plasma lipids nor the GBS were associated with clinical endpoints within 3 years of follow-up, with the notable exception of a negative association between HDL-C and composite cardiovascular endpoints. CONCLUSION: This study found no evidence that plasma lipid levels or their genetic determinants influence carotid plaque composition.

Genome-wide association study of lymphoblast cell viability after clozapine exposure.

Clozapine is an antipsychotic drug with proven efficacy in treatment-resistant schizophrenia but also known to induce potentially lethal agranulocytosis (CIA) in 1% of patients. Genetic factors are likely to play a role in the molecular basis of CIA. We explored an in vitro system to study the genetic susceptibility of CIA. Cell viability was measured in 90 lymphoblast cell lines exposed to a series of increasing concentrations of clozapine for 48 hr. Quantitative trait measures of cell viability as well as area under the survival curve were used in a linear mixed model for genome-wide association analyses. The estimated heritability of clozapine-induced cell viability reduction in these cell lines is h2=0.76. No genome-wide significant association was observed after correction for multiple testing. Two independent loci with nominal evidence of association were observed at 30× clinical clozapine concentration: rs2709505 (P=1.41×10(-8)) in an intron of MDFIC and rs10457252 (P=1.79×10(-8)) located in a gene desert at chromosome 6q21. We identified one locus (rs1293970) near PRG4 that was consistently associated for all separate concentration analyses at P<5×10(-5). PRG4 encodes hemangiopoietin, a growth stimulator for hematopoietic stem cells. No evidence was observed for involvement of the MHC region. Our results demonstrate that clozapine-induced viability reduction in lymphoblast cell lines is a heritable, polygenic trait. Thus, in vitro models of CIA might be a useful tool for future discovery of genetic risk factors, although larger sample sizes will be required to unambiguously identify these loci.

Cholesteryl ester transfer protein polymorphisms, statin use, and their impact on cholesterol levels and cardiovascular events.

The association of nonfunctional variants of the cholesteryl ester transfer protein (CETP) with efficacy of statins has been a subject of debate. We evaluated whether three functional CETP variants influence statin efficacy. The effect of CETP genotype on achieved levels of high-density lipoprotein cholesterol (HDLc), low-density lipoprotein cholesterol (LDLc), and total cholesterol during statin treatment was estimated by meta-analysis of the linear regression outcomes of three studies (11,021 individuals). The effect of these single-nucleotide polymorphisms (SNPs) on statin response in protecting against myocardial infarction (MI) was estimated by meta-analysis of statin × SNP interaction terms from logistic regression in five studies (16,570 individuals). The enhancer SNP rs3764261 significantly increased HDLc by 0.02 mmol/l per T allele (P = 6 × 10(-5)) and reduced protection against MI by statins (interaction odds ratio (OR) = 1.19 per T allele; P = 0.04). Focusing on functional CETP variants, we showed that in carriers of the rs3764261 T variant, HDLc increased more during statin treatment, and protection against MI by statins appeared to be reduced as compared with those in noncarriers.

Classical HLA-DRB1 and DPB1 alleles account for HLA associations with primary biliary cirrhosis.

Susceptibility to primary biliary cirrhosis (PBC) is strongly associated with human leukocyte antigen (HLA)-region polymorphisms. To determine if associations can be explained by classical HLA determinants, we studied Italian, 676 cases and 1440 controls, genotyped with dense single-nucleotide polymorphisms (SNPs) for which classical HLA alleles and amino acids were imputed. Although previous genome-wide association studies and our results show stronger SNP associations near DQB1, we demonstrate that the HLA signals can be attributed to classical DRB1 and DPB1 genes. Strong support for the predominant role of DRB1 is provided by our conditional analyses. We also demonstrate an independent association of DPB1. Specific HLA-DRB1 genes (*08, *11 and *14) account for most of the DRB1 association signal. Consistent with previous studies, DRB1*08 (P=1.59 × 10(-11)) was the strongest predisposing allele, whereas DRB1*11 (P=1.42 × 10(-10)) was protective. Additionally, DRB1*14 and the DPB1 association (DPB1*03:01; P=9.18 × 10(-7)) were predisposing risk alleles. No signal was observed in the HLA class 1 or class 3 regions. These findings better define the association of PBC with HLA and specifically support the role of classical HLA-DRB1 and DPB1 genes and alleles in susceptibility to PBC.

Amino acid position 11 of HLA-DRß1 is a major determinant of chromosome 6p association with ulcerative colitis.

The major histocompatibility complex (MHC) on chromosome 6p is an established risk locus for ulcerative colitis (UC) and Crohn's disease (CD). We aimed to better define MHC association signals in UC and CD by combining data from dense single-nucleotide polymorphism (SNP) genotyping and from imputation of classical human leukocyte antigen (HLA) types, their constituent SNPs and corresponding amino acids in 562 UC, 611 CD and 1428 control subjects. Univariate and multivariate association analyses were performed, controlling for ancestry. In univariate analyses, absence of the rs9269955 C allele was strongly associated with risk for UC (P = 2.67 × 10(-13)). rs9269955 is a SNP in the codon for amino acid position 11 of HLA-DRß1, located in the P6 pocket of the HLA-DR antigen binding cleft. This amino acid position was also the most significantly UC-associated amino acid in omnibus tests (P = 2.68 × 10(-13)). Multivariate modeling identified rs9269955-C and 13 other variants in best predicting UC vs control status. In contrast, there was only suggestive association evidence between the MHC and CD. Taken together, these data demonstrate that variation at HLA-DRß1, amino acid 11 in the P6 pocket of the HLA-DR complex antigen binding cleft is a major determinant of chromosome 6p association with UC.

Overview of the Rapid Response data.

The Type I Diabetes Genetics Consortium (T1DGC) Rapid Response Workshop was established to evaluate published candidate gene associations in a large collection of affected sib-pair (ASP) families. We report on our quality control (QC) and preliminary family-based association analyses. A random sample of blind duplicates was analyzed for QC. Quality checks, including examination of plate-panel yield, marker yield, Hardy-Weinberg equilibrium, mismatch error rate, Mendelian error rate, and allele distribution across plates, were performed. Genotypes from 2324 families within nine cohorts were obtained from a panel of 21 candidate genes, including 384 single-nucleotide polymorphisms on two genotyping platforms performed at the Broad Institute Center for Genotyping and Analysis (Cambridge, MA, USA). The T1DGC Rapid Response project, following rigorous QC procedures, resulted in a 2297 family, 9688 genotyped individual database on a single-candidate gene panel. The available data include 9005 individuals with genotype data from both platforms and 683 individuals genotyped (276 in Illumina; 407 in Sequenom) on only one platform.

Whole-genome association study of bipolar disorder.

We performed a genome-wide association scan in 1461 patients with bipolar (BP) 1 disorder, 2008 controls drawn from the Systematic Treatment Enhancement Program for Bipolar Disorder and the University College London sample collections with successful genotyping for 372,193 single nucleotide polymorphisms (SNPs). Our strongest single SNP results are found in myosin5B (MYO5B; P=1.66 x 10(-7)) and tetraspanin-8 (TSPAN8; P=6.11 x 10(-7)). Haplotype analysis further supported single SNP results highlighting MYO5B, TSPAN8 and the epidermal growth factor receptor (MYO5B; P=2.04 x 10(-8), TSPAN8; P=7.57 x 10(-7) and EGFR; P=8.36 x 10(-8)). For replication, we genotyped 304 SNPs in family-based NIMH samples (n=409 trios) and University of Edinburgh case-control samples (n=365 cases, 351 controls) that did not provide independent replication after correction for multiple testing. A comparison of our strongest associations with the genome-wide scan of 1868 patients with BP disorder and 2938 controls who completed the scan as part of the Wellcome Trust Case-Control Consortium indicates concordant signals for SNPs within the voltage-dependent calcium channel, L-type, alpha 1C subunit (CACNA1C) gene. Given the heritability of BP disorder, the lack of agreement between studies emphasizes that susceptibility alleles are likely to be modest in effect size and require even larger samples for detection.

HOMSTRAD: adding sequence information to structure-based alignments of homologous protein families.

summary: We describe an extension to the Homologous Structure Alignment Database (HOMSTRAD; Mizuguchi et al., Protein Sci., 7, 2469-2471, 1998a) to include homologous sequences derived from the protein families database Pfam (Bateman et al., Nucleic Acids Res., 28, 263-266, 2000). HOMSTRAD is integrated with the server FUGUE (Shi et al., submitted, 2001) for recognition and alignment of homologues, benefitting from the combination of abundant sequence information and accurate structure-based alignments. AVAILABILITY The HOMSTRAD database is available at: http://www-cryst.bioc.cam.ac.uk/homstrad/. Query sequences can be submitted to the homology recognition/alignment server FUGUE at: http://www-cryst.bioc.cam.ac.uk/fugue/.

Sequence-structure homology recognition by iterative alignment refinement and comparative modeling.

Our approach to fold recognition for the fourth critical assessment of techniques for protein structure prediction (CASP4) experiment involved the use of the FUGUE sequence-structure homology recognition program (http://www-cryst.bioc.cam.ac.uk/fugue), followed by model building. We treat models as hypotheses and examine these to determine whether they explain the available data. Our method depends heavily on environment-specific substitution tables derived from our database of structural alignments of homologous proteins (HOMSTRAD, http://www-cryst.bioc.cam.ac.uk/homstrad/). FUGUE uses these tables to incorporate structural information into profiles created from HOMSTRAD alignments that are matched against a profile created for the target from multiple sequence alignment. In addition, environment-specific substitution tables are used throughout the modeling procedure and as part of the model evaluation. Annotation of sequence alignments with JOY, to reflect local structural features, proved valuable, both for modifying hypotheses, and for rejecting predictions when the expected pattern of conservation is not observed. Our stringency in rejecting incorrect predictions led us to submit a relatively small number of models, including only a low number of false positives, resulting in a high average score.

MODBASE, a database of annotated comparative protein structure models.

MODBASE is a queryable database of annotated comparative protein structure models. The models are derived by MODPIPE, an automated modeling pipeline relying on the programs PSI-BLAST and MODELLER. The database currently contains 3D models for substantial portions of approximately 17 000 proteins from 10 complete genomes, including those of Caenorhabditis elegans, Saccharomyces cerevisiae and Escherichia coli, as well as all the available sequences from Arabidopsis thaliana and Homo sapiens. The database also includes fold assignments and alignments on which the models were based. In addition, special care is taken to assess the quality of the models. ModBase is accessible through a web interface at http://guitar.rockefeller.edu/modbase/

Molecular dynamics simulations of the hyperthermophilic protein sac7d from Sulfolobus acidocaldarius: contribution of salt bridges to thermostability.

Hyperthermophilic proteins often possess an increased number of surface salt bridges compared with their mesophilic homologues. However, salt bridges are generally thought to be of minor importance in protein stability at room temperature. In an effort to understand why this may no longer be true at elevated temperatures, we performed molecular dynamics simulations of the hyperthermophilic protein Sac7d at 300 K, 360 K, and 550 K. The three trajectories are stable on the nanosecond timescale, as evidenced by the analysis of several time-resolved properties. The simulations at 300 K and (to a lesser extent) 360 K are also compatible with nuclear Overhauser effect-derived distances. Raising the temperature from 300 K to 360 K results in a less favourable protein-solvent interaction energy, and a more favourable intraprotein interaction energy. Both effects are almost exclusively electrostatic in nature and dominated by contributions due to charged side-chains. The reduced solvation is due to a loss of spatial and orientational structure of water around charged side-chains, which is a consequence of the increased thermal motion in the solvent. The favourable change in the intraprotein Coulombic interaction energy is essentially due to the tightening of salt bridges. Assuming that charged side-chains are on average more distant from one another in the unfolded state than in the folded state, it follows that salt bridges may contribute to protein stability at elevated temperatures because (i) the solvation free energy of charged side-chains is more adversely affected in the unfolded state than in the folded state by an increase in temperature, and (ii) due to the tightening of salt bridges, unfolding implies a larger unfavourable increase in the intraprotein Coulombic energy at higher temperature. Possible causes for the unexpected stability of the protein at 550 K are also discussed.