Confident inference for SNP effects on treatment efficacy. (English) Zbl 1405.62170

Ding, Ying; Li, Ying Grace; Liu, Yushi; Ruberg, Stephen J.; Hsu, Jason C.

Ann. Appl. Stat. 12, No. 3, 1727-1748 (2018).

Summary: Our research is for finding SNPs that are predictive of treatment efficacy, to decide which subgroup (with enhanced treatment efficacy) to target in drug development. Testing SNPs for lack of association with treatment outcome is inherently challenging, because any linkage disequilibrium between a noncausal SNP with a causal SNP, however small, makes the zero-null (no association) hypothesis technically false. Control of Type I error rate in testing such null hypotheses are therefore difficult to interpret. We propose a completely different formulation to address this problem. For each SNP, we provide simultaneous confidence intervals directed toward detecting possible dominant, recessive, or additive effects. Across the SNPs, we control the expected number of SNPs with at least one false confidence interval coverage. Since our confidence intervals are constructed based on pivotal statistics, the false coverage control is guaranteed to be exact and unaffected by the true values of test quantities (whether zero or nonzero). Our method is applicable to the therapeutic areas of Diabetes and Alzheimer’s diseases, and perhaps more, as a step toward confidently targeting a patient subgroup in a tailored drug development process.

MSC:

62P10	Applications of statistics to biology and medical sciences; meta analysis
62J15	Paired and multiple comparisons; multiple testing
62H20	Measures of association (correlation, canonical correlation, etc.)
62G15	Nonparametric tolerance and confidence regions

Keywords:

multiple testing; simultaneous confidence intervals; SNP; tailored drug development; treatment efficacy

References:

[1]	Berger, R. L. and Hsu, J. C. (1996). Bioequivalence trials, intersection-union tests, and equivalence confidence sets. Statist. Sci.11 283–315. · Zbl 0955.62555 · doi:10.1214/ss/1032280304
[2]	de Bakker, P., McVean, G., Sabeti, P., Miretti, M., Green, T., Marchini, J., Ke, X., Monsuur, A., Whittaker, P., Delgado, M., Morrison, J., Richardson, A., Walsh, E., Gao, X., Galver, L., Hart, J., Hafler, D., Pericak-Vance, M., Todd, J., Daly, M., Trowsdale, J., Wijmenga, C., Vyse, T., Beck, S., Murray, S., Carrington, M., Gregory, S., Deloukas, P. and Rioux, J. (2006). A high-resolution HLA and SNP haplotype map for disease association studies in the extended human MHC. Nat. Genet.38 1166–1172.
[3]	Efron, B. (2007). Correlation and large-scale simultaneous significance testing. J. Amer. Statist. Assoc.102 93–103. · Zbl 1284.62340 · doi:10.1198/016214506000001211
[4]	Genz, A. and Bretz, F. (1999). Numerical computation of multivariate \(t\)-probabilities with application to power calculation of multiple contrasts. J. Stat. Comput. Simul.63 361–378. · Zbl 0934.62020 · doi:10.1080/00949659908811962
[5]	Hothorn, L. and Hothorn, T. (2009). Order-restricted scores test for the evaluation of population-based case-control studies when the genetic model is unknown. Biom. J.51 659–669. · Zbl 1442.62414
[6]	Hsu, J. C. (1996). Multiple Comparisons: Theory and Methods. Chapman & Hall, London. · Zbl 0898.62090
[7]	Lettre, G., Lange, C. and Hirschhorn, J. (2007). Genetic model testing and statistical power in population-based association studies of quantitative traits. Genet. Epidemiol.31 358–362.
[8]	Lipkovich, I., Dmeitrienko, A. and D’Agostino, R. B. (2017). Tutorial in biostatistics: Data-driven subgroup identification and analysis in clinical trials. Stat. Med.36 136–196.
[9]	Loh, W.-Y., He, X. and Man, M. (2015). A regression tree approach to identifying subgroups with differential treatment effects. Stat. Med.34 1818–1833. · doi:10.1002/sim.6454
[10]	Mallal, S., Nolan, D., Witt, C., Masel, G., Martin, A. M., Moore, C., Sayer, D., Castley, A., Mamotte, C., Maxwell, D., James, I. and Christiansen, F. T. (2002). Association between presence of HLA-B\(^{*}\)5701, HLA-DR7, and HLA-DQ3 and hypersensitivity to HIV-1 reverse-transcriptase inhibitor abacavir. Lancet359 727–732.
[11]	Mallal, S., Phillips, E., Carosi, G., Molina, J.-M., Workman, C., Tomažič, J., Jägel-Guedes, E., Rugina, S., Kozyrev, O., Cid, J. F., Hay, P., Nolan, D., Hughes, S., Hughes, A., Ryan, S., Fitch, N., Thorborn, D. and Benbow, A. (2008). HLA-B\(^{*}\)5701 screening for hypersensitivity to abacavir. N. Engl. J. Med.358 568–579.
[12]	So, H.-C. and Sham, P. C. (2011). Robust association tests under different genetic models, allowing for binary or quantitative traits and covariates. Behav. Genet.41 768–775.
[13]	The 1000 Genomes Project Consortium (2010). A map of human genome variation from population-scale sequencing. Nature467 1061–1073.
[14]	The 1000 Genomes Project Consortium (2012). An integrated map of genetic variation from 1092 human genomes. Nature491 56–65.
[15]	The 1000 Genomes Project Consortium (2015). A global reference for human genetic variation. Nature526 68–74.
[16]	Tukey, J. W. (1992). Where should multiple comparisons go next? In Multiple Comparisons, Selection, and Applications in Biometry: A Festschrift in Honor of Charles W. Dunnett (F. M. Hoppe, ed.) Chapter 12 187–208.