Computational and Mathematical Methods in Medicine
Volume 2017 (2017), Article ID 6742763, 17 pages
https://doi.org/10.1155/2017/6742763

Research Article

Pathway-Based Kernel Boosting for the Analysis of Genome-Wide Association Studies

Stefanie Friedrichs,¹ Juliane Manitz,2,3 Patricia Burger,¹ Christopher I. Amos,⁴ Angela Risch,5,6,7 Jenny Chang-Claude,⁸ Heinz-Erich Wichmann,9,10,11 Thomas Kneib,² Heike Bickeböller,¹ and Benjamin Hofner12,13

¹Institute of Genetic Epidemiology, University Medical Centre, Georg-August University Göttingen, Göttingen, Germany
²Department of Statistics and Econometrics, Georg-August University Göttingen, Göttingen, Germany
³Department of Mathematics and Statistics, Boston University, Boston, MA, USA
⁴Department of Community and Family Medicine, Geisel School of Medicine, Dartmouth College, Lebanon, NH, USA
⁵Division of Molecular Biology, University of Salzburg, Salzburg, Austria
⁶Translational Lung Research Center Heidelberg (TLRC-H), Member of the German Center for Lung Research (DZL), Heidelberg, Germany
⁷Division of Epigenomics and Cancer Risk Factors, German Cancer Research Center (DKFZ), Heidelberg, Germany
⁸Division of Cancer Epidemiology, German Cancer Research Center (DKFZ), Heidelberg, Germany
⁹Institute of Medical Informatics, Biometry and Epidemiology, Chair of Epidemiology, Ludwig-Maximilians University, Munich, Germany
¹⁰Helmholtz Center Munich, Institute of Epidemiology II, Munich, Germany
¹¹Institute of Medical Statistics and Epidemiology, Technical University Munich, Munich, Germany
¹²Department of Medical Informatics, Biometry and Epidemiology, Friedrich-Alexander-Universität Erlangen-Nürnberg, Erlangen, Germany
¹³Section Biostatistics, Paul-Ehrlich-Institut, Langen, Germany

Correspondence should be addressed to Stefanie Friedrichs; sfriedr2@gwdg.de

Received 10 February 2017; Revised 15 April 2017; Accepted 10 May 2017; Published 13 July 2017

Academic Editor: Angelo Facchiano

Copyright © 2017 Stefanie Friedrichs et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The analysis of genome-wide association studies (GWAS) benefits from the investigation of biologically meaningful gene sets, such as gene-interaction networks (pathways). We propose an extension to a successful kernel-based pathway analysis approach by integrating kernel functions into a powerful algorithmic framework for variable selection, to enable investigation of multiple pathways simultaneously. We employ genetic similarity kernels from the logistic kernel machine test (LKMT) as base-learners in a boosting algorithm. A model to explain case-control status is created iteratively by selecting pathways that improve its prediction ability. We evaluated our method in simulation studies adopting 50 pathways for different sample sizes and genetic effect strengths. Additionally, we included an exemplary application of kernel boosting to a rheumatoid arthritis and a lung cancer dataset. Simulations indicate that kernel boosting outperforms the LKMT in certain genetic scenarios. Applications to GWAS data on rheumatoid arthritis and lung cancer resulted in sparse models which were based on pathways interpretable in a clinical sense. Kernel boosting is highly flexible in terms of considered variables and overcomes the problem of multiple testing. Additionally, it enables the prediction of clinical outcomes. Thus, kernel boosting constitutes a new, powerful tool in the analysis of GWAS data and towards the understanding of biological processes involved in disease susceptibility.

1. Introduction

Many human diseases are complex in nature. They are caused by an interplay of several, often moderate genetic effects and environmental factors (i.e., demographic, clinical, and other nongenetic data [1]). Their genetic architecture is often analyzed in genome-wide association studies (GWAS). Herein, genetic information is represented by the genotypes of a multitude of single-nucleotide polymorphisms (SNPs) located across the whole genome. Numerous SNPs associated with various diseases have already been discovered in GWAS analyses; however they cannot account for the full heritability of the corresponding disease [2]. Different methods to approach this problem of missing heritability have been proposed, including the joint analysis of several SNPs representing a particular part of the genetic information, such as a gene or gene set.

Gene-set analysis methods facilitate the detection of associations between an individual’s genetic information and a phenotype of interest, for example, disease status. The joint analysis of several genes often leads to increased power, as it reduces the overall number of conducted tests and assists in the detection of moderate associations [3]. Furthermore, the results are usually more meaningful, as they are based on functional units rather than on single SNPs. One form of gene-set analysis is the investigation of pathways, such as networks of interacting genes responsible for a specific cell function or regulation [4]. The proteins coded by genes within a pathway can enhance or reduce the expression of other genes, to which we refer as activation or inhibition. Thus, genes interact directly as well as indirectly in a series of interconnected steps within pathways. Different types of biological pathway exist, for example, involved in metabolism or signal transduction. Faults in function can occur and such malfunction of biological pathways may lead to disease onset and development.

Large sample sizes are required to detect weak genetic effects influencing disease risk. Thanks to technical advances and the formation of data-sharing consortia in particular, larger GWAS datasets have become available over recent years. However, genotyping and participant recruitment are still cost and work intensive. Especially in rare diseases, taking as an example the analysis of histological subtypes of a disease, it is very challenging to achieve sample sizes that result in adequate power in analyses [5]. Another challenge we face is to understand the biological meaning of detected associations. It is often difficult to interpret the results of GWAS analysis in the elucidation of the precise biological processes and corresponding functional units influencing disease susceptibility. Single-pathway analysis methods are often successful in the identification of genetic effects influencing disease susceptibility. However, they usually can not discriminate causal biological processes from isolated effects included in pathways due to gene overlap [6, 7]. Another limitation of many pathway analysis approaches is the lacking ability to predict the disease state, or other outcomes of interest, based on the identified genetic effects.

Kernel methods in statistics have already been demonstrated as dealing well with the challenges faced when analyzing GWAS data [8, 9]. They are capable of handling high-dimensional data, without requiring any direct specification of the functional relationship between genetic effects. Furthermore, kernel methods are computationally efficient and allow the straightforward incorporation of environmental covariates [9–11]. Kernels are used to calculate a quantitative value from genotype data, which may be interpreted as reflecting the genetic similarity between each pair of individuals. Different kernels have been proposed in the analysis of pathways [9, 12, 13]. While some kernels only evaluate SNP membership in genes, others can also adjust for differing gene numbers and sizes or even include gene interaction structures or other information (please refer to Materials and Methods and [13] for an overview). We focussed on the network-based kernel, as it allows us to include interaction structures and has been demonstrated as being superior in performance for interconnected effects [13].

We extend kernel-based analysis of GWAS data by integrating a network-based kernel function into a boosting framework, in order to identify genetic variation modulating disease susceptibility. Boosting emerged from the field of machine learning and was later transferred to statistical modelling. It implements an ensemble of many weak learners (so-called base-learners, simple models that are slightly improved over random guessing) to optimize the predictive accuracy of a model [14]. Since it is able to combine the power from several predictors with weak signals into a strong prediction set [15, 16], it may prove to be a powerful tool in the analysis of GWAS. Component-wise boosting enforces variable selection and includes additional effect regularization, which makes it especially useful for high-dimensional data [17]. Model-based boosting can be seen as an extension of classic boosting approaches (see, e.g., [18, 19]). Diverse base-learners, which represent special effect types, may be chosen and combined arbitrarily [20]. Thus, boosting allows the simultaneous inclusion of genetic information and demographic or other environmental data. This joint investigation of multiple variables allows taking into account correlations between different pathways and will likely facilitate discrimination of causal biological processes from effects included in pathways only due to gene overlap. The derived models can be assessed and interpreted directly. Our kernel boosting approach overcomes the problem of multiple testing thanks to its inherent variable selection property [21]. Thereby the overall gain in power in the analysis of GWAS supports the analysis of smaller samples and moderate-to-weak genetic effects. Of note, the main focus of boosting (as well as of other machine learning methods) is not on hypothesis testing but on the development of a multivariable prediction model.

We applied our approach to two GWAS datasets, one on lung cancer and one on rheumatoid arthritis. Lung cancer is one of the most common forms of cancer, especially in industrialized nations. It is responsible for the greatest proportion of deaths caused by cancer worldwide [22]. Although the exposure to tobacco is known to be the major risk factor for lung cancer susceptibility, a number of genetic influences have been revealed by many studies [23]. The actual number of known genetic influences, excepting some specific lung cancer syndromes, is still limited, and each only accounts for a minor increase in disease risk. Rheumatoid arthritis is the most frequently occurring inflammatory disease of the joints, predominantly affecting the hands and feet. It is one of the major causes of disability and is strongly influenced by genetic factors in the human leukocyte antigen (HLA) region located on chromosome 6 [24, 25]. The investigation into these two diseases with different genetic architectures provides the ideal platform to evaluate the performance of our novel method.

In Section 2, we introduce the model structure utilized and describe the construction of network-based kernel functions. We provide a short introduction to boosting and derive the novel boosting algorithm with kernel-based base-learners. Section 3 comprises a description of the simulation study used to evaluate the method’s performance and an overview of the application to rheumatoid arthritis and lung cancer GWAS datasets. The results of the simulation study and GWAS analyses are summarized in Section 4. Finally, we end the paper with a discussion and an outlook.

1.1. Software

We used the statistical software environment R [26] to perform all analyses unless stated otherwise. The methodological developments were implemented in the R packages kangar00 [27] and mboost [28]. An exemplary application of the kernel boosting method to a simulated data set is given in Supplementary Material 2, available online at https://doi.org/10.1155/2017/6742763.

2. Materials and Methods

We aim to model the disease status of an individual, based on environmental covariates and genetic information obtained from GWAS. The genetic information given by the genotypes of different SNPs is mapped via genes to pathways. For each pathway, we compute a kernel matrix transforming the genotype vectors of each two individuals into a numeric value, which may be interpreted as the genetic similarity of the two individuals. Based on these matrices, we fit a kernel-based boosting model to identify relevant pathways and to find a prediction model for disease status. In the following paragraphs, we define all the relevant parts to this approach.

2.1. Model Definition and Notation

We assume an additive logistic regression model for the conditional probability of being a case for individual :

with additive predictor

where is the case-control indicator ( control; case), is the dimensional environmental covariate vector, and denotes the genotype vector of the SNPs of the th individual. Note that the non- or semiparametrically modelled genetic effects usually only depend on a pathway specific subset of SNPs, . However, for the sake of notational convenience we dropped the pathway index .

The vector represents the regression coefficients (including an intercept ) related to the environmental covariates. They typically include information on age, sex, or other traits relevant to the disease investigated. The genotype variables are coded as number of minor alleles, resulting in for any SNP and individual . The nonparametric functions , describe how the risk of being affected by the disease depends on the observed genotypes. Here, we aggregate the genotype information according to SNP membership in different gene interaction pathways.

2.2. Network-Based Kernels

Liu et al. [10] introduced the kernel machine framework to the field of pathway analysis. Since genes in pathways can include complex interactions, nonparametric approaches are advisable. The logistic kernel machine test (LKMT) can model the effect of a pathway on a binary outcome nonparametrically, while including parametrically modelled covariates. In the resulting logistic regression model, the genetic influence is incorporated by a function from the reproducing kernel Hilbert space generated by a positive definite kernel function .

In a genetic application, this kernel function is evaluated for the genotypes of each two individuals and , whereby the kernel matrix element is obtained. This value can be understood as the genetic similarity between the two individuals. To embed this definition into the mathematically well-defined framework of a reproducing kernel Hilbert space, the kernel matrix has to fulfill some requirements: it has to be quadratic, symmetric, and positive semidefinite. A variety of kernel functions are available. In the pathway-based analysis of GWAS data, a network-based kernel can be used, which is able to incorporate the pathway topology [13].

Assume denotes the pathway specific genotype matrix consisting of the genotype vectors , which include only the SNPs relevant for pathway , for all individuals. Then, the network-based kernel is defined by

where is an matrix mapping all SNPs to the investigated genes (including an adjustment to account for differing sizes of genes) and represents the (modified) matrix network adjacency matrix of gene interactions. To ensure positive semidefiniteness of the kernel, the network adjacency matrix is processed in a number of preparatory steps: if a gene is not represented by any SNPs in the investigated GWAS dataset, it cannot be considered in the analysis. To prevent loss of information about interactions in the network, genes which have previously been connected via the omitted gene will be linked directly. The new link’s weight is determined in a multiplicative fashion, based on the weights of the two omitted links. For a graphical representation refer to Figure 1. The resulting matrix is further mirrored along its diagonal and transformed to obtain positive semidefiniteness. The applied transformation is given by

where denotes the identity matrix and is a weight based on the smallest eigenvalue of . For more details, see [13].

Figure 1: Graphical representation of rewiring step in data preparation. Nodes are representing genes in the pathway, while edges indicate interactions between the corresponding genes. Assume the gene depicted in grey is not represented by any genetic markers in the considered study and thus cannot be analyzed. To retain information about the (indirect) interaction of the two genes previously linked to the omitted gene, a new direct link is established between them. Its interaction type is determined by multiplication of the weights inherent to the two dropped links.

2.3. Model-Based Boosting

Model fitting in general aims to minimize the loss when relating observed responses to an estimated model characterized by the additive predictor as defined in (2). Thus, boosting minimizes the empirical risk

where denotes a suitable loss function. Here, we use the negative binomial log-likelihood as loss function, which results in additive logistic regression models in analogy to the LKMT. In general, the loss function characterizes the model and can be defined in terms of a suitable negative log-likelihood or other appropriate loss functions, for example, the quadratic error loss for Gaussian regression or the absolute error loss for quantile regression. For an overview on loss functions see Hofner et al. [20]. Boosting solves this optimization problem via functional gradient descent by moving in the direction of the loss function’s steepest descent along the additive effects of predictor (2). This can be seen in the following (simplified) algorithm:

(1)Initialize the additive predictor with , and all function estimates with . Note that includes all kernels and possibly additional effects for environmental covariates.(2)For do the following:

(a)Compute the negative gradient of the loss function evaluated at the estimates of the previous iteration:

(b)Estimate the negative gradient vector separately for each effect in the additive predictor (2) by base-learners , with by fitting simple regression models via (penalized) least squares. Thus, each base-learner regresses the negative gradient vector separately on each of the predictors.(c)Choose the best-fitting base-learner with the minimal residual sum of squares.(d)Compute the update for the additive predictor by adding the best-fitting base-learner with a step-length factor :

 The corresponding update of function estimate is given by

 while

 for all .

Note that each base-learner usually depends on only one environmental covariate or one pathway based on a suitable subset of the genotypes of . However, other dependencies are also possible. For details on the algorithm, see [20]. A graphical display of the main features of the kernel boosting algorithm is given in Figure 2.

Figure 2: Graphical representation of the main features of the kernel boosting algorithm.

2.4. Model Tuning

The major tuning parameter of the functional gradient descent boosting algorithm is the number of iterations . We usually choose via cross-validation methods (such as bootstrap, -fold cross-validation, or subsampling) in order to avoid overfitting: one fits the model on the selected subset of the data and chooses such that it minimizes the empirical risk on the data that were not used to estimate the model. Subsampling is recommended to avoid overly complex models [29]. The step-length is another tuning parameter. In general it is of minor importance as long as it is relatively small. It determines the trade-off between speed of convergence and variable selection ability and is typically set to [30].

The current estimate of the additive predictor usually depends on only a subset of the possible predictors: as we select the best-fitting base-learner in each step and choose such that it maximizes prediction accuracy (i.e., usually relatively small so that not all base-learners are selected), boosting selects base-learners and thus variables. In our approach, we exploit this behaviour to identify genetic associations. Note that a base-learner can be selected multiple times. Hence, its function estimate , is the weighted sum with weights of the individual estimates over all iterations in which the base-learner was selected (see (8)).

2.5. Boosting with Network-Based Kernel as Base-Learner

To incorporate genotype data, aggregated to represent a particular pathway, we utilize kernel-based base-learners. Using a kernel function , we transform the definition of the genotypic information of all pairs of individuals to , as mentioned before, and collect them in the kernel matrix . With this matrix, we can estimate

The function is used to map the influence of SNP profiles to the clinical outcome (see (2)). As we expect patients with similar SNP profiles to have similar outcomes, we aim to discourage large differences in for genetically similar individuals. According to the standard penalization approaches in the boosting context, we thus introduce an additional smoothness constraint on the coefficient vector based on the kernel distances:

Thus, we define a separate kernel base-learner for each pathway in the boosting framework. Using the negative gradient vector from the th boosting iteration, we can estimate the coefficient vector of each base-learner (see step 2b of the algorithm) via penalized least squares

where we dropped the function index for the sake of notational convenience. Note that kernel matrix plays the role of design matrix as well as the role of penalty matrix with penalty parameter , which governs the smoothness of the estimate. Usually, the penalty parameter is chosen such that all base-learners have equal degrees of freedom to allow an unbiased selection. A common choice is four degrees of freedom if only smooth effects are used or one degree of freedom if linear effects are to be included; see Hofner et al. [21] for details.

In some rare cases, the derived kernel matrix is numerically not positive semidefinite (i.e., minimal deviations might occur), even though this should theoretically always be the case. To ensure a numerically positive semidefinite matrix , we apply transformation (4) not only to but also on the resulting kernel matrix . The proposed approach is very fast and results in smaller absolute differences in the matrix elements than alternatives such as the procedure suggested by Higham [31] (results not shown).

For numerical reasons, we reformulate the estimation problem from (12) by multiplying the design matrix with the inverse of the square root of the penalty matrix [32]. Thus, we obtain the design matrix

while the penalty matrix simplifies to the identity matrix . Now, we can equivalently write

A similar approach based on radial basis functions, which, for example, uses correlation functions to measure distances, was introduced to the boosting framework by Hofner [33].

2.6. Model Prediction Using Kernels

Boosting specifically aims to optimize prediction accuracy. As in all regression models, we can use the estimated coefficients to predict the outcome for new observations. However, some extra work is required to set up the kernel, that is, the design matrix, with new genotype data . In this context, the kernel can be understood to compute the similarity between genotype information of individuals to be predicted and the observations used to fit the model, the training data itself. Thus,

The resulting kernel has the dimension , with being new and previously used observations. Note that kernel matrix must no longer be of full rank nor be positive semidefinite. Using , we can predict the effect of a pathway on the outcome as

where is obtained as the weighted sum with weights over the estimates from (14) for all iterations in which the th base-learner was selected (see (8)).

2.7. Incorporation of Environmental Covariates

To incorporate environmental variables into the boosting model, we can choose different base-learners suited to different types of effect. Linear effect base-learners are suited to a continuous covariate such as patient age, while categorical effect base-learners facilitate the incorporation of categorical environmental variables such as gender. For details on inclusion of environmental variables, refer to [20].

With the inclusion of environmental variables as base-learners, these are also subject to the selection process inherent to boosting and compete with the pathway-based genetic effects. However, one usually wishes to consider only the added effect of genetic pathways. To ascertain that the model is corrected for environmental variables, one may include them as mandatory effects. This can be done by fitting a standard logistic regression model for the effect of the environmental variables on the clinical outcome and using the estimates as a start model (offset) for the boosting algorithm (see [34, 35]). This approach is very similar to the LKMT procedure, which tests if the logistic regression model can be improved via addition of a nonparametric effect incorporating a particular pathway.

3. Simulations and Applications

3.1. Simulation Study

To evaluate the performance of kernel boosting, we conducted a simulation study based on simulated SNP data in combination with gene networks from existing biological pathways. Pathway information was extracted from the Kyoto Encyclopedia of Genes and Genomes (KEGG) [36]. For simulation purposes, we considered a sample of networks, randomly chosen from the total of pathways available in January 2015. Please refer to Figure 3 for a list of these pathways and refer to Table 1 for their network topology characteristics. The primary aim of this study was to determine whether kernel boosting can detect associated pathways and is able to distinguish them from noninfluential pathways. Thus, we investigated the method’s performance on data without genetic effects (null case) including individuals and in six effect scenarios, differing in effect strengths (relative risk of and per allele) and sample sizes ( with a 1 : 1 ratio of cases to controls). Datasets for all scenarios were simulated for 100 replications. Note that these scenarios are small compared to typically available sample sizes nowadays. The reason can be found in the computational demands of the method for an insightful number of replications. Accordingly, comparably strong effects of markers were chosen to match the sample sizes used in our simulations.

Table 1: Description of network properties for pathway topology of pathways used in simulations, compared to the properties of the two effect pathways hsa04020 and hsa04022. Nodes equal the number of included genes, links give the number of interactions, inhibition links the count of interactions of inhibiting type, the average degree of a node is the mean number of adjacent edges, density is the ratio between numbers of existing links and possible links, diameter denotes the distance to the farthest node in the graph, transitivity (also called cluster coefficient) calculates the probability of adjacent vertices of a vertex being connected, and signed transitivity considers the type of interaction in this calculation.

Min
Mean
Median
Max
hsa04020
hsa04022

Nodes
29.00
103.60
86.5
398.00
180.00
167.00

Links
1.00
197.81
87.5
1493.00
297.00
372.00

Inhibition links
0.00
27.08
10.50
148.00
7.00
67.00

Average degree
0.07
3.18
2.36
15.62
3.30
4.46

Density
0.00
0.03
0.03
0.16
0.02
0.03

Inhibition degree
0.00
0.52
0.24
2.62
0.08
0.80

Diameter
1.00
7.36
7.00
18.00
6.00
7.00

Transitivity
0.00
0.02
0.00
0.14
0.00
0.03

Signed transitivity

0.01
0.00
0.10
0.00
0.03

Figure 3: Relative frequency of datasets in which a pathway was selected for pathways in the noninformative simulation scenario.

For each simulated dataset, we fitted a boosting model with pathway kernels. In order to tune the model, that is, to derive the optimal number of boosting steps , we used 20-fold subsampling for each model on each of the datasets with a maximum number of 200 iterations. Using the network-based kernel function in both methods, we compared the results from our kernel boosting approach on multiple pathways to those obtained from the single-pathway LKMT [9–11]. Additional simulations with cross-validated models and a maximum number of up to 1000 iterations were conducted to gain more insight into the proposed algorithm and are presented in Supplementary Material 1, Section A.

All genotypes were simulated with the help of a reference dataset from the International HapMap Consortium [37]. The reference data include 1,184 individuals of European descent (CEU) and a total of 1,440,616 SNPs, of which 116,565 are located on chromosome one. For each gene included in at least one of the selected pathways, we defined a pseudogene to represent the gene within our simulations. Such a pseudogene was a randomly selected DNA segment on chromosome one of the reference data including five different SNPs. Between each two sampled regions, we ensured a distance of at least 100 kilo base pairs to prevent distortive LD correlations between them [38]. The location of pseudogenes was left unchanged for all simulations, resulting in a realistic correlation structure for all simulation scenarios. In each of the 100 simulation runs, new genotype data for a total of SNPs in pseudogenes were simulated using the HAPGEN2 software. This software generates new haplotype data by combining a given set of reference haplotypes with previously simulated data. The detailed procedure is described in [39].

In the null case, noninformative genetic data were simulated for 1000 individuals. In each replication, new genotypes without association signals were generated for SNPs. The disease status was assigned at random with binomial probability of being a case, completely independent of genotype information. In each of the six effect scenarios, genotype data for a previously chosen equal number of cases and controls were simulated such that two pathways affected disease status. Association signals were included in three genes per causal pathway. In each of the resulting six genes, two randomly selected SNPs were chosen to be influential on the binary clinical outcome. Within one simulation scenario, all associated SNPs had the same effect strength and for each SNP the minor allele was influential. All effects were simulated as additive. To simplify the evaluation, we decided not to include environmental variables in these settings.

We chose two typical pathways (KEGG ids hsa04020 and hsa04022) to include causal genes. In accordance with the findings in [13], the influential genes in the two causal pathways were chosen to be interconnected within the corresponding pathway. Here, we additionally sampled one effect gene in each pathway, with the probability of being selected set to its betweenness centrality. Betweenness centrality measures the amount of shortest connections between each two genes in the network passing through the gene. Different studies have indicated that genes in topologically relevant positions of a pathway are more likely to be involved in disease association [40]. Two neighbouring genes of the sampled gene were randomly chosen to complete the connected scenario. In hsa04020, the genes GNA11, TACR1, and BDKRB2 were simulated to include SNPs influencing disease susceptibility. For hsa04022, genetic effects were placed on the genes PRKG2, ATP2B2, and KCNU1. For each of these genes, two SNPs were simulated as being influential on disease status. Note that existing biological pathways can have genes in common. Thus, beside our two pathways chosen to include influential effects, six additional pathways contain association signals. Refer to Table 2 for an overview of influential genes included in simulation pathways.

Table 2: Counts of included influential genes within pathways used for simulation purposes. Pathways without simulated causal genes are not displayed.

KEGG id
Name of pathway
Effect genes included

hsa04020
Calcium signaling pathway
4

hsa04022
cGMP-PKG signaling pathway
5

hsa04024
cAMP signaling pathway
1

hsa04080
Neuroactive ligand-receptor interaction
2

hsa04270
Vascular smooth muscle contraction
2

hsa04540
Gap junction
2

hsa04610
Complement and coagulation cascades
1

hsa05200
Pathways in cancer
2

Application: GWAS for Rheumatoid Arthritis and Lung Cancer. We considered the German Lung Cancer study (GLC) with 488 cases and 478 controls, based on the data of participants taken from the following three individual studies: Lung Cancer in the Young (LUCY), a population-based multicentre study run by the Helmholtz Zentrum Munich, and the University Medical Centre of the Georg-August-University in Goettingen. This study includes data of lung cancer patients under the age of 51 and family members recruited in German hospitals [41, 42]. The Heidelberg lung cancer case-control study, conducted by the German Cancer Research Centre (DKFZ) and the Thoraxklinik in Heidelberg, Germany, recruited cases and controls in a hospital-based study [43]. Additional controls were provided by Cooperative Health Research in the Augsburg Region (KORA), a population-based genome-wide study carried out by the Helmholtz Zentrum Munich [44]. A subset of the study participants of these three studies was chosen to form the German Lung Cancer GWAS. These individuals were genotyped on a HumanHap 550K SNP chip.

The second GWAS is a rheumatoid arthritis study of the North American Rheumatoid Arthritis Consortium (NARAC). It includes 868 cases from New York hospitals, in which rheumatoid arthritis was diagnosed based on the criteria of the American College of Rheumatology. Additionally, 1,194 controls matching in self-reported ethnic background were collected. All individuals were genotyped with the HumanHap500v1 array [45, 46].

For the rheumatoid arthritis study, we utilized gender as environmental covariate. In the lung cancer study, age and smoking exposure, measured in pack years, were also considered. To determine the pack year, one multiplies the number of packs of cigarettes smoked per day by the number of years an individual has smoked.

All GWAS data were subjected to strict quality control. Only individuals with a genotype call rate of at least 95% were considered. SNPs with more than 10% missing values or with a minor allele frequency (MAF) below 0.1% were excluded from further analysis. Missing values in remaining markers were imputed with BEAGLE [47]. No SNPs beyond the original chip were imputed. The base pair positions of all SNPs were updated to NCBI build 38 using the Ensembl database [48], which was accessed using the R package biomaRt [49, 50]. Gene start and end positions were extracted from the same database, also using NCBI build 38. SNPs with no unique position were excluded. Refer to Table 3 for an overview of the study characteristics. Note that, during analysis, only SNPs mapped to genes within pathways were considered. The assignment of SNPs to genes was based on their base pair location and gene boundaries. SNPs closely located to each other are often in linkage disequilibrium (LD). For SNP annotation, we specified gene regions including LD-blocks extending beyond gene boundaries, as recommended in [51].

Table 3: Characteristics of analyzed GWAS datasets. Numbers of case and control individuals after quality control and SNP numbers for several analysis stages are displayed. Preprocessing of SNPs included quality control of genotype data, as well as updating genomic SNP positions according to the latest information (genomic build 38). The last column indicates the total number of all SNPs annotated to a pathway under investigation.

Study
Cases/controls
SNPs genotyped
SNPs after preprocessing
SNPs in analysis

Lung cancer
467/468
561,466
533,062
148,938

Rheumatoid arthritis
866/1189
545,080
491,695
137,839

The KEGG database groups pathways in disjoint subsets according to their biological functionality. In the analysis of the rheumatoid arthritis and lung cancer data, we used a subgroup of 73 pathways connected to human diseases (see Table 4). The information on this group of pathways was downloaded in April 2016. An offset model containing only the environmental covariates was fitted for each of the studies to serve as start model for the kernel boosting of pathways.

Table 4: KEGG pathways in the human diseases class as downloaded in April 2016. Pathways are sorted according to value, derived from LKMT application on the rheumatoid arthritis dataset, in ascending order. values for pathways significantly associated after Bonferroni correction are listed. Pathways selected by kernel boosting on the same dataset are marked in italics. Pathways containing one or several genes belonging to the HLA complex are marked with an asterisk behind the id number.

KEGG id
Name of pathway
value

hsa05133
Pertussis

hsa05150
Staphylococcus aureus infection

hsa04933
AGE-RAGE signaling pathway in diabetic complications

hsa05169
Epstein-Barr virus infection

hsa05144
Malaria

hsa05206
MicroRNAs in cancer

hsa05330
Allograft rejection

hsa05200
Pathways in cancer

hsa05166
HTLV-I infection

hsa05030
Cocaine addiction

hsa05323
Rheumatoid arthritis

hsa05310
Asthma

hsa05134
Legionellosis

hsa04940
Type I diabetes mellitus

hsa05031
Amphetamine addiction

hsa05145
Toxoplasmosis

hsa05203
Viral carcinogenesis

hsa05332
Graft-versus-host disease

hsa05020
Prion diseases

hsa05143
African trypanosomiasis

hsa05222
Small-cell lung cancer

hsa05205
Proteoglycans in cancer

hsa05322
Systemic lupus erythematosus

hsa05161
Hepatitis B

hsa05410
Hypertrophic cardiomyopathy (HCM)

hsa05010
Alzheimer’s disease

hsa05142
Chagas disease (American trypanosomiasis)

hsa05168
Herpes simplex infection

hsa05012
Parkinson’s disease

hsa04932
Nonalcoholic fatty liver disease (NAFLD)

hsa05321
Inflammatory bowel disease (IBD)

hsa04931
Insulin resistance

hsa05219
Bladder cancer

hsa05215
Prostate cancer

hsa05202
Transcriptional misregulation in cancer

hsa05220
Chronic myeloid leukemia

hsa05146
Amoebiasis

hsa05414
Dilated cardiomyopathy

hsa05231
Choline metabolism in cancer

hsa05032
Morphine addiction

hsa05162
Measles

hsa05214
Glioma

hsa05164
Influenza A

hsa05416
Viral myocarditis

hsa05132
Salmonella infection

hsa05014
Amyotrophic lateral sclerosis (ALS)

hsa04930
Type II diabetes mellitus
Not significant

hsa05218
Melanoma
Not significant

hsa05140
Leishmaniasis
Not significant

hsa05213
Endometrial cancer
Not significant

hsa05211
Renal cell carcinoma
Not significant

hsa05340
Primary immunodeficiency
Not significant

hsa05160
Hepatitis C
Not significant

hsa05212
Pancreatic cancer
Not significant

hsa05016
Huntington’s disease
Not significant

hsa05221
Acute myeloid leukemia
Not significant

hsa04950
Maturity onset diabetes of the young
Not significant

hsa05412
Arrhythmogenic right ventricular cardiomyopathy (ARVC)
Not significant

hsa05223
Non-small-cell lung cancer
Not significant

hsa05034
Alcoholism
Not significant

hsa05130
Pathogenic Escherichia coli infection
Not significant

hsa05120
Epithelial cell signaling in Helicobacter pylori infection
Not significant

hsa05131
Shigellosis
Not significant

hsa05204
Chemical carcinogenesis
Not significant

hsa05100
Bacterial invasion of epithelial cells
Not significant

hsa05216
Thyroid cancer
Not significant

hsa05152
Tuberculosis
Not significant

hsa05210
Colorectal cancer
Not significant

hsa05230
Central carbon metabolism in cancer
Not significant

hsa05217
Basal cell carcinoma
Not significant

hsa05320
Autoimmune thyroid disease
Not significant

hsa05033
Nicotine addiction
Not significant

hsa05110
Vibrio cholerae infection
Not significant

For each pathway analyzed, the network-based kernel function with 4 degrees of freedom served as base-learner. The optimal number of iterations was derived via 20-fold subsampling and the default step length of was used. For the purpose of comparison, each of the pathways considered in GWAS data analysis was also tested individually on the corresponding data using the LKMT. The same environmental variables that were used in the offset model for boosting were also considered for the LKMT. Prediction accuracy was measured by the misclassification rate and the area under the ROC curve (AUC) for both datasets. Of note, prediction accuracy is influenced by the applied model but also by the dataset at hand, that is, the amount of information contained in the data. Additionally, we provided the cross-validation results, that is, the (average) negative binomial likelihood on the data that was not used for model fitting (see Supplementary Material , Section , for these results).

4. Results

4.1. Simulation Results

We compared the number of pathways each approach identified as associated with disease risk and considered the respective overlap in the results. The noninformative genetic data simulation comprised genotype data for 50 pathways and 1,000 individuals. Figure 3 displays the percentage of runs in which a pathway was selected. We can observe that the application of kernel boosting to these data does not lead to a high selection frequency for any pathway. Selection of pathways appears to be distributed randomly across all networks, not suggesting any clearly recognizable association with disease status. Note that, in kernel boosting, we do not conduct tests to evaluate the pathways’ influence but select pathways based on their predictive performance. Thus, we cannot calculate a type I error to evaluate our method’s performance. However, we can quantify the empirical type I error. Within simulation runs on pathways, a total number of false selections occurred. Thus, a pathway was falsely selected in % of all possible cases. In out of the simulation runs, no single pathway was chosen by the algorithm. Hence, we conclude that kernel boosting can be trusted to reliably avoid false positive selections in noninformative data.

Figures 4 and 5 compare the results of effect simulations with a relative risk of 1.5 per allele for 1,000 cases and 1,000 controls to those for 250 cases and 250 controls. (a) in each figure contains barplots indicating selection frequencies of the pathways across all simulation runs when applying kernel boosting to the corresponding simulation scenario. (b) compares these results with the selection frequencies using the LKMT. Here, both the percentages of results with a value below (lighter grey bars) and those with values below the Bonferroni-corrected significance level of (darker grey bars) are indicated. Pathways containing influential genes are additionally highlighted in italics.

(a)

(b)

Figure 4: Relative frequency of datasets in which a pathway was selected using (a) kernel boosting (, ) and (b) LKMT (, ) for a sample size of 2000 individuals. Pathways including effect genes are labeled in bold; numbers in brackets denote the count of included influential genes within the pathway. All effects were simulated with a relative risk of 1.5 per allele.

(a)

(b)

Figure 5: Relative frequency of datasets in which a pathway was selected using (a) kernel boosting (, ) and (b) LKMT (, ) for a sample size of 500 individuals. Pathways including effect genes are labeled in bold; numbers in brackets denote the count of included influential genes within the pathway. All effects were simulated with a relative risk of 1.5 per allele.

The results of kernel boosting in the sample of 2,000 individuals (Figure 4(a)) display three pathways clearly identified as influential on the clinical outcome, as their selection frequency is close to %. These are the pathways originally chosen to include genetic effects, hsa04020 and hsa04022, and the pathway hsa04610. It seems that the latter pathway is able to depict some of the information of the influential gene more effectively than the causal pathway for which it was originally simulated. This can be explained, as hsa04610 has the highest transitivity (), also known as global clustering coefficient, of all simulation pathways and contains an effect gene. As the network kernel was designed to work especially well in detection of interconnected genetic effects, the causal gene is identified very well in the pathway when using this base-learner. Note that the same pathway did not stand out in the noninformative simulation scenario. Thus, we conclude that high transitivity facilitates the detection of causal effects when using the network-based kernel but does not lead to false positives (i.e., here, pathways which do not contain any effect gene). Several other pathways were also selected, but only with very low frequencies. In the same simulation scenario, the LKMT had very high power to detect the two pathways simulated to affect disease risk, however, also detected other pathways including any of the causal genes on the Bonferroni-adjusted significance level (Figure 4(b)). Three of the six other effect-containing pathways were selected in almost % of the replications and two of the remaining ones in more than % and one other pathway which contained an effect gene was hardly selected.

Overall, this indicates that kernel boosting can identify the pathways with the most explanatory power with respect to disease status and is less likely than LKMT to select pathways due to overlapping effect genes (see [6] for a discussion). The reason can be found in the multivariate nature of the kernel boosting approach, in which pathways are not tested separately for their influence, but a multivariate model is fitted to incorporate multiple influential predictors at the same time.

Figure 5(a) reveals that the selection frequencies of associated pathways drop noticeably when sample size decreases. The same three pathways as in the larger sample reached the highest selection frequencies but here only between 20% and 60%. Simultaneously, the number of selections across nonassociated pathways increased slightly compared to the larger sample. This indicates that a reduction in sample size leads to less clear identification of the main influential pathways by kernel boosting. In Figure 5(b), we notice a similar behaviour of the selection frequency in LKMT analysis. Here again, the power to identify pathways, previously well detected in the larger sample, drops clearly with the smaller dataset. Regarding the percentage of detected pathways on the Bonferroni-corrected significance level, the drop is even more pronounced in the LKMT than for kernel boosting. This indicates that kernel boosting is less strongly influenced by sample size and may have greater potential in the identification of causal effects in smaller datasets for which the LKMT is underpowered.

Figures 6 and 7 compare the results of kernel boosting and the LKMT for differing effect sizes in equally sized samples of individuals. The graphics are structured as Figures 4 and 5, with kernel boosting selection frequencies plotted in (a) and LKMT selection frequencies in (b). Figure 6 contains a simulation scenario with relative risk of per causal allele and Figure 7 the results for a relative risk of per allele. Again, pathways containing influential genes are additionally highlighted.

(a)

(b)

Figure 6: Relative frequency of datasets in which a pathway was selected using (a) kernel boosting (, ) and (b) LKMT (, ) for sample sizes of 1000 individuals. Effect strength was set to relative risks of 1.5 per allele. Pathways including effect genes are labeled in bold; numbers in brackets denote the count of included influential genes within the pathway.

(a)

(b)

Figure 7: Relative frequency of datasets in which a pathway was selected using (a) kernel boosting (, ) and (b) LKMT (, ) for sample sizes of 1000 individuals. Effect strength was set to relative risks of 1.1 per allele. Pathways including effect genes are labeled in bold; numbers in brackets denote the count of included influential genes within the pathway.

In the kernel boosting plot in Figure 6(a), the three pathways standing out in Figure 4 again reached very high selection frequencies. All three bars decreased slightly in size compared to the scenario with individuals but still illustrate selections in more than % of simulation runs. Selection frequencies of the other effect pathways increased compared to the scenarios in Figure 4. However, as selections across noninfluential pathways occurred more frequently here, they cannot clearly be identified as influential based on their selection frequencies alone. In the LKMT analysis of this sample, the power to detect causal effects noticeably drops compared to the individuals’ sample illustrated in Figure 4(b). Comparing Figures 6 with 7, we can see a drop in selection frequencies as well as in power to detect associated pathways. In Figure 6, the two chosen effect pathways were detected in almost % and around % of simulation runs for both methods. In Figure 7, we observe that kernel boosting reaches selection frequencies of around % and %, while the LKMT with Bonferroni correction only achieves selection frequencies slightly greater than % and %, respectively. In a similar fashion to the results of the scenarios compared in Figures 4 and 5, both methods have higher power to detect associations for stronger effects; however the drop in power is less pronounced for kernel boosting. We conclude that kernel boosting firstly has no inferior performance in terms of power compared to the LKMT. It may even prove more likely to identify influential pathways with smaller genetic effects as it overcomes the multiple testing problem. Secondly, we infer that, in contrast to single-pathway testing approaches, kernel boosting has the ability to discriminate crucial biological processes associated with disease risk from effects included in pathways only due to overlapping genes.

4.2. GWAS Analysis Results

Kernel boosting on the human disease pathways in the lung cancer dataset resulted in selection of only the prion diseases pathway (KEGG id hsa05020). No other pathway was selected. The misclassification error of the tuned boosting model for lung cancer (evaluated at the optimal cut point as defined by the minimal Youden index) was and the AUC was . The ROC curve and the cross-validation results are presented in the Supplementary Material , Section . The LKMT with network-based kernel did not detect any associated pathway on the Bonferroni-corrected significance level. The prion diseases pathway appears ranked 20 out of 73 pathways, when sorting pathways according to ascending Bonferroni-corrected values. Prions are misfolded proteins capable of changing the structure of other, properly folded proteins into their own incorrect prion structure. They have mostly been reported in connection with neurodegenerative diseases [52]. Nevertheless, a connection with different forms of cancer has also previously been suspected [53, 54]. A full table of results from the analysis of the lung cancer dataset can be found in Supplementary Material , Section .

As expected, analysis of the rheumatoid arthritis dataset discovered a variety of pathways (compare results in [13]). Kernel boosting constructed an explanatory model for disease status based on selected pathways (see pathways written in italics in Table 4). It is well known that genes belonging to the human leukocyte antigen (HLA) complex are highly correlated with rheumatoid arthritis [55]. The HLA family, located on the short arm of chromosome 6, is a highly polymorphic genetic system mainly responsible for the regulation of the immune system [56]. In the human disease class, pathways contain at least one of the HLA genes. These pathways are marked with an asterisk in Table 4. Between the pathways containing HLA genes and the pathways selected by kernel boosting, there is an overlap of pathways. This may be explained by the multivariate nature of the method, in which only the pathway most clearly representing a particular genetic effect will be selected, conditionally on previously selected effects. Testing the human disease pathways’ influence on disease status with the LKMT resulted in a large number of 46 significantly associated pathways out of pathways after Bonferroni correction (see pathways with values in Table 4). These included almost all HLA pathways ( out of ). The more specific identification of influential pathways by kernel boosting provides a more complete basis to the understanding of the crucial biological processes involved in disease susceptibility. The misclassification error of the tuned boosting model for rheumatoid arthritis (evaluated at the optimal cut point as defined by the minimal Youden index) was and the AUC was . The ROC curve and the cross-validation results are presented in Supplementary Material , Section .

5. Discussion

We extend a successful method for single-pathway tests to a multivariate selection approach for simultaneous analysis of several pathways. The resulting kernel boosting method benefits from the advantages of a kernel-based analysis, while at the same time overcomes some of the limitations inherent to testing procedures.

Moreover, our multivariable approach to GWAS data analysis does not provide values, which only provide limited information on the relevance of a genetic effect. A more meaningful result would be an effect measure for the investigated trait or better still the ability to predict an outcome. Kernel boosting facilitates prediction, based on the selected influential variables, as was elucidated in the application where the overall prediction accuracy of each of the models was reported. Thus, it is also possible to interpret the influence of a specific genetic alteration by comparing the change in the predicted outcomes. A high degree of prediction accuracy for the model is ensured through the convenient evaluation of its performance on subsamples of the investigated dataset. This procedure usually results in good prediction accuracy and a sparse model.

Owing to the built-in shrinkage, our boosting approach is capable of dealing with correlated effects. Hence, correlated pathways, which partly include the same genes, can be handled within this framework. Thanks to the multivariable nature of the approach, only the best-fitting pathways, evaluated in terms of prediction accuracy, will be chosen to enter the model. Thus, only the pathway most clearly representing a particular genetic effect will be selected, depending on those pathways selected previously. Our observations support the statement by de Leeuw et al. [57] that competitive gene-set analysis methods (multivariate approach, pathways in competition), in contrast to self-contained approaches (univariate approach, one pathway at a time), can potentially differentiate widely spread heritability of polygenetic outcomes from causal biological processes. This property can be very helpful in the identification and understanding of specific biological functions involved in disease susceptibility.

We consider pathways as analysis units; however various other options exist. Single SNPs in transcribed or untranscribed regions, and SNP sets aggregated to represent a specific genomic region, environmental variables, or other variables, may be investigated and even combined arbitrarily within one model. For example, the application of our method to the genes comprising a pathway may help to identify key influential genes within the network (for gene boosting, see also the work of Ma et al. [58]; for good overviews of feature selection methods and machine learning tools in bioinformatics refer to [59, 60]). Known influential factors may be embedded in an initial model prior to the selection procedure to adjust for environmental or genetic effects. Furthermore, the considered effects can be incorporated into the model via a multitude of possible base-learners.

The choice of a base-learner can influence effect selections. We observed this behaviour during the simulations, in which the highly connected pathway containing only one effect gene was identified owing to the network-based kernel’s high power on interconnected effects. Thus, the well-considered selection of base-learners to be utilized is advisable. We account for the high complexity of possible gene interactions in pathways via the use of a kernel function, which accounts for additive and interaction effects. Such a kernel function will likely lead to a higher degree of prediction accuracy than a simple linear kernel. The application of our method to GWAS datasets on rheumatoid arthritis and lung cancer returned biologically plausible results. Particularly with the rheumatoid arthritis dataset, the number of identified pathways could be reduced considerably compared to single-pathway tests. While the LKMT resulted in significantly associated pathways, kernel boosting narrowed the selection down to pathways. Genes within the HLA region are known to have a strong influence on rheumatoid arthritis. Their effects can reach far across pathways, such that the LKMT detects many pathways including HLA genes as significantly associated. Boosting seems to help to pinpoint down signals even among those pathways and reduces the number of identified pathways to a more reasonable level.

Our results indicate that kernel boosting outperforms single kernel machine tests, as exemplified by the LKMT, in certain genetic scenarios. It may help to discriminate causal biological processes from isolated effects included in pathways only due to gene overlap and facilitate discovering weak signals, especially in studies of limited size. This is of particular interest in the investigation of rare diseases and disease subtypes, in which established methods often fail to find any significantly associated pathways owing to a lack of power.

Datasets of the size investigated here can be analyzed with kernel boosting quite efficiently on current high-performance cluster computing (HPCC) systems. However, such analysis of very large datasets places a rather high demand even on the most powerful HPCC systems to date. Usually, our kernel base-learners are based on the pairwise similarities of all observations. This leads to similarity matrices as design matrices and hence to parameter vectors of size . Instead of using all pairwise similarities, it is possible to compute the similarities only to a representative subset of the observations, or so-called knots. These knots can be chosen as subset of the observations which covers the complete observation space (space-filling algorithm; see [33, 61, 62]). Consequently, we obtain reduced-rank design matrices of dimension , where is the number of knots, and a parameter vector of size . This reduces the computational burden for the construction of the kernel base-learners and effect estimation and makes kernel-based methods even feasible in situations with many observations. The exact number of observations that can be processed depends, among others, on the considered number of individuals, SNPs, base-learners chosen, and the available hardware.

Kernel boosting constitutes a new and potentially powerful tool in the analysis of GWAS data. It offers a highly flexible and extensible framework, suitable for a wide range of application scenarios. We account for the high complexity of possible gene interactions via the use of kernel functions, while reducing the complexity of the resulting model with the built-in shrinkage of the boosting approach. The resulting model enables us to predict traits and returns more meaningful results than a testing procedure. We conclude that kernel boosting is a suitable methodological addition for the analysis of GWAS, which supports the detection and interpretation of genetic risk factors influencing disease susceptibility.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work was supported by the German Research Foundation, Research Training Group 1644 “Scaling Problems in Statistics.” The rheumatoid arthritis data considered in this article were provided by the National Institutes of Health (Grant AR44422). The analyzed lung cancer study was made available through TRICL Grant no. U19CA148127. The authors would like to thank Andrew Entwistle for his critical review of the manuscript.

References

J. Craig, “Complex diseases: Research and applications,” Nature Education, vol. 1, article 184, no. 1, 2008.

T. A. Manolio, F. S. Collins, N. J. Cox et al., “Finding the missing heritability of complex diseases,” Nature, vol. 461, no. 7265, pp. 747–753, 2009.

G. Fehringer, G. Liu, L. Briollais et al., “Comparison of pathway analysis approaches using lung cancer GWAS data sets,” PLoS ONE, vol. 7, no. 2, Article ID e31816, 2012.

R. M. Cantor, K. Lange, and J. S. Sinsheimer, “Prioritizing gwas results: a review of statistical methods and recommendations for their application,” The American Journal of Human Genetics, vol. 86, no. 1, pp. 6–22, 2010.

E. P. Hong and J. W. Park, “Sample size and statistical power calculation in genetic association studies,” Genomics & Informatics, vol. 10, no. 2, pp. 117–122, 2012.

P. Khatri, M. Sirota, and A. Butte, “Ten years of pathway analysis: current approaches and outstanding challenges,” PLoS Computational Biology, vol. 8, no. 2, Article ID e1002375, 2012.

M. García-Campos, J. Espinal-Enríquez, and E. Hernández-Lemus, “Pathway analysis: state of the art,” Frontiers in Physiology, vol. 6, 2015.

W. Pan, “Network-based model weighting to detect multiple loci influencing complex diseases,” Human Genetics, vol. 124, no. 3, pp. 225–234, 2008.

M. C. Wu, P. Kraft, M. P. Epstein et al., “Powerful snp-set analysis for case-control genome-wide association studies,” The American Journal of Human Genetics, vol. 86, no. 6, pp. 929–942, 2010.

D. Liu, D. Ghosh, and X. Lin, “Estimation and testing for the effect of a genetic pathway on a disease outcome using logistic kernel machine regression via logistic mixed models,” BMC Bioinformatics, vol. 9, article no. 292, 2008.

S. Basu, W. Pan, and W. S. Oetting, “A dimension reduction approach for modeling multi-locus interaction in case-control studies,” Human Heredity, vol. 71, no. 4, pp. 234–245, 2011.

S. Freytag, H. Bickeböller, C. I. Amos, T. Kneib, and M. Schlather, “A novel kernel for correcting size bias in the logistic kernel machine test with an application to rheumatoid arthritis,” Human Heredity, vol. 74, no. 2, pp. 97–108, 2013.

S. Freytag, J. Manitz, M. Schlather et al., “A network-based kernel machine test for the identification of risk pathways in genome-wide association studies,” Human Heredity, vol. 76, no. 2, pp. 64–75, 2014.

A. Mayr, H. Binder, O. Gefeller, and M. Schmid, “The evolution of boosting algorithms: From machine learning to statistical modelling,” Methods of Information in Medicine, vol. 53, no. 6, pp. 419–427, 2014.

J. H. Friedman, “Greedy function approximation: a gradient boosting machine,” The Annals of Statistics, vol. 29, no. 5, pp. 1189–1232, 2001.

J. Friedman, T. Hastie, and R. Tibshirani, “Additive logistic regression: a statistical view of boosting,” The Annals of Statistics, vol. 28, no. 2, pp. 337–407, 2000.

P. Bühlmann and B. Yu, “Boosting with the L 2 Loss: Regression and Classification,” Journal of the American Statistical Association, vol. 98, no. 462, pp. 324–339, 2003.

T. Kneib, T. Hothorn, and G. Tutz, “Variable selection and model choice in geoadditive regression models,” Biometrics. Journal of the International Biometric Society, vol. 65, no. 2, pp. 626–634, 2009.

T. Hothorn, P. Bühlmann, T. Kneib, M. Schmid, and B. Hofner, “Model-based boosting 2.0,” Journal of Machine Learning Research (JMLR), vol. 11, pp. 2109–2113, 2010.

B. Hofner, A. Mayr, N. Robinzonov, and M. Schmid, “Model-based boosting in R: a hands-on tutorial using the R package mboost,” Computational Statistics, vol. 29, no. 1-2, pp. 3–35, 2014.

B. Hofner, T. Hothorn, T. Kneib, and M. Schmid, “A framework for unbiased model selection based on boosting,” Journal of Computational and Graphical Statistics, vol. 20, no. 4, pp. 956–971, 2011.

J. Ferlay, I. Soerjomataram, R. Dikshit et al., “Cancer incidence and mortality worldwide: sources, methods and major patterns in GLOBOCAN 2012,” International Journal of Cancer, 2014.

P. Brennan, P. Hainaut, and P. Boffetta, “Genetics of lung-cancer susceptibility,” The Lancet Oncology, vol. 12, no. 4, pp. 399–408, 2011.

G. S. Firestein, “Evolving concepts of rheumatoid arthritis,” Nature, vol. 423, no. 6937, pp. 356–361, 2003.

S. Raychaudhuri, “Recent advances in the genetics of rheumatoid arthritis,” Current Opinion in Rheumatology, vol. 22, no. 2, pp. 109–118, 2010.

R Core Team, R: A Language and Environment for Statistical Computing, R Foundation for Statistical Computing, Vienna, Austria, 2011, https://www.R-project.org/.

J. Manitz, S. Friedrichs, P. Burger et al., “kangar00: Kernel Approaches for Nonlinear Genetic Association Regression,” R package version 1.0, 2017, https://CRAN.R-project.org/package=kangar00.

T. Hothorn, P. Buehlmann, T. Kneib, M. Schmid, and B. Hofner, “mboost: Model-Based Boosting,” R package version 2.8-0, 2017, http://CRAN.R-project.org/package=mboost.

A. Mayr, B. Hofner, and M. Schmid, “The importance of knowing when to stop—a sequential stopping rule for component-wise gradient boosting,” Methods of Information in Medicine, vol. 51, no. 2, pp. 178–186, 2012.

M. Schmid and T. Hothorn, “Boosting additive models using component-wise P-splines,” Computational Statistics and Data Analysis, vol. 53, no. 2, pp. 298–311, 2008.

N. J. Higham, “Computing the nearest correlation matrix—a problem from finance,” IMA Journal of Numerical Analysis, vol. 22, no. 3, pp. 329–343, 2002.

E. E. Kammann and M. P. Wand, “Geoadditive models,” Journal of the Royal Statistical Society. Series C. Applied Statistics, vol. 52, no. 1, pp. 1–18, 2003.

B. Hofner, Boosting in Structured Additive Models, LMU München, 2011, http://nbn-resolving.de/urn:nbn:de:bvb:19-138053.

A. L. Boulesteix and T. Hothorn, “Testing the additional predictive value of high-dimensional molecular data,” BMC Bioinformatics, vol. 11, article 78, 2010.

R. De Bin, “Boosting in Cox regression: a comparison between the likelihood-based and the model-based approaches with focus on the R-packages CoxBoost and mboost,” Computational Statistics, vol. 31, no. 2, pp. 513–531, 2016.

M. Kanehisa, S. Goto, Y. Sato, M. Kawashima, M. Furumichi, and M. Tanabe, “Data, information, knowledge and principle: back to metabolism in KEGG,” Nucleic Acids Research, vol. 42, no. 1, pp. D199–D205, 2014.

International HapMap Consortium, “The international hapmap project,” Nature, vol. 426, pp. 789–796, 2003.

D. E. Reich, M. Cargili, S. Boik et al., “Linkage disequilibrium in the human genome,” Nature, vol. 411, no. 6834, pp. 199–204, 2001.

Z. Su, J. Marchini, and P. Donnelly, “HAPGEN2: simulation of multiple disease SNPs,” Bioinformatics, vol. 27, no. 16, Article ID btr341, pp. 2304-2305, 2011.

Y. Lee, H. Li, J. Li et al., “Network models of genome-wide association studies uncover the topological centrality of protein interactions in complex diseases,” Nature, vol. 20, no. 4, pp. 619–629, 2013.

W. Sauter, A. Rosenberger, L. Beckmann et al., “Matrix metalloproteinase 1 (MMP1) is associated with early-onset lung cancer,” Cancer Epidemiology Biomarkers and Prevention, vol. 17, no. 5, pp. 1127–1135, 2008.

A. Rosenberger, T. Illig, K. Korb et al., “Do genetic factors protect for early onset lung cancer? A case control study before the age of 50 years,” BMC Cancer, vol. 8, no. 1, article 60, 2008.

H. Dally, K. Gassner, B. Jäger et al., “Myeloperoxidase (MPO) genotype and lung cancer histologic types: The MPO-463 a allele is associated with reduced risk for small cell lung cancer in smokers,” International Journal of Cancer, vol. 102, no. 5, pp. 530–535, 2002.

H.-E. Wichmann, C. Gieger, and T. Illig, “KORA-gen—resource for population genetics, controls and a broad spectrum of disease phenotypes,” Gesundheitswesen, vol. 67, Supplement 1, pp. S26–S30, 2005.

C. I. Amos, W. Chen, M. F. Seldin et al., “Data for genetic analysis workshop 16 problem 1, association analysis of rheumatoid arthritis data,” BMC Proceedings, vol. 3, Supplement 7, 2009.

R. M. Plenge, M. Seielstad, L. Padyukov et al., “TRAF1-C5 as a risk locus for rheumatoid arthritis—a genomewide study,” The New England Journal of Medicine, vol. 357, no. 12, pp. 1199–1209, 2007.

S. R. Browning and B. L. Browning, “Rapid and accurate haplotype phasing and missing data inference for whole genome association studies by use of localized haplotype clustering,” The American Journal of Human Genetics, vol. 81, no. 5, pp. 1084–1097, 2007.

A. Yates, W. Akanni, M. R. Amode et al., “Ensembl 2016,” Nucleic Acids Research, vol. 44, no. D1, pp. D710–D716, 2016.

S. Durinck, P. T. Spellman, E. Birney, and W. Huber, “Mapping identifiers for the integration of genomic datasets with the R/Bioconductor package biomaRt,” Nature Protocols, vol. 4, no. 8, pp. 1184–1191, 2009.

S. Durinck, Y. Moreau, A. Kasprzyk et al., “BioMart and bioconductor: a powerful link between biological databases and microarray data analysis,” Bioinformatics, vol. 21, no. 16, pp. 3439-3440, 2005.

D. Malzahn, S. Friedrichs, and H. Bickeböller, “Comparing strategies for combined testing of rare and common variants in whole sequence and genome-wide genotype data,” BMC Proceedings, vol. 10, Supplement 7, pp. 269–273, 2016.

P. L. A. Leighton and W. Ted Allison, “Protein misfolding in prion and prion-like diseases: Reconsidering a required role for protein loss-of-function,” Journal of Alzheimer's Disease, vol. 54, no. 1, pp. 3–29, 2016.

H. Antony, A. P. Wiegmans, M. Q. Wei, Y. O. Chernoff, K. K. Khanna, and A. L. Munn, “Potential roles for prions and protein-only inheritance in cancer.,” Cancer metastasis reviews, vol. 31, no. 1-2, pp. 1–19, 2012.

J. L. Silva, L. P. Rangel, D. C. F. Costa, Y. Cordeiro, and C. V. De Moura Gallo, “Expanding the prion concept to cancer biology: dominant-negative effect of aggregates of mutant p53 tumour suppressor,” Bioscience Reports, vol. 33, no. 4, pp. 593–603, 2013.

C. M. Weyand and J. J. Goronzy, “Association of MHC and rheumatoid arthritis. HLA polymorphisms in phenotypic variants of rheumatoid arthritis,” Arthritis Research, vol. 2, no. 3, pp. 212–216, 2000.

S. Y. Choo, “The HLA system: genetics, immunology, clinical testing, and clinical implications,” Yonsei Medical Journal, vol. 48, no. 1, pp. 11–23, 2007.

C. A. de Leeuw, B. M. Neale, T. Heskes, and D. Posthuma, “The statistical properties of gene-set analysis,” Nature Reviews Genetics, vol. 17, no. 6, pp. 353–364, 2016.

S. Ma, Y. Huang, J. Huang, and K. Fang, “Gene network-based cancer prognosis analysis with sparse boosting,” Genetics Research, vol. 94, no. 4, pp. 205–221, 2012.

Y. Saeys, I. Inza, and P. Larrañaga, “A review of feature selection techniques in bioinformatics,” Bioinformatics, vol. 23, no. 19, pp. 2507–2517, 2007.

P. Larrañaga, B. Calvo, R. Santana et al., “Machine learning in bioinformatics,” Briefings in Bioinformatics, vol. 7, no. 1, pp. 86–112, 2006.

M. E. Johnson, L. M. Moore, and D. Ylvisaker, “Minimax and maximin distance designs,” Journal of Statistical Planning and Inference, vol. 26, no. 2, pp. 131–148, 1990.

D. Nychka and N. Saltzman, “Design of air-quality monitoring networks,” in Case Studies in Environmental Statistics, D. Nychka, W. W. Piegorsch, and L. H. Cox, Eds., vol. 132, pp. 51–76, Springer US, New York, NY, USA, 1998.