Overview of PheWAS with Immune Variants.

<p>This flow chart provides an overview of the steps taken to perform PheWAS between immune variants and ICD-9 diagnosis codes. The final testing dataset (purple) was formed by selecting SNPs from our array data that also exist on Immunochip and/or are within immune-related genes (yellow) and removing samples with missing genotypic or phenotypic data (green). Comprehensive associations were calculated between all final dataset SNPs and ICD-9 code based case/control status using logistic regression, with all models adjusted for age, sex and first five principal components. Replication was sought following both an exact ICD-9 code and a category ICD-9 code approach following the specified criteria. Pooled analysis was performed for both approaches using METAL. See <a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0160573#pone.0160573.s001" target="_blank">S1 Fig</a> for the full workflow from imputation through quality control, association testing, and replication for this study.</p

Summary of Data Sets Used for the Study.

<p>Summary of Data Sets Used for the Study.</p

PheWAS View Plot of Meta-analysis Results with p<0.01 Replicating for the Same ICD-9 Category, Meeting Autoimmune and Immune-Related Diagnosis Criteria.

<p>The left track specifies the phenotype and ICD-9 Category code with which the SNP was associated. The next track indicates–log₁₀(p-value) from the meta-analysis performed on all replicating SNPs with p<0.01. The last track indicates the SNP that had the most significant p-value, and the direction of effect of the association (+, positive; -, negative). The total number of associations between the SNPs and diagnoses was 409.</p

Cytoscape Network Showing the Connections between Phenotypes, the Genes with SNPs, and Pathways.

<p>In this network, green squares represent phenotype; red triangles represent genes; and blue circles are KEGG pathways. The colored lines highlight the link between phenotype and pathway. For the gene <i>HLA-DRA</i> with SNPs associated with “<i>714</i>: <i>rheumatoid arthritis</i>” and “<i>250</i>: <i>type 1 diabetes</i>” is present in the KEGG pathway of “<i>rheumatoid arthritis</i>” (red line) and “<i>type 1 diabetes”</i> (green line) respectively. Also, the blue edge shows the connection between <i>“714</i>: <i>rheumatoid arthritis”</i>, <i>“716</i>: <i>other specified arthropathies”</i> and the KEGG “<i>JAK-STAT signaling pathway</i>” through two interleukin genes, <i>IL23R</i> and <i>IL6</i>.</p

Pleiotropy: SNPs Associated with more than One Phenotype and Replicating across more than One Study for the Same ICD-9 Category.

<p>This chromosomal ideogram has lines indicating the location of the SNP, with filled colored circles indicating different ICD-9 code diagnoses associated with that particular SNP. When there are multiple pairs of the same phenotypes in the same region, this indicates regions where several SNPs in close proximity were associated with the same pairs of phenotypes.</p

Immune-Related ICD-9 categories selected for further analysis.

<p>Immune-Related ICD-9 categories selected for further analysis.</p

Synthesis view plot showing PheWAS results replicating across MyCode^® and BioVU that have previously reported associations.

<p>The first track is the chromosomal location for each SNP. The next column lists the SNP identifier, the phenotype associated in our study, and the reported GWAS trait (p<10⁻⁵). Results representing exact matches with the NHGRI-EBI GWAS catalog and GRASP are annotated with a single asterisk and the closely related traits are represented with a double asterisk. Blue symbols represent results from MyCode^®, red symbols represent results from BioVU and green symbols are the pooled analysis results obtained using the program METAL.</p