Candidate Causal Single Nucleotide Polymorphisms Research Articles

THU0019 Genome-wide pathway analysis of genome-wide association studies on psoriasis and behcet’s disease

BackgroundGenome-wide association studies (GWASs) have been successfully used to identify novel common genetic variants that contribute to susceptibility to complex diseases, but individual GWASs are limited in terms of identifying...

Pathway analysis of a genome-wide association study in schizophrenia

ObjectiveThe aim of this study was to identify the candidate single nucleotide polymorphisms (SNPs) and candidate mechanisms that contribute to schizophrenia susceptibility and to generate a SNP to gene to pathway hypothesis using an analytical pathway-based approach. MethodsWe used schizophrenia GWAS data of the genotypes of 660,259 SNPs in 1378 controls and 1351 cases of European descent after quality control filtering. ICSNPathway (Identify candidate Causal SNPs and Pathways) analysis was applied to the schizophrenia GWAS dataset. The first stage involved the pre-selection of candidate SNPs by linkage disequilibrium analysis and the functional SNP annotation of the most significant SNPs found. The second stage involved the annotation of biological mechanisms for the pre-selected candidate SNPs using improved-gene set enrichment analysis. ResultsICSNPathway analysis identified fifteen candidate SNPs, ten candidate pathways, and nine hypothetical biological mechanisms. The most strongly associated potential pathways were as follows. First, rs1644731 and rs1644730 to RDH8 to estrogen biosynthetic process (p<0.001, FDR<0.001). The genes involved in this pathway are RDH8 and HSD3B1 (p<0.05). All-trans-retinol dehydrogenase (RDH8) is a visual cycle enzyme that reduces all-trans-retinal to all-trans-retinol in the presence of NADPH. The chemical reactions and pathways involved result in the formation of estrogens, which are C18 steroid hormones that can stimulate the development of female sexual characteristics. Second, rs1146031 to ACVR1 to mesoderm formation and activin binding (p<0.001, FDR=0.032, 0.034). Two of 15 candidate genes are known genes associated with schizophrenia: KCNQ2 and APOL2. One of the 10 candidate pathways, estrogen biosynthetic process, is known to be associated with schizophrenia (p<0.001, FDR<0.001). However, 13 of candidate genes (RDH8, ACVR1, PSMD9, KCNAB1, SLC17A3, ARCN1, COG7, STAB2, LRPAP1, STAB1, CXCL16, COL4A4, EXOSC3) and 9 of candidate pathways were novel. ConclusionBy applying ICSNPathway analysis to schizophrenia GWAS data, we identified candidate SNPs, genes like KCNQ2 and APOL2 and pathways involving the estrogen biosynthetic process may contribute to schizophrenia susceptibility. Further analyses are needed to validate the results of this analysis.

Genome-wide pathway analysis of a genome-wide association study on multiple sclerosis

The aims of this study were to identify candidate single nucleotide polymorphisms (SNPs) and mechanisms of multiple sclerosis (MS) and to generate SNP to gene to pathway hypotheses. A MS genome-wide association study (GWAS) dataset that included 505,763 SNPs in 500 cases and 500 controls of European descent was used in this study. Identify candidate Causal SNPs and Pathway (ICSNPathway) analysis was applied to the GWAS dataset. ICSNPathway analysis identified 9 candidate SNPs and 5 pathways, which provided 5 hypothetical biological mechanisms. The candidate SNPs, namely, rs1802127 (MSH5), rs9277471 (human leukocyte antigen [HLA]-DPB1), rs8084 (HLA-DRA), rs7192 (HLA-DRA), rs2072895 (HLA-F), rs2735059 (HLA-F), rs915669 (HLA-G), rs915668 (HLA-G), and rs1063320 (HLA-G) were all at HLA loci (-log10(P) = 3.301-4.000). The most strongly associated pathway was rs1802127 to MSH5 to meiotic recombination and meiotic cell cycle (nominal P < 0.001, false discovery rate [FDR] < 0.001). When HLA loci were excluded, ICSNPathway analysis identified seven candidate non-HLA SNPs (rs5896 [F2], rs8181979 [SHC1], rs9297605 [TAF2], rs669 [A2 M], rs2228043 [IL6ST], rs1061622 [TNFRSF1B], rs1801516 [ATM]) and ten candidate causal pathways, which provided seven hypothetical biological mechanisms (nominal P ≤ 0.001, FDR ≤ 0.047). The most strongly associated pathway was SNP rs5896 to F2 to the transcriptional activation DNA-binding protein B from mRNA (nominal P < 0.001, FDR = 0.006). The application of ICSNPathway analysis to the MS GWAS dataset resulted in the identification of candidate SNPs, pathways, and biological mechanisms that might contribute to MS susceptibility.

Pathway analysis of genome-wide association studies for Parkinson’s disease

The study was done to identify the candidate causal single nucleotide polymorphisms (SNPs) and candidate causal mechanisms that contribute to Parkinson's disease (PD) susceptibility and to generate a SNP to ene to pathway hypothesis using an analytical pathway-based approach. We used a PD genome-wide association study (GWAS) meta-analysis data of the genotypes of 2,525,705 SNPs in 4,238 PD cases and 4,239 controls. Identify candidate Causal SNPs and Pathways (ICSNPathway) analysis was applied to the PD GWAS dataset. The first stage involved the pre-selection of candidate causal SNPs by linkage disequilibrium analysis and the functional SNP annotation of the most significant SNPs found. The second stage involved the annotation of biological mechanisms for the pre-selected candidate causal SNPs using improved-gene set enrichment analysis. ICSNPathway analysis identified three candidate SNPs, two genes, twenty-one pathways, and three hypothetical biological mechanisms: (1) rs17651549 to microtubule-associated protein tau (MAPT) to protein domain specific binding (nominal p < 0.001, false discovery rate (FDR) < 0.001), neurogenesis (nominal p < 0.001, FDR < 0.001), regulation of neurogenesis (nominal p < 0.001, FDR = 0.001), positive regulation of axonogenesis (nominal p < 0.001, FDR = 0.001), regulation of protein polymerization (nominal p < 0.001, FDR = 0.004), negative regulation of organelle organization (nominal p < 0.001, FDR = 0.004), hsa01510 (nominal p < 0.001, FDR = 0.005), neuron differentiation (nominal p < 0.001, FDR = 0.009), and axonogenesis (nominal p < 0.001, FDR = 0.009); (2) rs10445337 to MAPT to protein domain specific binding (nominal p < 0.001, FDR < 0.001), neurogenesis (nominal p < 0.001, FDR < 0.001), regulation of neurogenesis (nominal p < 0.001, FDR = 0.001), and positive regulation of axonogenesis (nominal p < 0.001, FDR = 0.001); (3) rs9938550 to HSD3B7 to hsa00363 (nominal p < 0.001, FDR = 0.004), bile acid metabolic process (nominal p = 0.005, FDR = 0.019), and steroid metabolic process (nominal p = 0.010, FDR = 0.039). By applying the ICSNPathway analysis to PD GWAS meta-analysis data, three candidate SNPs, two genes (MAPT and HSD3B7), and 21 pathways involving protein domain specific binding and neurogenesis were identified, which may contribute to PD susceptibility.

Pathway analysis of genome-wide association study on asthma

The aims of this study were to identify the candidate causal single nucleotide polymorphisms (SNPs) and candidate causal mechanisms of asthma and to generate SNP to gene to pathway hypotheses. SNPs that met a threshold of p ≤ 0.001 in a genome-wide association study (GWAS) dataset of asthma, which included 292,443 SNPs in 473 asthma cases and 1892 controls, were used in the present study. Identify candidate causal SNPs and pathway (ICSNPathway) analysis was applied to this dataset. ICSNPathway analysis identified four candidate causal SNPs, four genes, and 21 candidate causal pathways, which in total provided four hypothetical biologic mechanisms: (1) rs7192 (nonsynonymous coding) to HLA-DRA to 21 pathways, such as, the role of eosinophils in the chemokine network of allergy, Th1/Th2 differentiation, and asthma (nominal p ≤ 0.001, FDR p ≤ 0.01); (2) rs20541 (nonsynonymous coding) to IL13 to asthma and cytokines and inflammatory response (nominal p<0.001, FDR p ≤ 0.008); (3) rs1058808 (frameshift coding) to ERBB2 to transmembrane receptor activity (nominal p=0.001, FDR p=0.01); (4) rs17350764 (nonsynonymous coding (deleterious)) to OR52J3 to transmembrane receptor activity (nominal p=0.001, FDR p=0.01). By applying ICSNPathway analysis to asthma GWAS data, we found four candidate causal SNPs, four genes involving HLA-DRA and IL-13, and four hypotheses, which may contribute to asthma susceptibility.

Genome-wide pathway analysis of genome-wide association studies on systemic lupus erythematosus and rheumatoid arthritis

The aim of this study was to explore candidate single nucleotide polymorphisms (SNPs) and candidate mechanisms of systemic lupus erythematosus (SLE) and rheumatoid arthritis (RA). Two SLE genome-wide association studies (GWASs) datasets were included in this study. Meta-analysis was conducted using 737,984 SNPs in 1,527 SLE cases and 3,421 controls of European ancestry, and 4,429 SNPs that met a threshold of p < 0.01 in a Korean RA GWAS dataset was used. ICSNPathway (identify candidate causal SNPs and pathways) analysis was applied to the meta-analysis results of the SLE GWAS datasets, and a RA GWAS dataset. The most significant result of SLE GWAS meta-analysis concerned rs2051549 in the human leukocyte antigen (HLA) region (p = 3.36E-22). In the non-HLA region, meta-analysis identified 6 SNPs associated with SLE with genome-wide significance (STAT4, TNPO3, BLK, FAM167A, and IRF5). ICSNPathway identified five candidate causal SNPs and 13 candidate causal pathways. This pathway-based analysis provides three hypotheses of the biological mechanism involved. First, rs8084 and rs7192 → HLA-DRA → bystander B cell activation. Second, rs1800629 → TNF → cytokine network. Third, rs1150752 and rs185819 → TNXB → collagen metabolic process. ICSNPathway analysis identified three candidate causal non-HLA SNPs and four candidate causal pathways involving the PADI4, MTR, PADI2, and TPH2 genes of RA. We identified five candidate SNPs and thirteen pathways, involving bystander B cell activation, cytokine network, and collagen metabolic processing, which may contribute to SLE susceptibility, and we revealed candidate causal non-HLA SNPs, genes, and pathways of RA.

Pathway Analysis of a Genome-Wide Association Study of Ileal Crohn's Disease

The aim of this study was to identify the candidate single-nucleotide polymorphisms (SNPs) and candidate mechanisms of Crohn's disease (CD) and to generate a SNP to gene to pathway hypothesis using a pathway-based approach, ICSNPathway (Identify candidate Causal SNPs and Pathways) of genome-wide association studies (GWAS). An ileal CD GWAS dataset downloaded from NCBI dbGap was used in this study, which was conducted using 297,857 SNPs in 968 CD cases and 995 controls after quality control filtering. ICSNPathway analysis was applied to the CD GWAS dataset. ICSNPathway analysis identified seven candidate SNPs, nine pathways, which provided five hypothetical biological mechanisms. First, rs4077515 to caspase recruitment domain-containing protein 9 (CARD9) to response to peptidoglycan, positive regulation of tumor necrosis factor production, response to exogenous dsRNA, positive regulation of interleukin-6 production, and positive regulation of innate immune response (nominal p≤0.002, false discovery rate [FDR]=0.027). Second, rs2066842, rs3135500, and rs5743291 to NOD2 to response to pepditoglycan, innate immune response activating signal transduction, and positive regulation of innate immune response (nominal p≤0.001, FDR=0.027). Third, rs8172678 to PPARGC1A to cellular glucose homeostasis (nominal p<0.001, FDR=0.031). Fourth, rs1050152 to SLC22A4 to cofactor transporter activity and vitamin transport (nominal p<0.001, FDR=0.044). Fifth, rs9621049 to TCN2 to vitamin transport (nominal p<0.001, FDR=0.046). In conclusion, by applying ICSNPathway analysis to CD GWAS, we identified candidate SNPs, genes involving CARD9, NOD2, PPARGC1A, SLC22A4, and TCN2, pathways, and hypothetical mechanisms, which may contribute to ileal CD susceptibility.

Pathway analysis of genome-wide association study for bone mineral density

The aim of this study was to identify the candidate causal single nucleotide polymorphisms (SNPs) and candidate causal mechanisms that contribute to bone mineral density (BMD) and to generate a SNP to gene to pathway hypothesis using an analytical pathway-based approach. We used hip BMD GWAS data of the genotypes of 301,019 SNPs in 5,715 Europeans. ICSNPathway (identify candidate causal SNPs and pathways) analysis was applied to the BMD GWAS dataset. The first stage involved the pre-selection of candidate causal SNPs by linkage disequilibrium analysis and the functional SNP annotation of the most significant SNPs found. The second stage involved the annotation of biological mechanisms for the pre-selected candidate causal SNPs using improved-gene set enrichment analysis. ICSNPathway analysis identified seven candidate SNPs, eight candidate pathways, and seven hypothetical biological mechanisms. Eight pathways are as follows; gamma-hexachlorocyclohexane degradation (nominal p-value < 0.001, false discovery rate (FDR) <0.001), regulation of the smoothened signaling pathway (nominal p-value < 0.001, FDR = 0.016), TACI and BCMA stimulation of B cell immune response (nominal p-value < 0.001, FDR = 0.021), endonuclease activity (nominal p-value = 0.001, FDR = 0,026), regulation of defense response to virus (nominal p-value = 0.001, FDR = 0.028), serine_type_endopeptidase_inhibitor_activity (nominal p-value = 0.001, FDR = 0.044), endoribonuclease activity (nominal p-value = 0.002, FDR = 0.045), and myeloid leukocyte differentiation (nominal p-value = 0.001, FDR = 0.050). The most significant causal pathway was gamma-hexachlorocyclohexane degradation. CYP3A5, PON2, PON3, CMBL, PON1, ALPL, CYP3A43, CYP3A7, ACP6, ACPP, and ALPI (p < 0.05) are involved in the pathway of gamma-hexachlorocyclohexane degradation. Further examination of the gene contents revealed that DBR1, DICER1, EXO1, FEN1, POP1, POP4, RPP30, and RPP38 were involved in 2 of the 8 pathways (p < 0.05). By applying ICSNPathway analysis to BMD GWAS data, we identified seven candidate SNPs and eight pathways involving gamma-hexachlorocyclohexane degradation, which may contribute to low BMD.

Genome-wide pathway analysis of a genome-wide association study on psoriasis and Behcet’s disease

The aim of this study was to identify candidate causal single nucleotide polymorphisms (SNPs) and candidate causal mechanisms of psoriasis and Behcets's disease (BD) and to generate an SNP→gene→pathway hypothesis. A psoriasis genome-wide association study (GWAS) dataset that included 436,192 SNPs in 1,409 psoriasis cases and 1,436 controls of European descent and a BD GWAS dataset that contained 310,324 SNPs in 1,215 BD cases and 1,278 controls were used in this study. Identify candidate causal SNPs and pathways (ICSNPathway) analysis was applied to the GWAS datasets. ICSNPathway analysis identified 15 candidate causal SNPs and 28 candidate causal pathways. The top five candidate causal SNPs were rs1063478 (P=1.45E-10), rs8084 (P=2.20E-08), rs7192 (P=5.18E-08), rs20541 (P=5.30E-06), and rs1130838 (P=5.65E-06), which with the exception of rs20541 [interleukin (IL)-13] are at human leukocyte antigen (HLA) loci. These candidate causal SNPs and pathways provided ten hypothetical biological mechanisms. The most strongly associated pathway concerned HLA. When HLA loci were excluded, ICSNPathway analysis provided one hypothetical biological mechanism. rs20541 (non_synonymous_coding)→IL-13→dendritic cell involvement in the regulation of Th1 and Th2 development, and the GATA3 pathway. ICSNPathway analysis identified four candidate causal SNPs, eleven candidate causal pathways, and three hypothetical biological mechanisms. One of them was as follows: rs2072895 (non_synonymous_coding & splice-site) and rs2735059 (non_synonymous_coding)→HLA-F→type I diabetes mellitus, antigen processing and presentation, and autoimmune thyroid disease. The application of ICSNPathway analysis to GWAS dataset of psoriasis and BD resulted in the identification of candidate causal SNPs and candidate pathways that might contribute to psoriasis susceptibility.

FGFR2 variants and breast cancer risk: fine-scale mapping using African American studies and analysis of chromatin conformation

Genome-wide association studies have identified FGFR2 as a breast cancer (BC) susceptibility gene in populations of European and Asian descent, but a causative variant has not yet been conclusively identified. We hypothesized that the weaker linkage disequilibrium across this associated region in populations of African ancestry might help refine the set of candidate-causal single nucleotide polymorphisms (SNPs) previously identified by our group. Eight candidate-causal SNPs were evaluated in 1253 African American invasive BC cases and 1245 controls. A significant association with BC risk was found with SNP rs2981578 (unadjusted per-allele odds ratio = 1.20, 95% confidence interval 1.03-1.41, P(trend) = 0.02), with the odds ratio estimate similar to that reported in European and Asian subjects. To extend the fine-mapping, genotype data from the African American studies were analyzed jointly with data from European (n = 7196 cases, 7275 controls) and Asian (n = 3901 cases, 3205 controls) studies. In the combined analysis, SNP rs2981578 was the most strongly associated. Five other SNPs were too strongly correlated to be excluded at a likelihood ratio of < 1/100 relative to rs2981578. Analysis of DNase I hypersensitive sites indicated that only two of these map to highly accessible chromatin, one of which, SNP rs2981578, has previously been implicated in up-regulating FGFR2 expression. Our results demonstrate that the association of SNPs in FGFR2 with BC risk extends to women of African American ethnicity, and illustrate the utility of combining association analysis in datasets of diverse ethnic groups with functional experiments to identify disease susceptibility variants.