Proteomic prediction of incident heart failure and its main subtypes. - Emilsson, Valur; Jonsson, Brynjolfur G; Austin, Thomas R; Gudmundsdottir, Valborg; Axelsson, Gisli T; Frick, Elisabet A; Jonmundsson, Thorarinn; Steindorsdottir, Anna E; Loureiro, Joseph; Brody, Jennifer A; Aspelund, Thor; Launer, Lenore J; et al. (Eur J Heart Fail, 2023)
AIM: To examine the ability of serum proteins in predicting future heart failure (HF) events, including HF with reduced or preserved ejection fraction (HFrEF or HFpEF), in relation to event time, and with or without considering established HF-associated clinical variables. METHODS AND RESULTS: In the prospective population-based Age, Gene/Environment Susceptibility Reykjavik Study (AGES-RS), 440 individuals developed HF after their first visit with a median follow-up of 5.45 years. Among them, 167 were diagnosed with HFrEF and 188 with HFpEF. A least absolute shrinkage and selection operator regression model with nonparametric bootstrap were used to select predictors from an analysis of 4782 serum proteins, and several pre-established clinical parameters linked to HF. A subset of 8-10 distinct or overlapping serum proteins predicted different future HF outcomes, and C-statistics were used to assess discrimination, revealing proteins combined with a C-index of 0.80 for all incident HF, 0.78 and 0.80 for incident HFpEF or HFrEF, respectively. In the AGES-RS, protein panels alone encompassed the risk contained in the clinical information and improved the performance characteristics of prediction models based on NT-proBNP and clinical risk factors. Finally, the protein predictors performed particularly well close to the time of an HF event, an outcome that was replicated in the Cardiovascular Health Study (CHS). CONCLUSION: A small number of circulating proteins accurately predicted future HF in the AGES-RS cohort of older adults, and they alone encompass the risk information found in a collection of clinical data. Incident HF events were predicted up to eight years, with predictor performance significantly improving for events occurring less than one year ahead, a finding replicated in an external cohort study. This article is protected by copyright. All rights reserved.

Rare variants in long non-coding RNAs are associated with blood lipid levels in the TOPMed Whole Genome Sequencing Study. - Wang, Yuxuan; Selvaraj, Margaret Sunitha; Li, Xihao; Li, Zilin; Holdcraft, Jacob A; Arnett, Donna K; Bis, Joshua C; Blangero, John; Boerwinkle, Eric; Bowden, Donald W; Cade, Brian E; Carlson, Jenna C; et al. (medRxiv, 2023)
Long non-coding RNAs (lncRNAs) are known to perform important regulatory functions. Large-scale whole genome sequencing (WGS) studies and new statistical methods for variant set tests now provide an opportunity to assess the 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 blood lipid levels (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 variations and rare protein coding variations 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 lncRNA, implicating new therapeutic opportunities.

Sleep Architecture, Obstructive Sleep Apnea, and Cognitive Function in Adults. - Pase, Matthew P; Harrison, Stephanie; Misialek, Jeffrey R; Kline, Christopher E; Cavuoto, Marina; Baril, Andree-Ann; Yiallourou, Stephanie; Bisson, Alycia; Himali, Dibya; Leng, Yue; Yang, Qiong; Seshadri, Sudha; et al. (JAMA Netw Open, 2023)
IMPORTANCE: Good sleep is essential for health, yet associations between sleep and dementia risk remain incompletely understood. The Sleep and Dementia Consortium was established to study associations between polysomnography (PSG)-derived sleep and the risk of dementia and related cognitive and brain magnetic resonance imaging endophenotypes. OBJECTIVE: To investigate association of sleep architecture and obstructive sleep apnea (OSA) with cognitive function in the Sleep and Dementia Consortium. DESIGN, SETTING, AND PARTICIPANTS: The Sleep and Dementia Consortium curated data from 5 population-based cohorts across the US with methodologically consistent, overnight, home-based type II PSG and neuropsychological assessments over 5 years of follow-up: the Atherosclerosis Risk in Communities study, Cardiovascular Health Study, Framingham Heart Study (FHS), Osteoporotic Fractures in Men Study, and Study of Osteoporotic Fractures. Sleep metrics were harmonized centrally and then distributed to participating cohorts for cohort-specific analysis using linear regression; study-level estimates were pooled in random effects meta-analyses. Results were adjusted for demographic variables, the time between PSG and neuropsychological assessment (0-5 years), body mass index, antidepressant use, and sedative use. There were 5946 participants included in the pooled analyses without stroke or dementia. Data were analyzed from March 2020 to June 2023. EXPOSURES: Measures of sleep architecture and OSA derived from in-home PSG. MAIN OUTCOMES AND MEASURES: The main outcomes were global cognitive composite z scores derived from principal component analysis, with cognitive domains investigated as secondary outcomes. Higher scores indicated better performance. RESULTS: Across cohorts, 5946 adults (1875 females [31.5%]; mean age range, 58-89 years) were included. The median (IQR) wake after sleep onset time ranged from 44 (27-73) to 101 (66-147) minutes, and the prevalence of moderate to severe OSA ranged from 16.9% to 28.9%. Across cohorts, higher sleep maintenance efficiency (pooled β per 1% increase, 0.08; 95% CI, 0.03 to 0.14; P < .01) and lower wake after sleep onset (pooled β per 1-min increase, -0.07; 95% CI, -0.13 to -0.01 per 1-min increase; P = .02) were associated with better global cognition. Mild to severe OSA (apnea-hypopnea index [AHI] ≥5) was associated with poorer global cognition (pooled β, -0.06; 95% CI, -0.11 to -0.01; P = .01) vs AHI less than 5; comparable results were found for moderate to severe OSA (pooled β, -0.06; 95% CI, -0.11 to -0.01; P = .02) vs AHI less than 5. Differences in sleep stages were not associated with cognition. CONCLUSIONS AND RELEVANCE: This study found that better sleep consolidation and the absence of OSA were associated with better global cognition over 5 years of follow-up. These findings suggest that the role of interventions to improve sleep for maintaining cognitive function requires investigation.

A statistical framework for powerful multi-trait rare variant analysis in large-scale whole-genome sequencing studies. - Li, Xihao; Chen, Han; Selvaraj, Margaret Sunitha; Van Buren, Eric; Zhou, Hufeng; Wang, Yuxuan; Sun, Ryan; McCaw, Zachary R; Yu, Zhi; Arnett, Donna K; Bis, Joshua C; Blangero, John; et al. (bioRxiv, 2023)
Large-scale whole-genome sequencing (WGS) studies have improved our understanding of the contributions of coding and noncoding rare variants to complex human traits. Leveraging association effect sizes across multiple traits in WGS rare variant association analysis can improve statistical power over single-trait analysis, and also detect pleiotropic genes and regions. Existing multi-trait methods have limited ability to perform rare variant analysis of large-scale WGS data. We propose MultiSTAAR, a statistical framework and computationally-scalable analytical pipeline for functionally-informed multi-trait rare variant analysis in large-scale WGS studies. MultiSTAAR accounts for relatedness, population structure and correlation among phenotypes by jointly analyzing multiple traits, and further empowers rare variant association analysis by incorporating multiple functional annotations. We applied MultiSTAAR to jointly analyze three lipid traits (low-density lipoprotein cholesterol, high-density lipoprotein cholesterol and triglycerides) in 61,861 multi-ethnic samples from the Trans-Omics for Precision Medicine (TOPMed) Program. We discovered new associations with lipid traits missed by single-trait analysis, including rare variants within an enhancer of and an intergenic region on chromosome 1.

Time-to-Event Genome-Wide Association Study for Incident Cardiovascular Disease in People with Type 2 Diabetes Mellitus. - Kwak, Soo Heon; Hernandez-Cancela, Ryan B; DiCorpo, Daniel A; Condon, David E; Merino, Jordi; Wu, Peitao; Brody, Jennifer A; Yao, Jie; Guo, Xiuqing; Ahmadizar, Fariba; Meyer, Mariah; Sincan, Murat; et al. (medRxiv, 2023)
BACKGROUND: Type 2 diabetes mellitus (T2D) confers a two- to three-fold increased risk of cardiovascular disease (CVD). However, the mechanisms underlying increased CVD risk among people with T2D are only partially understood. We hypothesized that a genetic association study among people with T2D at risk for developing incident cardiovascular complications could provide insights into molecular genetic aspects underlying CVD. METHODS: From 16 studies of the Cohorts for Heart & Aging Research in Genomic Epidemiology (CHARGE) Consortium, we conducted a multi-ancestry time-to-event genome-wide association study (GWAS) for incident CVD among people with T2D using Cox proportional hazards models. Incident CVD was defined based on a composite of coronary artery disease (CAD), stroke, and cardiovascular death that occurred at least one year after the diagnosis of T2D. Cohort-level estimated effect sizes were combined using inverse variance weighted fixed effects meta-analysis. We also tested 204 known CAD variants for association with incident CVD among patients with T2D. RESULTS: A total of 49,230 participants with T2D were included in the analyses (31,118 European ancestries and 18,112 non-European ancestries) which consisted of 8,956 incident CVD cases over a range of mean follow-up duration between 3.2 and 33.7 years (event rate 18.2%). We identified three novel, distinct genetic loci for incident CVD among individuals with T2D that reached the threshold for genome-wide significance ( <5.0×10 ): rs147138607 (intergenic variant between and ) with a hazard ratio (HR) 1.23, 95% confidence interval (CI) 1.15 - 1.32, =3.6×10 , rs11444867 (intergenic variant near ) with HR 1.89, 95% CI 1.52 - 2.35, =9.9×10 , and rs335407 (intergenic variant between and ) HR 1.25, 95% CI 1.16 - 1.35, =1.5×10 . Among 204 known CAD loci, 32 were associated with incident CVD in people with T2D with <0.05, and 5 were significant after Bonferroni correction ( <0.00024, 0.05/204). A polygenic score of these 204 variants was significantly associated with incident CVD with HR 1.14 (95% CI 1.12 - 1.16) per 1 standard deviation increase ( =1.0×10 ). CONCLUSIONS: The data point to novel and known genomic regions associated with incident CVD among individuals with T2D. CLINICAL PERSPECTIVE: We conducted a large-scale multi-ancestry time-to-event GWAS to identify genetic variants associated with CVD among people with T2D. Three variants were significantly associated with incident CVD in people with T2D: rs147138607 (intergenic variant between and ), rs11444867 (intergenic variant near ), and rs335407 (intergenic variant between and ). A polygenic score composed of known CAD variants identified in the general population was significantly associated with the risk of CVD in people with T2D. There are genetic risk factors specific to T2D that could at least partially explain the excess risk of CVD in people with T2D.In addition, we show that people with T2D have enrichment of known CAD association signals which could also explain the excess risk of CVD.

Traditional and novel risk factors for incident aortic stenosis in community-dwelling older adults. - Massera, Daniele; Bartz, Traci M; Biggs, Mary L; Sotoodehnia, Nona; Reiner, Alexander P; Semba, Richard D; Gottdiener, John S; Psaty, Bruce M; Owens, David S; Kizer, Jorge R (Heart, 2023)
OBJECTIVES: Calcific aortic stenosis (AS) is the most common valvular disease in older adults, yet its risk factors remain insufficiently studied in this population. Such studies are necessary to enhance understanding of mechanisms, disease management and therapeutics. METHODS: The Cardiovascular Health Study is a population-based investigation of older adults that completed adjudication of incident AS over long-term follow-up. We evaluated traditional cardiovascular risk factors or disease, as well as novel risk factors from lipid, inflammatory and mineral metabolism pathways, in relation to incident moderate or severe AS (including AS procedures) and clinically significant AS (severe AS, including procedures). RESULTS: Of 5390 participants (age 72.9±5.6 years, 57.6% female, 12.5% black), 287 developed moderate or severe AS, and 175 clinically significant AS, during median follow-up of 13.1 years. After full adjustment, age (HR=1.66 per SD (95% CI=1.45, 1.91)), male sex (HR=1.41 (1.06, 1.87)), diabetes (HR=1.53 (1.10, 2.13)), coronary heart disease (CHD, HR=1.36 (1.01, 1.84)), lipoprotein-associated phospholipase-A (LpPLA) activity (HR=1.21 per SD (1.07, 1.37)) and sCD14 (HR=1.16 per SD (1.01, 1.34)) were associated with incident moderate/severe AS, while black race demonstrated an inverse association (HR=0.40 (0.24, 0.65)), and creatinine-based estimated glomerular filtration rate (eGFR) showed a U-shaped relationship. Findings were similar for clinically significant AS, although CHD and sCD14 fell short of significance, but interleukin-(IL) 6 showed a positive association. CONCLUSION: This comprehensive evaluation of risk factors for long-term incidence of AS identified associations for diabetes and prevalent CHD, LpPLA activity, sCD14 and IL-6, and eGFR. These factors may hold clues to biology, preventive efforts and potential therapeutics for those at highest risk.

A Type 1 Diabetes Polygenic Score Is Not Associated With Prevalent Type 2 Diabetes in Large Population Studies. - Srinivasan, Shylaja; Wu, Peitao; Mercader, Josep M; Udler, Miriam S; Porneala, Bianca C; Bartz, Traci M; Floyd, James S; Sitlani, Colleen; Guo, Xiquing; Haessler, Jeffrey; Kooperberg, Charles; Liu, Jun; et al. (J Endocr Soc, 2023)
CONTEXT: Both type 1 diabetes (T1D) and type 2 diabetes (T2D) have significant genetic contributions to risk and understanding their overlap can offer clinical insight. OBJECTIVE: We examined whether a T1D polygenic score (PS) was associated with a diagnosis of T2D in the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) consortium. METHODS: We constructed a T1D PS using 79 known single nucleotide polymorphisms associated with T1D risk. We analyzed 13 792 T2D cases and 14 169 controls from CHARGE cohorts to determine the association between the T1D PS and T2D prevalence. We validated findings in an independent sample of 2256 T2D cases and 27 052 controls from the Mass General Brigham Biobank (MGB Biobank). As secondary analyses in 5228 T2D cases from CHARGE, we used multivariable regression models to assess the association of the T1D PS with clinical outcomes associated with T1D. RESULTS: The T1D PS was not associated with T2D both in CHARGE ( = .15) and in the MGB Biobank ( = .87). The partitioned human leukocyte antigens only PS was associated with T2D in CHARGE (OR 1.02 per 1 SD increase in PS, 95% CI 1.01-1.03, = .006) but not in the MGB Biobank. The T1D PS was weakly associated with insulin use (OR 1.007, 95% CI 1.001-1.012, = .03) in CHARGE T2D cases but not with other outcomes. CONCLUSION: In large biobank samples, a common variant PS for T1D was not consistently associated with prevalent T2D. However, possible heterogeneity in T2D cannot be ruled out and future studies are needed do subphenotyping.

Type 2 Diabetes Modifies the Association of CAD Genomic Risk Variants With Subclinical Atherosclerosis. - Hasbani, Natalie R; Westerman, Kenneth E; Heon Kwak, Soo; Chen, Han; Li, Xihao; DiCorpo, Daniel; Wessel, Jennifer; Bis, Joshua C; Sarnowski, Chloe; Wu, Peitao; Bielak, Lawrence F; Guo, Xiuqing; et al. (Circ Genom Precis Med, 2023)
BACKGROUND: Individuals with type 2 diabetes (T2D) have an increased risk of coronary artery disease (CAD), but questions remain about the underlying pathology. Identifying which CAD loci are modified by T2D in the development of subclinical atherosclerosis (coronary artery calcification [CAC], carotid intima-media thickness, or carotid plaque) may improve our understanding of the mechanisms leading to the increased CAD in T2D. METHODS: We compared the common and rare variant associations of known CAD loci from the literature on CAC, carotid intima-media thickness, and carotid plaque in up to 29 670 participants, including up to 24 157 normoglycemic controls and 5513 T2D cases leveraging whole-genome sequencing data from the Trans-Omics for Precision Medicine program. We included first-order T2D interaction terms in each model to determine whether CAD loci were modified by T2D. The genetic main and interaction effects were assessed using a joint test to determine whether a CAD variant, or gene-based rare variant set, was associated with the respective subclinical atherosclerosis measures and then further determined whether these loci had a significant interaction test. RESULTS: Using a Bonferroni-corrected significance threshold of <1.6×10, we identified 3 genes (, , and ) associated with CAC and 2 genes ( and ) associated with carotid intima-media thickness and carotid plaque, respectively, through gene-based rare variant set analysis. Both and also had significantly different associations for CAC in T2D cases versus controls. No significant interaction tests were identified through the candidate single-variant analysis. CONCLUSIONS: These results highlight T2D as an important modifier of rare variant associations in CAD loci with CAC.

Whole genome analysis of plasma fibrinogen reveals population-differentiated genetic regulators with putative liver roles. - Huffman, Jennifer E; Nicolas, Jayna; Hahn, Julie; Heath, Adam S; Raffield, Laura M; Yanek, Lisa R; Brody, Jennifer A; Thibord, Florian; Almasy, Laura; Bartz, Traci M; Bielak, Lawrence F; Bowler, Russell P; et al. (medRxiv, 2023)
UNLABELLED: Genetic studies have identified numerous regions associated with plasma fibrinogen levels in Europeans, yet missing heritability and limited inclusion of non-Europeans necessitates further studies with improved power and sensitivity. Compared with array-based genotyping, whole genome sequencing (WGS) data provides better coverage of the genome and better representation of non-European variants. To better understand the genetic landscape regulating plasma fibrinogen levels, we meta-analyzed WGS data from the NHLBI's Trans-Omics for Precision Medicine (TOPMed) program (n=32,572), with array-based genotype data from the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) Consortium (n=131,340) imputed to the TOPMed or Haplotype Reference Consortium panel. We identified 18 loci that have not been identified in prior genetic studies of fibrinogen. Of these, four are driven by common variants of small effect with reported MAF at least 10% higher in African populations. Three ( , and signals contain predicted deleterious missense variants. Two loci, and , each harbor two conditionally distinct, non-coding variants. The gene region encoding the protein chain subunits ( ), contains 7 distinct signals, including one novel signal driven by rs28577061, a variant common (MAF=0.180) in African reference panels but extremely rare (MAF=0.008) in Europeans. Through phenome-wide association studies in the VA Million Veteran Program, we found associations between fibrinogen polygenic risk scores and thrombotic and inflammatory disease phenotypes, including an association with gout. Our findings demonstrate the utility of WGS to augment genetic discovery in diverse populations and offer new insights for putative mechanisms of fibrinogen regulation. KEY POINTS: Largest and most diverse genetic study of plasma fibrinogen identifies 54 regions (18 novel), housing 69 conditionally distinct variants (20 novel).Sufficient power achieved to identify signal driven by African population variant.Links to (1) liver enzyme, blood cell and lipid genetic signals, (2) liver regulatory elements, and (3) thrombotic and inflammatory disease.

Whole Genome Analysis of Venous Thromboembolism: the Trans-Omics for Precision Medicine Program. - Seyerle, Amanda A; Laurie, Cecelia A; Coombes, Brandon J; Jain, Deepti; Conomos, Matthew P; Brody, Jennifer; Chen, Ming-Huei; Gogarten, Stephanie M; Beutel, Kathleen M; Gupta, Namrata; Heckbert, Susan R; Jackson, Rebecca D; et al. (Circ Genom Precis Med, 2023)
Background Risk for venous thromboembolism has a strong genetic component. Whole genome sequencingfrom the Trans-Omics for Precision Medicine program allowed us to look for new associations, particularly rare variants missed by standard genome-wide association studies. Methods The 3793 cases and 7834 controls (11.6% of cases were Black, Hispanic/Latino, or Asian American) were analyzed using a single variant approach and an aggregate gene-based approach using our primary filter (included only loss-of-function and missense variants predicted to be deleterious) and our secondary filter (included all missense variants). Results Single variant analyses identified associations at 5 known loci. Aggregate gene-based analyses identified only (odds ratio, 6.2 for carriers of rare variants; =7.4×10) when using our primary filter. Employing our secondary variant filter led to a smaller effect size at (odds ratio, 3.8; =1.6×10), while excluding variants found only in rare isoforms led to a larger one (odds ratio, 7.5). Different filtering strategies improved the signal for 2 other known genes: became significant (minimum =1.8×10 with the secondary filter), while did not (minimum =4.4×10 with minor allele frequency <0.0005). Results were largely the same when restricting the analyses to include only unprovoked cases; however, one novel gene, , became significant (=4.4×10 using all missense variants with minor allele frequency <0.0005). Conclusions Here, we have demonstrated the importance of using multiple variant filtering strategies, as we detected additional genes when filtering variants based on their predicted deleteriousness, frequency, and presence on the most expressed isoforms. Our primary analyses did not identify new candidate loci; thus larger follow-up studies are needed to replicate the novel locus and to identify additional rare variation associated with venous thromboembolism.

Whole genome sequence analysis of apparent treatment resistant hypertension status in participants from the Trans-Omics for Precision Medicine program. - Armstrong, Nicole D; Srinivasasainagendra, Vinodh; Ammous, Farah; Assimes, Themistocles L; Beitelshees, Amber L; Brody, Jennifer; Cade, Brian E; Ida Chen, Yii-Der; Chen, Han; de Vries, Paul S; Floyd, James S; Franceschini, Nora; et al. (Front Genet, 2023)
Apparent treatment-resistant hypertension (aTRH) is characterized by the use of four or more antihypertensive (AHT) classes to achieve blood pressure (BP) control. In the current study, we conducted single-variant and gene-based analyses of aTRH among individuals from 12 Trans-Omics for Precision Medicine cohorts with whole-genome sequencing data. Cases were defined as individuals treated for hypertension (HTN) taking three different AHT classes, with average systolic BP ≥ 140 or diastolic BP ≥ 90 mmHg, or four or more medications regardless of BP ( = 1,705). A normotensive control group was defined as individuals with BP < 140/90 mmHg ( = 22,079), not on AHT medication. A second control group comprised individuals who were treatment responsive on one AHT medication with BP < 140/ 90 mmHg ( = 5,424). Logistic regression with kinship adjustment using the Scalable and Accurate Implementation of Generalized mixed models (SAIGE) was performed, adjusting for age, sex, and genetic ancestry. We assessed variants using SKAT-O in rare-variant analyses. Single-variant and gene-based tests were conducted in a pooled multi-ethnicity stratum, as well as self-reported ethnic/racial strata (European and African American). One variant in the known HTN locus, , was a top finding in the multi-ethnic analysis ( = 8.23E-07) for the normotensive control group [rs12476527, odds ratio (95% confidence interval) = 0.80 (0.74-0.88)]. This variant was replicated in the Vanderbilt University Medical Center's DNA repository data. Aggregate gene-based signals included the genes and . Additional work validating these loci in larger, more diverse populations, is warranted to determine whether these regions influence the pathobiology of aTRH.

WHOLE GENOME SEQUENCING ANALYSIS OF BODY MASS INDEX IDENTIFIES NOVEL AFRICAN ANCESTRY-SPECIFIC RISK ALLELE. - Zhang, Xinruo; Brody, Jennifer A; Graff, Mariaelisa; Highland, Heather M; Chami, Nathalie; Xu, Hanfei; Wang, Zhe; Ferrier, Kendra; Chittoor, Geetha; Josyula, Navya S; Li, Xihao; Li, Zilin; et al. (medRxiv, 2023)
Obesity is a major public health crisis associated with high mortality rates. Previous genome-wide association studies (GWAS) investigating body mass index (BMI) have largely relied on imputed data from European individuals. This study leveraged whole-genome sequencing (WGS) data from 88,873 participants from the Trans-Omics for Precision Medicine (TOPMed) Program, of which 51% were of non-European population groups. We discovered 18 BMI-associated signals ( < 5 × 10 ). Notably, we identified and replicated a novel low frequency single nucleotide polymorphism (SNP) in that was common in individuals of African descent. Using a diverse study population, we further identified two novel secondary signals in known BMI loci and pinpointed two likely causal variants in the and loci. Our work demonstrates the benefits of combining WGS and diverse cohorts in expanding current catalog of variants and genes confer risk for obesity, bringing us one step closer to personalized medicine.

Whole Genome Sequencing Based Analysis of Inflammation Biomarkers in the Trans-Omics for Precision Medicine (TOPMed) Consortium. - Jiang, Min-Zhi; Gaynor, Sheila M; Li, Xihao; Van Buren, Eric; Stilp, Adrienne; Buth, Erin; Wang, Fei Fei; Manansala, Regina; Gogarten, Stephanie M; Li, Zilin; Polfus, Linda M; Salimi, Shabnam; et al. (bioRxiv, 2023)
Inflammation biomarkers can provide valuable insight into the role of inflammatory processes in many diseases and conditions. Sequencing based analyses of such biomarkers can also serve as an exemplar of the genetic architecture of quantitative traits. To evaluate the biological insight, which can be provided by a multi-ancestry, whole-genome based association study, we performed a comprehensive analysis of 21 inflammation biomarkers from up to 38,465 individuals with whole-genome sequencing from the Trans-Omics for Precision Medicine (TOPMed) program. We identified 22 distinct single-variant associations across 6 traits - E-selectin, intercellular adhesion molecule 1, interleukin-6, lipoprotein-associated phospholipase A2 activity and mass, and P-selectin - that remained significant after conditioning on previously identified associations for these inflammatory biomarkers. We further expanded upon known biomarker associations by pairing the single-variant analysis with a rare variant set-based analysis that further identified 19 significant rare variant set-based associations with 5 traits. These signals were distinct from both significant single variant association signals within TOPMed and genetic signals observed in prior studies, demonstrating the complementary value of performing both single and rare variant analyses when analyzing quantitative traits. We also confirm several previously reported signals from semi-quantitative proteomics platforms. Many of these signals demonstrate the extensive allelic heterogeneity and ancestry-differentiated variant-trait associations common for inflammation biomarkers, a characteristic we hypothesize will be increasingly observed with well-powered, large-scale analyses of complex traits.

Whole-Genome Sequencing Analysis of Human Metabolome in Multi-Ethnic Populations. - Feofanova, Elena V; Brown, Michael R; Alkis, Taryn; Manuel, Astrid M; Li, Xihao; Tahir, Usman A; Li, Zilin; Mendez, Kevin M; Kelly, Rachel S; Qi, Qibin; Chen, Han; Larson, Martin G; et al. (Nat Commun, 2023)
Circulating metabolite levels may reflect the state of the human organism in health and disease, however, the genetic architecture of metabolites is not fully understood. We have performed a whole-genome sequencing association analysis of both common and rare variants in up to 11,840 multi-ethnic participants from five studies with up to 1666 circulating metabolites. We have discovered 1985 novel variant-metabolite associations, and validated 761 locus-metabolite associations reported previously. Seventy-nine novel variant-metabolite associations have been replicated, including three genetic loci located on the X chromosome that have demonstrated its involvement in metabolic regulation. Gene-based analysis have provided further support for seven metabolite-replicated loci pairs and their biologically plausible genes. Among those novel replicated variant-metabolite pairs, follow-up analyses have revealed that 26 metabolites have colocalized with 21 tissues, seven metabolite-disease outcome associations have been putatively causal, and 7 metabolites might be regulated by plasma protein levels. Our results have depicted the genetic contribution to circulating metabolite levels, providing additional insights into understanding human disease.
Keywords: Ethnicity; Genome-Wide Association Study; Humans; Metabolome; Polymorphism, Single Nucleotide; Quantitative Trait Loci

Association analysis of mitochondrial DNA heteroplasmic variants: methods and application. - Sun, Xianbang; Bulekova, Katia; Yang, Jian; Lai, Meng; Pitsillides, Achilleas N; Liu, Xue; Zhang, Yuankai; Guo, Xiuqing; Yong, Qian; Raffield, Laura M; Rotter, Jerome I; Rich, Stephen S; et al. (medRxiv, 2024)
We rigorously assessed a comprehensive association testing framework for heteroplasmy, employing both simulated and real-world data. This framework employed a variant allele fraction (VAF) threshold and harnessed multiple gene-based tests for robust identification and association testing of heteroplasmy. Our simulation studies demonstrated that gene-based tests maintained an appropriate type I error rate at α=0.001. Notably, when 5% or more heteroplasmic variants within a target region were linked to an outcome, burden-extension tests (including the adaptive burden test, variable threshold burden test, and z-score weighting burden test) outperformed the sequence kernel association test (SKAT) and the original burden test. Applying this framework, we conducted association analyses on whole-blood derived heteroplasmy in 17,507 individuals of African and European ancestries (31% of African Ancestry, mean age of 62, with 58% women) with whole genome sequencing data. We performed both cohort- and ancestry-specific association analyses, followed by meta-analysis on bothpooled samples and within each ancestry group. Our results suggest that mtDNA-Enco ded genes/regions are likely to exhibit varying rates in somatic aging, with the notably strong associations observed between heteroplasmy in the and genes ( <0.001) and advance aging by the Original Burden test. In contrast, SKAT identified significant associations ( <0.001) between diabetes and the aggregated effects of heteroplasmy in several protein-coding genes. Further research is warranted to validate these findings. In summary, our proposed statistical framework represents a valuable tool for facilitating association testing of heteroplasmy with disease traits in large human populations.

Association of thyroid dysfunction in individuals ≥ 65 years of age with subclinical cardiac abnormalities. - Barasch, Eddy; Gottdiener, John; Bůzková, Petra; Cappola, Anne; Shah, Sanjiv; DeFilippi, Christopher; Gardin, Julius; Kizer, Jorge R (J Clin Endocrinol Metab, 2024)
CONTEXT: The relationship between thyroid dysfunction and measures of myocardial disease in older individuals remains to be defined. OBJECTIVE: To evaluate the impact of thyroid dysfunction on structure and function of the left-heart chambers and blood markers of cardiac disease. DESIGN: Cross-sectional analysis. SETTING: The Cardiovascular Health Study, a community-based cohort of older individuals recruited from four urban areas in the United States. PATIENTS: Of 3163 participants studied, 2477 were euthyroid, 465 had subclinical hypothyroidism (SCH), 47 overt hypothyroidism (OH), 45 endogenous (endo) subclinical hyperthyroidism (endo-SCT), and 129 had exogenous (exo) SCT due to thyroid hormone supplementation. INTERVENTIONS: Clinical evaluation, blood sampling and biomarker measurement, 2-dimensional and speckle-tracking echocardiography. MAIN OUTCOME MEASURE(S): Left heart myocardial deformation, circulating biomarkers of diastolic overload (NT-proBNP), fibrosis (sST2, gal-3), and cardiomyocyte injury (hs-cTnT). RESULTS: SCH was associated with higher NT-proBNP (beta = 0.17, p = 0.004), whereas OH was associated with higher hs-cTnT (beta = 0.29, p = 0.005). There were also suggestive associations of SCH with higher sST2, as well as endo-SCT with higher gal-3 and lower (worse) left atrial reservoir strain. Left ventricular longitudinal strain and end-diastolic strain rate did not differ significantly from euthyroid participants in SCH, OH, or exo-SCT. CONCLUSIONS: In this free-living elderly cohort, subclinical and overt hypothyroidism were associated with abnormalities of blood biomarkers consistent with diastolic overload and myocardial necrosis respectively, whereas subclinical hyperthyroidism tended to be associated with myocardial fibrosis and decreased left atrial strain. Our findings could represent stage B heart failure and illuminate distinct aspects of the pathobiology of heart disease related to thyroid gland dysfunction with potential clinical implications.

Development and Validation of the American Heart Association's PREVENT Equations. - Khan, Sadiya S; Matsushita, Kunihiro; Sang, Yingying; Ballew, Shoshana H; Grams, Morgan E; Surapaneni, Aditya; Blaha, Michael J; Carson, April P; Chang, Alexander R; Ciemins, Elizabeth; Go, Alan S; Gutierrez, Orlando M; et al. (Circulation, 2024)
BACKGROUND: Multivariable equations are recommended by primary prevention guidelines to assess absolute risk of cardiovascular disease (CVD). However, current equations have several limitations. Therefore, we developed and validated the American Heart Association Predicting Risk of CVD EVENTs (PREVENT) equations among US adults 30 to 79 years of age without known CVD. METHODS: The derivation sample included individual-level participant data from 25 data sets (N=3 281 919) between 1992 and 2017. The primary outcome was CVD (atherosclerotic CVD and heart failure). Predictors included traditional risk factors (smoking status, systolic blood pressure, cholesterol, antihypertensive or statin use, and diabetes) and estimated glomerular filtration rate. Models were sex-specific, race-free, developed on the age scale, and adjusted for competing risk of non-CVD death. Analyses were conducted in each data set and meta-analyzed. Discrimination was assessed using the Harrell C-statistic. Calibration was calculated as the slope of the observed versus predicted risk by decile. Additional equations to predict each CVD subtype (atherosclerotic CVD and heart failure) and include optional predictors (urine albumin-to-creatinine ratio and hemoglobin A1c), and social deprivation index were also developed. External validation was performed in 3 330 085 participants from 21 additional data sets. RESULTS: Among 6 612 004 adults included, mean±SD age was 53±12 years, and 56% were women. Over a mean±SD follow-up of 4.8±3.1 years, there were 211 515 incident total CVD events. The median C-statistics in external validation for CVD were 0.794 (interquartile interval, 0.763-0.809) in female and 0.757 (0.727-0.778) in male participants. The calibration slopes were 1.03 (interquartile interval, 0.81-1.16) and 0.94 (0.81-1.13) among female and male participants, respectively. Similar estimates for discrimination and calibration were observed for atherosclerotic CVD- and heart failure-specific models. The improvement in discrimination was small but statistically significant when urine albumin-to-creatinine ratio, hemoglobin A1c, and social deprivation index were added together to the base model to total CVD (ΔC-statistic [interquartile interval] 0.004 [0.004-0.005] and 0.005 [0.004-0.007] among female and male participants, respectively). Calibration improved significantly when the urine albumin-to-creatinine ratio was added to the base model among those with marked albuminuria (>300 mg/g; 1.05 [0.84-1.20] versus 1.39 [1.14-1.65]; =0.01). CONCLUSIONS: PREVENT equations accurately and precisely predicted risk for incident CVD and CVD subtypes in a large, diverse, and contemporary sample of US adults by using routinely available clinical variables.
Keywords: Adult; Aged; Albumins; American Heart Association; Atherosclerosis; Cardiovascular Diseases; Creatinine; Female; Glycated Hemoglobin; Heart Failure; Humans; Male; Middle Aged; Risk Assessment; Risk Factors

Familial Hypercholesterolemia Variant and Cardiovascular Risk in Individuals With Elevated Cholesterol. - Zhang, Yiyi; Dron, Jacqueline S; Bellows, Brandon K; Khera, Amit V; Liu, Junxiu; Balte, Pallavi P; Oelsner, Elizabeth C; Amr, Sami Samir; Lebo, Matthew S; Nagy, Anna; Peloso, Gina M; Natarajan, Pradeep; et al. (JAMA Cardiol, 2024)
IMPORTANCE: Familial hypercholesterolemia (FH) is a genetic disorder that often results in severely high low-density lipoprotein cholesterol (LDL-C) and high risk of premature coronary heart disease (CHD). However, the impact of FH variants on CHD risk among individuals with moderately elevated LDL-C is not well quantified. OBJECTIVE: To assess CHD risk associated with FH variants among individuals with moderately (130-189 mg/dL) and severely (≥190 mg/dL) elevated LDL-C and to quantify excess CHD deaths attributable to FH variants in US adults. DESIGN, SETTING, AND PARTICIPANTS: A total of 21 426 individuals without preexisting CHD from 6 US cohort studies (Atherosclerosis Risk in Communities study, Coronary Artery Risk Development in Young Adults study, Cardiovascular Health Study, Framingham Heart Study Offspring cohort, Jackson Heart Study, and Multi-Ethnic Study of Atherosclerosis) were included, 63 of whom had an FH variant. Data were collected from 1971 to 2018, and the median (IQR) follow-up was 18 (13-28) years. Data were analyzed from March to May 2023. EXPOSURES: LDL-C, cumulative past LDL-C, FH variant status. MAIN OUTCOMES AND MEASURES: Cox proportional hazards models estimated associations between FH variants and incident CHD. The Cardiovascular Disease Policy Model projected excess CHD deaths associated with FH variants in US adults. RESULTS: Of the 21 426 individuals without preexisting CHD (mean [SD] age 52.1 [15.5] years; 12 041 [56.2%] female), an FH variant was found in 22 individuals with moderately elevated LDL-C (0.3%) and in 33 individuals with severely elevated LDL-C (2.5%). The adjusted hazard ratios for incident CHD comparing those with and without FH variants were 2.9 (95% CI, 1.4-6.0) and 2.6 (95% CI, 1.4-4.9) among individuals with moderately and severely elevated LDL-C, respectively. The association between FH variants and CHD was slightly attenuated when further adjusting for baseline LDL-C level, whereas the association was no longer statistically significant after adjusting for cumulative past LDL-C exposure. Among US adults 20 years and older with no history of CHD and LDL-C 130 mg/dL or higher, more than 417 000 carry an FH variant and were projected to experience more than 12 000 excess CHD deaths in those with moderately elevated LDL-C and 15 000 in those with severely elevated LDL-C compared with individuals without an FH variant. CONCLUSIONS AND RELEVANCE: In this pooled cohort study, the presence of FH variants was associated with a 2-fold higher CHD risk, even when LDL-C was only moderately elevated. The increased CHD risk appeared to be largely explained by the higher cumulative LDL-C exposure in individuals with an FH variant compared to those without. Further research is needed to assess the value of adding genetic testing to traditional phenotypic FH screening.

{A genetic association study of circulating coagulation Factor VIII and von Willebrand Factor levels - de Vries, P. S.; Reventun, P.; Brown, M. R.; Heath, A. S.; Huffman, J. E.; Le, N. Q.; Bebo, A.; Brody, J. A.; Temprano-Sagrera, G.; Raffield, L. M.; Ozel, A. B.; Thibord, F.; et al. (Blood, 2024)
10-9) at seven new loci for FVIII (ST3GAL4, CLEC4M, B3GNT2, ASGR1, F12, KNG1, and TREM1/NCR2) and one for VWF (B3GNT2). VWF, ABO, and STAB2 were associated with FVIII and VWF in gene-based analyses. Multi-phenotype analysis of FVIII and VWF identified another three new loci, including PDIA3. Silencing of B3GNT2 and the previously reported CD36 gene decreased release of FVIII by HLECs, while 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 one 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.

{Genetic associations with dementia-related proteinopathy: Application of item response theory - Katsumata, Y.; Fardo, D. W.; Shade, L. M. P.; Wu, X.; Karanth, S. D.; Hohman, T. J.; Schneider, J. A.; Bennett, D. A.; Farfel, J. M.; Gauthreaux, K.; Mock, C.; Kukull, W. A.; et al. (Alzheimers Dement, 2024)
Although dementia-related proteinopathy has a strong negative impact on public health, and is highly heritable, understanding of the related genetic architecture is incomplete.\ -synuclein, and TDP-43.\ -synuclein/Lewy bodies. Novel suggestive proteinopathy-linked alleles were also discovered, including several (SDHAF1, TMEM68, and ARHGEF28) with colocalization analyses and/or high degrees of biologic credibility.\ A novel methodology using GPCM enabled insights into gene candidates for driving misfolded proteinopathies.\ Latent factor scores for proteinopathies were estimated using a generalized partial credit model. The three latent continuous scores corresponded well with proteinopathy severity. Novel genes associated with proteinopathies were identified. Several genes had high degrees of biologic credibility for dementia risk factors.

Genetic drivers of heterogeneity in type 2 diabetes pathophysiology. - Suzuki, Ken; Hatzikotoulas, Konstantinos; Southam, Lorraine; Taylor, Henry J; Yin, Xianyong; Lorenz, Kim M; Mandla, Ravi; Huerta-Chagoya, Alicia; Melloni, Giorgio E M; Kanoni, Stavroula; Rayner, Nigel W; Bocher, Ozvan; et al. (Nature, 2024)
Type 2 diabetes (T2D) is a heterogeneous disease that develops through diverse pathophysiological processes and molecular mechanisms that are often specific to cell type. 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) 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 scores 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.

{Genome-wide association study of preserved ratio impaired spirometry (PRISm) - Higbee, D. H.; Lirio, A.; Hamilton, F.; Granell, R.; Wyss, A. B.; London, S. J.; Bartz, T. M.; Gharib, S. A.; Cho, M. H.; Wan, E.; Silverman, E.; Crapo, J. D.; et al. (Eur Respir J, 2024)
{0.70. PRISm is associated with respiratory symptoms and comorbidities. Our objective was to discover novel genetic signals for PRISm and see if they provide insight into the pathogenesis of PRISm and associated comorbidities.\ We undertook a genome-wide association study (GWAS) of PRISm in UK Biobank participants (Stage 1), and selected single nucleotide polymorphisms (SNPs) reaching genome-wide significance for replication in 13 cohorts (Stage 2). A combined meta-analysis of Stage 1 and Stage 2 was done to determine top SNPs. We used cross-trait linkage disequilibrium score regression to estimate genome-wide genetic correlation between PRISm and pulmonary and extrapulmonary traits. Phenome-wide association studies of top SNPs were performed.\ =0.12

Human whole-exome genotype data for Alzheimer's disease. - Leung, Yuk Yee; Naj, Adam C; Chou, Yi-Fan; Valladares, Otto; Schmidt, Michael; Hamilton-Nelson, Kara; Wheeler, Nicholas; Lin, Honghuang; Gangadharan, Prabhakaran; Qu, Liming; Clark, Kaylyn; Kuzma, Amanda B; et al. (Nat Commun, 2024)
The heterogeneity of the whole-exome sequencing (WES) data generation methods present a challenge to a joint analysis. Here we present a bioinformatics strategy for joint-calling 20,504 WES samples collected across nine studies and sequenced using ten capture kits in fourteen sequencing centers in the Alzheimer's Disease Sequencing Project. The joint-genotype called variant-called format (VCF) file contains only positions within the union of capture kits. The VCF was then processed specifically to account for the batch effects arising from the use of different capture kits from different studies. We identified 8.2 million autosomal variants. 96.82% of the variants are high-quality, and are located in 28,579 Ensembl transcripts. 41% of the variants are intronic and 1.8% of the variants are with CADD > 30, indicating they are of high predicted pathogenicity. Here we show our new strategy can generate high-quality data from processing these diversely generated WES samples. The improved ability to combine data sequenced in different batches benefits the whole genomics research community.
Keywords: Alzheimer Disease; Computational Biology; Data Accuracy; Exome; Genotype; Humans

Iron Deficiency and Incident Heart Failure in Older Community-Dwelling Individuals. - Sharma, Shilpa; Katz, Ronit; Chaves, Paulo H M; Hoofnagle, Andrew N; Kizer, Jorge R; Bansal, Nisha; Ganz, Tomas; Ix, Joachim H (ESC Heart Fail, 2024)
AIMS: Among persons with prevalent heart failure (HF), iron deficiency has been linked to HF admissions, and intravenous iron replacement improves HF outcomes. Recent studies in persons with chronic kidney disease (CKD) demonstrate that iron deficiency is associated with incident HF. This study aimed to determine the relationship of iron status with incident HF in community-dwelling older adults irrespective of their kidney function. METHODS: In this case-cohort study, 1,006 Cardiovascular Health Study participants (785 from the random sub-cohort [including 193 HF cases] and 221 additional HF cases [N = 414 total HF cases]) aged ≥ 65 years without HF (41% with CKD), we used weighted Cox models to evaluate associations of iron status with incident HF. Participants were categorized based on quartiles of transferrin saturation and ferritin as "iron replete" (27.3%), "functional iron deficiency" (7.7%), "iron deficiency" (11.8%), "mixed iron deficiency" (5.6%), "high iron" (9.3%) and "non-classified" (38.1%), consistent with prior studies. RESULTS: Compared to older persons who were iron replete, those with iron deficiency were at higher risk of incident HF (HR 1.47; 1.02-2.11) in models adjusting for demographics, HF risk factors, and estimated glomerular filtration rate. Other iron categories did not associate with incident HF. The relationship of iron deficiency with incident HF did not differ by CKD status (interaction P value 0.2). CONCLUSIONS: Among community-dwelling elders, iron deficiency is independently associated with incident HF, an association that was similar irrespective of CKD status. Our findings support conduct of clinical trials of iron replacement for prevention of HF in older adults with iron deficiency.

Large-scale circulating proteome association study (CPAS) meta-analysis identifies circulating proteins and pathways predicting incident hip fractures. - Austin, Thomas R; Fink, Howard A; Jalal, Diana I; Törnqvist, Anna E; Bůzková, Petra; Barzilay, Joshua I; Lu, Tianyuan; Carbone, Laura; Gabrielsen, Maiken E; Grahnemo, Louise; Hveem, Kristian; Jonasson, Christian; et al. (J Bone Miner Res, 2024)
Hip fractures are associated with significant disability, high cost, and mortality. However, the exact biological mechanisms underlying susceptibility to hip fractures remain incompletely understood. In an exploratory search of the underlying biology as reflected through the circulating proteome, we performed a comprehensive Circulating Proteome Association Study (CPAS) meta-analysis for incident hip fractures. Analyses included 6430 subjects from two prospective cohort studies (Cardiovascular Health Study and Trøndelag Health Study) with circulating proteomics data (aptamer-based 5 K SomaScan version 4.0 assay; 4979 aptamers). Associations between circulating protein levels and incident hip fractures were estimated for each cohort using age and sex-adjusted Cox regression models. Participants experienced 643 incident hip fractures. Compared with the individual studies, inverse-variance weighted meta-analyses yielded more statistically significant associations, identifying 23 aptamers associated with incident hip fractures (conservative Bonferroni correction 0.05/4979, P < 1.0 × 10-5). The aptamers most strongly associated with hip fracture risk corresponded to two proteins of the growth hormone/insulin growth factor system (GHR and IGFBP2), as well as GDF15 and EGFR. High levels of several inflammation-related proteins (CD14, CXCL12, MMP12, ITIH3) were also associated with increased hip fracture risk. Ingenuity pathway analysis identified reduced LXR/RXR activation and increased acute phase response signaling to be overrepresented among those proteins associated with increased hip fracture risk. These analyses identified several circulating proteins and pathways consistently associated with incident hip fractures. These findings underscore the usefulness of the meta-analytic approach for comprehensive CPAS in a similar manner as has previously been observed for large-scale human genetic studies. Future studies should investigate the underlying biology of these potential novel drug targets.

A Large-Scale Genome-Wide Study of Gene-Sleep Duration Interactions for Blood Pressure in 811,405 Individuals from Diverse Populations. - Nagarajan, Pavithra; Winkler, Thomas W; Bentley, Amy R; Miller, Clint L; Kraja, Aldi T; Schwander, Karen; Lee, Songmi; Wang, Wenyi; Brown, Michael R; Morrison, John L; Giri, Ayush; O'Connell, Jeffrey R; et al. (medRxiv, 2024)
Although both short and long sleep duration are associated with elevated hypertension risk, our understanding of their interplay with biological pathways governing blood pressure remains limited. To address this, we carried out genome-wide cross-population gene-by-short-sleep and long-sleep duration interaction analyses for three blood pressure traits (systolic, diastolic, and pulse pressure) in 811,405 individuals from diverse population groups. We discover 22 novel gene-sleep duration interaction loci for blood pressure, mapped to genes involved in neurological, thyroidal, bone metabolism, and hematopoietic pathways. Non-overlap between short sleep (12) and long sleep (10) interactions underscores the plausibility of distinct influences of both sleep duration extremes in cardiovascular health. With several of our loci reflecting specificity towards population background or sex, our discovery sheds light on the importance of embracing granularity when addressing heterogeneity entangled in gene-environment interactions, and in therapeutic design approaches for blood pressure management.

Methylation patterns associated with C-reactive protein in racially and ethnically diverse populations. - Lundin, Jessica I; Peters, Ulrike; Hu, Yao; Ammous, Farah; Avery, Christy L; Benjamin, Emelia J; Bis, Joshua C; Brody, Jennifer A; Carlson, Chris; Cushman, Mary; Gignoux, Chris; Guo, Xiuqing; et al. (Epigenetics, 2024)
Systemic low-grade inflammation is a feature of chronic disease. C-reactive protein (CRP) is a common biomarker of inflammation and used as an indicator of disease risk; however, the role of inflammation in disease is not completely understood. Methylation is an epigenetic modification in the DNA which plays a pivotal role in gene expression. In this study we evaluated differential DNA methylation patterns associated with blood CRP level to elucidate biological pathways and genetic regulatory mechanisms to improve the understanding of chronic inflammation. The racially and ethnically diverse participants in this study were included as 50% White, 41% Black or African American, 7% Hispanic or Latino/a, and 2% Native Hawaiian, Asian American, American Indian, or Alaska Native (total  = 13,433) individuals. We replicated 113 CpG sites from 87 unique loci, of which five were novel (, , and cg03282312), across a discovery set of 1,150 CpG sites associated with CRP level ( < 1.2E-7). The downstream pathways affected by DNA methylation included the identification of and CpG-gene transcript pairs which contributed to the innate immune response gene enrichment pathway along with , , and . Gene enrichment analysis also identified the nuclear factor-kappaB transcription pathway. Using two-sample Mendelian randomization (MR) we inferred methylation at three CpG sites as causal for CRP levels using both White and Black or African American MR instrument variables. Overall, we identified novel CpG sites and gene transcripts that could be valuable in understanding the specific cellular processes and pathogenic mechanisms involved in inflammation.
Keywords: C-Reactive Protein; CpG Islands; DNA; DNA Methylation; Epigenesis, Genetic; Genome-Wide Association Study; Humans; Inflammation; Intracellular Signaling Peptides and Proteins

Multi-omics and pathway analyses of genome-wide associations implicate regulation and immunity in verbal declarative memory performance. - Mei, Hao; Simino, Jeannette; Li, Lianna; Jiang, Fan; Bis, Joshua C; Davies, Gail; Hill, W David; Xia, Charley; Gudnason, Vilmundur; Yang, Qiong; Lahti, Jari; Smith, Jennifer A; et al. (Alzheimers Res Ther, 2024)
BACKGROUND: Uncovering the functional relevance underlying verbal declarative memory (VDM) genome-wide association study (GWAS) results may facilitate the development of interventions to reduce age-related memory decline and dementia. METHODS: We performed multi-omics and pathway enrichment analyses of paragraph (PAR-dr) and word list (WL-dr) delayed recall GWAS from 29,076 older non-demented individuals of European descent. We assessed the relationship between single-variant associations and expression quantitative trait loci (eQTLs) in 44 tissues and methylation quantitative trait loci (meQTLs) in the hippocampus. We determined the relationship between gene associations and transcript levels in 53 tissues, annotation as immune genes, and regulation by transcription factors (TFs) and microRNAs. To identify significant pathways, gene set enrichment was tested in each cohort and meta-analyzed across cohorts. Analyses of differential expression in brain tissues were conducted for pathway component genes. RESULTS: The single-variant associations of VDM showed significant linkage disequilibrium (LD) with eQTLs across all tissues and meQTLs within the hippocampus. Stronger WL-dr gene associations correlated with reduced expression in four brain tissues, including the hippocampus. More robust PAR-dr and/or WL-dr gene associations were intricately linked with immunity and were influenced by 31 TFs and 2 microRNAs. Six pathways, including type I diabetes, exhibited significant associations with both PAR-dr and WL-dr. These pathways included fifteen MHC genes intricately linked to VDM performance, showing diverse expression patterns based on cognitive status in brain tissues. CONCLUSIONS: VDM genetic associations influence expression regulation via eQTLs and meQTLs. The involvement of TFs, microRNAs, MHC genes, and immune-related pathways contributes to VDM performance in older individuals.
Keywords: Aged; Cognition; Genome-Wide Association Study; Humans; Memory; MicroRNAs; Multiomics; Polymorphism, Single Nucleotide

{Multi-trait analysis characterizes the genetics of thyroid function and identifies causal associations with clinical implications - Sterenborg, R. B. T. M.; Steinbrenner, I.; Li, Y.; Bujnis, M. N.; Naito, T.; Marouli, E.; Galesloot, T. E.; Babajide, O.; Andreasen, L.; Astrup, A.; svold, B. O.; Bandinelli, S.; et al. (Nat Commun, 2024)
T3, whereas the FT4 associations represent the thyroid hormone metabolism. Polygenic risk score and Mendelian randomization analyses showed the effects of genetically determined variation in thyroid function on various clinical outcomes, including cardiovascular risk factors and diseases, autoimmune diseases, and cancer. In conclusion, our results improve the understanding of thyroid hormone physiology and highlight the pleiotropic effects of thyroid function on various diseases.

Plasma proteomic signature of human longevity. - Liu, Xiaojuan; Axelsson, Gisli Thor; Newman, Anne B; Psaty, Bruce M; Boudreau, Robert M; Wu, Chenkai; Arnold, Alice M; Aspelund, Thor; Austin, Thomas R; Gardin, Julius M; Siggeirsdottir, Kristin; Tracy, Russell P; et al. (Aging Cell, 2024)
The identification of protein targets that exhibit anti-aging clinical potential could inform interventions to lengthen the human health span. Most previous proteomics research has been focused on chronological age instead of longevity. We leveraged two large population-based prospective cohorts with long follow-ups to evaluate the proteomic signature of longevity defined by survival to 90 years of age. Plasma proteomics was measured using a SOMAscan assay in 3067 participants from the Cardiovascular Health Study (discovery cohort) and 4690 participants from the Age Gene/Environment Susceptibility-Reykjavik Study (replication cohort). Logistic regression identified 211 significant proteins in the CHS cohort using a Bonferroni-adjusted threshold, of which 168 were available in the replication cohort and 105 were replicated (corrected p value <0.05). The most significant proteins were GDF-15 and N-terminal pro-BNP in both cohorts. A parsimonious protein-based prediction model was built using 33 proteins selected by LASSO with 10-fold cross-validation and validated using 27 available proteins in the validation cohort. This protein model outperformed a basic model using traditional factors (demographics, height, weight, and smoking) by improving the AUC from 0.658 to 0.748 in the discovery cohort and from 0.755 to 0.802 in the validation cohort. We also found that the associations of 169 out of 211 proteins were partially mediated by physical and/or cognitive function. These findings could contribute to the identification of biomarkers and pathways of aging and potential therapeutic targets to delay aging and age-related diseases.

Role of Polyunsaturated Fat in Modifying Cardiovascular Risk Associated With Family History of Cardiovascular Disease: Pooled De Novo Results From 15 Observational Studies. - Laguzzi, Federica; Åkesson, Agneta; Marklund, Matti; Qian, Frank; Gigante, Bruna; Bartz, Traci M; Bassett, Julie K; Birukov, Anna; Campos, Hannia; Hirakawa, Yoichiro; Imamura, Fumiaki; Jäger, Susanne; et al. (Circulation, 2024)
BACKGROUND: It is unknown whether dietary intake of polyunsaturated fatty acids (PUFA) modifies the cardiovascular disease (CVD) risk associated with a family history of CVD. We assessed interactions between biomarkers of low PUFA intake and a family history in relation to long-term CVD risk in a large consortium. METHODS: Blood and tissue PUFA data from 40 885 CVD-free adults were assessed. PUFA levels ≤25th percentile were considered to reflect low intake of linoleic, alpha-linolenic, and eicosapentaenoic/docosahexaenoic acids (EPA/DHA). Family history was defined as having ≥1 first-degree relative who experienced a CVD event. Relative risks with 95% CI of CVD were estimated using Cox regression and meta-analyzed. Interactions were assessed by analyzing product terms and calculating relative excess risk due to interaction. RESULTS: After multivariable adjustments, a significant interaction between low EPA/DHA and family history was observed (product term pooled RR, 1.09 [95% CI, 1.02-1.16]; =0.01). The pooled relative risk of CVD associated with the combined exposure to low EPA/DHA, and family history was 1.41 (95% CI, 1.30-1.54), whereas it was 1.25 (95% CI, 1.16-1.33) for family history alone and 1.06 (95% CI, 0.98-1.14) for EPA/DHA alone, compared with those with neither exposure. The relative excess risk due to interaction results indicated no interactions. CONCLUSIONS: A significant interaction between biomarkers of low EPA/DHA intake, but not the other PUFA, and a family history was observed. This novel finding might suggest a need to emphasize the benefit of consuming oily fish for individuals with a family history of CVD.
Keywords: Animals; Biomarkers; Cardiovascular Diseases; Docosahexaenoic Acids; Fatty Acids, Omega-3; Risk Factors

Trajectory of Cognitive Function After Incident Heart Failure. - Shore, Supriya; Li, Hanyu; Zhang, Min; Whitney, Rachael; Gross, Alden L; Bhatt, Ankeet S; Nallamothu, Brahmajee K; Giordani, Bruno; Briceño, Emily M; Sussman, Jeremy B; Gutierrez, Jose; Yaffe, Kristine; et al. (medRxiv, 2024)
BACKGROUND: The size/magnitude of cognitive changes after incident heart failure (HF) are unclear. We assessed whether incident HF is associated with changes in cognitive function after accounting for pre-HF cognitive trajectories and known determinants of cognition. METHODS: This pooled cohort study included adults without HF, stroke, or dementia from six US population-based cohort studies from 1971-2019: Atherosclerosis Risk in Communities Study, Coronary Artery Risk Development in Young Adults Study, Cardiovascular Health Study, Framingham Offspring Study, Multi-Ethnic Study of Atherosclerosis, and Northern Manhattan Study. Linear mixed-effects models estimated changes in cognition at the time of HF (change in the intercept) and the rate of cognitive change over the years after HF (change in the slope), controlling for pre-HF cognitive trajectories and participant factors. Change in global cognition was the primary outcome. Change in executive function and memory were secondary outcomes. Cognitive outcomes were standardized to a -score metric (mean [SD], 50 [10]); a 1-point difference represented a 0.1-SD difference in cognition. RESULTS: The study included 29,614 adults (mean [SD] age was 61.1 [10.5] years, 55% female, 70.3% White, 22.2% Black 7.5% Hispanic). During a median follow-up of 6.6 (Q1-Q3: 5-19.8) years, 1,407 (4.7%) adults developed incident HF. Incident HF was associated with an acute decrease in global cognition (-1.08 points; 95% CI -1.36, -0.80) and executive function (-0.65 points; 95% CI -0.96, -0.34) but not memory (-0.51 points; 95% CI -1.37, 0.35) at the time of the event. Greater acute decreases in global cognition after HF were seen in those with older age, female sex and White race. Individuals with incident HF, compared to HF-free individuals, demonstrated faster declines in global cognition (-0.15 points per year; 95% CI, -0.21, -0.09) and executive function (-0.16 points per year; 95% CI -0.23, -0.09) but not memory ( -0.11 points per year; 95% CI -0.26, 0.04) compared with pre-HF slopes. CONCLUSIONS: In this pooled cohort study, incident HF was associated with an acute decrease in global cognition and executive function at the time of the event and faster declines in global cognition and executive function over the following years.

X-chromosome and kidney function: evidence from a multi-trait genetic analysis of 908,697 individuals reveals sex-specific and sex-differential findings in genes regulated by androgen response elements. - Scholz, Markus; Horn, Katrin; Pott, Janne; Wuttke, Matthias; Kühnapfel, Andreas; Nasr, M Kamal; Kirsten, Holger; Li, Yong; Hoppmann, Anselm; Gorski, Mathias; Ghasemi, Sahar; Li, Man; et al. (Nat Commun, 2024)
X-chromosomal genetic variants are understudied but can yield valuable insights into sexually dimorphic human traits and diseases. We performed a sex-stratified cross-ancestry X-chromosome-wide association meta-analysis of seven kidney-related traits (n = 908,697), identifying 23 loci genome-wide significantly associated with two of the traits: 7 for uric acid and 16 for estimated glomerular filtration rate (eGFR), including four novel eGFR loci containing the functionally plausible prioritized genes ACSL4, CLDN2, TSPAN6 and the female-specific DRP2. Further, we identified five novel sex-interactions, comprising male-specific effects at FAM9B and AR/EDA2R, and three sex-differential findings with larger genetic effect sizes in males at DCAF12L1 and MST4 and larger effect sizes in females at HPRT1. All prioritized genes in loci showing significant sex-interactions were located next to androgen response elements (ARE). Five ARE genes showed sex-differential expressions. This study contributes new insights into sex-dimorphisms of kidney traits along with new prioritized gene targets for further molecular research.
Keywords: Androgens; Chromosomes, Human, X; Female; Genetic Predisposition to Disease; Genome-Wide Association Study; Humans; Kidney; Male; Polymorphism, Single Nucleotide; Response Elements; Tetraspanins