A multi-cohort genome-wide association study in African ancestry individuals reveals risk loci for primary open-angle glaucoma.

Primary open-angle glaucoma (POAG), the leading cause of irreversible blindness worldwide, disproportionately affects individuals of African ancestry. We conducted a genome-wide association study (GWAS) for POAG in 11,275 individuals of African ancestry (6,003 cases; 5,272 controls). We detected 46 risk loci associated with POAG at genome-wide significance. Replication and post-GWAS analyses, including functionally informed fine-mapping, multiple trait co-localization, and in silico validation, implicated two previously undescribed variants (rs1666698 mapping to DBF4P2; rs34957764 mapping to ROCK1P1) and one previously associated variant (rs11824032 mapping to ARHGEF12) as likely causal. For individuals of African ancestry, a polygenic risk score (PRS) for POAG from our mega-analysis (African ancestry individuals) outperformed a PRS from summary statistics of a much larger GWAS derived from European ancestry individuals. This study quantifies the genetic architecture similarities and differences between African and non-African ancestry populations for this blinding disease.

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.

Aberrant activation of TCL1A promotes stem cell expansion in clonal haematopoiesis.

Mutations in a diverse set of driver genes increase the fitness of haematopoietic stem cells (HSCs), leading to clonal haematopoiesis1. These lesions are precursors for blood cancers2-6, but the basis of their fitness advantage remains largely unknown, partly owing to a paucity of large cohorts in which the clonal expansion rate has been assessed by longitudinal sampling. Here, to circumvent this limitation, we developed a method to infer the expansion rate from data from a single time point. We applied this method to 5,071 people with clonal haematopoiesis. A genome-wide association study revealed that a common inherited polymorphism in the TCL1A promoter was associated with a slower expansion rate in clonal haematopoiesis overall, but the effect varied by driver gene. Those carrying this protective allele exhibited markedly reduced growth rates or prevalence of clones with driver mutations in TET2, ASXL1, SF3B1 and SRSF2, but this effect was not seen in clones with driver mutations in DNMT3A. TCL1A was not expressed in normal or DNMT3A-mutated HSCs, but the introduction of mutations in TET2 or ASXL1 led to the expression of TCL1A protein and the expansion of HSCs in vitro. The protective allele restricted TCL1A expression and expansion of mutant HSCs, as did experimental knockdown of TCL1A expression. Forced expression of TCL1A promoted the expansion of human HSCs in vitro and mouse HSCs in vivo. Our results indicate that the fitness advantage of several commonly mutated driver genes in clonal haematopoiesis may be mediated by TCL1A activation.

Rare variants in long non-coding RNAs are associated with blood lipid levels in the TOPMed whole-genome sequencing study.

Long non-coding RNAs (lncRNAs) are known to perform important regulatory functions in lipid metabolism. Large-scale whole-genome sequencing (WGS) studies and new statistical methods for variant set tests now provide an opportunity to assess more associations between rare variants in lncRNA genes and complex traits across the genome. In this study, we used high-coverage WGS from 66,329 participants of diverse ancestries with measurement of blood lipids and lipoproteins (LDL-C, HDL-C, TC, and TG) in the National Heart, Lung, and Blood Institute (NHLBI) Trans-Omics for Precision Medicine (TOPMed) program to investigate the role of lncRNAs in lipid variability. We aggregated rare variants for 165,375 lncRNA genes based on their genomic locations and conducted rare-variant aggregate association tests using the STAAR (variant-set test for association using annotation information) framework. We performed STAAR conditional analysis adjusting for common variants in known lipid GWAS loci and rare-coding variants in nearby protein-coding genes. Our analyses revealed 83 rare lncRNA variant sets significantly associated with blood lipid levels, all of which were located in known lipid GWAS loci (in a ±500-kb window of a Global Lipids Genetics Consortium index variant). Notably, 61 out of 83 signals (73%) were conditionally independent of common regulatory variation and rare protein-coding variation at the same loci. We replicated 34 out of 61 (56%) conditionally independent associations using the independent UK Biobank WGS data. Our results expand the genetic architecture of blood lipids to rare variants in lncRNAs.

A genetic association study of circulating coagulation factor VIII and von Willebrand factor levels.

ABSTRACT: Coagulation factor VIII (FVIII) and its carrier protein von Willebrand factor (VWF) are critical to coagulation and platelet aggregation. We leveraged whole-genome sequence data from the Trans-Omics for Precision Medicine (TOPMed) program along with TOPMed-based imputation of genotypes in additional samples to identify genetic associations with circulating FVIII and VWF levels in a single-variant meta-analysis, including up to 45 289 participants. Gene-based aggregate tests were implemented in TOPMed. We identified 3 candidate causal genes and tested their functional effect on FVIII release from human liver endothelial cells (HLECs) and VWF release from human umbilical vein endothelial cells. Mendelian randomization was also performed to provide evidence for causal associations of FVIII and VWF with thrombotic outcomes. We identified associations (P < 5 × 10-9) at 7 new loci for FVIII (ST3GAL4, CLEC4M, B3GNT2, ASGR1, F12, KNG1, and TREM1/NCR2) and 1 for VWF (B3GNT2). VWF, ABO, and STAB2 were associated with FVIII and VWF in gene-based analyses. Multiphenotype analysis of FVIII and VWF identified another 3 new loci, including PDIA3. Silencing of B3GNT2 and the previously reported CD36 gene decreased release of FVIII by HLECs, whereas silencing of B3GNT2, CD36, and PDIA3 decreased release of VWF by HVECs. Mendelian randomization supports causal association of higher FVIII and VWF with increased risk of thrombotic outcomes. Seven new loci were identified for FVIII and 1 for VWF, with evidence supporting causal associations of FVIII and VWF with thrombotic outcomes. B3GNT2, CD36, and PDIA3 modulate the release of FVIII and/or VWF in vitro.

Inherited causes of clonal haematopoiesis in 97,691 whole genomes.

Age is the dominant risk factor for most chronic human diseases, but the mechanisms through which ageing confers this risk are largely unknown1. The age-related acquisition of somatic mutations that lead to clonal expansion in regenerating haematopoietic stem cell populations has recently been associated with both haematological cancer2-4 and coronary heart disease5-this phenomenon is termed clonal haematopoiesis of indeterminate potential (CHIP)6. Simultaneous analyses of germline and somatic whole-genome sequences provide the opportunity to identify root causes of CHIP. Here we analyse high-coverage whole-genome sequences from 97,691 participants of diverse ancestries in the National Heart, Lung, and Blood Institute Trans-omics for Precision Medicine (TOPMed) programme, and identify 4,229 individuals with CHIP. We identify associations with blood cell, lipid and inflammatory traits that are specific to different CHIP driver genes. Association of a genome-wide set of germline genetic variants enabled the identification of three genetic loci associated with CHIP status, including one locus at TET2 that was specific to individuals of African ancestry. In silico-informed in vitro evaluation of the TET2 germline locus enabled the identification of a causal variant that disrupts a TET2 distal enhancer, resulting in increased self-renewal of haematopoietic stem cells. Overall, we observe that germline genetic variation shapes haematopoietic stem cell function, leading to CHIP through mechanisms that are specific to clonal haematopoiesis as well as shared mechanisms that lead to somatic mutations across tissues.

Identification of type 2 diabetes loci in 433,540 East Asian individuals.

Meta-analyses of genome-wide association studies (GWAS) have identified more than 240 loci that are associated with type 2 diabetes (T2D)1,2; however, most of these loci have been identified in analyses of individuals with European ancestry. Here, to examine T2D risk in East Asian individuals, we carried out a meta-analysis of GWAS data from 77,418 individuals with T2D and 356,122 healthy control individuals. In the main analysis, we identified 301 distinct association signals at 183 loci, and across T2D association models with and without consideration of body mass index and sex, we identified 61 loci that are newly implicated in predisposition to T2D. Common variants associated with T2D in both East Asian and European populations exhibited strongly correlated effect sizes. Previously undescribed associations include signals in or near GDAP1, PTF1A, SIX3, ALDH2, a microRNA cluster, and genes that affect the differentiation of muscle and adipose cells3. At another locus, expression quantitative trait loci at two overlapping T2D signals affect two genes-NKX6-3 and ANK1-in different tissues4-6. Association studies in diverse populations identify additional loci and elucidate disease-associated genes, biology, and pathways.

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.

Whole-exome sequencing study identifies four novel gene loci associated with diabetic kidney disease.

Diabetic kidney disease (DKD) is recognized as an important public health challenge. However, its genomic mechanisms are poorly understood. To identify rare variants for DKD, we conducted a whole-exome sequencing (WES) study leveraging large cohorts well-phenotyped for chronic kidney disease and diabetes. Our two-stage WES study included 4372 European and African ancestry participants from the Chronic Renal Insufficiency Cohort and Atherosclerosis Risk in Communities studies (stage 1) and 11 487 multi-ancestry Trans-Omics for Precision Medicine participants (stage 2). Generalized linear mixed models, which accounted for genetic relatedness and adjusted for age, sex and ancestry, were used to test associations between single variants and DKD. Gene-based aggregate rare variant analyses were conducted using an optimized sequence kernel association test implemented within our mixed model framework. We identified four novel exome-wide significant DKD-related loci through initiating diabetes. In single-variant analyses, participants carrying a rare, in-frame insertion in the DIS3L2 gene (rs141560952) exhibited a 193-fold increased odds [95% confidence interval (CI): 33.6, 1105] of DKD compared with noncarriers (P = 3.59 × 10-9). Likewise, each copy of a low-frequency KRT6B splice-site variant (rs425827) conferred a 5.31-fold higher odds (95% CI: 3.06, 9.21) of DKD (P = 2.72 × 10-9). Aggregate gene-based analyses further identified ERAP2 (P = 4.03 × 10-8) and NPEPPS (P = 1.51 × 10-7), which are both expressed in the kidney and implicated in renin-angiotensin-aldosterone system modulated immune response. In the largest WES study of DKD, we identified novel rare variant loci attaining exome-wide significance. These findings provide new insights into the molecular mechanisms underlying DKD.

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.

Low MBOAT7 expression, a genetic risk for MASH, promotes a profibrotic pathway involving hepatocyte TAZ upregulation.

BACKGROUND AND AIMS: The common genetic variant rs641738 C>T is a risk factor for metabolic dysfunction-associated steatotic liver disease and metabolic dysfunction-associated steatohepatitis (MASH), including liver fibrosis, and is associated with decreased expression of the phospholipid-remodeling enzyme MBOAT7 (LPIAT1). However, whether restoring MBOAT7 expression in established metabolic dysfunction-associated steatotic liver disease dampens the progression to liver fibrosis and, importantly, the mechanism through which decreased MBOAT7 expression exacerbates MASH fibrosis remain unclear. APPROACH AND RESULTS: We first showed that hepatocyte MBOAT7 restoration in mice with diet-induced steatohepatitis slows the progression to liver fibrosis. Conversely, when hepatocyte-MBOAT7 was silenced in mice with established hepatosteatosis, liver fibrosis but not hepatosteatosis was exacerbated. Mechanistic studies revealed that hepatocyte-MBOAT7 restoration in MASH mice lowered hepatocyte-TAZ (WWTR1), which is known to promote MASH fibrosis. Conversely, hepatocyte-MBOAT7 silencing enhanced TAZ upregulation in MASH. Finally, we discovered that changes in hepatocyte phospholipids due to MBOAT7 loss-of-function promote a cholesterol trafficking pathway that upregulates TAZ and the TAZ-induced profibrotic factor Indian hedgehog (IHH). As evidence for relevance in humans, we found that the livers of individuals with MASH carrying the rs641738-T allele had higher hepatocyte nuclear TAZ, indicating higher TAZ activity and increased IHH mRNA. CONCLUSIONS: This study provides evidence for a novel mechanism linking MBOAT7-LoF to MASH fibrosis, adds new insight into an established genetic locus for MASH, and, given the druggability of hepatocyte TAZ for MASH fibrosis, suggests a personalized medicine approach for subjects at increased risk for MASH fibrosis due to inheritance of variants that lower MBOAT7.

Neurocognitive trajectory and proteomic signature of inherited risk for Alzheimer's disease.

For Alzheimer's disease-a leading cause of dementia and global morbidity-improved identification of presymptomatic high-risk individuals and identification of new circulating biomarkers are key public health needs. Here, we tested the hypothesis that a polygenic predictor of risk for Alzheimer's disease would identify a subset of the population with increased risk of clinically diagnosed dementia, subclinical neurocognitive dysfunction, and a differing circulating proteomic profile. Using summary association statistics from a recent genome-wide association study, we first developed a polygenic predictor of Alzheimer's disease comprised of 7.1 million common DNA variants. We noted a 7.3-fold (95% CI 4.8 to 11.0; p < 0.001) gradient in risk across deciles of the score among 288,289 middle-aged participants of the UK Biobank study. In cross-sectional analyses stratified by age, minimal differences in risk of Alzheimer's disease and performance on a digit recall test were present according to polygenic score decile at age 50 years, but significant gradients emerged by age 65. Similarly, among 30,541 participants of the Mass General Brigham Biobank, we again noted no significant differences in Alzheimer's disease diagnosis at younger ages across deciles of the score, but for those over 65 years we noted an odds ratio of 2.0 (95% CI 1.3 to 3.2; p = 0.002) in the top versus bottom decile of the polygenic score. To understand the proteomic signature of inherited risk, we performed aptamer-based profiling in 636 blood donors (mean age 43 years) with very high or low polygenic scores. In addition to the well-known apolipoprotein E biomarker, this analysis identified 27 additional proteins, several of which have known roles related to disease pathogenesis. Differences in protein concentrations were consistent even among the youngest subset of blood donors (mean age 33 years). Of these 28 proteins, 7 of the 8 proteins with concentrations available were similarly associated with the polygenic score in participants of the Multi-Ethnic Study of Atherosclerosis. These data highlight the potential for a DNA-based score to identify high-risk individuals during the prolonged presymptomatic phase of Alzheimer's disease and to enable biomarker discovery based on profiling of young individuals in the extremes of the score distribution.

Genetic and clinical determinants of abdominal aortic diameter: genome-wide association studies, exome array data and Mendelian randomization study.

Progressive dilation of the infrarenal aortic diameter is a consequence of the ageing process and is considered the main determinant of abdominal aortic aneurysm (AAA). We aimed to investigate the genetic and clinical determinants of abdominal aortic diameter (AAD). We conducted a meta-analysis of genome-wide association studies in 10 cohorts (n = 13 542) imputed to the 1000 Genome Project reference panel including 12 815 subjects in the discovery phase and 727 subjects [Partners Biobank cohort 1 (PBIO)] as replication. Maximum anterior-posterior diameter of the infrarenal aorta was used as AAD. We also included exome array data (n = 14 480) from seven epidemiologic studies. Single-variant and gene-based associations were done using SeqMeta package. A Mendelian randomization analysis was applied to investigate the causal effect of a number of clinical risk factors on AAD. In genome-wide association study (GWAS) on AAD, rs74448815 in the intronic region of LDLRAD4 reached genome-wide significance (beta = -0.02, SE = 0.004, P-value = 2.10 × 10-8). The association replicated in the PBIO1 cohort (P-value = 8.19 × 10-4). In exome-array single-variant analysis (P-value threshold = 9 × 10-7), the lowest P-value was found for rs239259 located in SLC22A20 (beta = 0.007, P-value = 1.2 × 10-5). In the gene-based analysis (P-value threshold = 1.85 × 10-6), PCSK5 showed an association with AAD (P-value = 8.03 × 10-7). Furthermore, in Mendelian randomization analyses, we found evidence for genetic association of pulse pressure (beta = -0.003, P-value = 0.02), triglycerides (beta = -0.16, P-value = 0.008) and height (beta = 0.03, P-value < 0.0001), known risk factors for AAA, consistent with a causal association with AAD. Our findings point to new biology as well as highlighting gene regions in mechanisms that have previously been implicated in the genetics of other vascular diseases.

Differences in Metabolomic Profiles Between Black and White Women and Risk of Coronary Heart Disease: an Observational Study of Women From Four US Cohorts.

BACKGROUND: Racial differences in metabolomic profiles may reflect underlying differences in social determinants of health by self-reported race and may be related to racial disparities in coronary heart disease (CHD) among women in the United States. However, the magnitude of differences in metabolomic profiles between Black and White women in the United States has not been well-described. It also remains unknown whether such differences are related to differences in CHD risk. METHODS: Plasma metabolomic profiles were analyzed using liquid chromatography-tandem mass spectrometry in the WHI-OS (Women's Health Initiative-Observational Study; 138 Black and 696 White women), WHI-HT trials (WHI-Hormone Therapy; 156 Black and 1138 White women), MESA (Multi-Ethnic Study of Atherosclerosis; 114 Black and 219 White women), JHS (Jackson Heart Study; 1465 Black women with 107 incident CHD cases), and NHS (Nurses' Health Study; 2506 White women with 136 incident CHD cases). First, linear regression models were used to estimate associations between self-reported race and 472 metabolites in WHI-OS (discovery); findings were replicated in WHI-HT and validated in MESA. Second, we used elastic net regression to construct a racial difference metabolomic pattern (RDMP) representing differences in the metabolomic patterns between Black and White women in the WHI-OS; the RDMP was validated in the WHI-HT and MESA. Third, using conditional logistic regressions in the WHI (717 CHD cases and 719 matched controls), we examined associations of metabolites with large differences in levels by race and the RDMP with risk of CHD, and the results were replicated in Black women from the JHS and White women from the NHS. RESULTS: Of the 472 tested metabolites, levels of 259 (54.9%) metabolites, mostly lipid metabolites and amino acids, significantly differed between Black and White women in both WHI-OS and WHI-HT after adjusting for baseline characteristics, socioeconomic status, lifestyle factors, baseline health conditions, and medication use (false discovery rate <0.05); similar trends were observed in MESA. The RDMP, composed of 152 metabolites, was identified in the WHI-OS and showed significantly different distributions between Black and White women in the WHI-HT and MESA. Higher RDMP quartiles were associated with an increased risk of incident CHD (odds ratio=1.51 [0.97-2.37] for the highest quartile comparing to the lowest; Ptrend=0.02), independent of self-reported race and known CHD risk factors. In race-stratified analyses, the RDMP-CHD associations were more pronounced in White women. Similar patterns were observed in Black women from the JHS and White women from the NHS. CONCLUSIONS: Metabolomic profiles significantly and substantially differ between Black and White women and may be associated with CHD risk and racial disparities in US women.

Determinants of Incident Atherosclerotic Cardiovascular Disease Events Among Those With Absent Coronary Artery Calcium: Multi-Ethnic Study of Atherosclerosis.

BACKGROUND: The 2018 American Heart Association/American College of Cardiology/Multisociety cholesterol guideline states that statin therapy may be withheld or delayed among intermediate-risk individuals in the absence of coronary artery calcium (CAC=0). We evaluated whether traditional cardiovascular risk factors are associated with incident atherosclerotic cardiovascular disease (ASCVD) events among individuals with CAC=0 over long-term follow-up. METHODS: We included participants with CAC=0 at baseline from the MESA (Multi-Ethnic Study of Atherosclerosis), a prospective cohort study of individuals free of clinical ASCVD at baseline. We used multivariable-adjusted Cox proportional hazards models to study the association between cardiovascular risk factors (cigarette smoking, diabetes, hypertension, preventive medication use [aspirin and statin], family history of premature ASCVD, chronic kidney disease, waist circumference, lipid and inflammatory markers) and adjudicated incident ASCVD outcomes. RESULTS: We studied 3416 individuals (mean [SD] age 58 [9] years; 63% were female, 33% White, 31% Black, 12% Chinese American, and 24% Hispanic). Over a median follow-up of 16 years, there were 189 ASCVD events (composite of coronary heart disease and stroke) of which 91 were coronary heart disease, 88 were stroke, and 10 were both coronary heart disease and stroke events. The unadjusted event rates of ASCVD were ≤5 per 1000 person-years among individuals with CAC=0 for most risk factors with the exception of current cigarette smoking (7.3), diabetes (8.9), hypertension (5.4), and chronic kidney disease (6.8). After multivariable adjustment, risk factors that were significantly associated with ASCVD included current cigarette smoking: hazard ratio, 2.12 (95% CI, 1.32-3.42); diabetes: hazard ratio, 1.68 (95% CI, 1.01-2.80); and hypertension: hazard ratio, 1.57 (95% CI, 1.06-2.33). CONCLUSIONS: Current cigarette smoking, diabetes, and hypertension are independently associated with incident ASCVD over a 16-year follow-up among those with CAC=0.

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.

Genome-wide association study of body fat distribution traits in Hispanics/Latinos from the HCHS/SOL.

Central obesity is a leading health concern with a great burden carried by ethnic minority populations, especially Hispanics/Latinos. Genetic factors contribute to the obesity burden overall and to inter-population differences. We aimed to identify the loci associated with central adiposity measured as waist-to-hip ratio (WHR), waist circumference (WC) and hip circumference (HIP) adjusted for body mass index (adjBMI) by using the Hispanic Community Health Study/Study of Latinos (HCHS/SOL); determine if differences in associations differ by background group within HCHS/SOL and determine whether previously reported associations generalize to HCHS/SOL. Our analyses included 7472 women and 5200 men of mainland (Mexican, Central and South American) and Caribbean (Puerto Rican, Cuban and Dominican) background residing in the USA. We performed genome-wide association analyses stratified and combined across sexes using linear mixed-model regression. We identified 16 variants for waist-to-hip ratio adjusted for body mass index (WHRadjBMI), 22 for waist circumference adjusted for body mass index (WCadjBMI) and 28 for hip circumference adjusted for body mass index (HIPadjBMI), which reached suggestive significance (P < 1 × 10-6). Many loci exhibited differences in strength of associations by ethnic background and sex. We brought a total of 66 variants forward for validation in cohorts (N = 34 161) with participants of Hispanic/Latino, African and European descent. We confirmed four novel loci (P < 0.05 and consistent direction of effect, and P < 5 × 10-8 after meta-analysis), including two for WHRadjBMI (rs13301996, rs79478137); one for WCadjBMI (rs3168072) and one for HIPadjBMI (rs28692724). Also, we generalized previously reported associations to HCHS/SOL, (8 for WHRadjBMI, 10 for WCadjBMI and 12 for HIPadjBMI). Our study highlights the importance of large-scale genomic studies in ancestrally diverse Hispanic/Latino populations for identifying and characterizing central obesity susceptibility that may be ancestry-specific.

Epigenome-wide association study of mitochondrial genome copy number.

We conducted cohort- and race-specific epigenome-wide association analyses of mitochondrial deoxyribonucleic acid (mtDNA) copy number (mtDNA CN) measured in whole blood from participants of African and European origins in five cohorts (n = 6182, mean age = 57-67 years, 65% women). In the meta-analysis of all the participants, we discovered 21 mtDNA CN-associated DNA methylation sites (CpG) (P < 1 × 10-7), with a 0.7-3.0 standard deviation increase (3 CpGs) or decrease (18 CpGs) in mtDNA CN corresponding to a 1% increase in DNA methylation. Several significant CpGs have been reported to be associated with at least two risk factors (e.g. chronological age or smoking) for cardiovascular disease (CVD). Five genes [PR/SET domain 16, nuclear receptor subfamily 1 group H member 3 (NR1H3), DNA repair protein, DNA polymerase kappa and decaprenyl-diphosphate synthase subunit 2], which harbor nine significant CpGs, are known to be involved in mitochondrial biosynthesis and functions. For example, NR1H3 encodes a transcription factor that is differentially expressed during an adipose tissue transition. The methylation level of cg09548275 in NR1H3 was negatively associated with mtDNA CN (effect size = -1.71, P = 4 × 10-8) and was positively associated with the NR1H3 expression level (effect size = 0.43, P = 0.0003), which indicates that the methylation level in NR1H3 may underlie the relationship between mtDNA CN, the NR1H3 transcription factor and energy expenditure. In summary, the study results suggest that mtDNA CN variation in whole blood is associated with DNA methylation levels in genes that are involved in a wide range of mitochondrial activities. These findings will help reveal molecular mechanisms between mtDNA CN and CVD.

Allele-specific variation at APOE increases nonalcoholic fatty liver disease and obesity but decreases risk of Alzheimer's disease and myocardial infarction.

Nonalcoholic fatty liver disease (NAFLD) is a leading cause of chronic liver disease and is highly correlated with metabolic disease. NAFLD results from environmental exposures acting on a susceptible polygenic background. This study performed the largest multiethnic investigation of exonic variation associated with NAFLD and correlated metabolic traits and diseases. An exome array meta-analysis was carried out among eight multiethnic population-based cohorts (n = 16 492) with computed tomography (CT) measured hepatic steatosis. A fixed effects meta-analysis identified five exome-wide significant loci (P < 5.30 × 10-7); including a novel signal near TOMM40/APOE. Joint analysis of TOMM40/APOE variants revealed the TOMM40 signal was attributed to APOE rs429358-T; APOE rs7412 was not associated with liver attenuation. Moreover, rs429358-T was associated with higher serum alanine aminotransferase, liver steatosis, cirrhosis, triglycerides and obesity; as well as, lower cholesterol and decreased risk of myocardial infarction and Alzheimer's disease (AD) in phenome-wide association analyses in the Michigan Genomics Initiative, United Kingdom Biobank and/or public datasets. These results implicate APOE in imaging-based identification of NAFLD. This association may or may not translate to nonalcoholic steatohepatitis; however, these results indicate a significant association with advanced liver disease and hepatic cirrhosis. These findings highlight allelic heterogeneity at the APOE locus and demonstrate an inverse link between NAFLD and AD at the exome level in the largest analysis to date.

Circulating metabolites may illustrate relationship of alcohol consumption with cardiovascular disease.

BACKGROUND: Metabolite signatures of long-term alcohol consumption are lacking. To better understand the molecular basis linking alcohol drinking and cardiovascular disease (CVD), we investigated circulating metabolites associated with long-term alcohol consumption and examined whether these metabolites were associated with incident CVD. METHODS: Cumulative average alcohol consumption (g/day) was derived from the total consumption of beer, wine, and liquor on average of 19 years in 2428 Framingham Heart Study Offspring participants (mean age 56 years, 52% women). We used linear mixed models to investigate the associations of alcohol consumption with 211 log-transformed plasma metabolites, adjusting for age, sex, batch, smoking, diet, physical activity, BMI, and familial relationship. Cox models were used to test the association of alcohol-related metabolite scores with fatal and nonfatal incident CVD (myocardial infarction, coronary heart disease, stroke, and heart failure). RESULTS: We identified 60 metabolites associated with cumulative average alcohol consumption (p < 0.05/211 ≈ 0.00024). For example, 1 g/day increase of alcohol consumption was associated with higher levels of cholesteryl esters (e.g., CE 16:1, beta = 0.023 ± 0.002, p = 6.3e - 45) and phosphatidylcholine (e.g., PC 32:1, beta = 0.021 ± 0.002, p = 3.1e - 38). Survival analysis identified that 10 alcohol-associated metabolites were also associated with a differential CVD risk after adjusting for age, sex, and batch. Further, we built two alcohol consumption weighted metabolite scores using these 10 metabolites and showed that, with adjustment age, sex, batch, and common CVD risk factors, the two scores had comparable but opposite associations with incident CVD, hazard ratio 1.11 (95% CI = [1.02, 1.21], p = 0.02) vs 0.88 (95% CI = [0.78, 0.98], p = 0.02). CONCLUSIONS: We identified 60 long-term alcohol consumption-associated metabolites. The association analysis with incident CVD suggests a complex metabolic basis between alcohol consumption and CVD.