Logistic Bayesian LASSO for detecting association combining family and case-control data

Functional coding haplotypes and machine-learning feature elimination identifies predictors of Methotrexate Response in Rheumatoid Arthritis patients

BACKGROUND

Major challenges in large scale genetic association studies include not only the identification of causative single nucleotide polymorphisms (SNPs), but also accounting for SNP-SNP interactions. This study thus proposes a novel feature engineering approach integrating potentially functional coding haplotypes (pfcHap) with machine-learning (ML) feature selection to identify biologically meaningful, possibly causative genetic factors, that take into consideration potential SNP-SNP interactions within the pfcHap, to best predict for methotrexate (MTX) response in rheumatoid arthritis (RA) patients.

METHODS

Exome sequencing from 349 RA patients were analysed, of which they were split into training and unseen test set. Inferred pfcHaps were combined with 30 non-genetic features to undergo ML recursive feature elimination with cross-validation using the training set. Predictive capacity and robustness of the selected features were assessed using six popular machine learning models through a train set cross-validation and evaluated in an unseen test set.

FINDINGS

Significantly, 100 features (95 pfcHaps, 5 non-genetic factors) were identified to have good predictive performance (AUC: 0.776-0.828; Sensitivity: 0.656-0.813; Specificity: 0.684-0.868) across all six ML models in an unseen test dataset for the prediction of MTX response in RA patients.

INTERPRETATION

Majority of the predictive pfcHap SNPs were predicted to be potentially functional and some of the genes in which the pfcHap resides in were identified to be associated with previously reported MTX/RA pathways.

FUNDING

Singapore Ministry of Health's National Medical Research Council (NMRC) [NMRC/CBRG/0095/2015; CG12Aug17; CGAug16M012; NMRC/CG/017/2013]; National Cancer Center Research Fund and block funding Duke-NUS Medical School.; Singapore Ministry of Education Academic Research Fund Tier 2 grant MOE2019-T2-1-138.

Detecting rare haplotype association with two correlated phenotypes of binary and continuous types

Multiple correlated traits/phenotypes are often collected in genetic association studies and they may share a common genetic mechanism. Joint analysis of correlated phenotypes has well-known advantages over one-at-a-time analysis including gain in power and better understanding of genetic etiology. However, when the phenotypes are of discordant types such as binary and continuous, the joint modeling is more challenging. Another research area of current interest is discovery of rare genetic variants. Currently there is no method available for detecting association of rare (or common) haplotypes with multiple discordant phenotypes jointly. Our goal is to fill this gap specifically for two discordant phenotypes. We consider a rare haplotype association method for a binary phenotype, logistic Bayesian LASSO (univariate LBL) and its extension for two correlated binary phenotypes (bivariate LBL-2B). Under this framework, we propose a haplotype association test with binary and continuous phenotypes jointly (bivariate LBL-BC). Specifically, we use a latent variable to induce correlation between the two phenotypes. We carry out extensive simulations to investigate bivariate LBL-BC and compare it with univariate LBL and bivariate LBL-2B. In most settings, bivariate LBL-BC performs the best. In only two situations, bivariate LBL-BC has similar performance-when the two phenotypes are (1) weakly or not correlated and the target haplotype affects the binary phenotype only and (2) strongly positively correlated and the target haplotype affects both phenotypes in positive direction. Finally, we apply the method to a data set on lung cancer and nicotine dependence and detect several haplotypes including a rare one.

The Association of ARMC5 with the Renin-Angiotensin-Aldosterone System, Blood Pressure, and Glycemia in African Americans

CONTEXT

Armadillo repeat containing 5 (ARMC5) on chromosome 16 is an adrenal gland tumor suppressor gene associated with primary aldosteronism, especially among African Americans (AAs). We examined the association of ARMC5 variants with aldosterone, plasma renin activity (PRA), blood pressure, glucose, and glycosylated hemoglobin A1c (HbA1c) in community-dwelling AAs.

METHODS

The Jackson Heart Study is a prospective cardiovascular cohort study in AAs with baseline data collection from 2000 to 2004. Kernel machine method was used to perform a single joint test to analyze for an overall association between the phenotypes of interest (aldosterone, PRA, systolic and diastolic blood pressure [SBP, DBP], glucose, and HbA1c) and the ARMC5 single nucleotide variants (SNVs) adjusted for age, sex, BMI, and medications; followed by Baysian Lasso methodology to identify sets of SNVs in terms of associated haplotypes with specific phenotypes.

RESULTS

Among 3223 participants (62% female; mean age 55.6 (SD ± 12.8) years), the average SBP and DBP were 127 and 76 mmHg, respectively. The average fasting plasma glucose and HbA1c were 101 mg/dL and 6.0%, respectively. ARMC5 variants were associated with all 6 phenotypes. Haplotype TCGCC (ch16:31476015-31476093) was negatively associated, whereas haplotype CCCCTTGCG (ch16:31477195-31477460) was positively associated with SBP, DBP, and glucose. Haplotypes GGACG (ch16:31477790-31478013) and ACGCG (ch16:31477834-31478113) were negatively associated with aldosterone and positively associated with HbA1c and glucose, respectively. Haplotype GCGCGAGC (ch16:31471193-ch16:31473597(rs114871627) was positively associated with PRA and negatively associated with HbA1c.

CONCLUSIONS

ARMC5 variants are associated with aldosterone, PRA, blood pressure, fasting glucose, and HbA1c in community-dwelling AAs, suggesting that germline mutations in ARMC5 may underlie cardiometabolic disease in AAs.

A Family-Based Rare Haplotype Association Method for Quantitative Traits

BACKGROUND

The variants identified in genome-wide association studies account for only a small fraction of disease heritability. A key to this "missing heritability" is believed to be rare variants. Specifically, we focus on rare haplotype variant (rHTV). The existing methods for detecting rHTV are mostly population-based, and as such, are susceptible to population stratification and admixture, leading to an inflated false-positive rate. Family-based methods are more robust in this respect.

METHODS

We propose a method for detecting rHTVs associated with quantitative traits called family-based quantitative Bayesian LASSO (famQBL). FamQBL can analyze any type of pedigree and is based on a mixed model framework. We regularize the haplotype effects using Bayesian LASSO and estimate the posterior distributions using Markov chain Monte Carlo methods.

RESULTS

We conduct simulation studies, including analyses of Genetic Analysis Workshop 18 simulated data, to study the properties of famQBL and compare with a standard family-based haplotype association test implemented in FBAT (family-based association test) software. We find famQBL to be more powerful than FBAT with well-controlled false-positive rates. We also apply famQBL to the Framingham Heart Study data and detect an rHTV associated with diastolic blood pressure.

CONCLUSION

FamQBL can help uncover rHTVs associated with quantitative traits.

Methods and results from the genome-wide association group at GAW20

BACKGROUND

This paper summarizes the contributions from the Genome-wide Association Study group (GWAS group) of the GAW20. The GWAS group contributions focused on topics such as association tests, phenotype imputation, and application of empirical kinships. The goals of the GWAS group contributions were varied. A real or a simulated data set based on the Genetics of Lipid Lowering Drugs and Diet Network (GOLDN) study was employed by different methods. Different outcomes and covariates were considered, and quality control procedures varied throughout the contributions.

RESULTS

The consideration of heritability and family structure played a major role in some contributions. The inclusion of family information and adaptive weights based on data were found to improve power in genome-wide association studies. It was proven that gene-level approaches are more powerful than single-marker analysis. Other contributions focused on the comparison between pedigree-based kinship and empirical kinship matrices, and investigated similar results in heritability estimation, association mapping, and genomic prediction. A new approach for linkage mapping of triglyceride levels was able to identify a novel linkage signal.

CONCLUSIONS

This summary paper reports on promising statistical approaches and findings of the members of the GWAS group applied on real and simulated data which encompass the current topics of epigenetic and pharmacogenomics.