Carriers of rare damaging genetic variants are at lower risk of atherosclerotic disease. - Georgakis, Marios K; Malik, Rainer; Hasbani, Natalie R; Shakt, Gabrielle; Morrison, Alanna C; Tsao, Noah L; Judy, Renae; Mitchell, Braxton D; Xu, Huichun; Montasser, May E; Do, Ron; Kenny, Eimear E; et al. (medRxiv, 2023)
BACKGROUND: The CCL2/CCR2 axis governs monocyte trafficking and recruitment to atherosclerotic lesions. Human genetic analyses and population-based studies support an association between circulating CCL2 levels and atherosclerosis. Still, it remains unknown whether pharmacological targeting of CCR2, the main CCL2 receptor, would provide protection against human atherosclerotic disease. METHODS: In whole-exome sequencing data from 454,775 UK Biobank participants (40-69 years), we identified predicted loss-of-function (LoF) or damaging missense (REVEL score >0.5) variants within the gene. We prioritized variants associated with lower monocyte count (p<0.05) and tested associations with vascular risk factors and risk of atherosclerotic disease over a mean follow-up of 14 years. The results were replicated in a pooled cohort of three independent datasets (TOPMed, deCODE and Penn Medicine BioBank; total n=441,445) and the effect of the most frequent damaging variant was experimentally validated. RESULTS: A total of 45 predicted LoF or damaging missense variants were identified in the gene, 4 of which were also significantly associated with lower monocyte count, but not with other white blood cell counts. Heterozygous carriers of these variants were at a lower risk of a combined atherosclerosis outcome, showed a lower burden of atherosclerosis across four vascular beds, and were at a lower lifetime risk of coronary artery disease and myocardial infarction. There was no evidence of association with vascular risk factors including LDL-cholesterol, blood pressure, glycemic status, or C-reactive protein. Using a cAMP assay, we found that cells transfected with the most frequent damaging variant (3:46358273:T:A, M249K, 547 carriers, frequency: 0.14%) show a decrease in signaling in response to CCL2. The associations of the M249K variant with myocardial infarction were consistent across cohorts (OR : 0.62 95%CI: 0.39-0.96; OR : 0.64 95%CI: 0.34-1.19; OR : 0.64 95%CI: 0.45-0.90). In a phenome-wide association study, we found no evidence for higher risk of common infections or mortality among carriers of damaging variants. CONCLUSIONS: Heterozygous carriers of damaging variants have a lower burden of atherosclerosis and lower lifetime risk of myocardial infarction. In conjunction with previous evidence from experimental and epidemiological studies, our findings highlight the translational potential of CCR2-targeting as an atheroprotective approach.

Clonal hematopoiesis is associated with protection from Alzheimer's disease. - Bouzid, Hind; Belk, Julia A; Jan, Max; Qi, Yanyan; Sarnowski, Chloe; Wirth, Sara; Ma, Lisa; Chrostek, Matthew R; Ahmad, Herra; Nachun, Daniel; Yao, Winnie; Beiser, Alexa; et al. (Nat Med, 2023)
Clonal hematopoiesis of indeterminate potential (CHIP) is a premalignant expansion of mutated hematopoietic stem cells. As CHIP-associated mutations are known to alter the development and function of myeloid cells, we hypothesized that CHIP may also be associated with the risk of Alzheimer's disease (AD), a disease in which brain-resident myeloid cells are thought to have a major role. To perform association tests between CHIP and AD dementia, we analyzed blood DNA sequencing data from 1,362 individuals with AD and 4,368 individuals without AD. Individuals with CHIP had a lower risk of AD dementia (meta-analysis odds ratio (OR) = 0.64, P = 3.8 × 10), and Mendelian randomization analyses supported a potential causal association. We observed that the same mutations found in blood were also detected in microglia-enriched fraction of the brain in seven of eight CHIP carriers. Single-nucleus chromatin accessibility profiling of brain-derived nuclei in six CHIP carriers revealed that the mutated cells comprised a large proportion of the microglial pool in the samples examined. While additional studies are required to validate the mechanistic findings, these results suggest that CHIP may have a role in attenuating the risk of AD.

Clonal Hematopoiesis of Indeterminate Potential (CHIP) and Incident Type 2 Diabetes Risk. - Tobias, Deirdre K; Manning, Alisa K; Wessel, Jennifer; Raghavan, Sridharan; Westerman, Kenneth E; Bick, Alexander G; DiCorpo, Daniel; Whitsel, Eric A; Collins, Jason; Correa, Adolfo; Cupples, L Adrienne; Dupuis, Josée; et al. (Diabetes Care, 2023)
OBJECTIVE: Clonal hematopoiesis of indeterminate potential (CHIP) is an aging-related accumulation of somatic mutations in hematopoietic stem cells, leading to clonal expansion. CHIP presence has been implicated in atherosclerotic coronary heart disease (CHD) and all-cause mortality, but its association with incident type 2 diabetes (T2D) is unknown. We hypothesized that CHIP is associated with elevated risk of T2D. RESEARCH DESIGN AND METHODS: CHIP was derived from whole-genome sequencing of blood DNA in the National Heart, Lung, and Blood Institute Trans-Omics for Precision Medicine (TOPMed) prospective cohorts. We performed analysis for 17,637 participants from six cohorts, without prior T2D, cardiovascular disease, or cancer. We evaluated baseline CHIP versus no CHIP prevalence with incident T2D, including associations with DNMT3A, TET2, ASXL1, JAK2, and TP53 variants. We estimated multivariable-adjusted hazard ratios (HRs) and 95% confidence intervals (CIs) with adjustment for age, sex, BMI, smoking, alcohol, education, self-reported race/ethnicity, and combined cohorts' estimates via fixed-effects meta-analysis. RESULTS: Mean (SD) age was 63.4 (11.5) years, 76% were female, and CHIP prevalence was 6.0% (n = 1,055) at baseline. T2D was diagnosed in n = 2,467 over mean follow-up of 9.8 years. Participants with CHIP had 23% (CI = 1.04, 1.45) higher risk of T2D than those with no CHIP. Specifically, higher risk was for TET2 (HR 1.48; CI = 1.05, 2.08) and ASXL1 (HR 1.76; CI = 1.03, 2.99) mutations; DNMT3A was nonsignificant (HR 1.15; CI = 0.93, 1.43). Statistical power was limited for JAK2 and TP53 analyses. CONCLUSIONS: CHIP was associated with higher incidence of T2D. CHIP mutations located on genes implicated in CHD and mortality were also related to T2D, suggesting shared aging-related pathology.

Clonal Hematopoiesis of Indeterminate Potential is Associated with Acute Kidney Injury. - Vlasschaert, Caitlyn; Robinson-Cohen, Cassianne; Kestenbaum, Bryan; Silver, Samuel A; Chen, Jian-Chun; Akwo, Elvis; Bhatraju, Pavan K; Zhang, Ming-Zhi; Cao, Shirong; Jiang, Ming; Wang, Yinqiu; Niu, Aolei; et al. (medRxiv, 2023)
Age is a predominant risk factor for acute kidney injury (AKI), yet the biological mechanisms underlying this risk are largely unknown and to date no genetic mechanisms for AKI have been established. Clonal hematopoiesis of indeterminate potential (CHIP) is a recently recognized biological mechanism conferring risk of several chronic aging diseases including cardiovascular disease, pulmonary disease and liver disease. In CHIP, blood stem cells acquire mutations in myeloid cancer driver genes such as and and the myeloid progeny of these mutated cells contribute to end-organ damage through inflammatory dysregulation. We sought to establish whether CHIP causes acute kidney injury (AKI). To address this question, we first evaluated associations with incident AKI events in three population-based epidemiology cohorts (N = 442,153). We found that CHIP was associated with a greater risk of AKI (adjusted HR 1.26, 95% CI: 1.19-1.34, p<0.0001), which was more pronounced in patients with AKI requiring dialysis (adjusted HR 1.65, 95% CI: 1.24-2.20, p=0.001). The risk was particularly high in the subset of individuals where CHIP was driven by mutations in genes other than (HR: 1.49, 95% CI: 1.37-1.61, p<0.0001). We then examined the association between CHIP and recovery from AKI in the ASSESS-AKI cohort and identified that non- CHIP was more common among those with a non-resolving pattern of injury (HR 2.3, 95% CI: 1.14-4.64, p = 0.03). To gain mechanistic insight, we evaluated the role of -CHIP to AKI in ischemia-reperfusion injury (IRI) and unilateral ureteral obstruction (UUO) mouse models. In both models, we observed more severe AKI and greater post-AKI kidney fibrosis in -CHIP mice. Kidney macrophage infiltration was markedly increased in -CHIP mice and -CHIP mutant renal macrophages displayed greater pro-inflammatory responses. In summary, this work establishes CHIP as a genetic mechanism conferring risk of AKI and impaired kidney function recovery following AKI via an aberrant inflammatory response in CHIP derived renal macrophages.

Determinants of mosaic chromosomal alteration fitness. - Pershad, Yash; Mack, Taralynn; Poisner, Hannah; Jakubek, Yasminka A; Stilp, Adrienne M; Mitchell, Braxton D; Lewis, Joshua P; Boerwinkle, Eric; Loos, Ruth J; Chami, Nathalie; Wang, Zhe; Barnes, Kathleen; et al. (medRxiv, 2023)
Clonal hematopoiesis (CH) is characterized by the acquisition of a somatic mutation in a hematopoietic stem cell that results in a clonal expansion. These driver mutations can be single nucleotide variants in cancer driver genes or larger structural rearrangements called mosaic chromosomal alterations (mCAs). The factors that influence the variations in mCA fitness and ultimately result in different clonal expansion rates are not well-understood. We used the Passenger-Approximated Clonal Expansion Rate (PACER) method to estimate clonal expansion rate for 6,381 individuals in the NHLBI TOPMed cohort with gain, loss, and copy-neutral loss of heterozygosity mCAs. Our estimates of mCA fitness were correlated (R = 0.49) with an alternative approach that estimated fitness of mCAs in the UK Biobank using a theoretical probability distribution. Individuals with lymphoid-associated mCAs had a significantly higher white blood cell count and faster clonal expansion rate. In a cross-sectional analysis, genome-wide association study of estimates of mCA expansion rate identified , , and locus variants as modulators of mCA clonal expansion rate.

Evaluating the use of blood pressure polygenic risk scores across race/ethnic background groups. - Kurniansyah, Nuzulul; Goodman, Matthew O; Khan, Alyna T; Wang, Jiongming; Feofanova, Elena; Bis, Joshua C; Wiggins, Kerri L; Huffman, Jennifer E; Kelly, Tanika; Elfassy, Tali; Guo, Xiuqing; Palmas, Walter; et al. (Nat Commun, 2023)
We assess performance and limitations of polygenic risk scores (PRSs) for multiple blood pressure (BP) phenotypes in diverse population groups. We compare "clumping-and-thresholding" (PRSice2) and LD-based (LDPred2) methods to construct PRSs from each of multiple GWAS, as well as multi-PRS approaches that sum PRSs with and without weights, including PRS-CSx. We use datasets from the MGB Biobank, TOPMed study, UK biobank, and from All of Us to train, assess, and validate PRSs in groups defined by self-reported race/ethnic background (Asian, Black, Hispanic/Latino, and White). For both SBP and DBP, the PRS-CSx based PRS, constructed as a weighted sum of PRSs developed from multiple independent GWAS, perform best across all race/ethnic backgrounds. Stratified analysis in All of Us shows that PRSs are better predictive of BP in females compared to males, individuals without obesity, and middle-aged (40-60 years) compared to older and younger individuals.
Keywords: Blood Pressure; Ethnicity; Female; Genetic Predisposition to Disease; Genome-Wide Association Study; Humans; Male; Multifactorial Inheritance; Population Health; Risk Factors

Evaluation of plasma sphingolipids as mediators of the relationship between kidney disease and cardiovascular events. - Lidgard, Benjamin; Bansal, Nisha; Zelnick, Leila R; Hoofnagle, Andrew N; Fretts, Amanda M; Longstreth, William T; Shlipak, Michael G; Siscovick, David S; Umans, Jason G; Lemaitre, Rozenn N (EBioMedicine, 2023)
BACKGROUND: Sphingolipids are a family of circulating lipids with regulatory and signaling roles that are strongly associated with both eGFR and cardiovascular disease. Patients with chronic kidney disease (CKD) are at high risk for cardiovascular events, and have different plasma concentrations of certain plasma sphingolipids compared to patients with normal kidney function. We hypothesize that circulating sphingolipids partially mediate the associations between eGFR and cardiovascular events. METHODS: We measured the circulating concentrations of 8 sphingolipids, including 4 ceramides and 4 sphingomyelins with the fatty acids 16:0, 20:0, 22:0, and 24:0, in plasma from 3,463 participants in a population-based cohort (Cardiovascular Health Study) without prevalent cardiovascular disease. We tested the adjusted mediation effects by these sphingolipids of the associations between eGFR and incident cardiovascular disease via quasi-Bayesian Monte Carlo method with 2,000 simulations, using a Bonferroni correction for significance. FINDINGS: The mean (±SD) eGFR was 70 (±16) mL/min/1.73 m; 62% of participants were women. Lower eGFR was associated with higher plasma ceramide-16:0 and sphingomyelin-16:0, and lower ceramides and sphingomyelins-20:0 and -22:0. Lower eGFR was associated with risk of incident heart failure and ischemic stroke, but not myocardial infarction. Five of eight sphingolipids partially mediated the association between eGFR and heart failure. The sphingolipids associated with the greatest proportion mediated were ceramide-16:0 (proportion mediated 13%, 95% CI 8-22%) and sphingomyelin-16:0 (proportion mediated 10%, 95% CI 5-17%). No sphingolipids mediated the association between eGFR and ischemic stroke. INTERPRETATION: Plasma sphingolipids partially mediated the association between lower eGFR and incident heart failure. Altered sphingolipids metabolism may be a novel mechanism for heart failure in patients with CKD. FUNDING: This study was supported by T32 DK007467 and a KidneyCure Ben J. Lipps Research Fellowship (Dr. Lidgard). Sphingolipid measurements were supported by R01 HL128575 (Dr. Lemaitre) and R01 HL111375 (Dr. Hoofnagle) from the National Heart, Lung, and Blood Institute (NHLBI).

Factors Associated With Circulating Sex Hormones in Men : Individual Participant Data Meta-analyses. - Marriott, Ross J; Murray, Kevin; Adams, Robert J; Antonio, Leen; Ballantyne, Christie M; Bauer, Douglas C; Bhasin, Shalender; Biggs, Mary L; Cawthon, Peggy M; Couper, David J; Dobs, Adrian S; Flicker, Leon; et al. (Ann Intern Med, 2023)
BACKGROUND: Various factors modulate circulating testosterone in men, affecting interpretation of testosterone measurements. PURPOSE: To clarify factors associated with variations in sex hormone concentrations. DATA SOURCES: Systematic literature searches (to July 2019). STUDY SELECTION: Prospective cohort studies of community-dwelling men with total testosterone measured using mass spectrometry. DATA EXTRACTION: Individual participant data (IPD) (9 studies;  = 21 074) and aggregate data (2 studies;  = 4075). Sociodemographic, lifestyle, and health factors and concentrations of total testosterone, sex hormone-binding globulin (SHBG), luteinizing hormone (LH), dihydrotestosterone, and estradiol were extracted. DATA SYNTHESIS: Two-stage random-effects IPD meta-analyses found a nonlinear association of testosterone with age, with negligible change among men aged 17 to 70 years (change per SD increase about the midpoint, -0.27 nmol/L [-7.8 ng/dL] [CI, -0.71 to 0.18 nmol/L {-20.5 to 5.2 ng/dL}]) and decreasing testosterone levels with age for men older than 70 years (-1.55 nmol/L [-44.7 ng/dL] [CI, -2.05 to -1.06 nmol/L {-59.1 to -30.6 ng/dL}]). Testosterone was inversely associated with body mass index (BMI) (change per SD increase, -2.42 nmol/L [-69.7 ng/dL] [CI, -2.70 to -2.13 nmol/L {-77.8 to -61.4 ng/dL}]). Testosterone concentrations were lower for men who were married (mean difference, -0.57 nmol/L [-16.4 ng/dL] [CI, -0.89 to -0.26 nmol/L {-25.6 to -7.5 ng/dL}]); undertook at most 75 minutes of vigorous physical activity per week (-0.51 nmol/L [-14.7 ng/dL] [CI, -0.90 to -0.13 nmol/L {-25.9 to -3.7 ng/dL}]); were former smokers (-0.34 nmol/L [-9.8 ng/dL] [CI, -0.55 to -0.12 nmol/L {-15.9 to -3.5 ng/dL}]); or had hypertension (-0.53 nmol/L [-15.3 ng/dL] [CI, -0.82 to -0.24 nmol/L {-23.6 to -6.9 ng/dL}]), cardiovascular disease (-0.35 nmol/L [-10.1 ng/dL] [CI, -0.55 to -0.15 nmol/L {-15.9 to -4.3 ng/dL}]), cancer (-1.39 nmol/L [-40.1 ng/dL] [CI, -1.79 to -0.99 nmol/L {-51.6 to -28.5 ng/dL}]), or diabetes (-1.43 nmol/L [-41.2 ng/dL] [CI, -1.65 to -1.22 nmol/L {-47.6 to -35.2 ng/dL}]). Sex hormone-binding globulin was directly associated with age and inversely associated with BMI. Luteinizing hormone was directly associated with age in men older than 70 years. LIMITATION: Cross-sectional analysis, heterogeneity between studies and in timing of blood sampling, and imputation for missing data. CONCLUSION: Multiple factors are associated with variation in male testosterone, SHBG, and LH concentrations. Reduced testosterone and increased LH concentrations may indicate impaired testicular function after age 70 years. Interpretation of individual testosterone measurements should account particularly for age older than 70 years, obesity, diabetes, and cancer. PRIMARY FUNDING SOURCE: Medical Research Future Fund, Government of Western Australia, and Lawley Pharmaceuticals. (PROSPERO: CRD42019139668).

Gene-educational attainment interactions in a multi-population genome-wide meta-analysis identify novel lipid loci. - de Las Fuentes, Lisa; Schwander, Karen L; Brown, Michael R; Bentley, Amy R; Winkler, Thomas W; Sung, Yun Ju; Munroe, Patricia B; Miller, Clint L; Aschard, Hugo; Aslibekyan, Stella; Bartz, Traci M; Bielak, Lawrence F; et al. (Front Genet, 2023)
Educational attainment, widely used in epidemiologic studies as a surrogate for socioeconomic status, is a predictor of cardiovascular health outcomes. A two-stage genome-wide meta-analysis of low-density lipoprotein cholesterol (LDL), high-density lipoprotein cholesterol (HDL), and triglyceride (TG) levels was performed while accounting for gene-educational attainment interactions in up to 226,315 individuals from five population groups. We considered two educational attainment variables: "Some College" (yes/no, for any education beyond high school) and "Graduated College" (yes/no, for completing a 4-year college degree). Genome-wide significant ( < 5 × 10) and suggestive ( < 1 × 10) variants were identified in Stage 1 (in up to 108,784 individuals) through genome-wide analysis, and those variants were followed up in Stage 2 studies (in up to 117,531 individuals). In combined analysis of Stages 1 and 2, we identified 18 novel lipid loci (nine for LDL, seven for HDL, and two for TG) by two degree-of-freedom (2 DF) joint tests of main and interaction effects. Four loci showed significant interaction with educational attainment. Two loci were significant only in cross-population analyses. Several loci include genes with known or suggested roles in adipose (), brain (), and liver () biology, highlighting the potential importance of brain-adipose-liver communication in the regulation of lipid metabolism. An investigation of the potential druggability of genes in identified loci resulted in five gene targets shown to interact with drugs approved by the Food and Drug Administration, including genes with roles in adipose and brain tissue. Genome-wide interaction analysis of educational attainment identified novel lipid loci not previously detected by analyses limited to main genetic effects.

Genetic architecture of spatial electrical biomarkers for cardiac arrhythmia and relationship with cardiovascular disease. - Young, William J; Haessler, Jeffrey; Benjamins, Jan-Walter; Repetto, Linda; Yao, Jie; Isaacs, Aaron; Harper, Andrew R; Ramirez, Julia; Garnier, Sophie; Van Duijvenboden, Stefan; Baldassari, Antoine R; Concas, Maria Pina; et al. (Nat Commun, 2023)
The 3-dimensional spatial and 2-dimensional frontal QRS-T angles are measures derived from the vectorcardiogram. They are independent risk predictors for arrhythmia, but the underlying biology is unknown. Using multi-ancestry genome-wide association studies we identify 61 (58 previously unreported) loci for the spatial QRS-T angle (N = 118,780) and 11 for the frontal QRS-T angle (N = 159,715). Seven out of the 61 spatial QRS-T angle loci have not been reported for other electrocardiographic measures. Enrichments are observed in pathways related to cardiac and vascular development, muscle contraction, and hypertrophy. Pairwise genome-wide association studies with classical ECG traits identify shared genetic influences with PR interval and QRS duration. Phenome-wide scanning indicate associations with atrial fibrillation, atrioventricular block and arterial embolism and genetically determined QRS-T angle measures are associated with fascicular and bundle branch block (and also atrioventricular block for the frontal QRS-T angle). We identify potential biology involved in the QRS-T angle and their genetic relationships with cardiovascular traits and diseases, may inform future research and risk prediction.
Keywords: Arrhythmias, Cardiac; Atrioventricular Block; Biomarkers; Cardiovascular Diseases; Electrocardiography; Genome-Wide Association Study; Humans; Risk Factors

Genetic control of mRNA splicing as a potential mechanism for incomplete penetrance of rare coding variants. - Einson, Jonah; Glinos, Dafni; Boerwinkle, Eric; Castaldi, Peter; Darbar, Dawood; de Andrade, Mariza; Ellinor, Patrick; Fornage, Myriam; Gabriel, Stacey; Germer, Soren; Gibbs, Richard; Hersh, Craig P; et al. (bioRxiv, 2023)
Exonic variants present some of the strongest links between genotype and phenotype. However, these variants can have significant inter-individual pathogenicity differences, known as variable penetrance. In this study, we propose a model where genetically controlled mRNA splicing modulates the pathogenicity of exonic variants. By first cataloging exonic inclusion from RNA-seq data in GTEx v8, we find that pathogenic alleles are depleted on highly included exons. Using a large-scale phased WGS data from the TOPMed consortium, we observe that this effect may be driven by common splice-regulatory genetic variants, and that natural selection acts on haplotype configurations that reduce the transcript inclusion of putatively pathogenic variants, especially when limiting to haploinsufficient genes. Finally, we test if this effect may be relevant for autism risk using families from the Simons Simplex Collection, but find that splicing of pathogenic alleles has a penetrance reducing effect here as well. Overall, our results indicate that common splice-regulatory variants may play a role in reducing the damaging effects of rare exonic variants.

The genetic determinants of recurrent somatic mutations in 43,693 blood genomes. - Weinstock, Joshua S; Laurie, Cecelia A; Broome, Jai G; Taylor, Kent D; Guo, Xiuqing; Shuldiner, Alan R; O'Connell, Jeffrey R; Lewis, Joshua P; Boerwinkle, Eric; Barnes, Kathleen C; Chami, Nathalie; Kenny, Eimear E; et al. (Sci Adv, 2023)
Nononcogenic somatic mutations are thought to be uncommon and inconsequential. To test this, we analyzed 43,693 National Heart, Lung and Blood Institute Trans-Omics for Precision Medicine blood whole genomes from 37 cohorts and identified 7131 non-missense somatic mutations that are recurrently mutated in at least 50 individuals. These recurrent non-missense somatic mutations (RNMSMs) are not clearly explained by other clonal phenomena such as clonal hematopoiesis. RNMSM prevalence increased with age, with an average 50-year-old having 27 RNMSMs. Inherited germline variation associated with RNMSM acquisition. These variants were found in genes involved in adaptive immune function, proinflammatory cytokine production, and lymphoid lineage commitment. In addition, the presence of eight specific RNMSMs associated with blood cell traits at effect sizes comparable to Mendelian genetic mutations. Overall, we found that somatic mutations in blood are an unexpectedly common phenomenon with ancestry-specific determinants and human health consequences.
Keywords: Germ-Line Mutation; Hematopoiesis; Humans; Middle Aged; Mutation; Mutation, Missense; Phenotype

Genome-Wide Interaction Analysis with DASH Diet Score Identified Novel Loci for Systolic Blood Pressure. - Guirette, Melanie; Lan, Jessie; McKeown, Nicola; Brown, Michael R; Chen, Han; de Vries, Paul S; Kim, Hyunju; Rebholz, Casey M; Morrison, Alanna C; Bartz, Traci M; Fretts, Amanda M; Guo, Xiuqing; et al. (medRxiv, 2023)
OBJECTIVE: We examined interactions between genotype and a Dietary Approaches to Stop Hypertension (DASH) diet score in relation to systolic blood pressure (SBP). METHODS: We analyzed up to 9,420,585 biallelic imputed single nucleotide polymorphisms (SNPs) in up to 127,282 individuals of six population groups (91% of European population) from the Cohorts for Heart and Aging Research in Genomic Epidemiology consortium (CHARGE; n=35,660) and UK Biobank (n=91,622) and performed European population-specific and cross-population meta-analyses. RESULTS: We identified three loci in European-specific analyses and an additional four loci in cross-population analyses at P for interaction < 5e-8. We observed a consistent interaction between rs117878928 at 15q25.1 (minor allele frequency = 0.03) and the DASH diet score (P for interaction = 4e-8; P for heterogeneity = 0.35) in European population, where the interaction effect size was 0.42±0.09 mm Hg (P for interaction = 9.4e-7) and 0.20±0.06 mm Hg (P for interaction = 0.001) in CHARGE and the UK Biobank, respectively. The 1 Mb region surrounding rs117878928 was enriched with -expression quantitative trait loci (eQTL) variants (P = 4e-273) and -DNA methylation quantitative trait loci (mQTL) variants (P = 1e-300). While the closest gene for rs117878928 is , the highest narrow sense heritability accounted by SNPs potentially interacting with the DASH diet score in this locus was for gene at 15q25.1. CONCLUSION: We demonstrated gene-DASH diet score interaction effects on SBP in several loci. Studies with larger diverse populations are needed to validate our findings.

{Identification of circulating proteins associated with general cognitive function among middle-aged and older adults - Tin, A.; Fohner, A. E.; Yang, Q.; Brody, J. A.; Davies, G.; Yao, J.; Liu, D.; Caro, I.; Lindbohm, J. V.; Duggan, M. R.; Meirelles, O.; Harris, S. E.; et al. (Commun Biol, 2023)
2.0E-4). Proteins implicated as causes or consequences of AD susceptibility may provide new insight into the potential relationship between immunity and AD susceptibility as well as potential therapeutic targets.

Inflammation and Incident Conduction Disease. - Frimodt-Møller, Emilie K; Gottdiener, John S; Soliman, Elsayed Z; Kizer, Jorge R; Vittinghoff, Eric; Psaty, Bruce M; Biering-Sørensen, Tor; Marcus, Gregory M (J Am Heart Assoc, 2023)
Keywords: Cardiac Conduction System Disease; Electrocardiography; Heart Block; Humans; Inflammation

Investigation of the independent role of a rare APOE variant (L28P; APOE*4Pittsburgh) in late-onset Alzheimer disease. - Fan, KangHsien; Francis, Lily; Aslam, M Muaaz; Bedison, Margret A; Lawrence, Elizabeth; Acharya, Vibha; Snitz, Beth E; Ganguli, Mary; DeKosky, Steven T; Lopez, Oscar L; Feingold, Eleanor; Kamboh, M Ilyas (Neurobiol Aging, 2023)
A rare missense APOE variant (L28P; APOE*4Pittsburgh), which is present only in populations with European ancestry, has been reported to be a risk factor for late-onset Alzheimer's disease (LOAD). However, due to the complete linkage disequilibrium of L28P with APOE*4 (C112R), its independent genetic association is uncertain. The original association study implicating L28P with LOAD risk was carried out in a relatively small sample size. In the current study, we have re-evaluated this association in a large case-control sample of 15,762 White U.S. subjects and investigated its independent effect in APOE 3/4 subjects, as L28P has been observed only in the heterozygous state of APOE*4 carriers and 3/4 is the most common genotype containing the APOE*4 allele. The heterozygous carrier frequency of L28P, all with APOE*4, was about 3-fold higher in AD cases than in cognitively intact controls (0.845% vs. 0.277%). The age- and sex-adjusted meta-analysis odds ratio (OR) was 2.87 (95% CI: 1.34 - 6.13; = 0.0066). Among APOE 3/4 subjects, age- and sex-adjusted meta-analysis OR was 1.53 (95% CI: 0.70 - 3.36; p = 0.28), indicating its effect was independent of APOE*4. The lack of statistical significance appears mainly due to the low power of 4138 subjects with the 3/4 genotype (12% power at α= 0.05) compared to the required sample of 139,088 subjects with the 3/4 genotype to detect an OR of 1.5 at α= 0.05 and 80% power. Our data suggesting that L28P has an independent genetic effect on AD risk is reinforced by earlier experimental findings showing that this mutation leads to significant structural and conformational changes in the ApoE4 molecule and can induce functional defects associated with neuronal Aβ42 accumulation and oxidative stress. Additional functional studies in cell-based systems and animal models will help to delineate its functional significance in the etiology of AD.
Keywords: Alleles; Alzheimer Disease; Animals; Apolipoprotein E3; Apolipoprotein E4; Apolipoproteins E; Genotype; Heterozygote

Lifestyle habits associated with cardiac conduction disease. - Frimodt-Møller, Emilie K; Soliman, Elsayed Z; Kizer, Jorge R; Vittinghoff, Eric; Psaty, Bruce M; Biering-Sørensen, Tor; Gottdiener, John S; Marcus, Gregory M (Eur Heart J, 2023)
AIMS: Cardiac conduction disease can lead to syncope, heart failure, and death. The only available therapy is pacemaker implantation, with no established prevention strategies. Research to identify modifiable risk factors has been scant. METHODS AND RESULTS: Data from the Cardiovascular Health Study, a population-based cohort study of adults ≥ 65 years with annual 12-lead electrocardiograms obtained over 10 years, were utilized to examine relationships between baseline characteristics, including lifestyle habits, and conduction disease. Of 5050 participants (mean age 73 ± 6 years; 52% women), prevalent conduction disease included 257 with first-degree atrioventricular block, 99 with left anterior fascicular block, 9 with left posterior fascicular block, 193 with right bundle branch block (BBB), 76 with left BBB, and 102 with intraventricular block at baseline. After multivariable adjustment, older age, male sex, a larger body mass index, hypertension, and coronary heart disease were associated with a higher prevalence of conduction disease, whereas White race and more physical activity were associated with a lower prevalence. Over a median follow-up on 7 (interquartile range 1-9) years, 1036 developed incident conduction disease. Older age, male sex, a larger BMI, and diabetes were each associated with incident conduction disease. Of lifestyle habits, more physical activity (hazard ratio 0.91, 95% confidence interval 0.84-0.98, P = 0.017) was associated with a reduced risk, while smoking and alcohol did not exhibit a significant association. CONCLUSION: While some difficult to control comorbidities were associated with conduction disease as expected, a readily modifiable lifestyle factor, physical activity, was associated with a lower risk.

Machine learning models for blood pressure phenotypes combining multiple polygenic risk scores. - Hrytsenko, Yana; Shea, Benjamin; Elgart, Michael; Kurniansyah, Nuzulul; Lyons, Genevieve; Morrison, Alanna C; Carson, April P; Haring, Bernhard; Mitchel, Braxton D; Psaty, Bruce M; Jaeger, Byron C; Gu, C Charles; et al. (medRxiv, 2023)
We construct non-linear machine learning (ML) prediction models for systolic and diastolic blood pressure (SBP, DBP) using demographic and clinical variables and polygenic risk scores (PRSs). We developed a two-model ensemble, consisting of a baseline model, where prediction is based on demographic and clinical variables only, and a genetic model, where we also include PRSs. We evaluate the use of a linear versus a non-linear model at both the baseline and the genetic model levels and assess the improvement in performance when incorporating multiple PRSs. We report the ensemble model's performance as percentage variance explained (PVE) on a held-out test dataset. A non-linear baseline model improved the PVEs from 28.1% to 30.1% (SBP) and 14.3% to 17.4% (DBP) compared with a linear baseline model. Including seven PRSs in the genetic model computed based on the largest available GWAS of SBP/DBP improved the genetic model PVE from 4.8% to 5.1% (SBP) and 4.7% to 5% (DBP) compared to using a single PRS. Adding additional 14 PRSs computed based on two independent GWASs further increased the genetic model PVE to 6.3% (SBP) and 5.7% (DBP). PVE differed across self-reported race/ethnicity groups, with primarily all non-White groups benefitting from the inclusion of additional PRSs.

Mortality Following Hip Fracture in Older Adults With and Without Coronary Heart Disease. - Robbins, John A; Bůzková, Petra; Barzilay, Joshua I; Cauley, Jane A; Fink, Howard A; Carbone, Laura D; Chen, Zhao; Stein, Phyllis K; Elam, Rachel; Sheets, Kerry; Mukamal, Kenneth J (Am J Med, 2023)
BACKGROUND: Comorbidities like coronary heart disease are common among older people who sustain an osteoporotic hip fracture. However, their impact on short- and long-term mortality post-hip fracture is not well quantified. METHODS: We examined 4092 and 1173 older adults without and with prevalent coronary heart disease, respectively. Post-hip fracture mortality rates were computed with Poisson models and hazard ratios with Cox regression. For perspective, we compared mortality rates among participants with prevalent coronary heart disease who had either a hip fracture or incident heart failure (but no hip fracture). RESULTS: Among participants without prevalent coronary heart disease, the mortality rate post-hip fracture was 21.83 per 100 participant years, including 49.27 per 100 participant years in the first 6 months following hip fracture. Among participants with prevalent coronary heart disease, the corresponding mortality rates were 32.52 and 79.44 per 100 participant years, respectively. Participants with prevalent coronary heart disease and incident heart failure (but no hip fracture) had corresponding post-incident heart failure mortality rates per 100 participant years of 25.62 overall and 46.4 in the first 6 months. In all 3 groups, the hazard ratio for mortality was similarly elevated: 5- to 7-fold at 6 months and 1.7- to 2.5-fold beyond 5 years. CONCLUSION: As a case study in the absolute effects of a comorbidity on post-hip fracture mortality, hip fracture in a person with coronary heart disease carries an exceedingly high mortality rate, even higher than that following incident heart failure in individuals with coronary heart disease.

Mosaic chromosomal alterations in blood across ancestries using whole-genome sequencing. - Jakubek, Yasminka A; Zhou, Ying; Stilp, Adrienne; Bacon, Jason; Wong, Justin W; Ozcan, Zuhal; Arnett, Donna; Barnes, Kathleen; Bis, Joshua C; Boerwinkle, Eric; Brody, Jennifer A; Carson, April P; et al. (Nat Genet, 2023)
Megabase-scale mosaic chromosomal alterations (mCAs) in blood are prognostic markers for a host of human diseases. Here, to gain a better understanding of mCA rates in genetically diverse populations, we analyzed whole-genome sequencing data from 67,390 individuals from the National Heart, Lung, and Blood Institute Trans-Omics for Precision Medicine program. We observed higher sensitivity with whole-genome sequencing data, compared with array-based data, in uncovering mCAs at low mutant cell fractions and found that individuals of European ancestry have the highest rates of autosomal mCAs and the lowest rates of chromosome X mCAs, compared with individuals of African or Hispanic ancestry. Although further studies in diverse populations will be needed to replicate our findings, we report three loci associated with loss of chromosome X, associations between autosomal mCAs and rare variants in DCPS, ADM17, PPP1R16B and TET2 and ancestry-specific variants in ATM and MPL with mCAs in cis.
Keywords: Black People; Genome, Human; Genome-Wide Association Study; Hispanic or Latino; Humans; Mosaicism; Precision Medicine

Multi-ancestry genome-wide study identifies effector genes and druggable pathways for coronary artery calcification. - Kavousi, Maryam; Bos, Maxime M; Barnes, Hanna J; Lino Cardenas, Christian L; Wong, Doris; Lu, Haojie; Hodonsky, Chani J; Landsmeer, Lennart P L; Turner, Adam W; Kho, Minjung; Hasbani, Natalie R; de Vries, Paul S; et al. (Nat Genet, 2023)
Coronary artery calcification (CAC), a measure of subclinical atherosclerosis, predicts future symptomatic coronary artery disease (CAD). Identifying genetic risk factors for CAC may point to new therapeutic avenues for prevention. Currently, there are only four known risk loci for CAC identified from genome-wide association studies (GWAS) in the general population. Here we conducted the largest multi-ancestry GWAS meta-analysis of CAC to date, which comprised 26,909 individuals of European ancestry and 8,867 individuals of African ancestry. We identified 11 independent risk loci, of which eight were new for CAC and five had not been reported for CAD. These new CAC loci are related to bone mineralization, phosphate catabolism and hormone metabolic pathways. Several new loci harbor candidate causal genes supported by multiple lines of functional evidence and are regulators of smooth muscle cell-mediated calcification ex vivo and in vitro. Together, these findings help refine the genetic architecture of CAC and extend our understanding of the biological and potential druggable pathways underlying CAC.

Multi-ancestry genome-wide study in >2.5 million individuals reveals heterogeneity in mechanistic pathways of type 2 diabetes and complications. - Suzuki, Ken; Hatzikotoulas, Konstantinos; Southam, Lorraine; Taylor, Henry J; Yin, Xianyong; Lorenz, Kim M; Mandla, Ravi; Huerta-Chagoya, Alicia; Rayner, Nigel W; Bocher, Ozvan; Ana Luiza de, S V Arruda; Sonehara, Kyuto; et al. (medRxiv, 2023)
Type 2 diabetes (T2D) is a heterogeneous disease that develops through diverse pathophysiological processes. To characterise the genetic contribution to these processes across ancestry groups, we aggregate genome-wide association study (GWAS) data from 2,535,601 individuals (39.7% non-European ancestry), including 428,452 T2D cases. We identify 1,289 independent association signals at genome-wide significance (P<5×10 ) that map to 611 loci, of which 145 loci are previously unreported. We define eight non-overlapping clusters of T2D signals characterised by distinct profiles of cardiometabolic trait associations. These clusters are differentially enriched for cell-type specific regions of open chromatin, including pancreatic islets, adipocytes, endothelial, and enteroendocrine cells. We build cluster-specific partitioned genetic risk scores (GRS) in an additional 137,559 individuals of diverse ancestry, including 10,159 T2D cases, and test their association with T2D-related vascular outcomes. Cluster-specific partitioned GRS are more strongly associated with coronary artery disease and end-stage diabetic nephropathy than an overall T2D GRS across ancestry groups, highlighting the importance of obesity-related processes in the development of vascular outcomes. Our findings demonstrate the value of integrating multi-ancestry GWAS with single-cell epigenomics to disentangle the aetiological heterogeneity driving the development and progression of T2D, which may offer a route to optimise global access to genetically-informed diabetes care.

Multi-ancestry transcriptome-wide association analyses yield insights into tobacco use biology and drug repurposing. - Chen, Fang; Wang, Xingyan; Jang, Seon-Kyeong; Quach, Bryan C; Weissenkampen, J Dylan; Khunsriraksakul, Chachrit; Yang, Lina; Sauteraud, Renan; Albert, Christine M; Allred, Nicholette D D; Arnett, Donna K; Ashley-Koch, Allison E; et al. (Nat Genet, 2023)
Most transcriptome-wide association studies (TWASs) so far focus on European ancestry and lack diversity. To overcome this limitation, we aggregated genome-wide association study (GWAS) summary statistics, whole-genome sequences and expression quantitative trait locus (eQTL) data from diverse ancestries. We developed a new approach, TESLA (multi-ancestry integrative study using an optimal linear combination of association statistics), to integrate an eQTL dataset with a multi-ancestry GWAS. By exploiting shared phenotypic effects between ancestries and accommodating potential effect heterogeneities, TESLA improves power over other TWAS methods. When applied to tobacco use phenotypes, TESLA identified 273 new genes, up to 55% more compared with alternative TWAS methods. These hits and subsequent fine mapping using TESLA point to target genes with biological relevance. In silico drug-repurposing analyses highlight several drugs with known efficacy, including dextromethorphan and galantamine, and new drugs such as muscle relaxants that may be repurposed for treating nicotine addiction.
Keywords: Biology; Drug Repositioning; Genetic Predisposition to Disease; Genome-Wide Association Study; Humans; Polymorphism, Single Nucleotide; Tobacco Use; Transcriptome

Neighborhood greenspace and cognition: The cardiovascular health study. - Godina, Sara L; Rosso, Andrea L; Hirsch, Jana A; Besser, Lilah M; Lovasi, Gina S; Donovan, Geoffrey H; Garg, Parveen K; Platt, Jonathan M; Fitzpatrick, Annette L; Lopez, Oscar L; Carlson, Michelle C; Michael, Yvonne L (Health Place, 2023)
OBJECTIVES: We examined whether greenspace measures (overall percent greenspace and forest, and number of greenspace types) were associated with clinically adjudicated dementia status. METHODS: In a sample of non-demented older adults (n = 2141, average age = 75.3 years) from the Cardiovascular Health and Cognition Study, Cox proportional hazard and logistic regression analyses were used to estimate associations of baseline greenspace with risks of incident dementia and MCI, respectively, while adjusting for demographics, co-morbidities, and other neighborhood factors. We derived quartiles of percent greenness (greenspace), forest (percent tree canopy cover), and tertiles of greenspace diversity (number of greenspace types) for 5-km radial buffers around participant's residences at study entry (1989-1990) from the 1992 National Land Cover Dataset. Dementia status and mild cognitive impairment (MCI) over 10 years was clinically adjudicated. RESULTS: We observed no significant association between overall percent greenspace and risk of mild cognitive impairment or dementia and mostly null results for forest and greenspace diversity. Forest greenspace was associated with lower odds of MCI (OR quartile 4 versus 1: 0.54, 95% CI: 0.29-0.98) and greenspace diversity was associated with lower hazard of incident dementia (HR tertile 2 versus 1: 0.70, 95% CI = 0.50-0.99). DISCUSSION: We found divergent results for different types of greenspace and mild cognitive impairment or dementia. Improved greenspace type and diversity measurement could better characterize the association between greenspace and cognition.

Plasma Ceramides and Sphingomyelins and Sudden Cardiac Death in the Cardiovascular Health Study. - Bockus, Lee B; Jensen, Paul N; Fretts, Amanda M; Hoofnagle, Andrew N; McKnight, Barbara; Sitlani, Colleen M; Siscovick, David S; King, Irena B; Psaty, Bruce M; Sotoodehnia, Nona; Lemaitre, Rozenn N (JAMA Netw Open, 2023)
IMPORTANCE: Sphingolipids, including ceramides and sphingomyelins, may influence the pathophysiology and risk of sudden cardiac death (SCD) through multiple biological activities. Whether the length of the fatty acid acylated to plasma sphingolipid species is associated with SCD risk is not known. OBJECTIVE: To determine whether the saturated fatty acid length of plasma ceramides and sphingomyelins influences the association with SCD risk. DESIGN, SETTING, AND PARTICIPANTS: In this cohort study, multivariable Cox proportional hazards regression models were used to examine the association of sphingolipid species with SCD risk. The study population included 4612 participants in the Cardiovascular Health Study followed up prospectively for a median of 10.2 (IQR, 5.5-11.6) years. Baseline data were collected from January 1992 to December 1995 during annual examinations. Data were analyzed from February 11, 2020, to September 9, 2023. EXPOSURES: Eight plasma sphingolipid species (4 ceramides and 4 sphingomyelins) with saturated fatty acids of 16, 20, 22, and 24 carbons. MAIN OUTCOME AND MEASURE: Association of plasma ceramides and sphingomyelins with saturated fatty acids of different lengths with SCD risk. RESULTS: Among the 4612 CHS participants included in the analysis (mean [SD] age, 77 [5] years; 2724 [59.1%] women; 6 [0.1%] American Indian; 4 [0.1%] Asian; 718 [15.6%] Black; 3869 [83.9%] White, and 15 [0.3%] Other), 215 SCD cases were identified. In adjusted Cox proportional hazards regression analyses, plasma ceramides and sphingomyelins with palmitic acid (Cer-16 and SM-16) were associated with higher SCD risk per higher SD of log sphingolipid levels (hazard ratio [HR] for Cer-16, 1.34 [95% CI, 1.12-1.59]; HR for SM-16, 1.37 [95% CI, 1.12-1.67]). Associations did not differ by baseline age, sex, race, or body mass index. No significant association of SCD with sphingolipids with very-long-chain saturated fatty acids was observed after correction for multiple testing (HR for ceramide with arachidic acid, 1.06 [95% CI, 0.90-1.24]; HR for ceramide with behenic acid, 0.92 [95% CI, 0.77-1.10]; HR for ceramide with lignoceric acid, 0.92 [95% CI, 0.77-1.09]; HR for sphingomyelin with arachidic acid, 0.83 [95% CI, 0.71-0.98]; HR for sphingomyelin with behenic acid, 0.84 [95% CI, 0.70-1.00]; HR for sphingomyelin with lignoceric acid, 0.86 [95% CI, 0.72-1.03]). CONCLUSIONS AND RELEVANCE: The findings of this large, population-based cohort study of SCD identified that higher plasma levels of Cer-16 and SM-16 were associated with higher risk of SCD. Future studies are needed to examine the underlying mechanism of these associations.
Keywords: Aged; Ceramides; Cohort Studies; Death, Sudden, Cardiac; Eicosanoic Acids; Fatty Acids; Female; Humans; Male; Sphingolipids; Sphingomyelins

Plasma Levels of Branched Chain Amino Acids, Incident Hip Fractures and Bone Mineral Density of the Hip and Spine. - Carbone, Laura; Bůzková, Petra; Fink, Howard A; Robbins, John A; Barzilay, Joshua I; Elam, Rachel E; Isales, Carlos; Connelly, Margery A; Mukamal, Kenneth J (J Clin Endocrinol Metab, 2023)
OBJECTIVE: Branched chain amino acids (BCAA) are building blocks for protein, an essential component of bone. However, the association of plasma levels of BCAA with fractures in populations outside of Hong Kong or with hip fractures in particular is not known. The purpose of these analyses was to determine the relationship of BCAA including valine, leucine and isoleucine and total BCAA (standard deviation of the sum of Z-scores for each BCAA) with incident hip fractures and bone mineral density (BMD) of the hip and lumbar spine in older African American and Caucasian men and women in the Cardiovascular Health Study (CHS). DESIGN: Longitudinal analyses of association of plasma levels of BCAA with incident hip fractures and cross-sectional BMD of the hip and lumbar spine from the CHS. SETTING: Community. PARTICIPANTS: 1850 men (38% of cohort) and women; mean age 73. MAIN OUTCOME MEASURES: Incident hip fractures and cross-sectional BMD of the total hip, femoral neck and lumbar spine. RESULTS: In fully adjusted models, over 12 years of follow-up, we observed no significant association between incident hip fracture and plasma values of valine, leucine, isoleucine or total BCAA per 1 standard deviation higher of each BCAA. Plasma values of leucine but not valine, isoleucine or total BCAA, were positively and significantly associated with BMD of the total hip (p = 0.03) and femoral neck (p = 0.02), but not the lumbar spine (p = 0.07). CONCLUSIONS: Plasma levels of the BCAA leucine may be associated with higher BMD in older men and women. However, given lack of a significant association with hip fracture risk, further information is needed to determine whether BCAAs would be novel targets for osteoporosis therapies.

Plasma Proteomic Associations With Incident Ischemic Stroke in Older Adults: The Cardiovascular Health Study. - Kalani, Rizwan; Bartz, Traci M; Psaty, Bruce M; Elkind, Mitchell S V; Floyd, James S; Gerszten, Robert E; Shojaie, Ali; Heckbert, Susan R; Bis, Joshua C; Austin, Thomas R; Tirschwell, David L; Delaney, Joseph A C; et al. (Neurology, 2023)
BACKGROUND: Plasma proteomics may elucidate novel insights into the pathophysiology of ischemic stroke (IS), identify biomarkers of IS risk, and guide development of nascent prevention strategies. We evaluated the relationship between the plasma proteome and IS risk in the population-based Cardiovascular Health Study (CHS). METHODS: Eligible CHS participants were free of prevalent stroke and underwent quantification of 1298 plasma proteins using the aptamer-based SOMAScan assay platform from the 1992-1993 study visit. Multivariable Cox proportional hazards regression was used to evaluate associations between a 1-standard deviation increase in the log-2 transformed estimated plasma protein concentrations and incident IS, adjusting for demographics, IS risk factors, and estimated glomerular filtration rate. For proteins independently associated with incident IS, a secondary stratified analysis evaluated associations in subgroups defined by sex and race. Exploratory analyses evaluated plasma proteomic associations with cardioembolic and non-cardioembolic IS as well as proteins associated with IS risk in participants with left atrial dysfunction but without atrial fibrillation. RESULTS: Of 2983 eligible participants, the mean age was 74.3 (± 4.8) years, 61.2% were women, and 15.4% were Black. Over a median follow-up of 12.6 years, 450 participants experienced an incident IS. N-terminal pro-brain natriuretic peptide (NTproBNP, adjusted HR 1.37, 95% CI 1.23-1.53, P=2.08x10) and macrophage metalloelastase (MMP12, adjusted HR 1.30, 95% CI 1.16-1.45, P=4.55x10) were independently associated with IS risk. These two associations were similar in men and women and in Black and non-Black participants. In exploratory analyses, NTproBNP was independently associated with incident cardioembolic IS, E-selectin with incident non-cardioembolic IS, and secreted frizzled-related protein 1 with IS risk in participants with left atrial dysfunction. CONCLUSIONS: In a cohort of older adults, NTproBNP and MMP12 were independently associated with IS risk. We identified plasma proteomic determinants of incident cardioembolic and non-cardioembolic IS and found a novel protein associated with IS risk in those with left atrial dysfunction.

Plasma Trimethylamine--Oxide and Incident Ischemic Stroke: The Cardiovascular Health Study and the Multi-Ethnic Study of Atherosclerosis. - Lemaitre, Rozenn N; Jensen, Paul N; Wang, Zeneng; Fretts, Amanda M; Sitlani, Colleen M; Nemet, Ina; Sotoodehnia, Nona; de Oliveira Otto, Marcia C; Zhu, Weifei; Budoff, Matt; Longstreth, W T; Psaty, Bruce M; et al. (J Am Heart Assoc, 2023)
Background The association of circulating trimethylamine--oxide (TMAO) with stroke has received limited attention. To address this gap, we examined the associations of serial measures of plasma TMAO with incident ischemic stroke. Methods and Results We used a prospective cohort design with data pooled from 2 cohorts. The settings were the CHS (Cardiovascular Health Study), a cohort of older adults, and the MESA (Multi-Ethnic Study of Atherosclerosis), both in the United States. We measured plasma concentrations of TMAO at baseline and again during the follow-up using high-performance liquid chromatography and mass spectrometry. We assessed the association of plasma TMAO with incident ischemic stroke using proportional hazards regression adjusted for risk factors. The combined cohorts included 11 785 participants without a history of stroke, on average 73 (CHS) and 62 (MESA) years old at baseline, including 60% (CHS) and 53% (MESA) women. We identified 1031 total incident ischemic strokes during a median 15-year follow-up in the combined cohorts. In multivariable analyses, TMAO was significantly associated with incident ischemic stroke risk (hazard ratios comparing a doubling of TMAO: 1.11 [1.03-1.18], =0.004). The association was linear over the range of TMAO concentrations and appeared restricted to those without diagnosed coronary heart disease. An association with hemorrhagic stroke was not found. Conclusions Plasma TMAO levels are associated with incident ischemic stroke in a diverse population. Registration URL: https://www.clinicaltrials.gov. Unique identifier: NCT00005133.
Keywords: Aged; Atherosclerosis; Female; Humans; Ischemic Stroke; Methylamines; Oxides; Prospective Studies; Risk Factors; Stroke; United States

Posttranslational modifications induce autoantibodies with risk prediction capability in patients with small cell lung cancer. - Lastwika, Kristin J; Kunihiro, Andrew; Solan, Joell L; Zhang, Yuzheng; Taverne, Lydia R; Shelley, David; Rho, Jung-Hyun; Randolph, Timothy W; Li, Christopher I; Grogan, Eric L; Massion, Pierre P; Fitzpatrick, Annette L; et al. (Sci Transl Med, 2023)
Small cell lung cancer (SCLC) elicits the generation of autoantibodies that result in unique paraneoplastic neurological syndromes. The mechanistic basis for the formation of such autoantibodies is largely unknown but is key to understanding their etiology. We developed a high-dimensional technique that enables detection of autoantibodies in complex with native antigens directly from patient plasma. Here, we used our platform to screen 1009 human plasma samples for 3600 autoantibody-antigen complexes, finding that plasma from patients with SCLC harbors, on average, fourfold higher disease-specific autoantibody signals compared with plasma from patients with other cancers. Across three independent SCLC cohorts, we identified a set of common but previously unknown autoantibodies that are produced in response to both intracellular and extracellular tumor antigens. We further characterized several disease-specific posttranslational modifications within extracellular proteins targeted by these autoantibodies, including citrullination, isoaspartylation, and cancer-specific glycosylation. Because most patients with SCLC have metastatic disease at diagnosis, we queried whether these autoantibodies could be used for SCLC early detection. We created a risk prediction model using five autoantibodies with an average area under the curve of 0.84 for the three cohorts that improved to 0.96 by incorporating cigarette smoke consumption in pack years. Together, our findings provide an innovative approach to identify circulating autoantibodies in SCLC with mechanistic insight into disease-specific immunogenicity and clinical utility.
Keywords: Autoantibodies; Humans; Lung Neoplasms; Protein Processing, Post-Translational; Small Cell Lung Carcinoma

Powerful, scalable and resource-efficient meta-analysis of rare variant associations in large whole genome sequencing studies. - Li, Xihao; Quick, Corbin; Zhou, Hufeng; Gaynor, Sheila M; Liu, Yaowu; Chen, Han; Selvaraj, Margaret Sunitha; Sun, Ryan; Dey, Rounak; Arnett, Donna K; Bielak, Lawrence F; Bis, Joshua C; et al. (Nat Genet, 2023)
Meta-analysis of whole genome sequencing/whole exome sequencing (WGS/WES) studies provides an attractive solution to the problem of collecting large sample sizes for discovering rare variants associated with complex phenotypes. Existing rare variant meta-analysis approaches are not scalable to biobank-scale WGS data. Here we present MetaSTAAR, a powerful and resource-efficient rare variant meta-analysis framework for large-scale WGS/WES studies. MetaSTAAR accounts for relatedness and population structure, can analyze both quantitative and dichotomous traits and boosts the power of rare variant tests by incorporating multiple variant functional annotations. Through meta-analysis of four lipid traits in 30,138 ancestrally diverse samples from 14 studies of the Trans Omics for Precision Medicine (TOPMed) Program, we show that MetaSTAAR performs rare variant meta-analysis at scale and produces results comparable to using pooled data. Additionally, we identified several conditionally significant rare variant associations with lipid traits. We further demonstrate that MetaSTAAR is scalable to biobank-scale cohorts through meta-analysis of TOPMed WGS data and UK Biobank WES data of ~200,000 samples.
Keywords: Exome Sequencing; Genome-Wide Association Study; Lipids; Phenotype; Whole Genome Sequencing

A proteomic analysis of atrial fibrillation in a prospective longitudinal cohort (AGES-Reykjavik study). - Jonmundsson, Thorarinn; Steindorsdottir, Anna E; Austin, Thomas R; Frick, Elisabet A; Axelsson, Gisli T; Launer, Lenore; Psaty, Bruce M; Loureiro, Joseph; Orth, Anthony P; Aspelund, Thor; Emilsson, Valur; Floyd, James S; et al. (Europace, 2023)
AIMS: Atrial fibrillation (AF) is associated with high risk of comorbidities and mortality. Our aim was to examine causal and predictive relationships between 4137 serum proteins and incident AF in the prospective population-based Age, Gene/Environment Susceptibility-Reykjavik (AGES-Reykjavik) study. METHODS AND RESULTS: The study included 4765 participants, of whom 1172 developed AF. Cox proportional hazards regression models were fitted for 4137 baseline protein measurements adjusting for known risk factors. Protein associations were tested for replication in the Cardiovascular Health Study (CHS). Causal relationships were examined in a bidirectional, two-sample Mendelian randomization analysis. The time-dependent area under the receiver operating characteristic curve (AUC)-statistic was examined as protein levels and an AF-polygenic risk score (PRS) were added to clinical risk models. The proteomic signature of incident AF consisted of 76 proteins, of which 63 (83%) were novel and 29 (38%) were replicated in CHS. The signature included both N-terminal prohormone of brain natriuretic peptide (NT-proBNP)-dependent (e.g. CHST15, ATP1B1, and SVEP1) and independent components (e.g. ASPN, AKR1B, and LAMA1/LAMB1/LAMC1). Nine causal candidates were identified (TAGLN, WARS, CHST15, CHMP3, COL15A1, DUSP13, MANBA, QSOX2, and SRL). The reverse causal analysis suggested that most AF-associated proteins were affected by the genetic liability to AF. N-terminal prohormone of brain natriuretic peptide improved the prediction of incident AF events close to baseline with further improvements gained by the AF-PRS at all time points. CONCLUSION: The AF proteomic signature includes biologically relevant proteins, some of which may be causal. It mainly reflects an NT-proBNP-dependent consequence of the genetic liability to AF. N-terminal prohormone of brain natriuretic peptide is a promising marker for incident AF in the short term, but risk assessment incorporating a PRS may improve long-term risk assessment.
Keywords: Atrial Fibrillation; Biomarkers; Endosomal Sorting Complexes Required for Transport; Humans; Natriuretic Peptide, Brain; Oxidoreductases Acting on Sulfur Group Donors; Peptide Fragments; Prognosis; Prospective Studies; Proteomics; Risk Factors

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.

{Serum NfL and GFAP are associated with incident dementia and dementia mortality in older adults: The cardiovascular health study - é, H. T.; Liu, X.; Odden, M. C.; Moseholm, K. F.; Seshadri, S.; Satizabal, C. L.; Lopez, O. L.; Bis, J. C.; é, L.; Fohner, A. E.; Psaty, B. M.; Tracy, R. P.; et al. (Alzheimers Dement, 2023)
Circulating neurofilament light chain (NfL) and glial fibrillary acidic protein (GFAP) have been independently associated with dementia risk. Their additive association, and their associations with dementia-specific mortality, have not been investigated.\ We associated serum NfL, GFAP, total tau ,and ubiquitin carboxyl-terminal hydrolase-L1, measured in 1712 dementia-free adults, with 19-year incident dementia and dementia-specific mortality risk, and with 3-year cognitive decline.\ 2.06 (1.60-2.67) and 9.22 (4.48-18.9). NfL was independently associated with accelerated cognitive decline.\ Circulating NfL and GFAP may, independently and jointly, provide useful clinical insight regarding dementia risk and prognosis.

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.

{Structural Variation Detection and Association Analysis of Whole-Genome-Sequence Data from 16,905 Alzheimer's Diseases Sequencing Project Subjects - Wang, H.; Dombroski, B. A.; Cheng, P. L.; Tucci, A.; Si, Y. Q.; Farrell, J. J.; Tzeng, J. Y.; Leung, Y. Y.; Malamon, J. S.; Wang, L. S.; Vardarajan, B. N.; Farrer, L. A.; et al. (medRxiv, 2023)
. We also identified 16 SVs associated with AD and 13 SVs associated with AD-related pathological/cognitive endophenotypes. Our findings demonstrate the broad impact of SVs on AD genetics.\ Alzheimer's disease, Structural variation, Copy number variation.

{Structural Variation Detection and Association Analysis of Whole-Genome-Sequence Data from 16,905 Alzheimer's Diseases Sequencing Project Subjects - Lee, W. P.; Wang, H.; Dombroski, B.; Cheng, P. L.; Tucci, A.; Si, Y. Q.; Farrell, J.; Tzeng, J. Y.; Leung, Y. Y.; Malamon, J.; Wang, L. S.; Vardarajan, B.; et al. (Res Sq, 2023)
. We also identified 16 SVs associated with AD and 13 SVs linked to AD-related pathological/cognitive endophenotypes. This study highlights the pivotal role of SVs in shaping our understanding of AD genetics.

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.

Time-Varying Food Retail and Incident Disease in the Cardiovascular Health Study. - Lovasi, Gina S; Boise, Sarah; Jogi, Siddharth; Hurvitz, Philip M; Rundle, Andrew G; Diez, Julia; Hirsch, Jana A; Fitzpatrick, Annette; Biggs, Mary L; Siscovick, David S (Am J Prev Med, 2023)
INTRODUCTION: Natural experiments can strengthen evidence linking neighborhood food retail presence to dietary intake patterns and cardiometabolic health outcomes, yet sample size and follow-up duration are typically not extensive. To complement natural experiment evidence, longitudinal data were used to estimate the impacts of neighborhood food retail presence on incident disease. METHODS: The Cardiovascular Health Study recruited adults aged 65+ years in 1989-1993. Analyses conducted in 2021-2022 included those in good baseline health, with addresses updated annually through the year of death (restricted to 91% who died during >2 decades of cohort follow-up). Baseline and annually updated presence of 2 combined food retail categories (supermarkets/produce markets and convenience/snack focused) was characterized using establishment-level data for 1-km and 5-km Euclidean buffers. Cox proportional hazards models estimated associations with time to each incident outcome (cardiovascular disease, diabetes), adjusting for individual and area-based confounders. RESULTS: Among 2,939 participants, 36% with baseline supermarket/produce market presence within 1 km had excess incident cardiovascular disease (hazard ratio=1.12; 95% CI=1.01, 1.24); the association was attenuated and no longer statistically significant after adjustment for sociodemographic characteristics. Adjusted associations were robustly null for time-varying supermarket/produce market or convenience/fast food retail presence across analyses with outcomes of cardiovascular disease or diabetes incidence. CONCLUSIONS: Food environment changes continue to be studied to provide an evidence base for policy decisions, and null findings in this longitudinal analysis add literature that casts doubt on the sufficiency of strategies targeting food retail presence alone of an elderly cohort for curtailing incident events of clinical importance.

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.

Validation of the CogDrisk Instrument as Predictive of Dementia in Four General Community-Dwelling Populations. - Kootar, S; Huque, M H; Eramudugolla, R; Rizzuto, D; Carlson, M C; Odden, M C; Lopez, O L; Qiu, C; Fratiglioni, L; Han, S D; Bennett, D A; Peters, R; et al. (J Prev Alzheimers Dis, 2023)
BACKGROUND: Lack of external validation of dementia risk tools is a major limitation for generalizability and translatability of prediction scores in clinical practice and research. OBJECTIVES: We aimed to validate a new dementia prediction risk tool called CogDrisk and a version, CogDrisk-AD for predicting Alzheimer's disease (AD) using cohort studies. DESIGN, SETTING, PARTICIPANTS AND MEASUREMENTS: Four cohort studies were identified that included majority of the dementia risk factors from the CogDrisk tool. Participants who were free of dementia at baseline were included. The predictors were component variables in the CogDrisk tool that include self-reported demographics, medical risk factors and lifestyle habits. Risk scores for Any Dementia and AD were computed and Area Under the Curve (AUC) was assessed. To examine modifiable risk factors for dementia, the CogDrisk tool was tested by excluding age and sex estimates from the model. RESULTS: The performance of the tool varied between studies. The overall AUC and 95% CI for predicting dementia was 0.77 (0.57, 0.97) for the Swedish National study on Aging and Care in Kungsholmen, 0.76 (0.70, 0.83) for the Health and Retirement Study - Aging, Demographics and Memory Study, 0.70 (0.67,0.72) for the Cardiovascular Health Study Cognition Study, and 0.66 (0.62,0.70) for the Rush Memory and Aging Project. CONCLUSIONS: The CogDrisk and CogDrisk-AD performed well in the four studies. Overall, this tool can be used to assess individualized risk factors of dementia and AD in various population settings.
Keywords: Aging; Alzheimer Disease; Cohort Studies; Dementia; Humans; Independent Living

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 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.

A comparative study of structural variant calling in WGS from Alzheimer's disease families. - Malamon, John S; Farrell, John J; Xia, Li Charlie; Dombroski, Beth A; Das, Rueben G; Way, Jessica; Kuzma, Amanda B; Valladares, Otto; Leung, Yuk Yee; Scanlon, Allison J; Lopez, Irving Antonio Barrera; Brehony, Jack; et al. (Life Sci Alliance, 2024)
Detecting structural variants (SVs) in whole-genome sequencing poses significant challenges. We present a protocol for variant calling, merging, genotyping, sensitivity analysis, and laboratory validation for generating a high-quality SV call set in whole-genome sequencing from the Alzheimer's Disease Sequencing Project comprising 578 individuals from 111 families. Employing two complementary pipelines, Scalpel and Parliament, for SV/indel calling, we assessed sensitivity through sample replicates (N = 9) with in silico variant spike-ins. We developed a novel metric, D-score, to evaluate caller specificity for deletions. The accuracy of deletions was evaluated by Sanger sequencing. We generated a high-quality call set of 152,301 deletions of diverse sizes. Sanger sequencing validated 114 of 146 detected deletions (78.1%). Scalpel excelled in accuracy for deletions ≤100 bp, whereas Parliament was optimal for deletions >900 bp. Overall, 83.0% and 72.5% of calls by Scalpel and Parliament were validated, respectively, including all 11 deletions called by both Parliament and Scalpel between 101 and 900 bp. Our flexible protocol successfully generated a high-quality deletion call set and a truth set of Sanger sequencing-validated deletions with precise breakpoints spanning 1-17,000 bp.
Keywords: Alzheimer Disease; Humans; Whole Genome Sequencing

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.

The Gut Microbial Metabolite Trimethylamine N-oxide, Incident CKD, and Kidney Function Decline. - Wang, Meng; Tang, W H Wilson; Li, Xinmin S; de Oliveira Otto, Marcia C; Lee, Yujin; Lemaitre, Rozenn N; Fretts, Amanda; Nemet, Ina; Sotoodehnia, Nona; Sitlani, Colleen M; Budoff, Matthew; DiDonato, Joseph A; et al. (J Am Soc Nephrol, 2024)
BACKGROUND: Trimethylamine N-oxide (TMAO) is a gut microbiota-derived metabolite of dietary phosphatidylcholine and carnitine. Experimentally, TMAO causes kidney injury and tubulointerstitial fibrosis. Little is known about prospective associations between TMAO and kidney outcomes, especially incident CKD. We hypothesized that higher plasma TMAO levels would be associated with higher risk of incident CKD and greater rate of kidney function decline. METHODS: We included 10,564 participants from two community-based, prospective cohorts with estimated glomerular filtration rate (eGFR) ≥ 60 mL/min/1.73m2 to assess incident CKD. TMAO was measured using targeted mass spectrometry at baseline and one follow-up visit. Creatinine and Cystatin C were measured up to 4 times during follow-up and used to compute eGFR. Incident CKD was defined as an eGFR decline ≥ 30% from baseline and a resulting eGFR<60 ml/min/1.73 m2. Time-varying Cox models assessed the association of serial TMAO measures with incident CKD, adjusting for sociodemographic, lifestyle, diet, and cardiovascular disease risk factors. Linear mixed models assessed the association with annualized eGFR change in 10,009 participants with at least one follow-up eGFR measure without exclusions for baseline eGFR levels. RESULTS: During a median follow-up of 9.4 years (interquartile range: 9.1-11.6 years), 979 incident CKD events occurred. Higher TMAO levels associated with higher risk of incident CKD (2nd to 5th vs. 1st quintile HR[95%CI]= 1.65 [1.22-2.23], 1.68 [1.26-2.25], 2.28 [1.72-3.02], and 2.24[1.68-2.98], respectively) and greater annualized eGFR decline ( 2nd to 5th vs. 1st quintile annualized eGFR change= -0.21 [-0.32, -0.09], -0.17 [-0.29, -0.05], -0.35 [-0.47, -0.22], and -0.43[-0.56, -0.30], respectively) with monotonic dose-response relationships. These associations were consistent across different racial/ethnic groups examined. The association with eGFR decline was similar to or larger than that seen for established CKD risk factors including diabetes, per 10 mmHg of higher systolic blood pressure, per 10 years of older age, and Black race. CONCLUSIONS: In community-based US adults, higher serial measures of plasma TMAO were associated with higher risk of incident CKD and greater annualized kidney function decline.

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

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.

Large-Scale Proteomics Identifies Novel Biomarkers and Circulating Risk Factors for Aortic Stenosis. - Shelbaya, Khaled; Arthur, Victoria; Yang, Yimin; Dorbala, Pranav; Buckley, Leo; Claggett, Brian; Skali, Hicham; Dufresne, Line; Yang, Ta-Yu; Engert, James C; Thanassoulis, George; Floyd, James; et al. (J Am Coll Cardiol, 2024)
BACKGROUND: Limited data exist regarding risk factors for aortic stenosis (AS). The plasma proteome is a promising phenotype for discovery of novel biomarkers and potentially causative mechanisms. OBJECTIVES: The aim of this study was to discover novel biomarkers with potentially causal associations with AS. METHODS: We measured 4,877 plasma proteins (SomaScan aptamer-affinity assay) among ARIC (Atherosclerosis Risk In Communities) study participants in mid-life (visit 3 [V3]; n = 11,430; age 60 ± 6 years) and in late-life (V5; n = 4,899; age 76 ± 5 years). We identified proteins cross-sectionally associated with aortic valve (AV) peak velocity (AVmax) and dimensionless index by echocardiography at V5 and with incident AV-related hospitalization after V3 with the use of multivariable linear and Cox proportional hazard regression. We assessed associations of candidate proteins with changes in AVmax over 6 years and with AV calcification with the use of cardiac computed tomography, replicated analysis in an independent sample, performed Mendelian randomization, and evaluated gene expression in explanted human AV tissue. RESULTS: Fifty-two proteins cross-sectionally were associated with AVmax and dimensionless index at V5 and with risk of incident AV-related hospitalization after V3. Among 3,413 participants in the Cardiovascular Health Study, 6 of those proteins were significantly associated with adjudicated moderate or severe AS, including matrix metalloproteinase 12 (MMP12), complement C1q tumor necrosis factor-related protein 1 (C1QTNF1), and growth differentiation factor-15. MMP12 was also associated with greater increase in AVmax over 6 years, greater degree of AV calcification, and greater expression in calcific compared with normal or fibrotic AV tissue. C1QTNF1 had consistent potential causal effects on both AS and AVmax according to Mendelian randomization analysis. CONCLUSIONS: These findings identify MMP12 as a potential novel circulating biomarker of AS risk and C1QTNF1 as a new putative target to prevent AS progression.
Keywords: Aged; Aged, 80 and over; Aortic Valve; Aortic Valve Stenosis; Biomarkers; Calcinosis; Humans; Matrix Metalloproteinase 12; Middle Aged; Proteomics; Risk Factors

Levels of procoagulant factors and peak thrombin generation in relation to dementia risk in older adults: The Cardiovascular Health Study. - Harrington, Laura B; Ehlert, Alexa N; Thacker, Evan L; Jenny, Nancy S; Lopez, Oscar; Cushman, Mary; Olson, Nels C; Fitzpatrick, Annette; Mukamal, Kenneth J; Jensen, Majken K (Thromb Res, 2024)
INTRODUCTION: Markers of hemostasis such as procoagulant factors and peak thrombin generation are associated with cardiovascular outcomes, but their associations with dementia risk are unclear. We aimed to evaluate prospective associations of selected procoagulant factors and peak thrombin generation with dementia risk. METHODS: We measured levels of 7 hemostatic factors (fibrinogen, factor VII coagulant activity [FVIIc], activated factor VII [FVIIa], factor VIIa-antithrombin [FVIIa-AT], factor XI antigen [FXI], peak thrombin generation, and platelet count) among participants in the Cardiovascular Health Study, a cohort of older adults free of dementia in 1992/1993 (n = 3185). Dementia was adjudicated and classified by DSM-IV criteria through 1998/1999. Cox proportional hazards models estimated hazard ratios (HRs) for any dementia associated with 1-standard deviation (SD) differences, adjusting for sociodemographic and clinical factors and APOE genotype. Secondary analyses separately evaluated the risk of vascular dementia, Alzheimer's disease, and mixed dementia. RESULTS: At baseline, participants had a median age of 73 years. Over 5.4 years of follow-up, we identified 448 dementia cases. There was no evidence of linear associations between levels of these hemostatic factors with any dementia risk (HRs per 1-SD difference ranged from 1.0 to 1.1; 95 % confidence intervals included 1.0). Results of secondary analyses by dementia subtype were similar. CONCLUSIONS: In this prospective study, there was no strong evidence of linear associations between levels of fibrinogen, FVIIc, FVIIa, FVIIa-AT, FXI, peak thrombin generation, or platelet count with dementia risk. Despite their associations with cardiovascular disease, higher levels of these biomarkers measured among older adults may not reflect dementia risk.
Keywords: Aged; Anticoagulants; Antithrombin III; Antithrombins; Dementia; Factor VIIa; Fibrinogen; Hemostatics; Humans; Prospective Studies; Thrombin

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.

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

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