Faculty Bibliography

A large portion of common variant loci associated with genetic risk for schizophrenia reside within noncoding sequence of unknown function. Here, we demonstrate promoter and enhancer enrichment in schizophrenia variants associated with expression quantitative trait loci (eQTL). The enrichment is greater when functional annotations derived from the human brain are used relative to peripheral tissues. Regulatory trait concordance analysis ranked genes within schizophrenia genome-wide significant loci for a potential functional role, based on colocalization of a risk SNP, eQTL, and regulatory element sequence. We identified potential physical interactions of noncontiguous proximal and distal regulatory elements. This was verified in prefrontal cortex and -induced pluripotent stem cell-derived neurons for the L-type calcium channel (CACNA1C) risk locus. Our findings point to a functional link between schizophrenia-associated noncoding SNPs and 3D genome architecture associated with chromosomal loopings and transcriptional regulation in the brain.

BACKGROUND: Joint destruction is a hallmark of autoantibody-positive rheumatoid arthritis (RA), though the severity is highly variable between patients. The processes underlying these interindividual differences are incompletely understood. METHODS: We performed a genome-wide association study on the radiological progression rate in 384 autoantibody-positive patients with RA. In stage-II 1557 X-rays of 301 Dutch autoantibody-positive patients with RA were studied and in stage-III 861 X-rays of 742 North American autoantibody-positive patients with RA. Sperm-Associated Antigen 16 (SPAG16) expression in RA synovium and fibroblast-like synoviocytes (FLS) was examined using Real-Time Quantitative Polymerase Chain Reaction (RT-qPCR) and immunohistochemistry. FLS secrete metalloproteinases that degrade cartilage and bone. SPAG16 genotypes were related to matrix metalloproteinase (MMP)-3 and MMP-1 expression by FLS in vitro and MMP-3 production ex vivo. RESULTS: A cluster of single nucleotide polymorphisms (SNPs) at 2q34, located at SPAG16, associated with the radiological progression rate; rs7607479 reached genome-wide significance. A protective role of rs7607479 was replicated in European and North American patients with RA. Per minor allele, patients had a 0.78-fold (95% CI 0.67 to 0.91) progression rate over 7 years. mRNA and protein expression of SPAG16 in RA synovium and FLS was verified. FLS carrying the minor allele secreted less MMP-3 (p=1.60x10-2). Furthermore, patients with RA carrying the minor allele had lower serum levels of MMP-3 (p=4.28x10-2). In a multivariate analysis on rs7607479 and MMP-3, only MMP-3 associated with progression (p=2.77x10-4), suggesting that the association between SPAG16-rs7607479 and joint damage is mediated via an effect on MMP-3 secretion. CONCLUSIONS: Genetic and functional analyses indicate that SPAG16 influences MMP-3 regulation and protects against joint destruction in autoantibody-positive RA. These findings could enhance risk stratification in autoantibody-positive RA.

The aim of this study was to identify genetic variants associated with rheumatoid arthritis (RA) risk in black South Africans. Black South African RA patients (n=263) were compared to healthy controls (n=374). Genotyping was performed using the Immunochip and four digit high resolution HLA typing was performed by DNA sequencing of exon 2. Standard quality control measures were implemented on the data. The strongest associations were in the intergenic region between the HLA DRB1 and HLA DQA1 loci. After conditioning on HLA DRB1 alleles the effect in the rest of the extended MHC diminished. Non-HLA single nucleotide polymorphisms (SNPs) in the intergenic regions LOC389203|RBPJ; LOC100131131|IL1R1; KIAA1919|REV3L; LOC643749|TRAF3IP2; and SNPs in the intron and UTR of IRF1 and the intronic region of ICOS and KIAA1542 showed association with RA (p < 5x10-5). Of the SNPs previously associated with RA in Caucasians, one SNP, rs874040, locating to the intergenic region LOC389203|RBPJ was replicated in this study. None of the variants in the PTPN22 gene was significantly associated. The sero-positive subgroups showed similar results to the overall cohort. The effects observed across the HLA region are most likely due to HLA-DRB1, and secondary effects in the extended MHC cannot be detected. Seven non-HLA loci are associated with RA in black South Africans. Similar to Caucasians, the intergenic region between LOC38920 and RBPJ is associated with RA in this population. The strong association of the R620W variant of the PTPN22 gene with RA in Caucasians was not replicated since this variant was monomorphic in our study, but other SNP variants of the PTPN22 gene were also not associated with RA in black South Africans, suggesting that this locus does not play a major role in RA in this population.

Psoriasis vulgaris (PsV) risk is strongly associated with variation within the major histocompatibility complex (MHC) region, but its genetic architecture has yet to be fully elucidated. Here, we conducted a large-scale fine-mapping study of PsV risk in the MHC region in 9,247 PsV-affected individuals and 13,589 controls of European descent by imputing class I and II human leukocyte antigen (HLA) genes from SNP genotype data. In addition, we imputed sequence variants for MICA, an MHC HLA-like gene that has been associated with PsV, to evaluate association at that locus as well. We observed that HLA-C *06:02 demonstrated the lowest p value for overall PsV risk (p = 1.7 x 10-364). Stepwise analysis revealed multiple HLA-C *06:02-independent risk variants in both class I and class II HLA genes for PsV susceptibility (HLA-C *12:03, HLA-B amino acid positions 67 and 9, HLA-A amino acid position 95, and HLA-DQalpha1 amino acid position 53; p < 5.0 x 10-8), but no apparent risk conferred by MICA. We further evaluated risk of two major clinical subtypes of PsV, psoriatic arthritis (PsA; n = 3,038) and cutaneous psoriasis (PsC; n = 3,098). We found that risk heterogeneity between PsA and PsC might be driven by HLA-B amino acid position 45 (pomnibus = 2.2 x 10-11), indicating that different genetic factors underlie the overall risk of PsV and the risk of specific PsV subphenotypes. Our study illustrates the value of high-resolution HLA and MICA imputation for fine mapping causal variants in the MHC.

Despite progress in defining human leukocyte antigen (HLA) alleles for anti-citrullinated-protein-autoantibody-positive (ACPA(+)) rheumatoid arthritis (RA), identifying HLA alleles for ACPA-negative (ACPA(-)) RA has been challenging because of clinical heterogeneity within clinical cohorts. We imputed 8,961 classical HLA alleles, amino acids, and SNPs from Immunochip data in a discovery set of 2,406 ACPA(-) RA case and 13,930 control individuals. We developed a statistical approach to identify and adjust for clinical heterogeneity within ACPA(-) RA and observed independent associations for serine and leucine at position 11 in HLA-DRbeta1 (p = 1.4 x 10(-13), odds ratio [OR] = 1.30) and for aspartate at position 9 in HLA-B (p = 2.7 x 10(-12), OR = 1.39) within the peptide binding grooves. These amino acid positions induced associations at HLA-DRB1( *)03 (encoding serine at 11) and HLA-B( *)08 (encoding aspartate at 9). We validated these findings in an independent set of 427 ACPA(-) case subjects, carefully phenotyped with a highly sensitive ACPA assay, and 1,691 control subjects (HLA-DRbeta1 Ser11+Leu11: p = 5.8 x 10(-4), OR = 1.28; HLA-B Asp9: p = 2.6 x 10(-3), OR = 1.34). Although both amino acid sites drove risk of ACPA(+) and ACPA(-) disease, the effects of individual residues at HLA-DRbeta1 position 11 were distinct (p < 2.9 x 10(-107)). We also identified an association with ACPA(+) RA at HLA-A position 77 (p = 2.7 x 10(-8), OR = 0.85) in 7,279 ACPA(+) RA case and 15,870 control subjects. These results contribute to mounting evidence that ACPA(+) and ACPA(-) RA are genetically distinct and potentially have separate autoantigens contributing to pathogenesis. We expect that our approach might have broad applications in analyzing clinical conditions with heterogeneity at both major histocompatibility complex (MHC) and non-MHC regions.