A parametric bootstrap approach for computing confidence intervals for genetic correlations with application to genetically determined protein-protein networks

Genetic correlation refers to the correlation between genetic determinants of a pair of traits. When using individual-level data, it is typically estimated based on a bivariate model specification where the correlation between the two variables is identifiable and can be estimated from a covariance model that incorporates the genetic relationship between individuals, e.g., using a pre-specified kinship matrix. Inference relying on asymptotic normality of the genetic correlation parameter estimates may be inaccurate when the sample size is low, when the genetic correlation is close to the boundary of the parameter space, and when the heritability of at least one of the traits is low. We address this problem by developing a parametric bootstrap procedure to construct confidence intervals for genetic correlation estimates. The procedure simulates paired traits under a range of heritability and genetic correlation parameters, and it uses the population structure encapsulated by the kinship matrix. Heritabilities and genetic correlations are estimated using the close-form, method of moment, Haseman-Elston regression estimators. The proposed parametric bootstrap procedure is especially useful when genetic correlations are computed on pairs of thousands of traits measured on the same exact set of individuals. We demonstrate the parametric bootstrap approach on a proteomics dataset from the Jackson Heart Study.

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.

Clonal Hematopoiesis of Indeterminate Potential (CHIP) and Incident Type 2 Diabetes Risk

OBJECTIVE

Clonal hematopoiesis of indeterminate potential (CHIP) is an aging-related accumulation of somatic mutations in hematopoietic stem cells, leading to clonal expansion. CHIP presence has been implicated in atherosclerotic coronary heart disease (CHD) and all-cause mortality, but its association with incident type 2 diabetes (T2D) is unknown. We hypothesized that CHIP is associated with elevated risk of T2D.

RESEARCH DESIGN AND METHODS

CHIP was derived from whole-genome sequencing of blood DNA in the National Heart, Lung, and Blood Institute Trans-Omics for Precision Medicine (TOPMed) prospective cohorts. We performed analysis for 17,637 participants from six cohorts, without prior T2D, cardiovascular disease, or cancer. We evaluated baseline CHIP versus no CHIP prevalence with incident T2D, including associations with DNMT3A, TET2, ASXL1, JAK2, and TP53 variants. We estimated multivariable-adjusted hazard ratios (HRs) and 95% CIs with adjustment for age, sex, BMI, smoking, alcohol, education, self-reported race/ethnicity, and combined cohorts' estimates via fixed-effects meta-analysis.

RESULTS

Mean (SD) age was 63.4 (11.5) years, 76% were female, and CHIP prevalence was 6.0% (n = 1,055) at baseline. T2D was diagnosed in n = 2,467 over mean follow-up of 9.8 years. Participants with CHIP had 23% (CI 1.04, 1.45) higher risk of T2D than those with no CHIP. Specifically, higher risk was for TET2 (HR 1.48; CI 1.05, 2.08) and ASXL1 (HR 1.76; CI 1.03, 2.99) mutations; DNMT3A was nonsignificant (HR 1.15; CI 0.93, 1.43). Statistical power was limited for JAK2 and TP53 analyses.

CONCLUSIONS

CHIP was associated with higher incidence of T2D. CHIP mutations located on genes implicated in CHD and mortality were also related to T2D, suggesting shared aging-related pathology.

A parametric bootstrap approach for computing confidence intervals for genetic correlations with application to genetically-determined protein-protein networks

Genetic correlation refers to the correlation between genetic determinants of a pair of traits. When using individual-level data, it is typically estimated based on a bivariate model specification where the correlation between the two variables is identifiable and can be estimated from a covariance model that incorporates the genetic relationship between individuals, e.g., using a pre-specified kinship matrix. Inference relying on asymptotic normality of the genetic correlation parameter estimates may be inaccurate when the sample size is low, when the genetic correlation is close to the boundary of the parameter space, and when the heritability of at least one of the traits is low. We address this problem by developing a parametric bootstrap procedure to construct confidence intervals for genetic correlation estimates. The procedure simulates paired traits under a range of heritability and genetic correlation parameters, and it uses the population structure encapsulated by the kinship matrix. Heritabilities and genetic correlations are estimated using the close-form, method of moment, Haseman-Elston regression estimators. The proposed parametric bootstrap procedure is especially useful when genetic correlations are computed on pairs of thousands of traits measured on the same exact set of individuals. We demonstrate the parametric bootstrap approach on a proteomics dataset from the Jackson Heart Study.

Association Between Whole Blood-Derived Mitochondrial DNA Copy Number, Low-Density Lipoprotein Cholesterol, and Cardiovascular Disease Risk

Background The relationship between mitochondrial DNA copy number (mtDNA CN) and cardiovascular disease remains elusive. Methods and Results We performed cross-sectional and prospective association analyses of blood-derived mtDNA CN and cardiovascular disease outcomes in 27 316 participants in 8 cohorts of multiple racial and ethnic groups with whole-genome sequencing. We also performed Mendelian randomization to explore causal relationships of mtDNA CN with coronary heart disease (CHD) and cardiometabolic risk factors (obesity, diabetes, hypertension, and hyperlipidemia). P<0.01 was used for significance. We validated most of the previously reported associations between mtDNA CN and cardiovascular disease outcomes. For example, 1-SD unit lower level of mtDNA CN was associated with 1.08 (95% CI, 1.04-1.12; P<0.001) times the hazard for developing incident CHD, adjusting for covariates. Mendelian randomization analyses showed no causal effect from a lower level of mtDNA CN to a higher CHD risk (β=0.091; P=0.11) or in the reverse direction (β=-0.012; P=0.076). Additional bidirectional Mendelian randomization analyses revealed that low-density lipoprotein cholesterol had a causal effect on mtDNA CN (β=-0.084; P<0.001), but the reverse direction was not significant (P=0.059). No causal associations were observed between mtDNA CN and obesity, diabetes, and hypertension, in either direction. Multivariable Mendelian randomization analyses showed no causal effect of CHD on mtDNA CN, controlling for low-density lipoprotein cholesterol level (P=0.52), whereas there was a strong direct causal effect of higher low-density lipoprotein cholesterol on lower mtDNA CN, adjusting for CHD status (β=-0.092; P<0.001). Conclusions Our findings indicate that high low-density lipoprotein cholesterol may underlie the complex relationships between mtDNA CN and vascular atherosclerosis.

Protein-metabolite association studies identify novel proteomic determinants of metabolite levels in human plasma

Although many novel gene-metabolite and gene-protein associations have been identified using high-throughput biochemical profiling, systematic studies that leverage human genetics to illuminate causal relationships between circulating proteins and metabolites are lacking. Here, we performed protein-metabolite association studies in 3,626 plasma samples from three human cohorts. We detected 171,800 significant protein-metabolite pairwise correlations between 1,265 proteins and 365 metabolites, including established relationships in metabolic and signaling pathways such as the protein thyroxine-binding globulin and the metabolite thyroxine, as well as thousands of new findings. In Mendelian randomization (MR) analyses, we identified putative causal protein-to-metabolite associations. We experimentally validated top MR associations in proof-of-concept plasma metabolomics studies in three murine knockout strains of key protein regulators. These analyses identified previously unrecognized associations between bioactive proteins and metabolites in human plasma. We provide publicly available data to be leveraged for studies in human metabolism and disease.

The Relationship of Duffy Gene Polymorphism, High Sensitivity C-Reactive Protein, and Long-term Outcomes

Background

Black adults have higher incidence of all-cause death and worse cardiovascular outcomes when compared to other populations. The Duffy chemokine receptor is not expressed in a large majority of Black adults and the clinical implications of this are unclear.

Methods

Here, we investigated the relationship of Duffy receptor status, high-sensitivity C-reactive protein (hs-CRP), and long-term cardiovascular outcomes in Black members of two contemporary, longitudinal cohort studies (the Jackson Heart Study and Multi-Ethnic Study of Atherosclerosis). Data on 4,307 Black participants (2,942 Duffy null and 1,365 Duffy receptor positive, as defined using Single Nucleotide Polymorphism (SNP) rs2814778) were included in this analysis.

Results

Duffy null was not independently associated with elevated levels of serum hs-CRP levels once conditioning for known CRP locus alleles in linkage disequilibrium with the Duffy gene. Duffy null status was not found to be independently associated with higher incidence of all-cause mortality or secondary outcomes after adjusting for possible confounders in Black participants.

Conclusions

These findings suggest that increased levels of hs-CRP found in Duffy null individuals is due to co-inheritance of CRP alleles known to influence circulating levels hs-CRP and that Duffy null status was not associated with worse adverse outcomes over the follow-up period in this cohort of well-balanced Black participants.

Systemic Markers of Lung Function and Forced Expiratory Volume in 1 Second Decline across Diverse Cohorts

Rationale: Chronic obstructive pulmonary disease (COPD) is a complex disease characterized by airway obstruction and accelerated lung function decline. Our understanding of systemic protein biomarkers associated with COPD remains incomplete. Objectives: To determine what proteins and pathways are associated with impaired pulmonary function in a diverse population. Methods: We studied 6,722 participants across six cohort studies with both aptamer-based proteomic and spirometry data (4,566 predominantly White participants in a discovery analysis and 2,156 African American cohort participants in a validation). In linear regression models, we examined protein associations with baseline forced expiratory volume in 1 second (FEV1) and FEV1/forced vital capacity (FVC). In linear mixed effects models, we investigated the associations of baseline protein levels with rate of FEV1 decline (ml/yr) in 2,777 participants with up to 7 years of follow-up spirometry. Results: We identified 254 proteins associated with FEV1 in our discovery analyses, with 80 proteins validated in the Jackson Heart Study. Novel validated protein associations include kallistatin serine protease inhibitor, growth differentiation factor 2, and tumor necrosis factor-like weak inducer of apoptosis (discovery β = 0.0561, Q = 4.05 × 10-10; β  = 0.0421, Q = 1.12 × 10-3; and β = 0.0358, Q = 1.67 × 10-3, respectively). In longitudinal analyses within cohorts with follow-up spirometry, we identified 15 proteins associated with FEV1 decline (Q < 0.05), including elafin leukocyte elastase inhibitor and mucin-associated TFF2 (trefoil factor 2; β = -4.3 ml/yr, Q = 0.049; β = -6.1 ml/yr, Q = 0.032, respectively). Pathways and processes highlighted by our study include aberrant extracellular matrix remodeling, enhanced innate immune response, dysregulation of angiogenesis, and coagulation. Conclusions: In this study, we identify and validate novel biomarkers and pathways associated with lung function traits in a racially diverse population. In addition, we identify novel protein markers associated with FEV1 decline. Several protein findings are supported by previously reported genetic signals, highlighting the plausibility of certain biologic pathways. These novel proteins might represent markers for risk stratification, as well as novel molecular targets for treatment of COPD.

Investigating Gene-Diet Interactions Impacting the Association Between Macronutrient Intake and Glycemic Traits

Few studies have demonstrated reproducible gene-diet interactions (GDIs) impacting metabolic disease risk factors, likely due in part to measurement error in dietary intake estimation and insufficient capture of rare genetic variation. We aimed to identify GDIs across the genetic frequency spectrum impacting the macronutrient-glycemia relationship in genetically and culturally diverse cohorts. We analyzed 33,187 participants free of diabetes from 10 National Heart, Lung, and Blood Institute Trans-Omics for Precision Medicine program cohorts with whole-genome sequencing, self-reported diet, and glycemic trait data. We fit cohort-specific, multivariable-adjusted linear mixed models for the effect of diet, modeled as an isocaloric substitution of carbohydrate for fat, and its interactions with common and rare variants genome-wide. In main effect meta-analyses, participants consuming more carbohydrate had modestly lower glycemic trait values (e.g., for glycated hemoglobin [HbA1c], -0.013% HbA1c/250 kcal substitution). In GDI meta-analyses, a common African ancestry-enriched variant (rs79762542) reached study-wide significance and replicated in the UK Biobank cohort, indicating a negative carbohydrate-HbA1c association among major allele homozygotes only. Simulations revealed that >150,000 samples may be necessary to identify similar macronutrient GDIs under realistic assumptions about effect size and measurement error. These results generate hypotheses for further exploration of modifiable metabolic disease risk in additional cohorts with African ancestry.

ARTICLE HIGHLIGHTS

We aimed to identify genetic modifiers of the dietary macronutrient-glycemia relationship using whole-genome sequence data from 10 Trans-Omics for Precision Medicine program cohorts. Substitution models indicated a modest reduction in glycemia associated with an increase in dietary carbohydrate at the expense of fat. Genome-wide interaction analysis identified one African ancestry-enriched variant near the FRAS1 gene that may interact with macronutrient intake to influence hemoglobin A1c. Simulation-based power calculations accounting for measurement error suggested that substantially larger sample sizes may be necessary to discover further gene-macronutrient interactions.

Canonical correlation analysis for multi-omics: Application to cross-cohort analysis

Integrative approaches that simultaneously model multi-omics data have gained increasing popularity because they provide holistic system biology views of multiple or all components in a biological system of interest. Canonical correlation analysis (CCA) is a correlation-based integrative method designed to extract latent features shared between multiple assays by finding the linear combinations of features-referred to as canonical variables (CVs)-within each assay that achieve maximal across-assay correlation. Although widely acknowledged as a powerful approach for multi-omics data, CCA has not been systematically applied to multi-omics data in large cohort studies, which has only recently become available. Here, we adapted sparse multiple CCA (SMCCA), a widely-used derivative of CCA, to proteomics and methylomics data from the Multi-Ethnic Study of Atherosclerosis (MESA) and Jackson Heart Study (JHS). To tackle challenges encountered when applying SMCCA to MESA and JHS, our adaptations include the incorporation of the Gram-Schmidt (GS) algorithm with SMCCA to improve orthogonality among CVs, and the development of Sparse Supervised Multiple CCA (SSMCCA) to allow supervised integration analysis for more than two assays. Effective application of SMCCA to the two real datasets reveals important findings. Applying our SMCCA-GS to MESA and JHS, we identified strong associations between blood cell counts and protein abundance, suggesting that adjustment of blood cell composition should be considered in protein-based association studies. Importantly, CVs obtained from two independent cohorts also demonstrate transferability across the cohorts. For example, proteomic CVs learned from JHS, when transferred to MESA, explain similar amounts of blood cell count phenotypic variance in MESA, explaining 39.0% ~ 50.0% variation in JHS and 38.9% ~ 49.1% in MESA. Similar transferability was observed for other omics-CV-trait pairs. This suggests that biologically meaningful and cohort-agnostic variation is captured by CVs. We anticipate that applying our SMCCA-GS and SSMCCA on various cohorts would help identify cohort-agnostic biologically meaningful relationships between multi-omics data and phenotypic traits.

CYP4F2 is a human-specific determinant of circulating N-acyl amino acid levels

N-acyl amino acids are a large family of circulating lipid metabolites that modulate energy expenditure and fat mass in rodents. However, little is known about the regulation and potential cardiometabolic functions of N-acyl amino acids in humans. Here, we analyze the cardiometabolic phenotype associations and genomic associations of four plasma N-acyl amino acids (N-oleoyl-leucine, N-oleoyl-phenylalanine, N-oleoyl-serine, and N-oleoyl-glycine) in 2,351 individuals from the Jackson Heart Study. We find that plasma levels of specific N-acyl amino acids are associated with cardiometabolic disease endpoints independent of free amino acid plasma levels and in patterns according to the amino acid head group. By integrating whole genome sequencing data with N-acyl amino acid levels, we identify that the genetic determinants of N-acyl amino acid levels also cluster according to amino acid head group. Furthermore, we identify the CYP4F2 locus as a genetic determinant of plasma N-oleoyl-leucine and N-oleoyl-phenylalanine levels in human plasma. In experimental studies, we demonstrate that CYP4F2-mediated hydroxylation of N-oleoyl-leucine and N-oleoyl-phenylalanine results in metabolic diversification and production of many previously unknown lipid metabolites with varying characteristics of the fatty acid tail group, including several that structurally resemble fatty acid hydroxy fatty acids (FAHFAs). These studies provide a structural framework for understanding the regulation and disease-associations of N-acyl amino acids in humans and identify that the diversity of this lipid signaling family can be significantly expanded through CYP4F-mediated ω-hydroxylation.

Plasma proteomic changes in response to exercise training are associated with cardiorespiratory fitness adaptations

Regular exercise leads to widespread salutary effects, and there is increasing recognition that exercise-stimulated circulating proteins can impart health benefits. Despite this, limited data exist regarding the plasma proteomic changes that occur in response to regular exercise. Here, we perform large-scale plasma proteomic profiling in 654 healthy human study participants before and after a supervised, 20-week endurance exercise training intervention. We identify hundreds of circulating proteins that are modulated, many of which are known to be secreted. We highlight proteins involved in angiogenesis, iron homeostasis, and the extracellular matrix, many of which are novel, including training-induced increases in fibroblast activation protein (FAP), a membrane-bound and circulating protein relevant in body-composition homeostasis. We relate protein changes to training-induced maximal oxygen uptake adaptations and validate our top findings in an external exercise cohort. Furthermore, we show that FAP is positively associated with survival in 3 separate, population-based cohorts.

Protein Markers of Diabetes Discovered in an African American Cohort

Proteomics has been used to study type 2 diabetes, but the majority of available data are from White participants. Here, we extend prior work by analyzing a large cohort of self-identified African Americans in the Jackson Heart Study (n = 1,313). We found 325 proteins associated with incident diabetes after adjusting for age, sex, and sample batch (false discovery rate q < 0.05) measured using a single-stranded DNA aptamer affinity-based method on fasting plasma samples. A subset was independent of established markers of diabetes development pathways, such as adiposity, glycemia, and/or insulin resistance, suggesting potential novel biological processes associated with disease development. Thirty-six associations remained significant after additional adjustments for BMI, fasting plasma glucose, cholesterol levels, hypertension, statin use, and renal function. Twelve associations, including the top associations of complement factor H, formimidoyltransferase cyclodeaminase, serine/threonine-protein kinase 17B, and high-mobility group protein B1, were replicated in a meta-analysis of two self-identified White cohorts-the Framingham Heart Study and the Malmö Diet and Cancer Study-supporting the generalizability of these biomarkers. A selection of these diabetes-associated proteins also improved risk prediction. Thus, we uncovered both novel and broadly generalizable associations by studying a diverse population, providing a more complete understanding of the diabetes-associated proteome.

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.

Genetic architecture of spatial electrical biomarkers for cardiac arrhythmia and relationship with cardiovascular disease

The 3-dimensional spatial and 2-dimensional frontal QRS-T angles are measures derived from the vectorcardiogram. They are independent risk predictors for arrhythmia, but the underlying biology is unknown. Using multi-ancestry genome-wide association studies we identify 61 (58 previously unreported) loci for the spatial QRS-T angle (N = 118,780) and 11 for the frontal QRS-T angle (N = 159,715). Seven out of the 61 spatial QRS-T angle loci have not been reported for other electrocardiographic measures. Enrichments are observed in pathways related to cardiac and vascular development, muscle contraction, and hypertrophy. Pairwise genome-wide association studies with classical ECG traits identify shared genetic influences with PR interval and QRS duration. Phenome-wide scanning indicate associations with atrial fibrillation, atrioventricular block and arterial embolism and genetically determined QRS-T angle measures are associated with fascicular and bundle branch block (and also atrioventricular block for the frontal QRS-T angle). We identify potential biology involved in the QRS-T angle and their genetic relationships with cardiovascular traits and diseases, may inform future research and risk prediction.

CYP4F2 is a human-specific determinant of circulating N-acyl amino acid levels

N-acyl amino acids are a large family of circulating lipid metabolites that modulate energy expenditure and fat mass in rodents. However, little is known about the regulation and potential cardiometabolic functions of N-acyl amino acids in humans. Here, we analyze the cardiometabolic phenotype associations and genetic regulation of four plasma N-fatty acyl amino acids (N-oleoyl-leucine, N-oleoyl-phenylalanine, N-oleoyl-serine, and N-oleoyl-glycine) in 2,351 individuals from the Jackson Heart Study. N-oleoyl-leucine and N-oleoyl-phenylalanine were positively associated with traits related to energy balance, including body mass index, waist circumference, and subcutaneous adipose tissue. In addition, we identify the CYP4F2 locus as a human-specific genetic determinant of plasma N-oleoyl-leucine and N-oleoyl-phenylalanine levels. In vitro, CYP4F2-mediated hydroxylation of N-oleoyl-leucine and N-oleoyl-phenylalanine results in metabolic diversification and production of many previously unknown lipid metabolites with varying characteristics of the fatty acid tail group, including several that structurally resemble fatty acid hydroxy fatty acids (FAHFAs). By contrast, FAAH-regulated N-oleoyl-glycine and N-oleoyl-serine were inversely associated with traits related to glucose and lipid homeostasis. These data uncover a human-specific enzymatic node for the metabolism of a subset of N-fatty acyl amino acids and establish a framework for understanding the cardiometabolic roles of individual N-fatty acyl amino acids in humans.

Powerful, scalable and resource-efficient meta-analysis of rare variant associations in large whole genome sequencing studies

Meta-analysis of whole genome sequencing/whole exome sequencing (WGS/WES) studies provides an attractive solution to the problem of collecting large sample sizes for discovering rare variants associated with complex phenotypes. Existing rare variant meta-analysis approaches are not scalable to biobank-scale WGS data. Here we present MetaSTAAR, a powerful and resource-efficient rare variant meta-analysis framework for large-scale WGS/WES studies. MetaSTAAR accounts for relatedness and population structure, can analyze both quantitative and dichotomous traits and boosts the power of rare variant tests by incorporating multiple variant functional annotations. Through meta-analysis of four lipid traits in 30,138 ancestrally diverse samples from 14 studies of the Trans Omics for Precision Medicine (TOPMed) Program, we show that MetaSTAAR performs rare variant meta-analysis at scale and produces results comparable to using pooled data. Additionally, we identified several conditionally significant rare variant associations with lipid traits. We further demonstrate that MetaSTAAR is scalable to biobank-scale cohorts through meta-analysis of TOPMed WGS data and UK Biobank WES data of ~200,000 samples.

Implicating genes, pleiotropy, and sexual dimorphism at blood lipid loci through multi-ancestry meta-analysis

BACKGROUND

Genetic variants within nearly 1000 loci are known to contribute to modulation of blood lipid levels. However, the biological pathways underlying these associations are frequently unknown, limiting understanding of these findings and hindering downstream translational efforts such as drug target discovery.

RESULTS

To expand our understanding of the underlying biological pathways and mechanisms controlling blood lipid levels, we leverage a large multi-ancestry meta-analysis (N = 1,654,960) of blood lipids to prioritize putative causal genes for 2286 lipid associations using six gene prediction approaches. Using phenome-wide association (PheWAS) scans, we identify relationships of genetically predicted lipid levels to other diseases and conditions. We confirm known pleiotropic associations with cardiovascular phenotypes and determine novel associations, notably with cholelithiasis risk. We perform sex-stratified GWAS meta-analysis of lipid levels and show that 3-5% of autosomal lipid-associated loci demonstrate sex-biased effects. Finally, we report 21 novel lipid loci identified on the X chromosome. Many of the sex-biased autosomal and X chromosome lipid loci show pleiotropic associations with sex hormones, emphasizing the role of hormone regulation in lipid metabolism.

CONCLUSIONS

Taken together, our findings provide insights into the biological mechanisms through which associated variants lead to altered lipid levels and potentially cardiovascular disease risk.

A framework for detecting noncoding rare-variant associations of large-scale whole-genome sequencing studies

Large-scale whole-genome sequencing studies have enabled analysis of noncoding rare-variant (RV) associations with complex human diseases and traits. Variant-set analysis is a powerful approach to study RV association. However, existing methods have limited ability in analyzing the noncoding genome. We propose a computationally efficient and robust noncoding RV association detection framework, STAARpipeline, to automatically annotate a whole-genome sequencing study and perform flexible noncoding RV association analysis, including gene-centric analysis and fixed window-based and dynamic window-based non-gene-centric analysis by incorporating variant functional annotations. In gene-centric analysis, STAARpipeline uses STAAR to group noncoding variants based on functional categories of genes and incorporate multiple functional annotations. In non-gene-centric analysis, STAARpipeline uses SCANG-STAAR to incorporate dynamic window sizes and multiple functional annotations. We apply STAARpipeline to identify noncoding RV sets associated with four lipid traits in 21,015 discovery samples from the Trans-Omics for Precision Medicine (TOPMed) program and replicate several of them in an additional 9,123 TOPMed samples. We also analyze five non-lipid TOPMed traits.

Protein biomarkers of cardiac remodeling and inflammation associated with HFpEF and incident events

There is increasing evidence that HFpEF is a heterogeneous clinical entity and distinct molecular pathways may contribute to pathophysiology. Leveraging unbiased proteomics to identify novel biomarkers, this study seeks to understand the underlying molecular mechanisms of HFpEF. The discovery cohort consisted of HFpEF cases and non-HF controls from the CATHGEN study (N = 176); the validation cohort consisted of participants from the TECOS trial of patients with diabetes (N = 109). Proteins associated with HFpEF were included in a LASSO model to create a discriminative multi-protein model and assessed in the validation cohort. Survival models and meta-analysis were used to test the association of proteins with incident clinical outcomes, including HF hospitalization, mortality and HFpEF hospitalization in CATHGEN, TECOS and the Jackson Heart Study. In the derivation set, 190 proteins were associated with HFpEF in univariate analysis, of which 65 remained significant in the multivariate model. Twenty (30.8%) of these proteins validated in TECOS, including LCN2, U-PAR, IL-1ra, KIM1, CSTB and Gal-9 (OR 1.93-2.77, p < 0.01). LASSO regression yielded a 13-protein model which, when added to a clinical model inclusive of NT-proBNP, improved the AUC from 0.82 to 0.92 (p = 1.5 × 10-4). Five proteins were associated with incident HF hospitalization, four with HFpEF hospitalization and eleven with mortality (p < 0.05). We identified and validated multiple circulating biomarkers associated with HFpEF as well as HF outcomes. These biomarkers added incremental discriminative capabilities beyond clinical factors and NT-proBNP.

Nontargeted and Targeted Metabolomic Profiling Reveals Novel Metabolite Biomarkers of Incident Diabetes in African Americans

Nontargeted metabolomics methods have increased potential to identify new disease biomarkers, but assessments of the additive information provided in large human cohorts by these less biased techniques are limited. To diversify our knowledge of diabetes-associated metabolites, we leveraged a method that measures 305 targeted or "known" and 2,342 nontargeted or "unknown" compounds in fasting plasma samples from 2,750 participants (315 incident cases) in the Jackson Heart Study (JHS)-a community cohort of self-identified African Americans-who are underrepresented in omics studies. We found 307 unique compounds (82 known) associated with diabetes after adjusting for age and sex at a false discovery rate of <0.05 and 124 compounds (35 known, including 11 not previously associated) after further adjustments for BMI and fasting plasma glucose. Of these, 144 and 68 associations, respectively, replicated in a multiethnic cohort. Among these is an apparently novel isomer of the 1-deoxyceramide Cer(m18:1/24:0) with functional geonomics and high-resolution mass spectrometry. Overall, known and unknown metabolites provided complementary information (median correlation ρ = 0.29), and their inclusion with clinical risk factors improved diabetes prediction modeling. Our findings highlight the importance of including nontargeted metabolomics methods to provide new insights into diabetes development in ethnically diverse cohorts.

Genetic analyses of the electrocardiographic QT interval and its components identify additional loci and pathways

The QT interval is an electrocardiographic measure representing the sum of ventricular depolarization and repolarization, estimated by QRS duration and JT interval, respectively. QT interval abnormalities are associated with potentially fatal ventricular arrhythmia. Using genome-wide multi-ancestry analyses (>250,000 individuals) we identify 177, 156 and 121 independent loci for QT, JT and QRS, respectively, including a male-specific X-chromosome locus. Using gene-based rare-variant methods, we identify associations with Mendelian disease genes. Enrichments are observed in established pathways for QT and JT, and previously unreported genes indicated in insulin-receptor signalling and cardiac energy metabolism. In contrast for QRS, connective tissue components and processes for cell growth and extracellular matrix interactions are significantly enriched. We demonstrate polygenic risk score associations with atrial fibrillation, conduction disease and sudden cardiac death. Prioritization of druggable genes highlight potential therapeutic targets for arrhythmia. Together, these results substantially advance our understanding of the genetic architecture of ventricular depolarization and repolarization.

Whole Genome Association Study of the Plasma Metabolome Identifies Metabolites Linked to Cardiometabolic Disease in Black Individuals

Integrating genetic information with metabolomics has provided new insights into genes affecting human metabolism. However, gene-metabolite integration has been primarily studied in individuals of European Ancestry, limiting the opportunity to leverage genomic diversity for discovery. In addition, these analyses have principally involved known metabolites, with the majority of the profiled peaks left unannotated. Here, we perform a whole genome association study of 2,291 metabolite peaks (known and unknown features) in 2,466 Black individuals from the Jackson Heart Study. We identify 519 locus-metabolite associations for 427 metabolite peaks and validate our findings in two multi-ethnic cohorts. A significant proportion of these associations are in ancestry specific alleles including findings in APOE, TTR and CD36. We leverage tandem mass spectrometry to annotate unknown metabolites, providing new insight into hereditary diseases including transthyretin amyloidosis and sickle cell disease. Our integrative omics approach leverages genomic diversity to provide novel insights into diverse cardiometabolic diseases.

Proteomic profiling platforms head to head: Leveraging genetics and clinical traits to compare aptamer- and antibody-based methods

High-throughput proteomic profiling using antibody or aptamer-based affinity reagents is used increasingly in human studies. However, direct analyses to address the relative strengths and weaknesses of these platforms are lacking. We assessed findings from the SomaScan1.3K (N = 1301 reagents), the SomaScan5K platform (N = 4979 reagents), and the Olink Explore (N = 1472 reagents) profiling techniques in 568 adults from the Jackson Heart Study and 219 participants in the HERITAGE Family Study across four performance domains: precision, accuracy, analytic breadth, and phenotypic associations leveraging detailed clinical phenotyping and genetic data. Across these studies, we show evidence supporting more reliable protein target specificity and a higher number of phenotypic associations for the Olink platform, while the Soma platforms benefit from greater measurement precision and analytic breadth across the proteome.

A multi-layer functional genomic analysis to understand noncoding genetic variation in lipids

A major challenge of genome-wide association studies (GWASs) is to translate phenotypic associations into biological insights. Here, we integrate a large GWAS on blood lipids involving 1.6 million individuals from five ancestries with a wide array of functional genomic datasets to discover regulatory mechanisms underlying lipid associations. We first prioritize lipid-associated genes with expression quantitative trait locus (eQTL) colocalizations and then add chromatin interaction data to narrow the search for functional genes. Polygenic enrichment analysis across 697 annotations from a host of tissues and cell types confirms the central role of the liver in lipid levels and highlights the selective enrichment of adipose-specific chromatin marks in high-density lipoprotein cholesterol and triglycerides. Overlapping transcription factor (TF) binding sites with lipid-associated loci identifies TFs relevant in lipid biology. In addition, we present an integrative framework to prioritize causal variants at GWAS loci, producing a comprehensive list of candidate causal genes and variants with multiple layers of functional evidence. We highlight two of the prioritized genes, CREBRF and RRBP1, which show convergent evidence across functional datasets supporting their roles in lipid biology.

Whole genome sequence association analysis of fasting glucose and fasting insulin levels in diverse cohorts from the NHLBI TOPMed program

The genetic determinants of fasting glucose (FG) and fasting insulin (FI) have been studied mostly through genome arrays, resulting in over 100 associated variants. We extended this work with high-coverage whole genome sequencing analyses from fifteen cohorts in NHLBI's Trans-Omics for Precision Medicine (TOPMed) program. Over 23,000 non-diabetic individuals from five race-ethnicities/populations (African, Asian, European, Hispanic and Samoan) were included. Eight variants were significantly associated with FG or FI across previously identified regions MTNR1B, G6PC2, GCK, GCKR and FOXA2. We additionally characterize suggestive associations with FG or FI near previously identified SLC30A8, TCF7L2, and ADCY5 regions as well as APOB, PTPRT, and ROBO1. Functional annotation resources including the Diabetes Epigenome Atlas were compiled for each signal (chromatin states, annotation principal components, and others) to elucidate variant-to-function hypotheses. We provide a catalog of nucleotide-resolution genomic variation spanning intergenic and intronic regions creating a foundation for future sequencing-based investigations of glycemic traits.

Polygenic transcriptome risk scores for COPD and lung function improve cross-ethnic portability of prediction in the NHLBI TOPMed program

While polygenic risk scores (PRSs) enable early identification of genetic risk for chronic obstructive pulmonary disease (COPD), predictive performance is limited when the discovery and target populations are not well matched. Hypothesizing that the biological mechanisms of disease are shared across ancestry groups, we introduce a PrediXcan-derived polygenic transcriptome risk score (PTRS) to improve cross-ethnic portability of risk prediction. We constructed the PTRS using summary statistics from application of PrediXcan on large-scale GWASs of lung function (forced expiratory volume in 1 s [FEV₁] and its ratio to forced vital capacity [FEV₁/FVC]) in the UK Biobank. We examined prediction performance and cross-ethnic portability of PTRS through smoking-stratified analyses both on 29,381 multi-ethnic participants from TOPMed population/family-based cohorts and on 11,771 multi-ethnic participants from TOPMed COPD-enriched studies. Analyses were carried out for two dichotomous COPD traits (moderate-to-severe and severe COPD) and two quantitative lung function traits (FEV₁ and FEV₁/FVC). While the proposed PTRS showed weaker associations with disease than PRS for European ancestry, the PTRS showed stronger association with COPD than PRS for African Americans (e.g., odds ratio [OR] = 1.24 [95% confidence interval [CI]: 1.08-1.43] for PTRS versus 1.10 [0.96-1.26] for PRS among heavy smokers with ≥ 40 pack-years of smoking) for moderate-to-severe COPD. Cross-ethnic portability of the PTRS was significantly higher than the PRS (paired t test p < 2.2 × 10^-16 with portability gains ranging from 5% to 28%) for both dichotomous COPD traits and across all smoking strata. Our study demonstrates the value of PTRS for improved cross-ethnic portability compared to PRS in predicting COPD risk.

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.

Whole Genome Sequence Analysis of the Plasma Proteome in Black Adults Provides Novel Insights Into Cardiovascular Disease

BACKGROUND

Plasma proteins are critical mediators of cardiovascular processes and are the targets of many drugs. Previous efforts to characterize the genetic architecture of the plasma proteome have been limited by a focus on individuals of European descent and leveraged genotyping arrays and imputation. Here we describe whole genome sequence analysis of the plasma proteome in individuals with greater African ancestry, increasing our power to identify novel genetic determinants.

METHODS

Proteomic profiling of 1301 proteins was performed in 1852 Black adults from the Jackson Heart Study using aptamer-based proteomics (SomaScan). Whole genome sequencing association analysis was ascertained for all variants with minor allele count ≥5. Results were validated using an alternative, antibody-based, proteomic platform (Olink) as well as replicated in the Multi-Ethnic Study of Atherosclerosis and the HERITAGE Family Study (Health, Risk Factors, Exercise Training and Genetics).

RESULTS

We identify 569 genetic associations between 479 proteins and 438 unique genetic regions at a Bonferroni-adjusted significance level of 3.8×10-11. These associations include 114 novel locus-protein relationships and an additional 217 novel sentinel variant-protein relationships. Novel cardiovascular findings include new protein associations at the APOE gene locus including ZAP70 (sentinel single nucleotide polymorphism [SNP] rs7412-T, β=0.61±0.05, P=3.27×10-30) and MMP-3 (β=-0.60±0.05, P=1.67×10-32), as well as a completely novel pleiotropic locus at the HPX gene, associated with 9 proteins. Further, the associations suggest new mechanisms of genetically mediated cardiovascular disease linked to African ancestry; we identify a novel association between variants linked to APOL1-associated chronic kidney and heart disease and the protein CKAP2 (rs73885319-G, β=0.34±0.04, P=1.34×10-17) as well as an association between ATTR amyloidosis and RBP4 levels in community-dwelling individuals without heart failure.

CONCLUSIONS

Taken together, these results provide evidence for the functional importance of variants in non-European populations, and suggest new biological mechanisms for ancestry-specific determinants of lipids, coagulation, and myocardial function.

Genetic determinants of telomere length from 109,122 ancestrally diverse whole-genome sequences in TOPMed

Genetic studies on telomere length are important for understanding age-related diseases. Prior GWAS for leukocyte TL have been limited to European and Asian populations. Here, we report the first sequencing-based association study for TL across ancestrally-diverse individuals (European, African, Asian and Hispanic/Latino) from the NHLBI Trans-Omics for Precision Medicine (TOPMed) program. We used whole genome sequencing (WGS) of whole blood for variant genotype calling and the bioinformatic estimation of telomere length in n=109,122 individuals. We identified 59 sentinel variants (p-value <5×10-9) in 36 loci associated with telomere length, including 20 newly associated loci (13 were replicated in external datasets). There was little evidence of effect size heterogeneity across populations. Fine-mapping at OBFC1 indicated the independent signals colocalized with cell-type specific eQTLs for OBFC1 (STN1). Using a multi-variant gene-based approach, we identified two genes newly implicated in telomere length, DCLRE1B (SNM1B) and PARN. In PheWAS, we demonstrated our TL polygenic trait scores (PTS) were associated with increased risk of cancer-related phenotypes.

Rare coding variants in 35 genes associate with circulating lipid levels-A multi-ancestry analysis of 170,000 exomes

Large-scale gene sequencing studies for complex traits have the potential to identify causal genes with therapeutic implications. We performed gene-based association testing of blood lipid levels with rare (minor allele frequency < 1%) predicted damaging coding variation by using sequence data from >170,000 individuals from multiple ancestries: 97,493 European, 30,025 South Asian, 16,507 African, 16,440 Hispanic/Latino, 10,420 East Asian, and 1,182 Samoan. We identified 35 genes associated with circulating lipid levels; some of these genes have not been previously associated with lipid levels when using rare coding variation from population-based samples. We prioritize 32 genes in array-based genome-wide association study (GWAS) loci based on aggregations of rare coding variants; three (EVI5, SH2B3, and PLIN1) had no prior association of rare coding variants with lipid levels. Most of our associated genes showed evidence of association among multiple ancestries. Finally, we observed an enrichment of gene-based associations for low-density lipoprotein cholesterol drug target genes and for genes closest to GWAS index single-nucleotide polymorphisms (SNPs). Our results demonstrate that gene-based associations can be beneficial for drug target development and provide evidence that the gene closest to the array-based GWAS index SNP is often the functional gene for blood lipid levels.

APOL1, Sickle Cell Trait, and CKD in the Jackson Heart Study

Rationale & Objective

Apolipoprotein L1 (APOL1) high-risk variants are associated with an increased risk for chronic kidney disease (CKD) among African Americans. Less is known regarding the risk for the development of CKD and kidney failure (end-stage kidney disease [ESKD]) among African Americans with only 1 APOL1 risk variant or whether the risk is modified by sickle cell trait.

Study Design

The Jackson Heart Study is a community-based longitudinal cohort study.

Setting & Participants

Self-reported African Americans in the Jackson Heart Study (n = 5,306).

Exposures

APOL1 G1 and G2 genotypes and sickle cell trait.

Outcomes

Incident CKD (estimated glomerular filtration rate < 60 mL/min/1.73 m²), albuminuria (urinary albumin-creatinine ratio ≥ 30 mg/g), continuous and rapid kidney function decline (≥30% decline), and incident ESKD.

Analytical Approach

Multivariable linear and logistic regression, and Cox proportional hazards models adjusted for age, sex, hypertension, diabetes, ancestry informative markers, and sickle cell trait.

Results

Of 2,300 participants, 41.3% had zero, 45.1% had 1, and 13.6% had 2 APOL1 risk variants. Sickle cell trait was present in 8.5%. Compared with participants with zero APOL1 risk variants, those with 2 alleles had an increased risk for incident albuminuria (adjusted HR [aHR], 1.88; 95% CI, 1.04 to 3.40), ESKD (aHR, 9.05; 95% CI, 1.79 to 45.85), incident CKD (aHR, 1.65; 95% CI, 1.06 to 2.57), continuous decline (β = −1.90; 95% CI, −3.35 to −0.45), and rapid kidney function decline (OR, 2.21; 95% CI, 1.22 to 4.00) after adjustment for sickle cell trait, with similar results after adjustment for ancestry informative markers. Having 1 APOL1 risk variant was not associated with CKD outcomes and there was no interaction of APOL1 with sickle cell trait.

Limitations

Single-site recruitment of African American individuals with APOL1 and sickle cell trait.

Conclusions

The presence of 1 APOL1 risk allele was not associated with increased risk for CKD outcomes, whereas 2 risk alleles were associated with incident albuminuria, CKD, ESKD, and rapid and continuous kidney function decline. Additional studies are needed to determine factors that might alter the risk for adverse kidney outcomes among individuals with high-risk APOL1 genotypes.

Associations between everyday discrimination and sleep quality and duration among African-Americans over time in the Jackson Heart Study

STUDY OBJECTIVES

African-Americans have a high burden of poor sleep, yet, psychosocial determinants (e.g. discrimination) are understudied. We investigated longitudinal associations between everyday discrimination and sleep quality and duration among African-Americans (N = 3404) in the Jackson Heart Study.

METHODS

At Exam 1 (2000-2004) and Exam 3 (2008-2013), participants completed the Everyday Discrimination Scale, rated their sleep quality (1 = poor to 5 = excellent), and self-reported hours of sleep. A subset of participants (N = 762) underwent 7-day actigraphy to objectively measure sleep duration and sleep quality (Sleep Exam 2012-2016). Changes in discrimination were defined as low stable (reference), increasing, decreasing, and high stable. Within-person changes in sleep from Exam 1 to Exam 3 were regressed on change in discrimination from Exam 1 to Exam 3 while adjusting for age, sex, education, income, employment, physical activity, smoking, body mass index, social support, and stress.

RESULTS

At Exam 1, the mean age was 54.1 (12.0) years; 64% were female, mean sleep quality was 3.0 (1.1) and 54% were short sleepers. The distribution of the discrimination change trajectories were 54.1% low stable, 13.5% increasing, 14.6% decreasing, and 17.7% were high stable. Participants who were in the increasing (vs. low stable) discrimination group had greater decrease in sleep quality. There was no association between change in discrimination and change in sleep duration. Among Sleep Exam participants, higher discrimination was cross-sectionally associated with shorter self-reported sleep duration, independent of stress.

CONCLUSION

Discrimination is a unique stressor for African-Americans; thus, future research should identify interventions to reduce the burden of discrimination on sleep quality.

Metabolomic Analysis of Coronary Heart Disease in an African American Cohort From the Jackson Heart Study

Importance

African American individuals have disproportionate rates of coronary heart disease (CHD) but lower levels of coronary artery calcium (CAC), a marker of subclinical CHD, than non-Hispanic White individuals. African American individuals may have distinct metabolite profiles associated with incident CHD risk compared with non-Hispanic White individuals, and examination of these differences could highlight important processes that differ between them.

Objectives

To identify novel biomarkers of incident CHD and CAC among African American individuals and to replicate incident CHD findings in a multiethnic cohort.

Design, Setting, and Participants

This analysis targeted plasma metabolomic profiling of 2346 participants in the Jackson Heart Study (JHS), a prospective population-based cohort study that included 5306 African American participants who were examined at baseline (2000-2004) and 2 follow-up visits. Replication of CHD-associated metabolites was sought among 1588 multiethnic participants from the Women's Health Initiative (WHI), a prospective population-based multiethnic cohort study of 161 808 postmenopausal women who were examined at baseline (1991-1995) and ongoing follow-up visits. Regression analyses were performed for each metabolite to examine the associations with incident CHD and CAC scores. Data were collected from the WHI between 1994 and 2009 and from the JHS between 2000 and 2015. All data were analyzed from November 2020 to August 2021.

Exposures

Plasma metabolites.

Main Outcomes and Measures

Incident CHD was defined as definite or probable myocardial infarction or definite fatal CHD in both the JHS and WHI cohorts. In the JHS cohort, silent myocardial infarction between examinations (as determined by electrocardiography) and coronary revascularization were included in the incident CHD analysis. Coronary artery calcium was measured using a 16-channel computed tomographic system and reported as an Agatston score.

Results

Among 2346 African American individuals in the JHS cohort, the mean (SD) age was 56 (13) years, and 1468 individuals (62.6%) were female. Among 1588 postmenopausal women in the WHI cohort, the mean (SD) age was 67 (7) years; 217 individuals (13.7%) self-identified as African American, 1219 (76.8%) as non-Hispanic White, and 152 (9.6%) as other races or ethnicities. In the fully adjusted model including 1876 individuals, 46 of 303 targeted metabolites were associated with incident CHD (false discovery rate q <0.100). Data for 32 of the 46 metabolites were available in the WHI cohort, and 13 incident CHD-associated metabolites from the JHS cohort were replicated in the WHI cohort. A total of 1439 participants from the JHS cohort with available CAC scores received metabolomic profiling. Nine metabolites were associated with CAC scores. Minimal overlap was found between the results from the incident CHD and CAC analyses, with only 3 metabolites shared between the 2 analyses.

Conclusions and Relevance

This cohort study identified metabolites that were associated with incident CHD among African American individuals, including 13 incident CHD-associated metabolites that were replicated in a multiethnic population and 9 novel metabolites that included N-acylamides, leucine, and lipid species. These findings may help to elucidate common and distinct metabolic processes that may be associated with CHD among individuals with different self-identified race.

The power of genetic diversity in genome-wide association studies of lipids

Increased blood lipid levels are heritable risk factors of cardiovascular disease with varied prevalence worldwide owing to different dietary patterns and medication use1. Despite advances in prevention and treatment, in particular through reducing low-density lipoprotein cholesterol levels2, heart disease remains the leading cause of death worldwide3. Genome-wideassociation studies (GWAS) of blood lipid levels have led to important biological and clinical insights, as well as new drug targets, for cardiovascular disease. However, most previous GWAS4-23 have been conducted in European ancestry populations and may have missed genetic variants that contribute to lipid-level variation in other ancestry groups. These include differences in allele frequencies, effect sizes and linkage-disequilibrium patterns24. Here we conduct a multi-ancestry, genome-wide genetic discovery meta-analysis of lipid levels in approximately 1.65 million individuals, including 350,000 of non-European ancestries. We quantify the gain in studying non-European ancestries and provide evidence to support the expansion of recruitment of additional ancestries, even with relatively small sample sizes. We find that increasing diversity rather than studying additional individuals of European ancestry results in substantial improvements in fine-mapping functional variants and portability of polygenic prediction (evaluated in approximately 295,000 individuals from 7 ancestry groupings). Modest gains in the number of discovered loci and ancestry-specific variants were also achieved. As GWAS expand emphasis beyond the identification of genes and fundamental biology towards the use of genetic variants for preventive and precision medicine25, we anticipate that increased diversity of participants will lead to more accurate and equitable26 application of polygenic scores in clinical practice.

APOL1 renal risk variants are associated with obesity and body composition in African ancestry adults: An observational genotype-phenotype association study

ABSTRACT

While increased obesity prevalence among persons of African ancestry (AAs) compared to persons of European ancestry (EAs) is linked to social, environmental and behavioral factors, there are no gene variants that are common and significantly associated with obesity in AA populations. We sought to explore the association between ancestry specific renal risk variants in the apolipoprotein L1 (APOL1) gene with obesity related traits in AAs.We conducted a genotype-phenotype association study from 3 electronic medical record linked cohorts (BioMe Biobank, BioVU, nuGENE); randomized controlled trials (genetic testing to understand and address renal disease disparities) and prospective cohort study (Jackson Heart Study). We analyzed association of APOL1 renal risk variants with cross-sectional measures of obesity (average body mass index (BMI), and proportion of overweight and obesity) and with measures of body composition (in Jackson Heart Study).We had data on 11,930 self-reported AA adults. Across cohorts, mean age was from 42 to 49 years and percentage female from 58% to 75.3%. Individuals who have 2 APOL1 risk alleles (14% of AAs) have 30% higher obesity odds compared to others (recessive model adjusted odds ratio 1.30; 95% confidence interval 1.16-1.41; P = 2.75 × 10-6). An additive model better fit the association, in which each allele (47% of AAs) increases obesity odds by 1.13-fold (adjusted odds ratio 1.13; 95% confidence interval 1.07-1.19; P = 3.07 × 10-6) and increases BMI by 0.36 kg/m2 (∼1 kg, for 1.7 m height; P = 2 × 10-4). APOL1 alleles are not associated with refined body composition traits overall but are significantly associated with fat free mass index in women [0.30 kg/m2 increment per allele; P = .03].Thus, renal risk variants in the APOL1 gene, found in nearly half of AAs, are associated with BMI and obesity in an additive manner. These variants could, either on their own or interacting with environmental factors, explain a proportion of ethnic disparities in obesity.

A high-resolution HLA reference panel capturing global population diversity enables multi-ancestry fine-mapping in HIV host response

Fine-mapping to plausible causal variation may be more effective in multi-ancestry cohorts, particularly in the MHC, which has population-specific structure. To enable such studies, we constructed a large (n = 21,546) HLA reference panel spanning five global populations based on whole-genome sequences. Despite population-specific long-range haplotypes, we demonstrated accurate imputation at G-group resolution (94.2%, 93.7%, 97.8% and 93.7% in admixed African (AA), East Asian (EAS), European (EUR) and Latino (LAT) populations). Applying HLA imputation to genome-wide association study data for HIV-1 viral load in three populations (EUR, AA and LAT), we obviated effects of previously reported associations from population-specific HIV studies and discovered a novel association at position 156 in HLA-B. We pinpointed the MHC association to three amino acid positions (97, 67 and 156) marking three consecutive pockets (C, B and D) within the HLA-B peptide-binding groove, explaining 12.9% of trait variance.

Presence and transmission of mitochondrial heteroplasmic mutations in human populations of European and African ancestry

We investigated the concordance of mitochondrial DNA heteroplasmic mutations (heteroplasmies) in 6745 maternal pairs of European (EA, n = 4718 pairs) and African (AA, n = 2027 pairs) Americans in whole blood. Mother-offspring pairs displayed the highest concordance rate, followed by sibling-sibling and more distantly-related maternal pairs. The allele fractions of concordant heteroplasmies exhibited high correlation (R2 = 0.8) between paired individuals. Discordant heteroplasmies were more likely to be in coding regions, be nonsynonymous or nonsynonymous-deleterious (p < 0.001). The number of deleterious heteroplasmies was significantly correlated with advancing age (20-44, 45-64, and ≥65 years, p-trend = 0.01). One standard deviation increase in heteroplasmic burden (i.e., the number of heteroplasmies carried by an individual) was associated with 0.17 to 0.26 (p < 1e - 23) standard deviation decrease in mtDNA copy number, independent of age. White blood cell count and differential count jointly explained 0.5% to 1.3% (p ≤ 0.001) variance in heteroplasmic burden. A genome-wide association and meta-analysis identified a region at 11p11.12 (top signal rs779031139, p = 2.0e - 18, minor allele frequency = 0.38) associated with the heteroplasmic burden. However, the 11p11.12 region is adjacent to a nuclear mitochondrial DNA (NUMT) corresponding to a 542 bp area of the D-loop. This region was no longer significant after excluding heteroplasmies within the 542 bp from the heteroplasmic burden. The discovery that blood mtDNA heteroplasmies were both inherited and somatic origins and that an increase in heteroplasmic burden was strongly associated with a decrease in average number of mtDNA copy number in blood are important findings to be considered in association studies of mtDNA with disease traits.

Whole-genome association analyses of sleep-disordered breathing phenotypes in the NHLBI TOPMed program

BACKGROUND

Sleep-disordered breathing is a common disorder associated with significant morbidity. The genetic architecture of sleep-disordered breathing remains poorly understood. Through the NHLBI Trans-Omics for Precision Medicine (TOPMed) program, we performed the first whole-genome sequence analysis of sleep-disordered breathing.

METHODS

The study sample was comprised of 7988 individuals of diverse ancestry. Common-variant and pathway analyses included an additional 13,257 individuals. We examined five complementary traits describing different aspects of sleep-disordered breathing: the apnea-hypopnea index, average oxyhemoglobin desaturation per event, average and minimum oxyhemoglobin saturation across the sleep episode, and the percentage of sleep with oxyhemoglobin saturation < 90%. We adjusted for age, sex, BMI, study, and family structure using MMSKAT and EMMAX mixed linear model approaches. Additional bioinformatics analyses were performed with MetaXcan, GIGSEA, and ReMap.

RESULTS

We identified a multi-ethnic set-based rare-variant association (p = 3.48 × 10-8) on chromosome X with ARMCX3. Additional rare-variant associations include ARMCX3-AS1, MRPS33, and C16orf90. Novel common-variant loci were identified in the NRG1 and SLC45A2 regions, and previously associated loci in the IL18RAP and ATP2B4 regions were associated with novel phenotypes. Transcription factor binding site enrichment identified associations with genes implicated with respiratory and craniofacial traits. Additional analyses identified significantly associated pathways.

CONCLUSIONS

We have identified the first gene-based rare-variant associations with objectively measured sleep-disordered breathing traits. Our results increase the understanding of the genetic architecture of sleep-disordered breathing and highlight associations in genes that modulate lung development, inflammation, respiratory rhythmogenesis, and HIF1A-mediated hypoxic response.

Multiethnic Genome-Wide Association Study of Subclinical Atherosclerosis in Individuals With Type 2 Diabetes

BACKGROUND

Coronary artery calcification (CAC) and carotid artery intima-media thickness (cIMT) are measures of subclinical atherosclerosis in asymptomatic individuals and strong risk factors for cardiovascular disease. Type 2 diabetes (T2D) is an independent cardiovascular disease risk factor that accelerates atherosclerosis.

METHODS

We performed meta-analyses of genome-wide association studies in up to 2500 T2D individuals of European ancestry (EA) and 1590 T2D individuals of African ancestry with or without exclusion of prevalent cardiovascular disease, for CAC measured by cardiac computed tomography, and 3608 individuals of EA and 838 individuals of African ancestry with T2D for cIMT measured by ultrasonography within the CHARGE (Cohorts for Heart and Aging Research in Genomic Epidemiology) Consortium.

RESULTS

We replicated 2 loci (rs9369640 and rs9349379 near PHACTR1 and rs10757278 near CDKN2B) for CAC and one locus for cIMT (rs7412 and rs445925 near APOE-APOC1) that were previously reported in the general EA populations. We identified one novel CAC locus (rs8000449 near CSNK1A1L/LINC00547/POSTN at 13q13.3) at P=2.0×10-8 in EA. No additional loci were identified with the meta-analyses of EA and African ancestry. The expression quantitative trait loci analysis with nearby expressed genes derived from arterial wall and metabolic tissues from the Genotype-Tissue Expression project pinpoints POSTN, encoding a matricellular protein involved in bone formation and bone matrix organization, as the potential candidate gene at this locus. In addition, we found significant associations (P<3.1×10-4) for 3 previously reported coronary artery disease loci for these subclinical atherosclerotic phenotypes (rs2891168 near CDKN2B-AS1 and rs11170820 near FLJ12825 for CAC, and rs7412 near APOE for cIMT).

CONCLUSIONS

Our results provide potential biological mechanisms that could link CAC and cIMT to increased cardiovascular disease risk in individuals with T2D.

Genome-wide association studies identify 137 genetic loci for DNA methylation biomarkers of aging

BACKGROUND

Biological aging estimators derived from DNA methylation data are heritable and correlate with morbidity and mortality. Consequently, identification of genetic and environmental contributors to the variation in these measures in populations has become a major goal in the field.

RESULTS

Leveraging DNA methylation and SNP data from more than 40,000 individuals, we identify 137 genome-wide significant loci, of which 113 are novel, from genome-wide association study (GWAS) meta-analyses of four epigenetic clocks and epigenetic surrogate markers for granulocyte proportions and plasminogen activator inhibitor 1 levels, respectively. We find evidence for shared genetic loci associated with the Horvath clock and expression of transcripts encoding genes linked to lipid metabolism and immune function. Notably, these loci are independent of those reported to regulate DNA methylation levels at constituent clock CpGs. A polygenic score for GrimAge acceleration showed strong associations with adiposity-related traits, educational attainment, parental longevity, and C-reactive protein levels.

CONCLUSION

This study illuminates the genetic architecture underlying epigenetic aging and its shared genetic contributions with lifestyle factors and longevity.

Determinants of penetrance and variable expressivity in monogenic metabolic conditions across 77,184 exomes

Hundreds of thousands of genetic variants have been reported to cause severe monogenic diseases, but the probability that a variant carrier develops the disease (termed penetrance) is unknown for virtually all of them. Additionally, the clinical utility of common polygenetic variation remains uncertain. Using exome sequencing from 77,184 adult individuals (38,618 multi-ancestral individuals from a type 2 diabetes case-control study and 38,566 participants from the UK Biobank, for whom genotype array data were also available), we apply clinical standard-of-care gene variant curation for eight monogenic metabolic conditions. Rare variants causing monogenic diabetes and dyslipidemias display effect sizes significantly larger than the top 1% of the corresponding polygenic scores. Nevertheless, penetrance estimates for monogenic variant carriers average 60% or lower for most conditions. We assess epidemiologic and genetic factors contributing to risk prediction in monogenic variant carriers, demonstrating that inclusion of polygenic variation significantly improves biomarker estimation for two monogenic dyslipidemias.

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 < 8.6 × 10⁻⁷) to 3.08 years (EEAA, p < 3.7 × 10⁻¹⁸). 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 >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 < 4.1 × 10⁻⁸) and coronary heart disease (CHD) (hazard ratio 3.24, p < 9.3 × 10⁻⁶) 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.

Multiomic Profiling in Black and White Populations Reveals Novel Candidate Pathways in Left Ventricular Hypertrophy and Incident Heart Failure Specific to Black Adults

BACKGROUND

Increased left ventricular (LV) mass is associated with adverse cardiovascular events including heart failure (HF). Both increased LV mass and HF disproportionately affect Black individuals. To understand the underlying mechanisms, we undertook a proteomic screen in a Black cohort and compared the findings to results from a White cohort.

METHODS

We measured 1305 plasma proteins using the SomaScan platform in 1772 Black participants (mean age, 56 years; 62% women) in JHS (Jackson Heart Study) with LV mass assessed by 2-dimensional echocardiography. Incident HF was assessed in 1600 participants. We then compared protein associations in JHS to those observed in White participants from FHS (Framingham Heart Study; mean age, 54 years; 56% women).

RESULTS

In JHS, there were 110 proteins associated with LV mass and 13 proteins associated with incident HF hospitalization with false discovery rate <5% after multivariable adjustment. Several proteins showed expected associations with both LV mass and HF, including NT-proBNP (N-terminal pro-B-type natriuretic peptide; β=0.04; P=2×10-8; hazard ratio, 1.48; P=0.0001). The strongest association with LV mass was novel: LKHA4 (leukotriene-A4 hydrolase; β=0.05; P=5×10-15). This association was confirmed on an alternate proteomics platform and further supported by related metabolomic data. Fractalkine/CX3CL1 (C-X3-C Motif Chemokine Ligand 1) showed a novel association with incident HF (hazard ratio, 1.32; P=0.0002). While established biomarkers such as cystatin C and NT-proBNP showed consistent associations in Black and White individuals, LKHA4 and fractalkine were significantly different between the two groups.

CONCLUSIONS

We identified several novel biological pathways specific to Black adults hypothesized to contribute to the pathophysiologic cascade of LV hypertrophy and incident HF including LKHA4 and fractalkine.

Chromosome Xq23 is associated with lower atherogenic lipid concentrations and favorable cardiometabolic indices

Autosomal genetic analyses of blood lipids have yielded key insights for coronary heart disease (CHD). However, X chromosome genetic variation is understudied for blood lipids in large sample sizes. We now analyze genetic and blood lipid data in a high-coverage whole X chromosome sequencing study of 65,322 multi-ancestry participants and perform replication among 456,893 European participants. Common alleles on chromosome Xq23 are strongly associated with reduced total cholesterol, LDL cholesterol, and triglycerides (min P = 8.5 × 10-72), with similar effects for males and females. Chromosome Xq23 lipid-lowering alleles are associated with reduced odds for CHD among 42,545 cases and 591,247 controls (P = 1.7 × 10-4), and reduced odds for diabetes mellitus type 2 among 54,095 cases and 573,885 controls (P = 1.4 × 10-5). Although we observe an association with increased BMI, waist-to-hip ratio adjusted for BMI is reduced, bioimpedance analyses indicate increased gluteofemoral fat, and abdominal MRI analyses indicate reduced visceral adiposity. Co-localization analyses strongly correlate increased CHRDL1 gene expression, particularly in adipose tissue, with reduced concentrations of blood lipids.

Discovery and fine-mapping of height loci via high-density imputation of GWASs in individuals of African ancestry

Although many loci have been associated with height in European ancestry populations, very few have been identified in African ancestry individuals. Furthermore, many of the known loci have yet to be generalized to and fine-mapped within a large-scale African ancestry sample. We performed sex-combined and sex-stratified meta-analyses in up to 52,764 individuals with height and genome-wide genotyping data from the African Ancestry Anthropometry Genetics Consortium (AAAGC). We additionally combined our African ancestry meta-analysis results with published European genome-wide association study (GWAS) data. In the African ancestry analyses, we identified three novel loci (SLC4A3, NCOA2, ECD/FAM149B1) in sex-combined results and two loci (CRB1, KLF6) in women only. In the African plus European sex-combined GWAS, we identified an additional three novel loci (RCCD1, G6PC3, CEP95) which were equally driven by AAAGC and European results. Among 39 genome-wide significant signals at known loci, conditioning index SNPs from European studies identified 20 secondary signals. Two of the 20 new secondary signals and none of the 8 novel loci had minor allele frequencies (MAF) < 5%. Of 802 known European height signals, 643 displayed directionally consistent associations with height, of which 205 were nominally significant (p < 0.05) in the African ancestry sex-combined sample. Furthermore, 148 of 241 loci contained ≤20 variants in the credible sets that jointly account for 99% of the posterior probability of driving the associations. In summary, trans-ethnic meta-analyses revealed novel signals and further improved fine-mapping of putative causal variants in loci shared between African and European ancestry populations.