A multi-population phenome-wide association study of genetically-predicted height in the Million Veteran Program.

BACKGROUND: Height has been associated with many clinical traits but whether such associations are causal versus secondary to confounding remains unclear in many cases. To systematically examine this question, we performed a Mendelian Randomization-Phenome-wide association study (MR-PheWAS) using clinical and genetic data from a national healthcare system biobank. METHODS AND FINDINGS: Analyses were performed using data from the US Veterans Affairs (VA) Million Veteran Program in non-Hispanic White (EA, n = 222,300) and non-Hispanic Black (AA, n = 58,151) adults in the US. We estimated height genetic risk based on 3290 height-associated variants from a recent European-ancestry genome-wide meta-analysis. We compared associations of measured and genetically-predicted height with phenome-wide traits derived from the VA electronic health record, adjusting for age, sex, and genetic principal components. We found 345 clinical traits associated with measured height in EA and an additional 17 in AA. Of these, 127 were associated with genetically-predicted height at phenome-wide significance in EA and 2 in AA. These associations were largely independent from body mass index. We confirmed several previously described MR associations between height and cardiovascular disease traits such as hypertension, hyperlipidemia, coronary heart disease (CHD), and atrial fibrillation, and further uncovered MR associations with venous circulatory disorders and peripheral neuropathy in the presence and absence of diabetes. As a number of traits associated with genetically-predicted height frequently co-occur with CHD, we evaluated effect modification by CHD status of genetically-predicted height associations with risk factors for and complications of CHD. We found modification of effects of MR associations by CHD status for atrial fibrillation/flutter but not for hypertension, hyperlipidemia, or venous circulatory disorders. CONCLUSIONS: We conclude that height may be an unrecognized but biologically plausible risk factor for several common conditions in adults. However, more studies are needed to reliably exclude horizontal pleiotropy as a driving force behind at least some of the MR associations observed in this study.

Mendelian randomization in the multivariate general linear model framework.

Mendelian randomization (MR) is an application of instrumental variable (IV) methods to observational data in which the IV is a genetic variant. MR methods applicable to the general exponential family of distributions are currently not well characterized. We adapt a general linear model framework to the IV setting and propose a general MR method applicable to any full-rank distribution from the exponential family. Empirical bias and coverage are estimated via simulations. The proposed method is compared to several existing MR methods. Real data analyses are performed using data from the REGARDS study to estimate the potential causal effect of smoking frequency on stroke risk in African Americans. In simulations with binary variates and very weak instruments the proposed method had the lowest median [Q1 , Q3 ] bias (0.10 [-3.68 to 3.62]); compared with 2SPS (0.27 [-3.74 to 4.26]) and the Wald method (-0.69 [-1.72 to 0.35]). Low bias was observed throughout other simulation scenarios; as well as more than 90% coverage for the proposed method. In simulations with count variates, the proposed method performed comparably to 2SPS; the Wald method maintained the most consistent low bias; and 2SRI was biased towards the null. Real data analyses find no evidence for a causal effect of smoking frequency on stroke risk. The proposed MR method has low bias and acceptable coverage across a wide range of distributional scenarios and instrument strengths; and provides a more parsimonious framework for asymptotic hypothesis testing compared to existing two-stage procedures.

Gene expression associations with body mass index in the Multi-Ethnic Study of Atherosclerosis.

BACKGROUND/OBJECTIVES: Obesity, defined as excessive fat accumulation that represents a health risk, is increasing in adults and children, reaching global epidemic proportions. Body mass index (BMI) correlates with body fat and future health risk, yet differs in prediction by fat distribution, across populations and by age. Nonetheless, few genetic studies of BMI have been conducted in ancestrally diverse populations. Gene expression association with BMI was assessed in the Multi-Ethnic Study of Atherosclerosis (MESA) in four self-identified race and ethnicity (SIRE) groups to identify genes associated with obesity. SUBJECTS/METHODS: RNA-sequencing was performed on 1096 MESA participants (37.8% white, 24.3% Hispanic, 28.4% African American, and 9.5% Chinese American) and linear models were used to assess the association of expression from each gene for its effect on BMI, adjusting for age, sex, sequencing center, study site, five expression and four genetic principal components in each self-identified race group. Sample-size-weighted meta-analysis was performed to identify genes with BMI-associated expression across ancestry groups. RESULTS: Within individual SIRE groups, there were zero to three genes whose expression is significantly (p < 1.97 × 10-6) associated with BMI. Across all groups, 45 genes were identified by meta-analysis whose expression was significantly associated with BMI, explaining 29.7% of BMI variation. The 45 genes are expressed in a variety of tissues and cell types and are enriched for obesity-related processes including erythrocyte function, oxygen binding and transport, and JAK-STAT signaling. CONCLUSIONS: We have identified genes whose expression is significantly associated with obesity in a multi-ethnic cohort. We have identified novel genes associated with BMI as well as confirmed previously identified genes from earlier genetic analyses. These novel genes and their biological pathways represent new targets for understanding the biology of obesity as well as new therapeutic intervention to reduce obesity and improve global public health.

Plasma neurofilament light chain concentrations are elevated in youth-onset type 2 diabetes and associated with neuropathy.

BACKGROUND AND AIMS: The lack of easily measurable biomarkers remains a challenge in executing clinical trials for diabetic neuropathy (DN). Plasma Neurofilament light chain (NFL) concentration is a promising biomarker in immune-mediated neuropathies. Longitudinal studies evaluating NFL in DN have not been performed. METHODS: A nested case-control study was performed on participants with youth-onset type 2 diabetes enrolled in the prospective Treatment Options for Type 2 Diabetes in Adolescents and Youth (TODAY) study. Plasma NFL concentrations were measured at 4-year intervals from 2008 to 2020 in 50 participants who developed DN and 50 participants with type 2 diabetes who did not develop DN. RESULTS: NFL concentrations were similar in the DN and no DN groups at the first assessment. Concentrations were higher in DN participants at all subsequent assessment periods (all p < .01). NFL concentrations increased over time in both groups, with higher degrees of change in DN participants (interaction p = .045). A doubling of the NFL value at Assessment 2 in those without DN increased the odds of ultimate DN outcome by an estimated ratio of 2.86 (95% CI: [1.30, 6.33], p = .0046). At the final study visit, positive Spearman correlations (controlled for age, sex, diabetes duration, and BMI) were observed between NFL and HbA1c (0.48, p < .0001), total cholesterol (0.25, p = .018), and low-density lipoprotein (LDL (0.30, p = .0037)). Negative correlations were observed with measures of heart rate variability (-0.42 to -0.46, p = <.0001). INTERPRETATION: The findings that NFL concentrations are elevated in individuals with youth-onset type 2 diabetes, and increase more rapidly in those who develop DN, suggest that NFL could be a valuable biomarker for DN.

Mendelian randomization study of diabetes and dementia in the Million Veteran Program.

INTRODUCTION: Diabetes and dementia are diseases of high health-care burden worldwide. Individuals with diabetes have 1.4 to 2.2 times higher risk of dementia. Our objective was to evaluate evidence of causality between these two common diseases. METHODS: We conducted a one-sample Mendelian randomization (MR) analysis in the US Department of Veterans Affairs Million Veteran program. The study included 334,672 participants ≥65 years of age with type 2 diabetes and dementia case-control status and genotype data. RESULTS: For each standard deviation increase in genetically predicted diabetes, we found increased odds of three dementia diagnoses in non-Hispanic White participants (all-cause: odds ratio [OR] = 1.07 [1.05-1.08], P = 3.40E-18; vascular: OR = 1.11 [1.07-1.15], P = 3.63E-09, Alzheimer's disease [AD]: OR = 1.06 [1.02-1.09], P = 6.84E-04) and non-Hispanic Black participants (all-cause: OR = 1.06 [1.02-1.10], P = 3.66E-03, vascular: OR = 1.11 [1.04-1.19], P = 2.20E-03, AD: OR = 1.12 [1.02-1.23], P = 1.60E-02) but not in Hispanic participants (all P > 0.05). DISCUSSION: We found evidence of causality between diabetes and dementia using a one-sample MR study, with access to individual level data, overcoming limitations of prior studies using two-sample MR techniques.

Soluble CD14 Levels in the Jackson Heart Study: Associations With Cardiovascular Disease Risk and Genetic Variants.

[Figure: see text].

Exome sequencing identifies rare LDLR and APOA5 alleles conferring risk for myocardial infarction.

Myocardial infarction (MI), a leading cause of death around the world, displays a complex pattern of inheritance. When MI occurs early in life, genetic inheritance is a major component to risk. Previously, rare mutations in low-density lipoprotein (LDL) genes have been shown to contribute to MI risk in individual families, whereas common variants at more than 45 loci have been associated with MI risk in the population. Here we evaluate how rare mutations contribute to early-onset MI risk in the population. We sequenced the protein-coding regions of 9,793 genomes from patients with MI at an early age (≤50 years in males and ≤60 years in females) along with MI-free controls. We identified two genes in which rare coding-sequence mutations were more frequent in MI cases versus controls at exome-wide significance. At low-density lipoprotein receptor (LDLR), carriers of rare non-synonymous mutations were at 4.2-fold increased risk for MI; carriers of null alleles at LDLR were at even higher risk (13-fold difference). Approximately 2% of early MI cases harbour a rare, damaging mutation in LDLR; this estimate is similar to one made more than 40 years ago using an analysis of total cholesterol. Among controls, about 1 in 217 carried an LDLR coding-sequence mutation and had plasma LDL cholesterol > 190 mg dl(-1). At apolipoprotein A-V (APOA5), carriers of rare non-synonymous mutations were at 2.2-fold increased risk for MI. When compared with non-carriers, LDLR mutation carriers had higher plasma LDL cholesterol, whereas APOA5 mutation carriers had higher plasma triglycerides. Recent evidence has connected MI risk with coding-sequence mutations at two genes functionally related to APOA5, namely lipoprotein lipase and apolipoprotein C-III (refs 18, 19). Combined, these observations suggest that, as well as LDL cholesterol, disordered metabolism of triglyceride-rich lipoproteins contributes to MI risk.

Identification and functional analysis of glycemic trait loci in the China Health and Nutrition Survey.

To identify genetic contributions to type 2 diabetes (T2D) and related glycemic traits (fasting glucose, fasting insulin, and HbA1c), we conducted genome-wide association analyses (GWAS) in up to 7,178 Chinese subjects from nine provinces in the China Health and Nutrition Survey (CHNS). We examined patterns of population structure within CHNS and found that allele frequencies differed across provinces, consistent with genetic drift and population substructure. We further validated 32 previously described T2D- and glycemic trait-loci, including G6PC2 and SIX3-SIX2 associated with fasting glucose. At G6PC2, we replicated a known fasting glucose-associated variant (rs34177044) and identified a second signal (rs2232326), a low-frequency (4%), probably damaging missense variant (S324P). A variant within the lead fasting glucose-associated signal at SIX3-SIX2 co-localized with pancreatic islet expression quantitative trait loci (eQTL) for SIX3, SIX2, and three noncoding transcripts. To identify variants functionally responsible for the fasting glucose association at SIX3-SIX2, we tested five candidate variants for allelic differences in regulatory function. The rs12712928-C allele, associated with higher fasting glucose and lower transcript expression level, showed lower transcriptional activity in reporter assays and increased binding to GABP compared to the rs12712928-G, suggesting that rs12712928-C contributes to elevated fasting glucose levels by disrupting an islet enhancer, resulting in reduced gene expression. Taken together, these analyses identified multiple loci associated with glycemic traits across China, and suggest a regulatory mechanism at the SIX3-SIX2 fasting glucose GWAS locus.

Common α-globin variants modify hematologic and other clinical phenotypes in sickle cell trait and disease.

Co-inheritance of α-thalassemia has a significant protective effect on the severity of complications of sickle cell disease (SCD), including stroke. However, little information exists on the association and interactions for the common African ancestral α-thalassemia mutation (-α3.7 deletion) and ß-globin traits (HbS trait [SCT] and HbC trait) on important clinical phenotypes such as red blood cell parameters, anemia, and chronic kidney disease (CKD). In a community-based cohort of 2,916 African Americans from the Jackson Heart Study, we confirmed the expected associations between SCT, HbC trait, and the -α3.7 deletion with lower mean corpuscular volume/mean corpuscular hemoglobin and higher red blood cell count and red cell distribution width. In addition to the recently recognized association of SCT with lower estimated glomerular filtration rate and glycated hemoglobin (HbA1c), we observed a novel association of the -α3.7 deletion with higher HbA1c levels. Co-inheritance of each additional copy of the -α3.7 deletion significantly lowered the risk of anemia and chronic kidney disease among individuals with SCT (P-interaction = 0.031 and 0.019, respectively). Furthermore, co-inheritance of a novel α-globin regulatory variant was associated with normalization of red cell parameters in individuals with the -α3.7 deletion and significantly negated the protective effect of α-thalassemia on stroke in 1,139 patients with sickle cell anemia from the Cooperative Study of Sickle Cell Disease (CSSCD) (P-interaction = 0.0049). Functional assays determined that rs11865131, located in the major alpha-globin enhancer MCS-R2, was the most likely causal variant. These findings suggest that common α- and ß-globin variants interact to influence hematologic and clinical phenotypes in African Americans, with potential implications for risk-stratification and counseling of individuals with SCD and SCT.

Genetic Risk for Hepatic Fat among an Ethnically Diverse Cohort of Youth: The Exploring Perinatal Outcomes among Children Study.

OBJECTIVE: To assess the importance of genetic and nongenetic risk factors contributing to hepatic fat accumulation in a multiethnic population of youth. STUDY DESIGN: We investigated the relationship between genetic factors and hepatic fat fraction (HFF) in 347 children aged 12.5-19.5 years. We examined 5 single nucleotide polymorphisms previously associated with HFF and a weighted genetic risk score (GRS) and examined how these associations varied with ethnicity (Hispanic vs non-Hispanic white) and body mass index (BMI) category. We also compared how much variation in HFF was explained by genetic factors vs cardiometabolic factors (BMI z-score and the Homeostasis Model of Insulin Resistance) or diet. RESULTS: PNPLA3 rs738409 and the GRS were each associated with HFF among Hispanic (ß = 0.39; 95% CI, 0.16-0.62; P = .001; and ß = 0.20; 95% CI, 0.05-0.34; P = .007, respectively) but not non-Hispanic white (ß = 0.04; 95% CI, -0.18 to 0.26; P = .696; and ß = 0.03; 95% CI, -0.09 to 0.14; P = .651, respectively) youth. Cardiometabolic risk factors explained more of the variation in HFF than genetic risk factors among non-lean Hispanic individuals (27.2% for cardiometabolic markers vs 6.4% for rs738409 and 4.3% for the GRS), and genetic risk factors were more important among lean individuals (2.7% for cardiometabolic markers vs 12.6% for rs738409 and 4.4% for the GRS). CONCLUSIONS: Poor cardiometabolic health may be more important than genetic factors when predicting HFF in overweight and obese young populations. Genetic risk is an important contributor to pediatric HFF among lean Hispanics, but further studies are necessary to elucidate the strength of the association between genetic risk and HFF in non-Hispanic white youth.

A robust and powerful two-step testing procedure for local ancestry adjusted allelic association analysis in admixed populations.

Genetic association studies in admixed populations allow us to gain deeper understanding of the genetic architecture of human diseases and traits. However, population stratification, complicated linkage disequilibrium (LD) patterns, and the complex interplay of allelic and ancestry effects on phenotypic traits pose challenges in such analyses. These issues may lead to detecting spurious associations and/or result in reduced statistical power. Fortunately, if handled appropriately, these same challenges provide unique opportunities for gene mapping. To address these challenges and to take these opportunities, we propose a robust and powerful two-step testing procedure Local Ancestry Adjusted Allelic (LAAA) association. In the first step, LAAA robustly captures associations due to allelic effect, ancestry effect, and interaction effect, allowing detection of effect heterogeneity across ancestral populations. In the second step, LAAA identifies the source of association, namely allelic, ancestry, or the combination. By jointly modeling allele, local ancestry, and ancestry-specific allelic effects, LAAA is highly powerful in capturing the presence of interaction between ancestry and allele effect. We evaluated the validity and statistical power of LAAA through simulations over a broad spectrum of scenarios. We further illustrated its usefulness by application to the Candidate Gene Association Resource (CARe) African American participants for association with hemoglobin levels. We were able to replicate independent groups' previously identified loci that would have been missed in CARe without joint testing. Moreover, the loci, for which LAAA detected potential effect heterogeneity, were replicated among African Americans from the Women's Health Initiative study. LAAA is freely available at https://yunliweb.its.unc.edu/LAAA.

REVEL: An Ensemble Method for Predicting the Pathogenicity of Rare Missense Variants.

The vast majority of coding variants are rare, and assessment of the contribution of rare variants to complex traits is hampered by low statistical power and limited functional data. Improved methods for predicting the pathogenicity of rare coding variants are needed to facilitate the discovery of disease variants from exome sequencing studies. We developed REVEL (rare exome variant ensemble learner), an ensemble method for predicting the pathogenicity of missense variants on the basis of individual tools: MutPred, FATHMM, VEST, PolyPhen, SIFT, PROVEAN, MutationAssessor, MutationTaster, LRT, GERP, SiPhy, phyloP, and phastCons. REVEL was trained with recently discovered pathogenic and rare neutral missense variants, excluding those previously used to train its constituent tools. When applied to two independent test sets, REVEL had the best overall performance (p < 10-12) as compared to any individual tool and seven ensemble methods: MetaSVM, MetaLR, KGGSeq, Condel, CADD, DANN, and Eigen. Importantly, REVEL also had the best performance for distinguishing pathogenic from rare neutral variants with allele frequencies <0.5%. The area under the receiver operating characteristic curve (AUC) for REVEL was 0.046-0.182 higher in an independent test set of 935 recent SwissVar disease variants and 123,935 putatively neutral exome sequencing variants and 0.027-0.143 higher in an independent test set of 1,953 pathogenic and 2,406 benign variants recently reported in ClinVar than the AUCs for other ensemble methods. We provide pre-computed REVEL scores for all possible human missense variants to facilitate the identification of pathogenic variants in the sea of rare variants discovered as sequencing studies expand in scale.

Exome Genotyping Identifies Pleiotropic Variants Associated with Red Blood Cell Traits.

Red blood cell (RBC) traits are important heritable clinical biomarkers and modifiers of disease severity. To identify coding genetic variants associated with these traits, we conducted meta-analyses of seven RBC phenotypes in 130,273 multi-ethnic individuals from studies genotyped on an exome array. After conditional analyses and replication in 27,480 independent individuals, we identified 16 new RBC variants. We found low-frequency missense variants in MAP1A (rs55707100, minor allele frequency [MAF] = 3.3%, p = 2 × 10(-10) for hemoglobin [HGB]) and HNF4A (rs1800961, MAF = 2.4%, p < 3 × 10(-8) for hematocrit [HCT] and HGB). In African Americans, we identified a nonsense variant in CD36 associated with higher RBC distribution width (rs3211938, MAF = 8.7%, p = 7 × 10(-11)) and showed that it is associated with lower CD36 expression and strong allelic imbalance in ex vivo differentiated human erythroblasts. We also identified a rare missense variant in ALAS2 (rs201062903, MAF = 0.2%) associated with lower mean corpuscular volume and mean corpuscular hemoglobin (p < 8 × 10(-9)). Mendelian mutations in ALAS2 are a cause of sideroblastic anemia and erythropoietic protoporphyria. Gene-based testing highlighted three rare missense variants in PKLR, a gene mutated in Mendelian non-spherocytic hemolytic anemia, associated with HGB and HCT (SKAT p < 8 × 10(-7)). These rare, low-frequency, and common RBC variants showed pleiotropy, being also associated with platelet, white blood cell, and lipid traits. Our association results and functional annotation suggest the involvement of new genes in human erythropoiesis. We also confirm that rare and low-frequency variants play a role in the architecture of complex human traits, although their phenotypic effect is generally smaller than originally anticipated.

Association of sickle cell trait with atrial fibrillation: The REGARDS cohort.

BACKGROUND: Sickle cell trait (SCT), sickle cell disease's (SCD) carrier status, has been recently associated with worse cardiovascular and renal outcomes. An increased prevalence of atrial fibrillation (AF) is documented in SCD patients; however, studies in individuals with SCT are lacking. OBJECTIVES: To determine the association of SCT with AF. METHODS: Among African-American participants in the REasons for Geographic and Racial Differences in Stroke (REGARDS) Study we assessed the association of SCT (by ECG or medical history) with prevalent AF using logistic regression adjusting for age, sex, income, education, history of stroke, myocardial infarction, diabetes, hypertension, and chronic kidney disease. A second evaluation was performed a mean of 9.2 years later among available participants, and the same model was used to test the association of SCT with incident AF. RESULTS: In 10,409 participants with baseline ECG data and genotyping, 778 (7.5%) had SCT and 811 (7.8%) had prevalent AF. After adjusting for age, sex, education and income, SCT was associated with AF, OR 1.32 (95% CI 1.03-1.70). The association with incident AF assessed at the second in-home visit with the same adjustments was similar; OR 1.25 (95% CI 0.77-2.03). CONCLUSIONS: SCT was associated with a higher prevalence of AF and a non-significantly higher incident AF over a 9.2 year period independent of AF risk factors. SCT remained associated with prevalent AF after adjusting for potential factors on the causal pathway such as hypertension and chronic kidney disease suggesting alternate mechanisms for the increased risk.

gsSKAT: Rapid gene set analysis and multiple testing correction for rare-variant association studies using weighted linear kernels.

Next-generation sequencing technologies have afforded unprecedented characterization of low-frequency and rare genetic variation. Due to low power for single-variant testing, aggregative methods are commonly used to combine observed rare variation within a single gene. Causal variation may also aggregate across multiple genes within relevant biomolecular pathways. Kernel-machine regression and adaptive testing methods for aggregative rare-variant association testing have been demonstrated to be powerful approaches for pathway-level analysis, although these methods tend to be computationally intensive at high-variant dimensionality and require access to complete data. An additional analytical issue in scans of large pathway definition sets is multiple testing correction. Gene set definitions may exhibit substantial genic overlap, and the impact of the resultant correlation in test statistics on Type I error rate control for large agnostic gene set scans has not been fully explored. Herein, we first outline a statistical strategy for aggregative rare-variant analysis using component gene-level linear kernel score test summary statistics as well as derive simple estimators of the effective number of tests for family-wise error rate control. We then conduct extensive simulation studies to characterize the behavior of our approach relative to direct application of kernel and adaptive methods under a variety of conditions. We also apply our method to two case-control studies, respectively, evaluating rare variation in hereditary prostate cancer and schizophrenia. Finally, we provide open-source R code for public use to facilitate easy application of our methods to existing rare-variant analysis results.

Genome-wide association study of homocysteine in African Americans from the Jackson Heart Study, the Multi-Ethnic Study of Atherosclerosis, and the Coronary Artery Risk in Young Adults study.

Homocysteine (Hcy) is a heritable biomarker for CVD, peripheral artery disease, stroke, and dementia. Little is known about genetic associations with Hcy in individuals of African ancestry. We performed a genome-wide association study for Hcy in 4927 AAs from the Jackson Heart Study (JHS), the Multi-Ethnic Study of Atherosclerosis (MESA), and the Coronary Artery Risk in Young Adults (CARDIA) study. Analyses were stratified by sex and results were meta-analyzed within and across sex. In the sex-combined meta-analysis, we observed genome-wide significant evidence (p < 5.0 × 10-8) for the NOX4 locus (lead variant rs2289125, ß = -0.15, p = 5.3 × 1011). While the NOX4 locus was previously reported as associated with Hcy in European-American populations, rs2289125 remained genome-wide significant when conditioned on the previously reported lead variants. Previously reported genome-wide significant associations at NOX4, MTR, CBS, and MMACHC were also nominally (p < 0.050) replicated in AAs. Associations at the CPS1 locus, previously reported in females only, also was replicated specifically in females in this analysis, supporting sex-specific effects for this locus. These results suggest that there may be a combination of cross-population and population-specific genetic effects, as well as differences in genetic effects between males and females, in the regulation of Hcy levels.

D-Dimer in African Americans: Whole Genome Sequence Analysis and Relationship to Cardiovascular Disease Risk in the Jackson Heart Study.

OBJECTIVE: Plasma levels of the fibrinogen degradation product D-dimer are higher among African Americans (AAs) compared with those of European ancestry and higher among women compared with men. Among AAs, little is known of the genetic architecture of D-dimer or the relationship of D-dimer to incident cardiovascular disease. APPROACH AND RESULTS: We measured baseline D-dimer in 4163 AAs aged 21 to 93 years from the prospective JHS (Jackson Heart Study) cohort and assessed association with incident cardiovascular disease events. In participants with whole genome sequencing data (n=2980), we evaluated common and rare genetic variants for association with D-dimer. Each standard deviation higher baseline D-dimer was associated with a 20% to 30% increased hazard for incident coronary heart disease, stroke, and all-cause mortality. Genetic variation near F3 was associated with higher D-dimer (rs2022030, ß=0.284, P=3.24×10-11). The rs2022030 effect size was nearly 3× larger among women (ß=0.373, P=9.06×10-13) than among men (ß=0.135, P=0.06; P interaction =0.009). The sex by rs2022030 interaction was replicated in an independent sample of 10 808 multiethnic men and women (P interaction =0.001). Finally, the African ancestral sickle cell variant (HBB rs334) was significantly associated with higher D-dimer in JHS (ß=0.507, P=1.41×10-14), and this association was successfully replicated in 1933 AAs (P=2.3×10-5). CONCLUSIONS: These results highlight D-dimer as an important predictor of cardiovascular disease risk in AAs and suggest that sex-specific and African ancestral genetic effects of the F3 and HBB loci contribute to the higher levels of D-dimer among women and AAs.

Post hoc Analysis for Detecting Individual Rare Variant Risk Associations Using Probit Regression Bayesian Variable Selection Methods in Case-Control Sequencing Studies.

Rare variants (RVs) have been shown to be significant contributors to complex disease risk. By definition, these variants have very low minor allele frequencies and traditional single-marker methods for statistical analysis are underpowered for typical sequencing study sample sizes. Multimarker burden-type approaches attempt to identify aggregation of RVs across case-control status by analyzing relatively small partitions of the genome, such as genes. However, it is generally the case that the aggregative measure would be a mixture of causal and neutral variants, and these omnibus tests do not directly provide any indication of which RVs may be driving a given association. Recently, Bayesian variable selection approaches have been proposed to identify RV associations from a large set of RVs under consideration. Although these approaches have been shown to be powerful at detecting associations at the RV level, there are often computational limitations on the total quantity of RVs under consideration and compromises are necessary for large-scale application. Here, we propose a computationally efficient alternative formulation of this method using a probit regression approach specifically capable of simultaneously analyzing hundreds to thousands of RVs. We evaluate our approach to detect causal variation on simulated data and examine sensitivity and specificity in instances of high RV dimensionality as well as apply it to pathway-level RV analysis results from a prostate cancer (PC) risk case-control sequencing study. Finally, we discuss potential extensions and future directions of this work.

Genome-wide admixture and association study of serum iron, ferritin, transferrin saturation and total iron binding capacity in African Americans.

Iron is an essential component of many important proteins and enzymes, including hemoglobin, which is responsible for carrying oxygen to the cells. African Americans (AAs) have a greater prevalence of iron deficiency compared with European Americans. We conducted genome-wide admixture-mapping and association studies for serum iron, serum ferritin, transferrin saturation (SAT) and total iron binding capacity (TIBC) in 2347 AAs participating in the Jackson Heart Study (JHS). Follow-up replication analyses for JHS iron-trait associated SNPs were conducted in 329 AA participants in the Healthy Aging in Neighborhoods of Diversity across the Life Span study (HANDLS). Higher estimated proportions of global African ancestry were significantly associated with lower levels of iron (P = 2.4 × 10(-5)), SAT (P = 0.0019) and TIBC (P = 0.042). We observed significant associations (P < 5 × 10(-8)) between serum TIBC levels and two independent SNPs around TF on chromosome 3, the first report of a genome-wide significant second independent signal in this region, and SNPs near two novel genes: HDGFL1 on chromosome 6 and MAF on chromosome 16. We also observed significant associations between ferritin levels and SNPs near GAB3 on chromosome X. We replicated our two independent associations at TF and our association at GAB3 in HANDLS. Our study provides evidence for both shared and unique genetic risk factors that are associated with iron-related measures in AAs. The top two variants in TF explain 11.2% of the total variation in TIBC levels in AAs after accounting for age, gender, body mass index and background ancestry.

Whole-exome sequencing identifies rare and low-frequency coding variants associated with LDL cholesterol.

Elevated low-density lipoprotein cholesterol (LDL-C) is a treatable, heritable risk factor for cardiovascular disease. Genome-wide association studies (GWASs) have identified 157 variants associated with lipid levels but are not well suited to assess the impact of rare and low-frequency variants. To determine whether rare or low-frequency coding variants are associated with LDL-C, we exome sequenced 2,005 individuals, including 554 individuals selected for extreme LDL-C (>98(th) or <2(nd) percentile). Follow-up analyses included sequencing of 1,302 additional individuals and genotype-based analysis of 52,221 individuals. We observed significant evidence of association between LDL-C and the burden of rare or low-frequency variants in PNPLA5, encoding a phospholipase-domain-containing protein, and both known and previously unidentified variants in PCSK9, LDLR and APOB, three known lipid-related genes. The effect sizes for the burden of rare variants for each associated gene were substantially higher than those observed for individual SNPs identified from GWASs. We replicated the PNPLA5 signal in an independent large-scale sequencing study of 2,084 individuals. In conclusion, this large whole-exome-sequencing study for LDL-C identified a gene not known to be implicated in LDL-C and provides unique insight into the design and analysis of similar experiments.