 Technical Advance
 Open Access
 Published:
Discovering weaker genetic associations guided by known associations
BMC Medical Genomics volume 13, Article number: 19 (2020)
Abstract
Background
The current understanding of the genetic basis of complex human diseases is that they are caused and affected by many common and rare genetic variants. A considerable number of the diseaseassociated variants have been identified by Genome Wide Association Studies, however, they can explain only a small proportion of heritability. One of the possible reasons for the missing heritability is that many undiscovered diseasecausing variants are weakly associated with the disease. This can pose serious challenges to many statistical methods, which seems to be only capable of identifying diseaseassociated variants with relatively stronger coefficients.
Results
In order to help identify weaker variants, we propose a novel statistical method, Constrained Sparse multilocus Linear Mixed Model (CSLMM) that aims to uncover genetic variants of weaker associations by incorporating known associations as a prior knowledge in the model. Moreover, CSLMM accounts for polygenic effects as well as corrects for complex relatednesses. Our simulation experiments show that CSLMM outperforms other competing existing methods in various settings when the combinations of MAFs and coefficients reflect different scenarios in complex human diseases.
Conclusions
We also apply our method to the GWAS data of alcoholism and Alzheimer’s disease and exploratively discover several SNPs. Many of these discoveries are supported through literature survey. Furthermore, our association results strengthen the belief in genetic links between alcoholism and Alzheimer’s disease.
Background
Genome Wide Association Studies (GWAS) have allowed people to address one of the most fundamental tasks in genetic research, which is to uncover associations between genetic variants and complex traits. Many efforts have been made which employ traditional statistical testing methods such as the Wald test to test the association of each individual SNP with a certain human disease, yet there are still a large amount of missing heritability to be discovered [1], which is due to the relatively low statistical power of these methods. In order to increase the power of the association mapping, many statistical approaches have been proposed.
For example, linear regression and the Lasso variants have been introduced to account for polygenic effects commonly seen in complex human diseases [2, 3]. Following the success of Lasso methods, the Adaptive Lasso with the oracle property under some regularity conditions [4], and the Precision Lasso that works with correlated and linearly dependent variables [3] were proposed.
However, a natural limitation of the Lassobased approaches is that they do not account for confounding effects raised by population structure and other complex relatedness in the GWAS data. In order to correct such effects, linear mixed models (LMMs) have been developed and received much attention in the recent years [5, 6]. Recently, Segural et al introduced a multilocus LMM that utilizes stepwise selection to model polygenetic effects [7]. Further Liu et al extended the multilocus LMM by dividing the model into fixed effect model and random effect model and use them iteratively [8]. On an alternative approach, recent studies also proposed a multilocus extension to the standard LMM to account for polygenic effects with the introduction of priors on coefficients [9, 10].
Despite the success of the aforementioned methods achieved, these methods are not effective in identifying genetic variants with weaker coefficients. Considering the current notion that many complex human diseases are likely to be caused and affected by many–rather than a few–genetic variants with small coefficients on a certain disease [11] and yet only a limited number of significant diseaseassociated variants have been identified from GWAS, we conjecture that the variants with small coefficients are difficult to identify given the presence of the variants with much larger coefficients, and that they will become easier to detect when conditioning on frequently reported SNPs which usually have larger coefficients. Following this belief, we propose a novel statistical method, Constrained Sparse Multilocus Linear Mixed Model (CSLMM), [12, 13] to uncover novel genetic variants of smaller coefficients by: 1) incorporating those frequently reported or known variants as a prior knowledge to the model, 2) accounting for polygenic association with a multivariate sparse regularized regression, and 3) correcting for population structure and complex relatedness (including family structure and other cypticx relatedness).
The performance of the CSLMM model is evaluated using extensive simulation experiments. We also apply our CSLMM model to an alcoholism and an Alzheimer’s Disease GWAS data, with the prior knowledge of the reported SNPs associated with each disease. We identify a set of SNPs having weak associations with each disease. Most of our findings are consistent with previously published results.
Methods
We formally introduce our model named Constrained Sparse Multilocus Linear Mixed Model (CSLMM) that aims to uncover genetic variants with weaker associations of a disease by incorporating variants of known associations as a prior knowledge.
Model
Given frequently reported or known variants (will be called known variants later for simplicity) with relatively larger coefficients, our model CSLMM aims to uncover novel variants of smaller coefficients. In order to achieve this, let X denote genotype data, Z denote population identification, y denote phenotype data (we first assume quantitative traits here, and discuss the casecontrol data or binary traits later), and let \(\mathcal {K}\) denote the set of the variants that are known or frequently reported. The “coefficient” is mathematically defined as the coefficient of linear regression [14]. With these settings, we have our CSLMM model formally presented as:
where β is the fixed genetic effects; u denotes the random population effects; ε is natural noise. We also introduce a constraint term β_{1}≤c with the belief that only a subset of the SNPs are associated with the phenotype, where c is a constant.
Algorithm
We proceed to introduce a threephase algorithm to estimate the parameter β, σ_{u}, and σ_{ε} in the CSLMM model.
Step I. Fitting known variants of larger coefficients: We first fit a linear regression model to determine the coefficients (magnitude of β_{i}) for the known SNPs, by solving the following equation:
$$\begin{array}{*{20}l} \hat{\beta_{i}} = \arg\min_{\beta_{i}} \mathbf{y}\sum_{i}\mathbf{X}^{i}\beta_{i}_{2}^{2}, \quad \forall i \in \mathcal{K} \end{array} $$(1)Step II. Correcting for population stratification and complex relatedness: Then, we consider to estimate σ_{u} and σ_{ε} for population stratification. Since y=Xβ+Zu+ε (u∼N(0,σ_{u}) and ε∼N(0,σ_{ε})) is equivalent to \(\mathbf {y} \sim N(\mathbf {X}\beta, \mathbf {Z}\mathbf {Z}^{T}\sigma _{u}^{2}+I\sigma _{\epsilon }^{2})\), we can estimate the variance term with a maximum likelihood estimation of Gaussian distribution by maximizing the following:
$$\begin{array}{*{20}l} l(\sigma_{u}, \sigma_{\epsilon}  \mathbf{y}', G) \propto N(\mathbf{y}'  \bar{\mathbf{y}'}  0, \sigma_{u}^{2}\mathbf{Z}\mathbf{Z}^{T} + \sigma_{\epsilon}^{2}\mathbf{I}) \end{array} $$(2)where \(\bar {\mathbf {y}'}\) is the empirical mean of y^{′} that is calculated by
$$\begin{array}{*{20}l} \mathbf{y}' = \mathbf{y}  \sum_{i} \mathbf{X}^{i}\hat{\beta_{i}} \end{array} $$(3)and ZZ^{T} is the genomic relationship matrix that is estimated as ZZ^{T}=(X^{j})(X^{j})^{T}, following the convention [15].
We then solve Eq. 2 for σ_{u} and σ_{ε}, where we can adopt the trick of introducing \(\delta = \frac {\sigma _{\epsilon }^{2}}{\sigma _{u}^{2}}\) to replace \(\sigma _{u}^{2}\) for more efficient optimization [16].
Finally, we can correct the population stratification by rotating the original data:
$$\begin{array}{*{20}l} \tilde{\mathbf{X}^{j}} &= (\textnormal{diag}(\mathbf{\Gamma})+\delta \mathbf{I})^{\frac{1}{2}}\mathbf{V}^{T}\mathbf{X}^{j} \\ \tilde{\mathbf{y}'} &= (\textnormal{diag}(\mathbf{\Gamma})+\delta \mathbf{I})^{\frac{1}{2}}\mathbf{V}^{T}\mathbf{y}' \end{array} $$where ZZ^{T}=UΓV^{T} is the singular value decomposition.
Step III. Fitting variants with smaller coefficients: Finally, we try to use the rest SNPs to explain the residual phenotypes, with solving the following:
$$\begin{array}{*{20}l} \hat{\beta_{j}} =& \arg\min_{\beta_{j}} \tilde{\mathbf{y}'}\sum_{j}\tilde{\mathbf{X}^{j}}\beta_{j}_{2}^{2} \\ & \textnormal{subject to}\quad \beta_{j} < \min\beta_{i}, \quad \forall j \quad \forall i \end{array} $$To solve this problem efficiently, we relax this constrain to a Lasso constrain as follows:
$$\begin{array}{*{20}l} \hat{\beta_{j}} =& \arg\min_{\beta_{j}} \tilde{\mathbf{y}'}\sum_{j}\tilde{\mathbf{X}^{j}}\beta_{j}_{2}^{2} + \sum_{j}\lambda\beta_{j}_{1} \end{array} $$(4)This new Lasso problem is solved via proximal gradient descent [17].
Stability Selection In Step III, to achieve a stable variable selection, we follow the regime of stability selection [18]: we run the algorithm 100 times, each time with half of the data points sampled without replacement from the original data. The final selected variables are the ones that are chosen more than 75% of chances over 100 runs.
Implementation
The implementation of CSLMM is available as a python software. Without installation, one can run the software with a single command line. It takes the Plink binary data as input. An extra file containing the known association variants is recommended. If this extra file is not available, CSLMM will first employ standard testing methods such as Wald test to select variants with the strongest signals. In order to identify a specific number (denoted as K) of SNPs associated with the disease, users can inquire the model with the number K or with a specific weight of the regularization term (λ in Eq. 4). If neither the number of SNPs nor the regularization weight is specified, the software will estimate the parameters using cross validation. The detailed instruction on how to use the software can be found in the Additional file 1. The implementation is available as a standalone software^{Footnote 1}. The computational complexity and scalability scales linearly with the number of samples and SNPs.
Results
Simulations
In order to evaluate the performance of CSLMM, we compare it with several existing association methods regarding their ability to uncover weaker associations. In particular, we compare CSLMM to the following methods:
Standard Wald test with the standard FDR control using the Benjamini–Hochberg (BH) procedure [19]: the most popular test used in GWA studies;
L1regularized linear regression (i.e. the Lasso);
Adaptive Lasso: an extension of Lasso that weighs the regularization term [4] (enabled by the method introduced in [20] for highdimensional data);
Precision Lasso: a novel improvement of Lasso that is more stable and consistent than Lasso [3];
Linear mixed model: the most popular method of population stratification;
Sparse linear mixed model (sparse LMM): a combination of sparse variable selection and population stratification [9, 21].
Multilocus linear mixed model (MLMM): an improvement of linear mixed model with stepwise selection to enable polygenetic modelling [7].
Fixed and random model Circulating Probability Unification (FarmCPU): a novel extension of MLMM that iteratively uses fixed effect model and random effect model [8]
Data generation
We generate the simulation data comprehensively to reflect real world scenarios of genetic data with population structure under different minor allele frequencies (MAFs) and coefficients. We use the SimuPop [22] software to simulate the real world genomic data with population structure. We simulate p SNPs for n individuals, denoted as X, and let X^{j} denote the j^{th} SNP. These individuals are from g populations and each population has f subpopulation.
In our simulation experiments, the SNPs come from two sets with two different MAFs: 20% of these SNPs are from one set (denoted as Set v) which has an MAF as m_{v} while the rest of the 80% SNPs are from the other set (denoted as Set u) which has a MAF as m_{u}. We assume there are k SNPs associated with the phenotype, of which, 20% are from set v and the rest are from set u.
In addition, the known SNPs in our simulation have higher MAFs and larger coefficients than the SNPs to be discovered. More specifically, for a SNP j, if j∈k and j∈v, it simulates the SNP that is already known to be associated with the trait and it has coefficient β_{j}=e_{v}c_{j}. On the other hand, if j∈k and j∈u, SNP j simulates the undiscovered associated SNP that has coefficient β_{j}=e_{u}c_{j}. If j∉k, SNP j simulates a SNP that is not associated with the phenotype and has the coefficient β_{j}=0c_{j}=0. c_{j} is the base coefficient, sampled from a uniform distribution U(0,1). This simulation process is showed in Fig. 1.
We generate the associated phenotype y as y=Xβ+ε, where ε∼N(0,1) is the natural noise. We further transform y into a binary phenotype with a Binomial sampling procedure with the probability of success achieved through feeding y into the inverse logit function.
Following [1], we conduct experiments with a variety of the settings with different combinations of MAFs (m_{u}=0.005,0.01), coefficients (e_{u}=5,10,25) of the SNPs to be discovered, and heritability (0.1,0.3,0.5,0.7) of the phenotype. For the known SNPs, we keep m_{v}=0.1 and e_{v}=50. We choose n=500, p=500000, and k=10 for the following experiments. For each configuration of the data, we repeat the experiments 10 times with different random seeds, and the reported result is based on the union of the results from all runs.
Evaluation
To conduct a fair comparison, we evaluate these models only regarding their ability to uncover the associated SNPs that are not already known to CSLMM, as CSLMM takes the known SNPs as a prior knowledge. For each method, we follow the convention to select the parameter λ (the weight of regularizer), which leads to the desired number of the selected variables (denoted as K) [3, 23]. This helps to avoid overly complex models, which tend to be selected by automatic measures such as cross validation, the Akaike information criterion (AIC), and the Bayesian information criterion (BIC) [24]. Moreover, it is known that the performance of parameter estimation and prediction are not directly coupled, e.g., as mentioned in [25] and the hyperparameter selected through crossvalidation tend to report more false positives [3]. In our experiments, we select exactly K=k variables.
Results
Figure 2 shows the precisionrecall curve of CSLMM compared to the Wald test, Lasso, Adaptive Lasso, Precision Lasso, LMM, sparse LMM, MLMM, and FarmCPU. The figure shows 24 experiments with three choices of coefficients (e_{u}) across two choices of MAFs m_{u} of the SNPs to be discovered, and four choices of heritability. In particular, plots in Figure 2 represent MAFs and coefficients correspond to heritability 0.1 (a), 0.3 (b), 0.5(c), and 0.7(d).
Figure 2a represents the most challenging case since the heratibility is as small as 0.1. All the methods do not behave well in this setting, and MLMM seems to have tiny advantages over other methods. Figure 2b and c illustrate the more realistic cases with heritabilities set as 0.3 and 0.5. Within this setup, we can see CSLMM has clear advantages over other methods. Sparse LMM and vanilla LMM are also behaving well, but still inferior to CSLMM. Figure 2d represents a simple scenario where the heritability is 0.7. In this setting, simpler univeriate testing methods, such as Wald and LMM, can also perform well, and CSLMM behave roughly slightly shy of these univariate testing methods. In general, CSLMM behave better than the competing methods in most settings of the experiments.
Other experiments
Other than the main experiment shown in Fig. 2, we have tested our methods in a larger range of choices of coefficients and MAF, tested the methods when we have different choices of k, and tested the methods under a larger number of samples. We also reported other other evaluation criteria including true positives, false positives and area under ROC (auROC) under a broader setting of the experiment. There more thorough tests are included in Additional file 1: Section 4–7.
Taken together, these results show that CSLMM outperforms other competing existing approaches in most cases, in particular, in the settings when the heratibility is at an intermediate level. Notably, these are also the settings that resemble real life scenarios for complex human diseases, and thus demonstrating the necessity and promising usages of CSLMM in the real life.
Application to real data
Alcoholism study
We apply our method CSLMM to the casecontrol GWAS data collected from subjects with and without alcoholism by The Center for Education and Drug Abuse Research (CEDAR) at the University of Pittsburgh. The data set consists of 383 individuals that include 305 subjects reported to be addicted to the consumption of alcohol through their lifetime. The data consists of 234 male subjects and 149 female subjects. The ages of these subjects range from 21 to 31. There are 519,138 genotyped SNPs in the data. The missing values are imputed as the mode of corresponding SNPs. To take the full advantage of our method, we collect the SNPs associated with alcoholism that are reported in GWAS Catalog [26] with pvalues smaller than 1e8 as the known SNPs to build in the CSLMM model. The four SNPs we collect include: rs1789891, rs7590720, rs2835872, and rs4478858. With these known alcoholismassociated SNPs fed into CSLMM, we run the model to uncover additional SNPs that have weaker associations with alcoholism.
We inquire 20 SNPs from the model, and CSLMM returns 21 predicted SNPs when converges, including the 4 known SNPs we feed into the model as a prior knowledge, and thus the model discovers 17 alcoholismassociated SNPs. Table 1 lists the SNPs associated with alcoholism that are identified by CSLMM. Since it is challenging to verify the reliability of these findings experimentally, we instead conduct a literature survey to find out whether the genes where these SNPs reside are linked to alcoholism or related disorders. Even though this type of “verification” may not provide conclusive evidence about the association between the identified SNPs and the disease, it can provide clues about whether the findings are worth further investigation.
Encouragingly, all the SNPs we discovered are linked to alcoholism, through the gene these SNPs reside in, in previously published results (shown in Table 1). For example, the 5^{th}, the 6^{th}, and the 17^{th} SNPs are within the region of the gene ADH1C, which encodes class I alcohol dehydrogenase, gamma subunit, a member of the alcohol dehydrogenase family. ADH1C has been shown to be associated with alcoholism in different populations [33]. Also, there are seven different SNPs residing within the region of KCNJ6, which encodes a member of the G proteincoupled inwardlyrectifying potassium channel. KCNJ6 is also reported to be associated with alcoholism previously [30]. The 9^{th} SNP resides within the region of PECR. Interestingly, previous evidence shows that PECR is not only associated with alcoholism [28], but also plays some role in Alzheimer’s disease [29]. A previous study reported that the protein level of PECR is significantly altered in the cortical lipid rafts of the murine model of AD, compared to the control mice [29]. This result is consistent with a previous study suggesting associations between daily alcohol users and Alzheimer’s patients [36].
The 10^{th} SNP is within the region of ADH1B, which is also known to be related with alcoholism. The 13^{th} SNP and the 20^{th} SNP are within in the region of gene ADH6, which is also known as an alcohol dependence gene [34].
Alzheimer’s disease study
Encouraged by our results from the alcoholism association mapping, we take a step further to investigate whether there is a genetic link between alcoholism and AD. We apply our method to a lateonset AD dataset provided by Harvard Brain Tissue Resource Center and Merck Research Laboratories [37]. The genotype data was generated from 540 subjects, and consists of the measurements for about 500,000 SNPs. There are 82 male subjects and 87 female subjects. The gender of the rest patients are unidentified. There are 366 subjects diagnosed with AD. The average age of these subjects is 56. The missing values are imputed as the mode of the corresponding SNPs. We use the two SNPs, rs2075650 (gene APOE) and rs157580 (gene TOMM40) as a prior knowledge to build into CSLMM. These two SNPs are reported to be associated with AD with pvalue less than 1e20 in GWAS Catalog [26]. We inquire the model for 20 SNPs that are associated with AD, and 22 SNPs are reported. The results are shown in Table 2. The reason that we use different thresholds (1e20 for Alzheimer’s disease and 1e8 for Alcoholism) to choose SNPs are prior knowledge is mainly due to the fact that Alzheimer’s disease is studied much more extensively than alcoholism in GWAS catalog, and pvalues for SNPs that are reported to be associated with Alzheimer’s disease tend to be smaller than those for alcoholism. We verify our findings following the same logic presented in the previous section.
Among the 19 SNPs associated with AD in Table 2, we found that the 6^{th} SNP within gene ABCA9 is previously reported associated with AD [41], confirming again that our method CSLMM can identify biologically meaningful variants. Also noticeably, the 15^{th} SNP resides within gene ESRRG, which encodes estrogen related receptor γ. Interestingly, evidence suggests that ERR γ plays key an role in alcoholinduced oxidative stress [42, 43]. This result also potentially verifies the existence of the pleiotropic effects between alcoholism and AD.
Since this short list of SNPs shows a promising application of CSLMM, we also apply CSLMM to identify a longer list of 200 SNPs for further studies. The longer list is reported in Additional file 1 (Section S2 and S3).
We also apply the competing existing methods to these two data sets, none of these methods identify a list of SNPs that are consistent with published results to the extent that CSLMM achieves.
Discussion
We developed a novel method: Constrained Sparse multilocus Linear Mixed Model (CSLMM) that conditions on the associations that have already been discovered to identify diseaseassociated SNPs with weaker signals. Our CSLMM model accounts for polygenic effects as well as corrects for complex relatedness such as population structure, family structure and cryptic relatedness. Our simulation experiments show that CSLMM outperforms other competing existing methods in terms of uncovering the variants with weaker signals in various settings which reflect real life scenarios for common and rare diseases. Interestingly, in the case of ’rare variants with weak coefficients’, which is categorized as the most challenging case in [1, 44], CSLMM is superior to other competing methods. Our simulations also show that CSLMM can particularly outperforms other methods consistently in terms of controlling false positives.
Furthermore, we apply CSLMM to alcoholism and AD studies. For about top 20 SNPs associated with either alcoholism or AD that CSLMM identifies, many of the SNPs reside within genes that were previously implicated in the corresponding diseases. Interestingly, our results further verify the pleiotropic effects between alcoholism and AD. The results indicate that two alcoholismassociated SNPs, rs7590720 (previously known) and rs1344694 (newly discovered), reside in PECR. The protein level of PECR was shown to be abnormally altered in a murine model of AD compared to the control mice, suggesting the involvement of PECR in the disease mechanism of AD. Similarly, our results also show that a novel ADassociated SNP, rs12563692, resides in ESRRG which encodes estrogen related receptor γ. Notably, ERR γ plays key an role in alcoholinduced oxidative stress and liver injury.
One interesting aspect regarding CSLMM is about the threephase learning algorithm we develop for estimating the parameters of the model. Two alternative strategies of learning the parameters are: 1) directly solving it as a convex optimization problem with explicit constrains; and 2) solving it as a standard Lasso with relaxation on the regularization on known associations. We tested these two algorithms in simulations, and our threephase learning algorithm outperforms these two alternative strategies.
To tailor CSLMM for casecontrol data or binary traits, a simple extension can be made that replaces the linear regression cost function with logistic regression cost function. Interestingly, our results indicates that CSLMM works well with casecontrol data as it is (data not shown), without any extensions required. In fact, extending CSLMM to logistic regression (or any other generalized linear models with a nontrivial link function) will affect the results adversely. For a generalized linear model, we believe CSLMM will only function as desire when the link function is identity.
Conclusions
In summary, we have proposed and developed a novel software tool, CSLMM, for disease association mapping which takes into account genetic variants of known associations, polygenic effects, as well as population structure and complex relatedness. The results from our simulation experiments and real data analysis demonstrate that CSLMM can be served as an effective tool for association studies for complex human diseases.
Availability of data and materials
The programs CSLMM is available at https://github.com/HaohanWang/CSLMM. The datasets used and analysed during the current study are available from the corresponding author on reasonable request.
Abbreviations
 AD:

Alzheimer’s disease
 CSLMM:

Constrained sparse multilocus linear mixed model
 GWAS:

Genome wide association studies
 LMM:

Linear mixed model
 MAF:

Minor allele frequency
 SNP:

Single nucleotide polymorphism
References
 1
Manolio TA, Collins FS, Cox NJ, Goldstein DB, Hindorff LA, Hunter DJ, McCarthy MI, Ramos EM, Cardon LR, Chakravarti A, et al.Finding the missing heritability of complex diseases. Nature. 2009; 461(7265):747–53.
 2
Ogutu JO, SchulzStreeck T, Piepho HP. Genomic selection using regularized linear regression models: ridge regression, lasso, elastic net and their extensions. BMC Proc. 2012; 6(S2). https://doi.org/10.1186/175365616s2s10.
 3
Wang H, Lengerich BJ, Aragam B, Xing EP, Stegle O. Precision Lasso: accounting for correlations and linear dependencies in highdimensional genomic data. Bioinformatics. 2019; 35(7):1181–7.
 4
Zou H. The adaptive lasso and its oracle properties. J Am Stat Assoc. 2006; 101(476):1418–29.
 5
Kang HM, Zaitlen NA, Wade CM, Kirby A, Heckerman D, Daly MJ, Eskin E. Efficient control of population structure in model organism association mapping. Genetics. 2008; 178(3):1709–23.
 6
Zhang Z, Ersoz E, Lai CQ, Todhunter RJ, Tiwari HK, Gore MA, Bradbury PJ, Yu J, Arnett DK, Ordovas JM, et al.Mixed linear model approach adapted for genomewide association studies. Nat Genet. 2010; 42(4):355–60.
 7
Segura V, Vilhjálmsson BJ, Platt A, Korte A, Seren Ü, Long Q, Nordborg M. An efficient multilocus mixedmodel approach for genomewide association studies in structured populations. Nat Genet. 2012; 44(7):825.
 8
Liu X, Huang M, Fan B, Buckler ES, Zhang Z. Iterative usage of fixed and random effect models for powerful and efficient genomewide association studies. PLoS Genet. 2016; 12(2):1005767.
 9
Rakitsch B, Lippert C, Stegle O, Borgwardt K. A lasso multimarker mixed model for association mapping with population structure correction. Bioinformatics. 2012; 29(2):206–14.
 10
Wang H, Aragam B, Xing EP. Variable selection in heterogeneous datasets: A truncatedrank sparse linear mixed model with applications to genomewide association studies. IEEE; 2017. https://doi.org/10.1109/bibm.2017.8217687.
 11
Eichler EE, Flint J, Gibson G, Kong A, Leal SM, Moore JH, Nadeau JH. Missing heritability and strategies for finding the underlying causes of complex disease. Nat Rev Genet. 2010; 11(6):446–50.
 12
Valdar W, Solberg LC, Gauguier D, Burnett S, Klenerman P, Cookson WO, Taylor MS, Rawlins JNP, Mott R, Flint J. Genomewide genetic association of complex traits in heterogeneous stock mice. Nat Genet. 2006; 38(8):879–87.
 13
Yang J, Ferreira T, Morris AP, Medland SE, Madden PA, Heath AC, Martin NG, Montgomery GW, Weedon MN, Loos RJ, et al.Conditional and joint multiplesnp analysis of gwas summary statistics identifies additional variants influencing complex traits. Nat Genet. 2012; 44(4):369–75.
 14
Park JH, Gail MH, Weinberg CR, Carroll RJ, Chung CC, Wang Z, Chanock SJ, Fraumeni JF, Chatterjee N. Distribution of allele frequencies and effect sizes and their interrelationships for common genetic susceptibility variants. Proc Nat Acad Sci. 2011; 108(44):18026–31.
 15
Kang HM, Sul JH, Service SK, Zaitlen NA, Kong Sy, Freimer NB, Sabatti C, Eskin E, et al.Variance component model to account for sample structure in genomewide association studies. Nat Genet. 2010; 42(4):348–54.
 16
Lippert C, Listgarten J, Liu Y, Kadie CM, Davidson RI, Heckerman D. Fast linear mixed models for genomewide association studies. Nat Methods. 2011; 8(10):833–5.
 17
Parikh N, Boyd S, et al.Proximal algorithms. Found Trends Ⓡ Optim. 2014; 1(3):127–239.
 18
Meinshausen N, Bühlmann P. Stability selection. J R Stat Soc Ser B Stat Methodol. 2010; 72(4):417–73.
 19
Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser B Stat Methodol. 1995; 57(1):289–300. https://doi.org/10.1111/j.25176161.1995.tb02031.x.
 20
Huang J, Ma S, Zhang CH. Adaptive Lasso for sparse highdimensional regression models. Stat Sin. 2008; Oct 1:1603–18.
 21
Wang H, Yang J. Multiple confounders correction with regularized linear mixed effect models, with application in biological processes. In: 2016 IEEE International Conference on Bioinformatics and Biomedicine (BIBM). IEEE: 2016. https://doi.org/10.1109/bibm.2016.7822753.
 22
Peng B, Kimmel M. simupop: a forwardtime population genetics simulation environment. Bioinformatics. 2005; 21(18):3686–3687.
 23
Wu TT, Chen YF, Hastie T, Sobel E, Lange K. Genomewide association analysis by lasso penalized logistic regression. Bioinformatics. 2009; 25(6):714–21.
 24
Meinshausen N, Bühlmann P. Highdimensional graphs and variable selection with the lasso. Ann Stat. 2006; 34(3):1436–62. https://doi.org/10.1214/009053606000000281.
 25
de los Campos G, Sorensen D, Gianola D. Genomic heritability: what is it?PLoS Genet. 2015; 11(5):1005048.
 26
MacArthur J, Bowler E, Cerezo M, Gil L, Hall P, Hastings E, Junkins H, McMahon A, Milano A, Morales J, et al.The new nhgriebi catalog of published genomewide association studies (gwas catalog). Nucleic Acids Res. 2017; 45(D1):896–901.
 27
Frank J, Cichon S, Treutlein J, Ridinger M, Mattheisen M, Hoffmann P, Herms S, Wodarz N, Soyka M, Zill P, et al.Genomewide significant association between alcohol dependence and a variant in the adh gene cluster. Addict Biol. 2012; 17(1):171–80.
 28
Treutlein J, Cichon S, Ridinger M, Wodarz N, Soyka M, Zill P, Maier W, Moessner R, Gaebel W, Dahmen N, et al.Genomewide association study of alcohol dependence. Arch Gen Psychiatr. 2009; 66(7):773–84.
 29
Chadwick W, Brenneman R, Martin B, Maudsley S. Complex and multidimensional lipid raft alterations in a murine model of alzheimer’s disease. Int J Alzheimers Dis. 2010; 2010:1–56. https://doi.org/10.4061/2010/604792.
 30
Kang SJ, Rangaswamy M, Manz N, Wang JC, Wetherill L, Hinrichs T, Almasy L, Brooks A, Chorlian DB, Dick D, et al.Familybased genomewide association study of frontal theta oscillations identifies potassium channel gene kcnj6. Genes Brain Behav. 2012; 11(6):712–9.
 31
Cooper A, Grigoryan G, GuyDavid L, Tsoory MM, Chen A, Reuveny E. Trisomy of the g proteincoupled k+ channel gene, kcnj6, affects reward mechanisms, cognitive functions, and synaptic plasticity in mice. Proc Nat Acad Sci. 2012; 109(7):2642–7.
 32
Zuo L, Wang K, Zhang XY, Krystal JH, Li CSR, Zhang F, Zhang H, Luo X. Nkain1–serinc2 is a functional, replicable and genomewide significant risk gene region specific for alcohol dependence in subjects of european descent. Drug Alcohol Depend. 2013; 129(3):254–64.
 33
Peng Q, Gizer IR, Wilhelmsen K, Ehlers C. Associations between genomic variants in alcohol dehydrogenase (adh) genes and alcohol symptomatology in american indians and european americans: Distinctions and convergence. Alcohol Clin Exp Res. 2017; 41(10):1695–704. https://doi.org/10.1111/acer.13480.
 34
Park BL, Kim JW, Cheong HS, Kim LH, Lee BC, Seo CH, Kang TC, Nam YW, Kim GB, Shin HD, et al.Extended genetic effects of adh cluster genes on the risk of alcohol dependence: from gwas to replication. Hum Genet. 2013; 132(6):657–68.
 35
Sherry ST, Ward MH, Kholodov M, Baker J, Phan L, Smigielski EM, Sirotkin K. dbsnp: the ncbi database of genetic variation. Nucleic Acids Res. 2001; 29(1):308–11.
 36
Zhou S, Zhou R, Zhong T, Li R, Tan J, Zhou H. Association of smoking and alcohol drinking with dementia risk among elderly men in china. Curr Alzheimer Res. 2014; 11(9):899–907.
 37
Zhang B, Gaiteri C, Bodea LG, Wang Z, McElwee J, Podtelezhnikov AA, Zhang C, Xie T, Tran L, Dobrin R, et al.Integrated systems approach identifies genetic nodes and networks in lateonset alzheimer’s disease. Cell. 2013; 153(3):707–20.
 38
Li H, Wetten S, Li L, Jean PLS, Upmanyu R, Surh L, Hosford D, Barnes MR, Briley JD, Borrie M, et al.Candidate singlenucleotide polymorphisms from a genomewide association study of alzheimer disease. Arch Neurol. 2008; 65(1):45–53.
 39
Naj AC, Beecham GW, Martin ER, Gallins PJ, Powell EH, Konidari I, Whitehead PL, Cai G, Haroutunian V, Scott WK, et al.Dementia revealed: novel chromosome 6 locus for lateonset alzheimer disease provides genetic evidence for folatepathway abnormalities. PLoS Genet. 2010; 6(9):1001130.
 40
Guipponi M, Santoni FA, Setola V, Gehrig C, Rotharmel M, Cuenca M, Guillin O, Dikeos D, Georgantopoulos G, Papadimitriou G, et al.Exome sequencing in 53 sporadic cases of schizophrenia identifies 18 putative candidate genes. PloS One. 2014; 9(11):112745.
 41
Piehler AP, Özcürümez M, Kaminski WE. Asubclass ATPbinding cassette proteins in brain lipid homeostasis and neurodegeneration. Front Psychiatry. 2012; 3:17. https://doi.org/10.3389/fpsyt.2012.00017.
 42
Kim DK, Kim YH, Jang HH, Park J, Kim JR, Koh M, Jeong WI, Koo SH, Park TS, Yun CH, et al.Estrogenrelated receptor γ controls hepatic cb1 receptormediated cyp2e1 expression and oxidative liver injury by alcohol. Gut. 2013; 62(7):1044–54. https://doi.org/10.1136/gutjnl2012303347.
 43
Han YH, Kim DK, Na TY, Ka NL, Choi HS, Lee MO. Ror α switches transcriptional mode of err γ that results in transcriptional repression of cyp2e1 under ethanolexposure. Nucleic Acids Res. 2016; 44(3):1095–104.
 44
McCarthy MI, Abecasis GR, Cardon LR, Goldstein DB, Little J, Ioannidis JP, Hirschhorn JN. Genomewide association studies for complex traits: consensus, uncertainty and challenges. Nature Rev Genet. 2008; 9(5):356–69.
Acknowledgements
The authors would like to thank Steven Knopf from University of Pittsburgh for instructions in using the Alcoholism data. The authors would also like to thank Dr. Bryon Aragam from Carnegie Mellon University for early stage discussion in the designing the method and experiments.
About this supplement
This article has been published as part of BMC Medical Genomics Volume 13 Supplement 3, 2020: Proceedings of the Joint International GIW & ABACBS2019 Conference: medical genomics (part 2). The full contents of the supplement are available online at https://bmcmedgenomics.biomedcentral.com/articles/supplements/volume13supplement3.
Funding
This work is funded and supported by the Department of Defense under Contract No. FA872105C0003 with Carnegie Mellon University for the operation of the Software Engineering Institute, a federally funded research and development center. This work is also supported by the National Institutes of Health grants R01GM093156 and P30DA035778.
Author information
Affiliations
Contributions
HW proposed and the idea, conducted the experiment and wrote the manuscript. MMV prepared the Alcoholism data. EPX read and wrote the manuscript. WW designed the experiment, read and wrote the manuscript. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Not applicable
Consent for publication
Not applicable
Competing interests
The authors declare that they have no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Additional file 1
Supplementary of Discovering Weaker Genetic Associations Guided by Known Associations, with Application to Alcoholism and Alzheimer’s Disease Studies. The file has instructions of using the software and extra experimental results.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Wang, H., Vanyukov, M.M., Xing, E.P. et al. Discovering weaker genetic associations guided by known associations. BMC Med Genomics 13, 19 (2020). https://doi.org/10.1186/s1292002006674
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s1292002006674
Keywords
 Weak association
 Linear mixed model
 GWAS