Genetic drivers of heterogeneity in type 2 diabetes pathophysiology

Type 2 diabetes (T2D) is a heterogeneous disease that develops through diverse pathophysiological processes1,2 and molecular mechanisms that are often specific to cell type3,4. Here, to characterize the genetic contribution to these processes across ancestry groups, we aggregate genome-wide association study data from 2,535,601 individuals (39.7% not of European ancestry), including 428,452 cases of T2D. We identify 1,289 independent association signals at genome-wide significance (P < 5 × 10-8) that map to 611 loci, of which 145 loci are, to our knowledge, previously unreported. We define eight non-overlapping clusters of T2D signals that are characterized by distinct profiles of cardiometabolic trait associations. These clusters are differentially enriched for cell-type-specific regions of open chromatin, including pancreatic islets, adipocytes, endothelial cells and enteroendocrine cells. We build cluster-specific partitioned polygenic scores5 in a further 279,552 individuals of diverse ancestry, including 30,288 cases of T2D, and test their association with T2D-related vascular outcomes. Cluster-specific partitioned polygenic scores are associated with coronary artery disease, peripheral artery disease and end-stage diabetic nephropathy across ancestry groups, highlighting the importance of obesity-related processes in the development of vascular outcomes. Our findings show the value of integrating multi-ancestry genome-wide association study data with single-cell epigenomics to disentangle the aetiological heterogeneity that drives the development and progression of T2D. This might offer a route to optimize global access to genetically informed diabetes care.

Gene-educational attainment interactions in a multi-population genome-wide meta-analysis identify novel lipid loci

Introduction: Educational attainment, widely used in epidemiologic studies as a surrogate for socioeconomic status, is a predictor of cardiovascular health outcomes. Methods: A two-stage genome-wide meta-analysis of low-density lipoprotein cholesterol (LDL), high-density lipoprotein cholesterol (HDL), and triglyceride (TG) levels was performed while accounting for gene-educational attainment interactions in up to 226,315 individuals from five population groups. We considered two educational attainment variables: "Some College" (yes/no, for any education beyond high school) and "Graduated College" (yes/no, for completing a 4-year college degree). Genome-wide significant (p < 5 × 10-8) and suggestive (p < 1 × 10-6) variants were identified in Stage 1 (in up to 108,784 individuals) through genome-wide analysis, and those variants were followed up in Stage 2 studies (in up to 117,531 individuals). Results: In combined analysis of Stages 1 and 2, we identified 18 novel lipid loci (nine for LDL, seven for HDL, and two for TG) by two degree-of-freedom (2 DF) joint tests of main and interaction effects. Four loci showed significant interaction with educational attainment. Two loci were significant only in cross-population analyses. Several loci include genes with known or suggested roles in adipose (FOXP1, MBOAT4, SKP2, STIM1, STX4), brain (BRI3, FILIP1, FOXP1, LINC00290, LMTK2, MBOAT4, MYO6, SENP6, SRGAP3, STIM1, TMEM167A, TMEM30A), and liver (BRI3, FOXP1) biology, highlighting the potential importance of brain-adipose-liver communication in the regulation of lipid metabolism. An investigation of the potential druggability of genes in identified loci resulted in five gene targets shown to interact with drugs approved by the Food and Drug Administration, including genes with roles in adipose and brain tissue. Discussion: Genome-wide interaction analysis of educational attainment identified novel lipid loci not previously detected by analyses limited to main genetic effects.

Multi-ancestry genome-wide study identifies effector genes and druggable pathways for coronary artery calcification

Coronary artery calcification (CAC), a measure of subclinical atherosclerosis, predicts future symptomatic coronary artery disease (CAD). Identifying genetic risk factors for CAC may point to new therapeutic avenues for prevention. Currently, there are only four known risk loci for CAC identified from genome-wide association studies (GWAS) in the general population. Here we conducted the largest multi-ancestry GWAS meta-analysis of CAC to date, which comprised 26,909 individuals of European ancestry and 8,867 individuals of African ancestry. We identified 11 independent risk loci, of which eight were new for CAC and five had not been reported for CAD. These new CAC loci are related to bone mineralization, phosphate catabolism and hormone metabolic pathways. Several new loci harbor candidate causal genes supported by multiple lines of functional evidence and are regulators of smooth muscle cell-mediated calcification ex vivo and in vitro. Together, these findings help refine the genetic architecture of CAC and extend our understanding of the biological and potential druggable pathways underlying CAC.

WHOLE GENOME SEQUENCING ANALYSIS OF BODY MASS INDEX IDENTIFIES NOVEL AFRICAN ANCESTRY-SPECIFIC RISK ALLELE

Obesity is a major public health crisis associated with high mortality rates. Previous genome-wide association studies (GWAS) investigating body mass index (BMI) have largely relied on imputed data from European individuals. This study leveraged whole-genome sequencing (WGS) data from 88,873 participants from the Trans-Omics for Precision Medicine (TOPMed) Program, of which 51% were of non-European population groups. We discovered 18 BMI-associated signals (P < 5 × 10-9). Notably, we identified and replicated a novel low frequency single nucleotide polymorphism (SNP) in MTMR3 that was common in individuals of African descent. Using a diverse study population, we further identified two novel secondary signals in known BMI loci and pinpointed two likely causal variants in the POC5 and DMD loci. Our work demonstrates the benefits of combining WGS and diverse cohorts in expanding current catalog of variants and genes confer risk for obesity, bringing us one step closer to personalized medicine.

Multi-ancestry genome-wide study in >2.5 million individuals reveals heterogeneity in mechanistic pathways of type 2 diabetes and complications

Type 2 diabetes (T2D) is a heterogeneous disease that develops through diverse pathophysiological processes. To characterise the genetic contribution to these processes across ancestry groups, we aggregate genome-wide association study (GWAS) data from 2,535,601 individuals (39.7% non-European ancestry), including 428,452 T2D cases. We identify 1,289 independent association signals at genome-wide significance (P<5×10-8) that map to 611 loci, of which 145 loci are previously unreported. We define eight non-overlapping clusters of T2D signals characterised by distinct profiles of cardiometabolic trait associations. These clusters are differentially enriched for cell-type specific regions of open chromatin, including pancreatic islets, adipocytes, endothelial, and enteroendocrine cells. We build cluster-specific partitioned genetic risk scores (GRS) in an additional 137,559 individuals of diverse ancestry, including 10,159 T2D cases, and test their association with T2D-related vascular outcomes. Cluster-specific partitioned GRS are more strongly associated with coronary artery disease and end-stage diabetic nephropathy than an overall T2D GRS across ancestry groups, highlighting the importance of obesity-related processes in the development of vascular outcomes. Our findings demonstrate the value of integrating multi-ancestry GWAS with single-cell epigenomics to disentangle the aetiological heterogeneity driving the development and progression of T2D, which may offer a route to optimise global access to genetically-informed diabetes care.

Multi-ancestry transcriptome-wide association analyses yield insights into tobacco use biology and drug repurposing

Most transcriptome-wide association studies (TWASs) so far focus on European ancestry and lack diversity. To overcome this limitation, we aggregated genome-wide association study (GWAS) summary statistics, whole-genome sequences and expression quantitative trait locus (eQTL) data from diverse ancestries. We developed a new approach, TESLA (multi-ancestry integrative study using an optimal linear combination of association statistics), to integrate an eQTL dataset with a multi-ancestry GWAS. By exploiting shared phenotypic effects between ancestries and accommodating potential effect heterogeneities, TESLA improves power over other TWAS methods. When applied to tobacco use phenotypes, TESLA identified 273 new genes, up to 55% more compared with alternative TWAS methods. These hits and subsequent fine mapping using TESLA point to target genes with biological relevance. In silico drug-repurposing analyses highlight several drugs with known efficacy, including dextromethorphan and galantamine, and new drugs such as muscle relaxants that may be repurposed for treating nicotine addiction.

Genetic diversity fuels gene discovery for tobacco and alcohol use

Tobacco and alcohol use are heritable behaviours associated with 15% and 5.3% of worldwide deaths, respectively, due largely to broad increased risk for disease and injury1-4. These substances are used across the globe, yet genome-wide association studies have focused largely on individuals of European ancestries5. Here we leveraged global genetic diversity across 3.4 million individuals from four major clines of global ancestry (approximately 21% non-European) to power the discovery and fine-mapping of genomic loci associated with tobacco and alcohol use, to inform function of these loci via ancestry-aware transcriptome-wide association studies, and to evaluate the genetic architecture and predictive power of polygenic risk within and across populations. We found that increases in sample size and genetic diversity improved locus identification and fine-mapping resolution, and that a large majority of the 3,823 associated variants (from 2,143 loci) showed consistent effect sizes across ancestry dimensions. However, polygenic risk scores developed in one ancestry performed poorly in others, highlighting the continued need to increase sample sizes of diverse ancestries to realize any potential benefit of polygenic prediction.

Rare genetic variants explain missing heritability in smoking

Common genetic variants explain less variation in complex phenotypes than inferred from family-based studies, and there is a debate on the source of this 'missing heritability'. We investigated the contribution of rare genetic variants to tobacco use with whole-genome sequences from up to 26,257 unrelated individuals of European ancestries and 11,743 individuals of African ancestries. Across four smoking traits, single-nucleotide-polymorphism-based heritability ([Formula: see text]) was estimated from 0.13 to 0.28 (s.e., 0.10-0.13) in European ancestries, with 35-74% of it attributable to rare variants with minor allele frequencies between 0.01% and 1%. These heritability estimates are 1.5-4 times higher than past estimates based on common variants alone and accounted for 60% to 100% of our pedigree-based estimates of narrow-sense heritability ([Formula: see text], 0.18-0.34). In the African ancestry samples, [Formula: see text] was estimated from 0.03 to 0.33 (s.e., 0.09-0.14) across the four smoking traits. These results suggest that rare variants are important contributors to the heritability of smoking.

Association of Vascular Properties With the Brain White Matter Hyperintensity in Middle-Aged Population

Background The periventricular white matter is more sensitive to the systemic hemodynamic alterations than the deep white matter because of differences in its vascular structure and systemic circulation relationship. We hypothesize that periventricular white matter hyperintensity (PVWMH) volume shows greater association than deep white matter hyperintensity (DWMH) volume with vascular properties (VPs) reflecting arterial stiffness and cardiovascular remodeling, indicators of the systemic circulation. Methods and Results A total of 426 participants (age, 59.0±6.1 years; 57.5% women; and 39.7% Black race) in the Genetic Study of Atherosclerosis Risk who were aged ≥50 years and had brain magnetic resonance imaging were studied. VPs included pulse pressure, hypertensive response to exercise, diastolic brachial artery diameter, diastolic common carotid artery diameter, common carotid artery distensibility coefficient, and left ventricular function. The relative associations of VPs with PVWMH and DWMH as multiple measures within the same individual were determined using multilevel linear models. We also determined if age modified the differences in VPs associations with PVWMH and DWMH. Our findings indicated that, within the same subject, PVWMH volume had greater association than DWMH volume with pulse pressure (P=0.002), hypertensive response to exercise (P=0.04), diastolic brachial artery diameter (P=0.012), and diastolic common carotid artery diameter (P=0.04), independent of age and cardiovascular risk factors. The differences of PVWMH versus DWMH associations with VPs did not differ at any age threshold. Conclusions We show, for the first time, that PVWMH has greater association than DWMH, independent of age, with vascular measurements of arterial stiffness and cardiovascular remodeling suggesting that changes in the systemic circulation affect the PVWMH and DWMH differently.

Multi-ancestry genetic study of type 2 diabetes highlights the power of diverse populations for discovery and translation

We assembled an ancestrally diverse collection of genome-wide association studies (GWAS) of type 2 diabetes (T2D) in 180,834 affected individuals and 1,159,055 controls (48.9% non-European descent) through the Diabetes Meta-Analysis of Trans-Ethnic association studies (DIAMANTE) Consortium. Multi-ancestry GWAS meta-analysis identified 237 loci attaining stringent genome-wide significance (P < 5 × 10-9), which were delineated to 338 distinct association signals. Fine-mapping of these signals was enhanced by the increased sample size and expanded population diversity of the multi-ancestry meta-analysis, which localized 54.4% of T2D associations to a single variant with >50% posterior probability. This improved fine-mapping enabled systematic assessment of candidate causal genes and molecular mechanisms through which T2D associations are mediated, laying the foundations for functional investigations. Multi-ancestry genetic risk scores enhanced transferability of T2D prediction across diverse populations. Our study provides a step toward more effective clinical translation of T2D GWAS to improve global health for all, irrespective of genetic background.

Abstract WMP20: Association Of Periventricular And Deep White Matter Hyperintensity With Cognitive Performance In A Middle-aged Population

Introduction: White matter hyperintensity (WMH) has been associated with manual dexterity performance in apparently healthy middle-aged population. However, the relative effect of WMH topographic distribution on cognition remains unknown. We hypothesized that periventricular WMH (PVWMH) is more strongly associated with manual dexterity and cognitive measures compared to deep WMH (DWMH) given its proximity to the white matter tracts involved in motor and executive functions.

Methods: Apparently healthy middle-aged participants in the Genetic Study of Atherosclerosis Risk with brain MRI and cognitive measurements were selected. Cognitive measures included delayed word recall short (DWRS) and long (DWRL), digit symbol substitution (DSS), word fluency (WF), digit span forward test (DSFT), backward (DSBT) and total (DSTS), and grooved peg board test (GPBT) including the dominant and non-dominant hand (GPBTD, GPBTnD). PVWMH and DWMH volumes were measured using semi-automated segmentation paradigm. The relative associations of these measurements with PVWMH and DWMH volumes as multiple measurements within the same individual further nested within families were determined using multilevel linear models adjusting for intracranial volume, age, sex and cardiovascular risk factors.

Results: 782 subjects were included (age 51.1±10.6, 58.4% female and 39% African American, PVWMH volume 437.4 [0-1646] mm 3 , DWMH volume 273.5 [76-665] mm 3 ). Within the same subject, PVWMH showed greater association than DWMH with worse performance on GPBT, GPBTD, GPBTnD, DWRS, and DWRL (P=0.006, 0.01, 0.01, 0.01 and 0.03, respectively) after adjusting for age and other cardiovascular risk factors. There was no difference in the association of PVWMH versus DWMH volumes with DSS, WF, DSFT, DSBT, DSTS (P=0.19, 0.52, 0.13, 0.6, and 0.6, respectively).

Conclusions: PVWMH shows greater association than DWMH with cognitive measures of manipulative manual dexterity and recall suggesting that PVWMH and DWMH may have different cognitive phenotypes in middle-aged population with cerebral small vessel disease.

Abstract WP14: Association Of Vascular Properties With The Periventricular And Deep White Matter Hyperintensity In Middle-age Population

Introduction: The periventricular white matter is more sensitive to alterations in the hemodynamics of the systemic circulation than the deep white matter due to differences in its vascular supply and relationship to systemic circulation. We hypothesize that periventricular white matter hyperintensity volume (PVWMH) shows greater association than deep white matter hyperintensity (DWMH) volume with vascular property measures of arterial stiffness (AS) and cardiovascular remodeling, important indicators of the systemic circulation.

Methods: Participants in the Genetic Study of Atherosclerosis Risk ≥ 50 years of age with brain MRI were selected. Vascular property measures included pulse pressure (PP), hypertension response to exercise (HRE), diastolic brachial artery (BA) diameter, diastolic common carotid artery (CCA) diameter, CCA distensibility coefficient, and left ventricular function. PVWMH and DWMH were measured using semi-automated segmentation approach. The relative associations of these parameters with PVWMH and DWMH volumes as multiple measures within the same individual further nested within families were determined using multilevel linear models. We also determined if differences between PVWMH and DWMH associations were modified by age.

Results: 426 participants (age 59.0±6.1 years, 57.5% females, and 39.7% African American) were included. Within the same subject, PVWMH volume showed greater association than DWMH volume with PP (P=0.002), HRE (p=0.04), diastolic BA diameter (p=0.012), and diastolic CCA diameter (p=0.04) independent of age and cardiovascular risk factors. This difference in PVWMH versus DWMH association with vascular properties did not differ at any age threshold.

Conclusions: We show for the first time that PVWMH volume has greater association than DWMH volume, independent of age, with vascular property measurements of AS and cardiovascular remodeling. Our findings suggest that changes in the systemic circulation may affect PVWMH and DWMH differently adding to a growing body of literature that PVWMH and DWMH may represent different phenotypes of cerebral small vessel disease.

FGL1 as a modulator of plasma D-dimer levels: Exome-wide marker analysis of plasma tPA, PAI-1, and D-dimer

BACKGROUND

Use of targeted exome-arrays with common, rare variants and functionally enriched variation has led to discovery of new genes contributing to population variation in risk factors. Plasminogen activator-inhibitor 1 (PAI-1), tissue plasminogen activator (tPA), and the plasma product D-dimer are important components of the fibrinolytic system. There have been few large-scale genome-wide or exome-wide studies of PAI-1, tPA, and D-dimer.

OBJECTIVES

We sought to discover new genetic loci contributing to variation in these traits using an exome-array approach.

METHODS

Cohort-level analyses and fixed effects meta-analyses of PAI-1 (n = 15 603), tPA (n = 6876,) and D-dimer (n = 19 306) from 12 cohorts of European ancestry with diverse study design were conducted, including single-variant analyses and gene-based burden testing.

RESULTS

Five variants located in NME7, FGL1, and the fibrinogen locus, all associated with D-dimer levels, achieved genome-wide significance (P < 5 × 10-8 ). Replication was sought for these 5 variants, as well as 45 well-imputed variants with P < 1 × 10-4 in the discovery using an independent cohort. Replication was observed for three out of the five significant associations, including a novel and uncommon (0.013 allele frequency) coding variant p.Trp256Leu in FGL1 (fibrinogen-like-1) with increased plasma D-dimer levels. Additionally, a candidate-gene approach revealed a suggestive association for a coding variant (rs143202684-C) in SERPINB2, and suggestive associations with consistent effect in the replication analysis include an intronic variant (rs11057830-A) in SCARB1 associated with increased D-dimer levels.

CONCLUSION

This work provides new evidence for a role of FGL1 in hemostasis.

Gene and protein expression in human megakaryocytes derived from induced pluripotent stem cells

BACKGROUND

There is interest in deriving megakaryocytes (MKs) from pluripotent stem cells (iPSC) for biological studies. We previously found that genomic structural integrity and genotype concordance is maintained in iPSC-derived MKs.

OBJECTIVE

To establish a comprehensive dataset of genes and proteins expressed in iPSC-derived MKs.

METHODS

iPSCs were reprogrammed from peripheral blood mononuclear cells (MNCs) and MKs were derived from the iPSCs in 194 healthy European American and African American subjects. mRNA was isolated and gene expression measured by RNA sequencing. Protein expression was measured in 62 of the subjects using mass spectrometry.

RESULTS AND CONCLUSIONS

MKs expressed genes and proteins known to be important in MK and platelet function and demonstrated good agreement with previous studies in human MKs derived from CD34+ progenitor cells. The percent of cells expressing the MK markers CD41 and CD42a was consistent in biological replicates, but variable across subjects, suggesting that unidentified subject-specific factors determine differentiation of MKs from iPSCs. Gene and protein sets important in platelet function were associated with increasing expression of CD41/42a, while those related to more basic cellular functions were associated with lower CD41/42a expression. There was differential gene expression by the sex and race (but not age) of the subject. Numerous genes and proteins were highly expressed in MKs but not known to play a role in MK or platelet function; these represent excellent candidates for future study of hematopoiesis, platelet formation, and/or platelet function.

Cerebral Small-Vessel Disease in Individuals with a Family History of Coronary Heart Disease: The Atherosclerosis Risk in Communities Study

INTRODUCTION

The degree to which a family history of coronary heart disease (FHCHD) is associated with silent cerebral small-vessel disease (cSVD) among healthy adults, independent of prevalent CHD and traditional risk factors, is unknown.

METHODS

The Atherosclerosis Risk in Communities (ARIC) study is a community-based cohort study with self-reported family history data and brain magnetic resonance imaging (ages 68-88). The association between markers of cSVD (lacunar infarcts and cerebral microbleeds), or log-transformed white matter hyperintensity (WMH) volume, and FHCHD, or the number of affected relatives was examined using separate adjusted logistic or linear regression models, respectively. Race interaction terms were evaluated.

RESULTS

Of 1,639 participants without prevalent CHD (76 ± 5 years, 62% female, 29% black), 686 (42%) had FHCHD. There were higher odds of lacunar infarct (OR 1.40, 95% CI 1.07-1.84) among those with parental FHCHD and higher odds of microhemorrhages (lobar OR 1.86, 95% CI 1.13-3.06; subcortical OR 1.47, 95% CI 1.01-2.15) among those with sibling FHCHD. A greater number of any relative affected was associated with higher odds of lacunar infarct (OR 1.24, 95% CI 1.04-1.47) and lobar microhemorrhages (OR 1.31, 95% CI 1.05-1.64) but not subcortical microhemorrhages (OR 1.09, 95% CI 0.92-1.28). Odds of having a lacunar infarct were higher among blacks (p-interaction 0.04) with paternal FHCHD (OR 2.20, CI 1.35-3.58) than whites with paternal FHCHD (OR 1.17, CI 0.87-1.56). There was no association with WMH.

DISCUSSION/CONCLUSION

Markers of cSVD, specifically lacunar infarcts and microhemorrhages, appear to be associated with FHCHD, potentially representing shared mechanisms in different vascular beds, and perhaps a genetic propensity for vascular disease.

The trans-ancestral genomic architecture of glycemic traits

Glycemic traits are used to diagnose and monitor type 2 diabetes and cardiometabolic health. To date, most genetic studies of glycemic traits have focused on individuals of European ancestry. Here we aggregated genome-wide association studies comprising up to 281,416 individuals without diabetes (30% non-European ancestry) for whom fasting glucose, 2-h glucose after an oral glucose challenge, glycated hemoglobin and fasting insulin data were available. Trans-ancestry and single-ancestry meta-analyses identified 242 loci (99 novel; P < 5 × 10-8), 80% of which had no significant evidence of between-ancestry heterogeneity. Analyses restricted to individuals of European ancestry with equivalent sample size would have led to 24 fewer new loci. Compared with single-ancestry analyses, equivalent-sized trans-ancestry fine-mapping reduced the number of estimated variants in 99% credible sets by a median of 37.5%. Genomic-feature, gene-expression and gene-set analyses revealed distinct biological signatures for each trait, highlighting different underlying biological pathways. Our results increase our understanding of diabetes pathophysiology by using trans-ancestry studies for improved power and resolution.

Association of Coronary Artery Atherosclerosis With Brain White Matter Hyperintensity

[Figure: see text].

Abstract P65: Arterial Stiffness Relationship With Periventricular and Deep White Matter Hyperintensities in a Middle-Aged Population

Introduction and hypothesis: Periventricular white matter (PVWM) is perfused by perforators originating from the proximal large caliber cerebral arteries with direct transmission of systemic hemodynamic changes as compared to the deep white matter (DWM) that is perfused by distal perforators with dampened transmission. We hypothesized that PVWM hyperintensities (PVWMH) are more strongly associated with the systemic vascular properties representing arterial stiffness (AS) due to vascular aging and hypertension compared to DWM hyperintensities (DWMH).

Methods: Apparently healthy first degree relatives of subjects with coronary artery disease who had brain MRI as participants in GeneSTAR were selected. PVWMH and DWMH volumes were computed by an automated analysis software. Arterial stiffness indicators included peak systolic blood pressure during treadmill test, pulse pressure (PP), brachial artery (BA) diameter, carotid distensibility index, and diastolic dysfunction estimated by filling and systolic emptying rates (SER) on sestamibi cardiac scan. Statistical analysis was completed using unadjusted and age-adjusted linear regression models in the cohort that was divided into 4 quartiles based on age distribution (1 st : 29-43, 2 nd : 44-51, 3 rd : 52-59, 4 th : 60-74).

Results: 782 participants were included (age 51.05±10.6 years, 58.4% female, 39% African American, PVWMH volume 1977.01±5195.81 ml and DWMH volume 625.93±1108.37 ml). In the unadjusted analysis, PVWMH showed association with PP in the 1 st and 3 rd age quartiles (p= 0.02 and 0.03, respectively) and SER in the 2 nd and 4 th age quartiles (p=0.04 and 0.03). This relationship remained significant in the adjusted analysis for SER (p=0.047, and 0.032, respectively) and in the 1 st age quartile for PP (P=0.02). DWMH relationship with BA diameter was significant in the 2 nd age quartile in the unadjusted and adjusted models (P=0.011, and 0.02, respectively). DWMH was also associated with ratio of cardiac filling/SER (P=0.039 and 0.047) in the 1 st age quartile.

Conclusions: In this exploratory analysis, AS measures show different relationships to PVWMH and DWMH across age groups in middle-aged population suggesting a potential mechanistic role of AS in WMH pathology.

Multi-Ancestry Genome-wide Association Study Accounting for Gene-Psychosocial Factor Interactions Identifies Novel Loci for Blood Pressure Traits

Psychological and social factors are known to influence blood pressure (BP) and risk of hypertension and associated cardiovascular diseases. To identify novel BP loci, we carried out genome-wide association meta-analyses of systolic, diastolic, pulse, and mean arterial BP taking into account the interaction effects of genetic variants with three psychosocial factors: depressive symptoms, anxiety symptoms, and social support. Analyses were performed using a two-stage design in a sample of up to 128,894 adults from 5 ancestry groups. In the combined meta-analyses of Stages 1 and 2, we identified 59 loci (p value <5e-8), including nine novel BP loci. The novel associations were observed mostly with pulse pressure, with fewer observed with mean arterial pressure. Five novel loci were identified in African ancestry, and all but one showed patterns of interaction with at least one psychosocial factor. Functional annotation of the novel loci supports a major role for genes implicated in the immune response (PLCL2), synaptic function and neurotransmission (LIN7A, PFIA2), as well as genes previously implicated in neuropsychiatric or stress-related disorders (FSTL5, CHODL). These findings underscore the importance of considering psychological and social factors in gene discovery for BP, especially in non-European populations.

White Matter Injury Is Associated with Reduced Manual Dexterity and Elevated Serum Ceramides in Subjects with Cerebral Small Vessel Disease

INTRODUCTION

We have demonstrated that asymptomatic cerebral small vessel disease (cSVD) measured by white matter hyperintensity volume is associated with reduced manipulative manual dexterity on the Grooved Peg Board Test (GPBT) in middle-aged healthy individuals with a family history of early coronary artery disease. In this current study, we aim to identify the association of subcortical white matter microstructural impairment measured by diffusion tensor imaging, manual dexterity measured by GPBT and circulating serums ceramide, another marker for white matter injury. We hypothesize that lower regional fractional anisotropy (rFA) is associated with worse performance on GPBT and elevated serum ceramides in the same study population.

METHODS

rFA of 48 regions representing the subcortical white matters were analyzed in GeneSTAR participants in addition to serum ceramides and GPBT scores. Unadjusted univariable analyses with Bonferroni correction for multiple comparisons were completed using Spearman correlation for testing the associations between ceramides, rFA of subcortical white matter, and GPBT performance. Subsequently, sensitivity analyses were performed after excluding the participants that had any physical limitation that may influence their performance on GPBT. Finally, in the adjusted analysis using generalized estimating equation, linear regression models were performed for the areas that met significance threshold in the unadjusted analyses.

RESULTS

112 subjects (age [49 ± 11], 51% female, 39.3% African American) were included. Adjusted analyses for the significant correlations that met the Bonferroni correction threshold in the unadjusted univariable analyses identified significant negative associations between rFA of the right fornix (RF) and log-GPBT score (β = -0.497, p = 0.037). In addition, rFA of RF negatively correlated with log serum ceramide levels (C18: β = -0.03, p = 0.003, C20: β = -0.0002, p = 0.004) and rFA of left genu of corpus callosum negatively correlated with log C18 level (β = -0.0103, p = 0.027).

CONCLUSIONS

These results demonstrate that subcortical microstructural white matter disruption is associated with elevated serum ceramides and reduced manual dexterity in a population with cSVD. These findings suggest that injury to white matter tracts undermines neural networks, with functional consequences in a middle-aged population with cardiovascular risk factors.

Genome-Wide Association Transethnic Meta-Analyses Identifies Novel Associations Regulating Coagulation Factor VIII and von Willebrand Factor Plasma Levels

BACKGROUND

Factor VIII (FVIII) and its carrier protein von Willebrand factor (VWF) are associated with risk of arterial and venous thrombosis and with hemorrhagic disorders. We aimed to identify and functionally test novel genetic associations regulating plasma FVIII and VWF.

METHODS

We meta-analyzed genome-wide association results from 46 354 individuals of European, African, East Asian, and Hispanic ancestry. All studies performed linear regression analysis using an additive genetic model and associated ≈35 million imputed variants with natural log-transformed phenotype levels. In vitro gene silencing in cultured endothelial cells was performed for candidate genes to provide additional evidence on association and function. Two-sample Mendelian randomization analyses were applied to test the causal role of FVIII and VWF plasma levels on the risk of arterial and venous thrombotic events.

RESULTS

We identified 13 novel genome-wide significant ( P≤2.5×10^-8) associations, 7 with FVIII levels ( FCHO2/TMEM171/TNPO1, HLA, SOX17/RP1, LINC00583/NFIB, RAB5C-KAT2A, RPL3/TAB1/SYNGR1, and ARSA) and 11 with VWF levels ( PDHB/PXK/KCTD6, SLC39A8, FCHO2/TMEM171/TNPO1, HLA, GIMAP7/GIMAP4, OR13C5/NIPSNAP, DAB2IP, C2CD4B, RAB5C-KAT2A, TAB1/SYNGR1, and ARSA), beyond 10 previously reported associations with these phenotypes. Functional validation provided further evidence of association for all loci on VWF except ARSA and DAB2IP. Mendelian randomization suggested causal effects of plasma FVIII activity levels on venous thrombosis and coronary artery disease risk and plasma VWF levels on ischemic stroke risk.

CONCLUSIONS

The meta-analysis identified 13 novel genetic loci regulating FVIII and VWF plasma levels, 10 of which we validated functionally. We provide some evidence for a causal role of these proteins in thrombotic events.

Multi-ancestry sleep-by-SNP interaction analysis in 126,926 individuals reveals lipid loci stratified by sleep duration

Both short and long sleep are associated with an adverse lipid profile, likely through different biological pathways. To elucidate the biology of sleep-associated adverse lipid profile, we conduct multi-ancestry genome-wide sleep-SNP interaction analyses on three lipid traits (HDL-c, LDL-c and triglycerides). In the total study sample (discovery + replication) of 126,926 individuals from 5 different ancestry groups, when considering either long or short total sleep time interactions in joint analyses, we identify 49 previously unreported lipid loci, and 10 additional previously unreported lipid loci in a restricted sample of European-ancestry cohorts. In addition, we identify new gene-sleep interactions for known lipid loci such as LPL and PCSK9. The previously unreported lipid loci have a modest explained variance in lipid levels: most notable, gene-short-sleep interactions explain 4.25% of the variance in triglyceride level. Collectively, these findings contribute to our understanding of the biological mechanisms involved in sleep-associated adverse lipid profiles.

Sleep duration is associated with an adverse lipid profile. Here, the authors perform genome-wide gene-by-sleep interaction analysis and find 49 previously unreported lipid loci when considering short or long total sleep time.

Associations of autozygosity with a broad range of human phenotypes

In many species, the offspring of related parents suffer reduced reproductive success, a phenomenon known as inbreeding depression. In humans, the importance of this effect has remained unclear, partly because reproduction between close relatives is both rare and frequently associated with confounding social factors. Here, using genomic inbreeding coefficients (FROH) for >1.4 million individuals, we show that FROH is significantly associated (p < 0.0005) with apparently deleterious changes in 32 out of 100 traits analysed. These changes are associated with runs of homozygosity (ROH), but not with common variant homozygosity, suggesting that genetic variants associated with inbreeding depression are predominantly rare. The effect on fertility is striking: FROH equivalent to the offspring of first cousins is associated with a 55% decrease [95% CI 44-66%] in the odds of having children. Finally, the effects of FROH are confirmed within full-sibling pairs, where the variation in FROH is independent of all environmental confounding.

Evolution of Hominin Polyunsaturated Fatty Acid Metabolism: From Africa to the New World

The metabolic conversion of dietary omega-3 and omega-6 18 carbon (18C) to long chain (>20 carbon) polyunsaturated fatty acids (LC-PUFAs) is vital for human life. The rate-limiting steps of this process are catalyzed by fatty acid desaturase (FADS) 1 and 2. Therefore, understanding the evolutionary history of the FADS genes is essential to our understanding of hominin evolution. The FADS genes have two haplogroups, ancestral and derived, with the derived haplogroup being associated with more efficient LC-PUFA biosynthesis than the ancestral haplogroup. In addition, there is a complex global distribution of these haplogroups that is suggestive of Neanderthal introgression. We confirm that Native American ancestry is nearly fixed for the ancestral haplogroup, and replicate a positive selection signal in Native Americans. This positive selection potentially continued after the founding of the Americas, although simulations suggest that the timing is dependent on the allele frequency of the ancestral Beringian population. We also find that the Neanderthal FADS haplotype is more closely related to the derived haplogroup and the Denisovan clusters closer to the ancestral haplogroup. Furthermore, the derived haplogroup has a time to the most recent common ancestor of 688,474 years before present. These results support an ancient polymorphism, as opposed to Neanderthal introgression, forming in the FADS region during the Pleistocene with possibly differential selection pressures on both haplogroups. The near fixation of the ancestral haplogroup in Native American ancestry calls for future studies to explore the potential health risk of associated low LC-PUFA levels in these populations.

Exome Chip Analysis Identifies Low-Frequency and Rare Variants in MRPL38 for White Matter Hyperintensities on Brain Magnetic Resonance Imaging

Background and Purpose- White matter hyperintensities (WMH) on brain magnetic resonance imaging are typical signs of cerebral small vessel disease and may indicate various preclinical, age-related neurological disorders, such as stroke. Though WMH are highly heritable, known common variants explain a small proportion of the WMH variance. The contribution of low-frequency/rare coding variants to WMH burden has not been explored. Methods- In the discovery sample we recruited 20 719 stroke/dementia-free adults from 13 population-based cohort studies within the Cohorts for Heart and Aging Research in Genomic Epidemiology consortium, among which 17 790 were of European ancestry and 2929 of African ancestry. We genotyped these participants at ≈250 000 mostly exonic variants with Illumina HumanExome BeadChip arrays. We performed ethnicity-specific linear regression on rank-normalized WMH in each study separately, which were then combined in meta-analyses to test for association with single variants and genes aggregating the effects of putatively functional low-frequency/rare variants. We then sought replication of the top findings in 1192 adults (European ancestry) with whole exome/genome sequencing data from 2 independent studies. Results- At 17q25, we confirmed the association of multiple common variants in TRIM65, FBF1, and ACOX1 ( P<6×10^-7). We also identified a novel association with 2 low-frequency nonsynonymous variants in MRPL38 (lead, rs34136221; P_EA=4.5×10^-8) partially independent of known common signal ( P_{EA(conditional)}=1.4×10^-3). We further identified a locus at 2q33 containing common variants in NBEAL1, CARF, and WDR12 (lead, rs2351524; P_all=1.9×10^-10). Although our novel findings were not replicated because of limited power and possible differences in study design, meta-analysis of the discovery and replication samples yielded stronger association for the 2 low-frequency MRPL38 variants ( P_rs34136221=2.8×10^-8). Conclusions- Both common and low-frequency/rare functional variants influence WMH. Larger replication and experimental follow-up are essential to confirm our findings and uncover the biological causal mechanisms of age-related WMH.

A genome-wide association study identifies genetic loci associated with specific lobar brain volumes

Brain lobar volumes are heritable but genetic studies are limited. We performed genome-wide association studies of frontal, occipital, parietal and temporal lobe volumes in 16,016 individuals, and replicated our findings in 8,789 individuals. We identified six genetic loci associated with specific lobar volumes independent of intracranial volume. Two loci, associated with occipital (6q22.32) and temporal lobe volume (12q14.3), were previously reported to associate with intracranial and hippocampal volume, respectively. We identified four loci previously unknown to affect brain volumes: 3q24 for parietal lobe volume, and 1q22, 4p16.3 and 14q23.1 for occipital lobe volume. The associated variants were located in regions enriched for histone modifications (DAAM1 and THBS3), or close to genes causing Mendelian brain-related diseases (ZIC4 and FGFRL1). No genetic overlap between lobar volumes and neurological or psychiatric diseases was observed. Our findings reveal part of the complex genetics underlying brain development and suggest a role for regulatory regions in determining brain volumes.

A multi-ancestry genome-wide study incorporating gene-smoking interactions identifies multiple new loci for pulse pressure and mean arterial pressure

Elevated blood pressure (BP), a leading cause of global morbidity and mortality, is influenced by both genetic and lifestyle factors. Cigarette smoking is one such lifestyle factor. Across five ancestries, we performed a genome-wide gene-smoking interaction study of mean arterial pressure (MAP) and pulse pressure (PP) in 129 913 individuals in stage 1 and follow-up analysis in 480 178 additional individuals in stage 2. We report here 136 loci significantly associated with MAP and/or PP. Of these, 61 were previously published through main-effect analysis of BP traits, 37 were recently reported by us for systolic BP and/or diastolic BP through gene-smoking interaction analysis and 38 were newly identified (P < 5 × 10-8, false discovery rate < 0.05). We also identified nine new signals near known loci. Of the 136 loci, 8 showed significant interaction with smoking status. They include CSMD1 previously reported for insulin resistance and BP in the spontaneously hypertensive rats. Many of the 38 new loci show biologic plausibility for a role in BP regulation. SLC26A7 encodes a chloride/bicarbonate exchanger expressed in the renal outer medullary collecting duct. AVPR1A is widely expressed, including in vascular smooth muscle cells, kidney, myocardium and brain. FHAD1 is a long non-coding RNA overexpressed in heart failure. TMEM51 was associated with contractile function in cardiomyocytes. CASP9 plays a central role in cardiomyocyte apoptosis. Identified only in African ancestry were 30 novel loci. Our findings highlight the value of multi-ancestry investigations, particularly in studies of interaction with lifestyle factors, where genomic and lifestyle differences may contribute to novel findings.

Genome-Wide Association Study of Heavy Smoking and Daily/Nondaily Smoking in the Hispanic Community Health Study/Study of Latinos (HCHS/SOL)

Introduction

Genetic variants associated with nicotine dependence have previously been identified, primarily in European-ancestry populations. No genome-wide association studies (GWAS) have been reported for smoking behaviors in Hispanics/Latinos in the United States and Latin America, who are of mixed ancestry with European, African, and American Indigenous components.

Methods

We examined genetic associations with smoking behaviors in the Hispanic Community Health Study/Study of Latinos (HCHS/SOL) (N = 12 741 with smoking data, 5119 ever-smokers), using ~2.3 million genotyped variants imputed to the 1000 Genomes Project phase 3. Mixed logistic regression models accounted for population structure, sampling, relatedness, sex, and age.

Results

The known region of CHRNA5, which encodes the α5 cholinergic nicotinic receptor subunit, was associated with heavy smoking at genome-wide significance (p ≤ 5 × 10-8) in a comparison of 1929 ever-smokers reporting cigarettes per day (CPD) > 10 versus 3156 reporting CPD ≤ 10. The functional variant rs16969968 in CHRNA5 had a p value of 2.20 × 10-7 and odds ratio (OR) of 1.32 for the minor allele (A); its minor allele frequency was 0.22 overall and similar across Hispanic/Latino background groups (Central American = 0.17; South American = 0.19; Mexican = 0.18; Puerto Rican = 0.22; Cuban = 0.29; Dominican = 0.19). CHRNA4 on chromosome 20 attained p < 10-4, supporting prior findings in non-Hispanics. For nondaily smoking, which is prevalent in Hispanic/Latino smokers, compared to daily smoking, loci on chromosomes 2 and 4 achieved genome-wide significance; replication attempts were limited by small Hispanic/Latino sample sizes.

Conclusions

Associations of nicotinic receptor gene variants with smoking, first reported in non-Hispanic European-ancestry populations, generalized to Hispanics/Latinos despite different patterns of smoking behavior.

Implications

We conducted the first large-scale genome-wide association study (GWAS) of smoking behavior in a US Hispanic/Latino cohort, and the first GWAS of daily/nondaily smoking in any population. Results show that the region of the nicotinic receptor subunit gene CHRNA5, which in non-Hispanic European-ancestry smokers has been associated with heavy smoking as well as cessation and treatment efficacy, is also significantly associated with heavy smoking in this Hispanic/Latino cohort. The results are an important addition to understanding the impact of genetic variants in understudied Hispanic/Latino smokers.

Multiancestry Genome-Wide Association Study of Lipid Levels Incorporating Gene-Alcohol Interactions

A person's lipid profile is influenced by genetic variants and alcohol consumption, but the contribution of interactions between these exposures has not been studied. We therefore incorporated gene-alcohol interactions into a multiancestry genome-wide association study of levels of high-density lipoprotein cholesterol, low-density lipoprotein cholesterol, and triglycerides. We included 45 studies in stage 1 (genome-wide discovery) and 66 studies in stage 2 (focused follow-up), for a total of 394,584 individuals from 5 ancestry groups. Analyses covered the period July 2014-November 2017. Genetic main effects and interaction effects were jointly assessed by means of a 2-degrees-of-freedom (df) test, and a 1-df test was used to assess the interaction effects alone. Variants at 495 loci were at least suggestively associated (P < 1 × 10-6) with lipid levels in stage 1 and were evaluated in stage 2, followed by combined analyses of stage 1 and stage 2. In the combined analysis of stages 1 and 2, a total of 147 independent loci were associated with lipid levels at P < 5 × 10-8 using 2-df tests, of which 18 were novel. No genome-wide-significant associations were found testing the interaction effect alone. The novel loci included several genes (proprotein convertase subtilisin/kexin type 5 (PCSK5), vascular endothelial growth factor B (VEGFB), and apolipoprotein B mRNA editing enzyme, catalytic polypeptide 1 (APOBEC1) complementation factor (A1CF)) that have a putative role in lipid metabolism on the basis of existing evidence from cellular and experimental models.

Multi-ancestry genome-wide gene-smoking interaction study of 387,272 individuals identifies new loci associated with serum lipids

The concentrations of high- and low-density-lipoprotein cholesterol and triglycerides are influenced by smoking, but it is unknown whether genetic associations with lipids may be modified by smoking. We conducted a multi-ancestry genome-wide gene-smoking interaction study in 133,805 individuals with follow-up in an additional 253,467 individuals. Combined meta-analyses identified 13 new loci associated with lipids, some of which were detected only because association differed by smoking status. Additionally, we demonstrate the importance of including diverse populations, particularly in studies of interactions with lifestyle factors, where genomic and lifestyle differences by ancestry may contribute to novel findings.

Associations of Mitochondrial and Nuclear Mitochondrial Variants and Genes with Seven Metabolic Traits

Mitochondria (MT), the major site of cellular energy production, are under dual genetic control by 37 mitochondrial DNA (mtDNA) genes and numerous nuclear genes (MT-nDNA). In the CHARGEmtDNA+ Consortium, we studied genetic associations of mtDNA and MT-nDNA associations with body mass index (BMI), waist-hip-ratio (WHR), glucose, insulin, HOMA-B, HOMA-IR, and HbA1c. This 45-cohort collaboration comprised 70,775 (insulin) to 170,202 (BMI) pan-ancestry individuals. Validation and imputation of mtDNA variants was followed by single-variant and gene-based association testing. We report two significant common variants, one in MT-ATP6 associated (p ≤ 5E-04) with WHR and one in the D-loop with glucose. Five rare variants in MT-ATP6, MT-ND5, and MT-ND6 associated with BMI, WHR, or insulin. Gene-based meta-analysis identified MT-ND3 associated with BMI (p ≤ 1E-03). We considered 2,282 MT-nDNA candidate gene associations compiled from online summary results for our traits (20 unique studies with 31 dataset consortia's genome-wide associations [GWASs]). Of these, 109 genes associated (p ≤ 1E-06) with at least 1 of our 7 traits. We assessed regulatory features of variants in the 109 genes, cis- and trans-gene expression regulation, and performed enrichment and protein-protein interactions analyses. Of the identified mtDNA and MT-nDNA genes, 79 associated with adipose measures, 49 with glucose/insulin, 13 with risk for type 2 diabetes, and 18 with cardiovascular disease, indicating for pleiotropic effects with health implications. Additionally, 21 genes related to cholesterol, suggesting additional important roles for the genes identified. Our results suggest that mtDNA and MT-nDNA genes and variants reported make important contributions to glucose and insulin metabolism, adipocyte regulation, diabetes, and cardiovascular disease.

Mendelian randomization evaluation of causal effects of fibrinogen on incident coronary heart disease

BACKGROUND

Fibrinogen is an essential hemostatic factor and cardiovascular disease risk factor. Early attempts at evaluating the causal effect of fibrinogen on coronary heart disease (CHD) and myocardial infraction (MI) using Mendelian randomization (MR) used single variant approaches, and did not take advantage of recent genome-wide association studies (GWAS) or multi-variant, pleiotropy robust MR methodologies.

METHODS AND FINDINGS

We evaluated evidence for a causal effect of fibrinogen on both CHD and MI using MR. We used both an allele score approach and pleiotropy robust MR models. The allele score was composed of 38 fibrinogen-associated variants from recent GWAS. Initial analyses using the allele score used a meta-analysis of 11 European-ancestry prospective cohorts, free of CHD and MI at baseline, to examine incidence CHD and MI. We also applied 2 sample MR methods with data from a prevalent CHD and MI GWAS. Results are given in terms of the hazard ratio (HR) or odds ratio (OR), depending on the study design, and associated 95% confidence interval (CI). In single variant analyses no causal effect of fibrinogen on CHD or MI was observed. In multi-variant analyses using incidence CHD cases and the allele score approach, the estimated causal effect (HR) of a 1 g/L higher fibrinogen concentration was 1.62 (CI = 1.12, 2.36) when using incident cases and the allele score approach. In 2 sample MR analyses that accounted for pleiotropy, the causal estimate (OR) was reduced to 1.18 (CI = 0.98, 1.42) and 1.09 (CI = 0.89, 1.33) in the 2 most precise (smallest CI) models, out of 4 models evaluated. In the 2 sample MR analyses for MI, there was only very weak evidence of a causal effect in only 1 out of 4 models.

CONCLUSIONS

A small causal effect of fibrinogen on CHD is observed using multi-variant MR approaches which account for pleiotropy, but not single variant MR approaches. Taken together, results indicate that even with large sample sizes and multi-variant approaches MR analyses still cannot exclude the null when estimating the causal effect of fibrinogen on CHD, but that any potential causal effect is likely to be much smaller than observed in epidemiological studies.

Diabetes and Platelet Response to Low-Dose Aspirin

CONTEXT

Previous studies have suggested less cardioprotective benefit of aspirin in adults with diabetes, raising concerns about "aspirin resistance" and potentially reduced effectiveness for prevention of cardiovascular disease (CVD).

OBJECTIVE

To examine differences in platelet response to aspirin by diabetes status.

DESIGN, SETTING, PARTICIPANTS

We examined platelet response before and after aspirin (81 mg/day for 14 days) in 2113 adults (175 with diabetes, 1,938 without diabetes), in the Genetic Study of Aspirin Responsiveness cohort, who had family history of early-onset CVD.

MAIN OUTCOME MEASURES

In vivo platelet activation (urinary thromboxane B2), in vitro platelet aggregation to agonists (arachidonic acid, adenosine diphosphate, collagen), and platelet function analyzer-100 closure time.

RESULTS

Although adults with diabetes had higher in vivo platelet activation before aspirin, the reduction in in vivo platelet activation after aspirin was similar in those with vs without diabetes. Likewise, the reduction in multiple in vitro platelet measures was similar after aspirin by diabetes status. In regression analyses adjusted for age, sex, race, BMI, smoking, platelet counts, and fibrinogen levels, in vivo platelet activation remained higher in adults with vs without diabetes after aspirin (P = 0.04), but this difference was attenuated after additional adjustment for preaspirin levels (P = 0.10). No differences by diabetes status were noted for any of the in vitro platelet measures after aspirin in fully adjusted models that also accounted for preaspirin levels.

CONCLUSIONS

In vitro platelet response to aspirin does not differ by diabetes status, suggesting no intrinsic differences in platelet response to aspirin. Instead, factors extrinsic to platelet function should be investigated to give further insights into aspirin use for primary prevention in diabetes.

Genome Analyses of >200,000 Individuals Identify 58 Loci for Chronic Inflammation and Highlight Pathways that Link Inflammation and Complex Disorders

C-reactive protein (CRP) is a sensitive biomarker of chronic low-grade inflammation and is associated with multiple complex diseases. The genetic determinants of chronic inflammation remain largely unknown, and the causal role of CRP in several clinical outcomes is debated. We performed two genome-wide association studies (GWASs), on HapMap and 1000 Genomes imputed data, of circulating amounts of CRP by using data from 88 studies comprising 204,402 European individuals. Additionally, we performed in silico functional analyses and Mendelian randomization analyses with several clinical outcomes. The GWAS meta-analyses of CRP revealed 58 distinct genetic loci (p < 5 × 10-8). After adjustment for body mass index in the regression analysis, the associations at all except three loci remained. The lead variants at the distinct loci explained up to 7.0% of the variance in circulating amounts of CRP. We identified 66 gene sets that were organized in two substantially correlated clusters, one mainly composed of immune pathways and the other characterized by metabolic pathways in the liver. Mendelian randomization analyses revealed a causal protective effect of CRP on schizophrenia and a risk-increasing effect on bipolar disorder. Our findings provide further insights into the biology of inflammation and could lead to interventions for treating inflammation and its clinical consequences.

Novel genetic associations for blood pressure identified via gene-alcohol interaction in up to 570K individuals across multiple ancestries

Heavy alcohol consumption is an established risk factor for hypertension; the mechanism by which alcohol consumption impact blood pressure (BP) regulation remains unknown. We hypothesized that a genome-wide association study accounting for gene-alcohol consumption interaction for BP might identify additional BP loci and contribute to the understanding of alcohol-related BP regulation. We conducted a large two-stage investigation incorporating joint testing of main genetic effects and single nucleotide variant (SNV)-alcohol consumption interactions. In Stage 1, genome-wide discovery meta-analyses in ≈131K individuals across several ancestry groups yielded 3,514 SNVs (245 loci) with suggestive evidence of association (P < 1.0 x 10-5). In Stage 2, these SNVs were tested for independent external replication in ≈440K individuals across multiple ancestries. We identified and replicated (at Bonferroni correction threshold) five novel BP loci (380 SNVs in 21 genes) and 49 previously reported BP loci (2,159 SNVs in 109 genes) in European ancestry, and in multi-ancestry meta-analyses (P < 5.0 x 10-8). For African ancestry samples, we detected 18 potentially novel BP loci (P < 5.0 x 10-8) in Stage 1 that warrant further replication. Additionally, correlated meta-analysis identified eight novel BP loci (11 genes). Several genes in these loci (e.g., PINX1, GATA4, BLK, FTO and GABBR2) have been previously reported to be associated with alcohol consumption. These findings provide insights into the role of alcohol consumption in the genetic architecture of hypertension.

The relationship of family history and risk of type 2 diabetes differs by ancestry

AIM

Type 2 diabetes (T2DM) in a first-degree relative is a risk factor for incident diabetes. Americans of African ancestry (AA) have higher rates of T2DM than Americans of European ancestry (EA). Thus, we aimed to determine whether the presence, number and kinship of affected relatives are associated with race-specific T2DM incidence in a prospective study of participants from the Genetic Study of Atherosclerosis Risk (GeneSTAR), who underwent baseline screening including a detailed family history.

METHODS

Nondiabetic healthy siblings (n=1405) of patients with early-onset coronary artery disease (18-59 years) were enrolled (861 EA and 544 AA) and followed for incident T2DM (mean 14±6 years).

RESULTS

Baseline age was 46.2±7.3 years and 56% were female. T2DM occurred in 12.3% of EA and 19.1% of AA. Among EA, 32.6% had ≥1 affected first-degree relatives versus 53.1% in AA, P<0.0001. In fully adjusted Cox proportional hazard analyses, any family history was related to incident T2DM in EA (HR=2.53, 95% CI: 1.58-4.06) but not in AA (HR=1.01, 0.67-1.53). The number of affected relatives conferred incremental risk of T2DM in EA with HR=1.82 (1.08-3.06), 4.83 (2.15-10.85) and 8.46 (3.09-23.91) for 1, 2, and ≥3 affected, respectively. In AA only ≥3 affected increased risk (HR=2.45, 1.44-4.19). Specific kinship patterns were associated with incident T2DM in EA but not in AA.

CONCLUSIONS

The presence of any first-degree relative with T2DM does not discriminate risk in AA given the high race-specific prevalence of diabetes. Accounting for the number of affected relatives may more appropriately estimate risk for incident diabetes in both races.

Targeted deep sequencing of the PEAR1 locus for platelet aggregation in European and African American families

Coronary artery disease (CAD) remains a major cause of mortality and morbidity worldwide. The aggregation of activated platelets on a ruptured atherosclerotic plaque is a critical step in most acute cardiovascular events like myocardial infarction. Platelet aggregation both at baseline and after aspirin is highly heritable. Genome-wide association studies (GWAS) have identified a common variant within the first intron of the platelet endothelial aggregation receptor1 (PEAR1), to be robustly associated with platelet aggregation. In this study, we used targeted deep sequencing to fine-map the prior GWAS peak and identify additional rare variants of PEAR1 that account for missing heritability in platelet aggregation within the GeneSTAR families. In this study, 1709 subjects (1043 European Americans, EA and 666 African Americans, AA) from families in the GeneSTAR study were included. In vitro platelet aggregation in response to collagen, ADP and epinephrine was measured at baseline and 14 days after aspirin therapy (81 mg/day). Targeted deep sequencing of PEAR1 in addition to 2kb of upstream and downstream of the gene was performed. Under an additive genetic model, the association of single variants of PEAR1 with platelet aggregation phenotypes were examined. Additionally, we examined the association between the burden of PEAR1 rare non-synonymous variants and platelet aggregation phenotypes. Of 532 variants identified through sequencing, the intron 1 variant, rs12041331, was significantly associated with all platelet aggregation phenotypes at baseline and after platelet inhibition with aspirin therapy. rs12566888, which is in linkage disequilibrium with rs12041331, was associated with platelet aggregation phenotypes but to a lesser extent. In the EA families, the burden of PEAR1 missense variants was associated with platelet aggregation after aspirin therapy when the platelets were stimulated with epinephrine (p = 0.0009) and collagen (p = 0.03). In AAs, the burden of PEAR1 missense variants was associated, to a lesser degree, with platelet aggregation in response to epinephrine (p = 0.02) and ADP (p = 0.04). Our study confirmed that the GWAS-identified variant, rs12041331, is the strongest variant associated with platelet aggregation both at baseline and after aspirin therapy in our GeneSTAR families in both races. We identified additional association of rare missense variants in PEAR1 with platelet aggregation following aspirin therapy. However, we observed a racial difference in the contribution of these rare variants to the platelet aggregation, most likely due to higher residual missing heritability of platelet aggregation after accounting for rs12041331 in the EAs compared to AAs.

A Large-Scale Multi-ancestry Genome-wide Study Accounting for Smoking Behavior Identifies Multiple Significant Loci for Blood Pressure

Genome-wide association analysis advanced understanding of blood pressure (BP), a major risk factor for vascular conditions such as coronary heart disease and stroke. Accounting for smoking behavior may help identify BP loci and extend our knowledge of its genetic architecture. We performed genome-wide association meta-analyses of systolic and diastolic BP incorporating gene-smoking interactions in 610,091 individuals. Stage 1 analysis examined ∼18.8 million SNPs and small insertion/deletion variants in 129,913 individuals from four ancestries (European, African, Asian, and Hispanic) with follow-up analysis of promising variants in 480,178 additional individuals from five ancestries. We identified 15 loci that were genome-wide significant (p < 5 × 10-8) in stage 1 and formally replicated in stage 2. A combined stage 1 and 2 meta-analysis identified 66 additional genome-wide significant loci (13, 35, and 18 loci in European, African, and trans-ancestry, respectively). A total of 56 known BP loci were also identified by our results (p < 5 × 10-8). Of the newly identified loci, ten showed significant interaction with smoking status, but none of them were replicated in stage 2. Several loci were identified in African ancestry, highlighting the importance of genetic studies in diverse populations. The identified loci show strong evidence for regulatory features and support shared pathophysiology with cardiometabolic and addiction traits. They also highlight a role in BP regulation for biological candidates such as modulators of vascular structure and function (CDKN1B, BCAR1-CFDP1, PXDN, EEA1), ciliopathies (SDCCAG8, RPGRIP1L), telomere maintenance (TNKS, PINX1, AKTIP), and central dopaminergic signaling (MSRA, EBF2).

Abstract TMP91: White Matter Lesion Progression and Cognitive Function Over 5-Years in a Young Susceptible Population

Introduction: Cerebrovascular disease (CVD) and vascular risk factors cluster in families of patients with early-onset coronary artery disease (CAD). In a prior cross-sectional study of 714 healthy middle-aged relatives, we found a high prevalence of early white matter hyperintensity lesions (WMHL) and associated lower cognitive function. Prior WMHL studies have included elderly or demented subjects. Rate of WMHL progression and cognitive decline in quiescent preclinical healthy middle-aged persons remains unknown.

Methods: We studied WMHL and cognitive function over 5 years in 20 healthy relatives, age 50-59, of patients with CAD <60 years of age (one subject per family). Persons with any physical or neurological abnormality were excluded. At baseline and follow-up, subjects underwent 3-Tesla brain magnetic resonance imaging. WMHL volumes were determined using automated segmented volumes measured from FLAIR and MPRAGE using a previously described automated pipeline. Cognitive function was measured with the: Delayed Recall, Digit Symbol Substitution, Verbal Fluency, Digit Span, and Grooved Pegboard tests.

Results: Our 20 participants were: 50% male, 45% African American (AA), average age was 54.1 (SD 1.57) years, 45% with HTN, 15% smokers, and 50% obese. WMHL prevalence was 92%. The average change in WMHL volume was 819 mm 3 (P<0.0008); this was significant for men (p=0.006), women (p= 0.062), AAs (p= 0.019), and European Americans (EA) (p= 0.018). Decline occurred in all cognitive function tests (average P across all tests was <0.001) in all races and sexes. Vascular risk factors did not change from baseline to follow-up.

Conclusion: These striking findings are among the first to demonstrate significant quantitative WMHL volume progression and cognitive decline in a healthy middle-aged population over just 5 years, in AA and EA and both sexes. This small pilot study suggests that early quiescent brain white matter lesion progression and early multidimensional cognitive decline can be detected in healthy susceptible populations using state-of-the-art quantitative volumetric imaging techniques.

Exome-chip meta-analysis identifies association between variation in ANKRD26 and platelet aggregation

Previous genome-wide association studies (GWAS) have identified several variants associated with platelet function phenotypes; however, the proportion of variance explained by the identified variants is mostly small. Rare coding variants, particularly those with high potential for impact on protein structure/function, may have substantial impact on phenotype but are difficult to detect by GWAS. The main purpose of this study was to identify low frequency or rare variants associated with platelet function using genotype data from the Illumina HumanExome Bead Chip. Three family-based cohorts of European ancestry, including ~4,000 total subjects, comprised the discovery cohort and two independent cohorts, one of European and one of African American ancestry, were used for replication. Optical aggregometry in platelet-rich plasma was performed in all the discovery cohorts in response to adenosine diphosphate (ADP), epinephrine, and collagen. Meta-analyses were performed using both gene-based and single nucleotide variant association methods. The gene-based meta-analysis identified a significant association (P = 7.13 × 10^-7) between rare genetic variants in ANKRD26 and ADP-induced platelet aggregation. One of the ANKRD26 SNVs - rs191015656, encoding a threonine to isoleucine substitution predicted to alter protein structure/function, was replicated in Europeans. Aggregation increases of ~20-50% were observed in heterozygotes in all cohorts. Novel genetic signals in ABCG1 and HCP5 were also associated with platelet aggregation to ADP in meta-analyses, although only results for HCP5 could be replicated. The SNV in HCP5 intersects epigenetic signatures in CD41+ megakaryocytes suggesting a new functional role in platelet biology for HCP5. This is the first study to use gene-based association methods from SNV array genotypes to identify rare variants related to platelet function. The molecular mechanisms and pathophysiological relevance for the identified genetic associations requires further study.

Genome-wide meta-analysis associates HLA-DQA1/DRB1 and LPA and lifestyle factors with human longevity

Genomic analysis of longevity offers the potential to illuminate the biology of human aging. Here, using genome-wide association meta-analysis of 606,059 parents’ survival, we discover two regions associated with longevity (HLA-DQA1/DRB1 and LPA). We also validate previous suggestions that APOE, CHRNA3/5, CDKN2A/B, SH2B3 and FOXO3A influence longevity. Next we show that giving up smoking, educational attainment, openness to new experience and high-density lipoprotein (HDL) cholesterol levels are most positively genetically correlated with lifespan while susceptibility to coronary artery disease (CAD), cigarettes smoked per day, lung cancer, insulin resistance and body fat are most negatively correlated. We suggest that the effect of education on lifespan is principally mediated through smoking while the effect of obesity appears to act via CAD. Using instrumental variables, we suggest that an increase of one body mass index unit reduces lifespan by 7 months while 1 year of education adds 11 months to expected lifespan.

Variability in human longevity is genetically influenced. Using genetic data of parental lifespan, the authors identify associations at HLA-DQA/DRB1 and LPA and find that genetic variants that increase educational attainment have a positive effect on lifespan whereas increasing BMI negatively affects lifespan.

Single-trait and multi-trait genome-wide association analyses identify novel loci for blood pressure in African-ancestry populations

Hypertension is a leading cause of global disease, mortality, and disability. While individuals of African descent suffer a disproportionate burden of hypertension and its complications, they have been underrepresented in genetic studies. To identify novel susceptibility loci for blood pressure and hypertension in people of African ancestry, we performed both single and multiple-trait genome-wide association analyses. We analyzed 21 genome-wide association studies comprised of 31,968 individuals of African ancestry, and validated our results with additional 54,395 individuals from multi-ethnic studies. These analyses identified nine loci with eleven independent variants which reached genome-wide significance (P < 1.25×10-8) for either systolic and diastolic blood pressure, hypertension, or for combined traits. Single-trait analyses identified two loci (TARID/TCF21 and LLPH/TMBIM4) and multiple-trait analyses identified one novel locus (FRMD3) for blood pressure. At these three loci, as well as at GRP20/CDH17, associated variants had alleles common only in African-ancestry populations. Functional annotation showed enrichment for genes expressed in immune and kidney cells, as well as in heart and vascular cells/tissues. Experiments driven by these findings and using angiotensin-II induced hypertension in mice showed altered kidney mRNA expression of six genes, suggesting their potential role in hypertension. Our study provides new evidence for genes related to hypertension susceptibility, and the need to study African-ancestry populations in order to identify biologic factors contributing to hypertension.

Discovery and fine-mapping of adiposity loci using high density imputation of genome-wide association studies in individuals of African ancestry: African Ancestry Anthropometry Genetics Consortium

Genome-wide association studies (GWAS) have identified >300 loci associated with measures of adiposity including body mass index (BMI) and waist-to-hip ratio (adjusted for BMI, WHR_adjBMI), but few have been identified through screening of the African ancestry genomes. We performed large scale meta-analyses and replications in up to 52,895 individuals for BMI and up to 23,095 individuals for WHR_adjBMI from the African Ancestry Anthropometry Genetics Consortium (AAAGC) using 1000 Genomes phase 1 imputed GWAS to improve coverage of both common and low frequency variants in the low linkage disequilibrium African ancestry genomes. In the sex-combined analyses, we identified one novel locus (TCF7L2/HABP2) for WHR_adjBMI and eight previously established loci at P < 5×10⁻⁸: seven for BMI, and one for WHR_adjBMI in African ancestry individuals. An additional novel locus (SPRYD7/DLEU2) was identified for WHR_adjBMI when combined with European GWAS. In the sex-stratified analyses, we identified three novel loci for BMI (INTS10/LPL and MLC1 in men, IRX4/IRX2 in women) and four for WHR_adjBMI (SSX2IP, CASC8, PDE3B and ZDHHC1/HSD11B2 in women) in individuals of African ancestry or both African and European ancestry. For four of the novel variants, the minor allele frequency was low (<5%). In the trans-ethnic fine mapping of 47 BMI loci and 27 WHR_adjBMI loci that were locus-wide significant (P < 0.05 adjusted for effective number of variants per locus) from the African ancestry sex-combined and sex-stratified analyses, 26 BMI loci and 17 WHR_adjBMI loci contained ≤ 20 variants in the credible sets that jointly account for 99% posterior probability of driving the associations. The lead variants in 13 of these loci had a high probability of being causal. As compared to our previous HapMap imputed GWAS for BMI and WHR_adjBMI including up to 71,412 and 27,350 African ancestry individuals, respectively, our results suggest that 1000 Genomes imputation showed modest improvement in identifying GWAS loci including low frequency variants. Trans-ethnic meta-analyses further improved fine mapping of putative causal variants in loci shared between the African and European ancestry populations.

Genome-wide association studies (GWAS) have identified >300 genetic regions that influence body size and shape as measured by body mass index (BMI) and waist-to-hip ratio (WHR), respectively, but few have been identified in populations of African ancestry. We conducted large scale high coverage GWAS and replication of these traits in 52,895 and 23,095 individuals of African ancestry, respectively, followed by additional replication in European populations. We identified 10 genome-wide significant loci in all individuals, and an additional seven loci by analyzing men and women separately. We combined African and European ancestry GWAS and were able to narrow down 43 out of 74 African ancestry associated genetic regions to contain small number of putative causal variants. Our results highlight the improvement of applying high density genome coverage and combining multiple ancestries in the identification and refinement of location of genetic regions associated with adiposity traits.

Whole exome sequencing in the Framingham Heart Study identifies rare variation in HYAL2 that influences platelet aggregation

Inhibition of platelet reactivity is a common therapeutic strategy in secondary prevention of cardiovascular disease. Genetic and environmental factors influence inter-individual variation in platelet reactivity. Identifying genes that contribute to platelet reactivity can reveal new biological mechanisms and possible therapeutic targets. Here, we examined rare coding variation to identify genes associated with platelet reactivity in a population-based cohort. To do so, we performed whole exome sequencing in the Framingham Heart Study and conducted single variant and gene-based association tests against platelet reactivity to collagen, adenosine diphosphate (ADP), and epinephrine agonists in up to 1,211 individuals. Single variant tests revealed no significant associations (p<1.44×10^-7), though we observed a suggestive association with previously implicated MRVI1 (rs11042902, p = 1.95×10^-7). Using gene-based association tests of rare and low-frequency variants, we found significant associations of HYAL2 with increased ADP-induced aggregation (p = 1.07×10^-7) and GSTZ1 with increased epinephrine-induced aggregation (p = 1.62×10^-6). HYAL2 also showed suggestive associations with epinephrine-induced aggregation (p = 2.64×10^-5). The rare variants in the HYAL2 gene-based association included a missense variant (N357S) at a known N-glycosylation site and a nonsense variant (Q406*) that removes a glycophosphatidylinositol (GPI) anchor from the resulting protein. These variants suggest that improper membrane trafficking of HYAL2 influences platelet reactivity. We also observed suggestive associations of AR (p = 7.39×10^-6) and MAPRE1 (p = 7.26×10^-6) with ADP-induced reactivity. Our study demonstrates that gene-based tests and other grouping strategies of rare variants are powerful approaches to detect associations in population-based analyses of complex traits not detected by single variant tests and possible new genetic influences on platelet reactivity.

Abstract TMP93: Hypertension is Associated with Disruption of White Matter Tracts in Healthy Middle-aged Persons at Risk for Vascular Disease

Background: Fractional anisotropy (FA) measured by Diffusion Tensor Imaging (DTI) is a sensitive biomarker of white matter (WM) integrity and offers a means to detect early axonal disruption and attendant cognitive impairment secondary to vascular disease.

Hypothesis: We hypothesized that hypertension (HTN), of all known risk factors, is most closely linked to early WM tract injury in apparently healthy high risk relatives of patients with early-onset coronary artery disease (CAD).

Methods: DTI was acquired on a Phillips 3T MRI scanner (slice thickness 2.2 mm, field of view 212X212 mm and b-value 700 mm 2 /s) for relatives of CAD patients. The brain was segmented into 181 regions using the Eve Atlas, Large Deformation Diffeomorphic Metric Mapping paradigm, and 48 regions representing the deep WM were selected for analysis. Mean regional FA (rFA) was assessed using t-tests in persons with and without smoking, HTN, diabetes. Continuous variable association with rFA was tested using Spearman correlations. Covariate-adjusted GEE linear regression (corrected for nonindependence of families) was performed for regions meeting the Bonferroni threshold of 0.001 (0.05/48).

Results: 116 subjects had DTI data (mean age 49.5±11, 57% males, 40.5% African Americans, 23% smokers, 36% hypertensives, 12% diabetics). rFA was not significantly associated with diabetes or current smoking. Moreover, rFA was not significantly correlated with glucose level, LDL, total cholesterol or body mass index. However, patients with HTN showed a significantly lower rFA compared to normotensives in the left (L) genu of corpus callosum, L external capsule, right (R) cingulum, L fornix/ stria terminalis, L and R column and body of fornix, and R sagittal stratum. After adjusting for age, sex and race, this relation remained significant in the R sagittal stratum (P=0.012) and R cingulum (P=0.001). Systolic blood pressure (SBP) correlated with rFA after adjustment in the L superior longitudinal fasciculus as well (unadjusted rho=-0.309, p=0.00073, adjusted beta =-1.08+-0.47, P=0.02).

Conclusions: Of all of the risk factors, as hypothesized, HTN and SBP were significantly associated with early WM tract injury in areas related to cognition in middle-aged persons at risk for vascular disease.

Integrity of Induced Pluripotent Stem Cell (iPSC) Derived Megakaryocytes as Assessed by Genetic and Transcriptomic Analysis

Previously, we have described our feeder-free, xeno-free approach to generate megakaryocytes (MKs) in culture from human induced pluripotent stem cells (iPSCs). Here, we focus specifically on the integrity of these MKs using: (1) genotype discordance between parent cell DNA to iPSC cell DNA and onward to the differentiated MK DNA; (2) genomic structural integrity using copy number variation (CNV); and (3) transcriptomic signatures of the derived MK lines compared to the iPSC lines. We detected a very low rate of genotype discordance; estimates were 0.0001%-0.01%, well below the genotyping error rate for our assay (0.37%). No CNVs were generated in the iPSCs that were subsequently passed on to the MKs. Finally, we observed highly biologically relevant gene sets as being upregulated in MKs relative to the iPSCs: platelet activation, blood coagulation, megakaryocyte development, platelet formation, platelet degranulation, and platelet aggregation. These data strongly support the integrity of the derived MK lines.

Novel genetic loci associated with hippocampal volume

The hippocampal formation is a brain structure integrally involved in episodic memory, spatial navigation, cognition and stress responsiveness. Structural abnormalities in hippocampal volume and shape are found in several common neuropsychiatric disorders. To identify the genetic underpinnings of hippocampal structure here we perform a genome-wide association study (GWAS) of 33,536 individuals and discover six independent loci significantly associated with hippocampal volume, four of them novel. Of the novel loci, three lie within genes (ASTN2, DPP4 and MAST4) and one is found 200 kb upstream of SHH. A hippocampal subfield analysis shows that a locus within the MSRB3 gene shows evidence of a localized effect along the dentate gyrus, subiculum, CA1 and fissure. Further, we show that genetic variants associated with decreased hippocampal volume are also associated with increased risk for Alzheimer's disease (rg=-0.155). Our findings suggest novel biological pathways through which human genetic variation influences hippocampal volume and risk for neuropsychiatric illness.

Genome-wide Trans-ethnic Meta-analysis Identifies Seven Genetic Loci Influencing Erythrocyte Traits and a Role for RBPMS in Erythropoiesis

Genome-wide association studies (GWASs) have identified loci for erythrocyte traits in primarily European ancestry populations. We conducted GWAS meta-analyses of six erythrocyte traits in 71,638 individuals from European, East Asian, and African ancestries using a Bayesian approach to account for heterogeneity in allelic effects and variation in the structure of linkage disequilibrium between ethnicities. We identified seven loci for erythrocyte traits including a locus (RBPMS/GTF2E2) associated with mean corpuscular hemoglobin and mean corpuscular volume. Statistical fine-mapping at this locus pointed to RBPMS at this locus and excluded nearby GTF2E2. Using zebrafish morpholino to evaluate loss of function, we observed a strong in vivo erythropoietic effect for RBPMS but not for GTF2E2, supporting the statistical fine-mapping at this locus and demonstrating that RBPMS is a regulator of erythropoiesis. Our findings show the utility of trans-ethnic GWASs for discovery and characterization of genetic loci influencing hematologic traits.