Confirmation of five novel susceptibility loci for systemic lupus erythematosus (SLE) and integrated network analysis of 82 SLE susceptibility loci

We recently identified ten novel SLE susceptibility loci in Asians and uncovered several additional suggestive loci requiring further validation. This study aimed to replicate five of these suggestive loci in a Han Chinese cohort from Hong Kong, followed by meta-analysis (11,656 cases and 23,968 controls) on previously reported Asian and European populations, and to perform bioinformatic analyses on all 82 reported SLE loci to identify shared regulatory signatures. We performed a battery of analyses for these five loci, as well as joint analyses on all 82 SLE loci. All five loci passed genome-wide significance: MYNN (rs10936599, Pmeta = 1.92 × 10-13, OR = 1.14), ATG16L2 (rs11235604, Pmeta = 8.87 × 10 -12, OR = 0.78), CCL22 (rs223881, Pmeta = 5.87 × 10-16, OR = 0.87), ANKS1A (rs2762340, Pmeta = 4.93 × 10-15, OR = 0.87) and RNASEH2C (rs1308020, Pmeta = 2.96 × 10-19, OR = 0.84) and co-located with annotated gene regulatory elements. The novel loci share genetic signatures with other reported SLE loci, including effects on gene expression, transcription factor binding, and epigenetic characteristics. Most (56%) of the correlated (r2 > 0.8) SNPs from the 82 SLE loci were implicated in differential expression (9.81 × 10-198 < P < 5 × 10-3) of cis-genes. Transcription factor binding sites for p53, MEF2A and E2F1 were significantly (P < 0.05) over-represented in SLE loci, consistent with apoptosis playing a critical role in SLE. Enrichment analysis revealed common pathways, gene ontology, protein domains, and cell type-specific expression. In summary, we provide evidence of five novel SLE susceptibility loci. Integrated bioinformatics using all 82 loci revealed that SLE susceptibility loci share many gene regulatory features, suggestive of conserved mechanisms of SLE etiopathogenesis.

[1]  J. Park,et al.  Replication of the results of genome-wide and candidate gene association studies on telomere length in a Korean population , 2015, The Korean journal of internal medicine.

[2]  Liangdan Sun,et al.  Whole-exome SNP array identifies 15 new susceptibility loci for psoriasis , 2015, Nature Communications.

[3]  Y. J. Kim,et al.  High-density genotyping of immune-related loci identifies new SLE risk variants in individuals with Asian ancestry , 2016, Nature Genetics.

[4]  R. Melamed,et al.  Gene expression microarrays: glimpses of the immunological genome , 2006, Nature Immunology.

[5]  Dmitri D. Pervouchine,et al.  The human transcriptome across tissues and individuals , 2015, Science.

[6]  Kai Wang,et al.  Enlight: web-based integration of GWAS results with biological annotations , 2015, Bioinform..

[7]  Hedi Peterson,et al.  g:Profiler—a web server for functional interpretation of gene lists (2016 update) , 2016, Nucleic Acids Res..

[8]  M. Peters,et al.  Systematic identification of trans eQTLs as putative drivers of known disease associations , 2013, Nature Genetics.

[9]  Yusun Jung,et al.  Genome-wide association study of Crohn's disease in Koreans revealed three new susceptibility loci and common attributes of genetic susceptibility across ethnic populations , 2013, Gut.

[10]  J. Kere,et al.  Genes identified in Asian SLE GWASs are also associated with SLE in Caucasian populations , 2012, European Journal of Human Genetics.

[11]  E. Papoutsakis,et al.  A global transcriptional view of apoptosis in human T-cell activation , 2008, BMC Medical Genomics.

[12]  Å. Lernmark,et al.  Identification of Non-HLA Genes Associated with Celiac Disease and Country-Specific Differences in a Large, International Pediatric Cohort , 2016, PloS one.

[13]  P. Visscher,et al.  Conditional and joint multiple-SNP analysis of GWAS summary statistics identifies additional variants influencing complex traits , 2012, Nature Genetics.

[14]  L. Lagae,et al.  Characterization of human disease phenotypes associated with mutations in TREX1, RNASEH2A, RNASEH2B, RNASEH2C, SAMHD1, ADAR, and IFIH1 , 2015, American journal of medical genetics. Part A.

[15]  K. Korotayev,et al.  Many pathways to apoptosis: E2F1 regulates splicing of apoptotic genes , 2008, Cell Death and Differentiation.

[16]  P. Deloukas,et al.  Multiple common variants for celiac disease influencing immune gene expression , 2010, Nature Genetics.

[17]  Cory Y. McLean,et al.  GREAT improves functional interpretation of cis-regulatory regions , 2010, Nature Biotechnology.

[18]  Troels Z. Kristiansen,et al.  Cloning of a novel phosphotyrosine binding domain containing molecule, Odin, involved in signaling by receptor tyrosine kinases , 2002, Oncogene.

[19]  C. Glass,et al.  Simple combinations of lineage-determining transcription factors prime cis-regulatory elements required for macrophage and B cell identities. , 2010, Molecular cell.

[20]  L. Looger,et al.  Combined protein- and nucleic acid-level effects of rs1143679 (R77H), a lupus-predisposing variant within ITGAM. , 2014, Human molecular genetics.

[21]  S. Meshaal,et al.  Signal transducer and activator of transcription 5 is implicated in disease activity in adult and juvenile onset systemic lupus erythematosus , 2016, Clinical Rheumatology.

[22]  Lu,et al.  Genome-wide association meta-analysis in Chinese and European individuals identifies ten new loci associated with systemic lupus erythematosus , 2016, Nature Genetics.

[23]  Simon G. Coetzee,et al.  motifbreakR: an R/Bioconductor package for predicting variant effects at transcription factor binding sites , 2015, Bioinform..

[24]  P. Sham,et al.  Genome-Wide Association Study in Asian Populations Identifies Variants in ETS1 and WDFY4 Associated with Systemic Lupus Erythematosus , 2010, PLoS genetics.

[25]  Anushya Muruganujan,et al.  PANTHER version 10: expanded protein families and functions, and analysis tools , 2015, Nucleic Acids Res..

[26]  A. Skol,et al.  Trait-stratified genome-wide association study identifies novel and diverse genetic associations with serologic and cytokine phenotypes in systemic lupus erythematosus , 2010, Arthritis research & therapy.

[27]  Avi Ma'ayan,et al.  Enrichr: interactive and collaborative HTML5 gene list enrichment analysis tool , 2013, BMC Bioinformatics.

[28]  Xuejun Zhang,et al.  Several Critical Cell Types, Tissues, and Pathways Are Implicated in Genome-Wide Association Studies for Systemic Lupus Erythematosus , 2016, G3: Genes, Genomes, Genetics.

[29]  Swneke D. Bailey,et al.  Breast cancer risk-associated SNPs modulate the affinity of chromatin for FOXA1 and alter gene expression , 2012, Nature Genetics.

[30]  Ralf Herwig,et al.  ConsensusPathDB—a database for integrating human functional interaction networks , 2008, Nucleic Acids Res..

[31]  C. Weston,et al.  Distinct Roles for CCR4 and CXCR3 in the Recruitment and Positioning of Regulatory T Cells in the Inflamed Human Liver , 2010, The Journal of Immunology.

[32]  F. Alt,et al.  Increased T-cell apoptosis and terminal B-cell differentiation induced by inactivation of the Ets-1 proto-oncogene , 1995, Nature.

[33]  Xinli Hu,et al.  SNPsea: an algorithm to identify cell types, tissues and pathways affected by risk loci , 2014, Bioinform..

[34]  J. Rioux,et al.  Genetic association analyses implicate aberrant regulation of innate and adaptive immunity genes in the pathogenesis of systemic lupus erythematosus , 2015, Nature Genetics.

[35]  M. Mohammadi,et al.  Higher circulating levels of chemokines CXCL10, CCL20 and CCL22 in patients with ischemic heart disease. , 2016, Cytokine.

[36]  N. Bresolin,et al.  MDC/CCL22 intrathecal levels in patients with multiple sclerosis , 2008, Multiple sclerosis.

[37]  U. Pati,et al.  Over‐expression of TATA binding protein (TBP) and p53 and autoantibodies to these antigens are features of systemic sclerosis, systemic lupus erythematosus and overlap syndromes , 2004, Clinical and experimental immunology.

[38]  Yun Li,et al.  METAL: fast and efficient meta-analysis of genomewide association scans , 2010, Bioinform..

[39]  S. Batalov,et al.  A gene atlas of the mouse and human protein-encoding transcriptomes. , 2004, Proceedings of the National Academy of Sciences of the United States of America.

[40]  N. Minato,et al.  Rap1 Functions as a Key Regulator of T-Cell and Antigen-Presenting Cell Interactions and Modulates T-Cell Responses , 2002, Molecular and Cellular Biology.

[41]  Ryan T. Phan,et al.  The BCL6 proto-oncogene suppresses p53 expression in germinal-centre B cells , 2004, Nature.

[42]  M. Hochberg,et al.  Updating the American College of Rheumatology revised criteria for the classification of systemic lupus erythematosus. , 1997, Arthritis and rheumatism.

[43]  Esra Bozgeyik,et al.  Gene expression profiles of autophagy-related genes in multiple sclerosis. , 2016, Gene.

[44]  G. Tsokos,et al.  Abnormal NF-kappa B activity in T lymphocytes from patients with systemic lupus erythematosus is associated with decreased p65-RelA protein expression. , 1999, Journal of immunology.

[45]  Scott M. Williams,et al.  Admixture Mapping in Lupus Identifies Multiple Functional Variants within IFIH1 Associated with Apoptosis, Inflammation, and Autoantibody Production , 2013, PLoS genetics.

[46]  C. Ponting,et al.  Mutations in genes encoding ribonuclease H2 subunits cause Aicardi-Goutières syndrome and mimic congenital viral brain infection , 2006, Nature Genetics.

[47]  T. Chan,et al.  The role of microRNA-1246 in the regulation of B cell activation and the pathogenesis of systemic lupus erythematosus , 2015, Clinical Epigenetics.

[48]  Takuro Nakamura,et al.  STAT5A Modulates Chemokine Receptor CCR6 Expression and Enhances Pre‐B Cell Growth in a CCL20‐Dependent Manner , 2016, Journal of cellular biochemistry.

[49]  Aneela Majid,et al.  A genome-wide association study identifies multiple susceptibility loci for chronic lymphocytic leukemia , 2013, Nature Genetics.

[50]  S. Feller,et al.  Odin (ANKS1A) is a Src family kinase target in colorectal cancer cells , 2008, Cell Communication and Signaling.

[51]  Ellen T. Gelfand,et al.  The Genotype-Tissue Expression (GTEx) project , 2013, Nature Genetics.

[52]  Marta E Alarcón-Riquelme,et al.  Genome-wide association scan in women with systemic lupus erythematosus identifies susceptibility variants in ITGAM, PXK, KIAA1542 and other loci , 2008, Nature Genetics.

[53]  C. Vestergaard,et al.  Expression of MDC/CCL22 and its receptor CCR4 in rheumatoid arthritis, psoriatic arthritis and osteoarthritis. , 2010, Cytokine.

[54]  Manuel A. R. Ferreira,et al.  PLINK: a tool set for whole-genome association and population-based linkage analyses. , 2007, American journal of human genetics.

[55]  A. Mahajan,et al.  Clearance Deficiency and Cell Death Pathways: A Model for the Pathogenesis of SLE , 2016, Front. Immunol..

[56]  Y. Okada,et al.  A Genome-Wide Association Study Identified AFF1 as a Susceptibility Locus for Systemic Lupus Eyrthematosus in Japanese , 2012, PLoS genetics.

[57]  T. Avčin,et al.  The Role of STAT Signaling Pathways in the Pathogenesis of Systemic Lupus Erythematosus , 2017, Clinical Reviews in Allergy & Immunology.

[58]  P. Alliel,et al.  Neuromuscular expression of the BTB/POZ and zinc finger protein myoneurin , 2004, Muscle & nerve.

[59]  William Stafford Noble,et al.  Sequence features and chromatin structure around the genomic regions bound by 119 human transcription factors , 2012, Genome research.

[60]  P. Sham,et al.  Evaluating the heritability explained by known susceptibility variants: a survey of ten complex diseases , 2011, Genetic epidemiology.

[61]  Yang Du,et al.  rSNPBase: a database for curated regulatory SNPs , 2013, Nucleic Acids Res..

[62]  Manolis Kellis,et al.  HaploReg: a resource for exploring chromatin states, conservation, and regulatory motif alterations within sets of genetically linked variants , 2011, Nucleic Acids Res..

[63]  P. Gaffney,et al.  Identification of a Systemic Lupus Erythematosus Risk Locus Spanning ATG16L2, FCHSD2, and P2RY2 in Koreans , 2015, Arthritis & rheumatology.

[64]  Vladimir B. Bajic,et al.  HOCOMOCO: a comprehensive collection of human transcription factor binding sites models , 2012, Nucleic Acids Res..

[65]  Midori Yoshida,et al.  Overexpression of PU.1 induces growth and differentiation inhibition and apoptotic cell death in murine erythroleukemia cells. , 1997, Blood.

[66]  M. Petri,et al.  Baseline characteristics of a multiethnic lupus cohort: PROFILE , 2002, Lupus.

[67]  Derek W Wright,et al.  Gateways to the FANTOM5 promoter level mammalian expression atlas , 2015, Genome Biology.

[68]  Kaanan P. Shah,et al.  A gene-based association method for mapping traits using reference transcriptome data , 2015, Nature Genetics.

[69]  P. Allavena,et al.  Human Macrophage–derived Chemokine (mdc), a Novel Chemoattractant for Monocytes, Monocyte-derived Dendritic Cells, and Natural Killer Cells , 1997 .

[70]  Kyong-Ah Yoon,et al.  Prognostic implications of genetic variants in advanced non-small cell lung cancer: a genome-wide association study. , 2013, Carcinogenesis.