Identification of two independent risk factors for lupus within the MHC in United Kingdom families

The association of the major histocompatibility complex (MHC) with SLE is well established yet the causal variants arising from this region remain to be identified, largely due to inadequate study design and the strong linkage disequilibrium demonstrated by genes across this locus. The majority of studies thus far have identified strong association with classical class II alleles, in particular HLA-DRB1*0301 and HLA-DRB1*1501. Additional associations have been reported with class III alleles; specifically, complement C4 null alleles and a tumor necrosis factor promoter SNP (TNF-308G/A). However, the relative effects of these class II and class III variants have not been determined. We have thus used a family-based approach to map association signals across the MHC class II and class III regions in a cohort of 314 complete United Kingdom Caucasian SLE trios by typing tagging SNPs together with classical typing of the HLA-DRB1 locus. Using TDT and conditional regression analyses, we have demonstrated the presence of two distinct and independent association signals in SLE: HLA-DRB1*0301 (nominal p = 4.9 x 10(-8), permuted p < 0.0001, OR = 2.3) and the T allele of SNP rs419788 (nominal p = 4.3 x 10(-8), permuted p < 0.0001, OR = 2.0) in intron 6 of the class III region gene SKIV2L. Assessment of genotypic risk demonstrates a likely dominant model of inheritance for HLA-DRB1*0301, while rs419788-T confers susceptibility in an additive manner. Furthermore, by comparing transmitted and untransmitted parental chromosomes, we have delimited our class II signal to a 180 kb region encompassing the alleles HLA-DRB1*0301-HLA-DQA1*0501-HLA-DQB1*0201 alone. Our class III signal importantly excludes independent association at the TNF promoter polymorphism, TNF-308G/A, in our SLE cohort and provides a potentially novel locus for future genetic and functional studies.     

Determination of the loss of function complement C4 exon 29 CT insertion using a novel paralog-specific assay in healthy UK and Spanish populations

Genetic variants resulting in non-expression of complement C4A and C4B genes are common in healthy European populations and have shown association with a number of diseases, most notably the autoimmune disease, systemic lupus erythematosus. The most frequent cause of a C4 "null" allele, following that of C4 gene copy number variation (CNV), is a non-sense mutation arising from a 2 bp CT insertion into codon 1232 of exon 29. Previous attempts to accurately genotype this polymorphism have not been amenable to high-throughput typing, and have been confounded by failure to account for CNV at this locus, as well as by inability to distinguish between paralogs. We have developed a novel, high-throughput, paralog-specific assay to detect the presence and copy number of this polymorphism. We have genotyped healthy cohorts from the United Kingdom (UK) and Spain. Overall, 30/719 (4.17%) individuals from the UK cohort and 8/449 (1.78%) individuals from the Spanish cohort harboured the CT insertion in a C4A gene. A single Spanish individual possessed a C4B CT insertion. There is weak correlation between the C4 CT insertion and flanking MHC polymorphism. Therefore it is important to note that, as with C4 gene CNV, disease-association due to this variant will be missed by current SNP-based genome-wide association strategies.     

Imputing Variants in HLA-DR Beta Genes Reveals That HLA-DRB1 Is Solely Associated with Rheumatoid Arthritis and Systemic Lupus Erythematosus

The genetic association of HLA-DRB1 with rheumatoid arthritis (RA) and systemic lupus erythematosus (SLE) is well documented, but association with other HLA-DR beta genes (HLA-DRB3, HLA-DRB4 and HLA-DRB5) has not been thoroughly studied, despite their similar functions and chromosomal positions. We examined variants in all functional HLA-DR beta genes in RA and SLE patients and controls, down to the amino-acid level, to better understand disease association with the HLA-DR locus. To this end, we improved an existing HLA reference panel to impute variants in all protein-coding HLA-DR beta genes. Using the reference panel, HLA variants were inferred from high-density SNP data of 9,271 RA-control subjects and 5,342 SLE-control subjects. Disease association tests were performed by logistic regression and log-likelihood ratio tests. After imputation using the newly constructed HLA reference panel and statistical analysis, we observed that HLA-DRB1 variants better accounted for the association between MHC and susceptibility to RA and SLE than did the other three HLA-DRB variants. Moreover, there were no secondary effects in HLA-DRB3, HLA-DRB4, or HLA-DRB5 in RA or SLE. Of all the HLA-DR beta chain paralogs, those encoded by HLA-DRB1 solely or dominantly influence susceptibility to RA and SLE.     

High-Density SNP Screening of the Major Histocompatibility Complex in Systemic Lupus Erythematosus Demonstrates Strong Evidence for Independent Susceptibility Regions

A substantial genetic contribution to systemic lupus erythematosus (SLE) risk is conferred by major histocompatibility complex (MHC) gene(s) on chromosome 6p21. Previous studies in SLE have lacked statistical power and genetic resolution to fully define MHC influences. We characterized 1,610 Caucasian SLE cases and 1,470 parents for 1,974 MHC SNPs, the highly polymorphic HLA-DRB1 locus, and a panel of ancestry informative markers. Single-marker analyses revealed strong signals for SNPs within several MHC regions, as well as with HLA-DRB1 (global p = 9.99×10⁻¹⁶). The most strongly associated DRB1 alleles were: *0301 (odds ratio, OR = 2.21, p = 2.53×10⁻¹²), *1401 (OR = 0.50, p = 0.0002), and *1501 (OR = 1.39, p = 0.0032). The MHC region SNP demonstrating the strongest evidence of association with SLE was rs3117103, with OR = 2.44 and p = 2.80×10⁻¹³. Conditional haplotype and stepwise logistic regression analyses identified strong evidence for association between SLE and the extended class I, class I, class III, class II, and the extended class II MHC regions. Sequential removal of SLE–associated DRB1 haplotypes revealed independent effects due to variation within OR2H2 (extended class I, rs362521, p = 0.006), CREBL1 (class III, rs8283, p = 0.01), and DQB2 (class II, rs7769979, p = 0.003, and rs10947345, p = 0.0004). Further, conditional haplotype analyses demonstrated that variation within MICB (class I, rs3828903, p = 0.006) also contributes to SLE risk independent of HLA-DRB1*0301. Our results for the first time delineate with high resolution several MHC regions with independent contributions to SLE risk. We provide a list of candidate variants based on biologic and functional considerations that may be causally related to SLE risk and warrant further investigation.     

Diversity of Extended HLA-DRB1 Haplotypes in the Finnish Population

The Major Histocompatibility Complex (MHC, 6p21) codes for traditional HLA and other host response related genes. The polymorphic HLA-DRB1 gene in MHC Class II has been associated with several complex diseases. In this study we focus on MHC haplotype structures in the Finnish population. We explore the variability of extended HLA-DRB1 haplotypes in relation to the other traditional HLA genes and a selected group of MHC class III genes. A total of 150 healthy Finnish individuals were included in the study. Subjects were genotyped for HLA alleles (HLA-A, -B, -DRB1, -DQB1, and -DPB1). The polymorphism of TNF, LTA, C4, BTNL2 and HLA-DRA genes was studied with 74 SNPs (single nucleotide polymorphism). The C4A and C4B gene copy numbers and a 2-bp silencing insertion at exon 29 in C4A gene were analysed with quantitative genomic realtime-PCR. The allele frequencies for each locus were calculated and haplotypes were constructed using both the traditional HLA alleles and SNP blocks. The most frequent Finnish A∼B∼DR -haplotype, uncommon in elsewhere in Europe, was A*03∼B*35∼DRB1*01∶01. The second most common haplotype was a common European ancestral haplotype AH 8.1 (A*01∼B*08∼DRB1*03∶01). Extended haplotypes containing HLA-B, TNF block, C4 and HLA-DPB1 strongly increased the number of HLA-DRB1 haplotypes showing variability in the extended HLA-DRB1 haplotype structures. On the contrary, BTNL2 block and HLA-DQB1 were more conserved showing linkage with the HLA-DRB1 alleles. We show that the use of HLA-DRB1 haplotypes rather than single HLA-DRB1 alleles is advantageous when studying the polymorphisms and LD patters of the MHC region. For disease association studies the HLA-DRB1 haplotypes with various MHC markers allows us to cluster haplotypes with functionally important gene variants such as C4 deficiency and cytokines TNF and LTA, and provides hypotheses for further assessment. Our study corroborates the importance of studying population-specific MHC haplotypes.     

MHC region and risk of systemic lupus erythematosus in African American women

The major histocompatibility complex (MHC) on chromosome 6p21 is a key contributor to the genetic basis of systemic lupus erythematosus (SLE). Although SLE affects African Americans disproportionately compared to European Americans, there has been no comprehensive analysis of the MHC region in relationship to SLE in African Americans. We conducted a screening of the MHC region for 1,536 single nucleotide polymorphisms (SNPs) and the deletion of the C4A gene in a SLE case–control study (380 cases, 765 age-matched controls) nested within the prospective Black Women’s Health Study. We also genotyped 1,509 ancestral informative markers throughout the genome to estimate European ancestry to control for population stratification due to population admixture. The most strongly associated SNP with SLE was the rs9271366 (odds ratio, OR = 1.70, p = 5.6 × 10⁻⁵) near the HLA-DRB1 gene. Conditional haplotype analysis revealed three other SNPs, rs204890 (OR = 1.86, p = 1.2 × 10⁻⁴), rs2071349 (OR = 1.53, p = 1.0 × 10⁻³), and rs2844580 (OR = 1.43, p = 1.3 × 10⁻³), to be associated with SLE independent of the rs9271366 SNP. In univariate analysis, the OR for the C4A deletion was 1.38, p = 0.075, but after simultaneous adjustment for the other four SNPs the odds ratio was 1.01, p = 0.98. A genotype score combining the four newly identified SNPs showed an additive risk according to the number of high-risk alleles (OR = 1.67 per high-risk allele, p < 0.0001). Our strongest signal, the rs9271366 SNP, was also associated with higher risk of SLE in a previous Chinese genome-wide association study (GWAS). In addition, two SNPs found in a GWAS of European ancestry women were confirmed in our study, indicating that African Americans share some genetic risk factors for SLE with European and Chinese subjects. In summary, we found four independent signals in the MHC region associated with risk of SLE in African American women.     

Several regions in the major histocompatibility complex confer risk for anti-CCP-antibody positive rheumatoid arthritis, independent of the DRB1 locus

Recent evidence suggests that additional risk loci for RA are present in the major histocompatibility complex (MHC), independent of the class II HLA-DRB1 locus. We have now tested a total of 1,769 SNPs across 7.5Mb of the MHC located from 6p22.2 (26.03 Mb) to 6p21.32 (33.59 Mb) derived from the Illumina 550K Beadchip (Illumina, San Diego, CA, USA). For an initial analysis in the whole dataset (869 RA CCP + cases, 1,193 controls), the strongest association signal was observed in markers near the HLA-DRB1 locus, with additional evidence for association extending out into the Class I HLA region. To avoid confounding that may arise due to linkage disequilibrium with DRB1 alleles, we analyzed a subset of the data by matching cases and controls by DRB1 genotype (both alleles matched 1:1), yielding a set of 372 cases with 372 controls. This analysis revealed the presence of at least two regions of association with RA in the Class I region, independent of DRB1 genotype. SNP alleles found on the conserved A1-B8-DR3 (8.1) haplotype show the strongest evidence of positive association (P ~ 0.00005) clustered in the region around the HLA-C locus. In addition, we identified risk alleles that are not present on the 8.1 haplotype, with maximal association signals (P ~ 0.001-0.0027) located near the ZNF311 locus. This latter association is enriched in DRB1*0404 individuals. Finally, several additional association signals were found in the extreme centromeric portion of the MHC, in regions containing the DOB1, TAP2, DPB1, and COL11A2 genes. These data emphasize that further analysis of the MHC is likely to reveal genetic risk factors for rheumatoid arthritis that are independent of the DRB1 shared epitope alleles.     

Deficiency of human complement protein C4 due to identical frameshift mutations in the C4A and C4B genes

The complement protein C4, encoded by two genes (C4A and C4B) on chromosome 6p, is the most polymorphic among the MHC III gene products. We investigated the molecular basis of C4 deficiency in a Finnish woman with systemic lupus erythematosus. C4-specific mRNA was present at low concentrations in C4-deficient (C4D) patient fibroblasts, but no pro-C4 protein was detected. This defect in C4 expression was specific in that synthesis of two other complement proteins was normal. Analysis of genomic DNA showed that the proposita had both deleted and nonexpressed C4 genes. Each of her nonexpressed genes, a C4A null gene inherited from the mother, a C4A null gene, and a C4B null gene inherited from the father, all contained an identical 2-bp insertion (TC) after nucleotide 5880 in exon 29, providing the first confirmatory proof of the C4B pseudogene. This mutation has been previously found only in C4A null genes. Although the exon 29/30 junction is spliced accurately, this frameshift mutation generates a premature stop at codon 3 in exon 30. These truncated C4A and C4B gene products were confirmed through RT-PCR and sequence analysis. Among the possible genetic mechanisms that produce identical mutations is both genes, the most likely is a mutation in C4A followed by a gene conversion to generate the mutated C4B allele.     

Gene copy-number variation and associated polymorphisms of complement component C4 in human systemic lupus erythematosus (SLE): low copy number is a risk factor for and high copy number is a protective factor against SLE susceptibility in European Americans

Interindividual gene copy-number variation (CNV) of complement component C4 and its associated polymorphisms in gene size (long and short) and protein isotypes (C4A and C4B) probably lead to different susceptibilities to autoimmune disease. We investigated the C4 gene CNV in 1,241 European Americans, including patients with systemic lupus erythematosus (SLE), their first-degree relatives, and unrelated healthy subjects, by definitive genotyping and phenotyping techniques. The gene copy number (GCN) varied from 2 to 6 for total C4, from 0 to 5 for C4A, and from 0 to 4 for C4B. Four copies of total C4, two copies of C4A, and two copies of C4B were the most common GCN counts, but each constituted only between one-half and three-quarters of the study populations. Long C4 genes were strongly correlated with C4A (R=0.695; P<.0001). Short C4 genes were correlated with C4B (R=0.437; P<.0001). In comparison with healthy subjects, patients with SLE clearly had the GCN of total C4 and C4A shifting to the lower side. The risk of SLE disease susceptibility significantly increased among subjects with only two copies of total C4 (patients 9.3%; unrelated controls 1.5%; odds ratio [OR] = 6.514; P=.00002) but decreased in those with > or =5 copies of C4 (patients 5.79%; controls 12%; OR=0.466; P=.016). Both zero copies (OR=5.267; P=.001) and one copy (OR=1.613; P=.022) of C4A were risk factors for SLE, whereas > or =3 copies of C4A appeared to be protective (OR=0.574; P=.012). Family-based association tests suggested that a specific haplotype with a single short C4B in tight linkage disequilibrium with the -308A allele of TNFA was more likely to be transmitted to patients with SLE. This work demonstrates how gene CNV and its related polymorphisms are associated with the susceptibility to a human complex disease.     

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).     

A comprehensive analysis of shared loci between systemic lupus erythematosus (SLE) and sixteen autoimmune diseases reveals limited genetic overlap

In spite of the well-known clustering of multiple autoimmune disorders in families, analyses of specific shared genes and polymorphisms between systemic lupus erythematosus (SLE) and other autoimmune diseases (ADs) have been limited. Therefore, we comprehensively tested autoimmune variants for association with SLE, aiming to identify pleiotropic genetic associations between these diseases. We compiled a list of 446 non–Major Histocompatibility Complex (MHC) variants identified in genome-wide association studies (GWAS) of populations of European ancestry across 17 ADs. We then tested these variants in our combined Caucasian SLE cohorts of 1,500 cases and 5,706 controls. We tested a subset of these polymorphisms in an independent Caucasian replication cohort of 2,085 SLE cases and 2,854 controls, allowing the computation of a meta-analysis between all cohorts. We have uncovered novel shared SLE loci that passed multiple comparisons adjustment, including the VTCN1 (rs12046117, P = 2.02×10⁻⁰⁶) region. We observed that the loci shared among the most ADs include IL23R, OLIG3/TNFAIP3, and IL2RA. Given the lack of a universal autoimmune risk locus outside of the MHC and variable specificities for different diseases, our data suggests partial pleiotropy among ADs. Hierarchical clustering of ADs suggested that the most genetically related ADs appear to be type 1 diabetes with rheumatoid arthritis and Crohn''s disease with ulcerative colitis. These findings support a relatively distinct genetic susceptibility for SLE. For many of the shared GWAS autoimmune loci, we found no evidence for association with SLE, including IL23R. Also, several established SLE loci are apparently not associated with other ADs, including the ITGAM-ITGAX and TNFSF4 regions. This study represents the most comprehensive evaluation of shared autoimmune loci to date, supports a relatively distinct non–MHC genetic susceptibility for SLE, provides further evidence for previously and newly identified shared genes in SLE, and highlights the value of studies of potentially pleiotropic genes in autoimmune diseases.     

Complete complement components C4A and C4B deficiencies in human kidney diseases and systemic lupus erythematosus

Although a heterozygous deficiency of either complement component C4A or C4B is common, and each has a frequency of approximately 20% in a Caucasian population, complete deficiencies of both C4A and C4B proteins are extremely rare. In this paper the clinical courses for seven complete C4 deficiency patients are described in detail, and the molecular defects for complete C4 deficiencies are elucidated. Three patients with homozygous HLA A24 Cw7 B38 DR13 had systemic lupus erythematosus, mesangial glomerulonephritis, and severe skin lesions or membranous nephropathy. Immunofixation, genomic restriction fragment length polymorphisms, and pulsed field gel electrophoresis experiments revealed the presence of monomodular RP-C4-CYP21-TNX (RCCX) modules, each containing a solitary, long C4A mutant gene. Sequencing of the mutant C4A genes revealed a 2-bp, GT deletion in exon 13 that leads to protein truncation. The other four patients with homozygous HLA A30 B18 DR7 had SLE, severe kidney disorders including mesangial or membranoproliferative glomerulonephritis, and/or Henoch Schoenlein purpura. Molecular genetic analyses revealed an unusual RCCX structure with two short C4B mutant genes, each followed by an intact gene for steroid 21-hydroxylase. Nine identical, intronic mutations were found in each mutant C4B. In particular, the 8127 g-->a mutation present at the donor site of intron 28 may cause an RNA splice defect. Analyses of 12 complete C4 deficiency patients revealed two hot spots of deleterious mutations: one is located at exon 13, the others within a 2.6-kb genomic region spanning exons 20-29. Screening of these mutations may facilitate epidemiologic studies of C4 in infectious, autoimmune, and kidney diseases.     

Alpha-defensin DEFA1A3 gene copy number variation in Asians and its genetic association study in Chinese systemic lupus erythematosus patients

Neutrophil extracellular traps (NETs) were closely associated with activation of type I interferon (IFN) pathway in systemic lupus erythematosus (SLE). We aimed to study the genetic basis of NETs-DEFA1A3 copy number variations (CNV) in SLE and HapMap CHB + JPT populations by quantitative real-time PCR and whole genome sequences data. DEFA1A3 CNs did not differ significantly between SLE patients and controls. DEFA1A3 CNs ranged from 3 to 11 in CHB and 4 to 16 in JPT. The median of DEFA1A3 CNV of CHB (6 copies) was significantly lower than that of JPT (9 copies). Associations of genotype of tag SNP rs2738113 with DEFA1A3 CNs and mRNA expression of IFNα were observed in CHB and JPT populations. Our data provided a genetic reference of DEFA1A3 CNV for further studies and suggested that the genetic pathogenesis of NETs, as well as DEFA1A3 in SLE should be further evaluated, specially in different populations.     

Copy number analysis of complement C4A, C4B and C4A silencing mutation by real-time quantitative polymerase chain reaction

Low protein levels and copy number variation (CNV) of the fourth component of human complement (C4A and C4B) have been associated with various diseases. High-throughput methods for analysing C4 CNV are available, but they commonly do not detect the most common C4A mutation, a silencing CT insertion (CTins) leading to low protein levels. We developed a SYBR? Green labelled real-time quantitative polymerase chain reaction (qPCR) with a novel concentration range approach to address C4 CNV and deficiencies due to CTins. This method was validated in three sample sets and applied to over 1600 patient samples. CTins caused C4A deficiency in more than 70% (76/105) of the carriers. Twenty per cent (76/381) of patients with a C4A deficiency would have been erroneously recorded as having none, if the CTins had not been assessed. C4A deficiency was more common in patients than a healthy reference population, (OR?=?1.60, 95%CI?=?1.02-2.52, p?=?0.039). The number of functional C4 genes can be straightforwardly analyzed by real-time qPCR, also with SYBR? Green labelling. Determination of CTins increases the frequency of C4A deficiency and thus helps to elucidate the genotypic versus phenotypic disease associations.     

Uncoupling the roles of HLA-DRB1 and HLA-DRB5 genes in multiple sclerosis

Genetic susceptibility to multiple sclerosis (MS) is associated with the MHC located on chromosome 6p21. This signal maps primarily to a 1-Mb region encompassing the HLA class II loci, and it segregates often with the HLA-DQB1*0602, -DQA1*0102, -DRB1*1501, -DRB5*0101 haplotype. However, the identification of the true predisposing gene or genes within the susceptibility haplotype has been handicapped by the strong linkage disequilibrium across the locus. African Americans have greater MHC haplotypic diversity and distinct patterns of linkage disequilibrium, which make this population particularly informative for fine mapping efforts. The purpose of this study was to establish the telomeric boundary of the HLA class II region affecting susceptibility to MS by assessing genetic association with the neighboring HLA-DRB5 gene as well as seven telomeric single nucleotide polymorphisms in a large, well-characterized African American dataset. Rare DRB5*null individuals were previously described in African populations. Although significant associations with both HLA-DRB1 and HLA-DRB5 loci were present, HLA-DRB1*1503 was associated with MS in the absence of HLA-DRB5, providing evidence for HLA-DRB1 as the primary susceptibility gene. Interestingly, the HLA-DRB5*null subjects appear to be at increased risk for developing secondary progressive MS. Thus, HLA-DRB5 attenuates MS severity, a finding consistent with HLA-DRB5's proposed role as a modifier in experimental autoimmune encephalomyelitis. Additionally, conditional haplotype analysis revealed a susceptibility signal at the class III AGER locus independent of DRB1. The data underscore the power of the African American MS dataset to identify disease genes by association in a region of high linkage disequilibrium.     

Risk alleles for systemic lupus erythematosus in a large case-control collection and associations with clinical subphenotypes

Systemic lupus erythematosus (SLE) is a genetically complex disease with heterogeneous clinical manifestations. Recent studies have greatly expanded the number of established SLE risk alleles, but the distribution of multiple risk alleles in cases versus controls and their relationship to subphenotypes have not been studied. We studied 22 SLE susceptibility polymorphisms with previous genome-wide evidence of association (p<5×10⁻⁸) in 1919 SLE cases from 9 independent Caucasian SLE case series and 4813 independent controls. The mean number of risk alleles in cases was 15.1 (SD 3.1) while the mean in controls was 13.1 (SD 2.8), with trend p = 4×10⁻¹²⁸. We defined a genetic risk score (GRS) for SLE as the number of risk alleles with each weighted by the SLE risk odds ratio (OR). The OR for high-low GRS tertiles, adjusted for intra-European ancestry, sex, and parent study, was 4.4 (95% CI 3.8–5.1). We studied associations of individual SNPs and the GRS with clinical manifestations for the cases: age at diagnosis, the 11 American College of Rheumatology classification criteria, and double-stranded DNA antibody (anti-dsDNA) production. Six subphenotypes were significantly associated with the GRS, most notably anti-dsDNA (OR_high-low = 2.36, p = 9e−9), the immunologic criterion (OR_high-low = 2.23, p = 3e−7), and age at diagnosis (OR_high-low = 1.45, p = 0.0060). Finally, we developed a subphenotype-specific GRS (sub-GRS) for each phenotype with more power to detect cumulative genetic associations. The sub-GRS was more strongly associated than any single SNP effect for 5 subphenotypes (the above plus hematologic disorder and oral ulcers), while single loci are more significantly associated with renal disease (HLA-DRB1, OR = 1.37, 95% CI 1.14–1.64) and arthritis (ITGAM, OR = 0.72, 95% CI 0.59–0.88). We did not observe significant associations for other subphenotypes, for individual loci or the sub-GRS. Thus our analysis categorizes SLE subphenotypes into three groups: those having cumulative, single, and no known genetic association with respect to the currently established SLE risk loci.     

Beyond C4: Analysis of the complement gene pathway shows enrichment for IQ in patients with psychotic disorders and healthy controls

Variation in cognitive performance, which strongly predicts functional outcome in schizophrenia (SZ), has been associated with multiple immune‐relevant genetic loci. These loci include complement component 4 (C4A), structural variation at which was recently associated with SZ risk and synaptic pruning during neurodevelopment and cognitive function. Here, we test whether this genetic association with cognition and SZ risk is specific to C4A, or extends more broadly to genes related to the complement system. Using a gene‐set with an identified role in “complement” function (excluding C4A), we used MAGMA to test if this gene‐set was enriched for genes associated with human intelligence and SZ risk, using genome‐wide association summary statistics (IQ; N = 269 867, SZ; N = 105 318). We followed up this gene‐set analysis with a complement gene‐set polygenic score (PGS) regression analysis in an independent data set of patients with psychotic disorders and healthy participants with cognitive and genomic data (N = 1000). Enrichment analysis suggested that genes within the complement pathway were significantly enriched for genes associated with IQ, but not SZ. In a gene‐based analysis of 90 genes, SERPING1 was the most enriched gene for the phenotype of IQ. In a PGS regression analysis, we found that a complement pathway PGS associated with IQ genome‐wide association studies statistics also predicted variation in IQ in our independent sample. This association (observed across both patients and controls) remained significant after controlling for the relationship between C4A and cognition. These results suggest a robust association between the complement system and cognitive function, extending beyond structural variation at C4A.     

The BTNL2 gene and sarcoidosis susceptibility in African Americans and Whites

The BTNL2 gene is a member of the B7 receptor family that probably functions as a T-cell costimulatory molecule. It resides in the class II major histocompatibility complex (MHC) region of chromosome 6p and has recently been associated with sarcoidosis susceptibility in a white German population. We sought to replicate the BTNL2 association in an African American family-based study population (n=219 nuclear families) and two case-control populations--one African American (n=295 pairs) and one white (n=366 pairs). Ten SNPs were detected within a 490-bp region spanning exon/intron 5 of BTNL2. Haplotype variation within this region was significantly associated with sarcoidosis in all three study populations but more so in whites (P=.0006) than in the African American case-control (P=.02) or family-based (P=.03) samples. The previously reported BTNL2 SNP with the strongest sarcoidosis association, rs2076530, was also the SNP with the strongest association in our white population (P<.0001). The A allele of rs2076530 results in a premature exon-splice site and increases risk for sarcoidosis (odds ratio=2.03; 95% confidence interval 1.32-3.12). Although rs2076530 was not associated with sarcoidosis in either African American sample, a three-locus haplotype that included rs2076530 was associated with sarcoidosis across all three study samples. Multivariable logistic regression analyses showed that BTNL2 effects are independent of human leukocyte antigen class II genes in whites but may interact antagonistically in African Americans. Our results underscore the complexity of genetic risk for sarcoidosis emanating from the MHC region.     

Molecular analysis of a novel hereditary C3 deficiency with systemic lupus erythematosus

A case of inherited homozygous complement C3 deficiency (C3D) in a patient with systemic lupus erythematosus (SLE) and the molecular basis for this deficiency are reported. A 22-year-old Japanese male was diagnosed as having SLE and his medical history revealed recurrent tonsillitis and pneumonia. He was diagnosed as having C3D because of undetectable serum C3 level. His parents were consanguineous. Sequence analysis of C3D cDNA revealed a homozygous deletion of exon 39 (84bp). A single base substitution (AG to GG) in the 3'-splice acceptor site of intron 38 was identified by sequencing the genomic DNA. Expression of C3Delta(ex39) cDNA, the C3cDNA lacking exon 39, in COS-7 cells revealed that C3Delta(ex39) was retained in endoplasmic reticulum-Golgi intermediate compartment because of defective secretion. These data indicate that a novel AG-->GG 3'-splice acceptor site mutation in intron 38 caused aberrant splicing of exon 39, resulting in defective secretion of C3.     

Identification of a novel risk locus for multiple sclerosis at 13q31.3 by a pooled genome-wide scan of 500,000 single nucleotide polymorphisms

Multiple sclerosis is a chronic inflammatory demyelinating disease of the central nervous system with an important genetic component and strongest association driven by the HLA genes. We performed a pooling-based genome-wide association study of 500,000 SNPs in order to find new loci associated with the disease. After applying several criteria, 320 SNPs were selected from the microarrays and individually genotyped in a first and independent Spanish Caucasian replication cohort. The 8 most significant SNPs validated in this cohort were also genotyped in a second US Caucasian replication cohort for confirmation. The most significant association was obtained for SNP rs3129934, which neighbors the HLA-DRB/DQA loci and validates our pooling-based strategy. The second strongest association signal was found for SNP rs1327328, which resides in an unannotated region of chromosome 13 but is in linkage disequilibrium with nearby functional elements that may play important roles in disease susceptibility. This region of chromosome 13 has not been previously identified in MS linkage genome screens and represents a novel risk locus for the disease.