Identification of a two-loci epistatic interaction associated with susceptibility to rheumatoid arthritis through reverse engineering and multifactor dimensionality reduction.

Altered synovial fibroblast (SF) transcriptional activity is a key factor in the disease progression of rheumatoid arthritis (RA). To determine the transcriptional regulatory network associated with SF response to an RA proinflammatory stimulus we applied a CARRIE reverse engineering approach to microarray gene expression data from SFs treated with RA synovial fluid. The association of the inferred gene network with RA susceptibility was further analyzed by a case-control study of promoter single-nucleotide polymorphisms, and the presence of epistatic interactions was determined using the multifactor dimensionality reduction methodology. Our findings suggest that a specific NF-kappaB transcriptional regulatory network of 13 genes is associated with SF response to RA proinflammatory stimulus and identify a significant epistatic association of two of its genes, IL6 and IL4I1, with RA susceptibility.

[1]  David B Allison,et al.  Novel tumor necrosis factor alpha-regulated genes in rheumatoid arthritis. , 2004, Arthritis and rheumatism.

[2]  R. Kastelein,et al.  Two novel IL-1 family members, IL-1 delta and IL-1 epsilon, function as an antagonist and agonist of NF-kappa B activation through the orphan IL-1 receptor-related protein 2. , 2001, Journal of immunology.

[3]  Z. Zhai,et al.  SINK Is a p65-interacting Negative Regulator of NF-κB-dependent Transcription* , 2003, Journal of Biological Chemistry.

[4]  M. Pierer,et al.  Chemokine Secretion of Rheumatoid Arthritis Synovial Fibroblasts Stimulated by Toll-Like Receptor 2 Ligands1 , 2004, The Journal of Immunology.

[5]  T. Mcclanahan,et al.  Divergent Pro- and Antiinflammatory Roles for IL-23 and IL-12 in Joint Autoimmune Inflammation , 2003, The Journal of experimental medicine.

[6]  A. Hoffmann,et al.  Regulation of IL-6 and IL-8 Expression in Rheumatoid Arthritis Synovial Fibroblasts: the Dominant Role for NF-κB But Not C/EBPβ or c-Jun1 , 2000, The Journal of Immunology.

[7]  R. Kastelein,et al.  Two Novel IL-1 Family Members, IL-1δ and IL-1ε, Function as an Antagonist and Agonist of NF-κB Activation Through the Orphan IL-1 Receptor-Related Protein 21 , 2001, The Journal of Immunology.

[8]  A. Julià,et al.  Lack of association between the corticotropin-releasing hormone locus and rheumatoid arthritis. , 2004, Arthritis and rheumatism.

[9]  C. Gabay,et al.  The new IL-1 family member IL-1F8 stimulates production of inflammatory mediators by synovial fibroblasts and articular chondrocytes , 2006, Arthritis research & therapy.

[10]  J. H. Moore,et al.  Multifactor-dimensionality reduction shows a two-locus interaction associated with Type 2 diabetes mellitus , 2004, Diabetologia.

[11]  K. Hiromura,et al.  Activin A induces cell proliferation of fibroblast-like synoviocytes in rheumatoid arthritis. , 2003, Arthritis and rheumatism.

[12]  Kutty Selva Nandakumar,et al.  Identification of epistasis through a partial advanced intercross reveals three arthritis loci within the Cia5 QTL in mice , 2005, Genes and Immunity.

[13]  Adam A. Margolin,et al.  Reverse engineering of regulatory networks in human B cells , 2005, Nature Genetics.

[14]  A. Balsa,et al.  IL-6 promoter polymorphisms in rheumatoid arthritis , 2000, Genes and Immunity.

[15]  Jason H. Moore,et al.  A global view of epistasis , 2005, Nature Genetics.

[16]  Ccggaattccatcacgtgttcgtggagact Tgtaaaacgacggccagtgc,et al.  Second Class Minors : Molecular Identification of the Autosomal H 46 Histocompatibility Locus as a Peptide Presented by Major Histocompatibility Complex Class II Molecules , 2003 .

[17]  T. Libermann,et al.  Activation of interleukin-6 gene expression through the NF-kappa B transcription factor , 1990, Molecular and cellular biology.

[18]  Jason H. Moore,et al.  Multifactor dimensionality reduction software for detecting gene-gene and gene-environment interactions , 2003, Bioinform..

[19]  E. G. de la Concha,et al.  Epistatic interaction between FCRL3 and NFκB1 genes in Spanish patients with rheumatoid arthritis , 2006, Annals of the rheumatic diseases.

[20]  S. Gabriel,et al.  The Structure of Haplotype Blocks in the Human Genome , 2002, Science.

[21]  G. Firestein.  NF-kappaB: Holy Grail for rheumatoid arthritis? , 2004, Arthritis and rheumatism.

[22]  C. Mackay,et al.  Overlapping gene expression profiles in rheumatoid fibroblast-like synoviocytes induced by the proinflammatory cytokines interleukin-1 β and tumor necrosis factor , 2004, Inflammation Research.

[23]  G. Firestein.  Evolving concepts of rheumatoid arthritis , 2003, Nature.

[24]  M. F. Shannon,et al.  Requirement for nuclear factor (NF)-kappa B p65 and NF-interleukin-6 binding elements in the tumor necrosis factor response region of the granulocyte colony-stimulating factor promoter. , 1994, Blood.

[25]  G. Firestein.  NF‐κB: Holy Grail for rheumatoid arthritis? , 2004 .

[26]  P. Sham,et al.  Model-Free Analysis and Permutation Tests for Allelic Associations , 1999, Human Heredity.

[27]  S. Haskill,et al.  Identification of three related human GRO genes encoding cytokine functions. , 1990, Proceedings of the National Academy of Sciences of the United States of America.

[28]  Todd Holden,et al.  A flexible computational framework for detecting, characterizing, and interpreting statistical patterns of epistasis in genetic studies of human disease susceptibility. , 2006, Journal of theoretical biology.

[29]  P. Tak,et al.  NF-kappaB: a key role in inflammatory diseases. , 2001, The Journal of clinical investigation.

[30]  A. Silman,et al.  Epistatic modeling in rheumatoid arthritis: An application of the Risch theory , 1993, Genetic epidemiology.

[31]  H. Makino,et al.  The synovial expression and serum levels of interleukin-6, interleukin-11, leukemia inhibitory factor, and oncostatin M in rheumatoid arthritis. , 1997, Arthritis and rheumatism.

[32]  T. Sumida,et al.  Direct evidence of high DNA binding activity of transcription factor AP-1 in rheumatoid arthritis synovium. , 1997, Arthritis and rheumatism.

[33]  P. Tak,et al.  NF-κB: a key role in inflammatory diseases , 2001 .

[34]  David B. Allison,et al.  Novel tumor necrosis factor α-regulated genes in rheumatoid arthritis , 2004 .

[35]  R. Holmdahl,et al.  Two-loci interaction confirms arthritis-regulating quantitative trait locus on rat chromosome 6. , 2003, Genomics.

[36]  J. Backman,et al.  Expression of cyclooxygenase 1 and cyclooxygenase 2 in human synovial tissue: differential elevation of cyclooxygenase 2 in inflammatory joint diseases. , 1998, Arthritis & Rheumatism.

[37]  N. Ito,et al.  Keratinocyte‐specific modulation of type VII collagen gene expression by pro‐inflammatory cytokines (tumor necrosis factor‐α and interleukin‐1β) , 2005, Experimental dermatology.

[38]  R. Gay,et al.  Synovial fibroblasts: key players in rheumatoid arthritis. , 2006, Rheumatology.

[39]  N E Morton,et al.  Tests and estimates of allelic association in complex inheritance. , 1998, Proceedings of the National Academy of Sciences of the United States of America.

[40]  P. Tak,et al.  Fibroblast-like synoviocytes derived from patients with rheumatoid arthritis show the imprint of synovial tissue heterogeneity: evidence of a link between an increased myofibroblast-like phenotype and high-inflammation synovitis. , 2005, Arthritis and rheumatism.

[41]  G. Church,et al.  Modular epistasis in yeast metabolism , 2005, Nature Genetics.

[42]  A. Manning,et al.  NF-κB Regulation by IκB Kinase-2 in Rheumatoid Arthritis Synoviocytes1 , 2001, The Journal of Immunology.

[43]  G. Rubin,et al.  Generation and initial analysis of more than 15,000 full-length human and mouse cDNA sequences , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[44]  L. Davis.  A question of transformation: the synovial fibroblast in rheumatoid arthritis. , 2003, The American journal of pathology.

[45]  S. Makarov,et al.  An essential role of NF-κB in the “tumor-like” phenotype of arthritic synoviocytes , 2006, Proceedings of the National Academy of Sciences.

[46]  S. Akira,et al.  The role of NF-kappa B and NF-IL6 transactivating factors in the synergistic activation of human serum amyloid A gene expression by interleukin-1 and interleukin-6. , 1993, The Journal of biological chemistry.

[47]  A. Manning,et al.  NF-kappa B regulation by I kappa B kinase-2 in rheumatoid arthritis synoviocytes. , 2001, Journal of immunology.

[48]  D. Green,et al.  Regional analysis of p53 mutations in rheumatoid arthritis synovium , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[49]  J. H. Moore,et al.  Multifactor-dimensionality reduction reveals high-order interactions among estrogen-metabolism genes in sporadic breast cancer. , 2001, American journal of human genetics.

[50]  D. Kwiatkowski,et al.  The effect of HLA-DR on susceptibility to rheumatoid arthritis is influenced by the associated lymphotoxin alpha-tumor necrosis factor haplotype. , 2003, Arthritis and rheumatism.

[51]  Z. Zhai,et al.  SINK is a p65-interacting negative regulator of NF-kappaB-dependent transcription. , 2003, The Journal of biological chemistry.

[52]  A. Hartemink.  Reverse engineering gene regulatory networks , 2005, Nature Biotechnology.

[53]  C. Copie-Bergman,et al.  Interleukin 4-induced gene 1 is activated in primary mediastinal large B-cell lymphoma. , 2003, Blood.

[54]  A. Hoffmann,et al.  Regulation of Il-6 and Il-8 Expression In , 2013 .

[55]  Zhiping Weng,et al.  Computational inference of transcriptional regulatory networks from expression profiling and transcription factor binding site identification. , 2004, Nucleic acids research.