F. Yesim Demirci

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.

Rare coding variants in the phospholipase D3 gene confer risk for Alzheimer's disease

Genome-wide association studies (GWAS) have identified several risk variants for late-onset Alzheimers disease (LOAD). These common variants have replicable but small effects on LOAD risk and generally do not have obvious functional effects. Low-frequency coding variants, not detected by GWAS, are predicted to include functional variants with larger effects on risk. To identify low-frequency coding variants with large effects on LOAD risk, we carried out whole-exome sequencing (WES) in 14 large LOAD families and follow-up analyses of the candidate variants in several large LOAD case–control data sets. A rare variant in PLD3 (phospholipase D3; Val232Met) segregated with disease status in two independent families and doubled risk for Alzheimer’s disease in seven independent case–control series with a total of more than 11,000 cases and controls of European descent. Gene-based burden analyses in 4,387 cases and controls of European descent and 302 African American cases and controls, with complete sequence data for PLD3, reveal that several variants in this gene increase risk for Alzheimer’s disease in both populations. PLD3 is highly expressed in brain regions that are vulnerable to Alzheimer’s disease pathology, including hippocampus and cortex, and is expressed at significantly lower levels in neurons from Alzheimer’s disease brains compared to control brains. Overexpression of PLD3 leads to a significant decrease in intracellular amyloid-β precursor protein (APP) and extracellular Aβ42 and Aβ40 (the 42- and 40-residue isoforms of the amyloid-β peptide), and knockdown of PLD3 leads to a significant increase in extracellular Aβ42 and Aβ40. Together, our genetic and functional data indicate that carriers of PLD3 coding variants have a twofold increased risk for LOAD and that PLD3 influences APP processing. This study provides an example of how densely affected families may help to identify rare variants with large effects on risk for disease or other complex traits.

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.

Differential Genetic Associations for Systemic Lupus Erythematosus Based on Anti–dsDNA Autoantibody Production

Systemic lupus erythematosus (SLE) is a clinically heterogeneous, systemic autoimmune disease characterized by autoantibody formation. Previously published genome-wide association studies (GWAS) have investigated SLE as a single phenotype. Therefore, we conducted a GWAS to identify genetic factors associated with anti–dsDNA autoantibody production, a SLE–related autoantibody with diagnostic and clinical importance. Using two independent datasets, over 400,000 single nucleotide polymorphisms (SNPs) were studied in a total of 1,717 SLE cases and 4,813 healthy controls. Anti–dsDNA autoantibody positive (anti–dsDNA +, n = 811) and anti–dsDNA autoantibody negative (anti–dsDNA –, n = 906) SLE cases were compared to healthy controls and to each other to identify SNPs associated specifically with these SLE subtypes. SNPs in the previously identified SLE susceptibility loci STAT4, IRF5, ITGAM, and the major histocompatibility complex were strongly associated with anti–dsDNA + SLE. Far fewer and weaker associations were observed for anti–dsDNA – SLE. For example, rs7574865 in STAT4 had an OR for anti–dsDNA + SLE of 1.77 (95% CI 1.57–1.99, p = 2.0E-20) compared to an OR for anti–dsDNA – SLE of 1.26 (95% CI 1.12–1.41, p = 2.4E-04), with pheterogeneity<0.0005. SNPs in the SLE susceptibility loci BANK1, KIAA1542, and UBE2L3 showed evidence of association with anti–dsDNA + SLE and were not associated with anti–dsDNA – SLE. In conclusion, we identified differential genetic associations with SLE based on anti–dsDNA autoantibody production. Many previously identified SLE susceptibility loci may confer disease risk through their role in autoantibody production and be more accurately described as autoantibody propensity loci. Lack of strong SNP associations may suggest that other types of genetic variation or non-genetic factors such as environmental exposures have a greater impact on susceptibility to anti–dsDNA – SLE.

X-linked cone-rod dystrophy (locus COD1): identification of mutations in RPGR exon ORF15.

X-linked cone-rod dystrophy (COD1) is a retinal disease that primarily affects the cone photoreceptors; the disease was originally mapped to a limited region of Xp11.4. We evaluated the three families from our original study with new markers and clinically reassessed all key recombinants; we determined that the critical intervals in families 2 and 3 overlapped the RP3 locus and that a status change (from affected to probably unaffected) of a key recombinant individual in family 1 also reassigned the disease locus to include RP3 as well. Mutation analysis of the entire RPGR coding region identified two different 2-nucleotide (nt) deletions in ORF15, in family 2 (delAG) and in families 1 and 3 (delGG), both of which result in a frameshift leading to altered amino acid structure and early termination. In addition, an independent individual with X-linked cone-rod dystrophy demonstrated a 1-nt insertion (insA) in ORF15. The presence of three distinct mutations associated with the same disease phenotype provides strong evidence that mutations in RPGR exon ORF15 are responsible for COD1. Genetic heterogeneity was observed in three other families, including the identification of an in-frame 12-nt deletion polymorphism in ORF15 that did not segregate with the disease in one of these families.

Association of CLU and PICALM variants with Alzheimer's disease

Two recent large genome-wide association studies have reported significant associations in the CLU (APOJ), CR1, and PICALM genes with the risk of Alzheimers disease (AD). In order to replicate these findings, we examined 7 single nucleotide polymorphisms (SNPs) most significantly implicated by these studies in a large case-control sample comprising 2707 individuals. Principle components analysis revealed no population substructure in our sample. While no association was observed with CR1 SNPs (p = 0.30-0.457), a trend of association was seen with the PICALM (p = 0.071-0.086) and CLU (p = 0.148-0.258) SNPs. A meta-analysis of 3 studies revealed significant associations with all 3 genes. Our data from an independent and large case-control sample suggest that these gene regions should be followed up by comprehensive resequencing to find functional variants.

Lupus Nephritis Susceptibility Loci in Women with Systemic Lupus Erythematosus

Lupus nephritis is a manifestation of SLE resulting from glomerular immune complex deposition and inflammation. Lupus nephritis demonstrates familial aggregation and accounts for significant morbidity and mortality. We completed a meta-analysis of three genome-wide association studies of SLE to identify lupus nephritis-predisposing loci. Through genotyping and imputation, >1.6 million markers were assessed in 2000 unrelated women of European descent with SLE (588 patients with lupus nephritis and 1412 patients with lupus without nephritis). Tests of association were computed using logistic regression adjusting for population substructure. The strongest evidence for association was observed outside the MHC and included markers localized to 4q11-q13 (PDGFRA, GSX2; P=4.5×10(-7)), 16p12 (SLC5A11; P=5.1×10(-7)), 6p22 (ID4; P=7.4×10(-7)), and 8q24.12 (HAS2, SNTB1; P=1.1×10(-6)). Both HLA-DR2 and HLA-DR3, two well established lupus susceptibility loci, showed evidence of association with lupus nephritis (P=0.06 and P=3.7×10(-5), respectively). Within the class I region, rs9263871 (C6orf15-HCG22) had the strongest evidence of association with lupus nephritis independent of HLA-DR2 and HLA-DR3 (P=8.5×10(-6)). Consistent with a functional role in lupus nephritis, intra-renal mRNA levels of PDGFRA and associated pathway members showed significant enrichment in patients with lupus nephritis (n=32) compared with controls (n=15). Results from this large-scale genome-wide investigation of lupus nephritis provide evidence of multiple biologically relevant lupus nephritis susceptibility loci.

No association between CALHM1 variation and risk of Alzheimer disease.

A polymorphism in the calcium homeostasis modulator 1 gene (CALHM1) has recently been associated with risk of late‐onset Alzheimer disease. We examined this variant (rs2986017) in 945 Caucasian Americans with late‐onset Alzheimer disease and 875 age‐matched Caucasian American controls. No association with risk of late‐onset Alzheimer disease (p=0.368 for genotypes; p=0.796 for alleles) was observed in our sample. However, a potential modest association of minor allele homozygosity (TT) with an earlier age‐at‐onset was seen (p=0.034).

Functional significance of lipoprotein lipase HindIII polymorphism associated with the risk of coronary artery disease

Lipoprotein lipase (LPL) plays a pivotal role in lipid metabolism by hydrolyzing triglyceride (TG)-rich lipoprotein particles. Abnormalities in normal LPL function are associated with the risk of coronary artery disease (CAD). A number of genetic variants have been identified in the LPL gene that affects different functions of the LPL protein. A common HindIII polymorphism in intron 8 (T/G) of the LPL gene has been found to be associated with altered plasma TG and HDL-cholesterol, and CAD risk in several studies, but its functional significance is unknown. It has been shown that certain intronic sequence contain regulatory elements that are important for transcription and translational regulation of a gene. In this study we tested the hypothesis that this polymorphism affects the binding site of a transcription factor that regulates the transcription of LPL gene. Electrophoretic mobility shift assays (EMSAs) revealed that the HindIII site binds to a transcription factor and that the mutant allele has lower binding affinity than the wild type allele. Transcription assays containing the entire intron 8 sequence along with full-length human LPL promoter were carried out in COS-1 and human vascular smooth muscle cells. The mutant allele was associated with significantly decreased luciferase expression level compared to the wild type allele in both the muscle (3.394+/-0.022 vs. 4.184+/-0.028; P=4.7 x 10(-6)) and COS-1 (11.603+/-0.409 vs. 14.373+/-1.096; P<0.0001) cells. In conclusion, this study demonstrates for the first time that the polymorphic HindIII site in the LPL gene is functional because it affects the binding of a transcription factor and it also has an impact on LPL expression.

Connecting the Dots: Potential of Data Integration to Identify Regulatory SNPs in Late-Onset Alzheimer's Disease GWAS Findings

Late-onset Alzheimers disease (LOAD) is a multifactorial disorder with over twenty loci associated with disease risk. Given the number of genome-wide significant variants that fall outside of coding regions, it is possible that some of these variants alter some function of gene expression rather than tagging coding variants that alter protein structure and/or function. RegulomeDB is a database that annotates regulatory functions of genetic variants. In this study, we utilized RegulomeDB to investigate potential regulatory functions of lead single nucleotide polymorphisms (SNPs) identified in five genome-wide association studies (GWAS) of risk and age-at onset (AAO) of LOAD, as well as SNPs in LD (r2≥0.80) with the lead GWAS SNPs. Of a total 614 SNPs examined, 394 returned RegulomeDB scores of 1–6. Of those 394 variants, 34 showed strong evidence of regulatory function (RegulomeDB score <3), and only 3 of them were genome-wide significant SNPs (ZCWPW1/rs1476679, CLU/rs1532278 and ABCA7/rs3764650). This study further supports the assumption that some of the non-coding GWAS SNPs are true associations rather than tagged associations and demonstrates the application of RegulomeDB to GWAS data.