Kernel machine SNP-set analysis for censored survival outcomes in genome-wide association studies

BayesKAT: bayesian optimal kernel-based test for genetic association studies reveals joint genetic effects in complex diseases

Genome-wide Association Studies (GWAS) methods have identified individual single-nucleotide polymorphisms (SNPs) significantly associated with specific phenotypes. Nonetheless, many complex diseases are polygenic and are controlled by multiple genetic variants that are usually non-linearly dependent. These genetic variants are marginally less effective and remain undetected in GWAS analysis. Kernel-based tests (KBT), which evaluate the joint effect of a group of genetic variants, are therefore critical for complex disease analysis. However, choosing different kernel functions in KBT can significantly influence the type I error control and power, and selecting the optimal kernel remains a statistically challenging task. A few existing methods suffer from inflated type 1 errors, limited scalability, inferior power or issues of ambiguous conclusions. Here, we present a new Bayesian framework, BayesKAT (https://github.com/wangjr03/BayesKAT), which overcomes these kernel specification issues by selecting the optimal composite kernel adaptively from the data while testing genetic associations simultaneously. Furthermore, BayesKAT implements a scalable computational strategy to boost its applicability, especially for high-dimensional cases where other methods become less effective. Based on a series of performance comparisons using both simulated and real large-scale genetics data, BayesKAT outperforms the available methods in detecting complex group-level associations and controlling type I errors simultaneously. Applied on a variety of groups of functionally related genetic variants based on biological pathways, co-expression gene modules and protein complexes, BayesKAT deciphers the complex genetic basis and provides mechanistic insights into human diseases.

BayesKAT: Bayesian Optimal Kernel-based Test for genetic association studies reveals joint genetic effects in complex diseases

GWAS methods have identified individual SNPs significantly associated with specific phenotypes. Nonetheless, many complex diseases are polygenic and are controlled by multiple genetic variants that are usually non-linearly dependent. These genetic variants are marginally less effective and remain undetected in GWAS analysis. Kernel-based tests (KBT), which evaluate the joint effect of a group of genetic variants, are therefore critical for complex disease analysis. However, choosing different kernel functions in KBT can significantly influence the type I error control and power, and selecting the optimal kernel remains a statistically challenging task. A few existing methods suffer from inflated type 1 errors, limited scalability, inferior power, or issues of ambiguous conclusions. Here, we present a new Bayesian framework, BayesKAT( https://github.com/wangjr03/BayesKAT ), which overcomes these kernel specification issues by selecting the optimal composite kernel adaptively from the data while testing genetic associations simultaneously. Furthermore, BayesKAT implements a scalable computational strategy to boost its applicability, especially for high-dimensional cases where other methods become less effective. Based on a series of performance comparisons using both simulated and real large-scale genetics data, BayesKAT outperforms the available methods in detecting complex group-level associations and controlling type I errors simultaneously. Applied on a variety of groups of functionally related genetic variants based on biological pathways, co-expression gene modules, and protein complexes, BayesKAT deciphers the complex genetic basis and provides mechanistic insights into human diseases.

Bayesian bi-level variable selection for genome-wide survival study

Mild cognitive impairment (MCI) is a clinical syndrome characterized by the onset and evolution of cognitive impairments, often considered a transitional stage to Alzheimer's disease (AD). The genetic traits of MCI patients who experience a rapid progression to AD can enhance early diagnosis capabilities and facilitate drug discovery for AD. While a genome-wide association study (GWAS) is a standard tool for identifying single nucleotide polymorphisms (SNPs) related to a disease, it fails to detect SNPs with small effect sizes due to stringent control for multiple testing. Additionally, the method does not consider the group structures of SNPs, such as genes or linkage disequilibrium blocks, which can provide valuable insights into the genetic architecture. To address the limitations, we propose a Bayesian bi-level variable selection method that detects SNPs associated with time of conversion from MCI to AD. Our approach integrates group inclusion indicators into an accelerated failure time model to identify important SNP groups. Additionally, we employ data augmentation techniques to impute censored time values using a predictive posterior. We adapt Dirichlet-Laplace shrinkage priors to incorporate the group structure for SNP-level variable selection. In the simulation study, our method outperformed other competing methods regarding variable selection. The analysis of Alzheimer's Disease Neuroimaging Initiative (ADNI) data revealed several genes directly or indirectly related to AD, whereas a classical GWAS did not identify any significant SNPs.

How can childhood maltreatment affect post-traumatic stress disorder in adult: Results from a composite null hypothesis perspective of mediation analysis

Background

A greatly growing body of literature has revealed the mediating role of DNA methylation in the influence path from childhood maltreatment to psychiatric disorders such as post-traumatic stress disorder (PTSD) in adult. However, the statistical method is challenging and powerful mediation analyses regarding this issue are lacking.

Methods

To study how the maltreatment in childhood alters long-lasting DNA methylation changes which further affect PTSD in adult, we here carried out a gene-based mediation analysis from a perspective of composite null hypothesis in the Grady Trauma Project (352 participants and 16,565 genes) with childhood maltreatment as exposure, multiple DNA methylation sites as mediators, and PTSD or its relevant scores as outcome. We effectively addressed the challenging issue of gene-based mediation analysis by taking its composite null hypothesis testing nature into consideration and fitting a weighted test statistic.

Results

We discovered that childhood maltreatment could substantially affected PTSD or PTSD-related scores, and that childhood maltreatment was associated with DNA methylation which further had significant roles in PTSD and these scores. Furthermore, using the proposed mediation method, we identified multiple genes within which DNA methylation sites exhibited mediating roles in the influence path from childhood maltreatment to PTSD-relevant scores in adult, with 13 for Beck Depression Inventory and 6 for modified PTSD Symptom Scale, respectively.

Conclusion

Our results have the potential to confer meaningful insights into the biological mechanism for the impact of early adverse experience on adult diseases; and our proposed mediation methods can be applied to other similar analysis settings.

Intravesical BCG in patients with non-muscle invasive bladder cancer induces trained immunity and decreases respiratory infections

BACKGROUND

BCG is recommended as intravesical immunotherapy to reduce the risk of tumor recurrence in patients with non-muscle invasive bladder cancer (NMIBC). Currently, it is unknown whether intravesical BCG application induces trained immunity.

METHODS

The aim of this research was to determine whether BCG immunotherapy induces trained immunity in NMIBC patients. We conducted a prospective observational cohort study in 17 NMIBC patients scheduled for BCG therapy and measured trained immunity parameters at 9 time points before and during a 1-year BCG maintenance regimen. Ex vivo cytokine production by peripheral blood mononuclear cells, epigenetic modifications, and changes in the monocyte transcriptome were measured. The frequency of respiratory infections was investigated in two larger cohorts of BCG-treated and non-BCG treated NMIBC patients as a surrogate measurement of trained immunity. Gene-based association analysis of genetic variants in candidate trained immunity genes and their association with recurrence-free survival and progression-free survival after BCG therapy was performed to investigate the hypothesized link between trained immunity and clinical response.

RESULTS

We found that intravesical BCG does induce trained immunity based on an increased production of TNF and IL-1β after heterologous ex vivo stimulation of circulating monocytes 6-12 weeks after intravesical BCG treatment; and a 37% decreased risk (OR 0.63 (95% CI 0.40 to 1.01)) for respiratory infections in BCG-treated versus non-BCG-treated NMIBC patients. An epigenomics approach combining chromatin immuno precipitation-sequencing and RNA-sequencing with in vitro trained immunity experiments identified enhanced inflammasome activity in BCG-treated individuals. Finally, germline variation in genes that affect trained immunity was associated with recurrence and progression after BCG therapy in NMIBC.

CONCLUSION

We conclude that BCG immunotherapy induces trained immunity in NMIBC patients and this may account for the protective effects against respiratory infections. The data of our gene-based association analysis suggest that a link between trained immunity and oncological outcome may exist. Future studies should further investigate how trained immunity affects the antitumor immune responses in BCG-treated NMIBC patients.

A comprehensive comparison of multilocus association methods with summary statistics in genome-wide association studies

BACKGROUND

Multilocus analysis on a set of single nucleotide polymorphisms (SNPs) pre-assigned within a gene constitutes a valuable complement to single-marker analysis by aggregating data on complex traits in a biologically meaningful way. However, despite the existence of a wide variety of SNP-set methods, few comprehensive comparison studies have been previously performed to evaluate the effectiveness of these methods.

RESULTS

We herein sought to fill this knowledge gap by conducting a comprehensive empirical comparison for 22 commonly-used summary-statistics based SNP-set methods. We showed that only seven methods could effectively control the type I error, and that these well-calibrated approaches had varying power performance under the simulation scenarios. Overall, we confirmed that the burden test was generally underpowered and score-based variance component tests (e.g., sequence kernel association test) were much powerful under the polygenic genetic architecture in both common and rare variant association analyses. We further revealed that two linkage-disequilibrium-free P value combination methods (e.g., harmonic mean P value method and aggregated Cauchy association test) behaved very well under the sparse genetic architecture in simulations and real-data applications to common and rare variant association analyses as well as in expression quantitative trait loci weighted integrative analysis. We also assessed the scalability of these approaches by recording computational time and found that all these methods can be scalable to biobank-scale data although some might be relatively slow.

CONCLUSION

In conclusion, we hope that our findings can offer an important guidance on how to choose appropriate multilocus association analysis methods in post-GWAS era. All the SNP-set methods are implemented in the R package called MCA, which is freely available at https://github.com/biostatpzeng/ .

Integrative eQTL-weighted hierarchical Cox models for SNP-set based time-to-event association studies

BACKGROUND

Integrating functional annotations into SNP-set association studies has been proven a powerful analysis strategy. Statistical methods for such integration have been developed for continuous and binary phenotypes; however, the SNP-set integrative approaches for time-to-event or survival outcomes are lacking.

METHODS

We here propose IEHC, an integrative eQTL (expression quantitative trait loci) hierarchical Cox regression, for SNP-set based survival association analysis by modeling effect sizes of genetic variants as a function of eQTL via a hierarchical manner. Three p-values combination tests are developed to examine the joint effects of eQTL and genetic variants after a novel decorrelated modification of statistics for the two components. An omnibus test (IEHC-ACAT) is further adapted to aggregate the strengths of all available tests.

RESULTS

Simulations demonstrated that the IEHC joint tests were more powerful if both eQTL and genetic variants contributed to association signal, while IEHC-ACAT was robust and often outperformed other approaches across various simulation scenarios. When applying IEHC to ten TCGA cancers by incorporating eQTL from relevant tissues of GTEx, we revealed that substantial correlations existed between the two types of effect sizes of genetic variants from TCGA and GTEx, and identified 21 (9 unique) cancer-associated genes which would otherwise be missed by approaches not incorporating eQTL.

CONCLUSION

IEHC represents a flexible, robust, and powerful approach to integrate functional omics information to enhance the power of identifying association signals for the survival risk of complex human cancers.

IUSMMT: Survival mediation analysis of gene expression with multiple DNA methylation exposures and its application to cancers of TCGA

Effective and powerful survival mediation models are currently lacking. To partly fill such knowledge gap, we particularly focus on the mediation analysis that includes multiple DNA methylations acting as exposures, one gene expression as the mediator and one survival time as the outcome. We proposed IUSMMT (intersection-union survival mixture-adjusted mediation test) to effectively examine the existence of mediation effect by fitting an empirical three-component mixture null distribution. With extensive simulation studies, we demonstrated the advantage of IUSMMT over existing methods. We applied IUSMMT to ten TCGA cancers and identified multiple genes that exhibited mediating effects. We further revealed that most of the identified regions, in which genes behaved as active mediators, were cancer type-specific and exhibited a full mediation from DNA methylation CpG sites to the survival risk of various types of cancers. Overall, IUSMMT represents an effective and powerful alternative for survival mediation analysis; our results also provide new insights into the functional role of DNA methylation and gene expression in cancer progression/prognosis and demonstrate potential therapeutic targets for future clinical practice.

Sequence kernel association test for survival outcomes in the presence of a non-susceptible fraction

In this work, we propose a single nucleotide polymorphism set association test for survival phenotypes in the presence of a non-susceptible fraction. We consider a mixture model with a logistic regression for the susceptibility indicator and a proportional hazards regression to model survival in the susceptible group. We propose a joint test to assess the significance of the genetic variant in both logistic and survival regressions simultaneously. We adopt the spirit of SKAT and conduct a variance-component test treating the genetic effects of multiple variants as random. We derive score-type test statistics, and we investigate several approaches to compute their $p$-values. The finite-sample properties of the proposed tests are assessed and compared to existing approaches by simulations and their use is illustrated through an application to ovarian cancer data from the Consortium of Investigators of Modifiers of BRCA1 and BRCA2.

A Fast and Accurate Method for Genome-Wide Time-to-Event Data Analysis and Its Application to UK Biobank

With increasing biobanking efforts connecting electronic health records and national registries to germline genetics, the time-to-event data analysis has attracted increasing attention in the genetics studies of human diseases. In time-to-event data analysis, the Cox proportional hazards (PH) regression model is one of the most used approaches. However, existing methods and tools are not scalable when analyzing a large biobank with hundreds of thousands of samples and endpoints, and they are not accurate when testing low-frequency and rare variants. Here, we propose a scalable and accurate method, SPACox (a saddlepoint approximation implementation based on the Cox PH regression model), that is applicable for genome-wide scale time-to-event data analysis. SPACox requires fitting a Cox PH regression model only once across the genome-wide analysis and then uses a saddlepoint approximation (SPA) to calibrate the test statistics. Simulation studies show that SPACox is 76-252 times faster than other existing alternatives, such as gwasurvivr, 185-511 times faster than the standard Wald test, and more than 6,000 times faster than the Firth correction and can control type I error rates at the genome-wide significance level regardless of minor allele frequencies. Through the analysis of UK Biobank inpatient data of 282,871 white British European ancestry samples, we show that SPACox can efficiently analyze large sample sizes and accurately control type I error rates. We identified 611 loci associated with time-to-event phenotypes of 12 common diseases, of which 38 loci would be missed within a logistic regression framework with a binary phenotype defined as event occurrence status during the follow-up period.

Identifying individual risk rare variants using protein structure guided local tests (POINT)

Rare variants are of increasing interest to genetic association studies because of their etiological contributions to human complex diseases. Due to the rarity of the mutant events, rare variants are routinely analyzed on an aggregate level. While aggregation analyses improve the detection of global-level signal, they are not able to pinpoint causal variants within a variant set. To perform inference on a localized level, additional information, e.g., biological annotation, is often needed to boost the information content of a rare variant. Following the observation that important variants are likely to cluster together on functional domains, we propose a protein structure guided local test (POINT) to provide variant-specific association information using structure-guided aggregation of signal. Constructed under a kernel machine framework, POINT performs local association testing by borrowing information from neighboring variants in the 3-dimensional protein space in a data-adaptive fashion. Besides merely providing a list of promising variants, POINT assigns each variant a p-value to permit variant ranking and prioritization. We assess the selection performance of POINT using simulations and illustrate how it can be used to prioritize individual rare variants in PCSK9, ANGPTL4 and CETP in the Action to Control Cardiovascular Risk in Diabetes (ACCORD) clinical trial data.

Family-based association tests for rare variants with censored traits

We propose a set of family-based burden and kernel tests for censored traits (FamBAC and FamKAC). Here, censored traits refer to time-to-event outcomes, for instance, age-at-onset of a disease. To model censored traits in family-based designs, we used the frailty model, which incorporated not only fixed genetic effects of rare variants in a region of interest but also random polygenic effects shared within families. We first partitioned genotype scores of rare variants into orthogonal between- and within-family components, and then derived their corresponding efficient score statistics from the frailty model. Finally, FamBAC and FamKAC were constructed by aggregating the weighted efficient scores of the within-family components across rare variants and subjects. FamBAC collapsed rare variants within subject first to form a burden test that followed a chi-squared distribution; whereas FamKAC was a variant component test following a mixture of chi-squared distributions. For FamKAC, p-values can be computed by permutation tests or for computational efficiency by approximation methods. Through simulation studies, we showed that type I error was correctly controlled by FamBAC for various variant weighting schemes (0.0371 to 0.0527). However, FamKAC type I error rates based on approximation methods were deflated (max 0.0376) but improved by permutation tests. Our simulations also demonstrated that burden test FamBAC had higher power than kernel test FamKAC when high proportion (e.g. ≥ 80%) of causal variants had effects in the same direction. In contrast, when the effects of causal variants on the censored trait were in mixed directions, FamKAC outperformed FamBAC and had comparable or higher power than an existing method, RVFam. Our proposed framework has the flexibility to accommodate general nuclear families, and can be used to analyze sequence data for censored traits such as age-at-onset of a complex disease of interest.

A review of kernel methods for genetic association studies

Evaluating the association of multiple genetic variants with a trait of interest by use of kernel-based methods has made a significant impact on how genetic association analyses are conducted. An advantage of kernel methods is that they tend to be robust when the genetic variants have effects that are a mixture of positive and negative effects, as well as when there is a small fraction of causal variants. Another advantage is that kernel methods fit within the framework of mixed models, providing flexible ways to adjust for additional covariates that influence traits. Herein, we review the basic ideas behind the use of kernel methods for genetic association analysis as well as recent methodological advancements for different types of traits, multivariate traits, pedigree data, and longitudinal data. Finally, we discuss opportunities for future research.

An adaptive microbiome α-diversity-based association analysis method

To relate microbial diversity with various host traits of interest (e.g., phenotypes, clinical interventions, environmental factors) is a critical step for generic assessments about the disparity in human microbiota among different populations. The performance of the current item-by-item α-diversity-based association tests is sensitive to the choice of α-diversity metric and unpredictable due to the unknown nature of the true association. The approach of cherry-picking a test for the smallest p-value or the largest effect size among multiple item-by-item analyses is not even statistically valid due to the inherent multiplicity issue. Investigators have recently introduced microbial community-level association tests while blustering statistical power increase of their proposed methods. However, they are purely a test for significance which does not provide any estimation facilities on the effect direction and size of a microbial community; hence, they are not in practical use. Here, I introduce a novel microbial diversity association test, namely, adaptive microbiome α-diversity-based association analysis (aMiAD). aMiAD simultaneously tests the significance and estimates the effect score of the microbial diversity on a host trait, while robustly maintaining high statistical power and accurate estimation with no issues in validity.

Variance components genetic association test for zero-inflated count outcomes

Commonly in biomedical research, studies collect data in which an outcome measure contains informative excess zeros; for example, when observing the burden of neuritic plaques (NPs) in brain pathology studies, those who show none contribute to our understanding of neurodegenerative disease. The outcome may be characterized by a mixture distribution with one component being the "structural zero" and the other component being a Poisson distribution. We propose a novel variance components score test of genetic association between a set of genetic markers and a zero-inflated count outcome from a mixture distribution. This test shares advantageous properties with single-nucleotide polymorphism (SNP)-set tests which have been previously devised for standard continuous or binary outcomes, such as the sequence kernel association test. In particular, our method has superior statistical power compared to competing methods, especially when there is correlation within the group of markers, and when the SNPs are associated with both the mixing proportion and the rate of the Poisson distribution. We apply the method to Alzheimer's data from the Rush University Religious Orders Study and Memory and Aging Project, where as proof of principle we find highly significant associations with the APOE gene, in both the "structural zero" and "count" parameters, when applied to a zero-inflated NPs count outcome.

Semiparametric Bayesian kernel survival model for evaluating pathway effects

Massive amounts of high-dimensional data have been accumulated over the past two decades, which has cultured increasing interests in identifying gene pathways related to certain biological processes. In particular, since pathway-based analysis has the ability to detect subtle changes of differentially expressed genes that could be missed when using gene-based analysis, detecting the gene pathways that regulate certain diseases can provide new strategies for medical procedures and new targets for drug discovery. Limited work has been carried out, primarily in regression settings, to study the effects of pathways on survival outcomes. Motivated by a breast cancer gene-pathway data set, which exhibits the "small n, large p" characteristics, we propose a semiparametric Bayesian kernel survival model (s-BKSurv) to study the effects of both clinical covariates and gene expression levels within a pathway on survival time. We model the unknown high-dimensional functions of pathways via Gaussian kernel machine to consider the possibility that genes within the same pathway interact with each other. To address the multiple comparisons problem under a full Bayesian setting, we propose a similarity-dependent procedure based on Bayes factor to control the family-wise error rate. We demonstrate the outperformance of our approach under various simulation settings and pathways data.

TSGSIS: a high-dimensional grouped variable selection approach for detection of whole-genome SNP-SNP interactions

Motivation

Identification of single nucleotide polymorphism (SNP) interactions is an important and challenging topic in genome-wide association studies (GWAS). Many approaches have been applied to detecting whole-genome interactions. However, these approaches to interaction analysis tend to miss causal interaction effects when the individual marginal effects are uncorrelated to trait, while their interaction effects are highly associated with the trait.

Results

A grouped variable selection technique, called two-stage grouped sure independence screening (TS-GSIS), is developed to study interactions that may not have marginal effects. The proposed TS-GSIS is shown to be very helpful in identifying not only causal SNP effects that are uncorrelated to trait but also their corresponding SNP-SNP interaction effects. The benefit of TS-GSIS are gaining detection of interaction effects by taking the joint information among the SNPs and determining the size of candidate sets in the model. Simulation studies under various scenarios are performed to compare performance of TS-GSIS and current approaches. We also apply our approach to a real rheumatoid arthritis (RA) dataset. Both the simulation and real data studies show that the TS-GSIS performs very well in detecting SNP-SNP interactions.

Availability and implementation

R-package is delivered through CRAN and is available at: https://cran.r-project.org/web/packages/TSGSIS/index.html.

Contact

hsiung@nhri.org.tw.

Supplementary information

Supplementary data are available at Bioinformatics online.

Conditional asymptotic inference for the kernel association test

Motivation

The kernel association test (KAT) is popular in biological studies for its ability to combine weak effects potentially of opposite direction. Its P-value is typically assessed via its (unconditional) asymptotic distribution. However, such an asymptotic distribution is known only for continuous traits and for dichotomous traits. Furthermore, the derived P-values are known to be conservative when sample size is small, especially for the important case of dichotomous traits. One alternative is the permutation test, a widely accepted approximation to the exact finite sample conditional inference. But it is time-consuming to use in practice due to stringent significance criteria commonly seen in these analyses.

Results

Based on a previous theoretical result a conditional asymptotic distribution for the KAT is introduced. This distribution provides an alternative approximation to the exact distribution of the KAT. An explicit expression of this distribution is provided from which P-values can be easily computed. This method applies to any type of traits. The usefulness of this approach is demonstrated via extensive simulation studies using real genotype data and an analysis of genetic data from the Ocular Hypertension Treatment Study. Numerical results showed that the new method can control the type I error rate and is a bit conservative when compared to the permutation method. Nevertheless the proposed method may be used as a fast screening method. A time-consuming permutation procedure may be conducted at locations that show signals of association.

Availability and implementation

An implementation of the proposed method is provided in the R package iGasso.

Contact

kai-wang@uiowa.edu.

Variance component score test for time-course gene set analysis of longitudinal RNA-seq data

As gene expression measurement technology is shifting from microarrays to sequencing, the statistical tools available for their analysis must be adapted since RNA-seq data are measured as counts. It has been proposed to model RNA-seq counts as continuous variables using nonparametric regression to account for their inherent heteroscedasticity. In this vein, we propose tcgsaseq, a principled, model-free, and efficient method for detecting longitudinal changes in RNA-seq gene sets defined a priori. The method identifies those gene sets whose expression varies over time, based on an original variance component score test accounting for both covariates and heteroscedasticity without assuming any specific parametric distribution for the (transformed) counts. We demonstrate that despite the presence of a nonparametric component, our test statistic has a simple form and limiting distribution, and both may be computed quickly. A permutation version of the test is additionally proposed for very small sample sizes. Applied to both simulated data and two real datasets, tcgsaseq is shown to exhibit very good statistical properties, with an increase in stability and power when compared to state-of-the-art methods ROAST (rotation gene set testing), edgeR, and DESeq2, which can fail to control the type I error under certain realistic settings. We have made the method available for the community in the R package tcgsaseq.

A highly adaptive microbiome-based association test for survival traits

BACKGROUND

There has been increasing interest in discovering microbial taxa that are associated with human health or disease, gathering momentum through the advances in next-generation sequencing technologies. Investigators have also increasingly employed prospective study designs to survey survival (i.e., time-to-event) outcomes, but current item-by-item statistical methods have limitations due to the unknown true association pattern. Here, we propose a new adaptive microbiome-based association test for survival outcomes, namely, optimal microbiome-based survival analysis (OMiSA). OMiSA approximates to the most powerful association test in two domains: 1) microbiome-based survival analysis using linear and non-linear bases of OTUs (MiSALN) which weighs rare, mid-abundant, and abundant OTUs, respectively, and 2) microbiome regression-based kernel association test for survival traits (MiRKAT-S) which incorporates different distance metrics (e.g., unique fraction (UniFrac) distance and Bray-Curtis dissimilarity), respectively.

RESULTS

We illustrate that OMiSA powerfully discovers microbial taxa whether their underlying associated lineages are rare or abundant and phylogenetically related or not. OMiSA is a semi-parametric method based on a variance-component score test and a re-sampling method; hence, it is free from any distributional assumption on the effect of microbial composition and advantageous to robustly control type I error rates. Our extensive simulations demonstrate the highly robust performance of OMiSA. We also present the use of OMiSA with real data applications.

CONCLUSIONS

OMiSA is attractive in practice as the true association pattern is unpredictable in advance and, for survival outcomes, no adaptive microbiome-based association test is currently available.

Germline polymorphisms in myeloid-associated genes are not associated with survival in glioma patients

Immune cells of myeloid origin, including microglia, macrophages, and myeloid-derived suppressor cells adopt immunosuppressive phenotypes that support gliomagenesis. Here, we tested an a priori hypothesis that single nucleotide polymorphisms (SNPs) in genes related to glioma-associated myeloid cell regulation and function are also associated with patient survival after glioma diagnosis. Subjects for this study were 992 glioma patients treated at The University of Texas MD Anderson Cancer Center in Houston, Texas between 1992 and 2008. Haplotype-tagging SNPs in 91 myeloid-associated genes were analyzed for association with survival by Cox regression. Individual SNP- and gene-based tests were performed separately in glioblastoma (WHO grade IV, n = 511) and lower-grade glioma (WHO grade II-III, n = 481) groups. After adjustment for multiple testing, no myeloid-associated gene variants were significantly associated with survival in glioblastoma. Two SNPs, rs147960238 in CD163 (p = 2.2 × 10-5) and rs17138945 in MET (p = 5.6 × 10-5) were significantly associated with survival of patients with lower-grade glioma. However, these associations were not confirmed in an independent analysis of 563 lower-grade glioma cases from the University of California at San Francisco Adult Glioma Study (p = 0.65 and p = 0.41, respectively). The results of this study do not support a role for inherited polymorphisms in myeloid-associated genes in affecting survival of patients diagnosed with glioblastoma or lower-grade glioma.

Assessing Rare Variation in Complex Traits

While genome-wide association studies have been very successful in identifying associations of common genetic variants with many different traits, the rarer frequency spectrum of the genome has not yet been comprehensively explored. Technological developments increasingly lift restrictions to access rare genetic variation. Dense reference panels enable improved genotype imputation for rarer variants in studies using DNA microarrays. Moreover, the decreasing cost of next generation sequencing makes whole exome and genome sequencing increasingly affordable for large samples. Large-scale efforts based on sequencing, such as ExAC, 100,000 Genomes, and TopMed, are likely to significantly advance this field.The main challenge in evaluating complex trait associations of rare variants is statistical power. The choice of population should be considered carefully because allele frequencies and linkage disequilibrium structure differ between populations. Genetically isolated populations can have favorable genomic characteristics for the study of rare variants.One strategy to increase power is to assess the combined effect of multiple rare variants within a region, known as aggregate testing. A  range of methods have been developed for this. Model performance depends on the genetic architecture of the region of interest.

Interactions Between Genome-Wide Significant Genetic Variants and Circulating Concentrations of 25-Hydroxyvitamin D in Relation to Prostate Cancer Risk in the National Cancer Institute BPC3

Genome-wide association studies (GWAS) have identified over 100 single nucleotide polymorphisms (SNPs) associated with prostate cancer. However, information on the mechanistic basis for some associations is limited. Recent research has been directed towards the potential association of vitamin D concentrations and prostate cancer, but little is known about whether the aforementioned genetic associations are modified by vitamin D. We investigated the associations of 46 GWAS-identified SNPs, circulating concentrations of 25-hydroxyvitamin D (25(OH)D), and prostate cancer (3,811 cases, 511 of whom died from the disease, compared with 2,980 controls-from 5 cohort studies that recruited participants over several periods beginning in the 1980s). We used logistic regression models with data from the National Cancer Institute Breast and Prostate Cancer Cohort Consortium (BPC3) to evaluate interactions on the multiplicative and additive scales. After allowing for multiple testing, none of the SNPs examined was significantly associated with 25(OH)D concentration, and the SNP-prostate cancer associations did not differ by these concentrations. A statistically significant interaction was observed for each of 2 SNPs in the 8q24 region (rs620861 and rs16902094), 25(OH)D concentration, and fatal prostate cancer on both multiplicative and additive scales (P ≤ 0.001). We did not find strong evidence that associations between GWAS-identified SNPs and prostate cancer are modified by circulating concentrations of 25(OH)D. The intriguing interactions between rs620861 and rs16902094, 25(OH)D concentration, and fatal prostate cancer warrant replication.

MiRKAT-S: a community-level test of association between the microbiota and survival times

BACKGROUND

Community-level analysis of the human microbiota has culminated in the discovery of relationships between overall shifts in the microbiota and a wide range of diseases and conditions. However, existing work has primarily focused on analysis of relatively simple dichotomous or quantitative outcomes, for example, disease status or biomarker levels. Recently, there is also considerable interest in the relationship between the microbiota and censored survival outcomes, such as in clinical trials. How to conduct community-level analysis with censored survival outcomes is unclear, since standard dissimilarity-based tests cannot accommodate censored survival times and no alternative methods exist.

METHODS

We develop a new approach, MiRKAT-S, for community-level analysis of microbiome data with censored survival times. MiRKAT-S uses ecologically informative distance metrics, such as the UniFrac distances, to generate matrices of pairwise distances between individuals' taxonomic profiles. The distance matrices are transformed into kernel (similarity) matrices, which are used to compare similarity in the microbiota to similarity in survival times between individuals.

RESULTS

Simulation studies using synthetic microbial communities demonstrate correct control of type I error and adequate power. We also apply MiRKAT-S to examine the relationship between the gut microbiota and survival after allogeneic blood or bone marrow transplant.

CONCLUSIONS

We present MiRKAT-S, a method that facilitates community-level analysis of the association between the microbiota and survival outcomes and therefore provides a new approach to analysis of microbiome data arising from clinical trials.

Exome-Wide Association Study Identifies Low-Frequency Coding Variants in 2p23.2 and 7p11.2 Associated with Survival of Non-Small Cell Lung Cancer Patients

INTRODUCTION

A growing body of evidence has suggested that low-frequency or rare coding variants might have strong effects on the development and prognosis of cancer. Here, we aim to assess the role of low-frequency and rare coding variants in the survival of NSCLC in Chinese populations.

METHODS

We performed an exome-wide scan of 247,870 variants in 1008 patients with NSCLC and replicated the promising variants by using imputed genotype data of The Cancer Genome Atlas (TCGA) with a Cox regression model. Gene-based and pathway-based analysis were also performed for nonsynonymous or splice site variants. Additionally, analysis of gene expression data in the TCGA was used to increase the reliability of candidate loci and genes.

RESULTS

A low-frequency missense variant in chaperonin containing TCP1 subunit 6A gene (CCT6A) (rs33922584: adjusted hazard ratio [HR_adjusted] = 1.75, p = 6.06 × 10^-4) was significantly related to the survival of patients with NSCLC, which was further replicated by the TCGA samples (HR_adjusted = 4.19, p = 0.015). Interestingly, the G allele of rs33922584 was significantly associated with high expression of CCT6A (p = 0.019) that might induce the worse survival in the TCGA samples (HR_adjusted = 1.15, p = 0.047). Besides, rs117512489 in gene phospholipase B1 gene (PLB1) (HR = 2.02, p = 7.28 × 10^-4) was also associated with survival of the patients with NSCLC in our samples, but it was supported only by gene expression analysis in the TCGA (HR_adjusted = 1.15, p = 0.023). Gene-based and pathway-based analysis revealed a total of 32 genes, including CCT6A and 34 potential pathways might account for the survival of NSCLC, respectively.

CONCLUSION

These results provided more evidence for the important role of low-frequency or rare variants in the survival of patients with NSCLC.

Prioritizing individual genetic variants after kernel machine testing using variable selection

Kernel machine learning methods, such as the SNP-set kernel association test (SKAT), have been widely used to test associations between traits and genetic polymorphisms. In contrast to traditional single-SNP analysis methods, these methods are designed to examine the joint effect of a set of related SNPs (such as a group of SNPs within a gene or a pathway) and are able to identify sets of SNPs that are associated with the trait of interest. However, as with many multi-SNP testing approaches, kernel machine testing can draw conclusion only at the SNP-set level, and does not directly inform on which one(s) of the identified SNP set is actually driving the associations. A recently proposed procedure, KerNel Iterative Feature Extraction (KNIFE), provides a general framework for incorporating variable selection into kernel machine methods. In this article, we focus on quantitative traits and relatively common SNPs, and adapt the KNIFE procedure to genetic association studies and propose an approach to identify driver SNPs after the application of SKAT to gene set analysis. Our approach accommodates several kernels that are widely used in SNP analysis, such as the linear kernel and the Identity by State (IBS) kernel. The proposed approach provides practically useful utilities to prioritize SNPs, and fills the gap between SNP set analysis and biological functional studies. Both simulation studies and real data application are used to demonstrate the proposed approach.

Boosting the Power of the Sequence Kernel Association Test by Properly Estimating Its Null Distribution

The sequence kernel association test (SKAT) is probably the most popular statistical test used in rare-variant association studies. Its null distribution involves unknown parameters that need to be estimated. The current estimation method has a valid type I error rate, but the power is compromised given that all subjects are used for estimation. I have developed an estimation method that uses only control subjects. Named SKAT+, this method uses the same test statistic as SKAT but differs in the way the null distribution is estimated. Extensive simulation studies and applications to data from the Genetic Analysis Workshop 17 and the Ocular Hypertension Treatment Study demonstrated that SKAT+ has superior power over SKAT while maintaining control over the type I error rate. This method is applicable to extensions of SKAT in the literature.

Predicting Alzheimer's Disease Using Combined Imaging-Whole Genome SNP Data

The growing public threat of Alzheimer's disease (AD) has raised the urgency to discover and validate prognostic biomarkers in order to predicting time to onset of AD. It is anticipated that both whole genome single nucleotide polymorphism (SNP) data and high dimensional whole brain imaging data offer predictive values to identify subjects at risk for progressing to AD. The aim of this paper is to test whether both whole genome SNP data and whole brain imaging data offer predictive values to identify subjects at risk for progressing to AD. In 343 subjects with mild cognitive impairment (MCI) enrolled in the Alzheimer's Disease Neuroimaging Initiative (ADNI-1), we extracted high dimensional MR imaging (volumetric data on 93 brain regions plus a surface fluid registration based hippocampal subregion and surface data), and whole genome data (504,095 SNPs from GWAS), as well as routine neurocognitive and clinical data at baseline. MCI patients were then followed over 48 months, with 150 participants progressing to AD. Combining information from whole brain MR imaging and whole genome data was substantially superior to the standard model for predicting time to onset of AD in a 48-month national study of subjects at risk. Our findings demonstrate the promise of combined imaging-whole genome prognostic markers in people with mild memory impairment.

Kernel machine score test for pathway analysis in the presence of semi-competing risks

In cancer studies, patients often experience two different types of events: a non-terminal event such as recurrence or metastasis, and a terminal event such as cancer-specific death. Identifying pathways and networks of genes associated with one or both of these events is an important step in understanding disease development and targeting new biological processes for potential intervention. These correlated outcomes are commonly dealt with by modeling progression-free survival, where the event time is the minimum between the times of recurrence and death. However, identifying pathways only associated with progression-free survival may miss out on pathways that affect time to recurrence but not death, or vice versa. We propose a combined testing procedure for a pathway's association with both the cause-specific hazard of recurrence and the marginal hazard of death. The dependency between the two outcomes is accounted for through perturbation resampling to approximate the test's null distribution, without any further assumption on the nature of the dependency. Even complex non-linear relationships between pathways and disease progression or death can be uncovered thanks to a flexible kernel machine framework. The superior statistical power of our approach is demonstrated in numerical studies and in a gene expression study of breast cancer.

On Sample Size and Power Calculation for Variant Set-Based Association Tests

Sample size and power calculations are an important part of designing new sequence-based association studies. The recently developed SEQPower and SPS programs adopted computationally intensive Monte Carlo simulations to empirically estimate power for a series of variant set association (VSA) test methods including the sequence kernel association test (SKAT). It is desirable to develop methods that can quickly and accurately compute power without intensive Monte Carlo simulations. We will show that the computed power for SKAT based on the existing analytical approach could be inflated especially for small significance levels, which are often of primary interest for large-scale whole genome and exome sequencing projects. We propose a new χ(2) -approximation-based approach to accurately and efficiently compute sample size and power. In addition, we propose and implement a more accurate "exact" method to compute power, which is more efficient than the Monte Carlo approach though generally involves more computations than the χ(2) approximation method. The exact approach could produce very accurate results and be used to verify alternative approximation approaches. We implement the proposed methods in publicly available R programs that can be readily adapted when planning sequencing projects.

Gene-Based Association Analysis for Censored Traits Via Fixed Effect Functional Regressions

Genetic studies of survival outcomes have been proposed and conducted recently, but statistical methods for identifying genetic variants that affect disease progression are rarely developed. Motivated by our ongoing real studies, here we develop Cox proportional hazard models using functional regression (FR) to perform gene-based association analysis of survival traits while adjusting for covariates. The proposed Cox models are fixed effect models where the genetic effects of multiple genetic variants are assumed to be fixed. We introduce likelihood ratio test (LRT) statistics to test for associations between the survival traits and multiple genetic variants in a genetic region. Extensive simulation studies demonstrate that the proposed Cox RF LRT statistics have well-controlled type I error rates. To evaluate power, we compare the Cox FR LRT with the previously developed burden test (BT) in a Cox model and sequence kernel association test (SKAT), which is based on mixed effect Cox models. The Cox FR LRT statistics have higher power than or similar power as Cox SKAT LRT except when 50%/50% causal variants had negative/positive effects and all causal variants are rare. In addition, the Cox FR LRT statistics have higher power than Cox BT LRT. The models and related test statistics can be useful in the whole genome and whole exome association studies. An age-related macular degeneration dataset was analyzed as an example.

Calcium intake, polymorphisms of the calcium-sensing receptor, and recurrent/aggressive prostate cancer

PURPOSE

To assess whether calcium intake and common genetic variants of the calcium-sensing receptor (CASR) are associated with either aggressive prostate cancer (PCa) or disease recurrence after prostatectomy.

METHODS

Calcium intake at diagnosis was assessed, and 65 common single-nucleotide polymorphisms (SNPs) in CASR were genotyped in 886 prostatectomy patients. We investigated the association between calcium intake and CASR variants with both PCa recurrence and aggressiveness (defined as Gleason score ≥4 + 3, stage ≥pT3, or nodal-positive disease).

RESULTS

A total of 285 men had aggressive disease and 91 experienced recurrence. A U-shaped relationship between calcium intake and both disease recurrence and aggressiveness was observed. Compared to the middle quintile, the HR for disease recurrence was 3.07 (95% CI 1.41-6.69) for the lowest quintile and 3.21 (95% CI 1.47-7.00) and 2.97 (95% CI 1.37-6.45) for the two upper quintiles, respectively. Compared to the middle quintile, the OR for aggressive disease was 1.80 (95% CI 1.11-2.91) for the lowest quintile and 1.75 (95% CI 1.08-2.85) for the highest quintile of calcium intake. The main effects of CASR variants were not associated with PCa recurrence or aggressiveness. In the subgroup of patients with moderate calcium intake, 31 SNPs in four distinct blocks of high linkage disequilibrium were associated with PCa recurrence.

CONCLUSIONS

We observed a protective effect of moderate calcium intake for PCa aggressiveness and recurrence. While CASR variants were not associated with these outcomes in the entire cohort, they may be associated with disease recurrence in men with moderate calcium intakes.

SNP Set Association Testing for Survival Outcomes in the Presence of Intrafamilial Correlation

In this work, we propose a single nucleotide polymorphism (SNP) set association test for censored phenotypes in the presence of a family-based design. The proposed test is valid for both common and rare variants. A proportional hazards Cox model is specified for the marginal distribution of the trait and the familial dependence is modeled via a Gaussian copula. Censored values are treated as partially missing data and a multiple imputation procedure is proposed in order to compute the test statistics. The P-value is then deduced analytically. The finite-sample empirical properties of the proposed method are evaluated and compared to existing competitors by simulations and its use is illustrated using a breast cancer data set from the Consortium of Investigators of Modifiers of BRCA1 and BRCA2.

A Fast Multiple-Kernel Method With Applications to Detect Gene-Environment Interaction

Kernel machine (KM) models are a powerful tool for exploring associations between sets of genetic variants and complex traits. Although most KM methods use a single kernel function to assess the marginal effect of a variable set, KM analyses involving multiple kernels have become increasingly popular. Multikernel analysis allows researchers to study more complex problems, such as assessing gene-gene or gene-environment interactions, incorporating variance-component based methods for population substructure into rare-variant association testing, and assessing the conditional effects of a variable set adjusting for other variable sets. The KM framework is robust, powerful, and provides efficient dimension reduction for multifactor analyses, but requires the estimation of high dimensional nuisance parameters. Traditional estimation techniques, including regularization and the "expectation-maximization (EM)" algorithm, have a large computational cost and are not scalable to large sample sizes needed for rare variant analysis. Therefore, under the context of gene-environment interaction, we propose a computationally efficient and statistically rigorous "fastKM" algorithm for multikernel analysis that is based on a low-rank approximation to the nuisance effect kernel matrices. Our algorithm is applicable to various trait types (e.g., continuous, binary, and survival traits) and can be implemented using any existing single-kernel analysis software. Through extensive simulation studies, we show that our algorithm has similar performance to an EM-based KM approach for quantitative traits while running much faster. We also apply our method to the Vitamin Intervention for Stroke Prevention (VISP) clinical trial, examining gene-by-vitamin effects on recurrent stroke risk and gene-by-age effects on change in homocysteine level.

Circulating vitamin D, vitamin D-related genetic variation, and risk of fatal prostate cancer in the National Cancer Institute Breast and Prostate Cancer Cohort Consortium

BACKGROUND

Evidence from experimental animal and cell line studies supports a beneficial role for vitamin D in prostate cancer (PCa). Although the results from human studies have been mainly null for overall PCa risk, there may be a benefit for survival. This study assessed the associations of circulating 25-hydroxyvitamin D (25(OH)D) and common variations in key vitamin D-related genes with fatal PCa.

METHODS

In a large cohort consortium, 518 fatal cases and 2986 controls with 25(OH)D data were identified. Genotyping information for 91 single-nucleotide polymorphisms (SNPs) in 7 vitamin D-related genes (vitamin D receptor, group-specific component, cytochrome P450 27A1 [CYP27A1], CYP27B1, CYP24A1, CYP2R1, and retinoid X receptor α) was available for 496 fatal cases and 3577 controls. Unconditional logistic regression was used to calculate odds ratios (ORs) and 95% confidence intervals (CIs) for the associations of 25(OH)D and SNPs with fatal PCa. The study also tested for 25(OH)D-SNP interactions among 264 fatal cases and 1169 controls.

RESULTS

No statistically significant relationship was observed between 25(OH)D and fatal PCa (OR for extreme quartiles, 0.86; 95% CI, 0.65-1.14; P for trend = .22) or the main effects of the SNPs and fatal PCa. There was evidence suggesting that associations of several SNPs, including 5 related to circulating 25(OH)D, with fatal PCa were modified by 25(OH)D. Individually, these associations did not remain significant after multiple testing; however, the P value for the set-based test for CYP2R1 was .002.

CONCLUSIONS

Statistically significant associations were not observed for either 25(OH)D or vitamin D-related SNPs with fatal PCa. The effect modification of 25(OH)D associations by biologically plausible genetic variation may deserve further exploration.

Circadian clock genes and risk of fatal prostate cancer

PURPOSE

Circadian genes may be involved in regulating cancer-related pathways, including cell proliferation, DNA damage response, and apoptosis. We aimed to assess the role of genetic variation in core circadian rhythm genes with the risk of fatal prostate cancer and first morning void urinary 6-sulfatoxymelatonin levels.

METHODS

We used unconditional logistic regression to evaluate the association of 96 single-nucleotide polymorphisms (SNPs) across 12 circadian-related genes with fatal prostate cancer in the AGES-Reykjavik cohort (n = 24 cases), the Health Professionals Follow-Up Study (HPFS) (n = 40 cases), and the Physicians' Health Study (PHS) (n = 105 cases). We used linear regression to evaluate the association between SNPs and first morning void urinary 6-sulfatoxymelatonin levels in AGES-Reykjavik. We used a kernel machine test to evaluate whether multimarker SNP sets in the pathway (gene based) were associated with our outcomes.

RESULTS

None of the individual SNPs were consistently associated with fatal prostate cancer across the three cohorts. In each cohort, gene-based analyses showed that variation in the CRY1 gene was nominally associated with fatal prostate cancer (p values = 0.01, 0.01, and 0.05 for AGES-Reykjavik, HPFS, and PHS, respectively). In AGES-Reykjavik, SNPs in TIMELESS (four SNPs), NPAS2 (six SNPs), PER3 (two SNPs) and CSNK1E (one SNP) were nominally associated with 6-sulfatoxymelatonin levels.

CONCLUSION

We did not find a strong and consistent association between variation in core circadian clock genes and fatal prostate cancer risk, but observed nominally significant gene-based associations with fatal prostate cancer and 6-sulfatoxymelatonin levels.

Global analysis of methylation profiles from high resolution CpG data

New high throughput technologies are now enabling simultaneous epigenetic profiling of DNA methylation at hundreds of thousands of CpGs across the genome. A problem of considerable practical interest is identification of large scale, global changes in methylation that are associated with environmental variables, clinical outcomes, or other experimental conditions. However, there has been little statistical research on methods for global methylation analysis using technologies with individual CpG resolution. To address this critical gap in the literature, we develop a new strategy for global analysis of methylation profiles using a functional regression approach wherein we approximate either the density or the cumulative distribution function (CDF) of the methylation values for each individual using B-spline basis functions. The spline coefficients for each individual are allowed to summarize the individual's overall methylation profile. We then test for association between the overall distribution and a continuous or dichotomous outcome variable using a variance component score test that naturally accommodates the correlation between spline coefficients. Simulations indicate that our proposed approach has desirable power while protecting type I error. The method was applied to detect methylation differences, both genome wide and at LINE1 elements, between the blood samples from rheumatoid arthritis patients and healthy controls and to detect the epigenetic changes of human hepatocarcinogenesis in the context of alcohol abuse and hepatitis C virus infection. A free implementation of our methods in the R language is available in the Global Analysis of Methylation Profiles (GAMP) package at http://research.fhcrc.org/wu/en.html.

A non-parametric approach for detecting gene-gene interactions associated with age-at-onset outcomes

Background

Cox-regression-based methods have been commonly used for the analyses of survival outcomes, such as age-at-disease-onset. These methods generally assume the hazard functions are proportional among various risk groups. However, such an assumption may not be valid in genetic association studies, especially when complex interactions are involved. In addition, genetic association studies commonly adopt case-control designs. Direct use of Cox regression to case-control data may yield biased estimators and incorrect statistical inference.

Results

We propose a non-parametric approach, the weighted Nelson-Aalen (WNA) approach, for detecting genetic variants that are associated with age-dependent outcomes. The proposed approach can be directly applied to prospective cohort studies, and can be easily extended for population-based case-control studies. Moreover, it does not rely on any assumptions of the disease inheritance models, and is able to capture high-order gene-gene interactions. Through simulations, we show the proposed approach outperforms Cox-regression-based methods in various scenarios. We also conduct an empirical study of progression of nicotine dependence by applying the WNA approach to three independent datasets from the Study of Addiction: Genetics and Environment. In the initial dataset, two SNPs, rs6570989 and rs2930357, located in genes GRIK2 and CSMD1, are found to be significantly associated with the progression of nicotine dependence (ND). The joint association is further replicated in two independent datasets. Further analysis suggests that these two genes may interact and be associated with the progression of ND.

Conclusions

As demonstrated by the simulation studies and real data analysis, the proposed approach provides an efficient tool for detecting genetic interactions associated with age-at-onset outcomes.

GENE-LEVEL PHARMACOGENETIC ANALYSIS ON SURVIVAL OUTCOMES USING GENE-TRAIT SIMILARITY REGRESSION

Gene/pathway-based methods are drawing significant attention due to their usefulness in detecting rare and common variants that affect disease susceptibility. The biological mechanism of drug responses indicates that a gene-based analysis has even greater potential in pharmacogenetics. Motivated by a study from the Vitamin Intervention for Stroke Prevention (VISP) trial, we develop a gene-trait similarity regression for survival analysis to assess the effect of a gene or pathway on time-to-event outcomes. The similarity regression has a general framework that covers a range of survival models, such as the proportional hazards model and the proportional odds model. The inference procedure developed under the proportional hazards model is robust against model misspecification. We derive the equivalence between the similarity survival regression and a random effects model, which further unifies the current variance-component based methods. We demonstrate the effectiveness of the proposed method through simulation studies. In addition, we apply the method to the VISP trial data to identify the genes that exhibit an association with the risk of a recurrent stroke. TCN2 gene was found to be associated with the recurrent stroke risk in the low-dose arm. This gene may impact recurrent stroke risk in response to cofactor therapy.

Prostate cancer (PCa) risk variants and risk of fatal PCa in the National Cancer Institute Breast and Prostate Cancer Cohort Consortium

BACKGROUND

Screening and diagnosis of prostate cancer (PCa) is hampered by an inability to predict who has the potential to develop fatal disease and who has indolent cancer. Studies have identified multiple genetic risk loci for PCa incidence, but it is unknown whether they could be used as biomarkers for PCa-specific mortality (PCSM).

OBJECTIVE

To examine the association of 47 established PCa risk single-nucleotide polymorphisms (SNPs) with PCSM.

DESIGN, SETTING, AND PARTICIPANTS

We included 10 487 men who had PCa and 11 024 controls, with a median follow-up of 8.3 yr, during which 1053 PCa deaths occurred.

OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS

The main outcome was PCSM. The risk allele was defined as the allele associated with an increased risk for PCa in the literature. We used Cox proportional hazards regression to calculate the hazard ratios of each SNP with time to progression to PCSM after diagnosis. We also used logistic regression to calculate odds ratios for each risk SNP, comparing fatal PCa cases to controls.

RESULTS AND LIMITATIONS

Among the cases, we found that 8 of the 47 SNPs were significantly associated (p<0.05) with time to PCSM. The risk allele of rs11672691 (intergenic) was associated with an increased risk for PCSM, while 7 SNPs had risk alleles inversely associated (rs13385191 [C2orf43], rs17021918 [PDLIM5], rs10486567 [JAZF1], rs6465657 [LMTK2], rs7127900 (intergenic), rs2735839 [KLK3], rs10993994 [MSMB], rs13385191 [C2orf43]). In the case-control analysis, 22 SNPs were associated (p<0.05) with the risk of fatal PCa, but most did not differentiate between fatal and nonfatal PCa. Rs11672691 and rs10993994 were associated with both fatal and nonfatal PCa, while rs6465657, rs7127900, rs2735839, and rs13385191 were associated with nonfatal PCa only.

CONCLUSIONS

Eight established risk loci were associated with progression to PCSM after diagnosis. Twenty-two SNPs were associated with fatal PCa incidence, but most did not differentiate between fatal and nonfatal PCa. The relatively small magnitudes of the associations do not translate well into risk prediction, but these findings merit further follow-up, because they may yield important clues about the complex biology of fatal PCa.

PATIENT SUMMARY

In this report, we assessed whether established PCa risk variants could predict PCSM. We found eight risk variants associated with PCSM: One predicted an increased risk of PCSM, while seven were associated with decreased risk. Larger studies that focus on fatal PCa are needed to identify more markers that could aid prediction.

Insulin-like growth factor pathway genetic polymorphisms, circulating IGF1 and IGFBP3, and prostate cancer survival

BACKGROUND

The insulin-like growth factor (IGF) signaling pathway has been implicated in prostate cancer (PCa) initiation, but its role in progression remains unknown.

METHODS

Among 5887 PCa patients (704 PCa deaths) of European ancestry from seven cohorts in the National Cancer Institute Breast and Prostate Cancer Cohort Consortium, we conducted Cox kernel machine pathway analysis to evaluate whether 530 tagging single nucleotide polymorphisms (SNPs) in 26 IGF pathway-related genes were collectively associated with PCa mortality. We also conducted SNP-specific analysis using stratified Cox models adjusting for multiple testing. In 2424 patients (313 PCa deaths), we evaluated the association of prediagnostic circulating IGF1 and IGFBP3 levels and PCa mortality. All statistical tests were two-sided.

RESULTS

The IGF signaling pathway was associated with PCa mortality (P = .03), and IGF2-AS and SSTR2 were the main contributors (both P = .04). In SNP-specific analysis, 36 SNPs were associated with PCa mortality with P trend less than .05, but only three SNPs in the IGF2-AS remained statistically significant after gene-based corrections. Two were in linkage disequilibrium (r 2 = 1 for rs1004446 and rs3741211), whereas the third, rs4366464, was independent (r 2 = 0.03). The hazard ratios (HRs) per each additional risk allele were 1.19 (95% confidence interval [CI] = 1.06 to 1.34; P trend = .003) for rs3741211 and 1.44 (95% CI = 1.20 to 1.73; P trend < .001) for rs4366464. rs4366464 remained statistically significant after correction for all SNPs (P trend.corr = .04). Prediagnostic IGF1 (HRhighest vs lowest quartile = 0.71; 95% CI = 0.48 to 1.04) and IGFBP3 (HR = 0.93; 95% CI = 0.65 to 1.34) levels were not associated with PCa mortality.

CONCLUSIONS

The IGF signaling pathway, primarily IGF2-AS and SSTR2 genes, may be important in PCa survival.

Insulin-like Growth Factor Pathway Genetic Polymorphisms, Circulating IGF1 and IGFBP3, and Prostate Cancer Survival

Background

The insulin-like growth factor (IGF) signaling pathway has been implicated in prostate cancer (PCa) initiation, but its role in progression remains unknown.

Methods

Among 5887 PCa patients (704 PCa deaths) of European ancestry from seven cohorts in the National Cancer Institute Breast and Prostate Cancer Cohort Consortium, we conducted Cox kernel machine pathway analysis to evaluate whether 530 tagging single nucleotide polymorphisms (SNPs) in 26 IGF pathway-related genes were collectively associated with PCa mortality. We also conducted SNP-specific analysis using stratified Cox models adjusting for multiple testing. In 2424 patients (313 PCa deaths), we evaluated the association of prediagnostic circulating IGF1 and IGFBP3 levels and PCa mortality. All statistical tests were two-sided.

Results

The IGF signaling pathway was associated with PCa mortality (P = .03), and IGF2-AS and SSTR2 were the main contributors (both P = .04). In SNP-specific analysis, 36 SNPs were associated with PCa mortality with P trend less than .05, but only three SNPs in the IGF2-AS remained statistically significant after gene-based corrections. Two were in linkage disequilibrium (r 2 = 1 for rs1004446 and rs3741211), whereas the third, rs4366464, was independent (r 2 = 0.03). The hazard ratios (HRs) per each additional risk allele were 1.19 (95% confidence interval [CI] = 1.06 to 1.34; P trend = .003) for rs3741211 and 1.44 (95% CI = 1.20 to 1.73; P trend < .001) for rs4366464. rs4366464 remained statistically significant after correction for all SNPs (P trend.corr = .04). Prediagnostic IGF1 (HRhighest vs lowest quartile = 0.71; 95% CI = 0.48 to 1.04) and IGFBP3 (HR = 0.93; 95% CI = 0.65 to 1.34) levels were not associated with PCa mortality.

Conclusions

The IGF signaling pathway, primarily IGF2-AS and SSTR2 genes, may be important in PCa survival.

Sequence kernel association test for survival traits

Rare variant tests have been of great interest in testing genetic associations with diseases and disease-related quantitative traits in recent years. Among these tests, the sequence kernel association test (SKAT) is an omnibus test for effects of rare genetic variants, in a linear or logistic regression framework. It is often described as a variance component test treating the genotypic effects as random. When the linear kernel is used, its test statistic can be expressed as a weighted sum of single-marker score test statistics. In this paper, we extend the test to survival phenotypes in a Cox regression framework. Because of the anticonservative small-sample performance of the score test in a Cox model, we substitute signed square-root likelihood ratio statistics for the score statistics, and confirm that the small-sample control of type I error is greatly improved. This test can also be applied in meta-analysis. We show in our simulation studies that this test has superior statistical power except in a few specific scenarios, as compared to burden tests in a Cox model. We also present results in an application to time-to-obesity using genotypes from Framingham Heart Study SNP Health Association Resource.

Adiponectin pathway polymorphisms and risk of breast cancer in African Americans and Hispanics in the Women's Health Initiative

Adiponectin, a protein secreted by the adipose tissue, is an endogenous insulin sensitizer with circulating levels that are decreased in obese and diabetic subjects. Recently, circulating levels of adiponectin have been correlated with breast cancer risk. Our previous work showed that polymorphisms of the adiponectin pathway are associated with breast cancer risk. We conducted the first study of adiponectin pathways in African Americans and Hispanics in the Women's Health Initiative SNP Health Association Resource cohort of 3,642 self-identified Hispanic women and 8,515 self-identified African American women who provided consent for DNA analysis. Single nucleotide polymorphisms (SNPs) from three genes were included in this analysis: ADIPOQ, ADIPOR1, and ADIPOR2. The genome-wide human SNP array 6.0 (909,622 SNPs) ( www.affymetrix.com ) was used. We found that rs1501299, a functional SNP of ADIPOQ that we previously reported was associated with breast cancer risk in a mostly Caucasian population, was also significantly associated with breast cancer incidence (HR for the GG/TG genotype: 1.23; 95 % CI 1.059-1.43) in African American women. We did not find any other SNPs in these genes to be associated with breast cancer incidence. This is the first study assessing the role of adiponectin pathway SNPs in breast cancer risk in African Americans and Hispanics. RS1501299 is significantly associated with breast cancer risk in African American women. As the rates of obesity and diabetes increase in African Americans and Hispanics, adiponectin and its functional SNPs may aid in breast cancer risk assessment.

Kernel machine SNP-set testing under multiple candidate kernels

Joint testing for the cumulative effect of multiple single-nucleotide polymorphisms grouped on the basis of prior biological knowledge has become a popular and powerful strategy for the analysis of large-scale genetic association studies. The kernel machine (KM)-testing framework is a useful approach that has been proposed for testing associations between multiple genetic variants and many different types of complex traits by comparing pairwise similarity in phenotype between subjects to pairwise similarity in genotype, with similarity in genotype defined via a kernel function. An advantage of the KM framework is its flexibility: choosing different kernel functions allows for different assumptions concerning the underlying model and can allow for improved power. In practice, it is difficult to know which kernel to use a priori because this depends on the unknown underlying trait architecture and selecting the kernel which gives the lowest P-value can lead to inflated type I error. Therefore, we propose practical strategies for KM testing when multiple candidate kernels are present based on constructing composite kernels and based on efficient perturbation procedures. We demonstrate through simulations and real data applications that the procedures protect the type I error rate and can lead to substantially improved power over poor choices of kernels and only modest differences in power vs. using the best candidate kernel.

Genetic susceptibility to lung cancer--light at the end of the tunnel?

Lung cancer is one of the most common and deadliest cancers in the world. The major socio-environmental risk factor involved in the development of lung cancer is cigarette smoking. Additionally, there are multiple genetic factors, which may also play a role in lung cancer risk. Early work focused on the presence of relatively prevalent but low-penetrance alterations in candidate genes leading to increased risk of lung cancer. Development of new technologies such as genomic profiling and genome-wide association studies has been helpful in the detection of new genetic variants likely involved in lung cancer risk. In this review, we discuss the role of multiple genetic variants and review their putative role in the risk of lung cancer. Identifying genetic biomarkers and patterns of genetic risk may be useful in the earlier detection and treatment of lung cancer patients.