Dataset Information

A towards-multidimensional screening approach to predict candidate genes of rheumatoid arthritis based on SNP, structural and functional annotations.

ABSTRACT:

Background

According to the Genetic Analysis Workshops (GAW), hundreds of thousands of SNPs have been tested for association with rheumatoid arthritis. Traditional genome-wide association studies (GWAS) have been developed to identify susceptibility genes using a "most significant SNPs/genes" model. However, many minor- or modest-risk genes are likely to be missed after adjustment of multiple testing. This screening process uses a strict selection of statistical thresholds that aim to identify susceptibility genes based only on statistical model, without considering multi-dimensional biological similarities in sequence arrangement, crystal structure, or functional categories/biological pathways between candidate and known disease genes.

Methods

Multidimensional screening approaches combined with traditional statistical genetics methods can consider multiple biological backgrounds of genetic mutation, structural, and functional annotations. Here we introduce a newly developed multidimensional screening approach for rheumatoid arthritis candidate genes that considers all SNPs with nominal evidence of Bayesian association (BFLn > 0), and structural and functional similarities of corresponding genes or proteins.

Results

Our multidimensional screening approach extracted all risk genes (BFLn > 0) by odd ratios of hypothesis H1 to H0, and determined whether a particular group of genes shared underlying biological similarities with known disease genes. Using this method, we found 6614 risk SNPs in our Bayesian screen result set. Finally, we identified 146 likely causal genes for rheumatoid arthritis, including CD4, FGFR1, and KDR, which have been reported as high risk factors by recent studies. We must denote that 790 (96.1%) of genes identified by GWAS could not easily be classified into related functional categories or biological processes associated with the disease, while our candidate genes shared underlying biological similarities (e.g. were in the same pathway or GO term) and contributed to disease etiology, but where common variations in each of these genes make modest contributions to disease risk. We also found 6141 risk SNPs that were too minor to be detected by conventional approaches, and associations between 58 candidate genes and rheumatoid arthritis were verified by literature retrieved from the NCBI PubMed module.

Conclusions

Our proposed approach to the analysis of GAW16 data for rheumatoid arthritis was based on an underlying biological similarities-based method applied to candidate and known disease genes. Application of our method could identify likely causal candidate disease genes of rheumatoid arthritis, and could yield biological insights that not detected when focusing only on genes that give the strongest evidence by multiple testing. We hope that our proposed method complements the "most significant SNPs/genes" model, and provides additional insights into the pathogenesis of rheumatoid arthritis and other diseases, when searching datasets for hundreds of genetic variances.

SUBMITTER: Zhang L

PROVIDER: S-EPMC2939610 | biostudies-literature | 2010 Aug

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

A towards-multidimensional screening approach to predict candidate genes of rheumatoid arthritis based on SNP, structural and functional annotations.

Zhang Liangcai L Li Wan W Song Leilei L Chen Lina L

BMC medical genomics 20100820

<h4>Background</h4>According to the Genetic Analysis Workshops (GAW), hundreds of thousands of SNPs have been tested for association with rheumatoid arthritis. Traditional genome-wide association studies (GWAS) have been developed to identify susceptibility genes using a "most significant SNPs/genes" model. However, many minor- or modest-risk genes are likely to be missed after adjustment of multiple testing. This screening process uses a strict selection of statistical thresholds that aim to id ...[more]

PMID: 20727150

Similar Datasets

Project description:Rheumatoid arthritis (RA), which normally manifests as a multi‑joint inflammatory reaction, is a common immunological disease in clinical practice. However, the pathogenesis of RA has not yet been fully elucidated. Rituximab (RTX) is an effective drug in the treatment of RA, however its therapeutic efficacy and mechanism of action require further investigation. Thus, the present study aimed to screen the candidate key regulatory genes and explain the potential mechanisms of RA. Gene chips of RA and normal joint tissues were analyzed and, gene chips of RTX before and after treatment were investigated. In the present study, strong evidence supporting the pathogenesis of RA and mechanism of action of RTX were also revealed. Differentially expressed genes (DEGs) were analyzed using the limma package of RStudio software. A total of 1,150 DEGs were detected in RA compared with normal joint tissues. The upregulated genes were enriched in 'interleukin‑12 production', 'I‑κB kinase/NF‑κB signaling', 'regulation of cytokine production involved in immune response' and 'cytokine metabolic process'. Functional enrichment analysis showed that RTX was primarily involved in the inhibition of 'adaptive immune response', 'B cell activation involved in immune response' and 'immune effector process'. Subsequently, leukocyte immunoglobulin‑like receptor subfamily B member 1 (LILRB1), a hub gene with high connectivity degree, was selected, and traditional Chinese medicine libraries were molecularly screened according to the structure of the LILRB1 protein. The results indicated that kaempferol 3‑O‑β‑D‑glucosyl‑(1→2)‑β‑D‑glucoside exhibited the highest docking score. In the present study, the DEGs and their biological functions in RA and the pharmacological mechanism of RTX action were determined. Taken together, the results suggested that LILRB1 may be used as a molecular target for RA treatment, and kaempferol 3‑O‑β‑D‑glucosyl‑(1→2)‑β‑D‑glucoside may inhibit the pathological process of RA.

Project description:Genotyping was performed using a custom-designed Illumina 384-SNP VeraCode microarray (Illumina) to determine possible associations of genes encoding PADI1, PADI2, PADI3, PADI4, PADI6, PRKD3, Gc, GLRX,CDSN, PSORS1C1, TXNDC5, CA1, bub3, SORBS1, VDR,SERPNA1, PCSK6, DNAH9 to rheumatoid arthrities (RA) and ankylosing spondylitis (AS). We previously found that encoding genes of PADI1, PADI2, PADI3, PADI4, PADI6, PRKD3, Gc, GLRX,CDSN, PSORS1C1, TXNDC5, CA1, bub3, SORBS1, VDR, SERPNA1, PCSK6 and DNAH9 could be associated with rheumatoid arthritis (RA) and ankylosing spondylitis (AS). We used illumina customer designed microarray to verify the above finding with an increased numbers of samples. Genotyping was performed using a custom-designed Illumina 384-SNP VeraCode microarray (Illumina) to determine possible associations of genes encoding PADI1, PADI2, PADI3, PADI4, PADI6, PRKD3, Gc, GLRX,CDSN, PSORS1C1, TXNDC5, CA1, bub3, SORBS1, VDR, SERPNA1, PCSK6, DNAH9 to RA and AS. A case-control study was performed with 51 AS patients, 266 rheumatoid arthritis (RA) and 163 health controls. Tag single nucleotide polymorphisms (tag SNPs) across the above genes were determined by searching the HapMap database. The tag SNPs were selected on the basis of linkage disequilibrium patterns observed in the Han Chinese in Beijing (CHB) samples that were genotyped as a part of the International HapMap Project. Only SNPs with a minor allele frequency of greater than 5% with a pair-wise r2 ≥ 0.8 in the HapMap were considered. SNPs that were in exons and untranslated regions (UTR), as well as promoters and introns within 500 bp of exons were also selected. Candidate SNPs were submitted to Illumina for a design score. The Illumina Assay Design Tool filtered out SNPs that were not suitable for the Illumina platform. Finally, 382 SNPs with a design score of “1” were selected. 382 SNPs within PADI1, PADI2, PADI3, PADI4, PADI6, PRKD3, Gc, GLRX,CDSN, PSORS1C1, TXNDC5, CA1, bub3, SORBS1, VDR,SERPNA1, PCSK6, DNAH9 encoding genes were selected. Genotyping was performed using Illumina VeraCode microarray in a case-control study including 51 AS patients, 266 rheumatoid arthritis (RA) and 163 health controls We performed genotyping using a custom-designed Illumina 384-SNP VeraCode microarray (Illumina). Peripheral blood samples were collected from patients with RA (n=266, 183 female) and AS (n=51, 10 female). RA patients had a mean age of 51.7 years, while AS patients had a mean age of 35.9 years. A total of 163 (60 female) healthy individuals with a mean age of 48.0 years were blood donors. Blood samples were put into Monovette tubes containing 3.8% sodium citrate. The genotyping was conducted with the BeadXpress Reader using the Illumina VeraCode GoldenGate Assay kit. 500 ng of sample DNA were used per assay. Genotype clustering and calling were performed using BeadStudio software (Genomestudio V2010.2 Genotyping module V1.6) (Illumina). This work was completed at the Beijing Institute of Genomics.

Project description:Rheumatoid arthritis (RA) is a chronic autoimmune, inflammatory disease, characterized by synovitis in small- and medium-sized joints and, if not treated early and efficiently, joint damage, and destruction. RA is a heterogeneous disease with a plethora of treatment options. The pro-inflammatory cytokine tumor necrosis factor (TNF) plays a central role in the pathogenesis of RA, and TNF inhibitors effectively repress inflammatory activity in RA. Currently, treatment decisions are primarily based on empirics and economic considerations. However, the considerable interpatient variability in response to treatment is a challenge. Markers for a more exact patient classification and stratification are lacking. The objective of this study was to identify markers in immune cell populations that distinguish RA patients from healthy donors with an emphasis on TNF signaling. We employed mass cytometry (CyTOF) with a panel of 13 phenotyping and 10 functional markers to explore signaling in unstimulated and TNF-stimulated peripheral blood mononuclear cells from 20 newly diagnosed, untreated RA patients and 20 healthy donors. The resulting high-dimensional data were analyzed in three independent analysis pipelines, characterized by differences in both data clean-up, identification of cell subsets/clustering and statistical approaches. All three analysis pipelines identified p-p38, IkBa, p-cJun, p-NFkB, and CD86 in cells of both the innate arm (myeloid dendritic cells and classical monocytes) and the adaptive arm (memory CD4+ T cells) of the immune system as markers for differentiation between RA patients and healthy donors. Inclusion of the markers p-Akt and CD120b resulted in the correct classification of 18 of 20 RA patients and 17 of 20 healthy donors in regression modeling based on a combined model of basal and TNF-induced signal. Expression patterns in a set of functional markers and specific immune cell subsets were distinct in RA patients compared to healthy individuals. These signatures may support studies of disease pathogenesis, provide candidate markers for response, and non-response to TNF inhibitor treatment, and aid the identification of future therapeutic targets.

Dataset Information

A towards-multidimensional screening approach to predict candidate genes of rheumatoid arthritis based on SNP, structural and functional annotations.

Background

Methods

Results

Conclusions

Publications

A towards-multidimensional screening approach to predict candidate genes of rheumatoid arthritis based on SNP, structural and functional annotations.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets