Human cytomegalovirus extensively re-organizes the human genome, diminishing TEAD1 transcription factor activity

Human cytomegalovirus (HCMV) infects up to 80% of the world's population. Here, we show that HCMV infection leads to widespread changes in human chromatin accessibility and chromatin looping, with hundreds of thousands of genomic regions affected 48 hours after infection. Integrative analyses reveal HCMV-induced perturbation of Hippo signaling through drastic reduction of TEAD1 transcription factor activity. We confirm extensive concordant loss of TEAD1 binding, active H3K27ac histone marks, and chromatin looping interactions upon infection. Our data position TEAD1 at the top of a hierarchy involving multiple altered important developmental pathways. HCMV infection reduces TEAD1 activity through four distinct mechanisms: closing of TEAD1-bound chromatin, reduction of YAP1 and phosphorylated YAP1 levels, reduction of TEAD1 transcript and protein levels, and alteration of TEAD1 exon-6 usage. Altered TEAD1-based mechanisms are highly enriched at genetic risk loci associated with eye and ear development, providing mechanistic insight into HCMV's established roles in these processes.

Association between systemic lupus erythematosus and myasthenia gravis: A population-based National Study

BACKGROUND

Systemic lupus erythematosus (SLE) and myasthenia gravis (MG) are autoimmune diseases. Previous case reports and case series suggest an association may exist between these diseases, as well as an increased risk of SLE after thymectomy for MG. We undertook this study to determine whether SLE and MG were associated in large cohorts.

METHODS

We searched the IBM Watson Health Explorys platform and the Department of Veterans Affairs Million Veteran Program (MVP) database for diagnoses of SLE and MG. In addition, we examined subjects enrolled in the Lupus Family Registry and Repository (LFRR) as well as controls for a diagnosis of MG.

RESULTS

Among 59,780,210 individuals captured in Explorys, there were 25,750 with MG and 65,370 with SLE. 370 subjects had both. Those with MG were >10 times more likely to have SLE than those without MG. Those with both diseases were more likely to be women, African American, and at a younger age than MG subjects without SLE. In addition, the MG patients who underwent thymectomy had an increased risk of SLE compared to MG patients who had not undergone thymectomy (OR 3.11, 95% CI: 2.12 to 4.55). Autoimmune diseases such as pernicious anemia and miscellaneous comorbidities such as chronic kidney disease were significantly more common in MG patients who developed SLE. In the MVP, SLE and MG were also significantly associated. Association of SLE and MG in a large SLE cohort with rigorous SLE classification confirmed the association of SLE with MG at a similar level.

CONCLUSION

While the number of patients with both MG and SLE is small, SLE and MG are strongly associated together in very large databases and a large SLE cohort.

Modeling of horizontal pleiotropy identifies possible causal gene expression in systemic lupus erythematosus

Background

Systemic lupus erythematosus (SLE) is a chronic autoimmune condition with complex causes involving genetic and environmental factors. While genome-wide association studies (GWASs) have identified genetic loci associated with SLE, the functional genomic elements responsible for disease development remain largely unknown. Mendelian Randomization (MR) is an instrumental variable approach to causal inference based on data from observational studies, where genetic variants are employed as instrumental variables (IVs).

Methods

This study utilized a two-step strategy to identify causal genes for SLE. In the first step, the classical MR method was employed, assuming the absence of horizontal pleiotropy, to estimate the causal effect of gene expression on SLE. In the second step, advanced probabilistic MR methods (PMR-Egger, MRAID, and MR-MtRobin) were applied to the genes identified in the first step, considering horizontal pleiotropy, to filter out false positives. PMR-Egger and MRAID analyses utilized whole blood expression quantitative trait loci (eQTL) and SLE GWAS summary data, while MR-MtRobin analysis used an independent eQTL dataset from multiple immune cell types along with the same SLE GWAS data.

Results

The initial MR analysis identified 142 genes, including 43 outside of chromosome 6. Subsequently, applying the advanced MR methods reduced the number of genes with significant causal effects on SLE to 66. PMR-Egger, MRAID, and MR-MtRobin, respectively, identified 13, 7, and 16 non-chromosome 6 genes with significant causal effects. All methods identified expression of PHRF1 gene as causal for SLE. A comprehensive literature review was conducted to enhance understanding of the functional roles and mechanisms of the identified genes in SLE development.

Conclusions

The findings from the three MR methods exhibited overlapping genes with causal effects on SLE, demonstrating consistent results. However, each method also uncovered unique genes due to different modelling assumptions and technical factors, highlighting the complementary nature of the approaches. Importantly, MRAID demonstrated a reduced percentage of causal genes from the Major Histocompatibility complex (MHC) region on chromosome 6, indicating its potential in minimizing false positive findings. This study contributes to unraveling the mechanisms underlying SLE by employing advanced probabilistic MR methods to identify causal genes, thereby enhancing our understanding of SLE pathogenesis.

Genome-wide association analyses define pathogenic signaling pathways and prioritize drug targets for IgA nephropathy

IgA nephropathy (IgAN) is a progressive form of kidney disease defined by glomerular deposition of IgA. Here we performed a genome-wide association study of 10,146 kidney-biopsy-diagnosed IgAN cases and 28,751 controls across 17 international cohorts. We defined 30 genome-wide significant risk loci explaining 11% of disease risk. A total of 16 loci were new, including TNFSF4/TNFSF18, REL, CD28, PF4V1, LY86, LYN, ANXA3, TNFSF8/TNFSF15, REEP3, ZMIZ1, OVOL1/RELA, ETS1, IGH, IRF8, TNFRSF13B and FCAR. The risk loci were enriched in gene orthologs causing abnormal IgA levels when genetically manipulated in mice. We also observed a positive genetic correlation between IgAN and serum IgA levels. High polygenic score for IgAN was associated with earlier onset of kidney failure. In a comprehensive functional annotation analysis of candidate causal genes, we observed convergence of biological candidates on a common set of inflammatory signaling pathways and cytokine ligand-receptor pairs, prioritizing potential new drug targets.

The US Department of Veterans Affairs Science and Health Initiative to Combat Infectious and Emerging Life-Threatening Diseases (VA SHIELD): A Biorepository Addressing National Health Threats

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.

Genetic regulation of serum IgA levels and susceptibility to common immune, infectious, kidney, and cardio-metabolic traits

Immunoglobulin A (IgA) mediates mucosal responses to food antigens and the intestinal microbiome and is involved in susceptibility to mucosal pathogens, celiac disease, inflammatory bowel disease, and IgA nephropathy. We performed a genome-wide association study of serum IgA levels in 41,263 individuals of diverse ancestries and identified 20 genome-wide significant loci, including 9 known and 11 novel loci. Co-localization analyses with expression QTLs prioritized candidate genes for 14 of 20 significant loci. Most loci encoded genes that produced immune defects and IgA abnormalities when genetically manipulated in mice. We also observed positive genetic correlations of serum IgA levels with IgA nephropathy, type 2 diabetes, and body mass index, and negative correlations with celiac disease, inflammatory bowel disease, and several infections. Mendelian randomization supported elevated serum IgA as a causal factor in IgA nephropathy. African ancestry was consistently associated with higher serum IgA levels and greater frequency of IgA-increasing alleles compared to other ancestries. Our findings provide novel insights into the genetic regulation of IgA levels and its potential role in human disease.

Haplotype-specific chromatin looping reveals genetic interactions of regulatory regions modulating gene expression in 8p23.1

A major goal of genetics research is to elucidate mechanisms explaining how genetic variation contributes to phenotypic variation. The genetic variants identified in genome-wide association studies (GWASs) generally explain only a small proportion of heritability of phenotypic traits, the so-called missing heritability problem. Recent evidence suggests that additional common variants beyond lead GWAS variants contribute to phenotypic variation; however, their mechanistic underpinnings generally remain unexplored. Herein, we undertake a study of haplotype-specific mechanisms of gene regulation at 8p23.1 in the human genome, a region associated with a number of complex diseases. The FAM167A-BLK locus in this region has been consistently found in the genome-wide association studies (GWASs) of systemic lupus erythematosus (SLE) in all major ancestries. Our haplotype-specific chromatin interaction (Hi-C) experiments, allele-specific enhancer activity measurements, genetic analyses, and epigenome editing experiments revealed that: 1) haplotype-specific long-range chromatin interactions are prevalent in 8p23.1; 2) BLK promoter and cis-regulatory elements cooperatively interact with haplotype-specificity; 3) genetic variants at distal regulatory elements are allele-specific modifiers of the promoter variants at FAM167A-BLK; 4) the BLK promoter interacts with and, as an enhancer-like promoter, regulates FAM167A expression and 5) local allele-specific enhancer activities are influenced by global haplotype structure due to chromatin looping. Although systemic lupus erythematosus causal variants at the FAM167A-BLK locus are thought to reside in the BLK promoter region, our results reveal that genetic variants at distal regulatory elements modulate promoter activity, changing BLK and FAM167A gene expression and disease risk. Our results suggest that global haplotype-specific 3-dimensional chromatin looping architecture has a strong influence on local allelic BLK and FAM167A gene expression, providing mechanistic details for how regional variants controlling the BLK promoter may influence disease risk.

Regulatory Architecture of the RCA Gene Cluster Captures an Intragenic TAD Boundary, CTCF-Mediated Chromatin Looping and a Long-Range Intergenic Enhancer

The Regulators of Complement Activation (RCA) gene cluster comprises several tandemly arranged genes with shared functions within the immune system. RCA members, such as complement receptor 2 (CR2), are well-established susceptibility genes in complex autoimmune diseases. Altered expression of RCA genes has been demonstrated at both the functional and genetic level, but the mechanisms underlying their regulation are not fully characterised. We aimed to investigate the structural organisation of the RCA gene cluster to identify key regulatory elements that influence the expression of CR2 and other genes in this immunomodulatory region. Using 4C, we captured extensive CTCF-mediated chromatin looping across the RCA gene cluster in B cells and showed these were organised into two topologically associated domains (TADs). Interestingly, an inter-TAD boundary was located within the CR1 gene at a well-characterised segmental duplication. Additionally, we mapped numerous gene-gene and gene-enhancer interactions across the region, revealing extensive co-regulation. Importantly, we identified an intergenic enhancer and functionally demonstrated this element upregulates two RCA members (CR2 and CD55) in B cells. We have uncovered novel, long-range mechanisms whereby autoimmune disease susceptibility may be influenced by genetic variants, thus highlighting the important contribution of chromatin topology to gene regulation and complex genetic disease.

Pre-Clinical Autoimmunity in Lupus Relatives: Self-Reported Questionnaires and Immune Dysregulation Distinguish Relatives Who Develop Incomplete or Classified Lupus From Clinically Unaffected Relatives and Unaffected, Unrelated Individuals

Systemic lupus erythematosus (SLE) is propelled by pathogenic autoantibody (AutoAb) and immune pathway dysregulation. Identifying populations at risk of reaching classified SLE is essential to curtail inflammatory damage. Lupus blood relatives (Rel) have an increased risk of developing SLE. We tested factors to identify Rel at risk of developing incomplete lupus (ILE) or classified SLE vs. clinically unaffected Rel and healthy controls (HC), drawing from two unique, well characterized lupus cohorts, the lupus autoimmunity in relatives (LAUREL) follow-up cohort, consisting of Rel meeting <4 ACR criteria at baseline, and the Lupus Family Registry and Repository (LFRR), made up of SLE patients, lupus Rel, and HC. Medical record review determined ACR SLE classification criteria; study participants completed the SLE portion of the connective tissue disease questionnaire (SLE-CSQ), type 2 symptom questions, and provided samples for assessment of serum SLE-associated AutoAb specificities and 52 plasma immune mediators. Elevated SLE-CSQ scores were associated with type 2 symptoms, ACR scores, and serology in both cohorts. Fatigue at BL was associated with transition to classified SLE in the LAUREL cohort (p≤0.01). Increased levels of BLyS and decreased levels of IL-10 were associated with type 2 symptoms (p<0.05). SLE-CSQ scores, ACR scores, and accumulated AutoAb specificities correlated with levels of multiple inflammatory immune mediators (p<0.05), including BLyS, IL-2Rα, stem cell factor (SCF), soluble TNF receptors, and Th-1 type mediators and chemokines. Transition to SLE was associated with increased levels of SCF (p<0.05). ILE Rel also had increased levels of TNF-α and IFN-γ, offset by increased levels of regulatory IL-10 and TGF-β (p<0.05). Clinically unaffected Rel (vs. HC) had higher SLE-CSQ scores (p<0.001), increased serology (p<0.05), and increased inflammatory mediator levels, offset by increased IL-10 and TGF-β (p<0.01). These findings suggest that Rel at highest risk of transitioning to classified SLE have increased inflammation coupled with decreased regulatory mediators. In contrast, clinically unaffected Rel and Rel with ILE demonstrate increased inflammation offset with increased immune regulation, intimating a window of opportunity for early intervention and enrollment in prevention trials.

A High Prevalence of Anti-EBNA1 Heteroantibodies in Systemic Lupus Erythematosus (SLE) Supports Anti-EBNA1 as an Origin for SLE Autoantibodies

Background

That Epstein–Barr virus (EBV) infection is associated with systemic lupus erythematosus (SLE) is established. The challenge is to explain mechanistic roles EBV has in SLE pathogenesis. Previous studies identify four examples of autoantibody cross-reactions between SLE autoantigens and Epstein–Barr nuclear antigen 1 (EBNA1). For two of these examples, the earliest detected autoantibody specifically cross-reacts with EBNA1; thereby, defined EBNA1 epitopes induce a robust autoantibody response in animals. These results suggest that the autoantibodies initiating the process leading to SLE may emerge from the anti-EBNA1 heteroimmune response. If this hypothesis is true, then anti-EBNA1 responses would be more frequent in EBV-infected SLE patients than in EBV-infected controls. We tested this prediction.

Methods

We evaluated published East Asian data by selecting those with a positive anti-viral capsid antigen (VCA) antibody immunoglobulin G (IgG) test and determining whether anti-EBNA1 was more common among the EBV-infected SLE cases than among matched EBV-infected controls with conditional logistic regression analysis.

Results

All the qualifying SLE patients (100%) in this dataset were EBV-infected compared to age- and sex-matched controls (92.2%) [odds ratio (OR) = 28.6, 95% CI 6.4–∞, p = 8.83 × 10^-8], confirming the known close association of EBV infection with SLE. Furthermore, virtually all the SLE cases have both anti-VCA IgG and anti-EBNA1 IgG antibodies [124 of 125 (99.2%)], which are more frequently present than in age- and sex-matched EBV-infected controls [232 of 250 (93.2%)] (OR = 9.7, 95% CI 1.5–414, p = 0.0078) for an 89.7% SLE attributable risk from anti-EBNA1, which is in addition to the 100% SLE risk attributable to EBV infection in these data.

Conclusions

The association of EBV infection with SLE is reconfirmed. The prediction that anti-EBNA1 is more frequent in these SLE cases than in EBV-infected controls is true, consistent with the hypothesis that anti-EBNA1 contributes to SLE. This second EBV-dependent risk factor is consistent with a molecular mimicry model for the generation of SLE, starting with EBV infection, progressing to anti-EBNA1 response; then molecular mimicry leads to anti-EBNA1 antibodies cross-reacting with an SLE autoantigen, causing autoantibody epitope spreading, and culminating in clinical SLE. These results support the anti-EBNA1 heteroimmune response being a foundation from which pathogenic SLE autoimmunity emerges.

Lupus enhancer risk variant causes dysregulation of IRF8 through cooperative lncRNA and DNA methylation machinery

Despite strong evidence that human genetic variants affect the expression of many key transcription factors involved in autoimmune diseases, establishing biological links between non-coding risk variants and the gene targets they regulate remains a considerable challenge. Here, we combine genetic, epigenomic, and CRISPR activation approaches to screen for functional variants that regulate IRF8 expression. We demonstrate that the locus containing rs2280381 is a cell-type-specific 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 regulating IRF8 expression by affecting methylation levels. The alleles of rs2280381 modulate PU.1 binding and chromatin state to regulate AC092723.1 and IRF8 expression differentially. Our work illustrates an integrative strategy to define functional genetic variants that regulate the expression of critical genes in autoimmune diseases and decipher the mechanisms underlying the dysregulation of IRF8 expression mediated by lupus risk variants.

Reanalysis of eMERGE phase III sequence variants in 10,500 participants and infrastructure to support the automated return of knowledge updates

PURPOSE

The clinical genomics knowledgebase is dynamic with variant classifications changing as newly identified cases, additional population data, and other evidence become available. This is a challenge for the clinical laboratory because of limited resource availability for variant reassessment.

METHODS

Throughout the Electronic Medical Records and Genomics phase III program, clinical sites associated with the Mass General Brigham/Broad sequencing center received automated, real-time notifications when reported variants were reclassified. In this study, we summarized the nature of these reclassifications and described the proactive reassessment framework we used for the Electronic Medical Records and Genomics program data set to identify variants most likely to undergo reclassification.

RESULTS

Reanalysis of 1855 variants led to the reclassification of 2% (n = 45) of variants, affecting 0.6% (n = 67) of participants. Of these reclassifications, 78% (n = 35) were high-impact changes affecting reportability, with 8 variants downgraded from likely pathogenic/pathogenic to variants of uncertain significance (VUS) and 27 variants upgraded from VUS to likely pathogenic/pathogenic. Most upgraded variants (67%) were initially classified as VUS-Favor Pathogenic, highlighting the benefit of VUS subcategorization. The most common reason for reclassification was new published case data and/or functional evidence.

CONCLUSION

Our results highlight the importance of periodic sequence variant reevaluation and the need for automated approaches to advance routine implementation of variant reevaluations in clinical practice.

Epstein-Barr virus nuclear antigen 2 extensively rewires the human chromatin landscape at autoimmune risk loci

The interplay between environmental and genetic factors plays a key role in the development of many autoimmune diseases. In particular, the Epstein-Barr virus (EBV) is an established contributor to multiple sclerosis, lupus, and other disorders. Previously, we showed that the EBV nuclear antigen 2 (EBNA2) transactivating protein occupies up to half of the risk loci for a set of seven autoimmune disorders. To further examine the mechanistic roles played by EBNA2 at these loci on a genome-wide scale, we globally examined gene expression, chromatin accessibility, chromatin looping, and EBNA2 binding in a B cell line that was (1) uninfected, (2) infected with a strain of EBV lacking EBNA2, or (3) infected with a strain that expresses EBNA2. We identified more than 400 EBNA2-dependent differentially expressed human genes and more than 5000 EBNA2 binding events in the human genome. ATAC-seq analysis revealed more than 2000 regions in the human genome with EBNA2-dependent chromatin accessibility, and HiChIP data revealed more than 1700 regions where EBNA2 altered chromatin looping interactions. Autoimmune genetic risk loci were highly enriched at the sites of these EBNA2-dependent chromatin-altering events. We present examples of autoimmune risk genotype-dependent EBNA2 events, nominating genetic risk mechanisms for autoimmune risk loci such as ZMIZ1 Taken together, our results reveal important interactions between host genetic variation and EBNA2-driven disease mechanisms. Further, our study highlights a critical role for EBNA2 in rewiring human gene regulatory programs through rearrangement of the chromatin landscape and nominates these interactions as components of genetic mechanisms that influence the risk of multiple autoimmune diseases.

Antibody Responses to Epstein-Barr Virus in the Preclinical Period of Rheumatoid Arthritis Suggest the Presence of Increased Viral Reactivation Cycles

OBJECTIVE

Patients with established rheumatoid arthritis (RA) demonstrate altered immune responses to Epstein-Barr virus (EBV), but the presence and role(s) of EBV have not been fully explored in the preclinical period. We hypothesized that EBV infection, as evidenced by an altered anti-EBV antibody response, could either play an important role in driving the development or be a result of expanded RA-related autoimmunity.

METHODS

83 subjects with RA based on 1987 ACR criteria and 83 matched controls were evaluated. Subject and matched control sera from the pre and post- RA diagnosis periods were tested for 5 anti-EBV antibodies (EBNA-1-IgG, VCA-IgG, EA-IgG, VCA-IgA, and EA-IgA), 7 RA-related autoantibodies (RF-neph, RF-IgA, RF-IgM, RF-IgG, CCP2, CCP3, CCP3.1), 22 cytokine/chemokine, 36 individual APCAs, and CMV-IgG antibodies. Random forest classification, mixed and joint mixed modelling were used to determine differences in anti-EBV antibodies between RA cases and controls.

RESULTS

Random Forest analysis identified preclinical EBV antibodies that differentiate RA subjects from controls. Specifically, EA-IgG antibody levels are higher in RA cases (0.82 ± 0.72) compared to controls (0.49 ± 0.28). Elevations in EA-IgG levels significantly correlated with increasing RF-IgM levels in future RA cases (p = 0.007) but not in controls (p = 0.150). CMV-IgG antibody levels did not differ between groups.

CONCLUSION

Subjects who eventually develop classified RA demonstrate elevated EA-IgG antibody levels in the preclinical period, which suggests the presence of increased EBV re-activation cycles. A combination of RF and EBV reactivation may play an important role in the development of RA.

Validation of low-coverage whole-genome sequencing for mitochondrial DNA variants suggests mitochondrial DNA as a genetic cause of preterm birth

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‐coverage 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.

Medical Records-Based Genetic Studies of the Complement System

BACKGROUND

Genetic variants in complement genes have been associated with a wide range of human disease states, but well-powered genetic association studies of complement activation have not been performed in large multiethnic cohorts.

METHODS

We performed medical records-based genome-wide and phenome-wide association studies for plasma C3 and C4 levels among participants of the Electronic Medical Records and Genomics (eMERGE) network.

RESULTS

In a GWAS for C3 levels in 3949 individuals, we detected two genome-wide significant loci: chr.1q31.3 (CFH locus; rs3753396-A; β=0.20; 95% CI, 0.14 to 0.25; P=1.52x10-11) and chr.19p13.3 (C3 locus; rs11569470-G; β=0.19; 95% CI, 0.13 to 0.24; P=1.29x10-8). These two loci explained approximately 2% of variance in C3 levels. GWAS for C4 levels involved 3998 individuals and revealed a genome-wide significant locus at chr.6p21.32 (C4 locus; rs3135353-C; β=0.40; 95% CI, 0.34 to 0.45; P=4.58x10-35). This locus explained approximately 13% of variance in C4 levels. The multiallelic copy number variant analysis defined two structural genomic C4 variants with large effect on blood C4 levels: C4-BS (β=-0.36; 95% CI, -0.42 to -0.30; P=2.98x10-22) and C4-AL-BS (β=0.25; 95% CI, 0.21 to 0.29; P=8.11x10-23). Overall, C4 levels were strongly correlated with copy numbers of C4A and C4B genes. In comprehensive phenome-wide association studies involving 102,138 eMERGE participants, we cataloged a full spectrum of autoimmune, cardiometabolic, and kidney diseases genetically related to systemic complement activation.

CONCLUSIONS

We discovered genetic determinants of plasma C3 and C4 levels using eMERGE genomic data linked to electronic medical records. Genetic variants regulating C3 and C4 levels have large effects and multiple clinical correlations across the spectrum of complement-related diseases in humans.

Pleiotropy of systemic lupus erythematosus risk alleles and cardiometabolic disorders: A phenome-wide association study and inverse-variance weighted meta-analysis

OBJECTIVES

To test the hypothesis that genetic predisposition to systemic lupus erythematosus (SLE) increases the risk of cardiometabolic disorders.

METHODS

Using 41 single nucleotide polymorphisms (SNPs) associated with SLE, we calculated a weighted genetic risk score (wGRS) for SLE. In a large biobank we tested the association between this wGRS and 9 cardiometabolic phenotypes previously associated with SLE: atrial fibrillation, ischemic stroke, coronary artery disease, type 1 and type 2 diabetes, obesity, chronic kidney disease, hypertension, and hypercholesterolemia. Additionally, we performed a phenome-wide association analysis (pheWAS) to discover novel clinical associations with a genetic predisposition to SLE. Findings were replicated in the Electronic Medical Records and Genomics (eMERGE) Network. To further define the association between SLE-related risk alleles and the selected cardiometabolic phenotypes, we performed an inverse variance weighted regression (IVWR) meta-analysis.

RESULTS

The wGRS for SLE was calculated in 74,759 individuals of European ancestry. Among the pre-selected phenotypes, the wGRS was significantly associated with type 1 diabetes (OR [95%CI] =1.11 [1.06, 1.17], P-value = 1.05x10-5). In the PheWAS, the wGRS was associated with several autoimmune phenotypes, kidney disorders, and skin neoplasm; but only the associations with autoimmune phenotypes were replicated. In the IVWR meta-analysis, SLE-related risk alleles were nominally associated with type 1 diabetes (P = 0.048) but the associations were heterogeneous and did not meet the adjusted significance threshold.

CONCLUSION

A weighted GRS for SLE was associated with an increased risk of several autoimmune-related phenotypes including type I diabetes but not with cardiometabolic disorders.

Novel EDGE encoding method enhances ability to identify genetic interactions

Assumptions are made about the genetic model of single nucleotide polymorphisms (SNPs) when choosing a traditional genetic encoding: additive, dominant, and recessive. Furthermore, SNPs across the genome are unlikely to demonstrate identical genetic models. However, running SNP-SNP interaction analyses with every combination of encodings raises the multiple testing burden. Here, we present a novel and flexible encoding for genetic interactions, the elastic data-driven genetic encoding (EDGE), in which SNPs are assigned a heterozygous value based on the genetic model they demonstrate in a dataset prior to interaction testing. We assessed the power of EDGE to detect genetic interactions using 29 combinations of simulated genetic models and found it outperformed the traditional encoding methods across 10%, 30%, and 50% minor allele frequencies (MAFs). Further, EDGE maintained a low false-positive rate, while additive and dominant encodings demonstrated inflation. We evaluated EDGE and the traditional encodings with genetic data from the Electronic Medical Records and Genomics (eMERGE) Network for five phenotypes: age-related macular degeneration (AMD), age-related cataract, glaucoma, type 2 diabetes (T2D), and resistant hypertension. A multi-encoding genome-wide association study (GWAS) for each phenotype was performed using the traditional encodings, and the top results of the multi-encoding GWAS were considered for SNP-SNP interaction using the traditional encodings and EDGE. EDGE identified a novel SNP-SNP interaction for age-related cataract that no other method identified: rs7787286 (MAF: 0.041; intergenic region of chromosome 7)–rs4695885 (MAF: 0.34; intergenic region of chromosome 4) with a Bonferroni LRT p of 0.018. A SNP-SNP interaction was found in data from the UK Biobank within 25 kb of these SNPs using the recessive encoding: rs60374751 (MAF: 0.030) and rs6843594 (MAF: 0.34) (Bonferroni LRT p: 0.026). We recommend using EDGE to flexibly detect interactions between SNPs exhibiting diverse action.

Although traditional genetic encodings are widely implemented in genetics research, including in genome-wide association studies (GWAS) and epistasis, each method makes assumptions that may not reflect the underlying etiology. Here, we introduce a novel encoding method that estimates and assigns an individualized data-driven encoding for each single nucleotide polymorphism (SNP): the elastic data-driven genetic encoding (EDGE). With simulations, we demonstrate that this novel method is more accurate and robust than traditional encoding methods in estimating heterozygous genotype values, reducing the type I error, and detecting SNP-SNP interactions. We further applied the traditional encodings and EDGE to biomedical data from the Electronic Medical Records and Genomics (eMERGE) Network for five phenotypes, and EDGE identified a novel interaction for age-related cataract not detected by traditional methods, which replicated in data from the UK Biobank. EDGE provides an alternative approach to understanding and modeling diverse SNP models and is recommended for studying complex genetics in common human phenotypes.

Lupus susceptibility region containing CDKN1B rs34330 mechanistically influences expression and function of multiple target genes, also linked to proliferation and apoptosis

OBJECTIVE

A recent genome-wide association study (GWAS) reported a significant genetic association between rs34330 of cyclin-dependent kinase inhibitor 1B (CDKN1B) and risk of systemic lupus erythematosus (SLE) in Han Chinese. This study aims to validate the reported association and elucidate the biochemical mechanisms underlying the variant's effect.

METHODS

We performed allelic association with SLE followed by meta-analysis across 11 independent cohorts (n=28,872). We applied in silico bioinformatics and experimental validation in SLE-relevant cell lines to determine the functional consequences of rs34330.

RESULTS

We replicated genetic association between SLE and rs34330 (P_meta =5.29x10^-22 , OR (95% CI)=0.84 (0.81-0.87)). Follow-up bioinformatics and eQTL analysis suggest that rs34330 is located in active chromatin and potentially regulates several target genes. Using luciferase and ChIP-qPCR, we demonstrated substantial allele-specific promoter and enhancer activity, and allele-specific binding of three histone marks (H3K27ac, H3K4me3, H3K4me1), RNA pol II, CTCF, and a critical immune transcription factor (IRF-1). Chromosome conformation capture (3C) detected long-range chromatin interactions between rs34330 and the promoters of neighboring genes APOLD1 and DDX47, and effects on CDKN1B and the other target genes were directly validated by CRISPR-based genome editing. Finally, CRISPR-dCas9-based epigenetic activation/silencing confirmed these results. Gene-edited cell lines also showed higher levels of proliferation and apoptosis.

CONCLUSION

Collectively, these findings suggest a mechanism whereby the rs34330 risk allele (C) influences the presence of histone marks, RNA pol II, and the IRF-1 transcription factor to regulate expression of several target genes linked to proliferation and apoptosis, which potentially underlie the association of rs34330 with SLE.

Meta-analysis of 208370 East Asians identifies 113 susceptibility loci for systemic lupus erythematosus

OBJECTIVE

Systemic lupus erythematosus (SLE), an autoimmune disorder, has been associated with nearly 100 susceptibility loci. Nevertheless, these loci only partially explain SLE heritability and their putative causal variants are rarely prioritised, which make challenging to elucidate disease biology. To detect new SLE loci and causal variants, we performed the largest genome-wide meta-analysis for SLE in East Asian populations.

METHODS

We newly genotyped 10 029 SLE cases and 180 167 controls and subsequently meta-analysed them jointly with 3348 SLE cases and 14 826 controls from published studies in East Asians. We further applied a Bayesian statistical approach to localise the putative causal variants for SLE associations.

RESULTS

We identified 113 genetic regions including 46 novel loci at genome-wide significance (p<5×10^-8). Conditional analysis detected 233 association signals within these loci, which suggest widespread allelic heterogeneity. We detected genome-wide associations at six new missense variants. Bayesian statistical fine-mapping analysis prioritised the putative causal variants to a small set of variants (95% credible set size ≤10) for 28 association signals. We identified 110 putative causal variants with posterior probabilities ≥0.1 for 57 SLE loci, among which we prioritised 10 most likely putative causal variants (posterior probability ≥0.8). Linkage disequilibrium score regression detected genetic correlations for SLE with albumin/globulin ratio (r_g=-0.242) and non-albumin protein (r_g=0.238).

CONCLUSION

This study reiterates the power of large-scale genome-wide meta-analysis for novel genetic discovery. These findings shed light on genetic and biological understandings of SLE.

Patients with Proliferative Lupus Nephritis Have Autoantibodies That React to Moesin and Demonstrate Increased Glomerular Moesin Expression

Kidney involvement in systemic lupus erythematosus (SLE)—termed lupus nephritis (LN)—is a severe manifestation of SLE that can lead to end-stage kidney disease (ESKD). LN is characterized by immune complex deposition and inflammation in the glomerulus. We tested the hypothesis that autoantibodies targeting podocyte and glomerular cell proteins contribute to the development of immune complex formation in LN. We used Western blotting with SLE sera from patients with and without LN to identify target antigens in human glomerular and cultured human-derived podocyte membrane proteins. Using liquid chromatography-tandem mass spectrometry (LC-MS/MS), we identified the proteins in the gel regions corresponding to reactive bands observed with sera from LN patients. We identified 102 proteins that were present in both the podocyte and glomerular samples. We identified 10 high-probability candidates, including moesin, using bioinformatic analysis. Confirmation of moesin as a target antigen was conducted using immunohistochemical analysis (IHC) of kidney biopsy tissue and enzyme-linked immunosorbent assay (ELISA) to detect circulating antibodies. By IHC, biopsies from patients with proliferative lupus nephritis (PLN, class III/IV) demonstrated significantly increased glomerular expression of moesin (p < 0.01). By ELISA, patients with proliferative LN demonstrated significantly increased antibodies against moesin (p < 0.01). This suggests that moesin is a target glomerular antigen in lupus nephritis.

SLE non-coding genetic risk variant determines the epigenetic dysfunction of an immune cell specific enhancer that controls disease-critical microRNA expression

Since most variants that impact polygenic disease phenotypes localize to non-coding genomic regions, understanding the consequences of regulatory element variants will advance understanding of human disease mechanisms. Here, we report that the systemic lupus erythematosus (SLE) risk variant rs2431697 as likely causal for SLE through disruption of a regulatory element, modulating miR-146a expression. Using epigenomic analysis, genome-editing and 3D chromatin structure analysis, we show that rs2431697 tags a cell-type dependent distal enhancer specific for miR-146a that physically interacts with the miR-146a promoter. NF-kB binds the disease protective allele in a sequence-specific manner, increasing expression of this immunoregulatory microRNA. Finally, CRISPR activation-based modulation of this enhancer in the PBMCs of SLE patients attenuates type I interferon pathway activation by increasing miR-146a expression. Our work provides a strategy to define non-coding RNA functional regulatory elements using disease-associated variants and provides mechanistic links between autoimmune disease risk genetic variation and disease etiology.

Comprehensive Review of Steroid-Sensitive Nephrotic Syndrome Genetic Risk Loci and Transcriptional Regulation as a Possible Mechanistic Link to Disease Risk

Introduction

The etiology of steroid-sensitive nephrotic syndrome (SSNS) is not well understood. Genetic studies have established common single nucleotide polymorphisms (SNPs) that are associated with increased SSNS disease risk. We review previous genetic association studies of SSNS and nominate particular transcriptional regulators and immune cells as potential key players in the etiology of this disease.

Methods

A list of SNPs associated with SSNS was compiled from published genome wide association and candidate gene studies. The Regulatory Element Locus Intersection (RELI) tool was used to calculate the enrichment of the overlap between disease risk SNPs and the genomic coordinates of data from a collection of >10,000 chromatin immunoprecipitation sequencing experiments.

Results

After linkage disequilibrium expansion of the previously reported tag associated SNPs, we identified 192 genetic variants at 8 independent risk loci. Using the Regulatory Element Locus Intersection algorithm, we identified transcriptional regulators with enriched binding at SSNS risk loci (10^-05 < P_corrected < 10^-124), including ZNF530, CIITA, CD74, RFX5, and ZNF425. Many of these regulators have well-described roles in the immune response. RNA polymerase II binding in B cells also demonstrated enriched binding at SSNS risk loci (10^-37<P_corrected<10^-5).

Conclusion

SSNS is a complex disease, and immune dysregulation has been previously implicated as a potential underlying cause. This assessment of established SSNS risk loci and analysis of possible function implicates transcriptional dysregulation, and specifically particular transcriptional regulators with known roles in the immune response, as important in the genetic etiology of SSNS.

Relationship Between Genetic Risk and Age of Diagnosis in Systemic Lupus Erythematosus

OBJECTIVE

Specific risk alleles for childhood-onset systemic lupus erythematosus SLE (cSLE) vs adult-onset SLE (aSLE) patients have not been identified. The aims of this study were to determine if there is an association (1) between non-HLA-related genetic risk score (GRS) and age of SLE diagnosis, and (2) between HLA-related GRS and age of SLE diagnosis.

METHODS

Genomic DNA was obtained from 2001 multiethnic patients and genotyped using the Immunochip. Following quality control, genetic risk counting (GRCS), weighted (GRWS), standardized counting (GRSCS), and standardized weighted (GRSWS) scores were calculated based on independent single-nucleotide polymorphisms from validated SLE loci. Scores were analyzed in a regression model and adjusted by sex and ancestral population.

RESULTS

The analyzed cohort consisted of 1540 patients: 1351 females and 189 males (675 cSLE and 865 aSLE). There were significant negative associations between all non-HLA GRS and age of SLE diagnosis: P = 0.011 and r² = 0.175 for GRWS; P = 0.008 and r² = 0.178 for GRSCS; P = 0.002 and r² = 0.176 for GRSWS (higher GRS correlated with lower age of diagnosis.) All HLA GRS showed significant positive associations with age of diagnosis: P = 0.049 and r² = 0.176 for GRCS; P = 0.022 and r² = 0.176 for GRWS; P = 0.022 and r² = 0.176 for GRSCS; P = 0.011 and r² = 0.177 for GRSWS (higher GRS correlated with higher age of diagnosis).

CONCLUSION

Our data suggest that there is a linear relationship between genetic risk and age of SLE diagnosis and that HLA and non-HLA GRS are associated with age of diagnosis in opposite directions.

Genome-wide Modeling of Polygenic Risk Score in Colorectal Cancer Risk

Accurate colorectal cancer (CRC) risk prediction models are critical for identifying individuals at low and high risk of developing CRC, as they can then be offered targeted screening and interventions to address their risks of developing disease (if they are in a high-risk group) and avoid unnecessary screening and interventions (if they are in a low-risk group). As it is likely that thousands of genetic variants contribute to CRC risk, it is clinically important to investigate whether these genetic variants can be used jointly for CRC risk prediction. In this paper, we derived and compared different approaches to generating predictive polygenic risk scores (PRS) from genome-wide association studies (GWASs) including 55,105 CRC-affected case subjects and 65,079 control subjects of European ancestry. We built the PRS in three ways, using (1) 140 previously identified and validated CRC loci; (2) SNP selection based on linkage disequilibrium (LD) clumping followed by machine-learning approaches; and (3) LDpred, a Bayesian approach for genome-wide risk prediction. We tested the PRS in an independent cohort of 101,987 individuals with 1,699 CRC-affected case subjects. The discriminatory accuracy, calculated by the age- and sex-adjusted area under the receiver operating characteristics curve (AUC), was highest for the LDpred-derived PRS (AUC = 0.654) including nearly 1.2 M genetic variants (the proportion of causal genetic variants for CRC assumed to be 0.003), whereas the PRS of the 140 known variants identified from GWASs had the lowest AUC (AUC = 0.629). Based on the LDpred-derived PRS, we are able to identify 30% of individuals without a family history as having risk for CRC similar to those with a family history of CRC, whereas the PRS based on known GWAS variants identified only top 10% as having a similar relative risk. About 90% of these individuals have no family history and would have been considered average risk under current screening guidelines, but might benefit from earlier screening. The developed PRS offers a way for risk-stratified CRC screening and other targeted interventions.

Pleiotropy in the Genetic Predisposition to Rheumatoid Arthritis: A Phenome-Wide Association Study and Inverse Variance-Weighted Meta-Analysis

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.

A Polygenic and Phenotypic Risk Prediction for Polycystic Ovary Syndrome Evaluated by Phenome-Wide Association Studies

CONTEXT

As many as 75% of patients with polycystic ovary syndrome (PCOS) are estimated to be unidentified in clinical practice.

OBJECTIVE

Utilizing polygenic risk prediction, we aim to identify the phenome-wide comorbidity patterns characteristic of PCOS to improve accurate diagnosis and preventive treatment.

DESIGN, PATIENTS, AND METHODS

Leveraging the electronic health records (EHRs) of 124 852 individuals, we developed a PCOS risk prediction algorithm by combining polygenic risk scores (PRS) with PCOS component phenotypes into a polygenic and phenotypic risk score (PPRS). We evaluated its predictive capability across different ancestries and perform a PRS-based phenome-wide association study (PheWAS) to assess the phenomic expression of the heightened risk of PCOS.

RESULTS

The integrated polygenic prediction improved the average performance (pseudo-R2) for PCOS detection by 0.228 (61.5-fold), 0.224 (58.8-fold), 0.211 (57.0-fold) over the null model across European, African, and multi-ancestry participants respectively. The subsequent PRS-powered PheWAS identified a high level of shared biology between PCOS and a range of metabolic and endocrine outcomes, especially with obesity and diabetes: "morbid obesity", "type 2 diabetes", "hypercholesterolemia", "disorders of lipid metabolism", "hypertension", and "sleep apnea" reaching phenome-wide significance.

CONCLUSIONS

Our study has expanded the methodological utility of PRS in patient stratification and risk prediction, especially in a multifactorial condition like PCOS, across different genetic origins. By utilizing the individual genome-phenome data available from the EHR, our approach also demonstrates that polygenic prediction by PRS can provide valuable opportunities to discover the pleiotropic phenomic network associated with PCOS pathogenesis.

Complement genes contribute sex-biased vulnerability in diverse disorders

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 men¹, whereas schizophrenia affects men more frequently and severely². 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 alleles of the human leukocyte antigen (HLA) genes at that locus^3–6. Here we show that the complement component 4 (C4) genes, which are also in the MHC locus and were recently found to increase risk for schizophrenia⁷, generate 7-fold variation in risk for lupus (95% CI: 5.88–8.61; p < 10⁻¹¹⁷ in total) and 16-fold variation in risk for Sjögren’s syndrome (95% CI: 8.59–30.89; p < 10⁻²³ in total) among individuals with common C4 genotypes, 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, 31-fold variation in risk for Sjögren’s syndrome, and 1.7-fold variation in schizophrenia risk among men (vs. 6-fold, 15-fold, and 1.26-fold among women respectively). At a protein level, both C4 and its effector C3 were present at greater levels in men than women in cerebrospinal fluid (p < 10⁻⁵ for both C4 and C3) and plasma^8,9 among adults ages 20–50, 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 implicate the complement system as a source of sexual dimorphism in vulnerability to diverse illnesses.

The Ikaros1, IKZF1, risk locus for Acute Lymphoblastic Leukemia (ALL) at chromosome 7p12.2 regulates the expression of the FIGNL1 gene

ALL, the most common pediatric cancer, originates from progenitor B-cells in ~85% of cases. IKAROS Family Zinc Finger 1 (IKZF1) is an important locus for ALL. Our goal was to identify the functional variants, thereby nominating a potential causal mechanism for ALL.

We evaluated all 27 genetic polymorphisms in disequilibrium at r2&gt;0.8 with the tag, the most highly associated risk variant at IKZF1 locus. They map to intron 7, 3′ UTR and 3′ flanking sequence of IKZF1. Fidgetin-like protein 1 (FIGNL1) is 45,524 bases from the tag variant and is nominated by eQTL and DNA looping data along with IKZF1 as possible targets. Epstein-Barr virus (EBV) transformed GM12878 cells were transduced with lentiviral vectors for stable expression of fusion proteins dead (d) Cas9-VP64 and activator MS2-p65-HSF1 (synergistic activation mediator (SAM) system). We designed 27 single guide (sg) RNAs, one for each of the 27 variants and evaluated expression perturbation in cells transfected with pools and individual sgRNA. Expression of FIGNL1, but not IKZF1, was convincingly altered in this system. The reduction of FIGNL1 expression by ~50% was isolated to two neighboring variants separated by 1 kb. The risk allele at these variants increases FIGNL1 expression in wild type cells, however in EBV transformed B-cells FIGNL1 expression is decreased. Other data show that the EBV encoded DNA binding co-factor EBNA3C binds at or very near this locus in EBV transformed B cells.

Our results nominate two closely located SNPs in the last intron of IKZF1 as potentially causal variants, controlling the expression of FIGNL1, which ordinarily participates in homologous recombination during DNA double-strand break repair, maintenance of genomic stability and prevention of cancer.

Conditional abrogation of Ddx3x reveals its essential role in lymphocyte development and partial compensation by Y chromosome

Ddx3x on the X-chromosome encodes a DEAD-box RNA helicase with putative roles in RNA metabolism, cell cycle progression, apoptosis, and viral immunity. We previously found that hematopoietic-specific hemizygosity of Ddx3x studies caused a marked reduction in the frequencies of natural killer (NK) and B cells in peripheral lymphoid organs of male mice (Liu et al. bioRxiv 2018). No homozygous Ddx3x-deficient female pups were ever obtained, suggesting that hematopoietic loss of this gene contributes to the embryonic lethality observed in germline Ddx3x knockout mice. To assess cell-lineage intrinsic roles for Ddx3x in B-cell and NK-cell development, we bred Ddx3x-floxed mice with Mb1-Cre and Ncr1-iCre mice, respectively. While hemizygous male Ddx3x-deficient mice display mild reductions in B-cell (Mb1-Cre) and NK-cell (Ncr1-iCre) frequencies, homozygous Ddx3x-deficient female mice from each were almost completely devoid of the relevant cell lineages. B-cell loss was evident from the pro-B cells stage, with B cells surviving beyond this stage showing evidence of having escaped deletion of Ddx3x. Ncr1-iCre-mediated deletion of Ddx3x completely abrogated development of both NK cells and other Ncr1-expressing innate lymphoid cells. Thus, DDX3 plays an essential role in the lymphocyte development, where Ddx3y may partially compensate for loss of Ddx3x in these lineages.

Latency III gene products of Epstein-Barr Virus (EBV) are associated with Systemic Lupus Erythematosus (SLE) genetic risk loci

SLE affects millions with increasing diagnostic prevalence. Genetic studies have identified &gt;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 &gt;100 risk alleles associated at p&lt;5E-8 for SLE were curated and reduced to 83 loci by linkage disequilibrium pruning. We evaluated the available 53 virally encoded and 13,051 human TF ChIP-seq (chromatin immunoprecipitation with DNA sequencing) datasets. In addition to EBNA2 (now OR=3.7, Pc=1.96E-18 after Bonferroni correction), two of the 53 viral TF data sets, EBNA3C and EBNALP, each of which are EBV Latency III gene products, were also strongly associated with the SLE risk loci (OR=6.0, Pc=4.19E-18; and OR=3.3, Pc=5.53E-18). Moreover, among the human TF datasets strongly associated with SLE risk loci (Pc&lt;E-6), those from EBV infected and transformed B cells are highly concentrated (OR≈56, P&lt;E-100) and, in addition, the human TFs that are known to form super-enhancers upon B cell transformation by EBV (e.g., NFκB subunits) also tend to be concentrated at the SLE risk loci (OR&gt;30, p&lt;E-25). Meanwhile, genetic associations for depression, anxiety, and schizophrenia show no such relationships.

These new results further the possible validity of the hypothesis that SLE is originally caused by the Latency III expression program of EBV.

Genetic influences on susceptibility to rheumatoid arthritis in African-Americans

Large meta-analyses of rheumatoid arthritis (RA) susceptibility in European (EUR) and East Asian (EAS) populations have identified >100 RA risk loci, but genome-wide studies of RA in African-Americans (AAs) are absent. To address this disparity, we performed an analysis of 916 AA RA patients and 1392 controls and aggregated our data with genotyping data from >100 000 EUR and Asian RA patients and controls. We identified two novel risk loci that appear to be specific to AAs: GPC5 and RBFOX1 (PAA < 5 × 10-9). Most RA risk loci are shared across different ethnicities, but among discordant loci, we observed strong enrichment of variants having large effect sizes. We found strong evidence of effect concordance for only 3 of the 21 largest effect index variants in EURs. We used the trans-ethnic fine-mapping algorithm PAINTOR3 to prioritize risk variants in >90 RA risk loci. Addition of AA data to those of EUR and EAS descent enabled identification of seven novel high-confidence candidate pathogenic variants (defined by posterior probability > 0.8). In summary, our trans-ethnic analyses are the first to include AAs, identified several new RA risk loci and point to candidate pathogenic variants that may underlie this common autoimmune disease. These findings may lead to better ways to diagnose or stratify treatment approaches in RA.

Epstein Barr virus nuclear antigen 1 (EBNA-1) peptides recognized by adult multiple sclerosis patient sera induce neurologic symptoms in a murine model

Multiple sclerosis (MS) is an autoimmune demyelinating disease with progressive neurodegeneration and complex etiology likely involving genetic and environmental factors. MS has been associated with Epstein Barr virus (EBV) infection, with patients often showing enhanced responses to EBV antigens. To determine whether abnormal EBV nuclear antigen-1 (EBNA-1) humoral immunity can serve as an initiator of autoimmune responses in MS, we investigated the fine specificities of the humoral immune response against EBNA-1 in MS patients using solid phase epitope mapping. Antibodies from MS patients recognized an EBNA-1 epitope spanning amino acids 411-426, previously unknown to be recognized specifically by untreated MS patients. Antibodies against this epitope cross-reacted to myelin basic protein (MBP). Furthermore, animals immunized with this EBNA-1 polypeptide mounted a response against MBP and developed signs of experimental autoimmune encephalitis (EAE). These data support a link between MS and EBV through antibodies that cross-react between EBV proteins and the MBP autoantigen.

Complete Tracheal Ring Deformity. A Translational Genomics Approach to Pathogenesis

Rationale: Complete tracheal ring deformity (CTRD) is a rare congenital abnormality of unknown etiology characterized by circumferentially continuous or nearly continuous cartilaginous tracheal rings, variable degrees of tracheal stenosis and/or shortening, and/or pulmonary arterial sling anomaly.Objectives: To test the hypothesis that CTRD is caused by inherited or de novo mutations in genes required for normal tracheal development.Methods: CTRD and normal tracheal tissues were examined microscopically to define the tracheal abnormalities present in CTRD. Whole-exome sequencing was performed in children with CTRD and their biological parents ("trio analysis") to identify gene variants in patients with CTRD. Mutations were confirmed by Sanger sequencing, and their potential impact on structure and/or function of encoded proteins was examined using human gene mutation databases. Relevance was further examined by comparison with the effects of targeted deletion of murine homologs important to tracheal development in mice.Measurements and Main Results: The trachealis muscle was absent in all of five patients with CTRD. Exome analysis identified six de novo, three recessive, and multiple compound-heterozygous or rare hemizygous variants in children with CTRD. De novo variants were identified in SHH (Sonic Hedgehog), and inherited variants were identified in HSPG2 (perlecan), ROR2 (receptor tyrosine kinase-like orphan receptor 2), and WLS (Wntless), genes involved in morphogenetic pathways known to mediate tracheoesophageal development in mice.Conclusions: The results of the present study demonstrate that absence of the trachealis muscle is associated with CTRD. Variants predicted to cause disease were identified in genes encoding Hedgehog and Wnt signaling pathway molecules, which are critical to cartilage formation and normal upper airway development in mice.

Characterization of cxorf21 Provides Molecular Insight Into Female-Bias Immune Response in SLE Pathogenesis

Background: Ninety percent of systemic lupus erythematosus (SLE) patients are women. X chromosome-dosage increases susceptibility to SLE and primary Sjögren's syndrome (pSS). Chromosome X open reading frame 21 (CXorf21) escapes X-inactivation and is an SLE risk gene of previously unknown function. We undertook the present study to delineate the function of CXorf21 in the immune system as well as investigate a potential role in the sex bias of SLE and pSS. Methods: Western blot protein analysis, qPCR, BioPlex cytokine immunoassay, pHrodo™ assays, as well as in vitro CRISPR-Cas9 knockdown experiments were employed to delineate the role of CXorf21 in relevant immunocytes. Results: Expressed in monocytes and B cells, CXorf21 basal Mrna, and protein expression levels are elevated in female primary monocytes, B cells, and EBV-transformed B cells compared to male cells. We also found CXorf21 mRNA and protein expression is higher in both male and female cells from SLE patients compared to control subjects. TLR7 ligation increased CXorf21 protein expression and CXorf21 knockdown abrogated TLR7-driven increased IFNA1 mRNA expression, and reduced secretion of both TNF-alpha and IL-6 in healthy female monocytes. Similarly, we found increased pH in the lysosomes of CXorf21-deficient female monocytes. Conclusion: CXorf21 is more highly expressed in female compared to male cells and is involved in a sexually dimorphic response to TLR7 activation. In addition, CXorf21 expression regulates lysosomal pH in a sexually dimorphic manner. Thus, sexually dimorphic expression of CXorf21 skews cellular immune responses in manner consistent with expected properties of a mediator of the X chromosome dose risk in SLE and pSS.

GWAS and enrichment analyses of non-alcoholic fatty liver disease identify new trait-associated genes and pathways across eMERGE Network

BACKGROUND

Non-alcoholic fatty liver disease (NAFLD) is a common chronic liver illness with a genetically heterogeneous background that can be accompanied by considerable morbidity and attendant health care costs. The pathogenesis and progression of NAFLD is complex with many unanswered questions. We conducted genome-wide association studies (GWASs) using both adult and pediatric participants from the Electronic Medical Records and Genomics (eMERGE) Network to identify novel genetic contributors to this condition.

METHODS

First, a natural language processing (NLP) algorithm was developed, tested, and deployed at each site to identify 1106 NAFLD cases and 8571 controls and histological data from liver tissue in 235 available participants. These include 1242 pediatric participants (396 cases, 846 controls). The algorithm included billing codes, text queries, laboratory values, and medication records. Next, GWASs were performed on NAFLD cases and controls and case-only analyses using histologic scores and liver function tests adjusting for age, sex, site, ancestry, PC, and body mass index (BMI).

RESULTS

Consistent with previous results, a robust association was detected for the PNPLA3 gene cluster in participants with European ancestry. At the PNPLA3-SAMM50 region, three SNPs, rs738409, rs738408, and rs3747207, showed strongest association (best SNP rs738409 p = 1.70 × 10- 20). This effect was consistent in both pediatric (p = 9.92 × 10- 6) and adult (p = 9.73 × 10- 15) cohorts. Additionally, this variant was also associated with disease severity and NAFLD Activity Score (NAS) (p = 3.94 × 10- 8, beta = 0.85). PheWAS analysis link this locus to a spectrum of liver diseases beyond NAFLD with a novel negative correlation with gout (p = 1.09 × 10- 4). We also identified novel loci for NAFLD disease severity, including one novel locus for NAS score near IL17RA (rs5748926, p = 3.80 × 10- 8), and another near ZFP90-CDH1 for fibrosis (rs698718, p = 2.74 × 10- 11). Post-GWAS and gene-based analyses identified more than 300 genes that were used for functional and pathway enrichment analyses.

CONCLUSIONS

In summary, this study demonstrates clear confirmation of a previously described NAFLD risk locus and several novel associations. Further collaborative studies including an ethnically diverse population with well-characterized liver histologic features of NAFLD are needed to further validate the novel findings.

Association of Epstein-Barr virus serological reactivation with transitioning to systemic lupus erythematosus in at-risk individuals

Objective.

Systemic lupus erythematosus (SLE) is a systemic autoimmune disease with unknown etiology. Epstein-Barr virus (EBV) is an environmental factor associated with SLE. EBV maintains latency in B cells with frequent reactivation measured by antibodies against viral capsid antigen (VCA) and early antigen (EA). In this study, we determined whether EBV reactivation and single nucleotide polymorphisms (SNPs) in EBV-associated host genes are associated with SLE transition.

Methods.

SLE patient relatives (n=436) who did not have SLE at baseline were re-contacted after 6.3 (±3.9) years and evaluated for interim transitioning to SLE (≥4 cumulative ACR criteria); 56 (13%) transitioned to SLE prior to the follow-up visit. At both visits, detailed demographic, environmental, clinical information, and blood samples were obtained. Antibodies against viral antigens were measured by ELISA. SNPs in IL10, CR2, TNFAIP3, and CD40 genes were typed by ImmunoChip™. Generalized estimating equations were used to test associations between viral antibody levels and transitioning to SLE.

Results.

Mean baseline VCA IgG (4.879±1.797 vs 3.866±1.795, p=0.0003) and EA IgG (1.192±1.113 vs 0.7774±0.8484, p=0.0236) levels were higher in transitioned compared to autoantibody negative non-transition relatives. Increased VCA IgG and EA IgG were associated with transitioning to SLE (OR 1.28 95%CI 1.07–1.53 p=0.007, OR 1.43 95%CI 1.06–1.93 p=0.02, respectively). Significant interactions were observed between CD40 variant rs48100485 and VCA IgG levels, and IL10 variant rs3024493 and VCA IgA levels in transitioning to SLE.

Conclusion.

Heightened serologic reactivation of EBV increases the probability of transitioning to SLE in unaffected SLE relatives.

New role of Epstein-Barr virus in pathogenesis of acute and chronic lymphocytic leukemia

Acute lymphoblastic leukemia (ALL) is the most common mortal cancer in children. Chronic lymphocytic leukemia (CLL) is the most prevalent form of adult leukemia in western countries. B lymphocytes dominate the origin of both diseases. We have found a possible role for Epstein-Barr virus (EBV) in the origins of both ALL and CLL.

We applied our strategy (Nat Genet 50:699, 2018) to determine whether the binding of EBV transcription factors (TFs) was concentrate at the 84 known risk loci for CLL and 16 loci for ALL. We evaluated 52 virally encoded TF ChIP-seq (chromatin immunoprecipitation with DNA sequencing) datasets and complemented this analysis with the results from 1535 human TF ChIP-seq datasets.

We found that Epstein-Barr nuclear antigen leader protein (EBNALP), EBNA3C and EBNA2 were concentrated in the CLL loci by 3.7, 3.7 and 3.5-fold with p=4.89*10−19, p=2.74*10−11 and p=1.07*10−8, respectively. The viral (n=3) and human TFs (n=40) cluster together in an optimal subset of ~15 of the 84 known loci in CLL at p&lt;10−6. Eighty percent of the most highly associated viral and human TF ChIP-seq datasets were collected from EBV transformed B cell lines in the Latency III program of viral expression, for which EBNALP, EBNA3C and EBNA2 are viral gene products.

In ALL 3 of 16 loci were occupied by EBNA3ABC by a 17.7-fold enrichment with p=2.0*10−10. Among human TFs only c-MYC reached statistical significance binding 2 non-overlapping risk loci with 28.4-fold enrichment with p=2.9*10−10.

These results nominate EBV for a role in the pathogenesis of CLL and ALL by a mechanism operating in transformed B cells through the EBV Latency III program of viral expression.

Amino acid signatures of HLA Class-I and II molecules are strongly associated with SLE susceptibility and autoantibody production in Eastern Asians

Human leukocyte antigen (HLA) is a key genetic factor conferring risk of systemic lupus erythematosus (SLE), but precise independent localization of HLA effects is extremely challenging. As a result, the contribution of specific HLA alleles and amino-acid residues to the overall risk of SLE and to risk of specific autoantibodies are far from completely understood. Here, we dissected (a) overall SLE association signals across HLA, (b) HLA-peptide interaction, and (c) residue-autoantibody association. Classical alleles, SNPs, and amino-acid residues of eight HLA genes were imputed across 4,915 SLE cases and 13,513 controls from Eastern Asia. We performed association followed by conditional analysis across HLA, assessing both overall SLE risk and risk of autoantibody production. DR15 alleles HLA-DRB1*15:01 (P = 1.4x10^-27, odds ratio (OR) = 1.57) and HLA-DQB1*06:02 (P = 7.4x10^-23, OR = 1.55) formed the most significant haplotype (OR = 2.33). Conditioned protein-residue signals were stronger than allele signals and mapped predominantly to HLA-DRB1 residue 13 (P = 2.2x10^-75) and its proxy position 11 (P = 1.1x10^-67), followed by HLA-DRB1-37 (P = 4.5x10^-24). After conditioning on HLA-DRB1, novel associations at HLA-A-70 (P = 1.4x10^-8), HLA-DPB1-35 (P = 9.0x10^-16), HLA-DQB1-37 (P = 2.7x10^-14), and HLA-B-9 (P = 6.5x10^-15) emerged. Together, these seven residues increased the proportion of explained heritability due to HLA to 2.6%. Risk residues for both overall disease and hallmark autoantibodies (i.e., nRNP: DRB1-11, P = 2.0x10^-14; DRB1-13, P = 2.9x10^-13; DRB1-30, P = 3.9x10^-14) localized to the peptide-binding groove of HLA-DRB1. Enrichment for specific amino-acid characteristics in the peptide-binding groove correlated with overall SLE risk and with autoantibody presence. Risk residues were in primarily negatively charged side-chains, in contrast with rheumatoid arthritis. We identified novel SLE signals in HLA Class I loci (HLA-A, HLA-B), and localized primary Class II signals to five residues in HLA-DRB1, HLA-DPB1, and HLA-DQB1. These findings provide insights about the mechanisms by which the risk residues interact with each other to produce autoantibodies and are involved in SLE pathophysiology.

The Human leukocyte antigen (HLA) region is a key genetic factor conferring risk of systemic lupus erythematosus (SLE). In spite of multiple SLE association signals identified in the HLA region, only amino-acid residues within HLA-DRB1 have been specifically described previously. In this study, we performed an imputation-based analysis on individuals with East Asian ancestry, and characterized SLE risk within the HLA region for all involved independent genes (HLA-DRB1, HLA-DPB1, HLA-DQB1, HLA-A, and HLA-B). Furthermore, we identified a characteristic SLE risk residue signature as well as a pattern of specific nRNP and Ro/La autoantibody residues located in the peptide-binding grooves, suggesting their key involvement in autoantibody production.

Genetic fine mapping of systemic lupus erythematosus MHC associations in Europeans and African Americans

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.

Latent autoimmunity across disease-specific boundaries in at-risk first-degree relatives of SLE and RA patients

BACKGROUND

Autoimmune disease prevention requires tools to assess an individual's risk of developing a specific disease. One tool is disease-associated autoantibodies, which accumulate in an asymptomatic preclinical period. However, patients sometimes exhibit autoantibodies associated with a different disease classification. When and how these alternative autoantibodies first appear remain unknown. This cross-sectional study characterizes alternative autoimmunity, and associated genetic and environmental factors, in unaffected first-degree relatives (FDRs) of patients, who exhibit increased future risk for the same disease.

METHODS

Samples (n = 1321) from disease-specific autoantibody-positive (aAb+) systemic lupus erythematosus (SLE), rheumatoid arthritis (RA), and type 1 diabetes (T1D) patients; and unaffected aAb+ and autoantibody-negative (aAb-) SLE and RA FDRs were tested for SLE, RA, and T1D aAbs, as well as anti-tissue transglutaminase, anti-cardiolipin and anti-thyroperoxidase. FDR SLE and RA genetic risk scores (GRS) were calculated.

FINDINGS

Alternative autoimmunity occurred in SLE patients (56%) and FDRs (57·4%), RA patients (32·6%) and FDRs (34·8%), and T1D patients (43%). Expanded autoimmunity, defined as autoantibodies spanning at least two other diseases, occurred in 18·5% of SLE patients, 16·4% of SLE FDRs, 7·8% of RA patients, 5·3% of RA FDRs, and 10·8% of T1D patients. SLE FDRs were more likely to have alternative (odds ratio [OR] 2·44) and expanded (OR 3·27) autoimmunity than RA FDRs. Alternative and expanded autoimmunity were associated with several environmental exposures. Alternative autoimmunity was associated with a higher RA GRS in RA FDRs (OR 1·41), and a higher SLE GRS in aAb+ RA FDRs (OR 1·87), but not in SLE FDRs.

INTERPRETATION

Autoimmunity commonly crosses disease-specific boundaries in systemic (RA, SLE) and organ-specific (T1D) autoimmune diseases. Alternative autoimmunity is more common in SLE FDRs than RA FDRs, and is influenced by genetic and environmental factors. These findings have substantial implications for preclinical disease pathogenesis and autoimmune disease prevention studies. FUND: NIH U01AI101981, R01AR051394, U19AI082714, P30AR053483, P30GM103510, U54GM104938, U01AI101934, R01AI024717, U01AI130830, I01BX001834, & U01HG008666.

Screening characteristics for enrichment of individuals at higher risk for transitioning to classified SLE

OBJECTIVE

Further prospective study is needed to elucidate the etiology and natural history of systemic lupus erythematosus development. The clinical complexity of this heterogeneous disease makes study design challenging. Our objective was to ascertain useful screening factors for identifying at-risk individuals for follow-up rheumatologic assessment or inclusion in prospective studies.

METHODS

We attempted to re-contact 3823 subjects with a family history of systemic lupus erythematosus, who did not meet American College of Rheumatology systemic lupus erythematosus classification at a baseline study visit; 436 agreed to follow-up participation an average of 6.3 years after baseline. In total, 56 of these individuals had transitioned to classified systemic lupus erythematosus (≥ 4 cumulative American College of Rheumatology criteria, verified by medical record review) by the time of follow up. Generalized estimating equations assessed associations between our dichotomous outcome of transitioning to systemic lupus erythematosus with baseline characteristics, including ANA positivity, Connective Tissue Disease Screening questionnaire systemic lupus erythematosus score, and number of American College of Rheumatology criteria. We analyzed predictive accuracy of characteristics on transitioning.

RESULTS

ANA positivity, Connective Tissue Disease Screening questionnaire systemic lupus erythematosus score categorization of possible or probable systemic lupus erythematosus, and greater number of American College of Rheumatology criteria at baseline were each associated with transitioning to systemic lupus erythematosus classification. Being ANA positive and having confirmed immunologic criteria at baseline had the highest positive predictive value and specificity for transitioning to systemic lupus erythematosus. American College of Rheumatology Connective Tissue Disease Screening questionnaire systemic lupus erythematosus score categorization of possible or probable systemic lupus erythematosus had a better positive predictive value, negative predictive value, sensitivity, and specificity than ANA positivity.

CONCLUSION

Given limited resources, identifying individuals for follow up based on the systemic lupus erythematosus portion of the Connective Tissue Disease Screening questionnaire could be an efficient way to identify family members at highest risk of disease transition.

Novel genetic associations with interferon in systemic lupus erythematosus identified by replication and fine-mapping of trait-stratified genome-wide screen

BACKGROUND/PURPOSE

High serum interferon alpha (IFN-α) is an important heritable phenotype in systemic lupus erythematosus (SLE) which is involved in primary disease pathogenesis. High vs. low levels of IFN-α are associated with disease severity and account for some of the biological heterogeneity between SLE patients. The aim of the study was to replicate and fine-map previously detected genetic associations with serum IFN-α in SLE.

METHODS

We previously undertook a case-case genome-wide association study of SLE patients stratified by ancestry and extremes of phenotype in serum IFN-α. Single nucleotide polymorphisms (SNPs) in seven loci identified in this screen were selected for follow up in a large independent cohort of 1370 SLE patients (703 European-ancestry, 432 African ancestry, and 235 Amerindian ancestry). Each ancestral background was analyzed separately, and ancestry-informative markers were used to control for ancestry and admixture.

RESULTS

We find a rare haplotype spanning the promoter region of EFNA5 that is strongly associated with serum IFN-α in both African-American and European-American SLE patients (OR = 3.0, p = 3.7 × 10^-6). We also find SNPs in the PPM1H, PTPRM, and NRGN regions associated with IFN-α levels in European-American, Amerindian, and African-American SLE patients respectively. Many of these associations are within regulatory regions of the gene, suggesting an impact on transcription.

CONCLUSION

This study demonstrates the power of molecular sub-phenotypes to reveal genetic factors involved in complex autoimmune disease. The distinct associations observed in different ancestral backgrounds emphasize the heterogeneity of molecular pathogenesis in SLE.

Trans-Ethnic Mapping of BANK1 Identifies Two Independent SLE-Risk Linkage Groups Enriched for Co-Transcriptional Splicing Marks

BANK1 is a susceptibility gene for several systemic autoimmune diseases in several populations. Using the genome-wide association study (GWAS) data from Europeans (EUR) and African Americans (AA), we performed an extensive fine mapping of ankyrin repeats 1 (BANK1). To increase the SNP density, we used imputation followed by univariate and conditional analysis, combined with a haplotypic and expression quantitative trait locus (eQTL) analysis. The data from Europeans showed that the associated region was restricted to a minimal and dependent set of SNPs covering introns two and three, and exon two. In AA, the signal found in the Europeans was split into two independent effects. All of the major risk associated SNPs were eQTLs, and the risks were associated with an increased BANK1 gene expression. Functional annotation analysis revealed the enrichment of repressive B cell epigenomic marks (EZH2 and H3K27me3) and a strong enrichment of splice junctions. Furthermore, one eQTL located in intron two, rs13106926, was found within the binding site for RUNX3, a transcriptional activator. These results connect the local genome topography, chromatin structure, and the regulatory landscape of BANK1 with co-transcriptional splicing of exon two. Our data defines a minimal set of risk associated eQTLs predicted to be involved in the expression of BANK1 modulated through epigenetic regulation and splicing. These findings allow us to suggest that the increased expression of BANK1 will have an impact on B-cell mediated disease pathways.

A plausibly causal functional lupus-associated risk variant in the STAT1-STAT4 locus

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.