Unraveling multiple MHC gene associations with systemic lupus erythematosus: model choice indicates a role for HLA alleles and non-HLA genes in Europeans.

We have performed a meta-analysis of the major-histocompatibility-complex (MHC) region in systemic lupus erythematosus (SLE) to determine the association with both SNPs and classical human-leukocyte-antigen (HLA) alleles. More specifically, we combined results from six studies and well-known out-of-study control data sets, providing us with 3,701 independent SLE cases and 12,110 independent controls of European ancestry. This study used genotypes for 7,199 SNPs within the MHC region and for classical HLA alleles (typed and imputed). Our results from conditional analysis and model choice with the use of the Bayesian information criterion show that the best model for SLE association includes both classical loci (HLA-DRB1(∗)03:01, HLA-DRB1(∗)08:01, and HLA-DQA1(∗)01:02) and two SNPs, rs8192591 (in class III and upstream of NOTCH4) and rs2246618 (MICB in class I). Our approach was to perform a stepwise search from multiple baseline models deduced from a priori evidence on HLA-DRB1 lupus-associated alleles, a stepwise regression on SNPs alone, and a stepwise regression on HLA alleles. With this approach, we were able to identify a model that was an overwhelmingly better fit to the data than one identified by simple stepwise regression either on SNPs alone (Bayes factor [BF] > 50) or on classical HLA alleles alone (BF > 1,000).

Association of systemic lupus erythematosus with C8orf13-BLK and ITGAM-ITGAX.

BACKGROUND: Systemic lupus erythematosus (SLE) is a clinically heterogeneous disease in which the risk of disease is influenced by complex genetic and environmental contributions. Alleles of HLA-DRB1, IRF5, and STAT4 are established susceptibility genes; there is strong evidence for the existence of additional risk loci. METHODS: We genotyped more than 500,000 single-nucleotide polymorphisms (SNPs) in DNA samples from 1311 case subjects with SLE and 1783 control subjects; all subjects were North Americans of European descent. Genotypes from 1557 additional control subjects were obtained from public data repositories. We measured the association between the SNPs and SLE after applying strict quality-control filters to reduce technical artifacts and to correct for the presence of population stratification. Replication of the top loci was performed in 793 case subjects and 857 control subjects from Sweden. RESULTS: Genetic variation in the region upstream from the transcription initiation site of the gene encoding B lymphoid tyrosine kinase (BLK) and C8orf13 (chromosome 8p23.1) was associated with disease risk in both the U.S. and Swedish case-control series (rs13277113; odds ratio, 1.39; P=1x10(-10)) and also with altered levels of messenger RNA in B-cell lines. In addition, variants on chromosome 16p11.22, near the genes encoding integrin alpha M (ITGAM, or CD11b) and integrin alpha X (ITGAX), were associated with SLE in the combined sample (rs11574637; odds ratio, 1.33; P=3x10(-11)). CONCLUSIONS: We identified and then confirmed through replication two new genetic loci for SLE: a promoter-region allele associated with reduced expression of BLK and increased expression of C8orf13 and variants in the ITGAM-ITGAX region.

Specificity of the STAT4 genetic association for severe disease manifestations of systemic lupus erythematosus.

Systemic lupus erythematosus (SLE) is a genetically complex disease with heterogeneous clinical manifestations. A polymorphism in the STAT4 gene has recently been established as a risk factor for SLE, but the relationship with specific SLE subphenotypes has not been studied. We studied 137 SNPs in the STAT4 region genotyped in 4 independent SLE case series (total n = 1398) and 2560 healthy controls, along with clinical data for the cases. Using conditional testing, we confirmed the most significant STAT4 haplotype for SLE risk. We then studied a SNP marking this haplotype for association with specific SLE subphenotypes, including autoantibody production, nephritis, arthritis, mucocutaneous manifestations, and age at diagnosis. To prevent possible type-I errors from population stratification, we reanalyzed the data using a subset of subjects determined to be most homogeneous based on principal components analysis of genome-wide data. We confirmed that four SNPs in very high LD (r(2) = 0.94 to 0.99) were most strongly associated with SLE, and there was no compelling evidence for additional SLE risk loci in the STAT4 region. SNP rs7574865 marking this haplotype had a minor allele frequency (MAF) = 31.1% in SLE cases compared with 22.5% in controls (OR = 1.56, p = 10(-16)). This SNP was more strongly associated with SLE characterized by double-stranded DNA autoantibodies (MAF = 35.1%, OR = 1.86, p<10(-19)), nephritis (MAF = 34.3%, OR = 1.80, p<10(-11)), and age at diagnosis<30 years (MAF = 33.8%, OR = 1.77, p<10(-13)). An association with severe nephritis was even more striking (MAF = 39.2%, OR = 2.35, p<10(-4) in the homogeneous subset of subjects). In contrast, STAT4 was less strongly associated with oral ulcers, a manifestation associated with milder disease. We conclude that this common polymorphism of STAT4 contributes to the phenotypic heterogeneity of SLE, predisposing specifically to more severe disease.

STAT4 and the risk of rheumatoid arthritis and systemic lupus erythematosus.

BACKGROUND: Rheumatoid arthritis is a chronic inflammatory disease with a substantial genetic component. Susceptibility to disease has been linked with a region on chromosome 2q. METHODS: We tested single-nucleotide polymorphisms (SNPs) in and around 13 candidate genes within the previously linked chromosome 2q region for association with rheumatoid arthritis. We then performed fine mapping of the STAT1-STAT4 region in a total of 1620 case patients with established rheumatoid arthritis and 2635 controls, all from North America. Implicated SNPs were further tested in an independent case-control series of 1529 patients with early rheumatoid arthritis and 881 controls, all from Sweden, and in a total of 1039 case patients and 1248 controls from three series of patients with systemic lupus erythematosus. RESULTS: A SNP haplotype in the third intron of STAT4 was associated with susceptibility to both rheumatoid arthritis and systemic lupus erythematosus. The minor alleles of the haplotype-defining SNPs were present in 27% of chromosomes of patients with established rheumatoid arthritis, as compared with 22% of those of controls (for the SNP rs7574865, P=2.81x10(-7); odds ratio for having the risk allele in chromosomes of patients vs. those of controls, 1.32). The association was replicated in Swedish patients with recent-onset rheumatoid arthritis (P=0.02) and matched controls. The haplotype marked by rs7574865 was strongly associated with lupus, being present on 31% of chromosomes of case patients and 22% of those of controls (P=1.87x10(-9); odds ratio for having the risk allele in chromosomes of patients vs. those of controls, 1.55). Homozygosity of the risk allele, as compared with absence of the allele, was associated with a more than doubled risk for lupus and a 60% increased risk for rheumatoid arthritis. CONCLUSIONS: A haplotype of STAT4 is associated with increased risk for both rheumatoid arthritis and systemic lupus erythematosus, suggesting a shared pathway for these illnesses.

Full-sequence computational design and solution structure of a thermostable protein variant.

Computational protein design procedures were applied to the redesign of the entire sequence of a 51 amino acid residue protein, Drosophila melanogaster engrailed homeodomain. Various sequence optimization algorithms were compared and two resulting designed sequences were experimentally evaluated. The two sequences differ by 11 mutations and share 22% and 24% sequence identity with the wild-type protein. Both computationally designed proteins were considerably more stable than the naturally occurring protein, with midpoints of thermal denaturation greater than 99 degrees C. The solution structure was determined for one of the two sequences using multidimensional heteronuclear NMR spectroscopy, and the structure was found to closely match the original design template scaffold.

Dioxane contributes to the altered conformation and oligomerization state of a designed engrailed homeodomain variant.

Our goal was to compute a stable, full-sequence design of the Drosophila melanogaster engrailed homeodomain. Thermal and chemical denaturation data indicated the design was significantly more stable than was the wild-type protein. The data were also nearly identical to those for a similar, later full-sequence design, which was shown by NMR to adopt the homeodomain fold: a three-helix, globular monomer. However, a 1.65 A crystal structure of the design described here turned out to be of a completely different fold: a four-helix, rodlike tetramer. The crystallization conditions included approximately 25% dioxane, and subsequent experiments by circular dichroism and sedimentation velocity analytical ultracentrifugation indicated that dioxane increases the helicity and oligomerization state of the designed protein. We attribute at least part of the discrepancy between the target fold and the crystal structure to the presence of a high concentration of dioxane.

Replication of the association between the C8orf13-BLK region and systemic lupus erythematosus in a Japanese population.

OBJECTIVE: Recent genome-wide association studies identified an association between single-nucleotide polymorphisms (SNPs) in the C8orf13 region of BLK, the B lymphoid tyrosine kinase gene, with systemic lupus erythematosus (SLE) in Caucasians. The purpose of this study was to evaluate the significance of this region in the genetic background of Japanese patients with SLE. METHODS: Fourteen tag SNPs in the C8orf13-BLK region were genotyped in 327 Japanese patients with SLE and 322 healthy Japanese controls. The population-attributable risk percentage (PAR%) of rs13277113 in Japanese was compared with that in Caucasians as well as with that of other SLE susceptibility genes in Japanese. RESULTS: As in Caucasians, rs13277113A demonstrated the strongest association in Japanese (P = 1.73 x 10(-6) for the genotype frequency, P = 4.75 x 10(-7) for the allele frequency, odds ratio [OR] 2.44 [95% confidence interval (95% CI) 1.43-4.16]). The association in Japanese was consistent with a recessive model (P = 2.74 x 10(-7), OR 2.27 [95% CI 1.66-3.11]). In contrast to the Caucasian population, this risk allele was the major allele in the Japanese population. Because both the risk allele frequency and the OR were higher in Japanese than in Caucasians, the PAR% of rs13277113 was estimated to be much higher in Japanese (35.4%) than in Caucasians (16.2%), and the second highest among the 6 confirmed SLE susceptibility genes in Japanese. CONCLUSION: The association of the C8orf13-BLK region with SLE was replicated in a Japanese population. Contribution of this region to the genetic predisposition to SLE appeared to be greater in Japanese than in Caucasians.

European population substructure is associated with mucocutaneous manifestations and autoantibody production in systemic lupus erythematosus.

OBJECTIVE: To determine whether genetic substructure in European-derived populations is associated with specific manifestations of systemic lupus erythematosus (SLE), including mucocutaneous phenotypes, autoantibody production, and renal disease. METHODS: SLE patients of European descent (n=1,754) from 8 case collections were genotyped for >1,400 ancestry informative markers that define a north-south gradient of European substructure. Using the Structure program, each SLE patient was characterized in terms of percent Northern (versus percent Southern) European ancestry based on these genetic markers. Nonparametric methods, including tests for trend, were used to identify associations between Northern European ancestry and specific SLE manifestations. RESULTS: In multivariate analyses, increasing levels of Northern European ancestry were significantly associated with photosensitivity (Ptrend=0.0021, odds ratio for highest quartile of Northern European ancestry versus lowest quartile [ORhigh-low] 1.64, 95% confidence interval [95% CI] 1.13-2.35) and discoid rash (Ptrend=0.014, ORhigh-low 1.93, 95% CI 0.98-3.83). In contrast, increasing levels of Northern European ancestry had a protective effect against the production of anticardiolipin autoantibodies (Ptrend=1.6x10(-4), ORhigh-low 0.46, 95% CI 0.30-0.69) and anti-double-stranded DNA autoantibodies (Ptrend=0.017, ORhigh-low 0.67, 95% CI 0.46-0.96). CONCLUSION: This study demonstrates that specific SLE manifestations vary according to Northern versus Southern European ancestry. Thus, genetic ancestry may contribute to the clinical heterogeneity and variation in disease outcomes among SLE patients of European descent. Moreover, these results suggest that genetic studies of SLE subphenotypes will need to carefully address issues of population substructure based on genetic ancestry.

Role of STAT4 polymorphisms in systemic lupus erythematosus in a Japanese population: a case-control association study of the STAT1-STAT4 region.

INTRODUCTION: Recent studies identified STAT4 (signal transducers and activators of transcription-4) as a susceptibility gene for systemic lupus erythematosus (SLE). STAT1 is encoded adjacently to STAT4 on 2q32.2-q32.3, upregulated in peripheral blood mononuclear cells from SLE patients, and functionally relevant to SLE. This study was conducted to test whether STAT4 is associated with SLE in a Japanese population also, to identify the risk haplotype, and to examine the potential genetic contribution of STAT1. To accomplish these aims, we carried out a comprehensive association analysis of 52 tag single nucleotide polymorphisms (SNPs) encompassing the STAT1-STAT4 region. METHODS: In the first screening, 52 tag SNPs were selected based on HapMap Phase II JPT (Japanese in Tokyo, Japan) data, and case-control association analysis was carried out on 105 Japanese female patients with SLE and 102 female controls. For associated SNPs, additional cases and controls were genotyped and association was analyzed using 308 SLE patients and 306 controls. Estimation of haplotype frequencies and an association study using the permutation test were performed with Haploview version 4.0 software. Population attributable risk percentage was estimated to compare the epidemiological significance of the risk genotype among populations. RESULTS: In the first screening, rs7574865, rs11889341, and rs10168266 in STAT4 were most significantly associated (P < 0.01). Significant association was not observed for STAT1. Subsequent association studies of the three SNPs using 308 SLE patients and 306 controls confirmed a strong association of the rs7574865T allele (SLE patients: 46.3%, controls: 33.5%, P = 4.9 x 10(-6), odds ratio 1.71) as well as TTT haplotype (rs10168266/rs11889341/rs7574865) (P = 1.5 x 10(-6)). The association was stronger in subgroups of SLE with nephritis and anti-double-stranded DNA antibodies. Population attributable risk percentage was estimated to be higher in the Japanese population (40.2%) than in Americans of European descent (19.5%). CONCLUSIONS: The same STAT4 risk allele is associated with SLE in Caucasian and Japanese populations. Evidence for a role of STAT1 in genetic susceptibility to SLE was not detected. The contribution of STAT4 for the genetic background of SLE may be greater in the Japanese population than in Americans of European descent.

Three functional variants of IFN regulatory factor 5 (IRF5) define risk and protective haplotypes for human lupus.

Systematic genome-wide studies to map genomic regions associated with human diseases are becoming more practical. Increasingly, efforts will be focused on the identification of the specific functional variants responsible for the disease. The challenges of identifying causal variants include the need for complete ascertainment of genetic variants and the need to consider the possibility of multiple causal alleles. We recently reported that risk of systemic lupus erythematosus (SLE) is strongly associated with a common SNP in IFN regulatory factor 5 (IRF5), and that this variant altered spicing in a way that might provide a functional explanation for the reproducible association to SLE risk. Here, by resequencing and genotyping in patients with SLE, we find evidence for three functional alleles of IRF5: the previously described exon 1B splice site variant, a 30-bp in-frame insertion/deletion variant of exon 6 that alters a proline-, glutamic acid-, serine- and threonine-rich domain region, and a variant in a conserved polyA+ signal sequence that alters the length of the 3' UTR and stability of IRF5 mRNAs. Haplotypes of these three variants define at least three distinct levels of risk to SLE. Understanding how combinations of variants influence IRF5 function may offer etiological and therapeutic insights in SLE; more generally, IRF5 and SLE illustrates how multiple common variants of the same gene can together influence risk of common disease.

A search algorithm for fixed-composition protein design.

We present a computational protein design algorithm for finding low-energy sequences of fixed amino acid composition. The search algorithms used in protein design typically do not restrict amino acid composition. However, the random energy model of Shakhnovich suggests that the use of fixed-composition sequences may circumvent defects in the modeling of the denatured state. Our algorithm, FC_FASTER, links fixed-composition versions of Monte Carlo and the FASTER algorithm. As proof of principle, FC_FASTER was tested on an experimentally validated, full-sequence design of the beta1 domain of protein G. For the wild-type composition, FC_FASTER found a lower energy sequence than the experimentally validated sequence. Also, for a different composition, FC_FASTER found the hypothetical lowest-energy sequence in 14 out of 32 trials.

Preprocessing of rotamers for protein design calculations.

We have developed a process that significantly reduces the number of rotamers in computational protein design calculations. This process, which we call Vegas, results in dramatic computational performance increases when used with algorithms based on the dead-end elimination (DEE) theorem. Vegas estimates the energy of each rotamer at each position by fixing each rotamer in turn and utilizing various search algorithms to optimize the remaining positions. Algorithms used for this context specific optimization can include Monte Carlo, self-consistent mean field, and the evaluation of an expression that generates a lower bound energy for the fixed rotamer. Rotamers with energies above a user-defined cutoff value are eliminated. We found that using Vegas to preprocess rotamers significantly reduced the calculation time of subsequent DEE-based algorithms while retaining the global minimum energy conformation. For a full boundary design of a 51 amino acid fragment of engrailed homeodomain, the total calculation time was reduced by 12-fold.

Exact rotamer optimization for protein design.

Computational methods play a central role in the rational design of novel proteins. The present work describes a new hybrid exact rotamer optimization (HERO) method that builds on previous dead-end elimination algorithms to yield dramatic performance enhancements. Measured on experimentally validated physical models, these improvements make it possible to perform previously intractable designs of entire protein core, surface, or boundary regions. Computational demonstrations include a full core design of the variable domains of the light and heavy chains of catalytic antibody 48G7 FAB with 74 residues and 10(128) conformations, a full core/boundary design of the beta1 domain of protein G with 25 residues and 10(53) conformations, and a full surface design of the beta1 domain of protein G with 27 residues and 10(60) conformations. In addition, a full sequence design of the beta1 domain of protein G is used to demonstrate the strong dependence of algorithm performance on the exact form of the potential function and the fidelity of the rotamer library. These results emphasize that search algorithm performance for protein design can only be meaningfully evaluated on physical models that have been subjected to experimental scrutiny. The new algorithm greatly facilitates ongoing efforts to engineer increasingly complex protein features.