Empirical Bayesian Approach to Testing Multiple Hypotheses with Separate Priors for Left and Right Alternatives

We consider a multiple hypotheses problem with directional alternatives in a decision theoretic framework. We obtain an empirical Bayes rule subject to a constraint on mixed directional false discovery rate (mdFDR≤α) under the semiparametric setting where the distribution of the test statistic is parametric, but the prior distribution is nonparametric. We proposed separate priors for the left tail and right tail alternatives as it may be required for many applications. The proposed Bayes rule is compared through simulation against rules proposed by Benjamini and Yekutieli and Efron. We illustrate the proposed methodology for two sets of data from biological experiments: HIV-transfected cell-line mRNA expression data, and a quantitative trait genome-wide SNP data set. We have developed a user-friendly web-based shiny App for the proposed method which is available through URL https://npseb.shinyapps.io/npseb/. The HIV and SNP data can be directly accessed, and the results presented in this paper can be executed

The Decay of Disease Association with Declining Linkage Disequilibrium: A Fine Mapping Theorem

Several important and fundamental aspects of disease genetics models have yet to be described. One such property is the relationship of disease association statistics at a marker site closely linked to a disease causing site. A complete description of this two-locus system is of particular importance to experimental efforts to fine map association signals for complex diseases. Here, we present a simple relationship between disease association statistics and the decline of linkage disequilibrium from a causal site. Specifically, the ratio of Chi-square disease association statistics at a marker site and causal site is equivalent to the standard measure of pairwise linkage disequilibrium, r2. A complete derivation of this relationship from a general disease model is shown. Quite interestingly, this relationship holds across all modes of inheritance. Extensive Monte Carlo simulations using a disease genetics model applied to chromosomes subjected to a standard model of recombination are employed to better understand the variation around this fine mapping theorem due to sampling effects. We also use this relationship to provide a framework for estimating properties of a non-interrogated causal site using data at closely linked markers. Lastly, we apply this way of examining association data from high-density genotyping in a large, publicly-available data set investigating extreme BMI. We anticipate that understanding the patterns of disease association decay with declining linkage disequilibrium from a causal site will enable more powerful fine mapping methods and provide new avenues for identifying causal sites/genes from fine-mapping studies

The 5q31 variants associated with psoriasis and Crohn's disease are distinct

Predisposition to psoriasis is known to be affected by genetic variation in HLA-C, IL12B and IL23R, but other genetic risk factors also exist. We recently reported three psoriasis-associated single nucleotide polymorphisms (SNPs) in the 5q31 locus, a region of high linkage disequilibrium laden with inflammatory pathway genes. The aim of this study was to assess whether other variants in the 5q31 region are causal to these SNPs or make independent contributions to psoriasis risk by genotyping a comprehensive set of tagging SNPs in a 725 kb region bounded by IL3 and IL4 and testing for disease association. Ninety SNPs, capturing 86.4% of the genetic diversity, were tested in one case–control sample set (467 cases/460 controls) and significant markers (Pallelic < 0.05) (n = 9) were then tested in two other sample sets (981 cases/925 controls). All nine SNPs were significant in a meta-analysis of the combined sample sets. Pair-wise conditional association tests showed rs1800925, an intergenic SNP located just upstream of IL13 (Mantel–Haenszel Pcombined = 1.5 × 10−4, OR = 0.77 [0.67–0.88]), could account for observed significant association of all but one other SNP, rs11568506 in SLC22A4 [Mantel–Haenszel Pcombined = 0.043, OR = 0.68 (0.47–0.99)]. Haplotype analysis of these two SNPs showed increased significance for the two common haplotypes (rs11568506–rs1800925: GC, Pcombined = 5.67 × 10−6, OR = 1.37; GT, Pcombined = 6.01 × 10−5, OR = 0.75; global haplotype P = 8.93 × 10−5). Several 5q31-region SNPs strongly associated with Crohn's disease (CD) in the recent WTCCC study were not significant in the psoriasis sample sets tested here. These results identify the most significant 5q31 risk variants for psoriasis and suggest that distinct 5q31 variants contribute to CD and psoriasis risk

A Candidate Gene Approach Identifies the TRAF1/C5 Region as a Risk Factor for Rheumatoid Arthritis

Using a candidate-gene approach, Rene Toes and colleagues identified a novel genetic risk factor for rheumatoid arthritis in theTRAF1/C5 region

A comprehensive review of Tripterygium wilfordii hook. f. in the treatment of rheumatic and autoimmune diseases: Bioactive compounds, mechanisms of action, and future directions

Rheumatic and autoimmune diseases are a group of immune system-related disorders wherein the immune system mistakenly attacks and damages the body’s tissues and organs. This excessive immune response leads to inflammation, tissue damage, and functional impairment. Therapeutic approaches typically involve medications that regulate immune responses, reduce inflammation, alleviate symptoms, and target specific damaged organs. Tripterygium wilfordii Hook. f., a traditional Chinese medicinal plant, has been widely studied in recent years for its application in the treatment of autoimmune diseases, including rheumatoid arthritis, systemic lupus erythematosus, and multiple sclerosis. Numerous studies have shown that preparations of Tripterygium wilfordii have anti-inflammatory, immunomodulatory, and immunosuppressive effects, which effectively improve the symptoms and quality of life of patients with autoimmune diseases, whereas the active metabolites of T. wilfordii have been demonstrated to inhibit immune cell activation, regulate the production of inflammatory factors, and modulate the immune system. However, although these effects contribute to reductions in inflammatory responses and the suppression of autoimmune reactions, as well as minimize tissue and organ damage, the underlying mechanisms of action require further investigation. Moreover, despite the efficacy of T. wilfordii in the treatment of autoimmune diseases, its toxicity and side effects, including its potential hepatotoxicity and nephrotoxicity, warrant a thorough assessment. Furthermore, to maximize the therapeutic benefits of this plant in the treatment of autoimmune diseases and enable more patients to utilize these benefits, efforts should be made to strengthen the regulation and standardized use of T. wilfordii

Neither Replication nor Simulation Supports a Role for the Axon Guidance Pathway in the Genetics of Parkinson's Disease

Susceptibility to sporadic Parkinson's disease (PD) is thought to be influenced by both genetic and environmental factors and their interaction with each other. Statistical models including multiple variants in axon guidance pathway genes have recently been purported to be capable of predicting PD risk, survival free of the disease and age at disease onset; however the specific models have not undergone independent validation. Here we tested the best proposed risk panel of 23 single nucleotide polymorphisms (SNPs) in two PD sample sets, with a total of 525 cases and 518 controls. By single marker analysis, only one marker was significantly associated with PD risk in one of our sample sets (rs6692804: P = 0.03). Multi-marker analysis using the reported model found a mild association in one sample set (two sided P = 0.049, odds ratio for each score change = 1.07) but no significance in the other (two sided P = 0.98, odds ratio = 1), a stark contrast to the reported strong association with PD risk (P = 4.64×10−38, odds ratio as high as 90.8). Following a procedure similar to that used to build the reported model, simulated multi-marker models containing SNPs from randomly chosen genes in a genome wide PD dataset produced P-values that were highly significant and indistinguishable from similar models where disease status was permuted (3.13×10−23 to 4.90×10−64), demonstrating the potential for overfitting in the model building process. Together, these results challenge the robustness of the reported panel of genetic markers to predict PD risk in particular and a role of the axon guidance pathway in PD genetics in general

A Large-Scale Rheumatoid Arthritis Genetic Study Identifies Association at Chromosome 9q33.2

Rheumatoid arthritis (RA) is a chronic, systemic autoimmune disease affecting both joints and extra-articular tissues. Although some genetic risk factors for RA are well-established, most notably HLA-DRB1 and PTPN22, these markers do not fully account for the observed heritability. To identify additional susceptibility loci, we carried out a multi-tiered, case-control association study, genotyping 25,966 putative functional SNPs in 475 white North American RA patients and 475 matched controls. Significant markers were genotyped in two additional, independent, white case-control sample sets (661 cases/1322 controls from North America and 596 cases/705 controls from The Netherlands) identifying a SNP, rs1953126, on chromosome 9q33.2 that was significantly associated with RA (ORcommon = 1.28, trend Pcomb = 1.45E-06). Through a comprehensive fine-scale-mapping SNP-selection procedure, 137 additional SNPs in a 668 kb region from MEGF9 to STOM on 9q33.2 were chosen for follow-up genotyping in a staged-approach. Significant single marker results (Pcomb<0.01) spanned a large 525 kb region from FBXW2 to GSN. However, a variety of analyses identified SNPs in a 70 kb region extending from the third intron of PHF19 across TRAF1 into the TRAF1-C5 intergenic region, but excluding the C5 coding region, as the most interesting (trend Pcomb: 1.45E-06 → 5.41E-09). The observed association patterns for these SNPs had heightened statistical significance and a higher degree of consistency across sample sets. In addition, the allele frequencies for these SNPs displayed reduced variability between control groups when compared to other SNPs. Lastly, in combination with the other two known genetic risk factors, HLA-DRB1 and PTPN22, the variants reported here generate more than a 45-fold RA-risk differential

The Impact of Diagnostic Code Misclassification on Optimizing the Experimental Design of Genetic Association Studies

Diagnostic codes within electronic health record systems can vary widely in accuracy. It has been noted that the number of instances of a particular diagnostic code monotonically increases with the accuracy of disease phenotype classification. As a growing number of health system databases become linked with genomic data, it is critically important to understand the effect of this misclassification on the power of genetic association studies. Here, I investigate the impact of this diagnostic code misclassification on the power of genetic association studies with the aim to better inform experimental designs using health informatics data. The trade-off between (i) reduced misclassification rates from utilizing additional instances of a diagnostic code per individual and (ii) the resulting smaller sample size is explored, and general rules are presented to improve experimental designs