Obesity Duration, Severity, and Distribution Trajectories and Cardiovascular Disease Risk in the Atherosclerosis Risk in Communities Study

BACKGROUND: Research examining the role of obesity in cardiovascular disease (CVD) often fails to adequately consider heterogeneity in obesity severity, distribution, and duration. METHODS AND RESULTS: We here use multivariate latent class mixed models in the biracial Atherosclerosis Risk in Communities study (N=14 514; mean age=54 years; 55% female) to associate obesity subclasses (derived from body mass index, waist circumference, self-reported weight at age 25, tricep skinfold, and calf circumference across up to four triennial visits) with total mortality, incident CVD, and CVD risk factors. We identified four obesity subclasses, summarized by their body mass index and waist circumference slope as decline (4.1%), stable/slow decline (67.8%), moderate increase (24.6%), and rapid increase (3.6%) subclasses. Compared with participants in the stable/slow decline subclass, the decline subclass was associated with elevated mortality (hazard ratio [HR] 1.45, 95% CI 1.31, 1.60, P<0.0001) and with heart failure (HR 1.41, 95% CI 1.22, 1.63, P<0.0001), stroke (HR 1.53, 95% CI 1.22, 1.92, P=0.0002), and coronary heart disease (HR 1.36, 95% CI 1.14, 1.63, P=0.0008), adjusting for baseline body mass index and CVD risk factor profile. The moderate increase latent class was not associated with any significant differences in CVD risk as compared to the stable/slow decline latent class and was associated with a lower overall risk of mortality (HR 0.85, 95% CI 0.80, 0.90, P<0.0001), despite higher body mass index at baseline. The rapid increase latent class was associated with a higher risk of heart failure versus the stable/slow decline latent class (HR 1.34, 95% CI 1.10, 1.62, P=0.004). CONCLUSIONS: Consideration of heterogeneity and longitudinal changes in obesity measures is needed in clinical care for a more precision-oriented view of CVD risk

Analyses of biomarker traits in diverse UK biobank participants identify associations missed by European-centric analysis strategies

Despite the dramatic underrepresentation of non-European populations in human genetics studies, researchers continue to exclude participants of non-European ancestry, as well as variants rare in European populations, even when these data are available. This practice perpetuates existing research disparities and can lead to important and large effect size associations being missed. Here, we conducted genome-wide association studies (GWAS) of 31 serum and urine biomarker quantitative traits in African (n = 9354), East Asian (n = 2559), and South Asian (n = 9823) ancestry UK Biobank (UKBB) participants. We adjusted for all known GWAS catalog variants for each trait, as well as novel signals identified in a recent European ancestry-focused analysis of UKBB participants. We identify 7 novel signals in African ancestry and 2 novel signals in South Asian ancestry participants (p < 1.61E−10). Many of these signals are highly plausible, including a cis pQTL for the gene encoding gamma-glutamyl transferase and PIEZO1 and G6PD variants with impacts on HbA1c through likely erythrocytic mechanisms. This work illustrates the importance of using the genetic data we already have in diverse populations, with novel discoveries possible in even modest sample sizes

Ancestry-specific associations identified in genome-wide combined-phenotype study of red blood cell traits emphasize benefits of diversity in genomics

Background: Quantitative red blood cell (RBC) traits are highly polygenic clinically relevant traits, with approximately 500 reported GWAS loci. The majority of RBC trait GWAS have been performed in European- or East Asian-ancestry populations, despite evidence that rare or ancestry-specific variation contributes substantially to RBC trait heritability. Recently developed combined-phenotype methods which leverage genetic trait correlation to improve statistical power have not yet been applied to these traits. Here we leveraged correlation of seven quantitative RBC traits in performing a combined-phenotype analysis in a multi-ethnic study population. Results: We used the adaptive sum of powered scores (aSPU) test to assess combined-phenotype associations between ~ 21 million SNPs and seven RBC traits in a multi-ethnic population (maximum n = 67,885 participants; 24% African American, 30% Hispanic/Latino, and 43% European American; 76% female). Thirty-nine loci in our multi-ethnic population contained at least one significant association signal (p 5%) across all ancestral populations. Nineteen additional independent association signals were identified at seven known loci (HFE, KIT, HBS1L/MYB, CITED2/FILNC1, ABO, HBA1/2, and PLIN4/5). For example, the HBA1/2 locus contained 14 conditionally independent association signals, 11 of which were previously unreported and are specific to African and Amerindian ancestries. One variant in this region was common in all ancestries, but exhibited a narrower LD block in African Americans than European Americans or Hispanics/Latinos. GTEx eQTL analysis of all independent lead SNPs yielded 31 significant associations in relevant tissues, over half of which were not at the gene immediately proximal to the lead SNP. Conclusion: This work identified seven loci containing multiple independent association signals for RBC traits using a combined-phenotype approach, which may improve discovery in genetically correlated traits. Highly complex genetic architecture at the HBA1/2 locus was only revealed by the inclusion of African Americans and Hispanics/Latinos, underscoring the continued importance of expanding large GWAS to include ancestrally diverse populations. © 2020 The Author(s)

Multi-ethnic GWAS and fine-mapping of glycaemic traits identify novel loci in the PAGE Study

Aims/hypothesis: Type 2 diabetes is a growing global public health challenge. Investigating quantitative traits, including fasting glucose, fasting insulin and HbA1c, that serve as early markers of type 2 diabetes progression may lead to a deeper understanding of the genetic aetiology of type 2 diabetes development. Previous genome-wide association studies (GWAS) have identified over 500 loci associated with type 2 diabetes, glycaemic traits and insulin-related traits. However, most of these findings were based only on populations of European ancestry. To address this research gap, we examined the genetic basis of fasting glucose, fasting insulin and HbA1c in participants of the diverse Population Architecture using Genomics and Epidemiology (PAGE) Study. Methods: We conducted a GWAS of fasting glucose (n = 52,267), fasting insulin (n = 48,395) and HbA1c (n = 23,357) in participants without diabetes from the diverse PAGE Study (23% self-reported African American, 46% Hispanic/Latino, 40% European, 4% Asian, 3% Native Hawaiian, 0.8% Native American), performing transethnic and population-specific GWAS meta-analyses, followed by fine-mapping to identify and characterise novel loci and independent secondary signals in known loci. Results: Four novel associations were identified (p < 5 × 10−9), including three loci associated with fasting insulin, and a novel, low-frequency African American-specific locus associated with fasting glucose. Additionally, seven secondary signals were identified, including novel independent secondary signals for fasting glucose at the known GCK locus and for fasting insulin at the known PPP1R3B locus in transethnic meta-analysis. Conclusions/interpretation: Our findings provide new insights into the genetic architecture of glycaemic traits and highlight the continued importance of conducting genetic studies in diverse populations. Data availability: Full summary statistics from each of the population-specific and transethnic results are available at NHGRI-EBI GWAS catalog (https://www.ebi.ac.uk/gwas/downloads/summary-statistics)

Multi-ethnic genome-wide association analyses of white blood cell and platelet traits in the Population Architecture using Genomics and Epidemiology (PAGE) study

Background: Circulating white blood cell and platelet traits are clinically linked to various disease outcomes and differ across individuals and ancestry groups. Genetic factors play an important role in determining these traits and many loci have been identified. However, most of these findings were identified in populations of European ancestry (EA), with African Americans (AA), Hispanics/Latinos (HL), and other races/ethnicities being severely underrepresented. Results: We performed ancestry-combined and ancestry-specific genome-wide association studies (GWAS) for white blood cell and platelet traits in the ancestrally diverse Population Architecture using Genomics and Epidemiology (PAGE) Study, including 16,201 AA, 21,347 HL, and 27,236 EA participants. We identified six novel findings at suggestive significance (P < 5E-8), which need confirmation, and independent signals at six previously established regions at genome-wide significance (P < 2E-9). We confirmed multiple previously reported genome-wide significant variants in the single variant association analysis and multiple genes using PrediXcan. Evaluation of loci reported from a Euro-centric GWAS indicated attenuation of effect estimates in AA and HL compared to EA populations. Conclusions: Our results highlighted the potential to identify ancestry-specific and ancestry-agnostic variants in participants with diverse backgrounds and advocate for continued efforts in improving inclusion of racially/ethnically diverse populations in genetic association studies for complex traits

Supplemental Association of Clonal Hematopoiesis With Incident Heart Failure

Background: Age-related clonal hematopoiesis of indeterminate potential (CHIP), defined as clonally expanded leukemogenic sequence variations (particularly in DNMT3A, TET2, ASXL1, and JAK2) in asymptomatic individuals, is associated with cardiovascular events, including recurrent heart failure (HF). Objectives: This study sought to evaluate whether CHIP is associated with incident HF. Methods: CHIP status was obtained from whole exome or genome sequencing of blood DNA in participants without prevalent HF or hematological malignancy from 5 cohorts. Cox proportional hazards models were performed within each cohort, adjusting for demographic and clinical risk factors, followed by fixed-effect meta-analyses. Large CHIP clones (defined as variant allele frequency >10%), HF with or without baseline coronary heart disease, and left ventricular ejection fraction were evaluated in secondary analyses. Results: Of 56,597 individuals (59% women, mean age 58 years at baseline), 3,406 (6%) had CHIP, and 4,694 developed HF (8.3%) over up to 20 years of follow-up. CHIP was prospectively associated with a 25% increased risk of HF in meta-analysis (hazard ratio: 1.25; 95% confidence interval: 1.13-1.38) with consistent associations across cohorts. ASXL1, TET2, and JAK2 sequence variations were each associated with an increased risk of HF, whereas DNMT3A sequence variations were not associated with HF. Secondary analyses suggested large CHIP was associated with a greater risk of HF (hazard ratio: 1.29; 95% confidence interval: 1.15-1.44), and the associations for CHIP on HF with and without prior coronary heart disease were homogenous. ASXL1 sequence variations were associated with reduced left ventricular ejection fraction. Conclusions: CHIP, particularly sequence variations in ASXL1, TET2, and JAK2, represents a new risk factor for HF

Genome-wide association study of heart rate and its variability in Hispanic/Latino cohorts

Background Although time–domain measures of heart rate variability (HRV) are used to estimate cardiac autonomic tone and disease risk in multiethnic populations, the genetic epidemiology of HRV in Hispanics/Latinos has not been characterized. Objective The purpose of this study was to conduct a genome-wide association study of heart rate (HR) and its variability in the Hispanic Community Health Study/Study of Latinos, Multi-Ethnic Study of Atherosclerosis, and Women's Health Initiative Hispanic SNP-Health Association Resource project (n = 13,767). Methods We estimated HR (bpm), standard deviation of normal-to-normal interbeat intervals (SDNN, ms), and root mean squared difference in successive, normal-to-normal interbeat intervals (RMSSD, ms) from resting, standard 12-lead ECGs. We estimated associations between each phenotype and 17 million genotyped or imputed single nucleotide polymorphisms (SNPs), accounting for relatedness and adjusting for age, sex, study site, and ancestry. Cohort-specific estimates were combined using fixed-effects, inverse-variance meta-analysis. We investigated replication for select SNPs exceeding genome-wide (P &lt;5 × 10–8) or suggestive (P &lt;10–6) significance thresholds. Results Two genome-wide significant SNPs replicated in a European ancestry cohort, 1 one for RMSSD (rs4963772; chromosome 12) and another for SDNN (rs12982903; chromosome 19). A suggestive SNP for HR (rs236352; chromosome 6) replicated in an African-American cohort. Functional annotation of replicated SNPs in cardiac and neuronal tissues identified potentially causal variants and mechanisms. Conclusion This first genome-wide association study of HRV and HR in Hispanics/Latinos underscores the potential for even modestly sized samples of non-European ancestry to inform the genetic epidemiology of complex traits

Transcriptome-wide association study of blood cell traits in african ancestry and hispanic/latino populations

Background: Thousands of genetic variants have been associated with hematological traits, though target genes remain unknown at most loci. Moreover, limited analyses have been conducted in African ancestry and Hispanic/Latino populations; hematological trait associated variants more common in these populations have likely been missed. Methods: To derive gene expression prediction models, we used ancestry-stratified datasets from the Multi-Ethnic Study of Atherosclerosis (MESA, including n = 229 African American and n = 381 Hispanic/Latino participants, monocytes) and the Depression Genes and Networks study (DGN, n = 922 European ancestry participants, whole blood). We then performed a transcriptome-wide association study (TWAS) for platelet count, hemoglobin, hematocrit, and white blood cell count in African (n = 27,955) and Hispanic/Latino (n = 28,324) ancestry participants. Results: Our results revealed 24 suggestive signals (p &lt; 1 × 10−4 ) that were conditionally distinct from known GWAS identified variants and successfully replicated these signals in European ancestry subjects from UK Biobank. We found modestly improved correlation of predicted and measured gene expression in an independent African American cohort (the Genetic Epidemiology Network of Arteriopathy (GENOA) study (n = 802), lymphoblastoid cell lines) using the larger DGN reference panel; however, some genes were well predicted using MESA but not DGN. Conclusions: These analyses demonstrate the importance of performing TWAS and other genetic analyses across diverse populations and of balancing sample size and ancestry background matching when selecting a TWAS reference panel

Clonal hematopoiesis associated with epigenetic aging and clinical outcomes

Clonal hematopoiesis of indeterminate potential (CHIP) is a common precursor state for blood cancers that most frequently occurs due to mutations in the DNA-methylation modifying enzymes DNMT3A or TET2. We used DNA-methylation array and whole-genome sequencing data from four cohorts together comprising 5522 persons to study the association between CHIP, epigenetic clocks, and health outcomes. CHIP was strongly associated with epigenetic age acceleration, defined as the residual after regressing epigenetic clock age on chronological age, in several clocks, ranging from 1.31 years (GrimAge, p &lt; 8.6 × 10−7) to 3.08 years (EEAA, p &lt; 3.7 × 10−18). Mutations in most CHIP genes except DNA-damage response genes were associated with increases in several measures of age acceleration. CHIP carriers with mutations in multiple genes had the largest increases in age acceleration and decrease in estimated telomere length. Finally, we found that ~40% of CHIP carriers had acceleration &gt;0 in both Hannum and GrimAge (referred to as AgeAccelHG+). This group was at high risk of all-cause mortality (hazard ratio 2.90, p &lt; 4.1 × 10−8) and coronary heart disease (CHD) (hazard ratio 3.24, p &lt; 9.3 × 10−6) compared to those who were CHIP−/AgeAccelHG−. In contrast, the other ~60% of CHIP carriers who were AgeAccelHG− were not at increased risk of these outcomes. In summary, CHIP is strongly linked to age acceleration in multiple clocks, and the combination of CHIP and epigenetic aging may be used to identify a population at high risk for adverse outcomes and who may be a target for clinical interventions

Comparison of Proteomic Assessment Methods in Multiple Cohort Studies

Novel proteomics platforms, such as the aptamer-based SOMAscan platform, can quantify large numbers of proteins efficiently and cost-effectively and are rapidly growing in popularity. However, comparisons to conventional immunoassays remain underexplored, leaving investigators unsure when cross-assay comparisons are appropriate. The correlation of results from immunoassays with relative protein quantification is explored by SOMAscan. For 63 proteins assessed in two chronic obstructive pulmonary disease (COPD) cohorts, subpopulations and intermediate outcome measures in COPD Study (SPIROMICS), and COPDGene, using myriad rules based medicine multiplex immunoassays and SOMAscan, Spearman correlation coefficients range from −0.13 to 0.97, with a median correlation coefficient of ≈0.5 and consistent results across cohorts. A similar range is observed for immunoassays in the population-based Multi-Ethnic Study of Atherosclerosis and for other assays in COPDGene and SPIROMICS. Comparisons of relative quantification from the antibody-based Olink platform and SOMAscan in a small cohort of myocardial infarction patients also show a wide correlation range. Finally, cis pQTL data, mass spectrometry aptamer confirmation, and other publicly available data are integrated to assess relationships with observed correlations. Correlation between proteomics assays shows a wide range and should be carefully considered when comparing and meta-analyzing proteomics data across assays and studies