Search

Background The coronavirus disease 2019 (COVID-19) pandemic, caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has demonstrated the need to share data and biospecimens broadly to optimize clinical outcomes for US military Veterans. Methods In response, the Veterans Health Administration established VA SHIELD (Science and Health Initiative to Combat Infectious and Emerging Life-threatening Diseases), a comprehensive biorepository of specimens and clinical data from affected Veterans to advance research and public health surveillance and to improve diagnostic and therapeutic capabilities. Results VA SHIELD now comprises 12 sites collecting de-identified biospecimens from US Veterans affected by SARS-CoV-2. In addition, 2 biorepository sites, a data processing center, and a coordinating center have been established under the direction of the Veterans Affairs Office of Research and Development. Phase 1 of VA SHIELD comprises 34 157 samples. Of these, 83.8% had positive tests for SARS-CoV-2, with the remainder serving as contemporaneous controls. The samples include nasopharyngeal swabs (57.9%), plasma (27.9%), and sera (12.5%). The associated clinical and demographic information available permits the evaluation of biological data in the context of patient demographics, clinical experience and management, vaccinations, and comorbidities. Conclusions VA SHIELD is representative of US national diversity with a significant potential to impact national healthcare. VA SHIELD will support future projects designed to better understand SARS-CoV-2 and other emergent healthcare crises. To the extent possible, VA SHIELD will facilitate the discovery of diagnostics and therapeutics intended to diminish COVID-19 morbidity and mortality and to reduce the impact of new emerging threats to the health of US Veterans and populations worldwide.

Preterm birth (PTB), or birth that occurs earlier than 37 weeks of gestational age, is a major contributor to infant mortality and neonatal hospitalization. Mutations in the mitochondrial genome (mtDNA) have been linked to various rare mitochondrial disorders, and may be a contributing factor in PTB given that maternal genetic factors have been strongly linked to PTB. However, to date, no study has found a conclusive connection between a particular mtDNA variant and PTB. Given the high mtDNA copy number per cell, an automated pipeline was developed for detecting mtDNA variants using low-pass whole genome sequencing (lcWGS) data. The pipeline was first validated against samples of known heteroplasmy, and then applied to 929 samples from a PTB cohort from diverse ethnic backgrounds with an average gestational age of 27.18 weeks (range: 21-30). Our new pipeline successfully identified haplogroups and a large number of mtDNA variants in this large PTB cohort, including 8 samples carrying known pathogenic variants and 47 samples carrying rare mtDNA variants. These results confirm that lcWGS can be utilized to reliably identify mtDNA variants. These mtDNA variants may make a contribution toward preterm birth in a small proportion of live births. This article is protected by copyright. All rights reserved.

Objective This study was undertaken to investigate the hypothesis that a genetic predisposition toward rheumatoid arthritis (RA) increases the risk of 10 cardiometabolic and autoimmune disorders previously associated with RA in epidemiologic studies, and to define new genetic pleiotropy present in RA. Methods Two approaches were used to test our hypothesis. First, we constructed a weighted genetic risk score (wGRS) and then examined its association with 10 prespecified disorders. Additionally, a phenome-wide association study (PheWAS) was carried out to identify potential new associations. Second, inverse variance-weighted regression (IVWR) meta-analysis was used to characterize the association between genetic susceptibility to RA and the prespecified disorders, with the results expressed as odds ratios (ORs) and 95% confidence intervals (95% CIs). Results The wGRS for RA was significantly associated with type 1 diabetes mellitus (DM) (OR 1.10 [95% CI 1.04-1.16]; P = 9.82 × 10-4 ) and multiple sclerosis (OR 0.82 [95% CI 0.77-0.88]; P = 1.73 × 10-8 ), but not with other cardiometabolic phenotypes. In the PheWAS, wGRS was also associated with an increased risk of several autoimmune phenotypes including RA, thyroiditis, and systemic sclerosis, and with a decreased risk of demyelinating disorders. In the IVWR meta-analyses, RA was significantly associated with an increased risk of type 1 DM (P = 1.15 × 10-14 ), with evidence of horizontal pleiotropy (Mendelian Randomization-Egger intercept estimate P = 0.001) likely driven by rs2476601, a PTPN22 variant. The association between type 1 DM and RA remained significant (P = 9.53 × 10-9 ) after excluding rs2476601, with no evidence of horizontal pleiotropy (intercept estimate P = 0.939). RA was also significantly associated with type 2 DM and C-reactive protein levels. These associations were driven by variation in the major histocompatibility complex region. Conclusion This study presents evidence of pleiotropy between the genetic predisposition to RA and associated phenotypes found in other autoimmune and cardiometabolic disorders, including type 1 DM.

Eosinophilic esophagitis (EoE) is a chronic allergic inflammatory disease with a complex underlying genetic etiology. Herein, we conduct whole-exome sequencing of a multigeneration EoE pedigree (discovery set) and 61 additional multiplex families with EoE (replication set). A series of rare, heterozygous, missense variants are identified in the genes encoding the desmosome-associated proteins DSP and PPL in 21% of the multiplex families. Esophageal biopsies from patients with these variants retain dilated intercellular spaces and decrease DSP and PPL expression even during disease remission. These variants affect barrier integrity, cell motility and RhoGTPase activity in esophageal epithelial cells and have increased susceptibility to calpain-14–mediated degradation. An acquired loss of esophageal DSP and PPL is present in non-familial EoE. Taken together, herein, we uncover a pathogenic role for desmosomal dysfunction in EoE, providing a deeper mechanistic understanding of tissue-specific allergic responses., Eosinophilic esophagitis (EoE) is a chronic allergic inflammatory disease with a complex underlying genetic etiology. Here, the authors identify a series of rare variants in DSP and PPL in multiplex families with EoE and uncover a pathogenic role for desmosomal dysfunction in EoE.

Rationale: Complete tracheal ring deformity (CTRD) is a rare congenital abnormality of unknown etiology characterized by circumferentially continuous or nearly-continuous cartilaginous tracheal rin...

Dysregulated transcription factors represent a major class of drug targets that mediate the abnormal expression of many critical genes involved in SLE and other autoimmune diseases. Although strong evidence suggests that natural human genetic variation affects basal and inducible gene expression, it is still a considerable challenge to establish a biological link between GWAS-identified non-coding genetic risk variants and their regulated gene targets. Here, we combine genetic data, epigenomic data, and CRISPR activation (CRISPRa) assays to screen for functional variants regulating IRF8 expression. Using CRISPR-mediated deletion and 3D chromatin structure analysis, we demonstrate that the locus containing rs2280381 is a cell-type-specific distal enhancer for IRF8 that spatially interacts with the IRF8 promoter. Further, rs2280381 mediates IRF8 expression through enhancer RNA AC092723.1, which recruits TET1 to the IRF8 promoter to modulate IRF8 expression by affecting methylation levels. The alleles of rs2280381 modulate PU.1 binding and chromatin state to differentially regulate AC092723.1 and IRF8 expression. Our work illustrates a strategy to define the functional genetic variants modulating transcription factor gene expression levels and identifies the biologic mechanism by which autoimmune diseases risk genetic variants contribute to the pathogenesis of disease.

Coronary artery disease (CAD) is a leading cause of death, yet its genetic determinants are not fully elucidated. We report a multi-ethnic genome-wide association study of CAD involving nearly a quarter of a million cases, incorporating the largest cohorts to date of Whites, Blacks, and Hispanics from the Million Veteran Program with existing studies including CARDIoGRAMplusC4D, UK Biobank, and Biobank Japan. We verify substantial and nearly equivalent heritability of CAD across multiple ancestral groups, discover 107 novel loci including the first nine on the X-chromosome, identify the first eight genome-wide significant loci among Blacks and Hispanics, and demonstrate that two common haplotypes are largely responsible for the risk stratification at the well-known 9p21 locus in most populations except those of African origin where both haplotypes are virtually absent. We identify 15 loci for angiographically derived burden of coronary atherosclerosis, which robustly overlap with the strongest and earliest loci reported to date for clinical CAD. Phenome-wide association analyses of novel loci and externally validated polygenic risk scores (PRS) augment signals from the insulin resistance cluster of risk factors and consequences, extend previously established pleiotropic associations of loci with traditional risk factors to include smoking and family history, and confirm a substantially reduced transferability of existing PRS to Blacks. Downstream integrative genomic analyses reinforce the critical role of endothelial, fibroblast, and smooth muscle cells within the coronary vessel wall in CAD susceptibility. Our study highlights the value of a multi-ethnic design in efficiently characterizing the genetic architecture of CAD across all human populations.

Genetic variation within the major histocompatibility complex (MHC) contributes substantial risk for systemic lupus erythematosus, but high gene density, extreme polymorphism and extensive linkage disequilibrium (LD) have made fine mapping challenging. To address the problem, we compared two association techniques in two ancestrally diverse populations, African Americans (AAs) and Europeans (EURs). We observed a greater number of Human Leucocyte Antigen (HLA) alleles in AA consistent with the elevated level of recombination in this population. In EUR we observed 50 different A—C—B—DRB1—DQA—DQB multilocus haplotype sequences per hundred individuals; in the AA sample, these multilocus haplotypes were twice as common compared to Europeans. We also observed a strong narrow class II signal in AA as opposed to the long-range LD observed in EUR that includes class I alleles. We performed a Bayesian model choice of the classical HLA alleles and a frequentist analysis that combined both single nucleotide polymorphisms (SNPs) and classical HLA alleles. Both analyses converged on a similar subset of risk HLA alleles: in EUR HLA– B*08:01 + B*18:01 + (DRB1*15:01 frequentist only) + DQA*01:02 + DQB*02:01 + DRB3*02 and in AA HLA–C*17:01 + B*08:01 + DRB1*15:03 + (DQA*01:02 frequentist only) + DQA*02:01 + DQA*05:01+ DQA*05:05 + DQB*03:19 + DQB*02:02. We observed two additional independent SNP associations in both populations: EUR rs146903072 and rs501480; AA rs389883 and rs114118665. The DR2 serotype was best explained by DRB1*15:03 + DQA*01:02 in AA and by DRB1*15:01 + DQA*01:02 in EUR. The DR3 serotype was best explained by DQA*05:01 in AA and by DQB*02:01 in EUR. Despite some differences in underlying HLA allele risk models in EUR and AA, SNP signals across the extended MHC showed remarkable similarity and significant concordance in direction of effect for risk-associated variants.

Explaining the genetics of many diseases is challenging because most associations localize to incompletely characterized regulatory regions. Using new computational methods, we show that transcription factors (TFs) occupy multiple loci associated with individual complex genetic disorders. Application to 213 phenotypes and 1,544 TF binding datasets identified 2,264 relationships between hundreds of TFs and 94 phenotypes, including androgen receptor in prostate cancer and GATA3 in breast cancer. Strikingly, nearly half of systemic lupus erythematosus risk loci are occupied by the Epstein–Barr virus EBNA2 protein and many coclustering human TFs, showing gene–environment interaction. Similar EBNA2-anchored associations exist in multiple sclerosis, rheumatoid arthritis, inflammatory bowel disease, type 1 diabetes, juvenile idiopathic arthritis and celiac disease. Instances of allele-dependent DNA binding and downstream effects on gene expression at plausibly causal variants support genetic mechanisms dependent on EBNA2. Our results nominate mechanisms that operate across risk loci within disease phenotypes, suggesting new models for disease origins.

Objective This study evaluated single nucleotide polymorphisms (SNPs) associated with progression free (PFS) and overall survival (OS) in patients with advanced stage serous EOC. Methods Patients enrolled in GOG-172 and 182 who provided specimens for translational research and consent were included. Germline DNA was evaluated with the Illumina's HumanOMNI1-Quad beadchips and scanned using Illumina's iScan optical imaging system. SNPs with allele frequency>0.05 and genotyping rate>0.98 were included. Analysis of SNPs for PFS and OS was done using Cox regression. Statistical significance was determined using Bonferroni corrected p -values with genomic control adjustment. Results The initial GWAS analysis included 1,124,677 markers in 396 patients. To obtain the final data set, quality control checks were performed and limited to serous tumors and self-identified Caucasian race. In total 636,555 SNPs and 289 patients passed all the filters. The pre-specified statistical level of significance was 7.855e −08 . No SNPs met this criteria for PFS or OS, however, two SNPs were close to significance (rs10899426 p -2.144e − 08 ) (rs6256 p -9.774e −07 ) for PFS and 2 different SNPs were identified (rs295315 p -7.536e − 07 ; rs17693104 p -7.734e − 07 ) which were close to significance for OS. Conclusions Using the pre-specified level of significance of 1×10 −08 , we did not identify any SNPs of statistical significance for OS or PFS, however several were close. The SNP's identified in this GWAS study will require validation and these preliminary findings may lead to identification of novel pathways and biomarkers.

Primary microcephaly is a congenital brain malformation characterized by a head circumference less than three standard deviations below the mean for age and sex and results in moderate to severe mental deficiencies and decreased lifespan. We recently studied two children with primary microcephaly in an otherwise unaffected family. Exome sequencing identified an autosomal recessive mutation leading to an amino acid substitution in a WD40 domain of the highly conserved Coatomer Protein Complex, Subunit Beta 2 (COPB2). To study the role of Copb2 in neural development, we utilized genome editing technology to generate an allelic series in the mouse. Two independent null alleles revealed that Copb2 is essential for early stages of embryogenesis. Mice homozygous for the patient variant (Copb2R254C/R254C) appear to have a grossly normal phenotype, likely due to differences in corticogenesis between the two species. Strikingly, mice heterozygous for the patient mutation and a null allele (Copb2R254C/Znf) show a severe perinatal phenotype including low neonatal weight, significantly increased apoptosis in the brain, and death within the first week of life. Immunostaining of the Copb2R254C/Znf brain revealed a reduction in layer V (CTIP2+) neurons, while the overall cell density of the cortex is unchanged. Moreover, disruption of Copb2 in mouse neurospheres resulted in reduced proliferation. These results identify a general requirement for COPB2 in embryogenesis and a specific role in corticogenesis. We further demonstrate the utility of CRISPR-Cas9 generated mouse models in the study of potential pathogenicity of variants of potential clinical interest.

Aim: Fatty acid amide hydrolase (FAAH) degrades anandamide, an endogenous cannabinoid. We hypothesized that FAAH variants will predict risk of morphine-related adverse outcomes due to opioid–endocannabinoid interactions. Patients & methods: In 101 postsurgical adolescents receiving morphine analgesia, we prospectively studied ventilatory response to 5% CO2(HCVR), respiratory depression (RD) and vomiting. Blood was collected for genotyping and morphine pharmacokinetics. Results: We found significant FAAH–morphine interaction for missense (rs324420) and several regulatory variants, with HCVR (p < 0.0001) and vomiting (p = 0.0339). HCVR was more depressed in patients who developed RD compared with those who did not (p = 0.0034), thus FAAH–HCVR association predicts risk of impending RD from morphine use. Conclusion: FAAH genotypes predict risk for morphine-related adverse outcomes.

Objective IBD is a chronic inflammatory disorder of the GI tract with complex etiology that involves both genetic variants and environmental factors. Several viruses and bacteria have been suggested as potential causes of the disease. The clinical course of IBD appears to be influenced by EBV tissue infection. With a prevalence of >90% in the adult human population, EBV infection is nearly ubiquitous. We reported that the EBV transcription co-factor EBNA2 was concentrated in the risk loci of IBD relative to the remainder of the genome (Pc=1.24E-13, (RR) = 3.33, Nature Genetics 50:699, 2018). The objective of this study was to determine whether new data and studies made available since the data used for our 2018 study further supported the possibility that EBV was related to some IBD as a potential etiologic agent. Methods Genomic approaches offer previously unavailable perspectives for understanding disease mechanisms. The previous analysis was based upon the 9 viral transcription factor and co-factor (TF) datasets available in human cells in 2015.Of the 52 now available there are 29 from EBV. The number of established IBD loci qualifying for analysis has expanded from 112 in 2015 to 324 available now. We applied simulation analysis to assess statistical significance of TF associations with IBD risk loci using RELI (Regulatory Element Locus Intersection) (Nat Genet 50:699, 2018). RELI tests the probability of the observed association by the count of intersections of the same variant from two sources. We extracted 700 SNPs in European ancestry from 93 published GWAS and 350 additional candidate gene studies for IBD. All variants with linkage disequilibrium r2>0.8 of the best variant were included, while loci with r2>0.2 with a more highly significant locus were excluded. Results Among the 52 viral TF datasets, 139 of the 324 IBD loci were occupied by EBNALP (RR=2.08793, Pc=3.7321E-23), 113 by EBNA3C (RR=2.35113, Pc=2.84678E-22), 88 loci by EBNA2 (RR=3.2, Pc=7.8E-32) and 28 loci by EBNA3ABC, which used an antibody that did not distinguish between EBNA3 subtypes. Analysis of 11,535 human TF ChIP-seq datasets produced 1578 with strong associations (p The locus intersections of the three EBV gene products optimally clustered at 150 IBD loci (46%) with a set of 26 human TFs (p The human TFs were enriched for (ChIP-seq) data obtained from EBV-infected B lymphocyte cell line (OR = 10.12425329, Chi-squared = 241.2925605) Moreover, the human TFs that participate in the super-enhancers that form in the latency III EBV infected and transformed B cell are also concentrated among the statistically significant associations (p Conclusion These preliminary findings nominate EBV for a role in the pathogenesis and etiology of IBD by mechanisms operating in transformed B cells through the latency III EBV program of viral gene expression.

Many common illnesses differentially affect men and women for unknown reasons. The autoimmune diseases lupus and Sjögren’s syndrome affect nine times more women than men1,2, whereas schizophrenia affects men more frequently and severely3–5. All three illnesses have their strongest common-genetic associations in the Major Histocompatibility Complex (MHC) locus, an association that in lupus and Sjögren’s syndrome has long been thought to arise from HLA alleles6–13. Here we show that the complement component 4 (C4) genes in the MHC locus, recently found to increase risk for schizophrenia14, generate 7-fold variation in risk for lupus (95% CI: 5.88-8.61; p < 10−117 in total) and 16-fold variation in risk for Sjögren’s syndrome (95% CI: 8.59-30.89; p < 10−23 in total), with C4A protecting more strongly than C4B in both illnesses. The same alleles that increase risk for schizophrenia, greatly reduced risk for lupus and Sjögren’s syndrome. In all three illnesses, C4 alleles acted more strongly in men than in women: common combinations of C4A and C4B generated 14-fold variation in risk for lupus and 31-fold variation in risk for Sjögren’s syndrome in men (vs. 6-fold and 15-fold among women respectively) and affected schizophrenia risk about twice as strongly in men as in women. At a protein level, both C4 and its effector (C3) were present at greater levels in men than women in cerebrospinal fluid (p < 10−5 for both C4 and C3) and plasma among adults ages 20-5015–17, corresponding to the ages of differential disease vulnerability. Sex differences in complement protein levels may help explain the larger effects of C4 alleles in men, women’s greater risk of SLE and Sjögren’s, and men’s greater vulnerability in schizophrenia. These results nominate the complement system as a source of sexual dimorphism in vulnerability to diverse illnesses.

Annals of the rheumatic diseases 76(5), 906-913 (2017). doi:10.1136/annrheumdis-2016-210324, Published by BMJ Publ. Group, London

Objective Systemic lupus erythematosus (SLE) is a chronic autoimmune disease with a strong genetic component. We undertook the present work to perform the first genome-wide association study on individuals from the Americas who are enriched for Native American heritage. Methods We analyzed 3,710 individuals from the US and 4 countries of Latin America who were diagnosed as having SLE, and healthy controls. Samples were genotyped with HumanOmni1 BeadChip. Data on out-of-study controls genotyped with HumanOmni2.5 were also included. Statistical analyses were performed using SNPtest and SNPGWA. Data were adjusted for genomic control and false discovery rate. Imputation was performed using Impute2 and, for classic HLA alleles, HiBag. Odds ratios (ORs) and 95% confidence intervals (95% CIs) were calculated. Results The IRF5–TNPO3 region showed the strongest association and largest OR for SLE (rs10488631: genomic control–adjusted P [Pgcadj] = 2.61 × 10−29, OR 2.12 [95% CI 1.88–2.39]), followed by HLA class II on the DQA2–DQB1 loci (rs9275572: Pgcadj = 1.11 × 10−16, OR 1.62 [95% CI 1.46–1.80] and rs9271366: Pgcadj = 6.46 × 10−12, OR 2.06 [95% CI 1.71–2.50]). Other known SLE loci found to be associated in this population were ITGAM, STAT4, TNIP1, NCF2, and IRAK1. We identified a novel locus on 10q24.33 (rs4917385: Pgcadj = 1.39 × 10−8) with an expression quantitative trait locus (eQTL) effect (Peqtl = 8.0 × 10−37 at USMG5/miR1307), and several new suggestive loci. SLE risk loci previously identified in Europeans and Asians were corroborated. Local ancestry estimation showed that the HLA allele risk contribution is of European ancestral origin. Imputation of HLA alleles suggested that autochthonous Native American haplotypes provide protection against development of SLE. Conclusion Our results demonstrate that studying admixed populations provides new insights in the delineation of the genetic architecture that underlies autoimmune and complex diseases.

Systemic lupus erythematosus (SLE) has a strong but incompletely understood genetic architecture. We conducted an association study with replication in 4,478 SLE cases and 12,656 controls from six East Asian cohorts to identify new SLE susceptibility loci and better localize known loci. We identified ten new loci and confirmed 20 known loci with genome-wide significance. Among the new loci, the most significant locus was GTF2IRD1-GTF2I at 7q11.23 (rs73366469, Pmeta = 3.75 × 10(-117), odds ratio (OR) = 2.38), followed by DEF6, IL12B, TCF7, TERT, CD226, PCNXL3, RASGRP1, SYNGR1 and SIGLEC6. We identified the most likely functional variants at each locus by analyzing epigenetic marks and gene expression data. Ten candidate variants are known to alter gene expression in cis or in trans. Enrichment analysis highlights the importance of these loci in B cell and T cell biology. The new loci, together with previously known loci, increase the explained heritability of SLE to 24%. The new loci share functional and ontological characteristics with previously reported loci and are possible drug targets for SLE therapeutics.

SLE affects millions with increasing diagnostic prevalence. Genetic studies have identified >600 mostly regulatory variant associations. EBV has been nominated as a causal factor for SLE from immunochemistry and epidemiologic studies. EBNA2 concentrates at SLE risk loci (Nat Genet 50:699 2018) suggesting mechanisms. Additional evidence supporting EBV as an SLE origin would advance this hypothesis of etiology. Previous simulation using our RELI algorithm revealed an astonishing association with Epstein-Barr virus nuclear antigen 2 (EBNA2) (OR=5.96, Pc=E-25) (Nat Genet 50:699, 2018). The >100 risk alleles associated at p30, p These new results further the possible validity of the hypothesis that SLE is originally caused by the Latency III expression program of EBV.

Background Explaining the genetics of many diseases is challenging because most associations localize to incompletely understood regulatory regions. Methods We show that transcription factors (TFs) occupy multiple loci of individual complex genetic disorders much more than expected by chance using novel computational methods. Results Application to 213 phenotypes and 1,544 TF binding datasets identifies 2,264 relationships between hundreds of TFs and 94 phenotypes, including AR in prostate cancer and GATA3 in breast cancer. Strikingly, nearly half of the systemic lupus erythematosus risk loci are occupied by the Epstein-Barr virus (EBV) Nuclear Antigen 2 (EBNA2) protein (OR=6, P Conclusions Our results nominate mechanisms that operate across risk loci within disease phenotypes; they suggest new paradigms for disease origin and strongly support a role for Epstein-Barr virus in the generation of systemic lupus erythematosus, as well as of particular other autoimmune diseases, apparently related to lupus by the genomic mechanisms that produce them.

Eosinophilic esophagitis (EoE) is an allergic inflammatory esophageal disorder with a complex underlying genetic etiology often associated with other comorbidities. Using whole-exome sequencing (WES) of 63 patients with EoE and 60 unaffected family members and family-based trio analysis, we sought to uncover rare coding variants. WES analysis identified 5 rare, damaging variants in dehydrogenase E1 and transketolase domain–containing 1 (DHTKD1). Rare variant burden analysis revealed an overabundance of putative, potentially damaging DHTKD1 mutations in EoE (P = 0.01). Interestingly, we also identified 7 variants in the DHTKD1 homolog oxoglutarate dehydrogenase-like (OGDHL). Using shRNA-transduced esophageal epithelial cells and/or patient fibroblasts, we further showed that disruption of normal DHTKD1 or OGDHL expression blunts mitochondrial function. Finally, we demonstrated that the loss of DHTKD1 expression increased ROS production and induced the expression of viperin, a gene previously shown to be involved in production of Th2 cytokines in T cells. Viperin had increased expression in esophageal biopsies of EoE patients compared with control individuals and was upregulated by IL-13 in esophageal epithelial cells. These data identify a series of rare genetic variants implicating DHTKD1 and OGDHL in the genetic etiology of EoE and underscore a potential pathogenic role for mitochondrial dysfunction in EoE.

Background Isocyanates are major causes of occupational asthma, but susceptibility and mechanisms of diisocyanate-induced asthma (DA) remain uncertain. Objective The aim of this study was to identify DA-associated functional genetic variants through next-generation sequencing (NGS), bioinformatics, and functional assays. Methods NGS was performed in 91 workers with DA. Fourteen loci with known DA-associated single nucleotide polymorphisms (SNPs) were sequenced and compared with data from 238 unexposed subjects. Ranking of DA-associated SNPs based on their likelihood to affect gene regulatory mechanisms in the lung yielded 21 prioritized SNPs. Risk and nonrisk oligonucleotides were tested for binding of nuclear extracts from A549, BEAS-2B, and IMR-90 lung cell lines by using electrophoretic mobility shift assays. DNA constructs were cloned into a pGL3 promoter vector for luciferase gene reporter assays. Results NGS detected 130 risk variants associated with DA (3.1 × 10−6 to 6.21 × 10−4), 129 of which were located in noncoding regions. The 21 SNPs prioritized by using functional genomic data sets were in or proximal to 5 genes: cadherin 17 (CDH17; n = 10), activating transcription factor 3 (ATF3; n = 7), family with sequence similarity, member A (FAM71A; n = 2), tachykinin receptor 1 (TACR1; n = 1), and zinc finger and BTB domain-containing protein 16 (ZBTB16; n = 1). Electrophoretic mobility shift assays detected allele-dependent nuclear protein binding in A549 cells for 8 of 21 variants. In the luciferase assay 4 of the 21 SNPs exhibited allele-dependent changes in gene expression. DNA affinity precipitation and mass spectroscopy of rs147978008 revealed allele-dependent binding of H1 histones, which was confirmed by using Western blotting. Conclusions We identified 5 DA-associated potential regulatory SNPs. Four variants exhibited effects on gene regulation (ATF rs11571537, CDH17 rs2446824 and rs2513789, and TACR1 rs2287231). A fifth variant (FAM71A rs147978008) showed nonrisk allele preferential binding to H1 histones. These results demonstrate that many DA-associated genetic variants likely act by modulating gene regulation.

Eosinophilic esophagitis (EoE) is a chronic inflammatory disease of the esophagus triggered by immune hypersensitivity to food. Herein, we tested whether genetic risk factors for known, non-allergic, immune-mediated diseases, particularly those involving autoimmunity, were associated with EoE risk. We used the high-density Immunochip platform, encoding 200,000 genetic variants for major auto-immune disease. Accordingly, 1,214 subjects with EoE of European ancestry and 3,734 population controls were genotyped and assessed using data directly generated or imputed from the previously published GWAS. We found lack of association of EoE with the genetic variants in the major histocompatibility complex (MHC) class I, II, and III genes and nearly all other loci using a highly powered study design with dense genotyping throughout the locus. Importantly, we identified an EoE risk locus at 16p13 with genome-wide significance (Pcombined = 2.05 × 10−9, odds ratio = 0.76–0.81). This region is known to encode for the genes CLEC16A, DEXI, and CIITI, which are expressed in immune cells and esophageal epithelial cells. Suggestive EoE risk were also seen 5q23 (intergenic) and 7p15 (JAZF1). Overall, we have identified an additional EoE risk locus at 16p13 and highlight a shared and unique genetic etiology of EoE with a spectrum of immune-associated diseases.

Systemic lupus erythematosus (SLE or lupus) (OMIM: 152700) is a chronic autoimmune disease with debilitating inflammation that affects multiple organ systems. The STAT1-STAT4 locus is one of the first and most highly replicated genetic loci associated with lupus risk. We performed a fine-mapping study to identify plausible causal variants within the STAT1-STAT4 locus associated with increased lupus disease risk. Using complementary frequentist and Bayesian approaches in trans-ancestral Discovery and Replication cohorts, we found one variant whose association with lupus risk is supported across ancestries in both the Discovery and Replication cohorts: rs11889341. In B cell lines from patients with lupus and healthy controls, the lupus risk allele of rs11889341 was associated with increased STAT1 expression. We demonstrated that the transcription factor HMGA1, a member of the HMG transcription factor family with an AT-hook DNA-binding domain, has enriched binding to the risk allele compared with the non-risk allele of rs11889341.We identified a genotype-dependent repressive element in the DNA within the intron of STAT4 surrounding rs11889341. Consistent with expression quantitative trait locus (eQTL) analysis, the lupus risk allele of rs11889341 decreased the activity of this putative repressor. Altogether, we present a plausible molecular mechanism for increased lupus risk at the STAT1-STAT4 locus in which the risk allele of rs11889341, the most probable causal variant, leads to elevated STAT1 expression in B cells due to decreased repressor activity mediated by increased binding of HMGA1. © The Author(s) 2018. Published by Oxford University Press. All rights reserved.

Leptin is abnormally elevated in the plasma of patients with systemic lupus erythematosus (SLE), where it is thought to promote and/or sustain pro-inflammatory responses. Whether this association could reflect an increased genetic susceptibility to develop SLE is not known, and studies of genetic associations with leptin-related polymorphisms in SLE patients have been so far inconclusive. Here we genotyped DNA samples from 15,706 SLE patients and healthy matched controls from four different ancestral groups, to correlate polymorphisms of genes of the leptin pathway to risk for SLE. It was found that although several SNPs showed weak associations, those associations did not remain significant after correction for multiple testing. These data do not support associations between defined leptin-related polymorphisms and increased susceptibility to develop SLE. © 2015 Elsevier Inc.

Diisocyanates, reactive chemicals used to produce polyurethane products, are the most common causes of occupational asthma. The aim of this study is to identify susceptibility gene variants that could contribute to the pathogenesis of diisocyanate asthma (DA) using a Genome-Wide Association Study (GWAS) approach. Genome-wide single nucleotide polymorphism (SNP) genotyping was performed in 74 diisocyanate-exposed workers with DA and 824 healthy controls using Omni-2.5 and Omni-5 SNP microarrays. We identified 11 SNPs that exceeded genome-wide significance; the strongest association was for the rs12913832 SNP located on chromosome 15, which has been mapped to the HERC2 gene (p = 6.94 × 10(-14)). Strong associations were also found for SNPs near the ODZ3 and CDH17 genes on chromosomes 4 and 8 (rs908084, p = 8.59 × 10(-9) and rs2514805, p = 1.22 × 10(-8), respectively). We also prioritized 38 SNPs with suggestive genome-wide significance (p < 1 × 10(-6)). Among them, 17 SNPs map to the PITPNC1, ACMSD, ZBTB16, ODZ3, and CDH17 gene loci. Functional genomics data indicate that 2 of the suggestive SNPs (rs2446823 and rs2446824) are located within putative binding sites for the CCAAT/Enhancer Binding Protein (CEBP) and Hepatocyte Nuclear Factor 4, Alpha transcription factors (TFs), respectively. This study identified SNPs mapping to the HERC2, CDH17, and ODZ3 genes as potential susceptibility loci for DA. Pathway analysis indicated that these genes are associated with antigen processing and presentation, and other immune pathways. Overlap of 2 suggestive SNPs with likely TF binding sites suggests possible roles in disruption of gene regulation. These results provide new insights into the genetic architecture of DA and serve as a basis for future functional and mechanistic studies.

Complement receptor 2 (CR2/CD21) is predominantly expressed on the surface of mature B cells where it forms part of a coreceptor complex that functions, in part, to modulate B-cell receptor signal strength. CR2/CD21 expression is tightly regulated throughout B-cell development such that CR2/CD21 cannot be detected on pre-B or terminally differentiated plasma cells. CR2/CD21 expression is upregulated at B-cell maturation and can be induced by IL-4 and CD40 signaling pathways. We have previously characterized elements in the proximal promoter and first intron of CR2/CD21 that are involved in regulating basal and tissue-specific expression. We now extend these analyses to the CR2/CD21 core promoter. We show that in mature B cells, CR2/CD21 transcription proceeds from a focused TSS regulated by a non-consensus TATA box, an initiator element and a downstream promoter element. Furthermore, occupancy of the general transcriptional machinery in pre-B versus mature B-cell lines correlate with CR2/CD21 expression level and indicate that promoter accessibility must switch from inactive to active during the transitional B-cell window.

Opioid effects are potentiated by cannabinoid agonists including anandamide, an endocannabinoid. Inter-individual variability in responses to opioids is a major clinical problem. Multiple deaths and anoxic brain injuries occur every year because of opioid-induced respiratory depression (RD) in surgical patients and drug abusers of opioids and cannabinoids. This study aimed to determine specific associations between genetic variants of fatty acid amide hydrolase (FAAH) and postoperative central opioid adverse effects in children undergoing tonsillectomy. This is a prospective genotype-blinded observational study in which 259 healthy children between 6 and 15 years of age who received standard perioperative care with a standard anesthetic and an intraoperative dose of morphine were enrolled. Associations between frequent polymorphisms of FAAH and central postoperative opioid adverse effects including, RD, postoperative nausea and vomiting (PONV) and prolonged stay in Post Anesthesia Recovery Room (postoperative anesthesia care unit, PACU) due to RD and PONV were analyzed. Five specific FAAH single nucleotide polymorphisms (SNPs) had significant associations with more than twofold increased risk for refractory PONV (adjusted P

Sjögren’s syndrome (SS) is a common, autoimmune exocrinopathy distinguished by keratoconjunctivitis sicca and xerostomia. Patients frequently develop serious complications including lymphoma, pulmonary dysfunction, neuropathy, vasculitis, and debilitating fatigue. Dysregulation of type I interferon (IFN) pathway is a prominent feature of SS and is correlated with increased autoantibody titers and disease severity. To identify genetic determinants of IFN pathway dysregulation in SS, we performed cis-expression quantitative trait locus (eQTL) analyses focusing on differentially expressed type I IFN-inducible transcripts identified through a transcriptome profiling study. Multiple cis-eQTLs were associated with transcript levels of 2'-5'-oligoadenylate synthetase 1 (OAS1) peaking at rs10774671 (PeQTL = 6.05 × 10−14). Association of rs10774671 with SS susceptibility was identified and confirmed through meta-analysis of two independent cohorts (Pmeta = 2.59 × 10−9; odds ratio = 0.75; 95% confidence interval = 0.66–0.86). The risk allele of rs10774671 shifts splicing of OAS1 from production of the p46 isoform to multiple alternative transcripts, including p42, p48, and p44. We found that the isoforms were differentially expressed within each genotype in controls and patients with and without autoantibodies. Furthermore, our results showed that the three alternatively spliced isoforms lacked translational response to type I IFN stimulation. The p48 and p44 isoforms also had impaired protein expression governed by the 3' end of the transcripts. The SS risk allele of rs10774671 has been shown by others to be associated with reduced OAS1 enzymatic activity and ability to clear viral infections, as well as reduced responsiveness to IFN treatment. Our results establish OAS1 as a risk locus for SS and support a potential role for defective viral clearance due to altered IFN response as a genetic pathophysiological basis of this complex autoimmune disease., Author summary Sjögren’s syndrome (SS) is a common autoimmune condition where immune cells infiltrate moisture-producing glands, leading to dryness typically in the eyes and mouth. SS patients also manifest debilitating fatigue as well as other diseases in liver, lung, kidney, and skin. The cause of this complex disease is still not fully understood; however, an environmental trigger, such as viral infections, in individuals with genetic risk factor(s) is thought to contribute to the development of SS. Type 1 interferons (IFNs) are one of the first defenders after viral infection and induce the expression of various virus-responding genes. Perpetual elevation of type 1 IFN signaling has been observed in SS patients. Here, we first replicated previously identified RNA transcripts contributing to the abnormal type 1 IFN signaling in SS patients. We then identified a disease-associated genetic variant in an IFN-inducible gene, OAS1. This variant governs splicing of OAS1, altering the transcript into multiple isoforms that lack protein expression and responsiveness to IFNs. The results of this study may provide insight into the genetic basis of SS, as well as other autoimmune disease with similar dysregulation in the type 1 IFN system.

Systemic lupus erythematosus (SLE) is an autoimmune disease with marked gender and ethnic disparities. We report a large transancestral association study of SLE using Immunochip genotype data from 27,574 individuals of European (EA), African (AA) and Hispanic Amerindian (HA) ancestry. We identify 58 distinct non-HLA regions in EA, 9 in AA and 16 in HA (∼50% of these regions have multiple independent associations); these include 24 novel SLE regions (P, Systemic lupus erythematosus (SLE) is an autoimmune disease with a strong ethnic and gender bias. In a transancestral genetic association study, Langefeld et al. identify 24 novel regions associated with risk to lupus and propose a cumulative hits hypothesis for loci conferring risk to SLE.

Objective Macrophage activation syndrome (MAS), a life-threatening complication of systemic juvenile idiopathic arthritis (JIA), resembles familial hemophagocytic lymphohistiocytosis (HLH), a constellation of autosomal-recessive immune disorders resulting from deficiency in cytolytic pathway proteins. We undertook this study to test our hypothesis that MAS predisposition in systemic JIA could be attributed to rare gene sequence variants affecting the cytotolytic pathway. Methods Whole-exome sequencing was used in 14 patients with systemic JIA and MAS and in their parents to identify protein-altering single-nucleotide polymorphisms/indels in known HLH-associated genes. To discover new candidate genes, the entire whole-exome sequencing data were filtered to identify protein-altering, rare recessive homozygous, compound heterozygous, and de novo variants with the potential to affect the cytolytic pathway. Results Heterozygous protein-altering rare variants in the known genes (LYST,MUNC13-4, and STXBP2) were found in 5 of 14 patients with systemic JIA and MAS (35.7%). This was in contrast to only 4 variants in 4 of 29 patients with systemic JIA without MAS (13.8%). Homozygosity and compound heterozygosity analysis applied to the entire whole-exome sequencing data in systemic JIA/MAS revealed 3 recessive pairs in 3 genes and compound heterozygotes in 73 genes. We also identified 20 heterozygous rare protein-altering variants that occurred in at least 2 patients. Many of the identified genes encoded proteins with a role in actin and microtubule reorganization and vesicle-mediated transport. "Cellular assembly and organization" was the top cellular function category based on Ingenuity Pathways Analysis (P < 3.10 × 10-5). Conclusion Whole-exome sequencing performed in patients with systemic JIA and MAS identified rare protein-altering variants in known HLH-associated genes as well as in new candidate genes.

Exploiting genotyping, DNA sequencing, imputation and trans-ancestral mapping, we used Bayesian and frequentist approaches to model the IRF5-TNPO3 locus association, now implicated in two immunotherapies and seven autoimmune diseases. Specifically, in systemic lupus erythematosus (SLE), we resolved separate associations in the IRF5 promoter (all ancestries) and with an extended European haplotype. We captured 3230 IRF5-TNPO3 high-quality, common variants across 5 ethnicities in 8395 SLE cases and 7367 controls. The genetic effect from the IRF5 promoter can be explained by any one of four variants in 5.7 kb (P-valuemeta = 6 × 10-49; OR = 1.38-1.97). The second genetic effect spanned an 85.5-kb, 24-variant haplotype that included the genes IRF5 and TNPO3(P-valuesEU = 10-27-10-32, OR = 1.7-1.81). Many variants at the IRF5 locus with previously assigned biological function are not members of either final credibleset of potential causal variants identified herein. In addition to the known biologically functional variants, we demonstrated that the risk allele of rs4728142, a variant in the promoter among the lowest frequentist probability and highest Bayesian posterior probability, was correlated with IRF5 expression and differentially binds the transcription factor ZBTB3. Our analytical strategy provides a novel framework for future studies aimed at dissecting etiological genetic effects. Finally, both SLE elements of the statistical model appear to operate in Sjögren's syndrome and systemic sclerosis whereas only the IRF5-TNPO3 gene-spanning haplotype is associated with primary biliary cirrhosis, demonstrating the nuance of similarity and difference in autoimmune disease risk mechanisms at IRF5-TNPO3. © The Author 2014.

The genetic factors underlying the pathogenesis of lupus nephritis associated with systemic lupus erythematosus are largely unknown, although animal studies indicate that nuclear factor (NF)-?B is involved. We reported previously that aknockin mouse expressinganin active form of ABIN1 (ABIN1[D485N]) develops lupus-like autoimmune disease and demonstrates enhanced activation of NF-?B and mitogen-activated protein kinases in immune cells after toll-like receptor stimulation. In the current study, we show that ABIN1[D485N] mice develop progressive GN similar to class III and IV lupus nephritis in humans. To investigate the clinical relevance of ABIN1 dysfunction, we genotyped five single-nucleotide polymorphisms in the gene encoding ABIN1, TNIP1, in samples from European-American, African American, Asian, Gullah, and Hispanic participants in the Large Lupus Association Study 2. Comparing cases of systemic lupus erythematosus with nephritis and cases ofsystemic lupus erythematosus without nephritis revealed strong associations with lupus nephritis at rs7708392 in European Americans and rs4958881 in African Americans. Comparing cases of systemic lupus erythematosus with nephritis and healthy controls revealed a stronger association at rs7708392 in European Americans but not at rs4958881 in African Americans. Our data suggest that variants in the TNIP1 gene are associated with the risk for lupus nephritis and could be mechanistically involved in disease development via aberrant regulation of NF-?B and mitogen-activated protein kinase activity. Copyright © 2013 by the American Society of Nephrology.

Objective.Little is known about the genetic etiology of systemic lupus erythematosus (SLE) in individuals of African ancestry, despite its higher prevalence and greater disease severity. Overproduction of nitric oxide (NO) and reactive oxygen species are implicated in the pathogenesis and severity of SLE, making NO synthases and other reactive intermediate-related genes biological candidates for disease susceptibility. We analyzed variation in reactive intermediate genes for association with SLE in 2 populations with African ancestry.Methods.A total of 244 single-nucleotide polymorphisms (SNP) from 53 regions were analyzed in non-Gullah African Americans (AA; 1432 cases and 1687 controls) and the genetically more homogeneous Gullah of the Sea Islands of South Carolina (133 cases and 112 controls). Single-marker, haplotype, and 2-locus interaction tests were computed for these populations.Results.The glutathione reductase geneGSR(rs2253409; p = 0.0014, OR 1.26, 95% CI 1.09–1.44) was the most significant single SNP association in AA. In the Gullah, the NADH dehydrogenaseNDUFS4(rs381575; p = 0.0065, OR 2.10, 95% CI 1.23–3.59) and NO synthase geneNOS1(rs561712; p = 0.0072, OR 0.62, 95% CI 0.44–0.88) were most strongly associated with SLE. When both populations were analyzed together,GSRremained the most significant effect (rs2253409; p = 0.00072, OR 1.26, 95% CI 1.10–1.44). Haplotype and 2-locus interaction analyses also uncovered different loci in each population.Conclusion.These results suggest distinct patterns of association with SLE in African-derived populations; specific loci may be more strongly associated within select population groups.

Background A key question in the allergy field is to understand how tissue-specific disease is manifested. Eosinophilic esophagitis (EoE) is an emerging tissue-specific allergic disease with an unclear pathogenesis. Objective Herein we tested the hypothesis that a defect in tissue-specific esophageal genes is an integral part of EoE pathogenesis. Methods We interrogated the pattern of expression of esophagus-specific signature genes derived from the Human Protein Atlas in the EoE transcriptome and in EPC2 esophageal epithelial cells. Western blotting and immunofluorescence were used for evaluating expression of esophageal proteins in biopsy specimens from control subjects and patients with active EoE. Whole-exome sequencing was performed to identify mutations in esophagus-specific genes. Results We found that approximately 39% of the esophagus-specific transcripts were altered in patients with EoE, with approximately 90% being downregulated. The majority of transcriptional changes observed in esophagus-specific genes were reproduced in vitro in esophageal epithelial cells differentiated in the presence of IL-13. Functional enrichment analysis revealed keratinization and differentiation as the most affected biological processes and identified IL-1 cytokines and serine peptidase inhibitors as the most dysregulated esophagus-specific protein families in patients with EoE. Accordingly, biopsy specimens from patients with EoE evidenced a profound loss of tissue differentiation, decreased expression of keratin 4 (KRT4) and cornulin (CRNN) , and increased expression of KRT5 and KRT14 . Whole-exome sequencing of 33 unrelated patients with EoE revealed 39 rare mutations in 18 esophagus-specific differentially expressed genes. Conclusions A tissue-centered analysis has revealed a profound loss of esophageal tissue differentiation (identity) as an integral and specific part of the pathophysiology of EoE and implicated protease- and IL-1–related activities as putative central pathways in disease pathogenesis.

The widespread availability of next-generation sequencing (NGS) has transformed our understanding of human genetic variation and its impact on human health. This chapter describes the most common DNA sequencing technologies available to research and clinical laboratories today, and resources for interpreting the functional impact of genetic variants identified with these technologies. Targeted genetic capture techniques were developed to dramatically decrease the cost of determining variant genotypes, although as prices decrease many of the advantages of targeted experiments diminish. Targeted whole exome sequencing is used to identify variants that alter the amino acid sequence of a protein. Sequencing can also be performed for the whole genome, enabling the identification of variants that fall within non-coding regions, which might alter the expression of a gene by disrupting regulatory sequences such as transcription factor binding sites. Careful consideration of the study design, particularly the use of prior knowledge about the phenotype in question, family history, and the availability of affected and unaffected family members, increases the chances that meaningful results are obtained. Identifying variants in the short, error-prone sequencing reads generated by modern technologies is challenging, although software packages exist that mitigate the most common types of errors. The most difficult task in analyzing results is interpreting the functional impact of putative variants and differentiating between clinically reportable variants and variants of unknown significance (VUS) or variants within genes of unknown significance (GUS). Additional information, such as population allele frequencies, genetic inheritance patterns, and functional genomics data, can help to identify the variants most likely to be involved in disease pathogenesis.

ObjectivesFollowing up the systemic lupus erythematosus (SLE) genome-wide association studies (GWAS) identification of NMNAT2 at rs2022013, we fine-mapped its 150 kb flanking regions containing NMNAT2 and SMG7 in a 15 292 case–control multi-ancestry population and tested functions of identified variants.MethodsWe performed genotyping using custom array, imputation by IMPUTE 2.1.2 and allele specific functions using quantitative real-time PCR and luciferase reporter transfections. SLE peripheral blood mononuclear cells (PBMCs) were cultured with small interfering RNAs to measure antinuclear antibody (ANA) and cyto/chemokine levels in supernatants using ELISA.ResultsWe confirmed association at NMNAT2 in European American (EA) and Amerindian/Hispanic ancestries, and identified independent signal at SMG7 tagged by rs2702178 in EA only (p=2.4×10−8, OR=1.23 (95% CI 1.14 to 1.32)). In complete linkage disequilibrium with rs2702178, rs2275675 in the promoter region robustly associated with SMG7 mRNA levels in multiple expression quantitative trait locus (eQTL) datasets. Its risk allele was dose-dependently associated with decreased SMG7 mRNA levels in PBMCs of 86 patients with SLE and 119 controls (p=1.1×10−3 and 6.8×10−8, respectively) and conferred reduced transcription activity in transfected HEK-293 (human embryonic kidney cell line) and Raji cells (p=0.0035 and 0.0037, respectively). As a critical component in the nonsense-mediated mRNA decay pathway, SMG7 could regulate autoantigens including ribonucleoprotein (RNP) and Smith (Sm). We showed SMG7 mRNA levels in PBMCs correlated inversely with ANA titres of patients with SLE (r=−0.31, p=0.01), and SMG7 knockdown increased levels of ANA IgG and chemokine (C-C motif) ligand 19 in SLE PBMCs (p=2.0×10−5 and 2.0×10−4, respectively).ConclusionWe confirmed NMNAT2 and identified independent SMG7 association with SLE. The inverse relationship between levels of the risk allele-associated SMG7 mRNAs and ANA suggested the novel contribution of mRNA surveillance pathway to SLE pathogenesis.

Objective American Indian-Europeans, Asians, and African Americans have an excess morbidity from systemic lupus erythematosus (SLE) and a higher prevalence of lupus nephritis than do Caucasians. The aim of this study was to analyze the relationship between genetic ancestry and sociodemographic characteristics and clinical features in a large cohort of American Indian-European SLE patients. Methods A total of 2,116 SLE patients of American Indian-European origin and 4,001 SLE patients of European descent for whom we had clinical data were included in the study. Genotyping of 253 continental ancestry-informative markers was performed on the Illumina platform. Structure and Admixture software were used to determine genetic ancestry proportions of each individual. Logistic regression was used to test the association between genetic ancestry and sociodemographic and clinical characteristics. Odds ratios (ORs) were calculated with 95% confidence intervals (95% CIs). Results The average American Indian genetic ancestry of 2,116 SLE patients was 40.7%. American Indian genetic ancestry conferred increased risks of renal involvement (P less than 0.0001, OR 3.50 [95% CI 2.63- 4.63]) and early age at onset (P less than 0.0001). American Indian ancestry protected against photosensitivity (P less than 0.0001, OR 0.58 [95% CI 0.44-0.76]), oral ulcers (P less than 0.0001, OR 0.55 [95% CI 0.42-0.72]), and serositis (P less than 0.0001, OR 0.56 [95% CI 0.41-0.75]) after adjustment for age, sex, and age at onset. However, age and sex had stronger effects than genetic ancestry on malar rash, discoid rash, arthritis, and neurologic involvement. Conclusion In general, American Indian genetic ancestry correlates with lower sociodemographic status and increases the risk of developing renal involvement and SLE at an earlier age. Copyright © 2012 by the American College of Rheumatology.

Objectives: Systemic lupus erythematosus (SLE) is a sexually dimorphic autoimmune disease which is more common in women, but affected men often experience a more severe disease. The genetic basis of sexual dimorphism in SLE is not clearly defined. A study was undertaken to examine sex-specific genetic effects among SLE susceptibility loci. Methods: A total of 18 autosomal genetic susceptibility loci for SLE were genotyped in a large set of patients with SLE and controls of European descent, consisting of 5932 female and 1495 male samples. Sex-specific genetic association analyses were performed. The sex-gene interaction was further validated using parametric and non-parametric methods. Aggregate differences in sex-specific genetic risk were examined by calculating a cumulative genetic risk score for SLE in each individual and comparing the average genetic risk between male and female patients. Results: A significantly higher cumulative genetic risk for SLE was observed in men than in women. (P=4.52x10 -8) A significant sex-gene interaction was seen primarily in the human leucocyte antigen (HLA) region but also in IRF5, whereby men with SLE possess a significantly higher frequency of risk alleles than women. The genetic effect observed in KIAA1542 is specific to women with SLE and does not seem to have a role in men. Conclusions: The data indicate that men require a higher cumulative genetic load than women to develop SLE. These observations suggest that sex bias in autoimmunity could be influenced by autosomal genetic susceptibility loci.