Novel Polygenic Risk Score and Established Clinical Risk Factors for Risk Estimation of Aortic Stenosis

Importance

Polygenic risk scores (PRSs) have proven to be as strong as or stronger than established clinical risk factors for many cardiovascular phenotypes. Whether this is true for aortic stenosis remains unknown.

Objective

To develop a novel aortic stenosis PRS and compare its aortic stenosis risk estimation to established clinical risk factors.

Design, Setting, and Participants

This was a longitudinal cohort study using data from the Million Veteran Program (MVP; 2011-2020), UK Biobank (2006-2010), and 6 Thrombolysis in Myocardial Infarction (TIMI) trials, including DECLARE-TIMI 58 (2013-2018), FOURIER (TIMI 59; 2013-2017), PEGASUS-TIMI 54 (2010-2014), SAVOR-TIMI 53 (2010-2013), SOLID-TIMI 52 (2009-2014), and ENGAGE AF-TIMI 48 (2008-2013), which were a mix of population-based and randomized clinical trials. Individuals from UK Biobank and the MVP meeting a previously validated case/control definition for aortic stenosis were included. All individuals from TIMI trials were included unless they had a documented preexisting aortic valve replacement. Analysis took place from January 2022 to December 2023.

Exposures

PRS for aortic stenosis (developed using data from MVP and validated in UK Biobank) and other previously validated cardiovascular PRSs, defined either as a continuous variable or as low (bottom 20%), intermediate, and high (top 20%), and clinical risk factors.

Main Outcomes

Aortic stenosis (defined using International Classification of Diseases or Current Procedural Terminology codes in UK Biobank and MVP or safety event data in the TIMI trials).

Results

The median (IQR) age in MVP was 67 (57-73) years, and 135 140 of 147 104 participants (92%) were male. The median (IQR) age in the TIMI trials was 66 (54-78) years, and 45 524 of 59 866 participants (71%) were male. The best aortic stenosis PRS incorporated 5 170 041 single-nucleotide variants and was associated with aortic stenosis in both the MVP testing sample (odds ratio, 1.41; 95% CI, 1.37-1.45 per 1 SD PRS; P = 4.6 × 10-116) and TIMI trials (hazard ratio, 1.44; 95% CI, 1.27-1.62 per 1 SD PRS; P = 3.2 × 10-9). Among genetic and clinical risk factors, the aortic stenosis PRS performed comparably to most risk factors besides age, and within a given age range, the combination of clinical and genetic risk factors was additive, providing a 3- to 4-fold increased gradient of risk of aortic stenosis. However, the addition of the aortic stenosis PRS to a model including clinical risk factors only improved risk discrimination of aortic stenosis by 0.01 to 0.02 (C index in MVP: 0.78 with clinical risk factors, 0.79 with risk factors and aortic stenosis PRS; C index in TIMI: 0.71 with clinical risk factors, 0.73 with risk factors and aortic stenosis PRS).

Conclusions

This study developed and validated 1 of the first aortic stenosis PRSs. While aortic stenosis genetic risk was independent from clinical risk factors and performed comparably to all other risk factors besides age, genetic risk resulted in only a small improvement in overall aortic stenosis risk discrimination beyond age and clinical risk factors. This work sets the stage for further development of an aortic stenosis PRS.

Lipoprotein(a), C-Reactive Protein, and Cardiovascular Risk in Primary and Secondary Prevention Populations

Importance

Elevated lipoprotein(a) (Lp[a]) is a putative causal risk factor for atherosclerotic cardiovascular disease (ASCVD). There are conflicting data as to whether Lp(a) may increase cardiovascular risk only in the presence of concomitant inflammation.

Objective

To investigate whether Lp(a) is associated with cardiovascular risk independent of high-sensitivity C-reactive protein (hs-CRP) in both primary and secondary prevention populations.

Design, Setting, and Participants

This cohort study uses data from 3 distinct cohorts, 1 population-based cohort and 2 randomized clinical trials. Participants included individuals from the UK Biobank (data from 2006-2010) without prevalent ASCVD, participants in the FOURIER (TIMI 59) trial (data from 2013-2017) who had baseline Lp(a) and hs-CRP data, and participants in the SAVOR-TIMI 53 trial (data from 2010-2013) who had prevalent ASCVD and baseline values for Lp(a) and hs-CRP. The data analysis took place from November 2022 to November 2023.

Exposure

Baseline plasma Lp(a), considered either as a continuous variable or dichotomized at 125 nmol/L.

Main Outcomes and Measures

Risk of major adverse cardiovascular events (MACE) (composite of cardiovascular death, myocardial infarction [MI], or ischemic stroke), the individual MACE components, and peripheral artery disease (PAD).

Results

Among 357 220 individuals in the UK Biobank without prevalent ASCVD, 232 699 (65%) had low hs-CRP (<2 mg/L), and 124 521 (35%) had high hs-CRP (≥2 mg/L) values. In a Cox proportional hazard model adjusted for ASCVD risk factors, higher Lp(a) was associated with increased cardiovascular risk regardless of baseline hs-CRP value for MACE (hs-CRP ≥2 mg/L: hazard ratio [HR] per 50-nmol/L higher Lp[a], 1.05; 95% CI, 1.04-1.07; P < .001; for hs-CRP <2 mg/L: HR, 1.05; 95% CI, 1.04-1.07; P < .001; P = .80 for interaction), as well as MI, ischemic stroke, and PAD individually. Among 34 020 individuals in the FOURIER and SAVOR trials with baseline cardiometabolic disease, there were 17 643 (52%) with low and 16 377 (48%) with high baseline hs-CRP values. In Cox proportional hazard models using aggregated data from FOURIER and SAVOR, higher baseline Lp(a) was associated with increased cardiovascular risk regardless of baseline hs-CRP for MACE (hs-CRP ≥2 mg/L: HR per 50-nmol/L higher Lp[a], 1.02; 95% CI, 1.00-1.05; P = .04; hs-CRP <2 mg/L: HR, 1.05; 95% CI, 1.02-1.08; P < .001; P = .16 for interaction), MI, and PAD.

Conclusions and Relevance

In this study, higher levels of Lp(a) were associated with MACE, MI, and PAD in both primary and secondary prevention populations regardless of baseline hs-CRP value.

Frequency of Variants in Mendelian Alzheimer's Disease Genes within the Alzheimer's Disease Sequencing Project (ADSP)

Background

Prior studies using the ADSP data examined variants within presenilin-2 ( PSEN2 ), presenilin-1 ( PSEN1 ), and amyloid precursor protein ( APP ) genes. However, previously-reported clinically-relevant variants and other predicted damaging missense (DM) variants have not been characterized in a newer release of the Alzheimer's Disease Sequencing Project (ADSP).

Objective

To characterize previously-reported clinically-relevant variants and DM variants in PSEN2, PSEN1, APP within the participants from the ADSP.

Methods

We identified rare variants (MAF <1%) previously-reported in PSEN2 , PSEN1, and APP in the available ADSP sample of 14,641 individuals with whole genome sequencing and 16,849 individuals with whole exome sequencing available for research-use (N total = 31,490). We additionally curated variants in these three genes from ClinVar, OMIM, and Alzforum and report carriers of variants in clinical databases as well as predicted DM variants in these genes.

Results

We detected 31 previously-reported clinically-relevant variants with alternate alleles observed within the ADSP: 4 variants in PSEN2 , 25 in PSEN1 , and 2 in APP . The overall variant carrier rate for the 31 clinically-relevant variants in the ADSP was 0.3%. We observed that 79.5% of the variant carriers were cases compared to 3.9% were controls. In those with AD, the mean age of onset of AD among carriers of these clinically-relevant variants was 19.6 ± 1.4 years earlier compared with non-carriers (p-value=7.8×10 -57 ).

Conclusion

A small proportion of individuals in the ADSP are carriers of a previously-reported clinically-relevant variant allele for AD and these participants have significantly earlier age of AD onset compared to non-carriers.

Key variants via the Alzheimer's Disease Sequencing Project whole genome sequence data

INTRODUCTION

Genome-wide association studies (GWAS) have identified loci associated with Alzheimer's disease (AD) but did not identify specific causal genes or variants within those loci. Analysis of whole genome sequence (WGS) data, which interrogates the entire genome and captures rare variations, may identify causal variants within GWAS loci.

METHODS

We performed single common variant association analysis and rare variant aggregate analyses in the pooled population (N cases = 2184, N controls = 2383) and targeted analyses in subpopulations using WGS data from the Alzheimer's Disease Sequencing Project (ADSP). The analyses were restricted to variants within 100 kb of 83 previously identified GWAS lead variants.

RESULTS

Seventeen variants were significantly associated with AD within five genomic regions implicating the genes OARD1/NFYA/TREML1, JAZF1, FERMT2, and SLC24A4. KAT8 was implicated by both single variant and rare variant aggregate analyses.

DISCUSSION

This study demonstrates the utility of leveraging WGS to gain insights into AD loci identified via GWAS.

Lipid levels and risk of acute pancreatitis using bidirectional Mendelian randomization

Previous studies found lipid levels, especially triglycerides (TG), are associated with acute pancreatitis, but their causalities and bi-directions were not fully examined. We determined whether abnormal levels of TG, high-density lipoprotein cholesterol (HDL-C), and low-density lipoprotein cholesterol (LDL-C) are precursors and/or consequences of acute pancreatitis using bidirectional two-sample Mendelian randomization (MR) with two non-overlapping genome-wide association study (GWAS) summary statistics for lipid levels and acute pancreatitis. We found phenotypic associations that both higher TG levels and lower HDL-C levels contributed to increased risk of acute pancreatitis. Our GWAS meta-analysis of acute pancreatitis identified seven independent signals. Genetically predicted TG was positively associated with acute pancreatitis when using the variants specifically associated with TG using univariable MR [Odds ratio (OR), 95% CI 2.02, 1.22-3.31], but the reversed direction from acute pancreatitis to TG was not observed (mean difference = 0.003, SE = 0.002, P-value = 0.138). However, a bidirectional relationship of HDL-C and acute pancreatitis was observed: A 1-SD increment of genetically predicted HDL-C was associated with lower risk of acute pancreatitis (OR, 95% CI 0.84, 0.76-0.92) and genetically predisposed individuals with acute pancreatitis have, on average, 0.005 SD lower HDL-C (mean difference = - 0.005, SE = 0.002, P-value = 0.004). Our MR analysis confirms the evidence of TG as a risk factor of acute pancreatitis but not a consequence. A potential bidirectional relationship of HDL-C and acute pancreatitis occurs and raises the prospect of HDL-C modulation in the acute pancreatitis prevention and treatment.

Familial Hypercholesterolemia Variant and Cardiovascular Risk in Individuals With Elevated Cholesterol

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.

Protein identification for stroke progression via Mendelian Randomization in Million Veteran Program and UK Biobank

Background

Individuals who have experienced a stroke, or transient ischemic attack, face a heightened risk of future cardiovascular events. Identification of genetic and molecular risk factors for subsequent cardiovascular outcomes may identify effective therapeutic targets to improve prognosis after an incident stroke.

Methods

We performed genome-wide association studies (GWAS) for subsequent major adverse cardiovascular events (MACE) (Ncases=51,929, Ncntrl=39,980) and subsequent arterial ischemic stroke (AIS) Ncases=45,120, Ncntrl=46,789) after first incident stroke within the Million Veteran Program and UK Biobank. We then used genetic variants associated with proteins (pQTLs) to determine the effect of 1,463 plasma protein abundances on subsequent MACE using Mendelian randomization (MR).

Results

Two variants were significantly associated with subsequent cardiovascular events: rs76472767 (OR=0.75, 95% CI = 0.64-0.85, p= 3.69×10-08) with subsequent AIS and rs13294166 (OR=1.52, 95% CI = 1.37-1.67, p=3.77×10-08) with subsequent MACE. Using MR, we identified 2 proteins with an effect on subsequent MACE after a stroke: CCL27 (effect OR= 0.77, 95% CI = 0.66-0.88, adj. p=0.05), and TNFRSF14 (effect OR=1.42, 95% CI = 1.24-1.60, adj. p=0.006). These proteins are not associated with incident AIS and are implicated to have a role in inflammation.

Conclusions

We found evidence that two proteins with little effect on incident stroke appear to influence subsequent MACE after incident AIS. These associations suggest that inflammation is a contributing factor to subsequent MACE outcomes after incident AIS and highlights potential novel targets.

A statistical framework for powerful multi-trait rare variant analysis in large-scale whole-genome sequencing studies

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 NIPSNAP3A and an intergenic region on chromosome 1.

Genetic and clinical factors underlying a self-reported family history of heart disease

AIMS

To estimate how much information conveyed by self-reported family history of heart disease (FHHD) is already explained by clinical and genetic risk factors.

METHODS AND RESULTS

Cross-sectional analysis of UK Biobank participants without pre-existing coronary artery disease using a multivariable model with self-reported FHHD as the outcome. Clinical (diabetes, hypertension, smoking, apolipoprotein B-to-apolipoprotein AI ratio, waist-to-hip ratio, high sensitivity C-reactive protein, lipoprotein(a), triglycerides) and genetic risk factors (polygenic risk score for coronary artery disease [PRSCAD], heterozygous familial hypercholesterolemia [HeFH]) were exposures. Models were adjusted for age, sex, and cholesterol-lowering medication use. Multiple logistic regression models were fitted to associate FHHD with risk factors, with continuous variables treated as quintiles. Population attributable risks (PAR) were subsequently calculated from the resultant odds ratios. Among 166 714 individuals, 72 052 (43.2%) participants reported an FHHD. In a multivariable model, genetic risk factors PRSCAD (OR 1.30, CI 1.27-1.33) and HeFH (OR 1.31, 1.11-1.54) were most strongly associated with FHHD. Clinical risk factors followed: hypertension (OR 1.18, CI 1.15-1.21), lipoprotein(a) (OR 1.17, CI 1.14-1.20), apolipoprotein B-to-apolipoprotein AI ratio (OR 1.13, 95% CI 1.10-1.16), and triglycerides (OR 1.07, CI 1.04-1.10). For the PAR analyses: 21.9% (CI 18.19-25.63) of the risk of reporting an FHHD is attributed to clinical factors, 22.2% (CI% 20.44-23.88) is attributed to genetic factors, and 36.0% (CI 33.31-38.68) is attributed to genetic and clinical factors combined.

CONCLUSIONS

A combined model of clinical and genetic risk factors explains only 36% of the likelihood of FHHD, implying additional value in the family history.

Bayesian multivariate genetic analysis improves translational insights

While lipid traits are known essential mediators of cardiovascular disease, few approaches have taken advantage of their shared genetic effects. We apply a Bayesian multivariate size estimator, mash, to GWAS of four lipid traits in the Million Veterans Program (MVP) and provide posterior mean and local false sign rates for all effects. These estimates borrow information across traits to improve effect size accuracy. We show that controlling local false sign rates accurately and powerfully identifies replicable genetic associations and that multivariate control furthers the ability to explain complex diseases. Our application yields high concordance between independent datasets, more accurately prioritizes causal genes, and significantly improves polygenic prediction beyond state-of-the-art methods by up to 59% for lipid traits. The use of Bayesian multivariate genetic shrinkage has yet to be applied to human quantitative trait GWAS results, and we present a staged approach to prediction on a polygenic scale.

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

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

Genetic modification of inflammation- and clonal hematopoiesis-associated cardiovascular risk

Clonal hematopoiesis of indeterminate potential (CHIP) is associated with an increased risk of cardiovascular diseases (CVDs), putatively via inflammasome activation. We pursued an inflammatory gene modifier scan for CHIP-associated CVD risk among 424,651 UK Biobank participants. We identified CHIP using whole-exome sequencing data of blood DNA and modeled as a composite, considering all driver genes together, as well as separately for common drivers (DNMT3A, TET2, ASXL1, and JAK2). We developed predicted gene expression scores for 26 inflammasome-related genes and assessed how they modify CHIP-associated CVD risk. We identified IL1RAP as a potential key molecule for CHIP-associated CVD risk across genes and increased AIM2 gene expression leading to heightened JAK2- and ASXL1-associated CVD risk. We show that CRISPR-induced Asxl1-mutated murine macrophages had a particularly heightened inflammatory response to AIM2 agonism, associated with an increased DNA damage response, as well as increased IL-10 secretion, mirroring a CVD-protective effect of IL10 expression in ASXL1 CHIP. Our study supports the role of inflammasomes in CHIP-associated CVD and provides evidence to support gene-specific strategies to address CHIP-associated CVD risk.

Association of Common and Rare Variants with Alzheimer's Disease in over 13,000 Diverse Individuals with Whole-Genome Sequencing from the Alzheimer's Disease Sequencing Project

Alzheimer's Disease (AD) is a common disorder of the elderly that is both highly heritable and genetically heterogeneous. Here, we investigated the association between AD and both common variants and aggregates of rare coding and noncoding variants in 13,371 individuals of diverse ancestry with whole genome sequence (WGS) data. Pooled-population analyses identified genetic variants in or near APOE, BIN1, and LINC00320 significantly associated with AD (p < 5×10-8). Population-specific analyses identified a haplotype on chromosome 14 including PSEN1 associated with AD in Hispanics, further supported by aggregate testing of rare coding and noncoding variants in this region. Finally, we observed suggestive associations (p < 5×10-5) of aggregates of rare coding rare variants in ABCA7 among non-Hispanic Whites (p=5.4×10-6), and rare noncoding variants in the promoter of TOMM40 distinct of APOE in pooled-population analyses (p=7.2×10-8). Complementary pooled-population and population-specific analyses offered unique insights into the genetic architecture of AD.

Key variants via Alzheimer's Disease Sequencing Project whole genome sequence data

INTRODUCTION

Genome-wide association studies (GWAS) have identified loci associated with Alzheimer's disease (AD) but did not identify specific causal genes or variants within those loci. Analysis of whole genome sequence (WGS) data, which interrogates the entire genome and captures rare variations, may identify causal variants within GWAS loci.

METHODS

We performed single common variant association analysis and rare variant aggregate analyses in the pooled population (N cases=2,184, N controls=2,383) and targeted analyses in sub-populations using WGS data from the Alzheimer's Disease Sequencing Project (ADSP). The analyses were restricted to variants within 100 kb of 83 previously identified GWAS lead variants.

RESULTS

Seventeen variants were significantly associated with AD within five genomic regions implicating the genes OARD1/NFYA/TREML1, JAZF1, FERMT2, and SLC24A4. KAT8 was implicated by both single variant and rare variant aggregate analyses.

DISCUSSION

This study demonstrates the utility of leveraging WGS to gain insights into AD loci identified via GWAS.

Proteome-Wide Association Studies for Blood Lipids and Comparison with Transcriptome-Wide Association Studies

Blood lipid traits are treatable and heritable risk factors for heart disease, a leading cause of mortality worldwide. Although genome-wide association studies (GWAS) have discovered hundreds of variants associated with lipids in humans, most of the causal mechanisms of lipids remain unknown. To better understand the biological processes underlying lipid metabolism, we investigated the associations of plasma protein levels with total cholesterol (TC), triglycerides (TG), high-density lipoprotein cholesterol (HDL), and low-density lipoprotein cholesterol (LDL) in blood. We trained protein prediction models based on samples in the Multi-Ethnic Study of Atherosclerosis (MESA) and applied them to conduct proteome-wide association studies (PWAS) for lipids using the Global Lipids Genetics Consortium (GLGC) data. Of the 749 proteins tested, 42 were significantly associated with at least one lipid trait. Furthermore, we performed transcriptome-wide association studies (TWAS) for lipids using 9,714 gene expression prediction models trained on samples from peripheral blood mononuclear cells (PBMCs) in MESA and 49 tissues in the Genotype-Tissue Expression (GTEx) project. We found that although PWAS and TWAS can show different directions of associations in an individual gene, 40 out of 49 tissues showed a positive correlation between PWAS and TWAS signed p-values across all the genes, which suggests a high-level consistency between proteome-lipid associations and transcriptome-lipid associations.

Expression quantitative trait methylation analysis elucidates gene regulatory effects of DNA methylation: the Framingham Heart Study

Expression quantitative trait methylation (eQTM) analysis identifies DNA CpG sites at which methylation is associated with gene expression. The present study describes an eQTM resource of CpG-transcript pairs derived from whole blood DNA methylation and RNA sequencing gene expression data in 2115 Framingham Heart Study participants. We identified 70,047 significant cis CpG-transcript pairs at p < 1E-7 where the top most significant eGenes (i.e., gene transcripts associated with a CpG) were enriched in biological pathways related to cell signaling, and for 1208 clinical traits (enrichment false discovery rate [FDR] ≤ 0.05). We also identified 246,667 significant trans CpG-transcript pairs at p < 1E-14 where the top most significant eGenes were enriched in biological pathways related to activation of the immune response, and for 1191 clinical traits (enrichment FDR ≤ 0.05). Independent and external replication of the top 1000 significant cis and trans CpG-transcript pairs was completed in the Women's Health Initiative and Jackson Heart Study cohorts. Using significant cis CpG-transcript pairs, we identified significant mediation of the association between CpG sites and cardiometabolic traits through gene expression and identified shared genetic regulation between CpGs and transcripts associated with cardiometabolic traits. In conclusion, we developed a robust and powerful resource of whole blood eQTM CpG-transcript pairs that can help inform future functional studies that seek to understand the molecular basis of disease.

Novel Association of the NOTCH Pathway Regulator MIB1 Gene With the Development of Bicuspid Aortic Valve

Importance

Nonsyndromic bicuspid aortic valve (nsBAV) is the most common congenital heart valve malformation. BAV has a heritable component, yet only a few causative genes have been identified; understanding BAV genetics is a key point in developing personalized medicine.

Objective

To identify a new gene for nsBAV.

Design, Setting, and Participants

This was a comprehensive, multicenter, genetic association study based on candidate gene prioritization in a familial cohort followed by rare and common association studies in replication cohorts. Further validation was done using in vivo mice models. Study data were analyzed from October 2019 to October 2022. Three cohorts of patients with BAV were included in the study: (1) the discovery cohort was a large cohort of inherited cases from 29 pedigrees of French and Israeli origin; (2) the replication cohort 1 for rare variants included unrelated sporadic cases from various European ancestries; and (3) replication cohort 2 was a second validation cohort for common variants in unrelated sporadic cases from Europe and the US.

Main Outcomes and Measures

To identify a candidate gene for nsBAV through analysis of familial cases exome sequencing and gene prioritization tools. Replication cohort 1 was searched for rare and predicted deleterious variants and genetic association. Replication cohort 2 was used to investigate the association of common variants with BAV.

Results

A total of 938 patients with BAV were included in this study: 69 (7.4%) in the discovery cohort, 417 (44.5%) in replication cohort 1, and 452 (48.2%) in replication cohort 2. A novel human nsBAV gene, MINDBOMB1 homologue MIB1, was identified. MINDBOMB1 homologue (MIB1) is an E3-ubiquitin ligase essential for NOTCH-signal activation during heart development. In approximately 2% of nsBAV index cases from the discovery and replication 1 cohorts, rare MIB1 variants were detected, predicted to be damaging, and were significantly enriched compared with population-based controls (2% cases vs 0.9% controls; P = .03). In replication cohort 2, MIB1 risk haplotypes significantly associated with nsBAV were identified (permutation test, 1000 repeats; P = .02). Two genetically modified mice models carrying Mib1 variants identified in our cohort showed BAV on a NOTCH1-sensitized genetic background.

Conclusions and Relevance

This genetic association study identified the MIB1 gene as associated with nsBAV. This underscores the crucial role of the NOTCH pathway in the pathophysiology of BAV and its potential as a target for future diagnostic and therapeutic intervention.

Human gain-of-function variants in HNF1A confer protection from diabetes but independently increase hepatic secretion of atherogenic lipoproteins

Loss-of-function mutations in hepatocyte nuclear factor 1A (HNF1A) are known to cause rare forms of diabetes and alter hepatic physiology through unclear mechanisms. In the general population, 1:100 individuals carry a rare, protein-coding HNF1A variant, most of unknown functional consequence. To characterize the full allelic series, we performed deep mutational scanning of 11,970 protein-coding HNF1A variants in human hepatocytes and clinical correlation with 553,246 exome-sequenced individuals. Surprisingly, we found that ∼1:5 rare protein-coding HNF1A variants in the general population cause molecular gain of function (GOF), increasing the transcriptional activity of HNF1A by up to 50% and conferring protection from type 2 diabetes (odds ratio [OR] = 0.77, p = 0.007). Increased hepatic expression of HNF1A promoted a pro-atherogenic serum profile mediated in part by enhanced transcription of risk genes including ANGPTL3 and PCSK9. In summary, ∼1:300 individuals carry a GOF variant in HNF1A that protects carriers from diabetes but enhances hepatic secretion of atherogenic lipoproteins.

Genome-wide association analysis and Mendelian randomization proteomics identify drug targets for heart failure

We conduct a large-scale meta-analysis of heart failure genome-wide association studies (GWAS) consisting of over 90,000 heart failure cases and more than 1 million control individuals of European ancestry to uncover novel genetic determinants for heart failure. Using the GWAS results and blood protein quantitative loci, we perform Mendelian randomization and colocalization analyses on human proteins to provide putative causal evidence for the role of druggable proteins in the genesis of heart failure. We identify 39 genome-wide significant heart failure risk variants, of which 18 are previously unreported. Using a combination of Mendelian randomization proteomics and genetic cis-only colocalization analyses, we identify 10 additional putatively causal genes for heart failure. Findings from GWAS and Mendelian randomization-proteomics identify seven (CAMK2D, PRKD1, PRKD3, MAPK3, TNFSF12, APOC3 and NAE1) proteins as potential targets for interventions to be used in primary prevention of heart failure.

Rare variants in long non-coding RNAs are associated with blood lipid levels in the TOPMed Whole Genome Sequencing Study

Long non-coding RNAs (lncRNAs) are known to perform important regulatory functions. 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.

Integrative metabolomics differentiate coronary artery disease, peripheral artery disease, and venous thromboembolism risks

Rationale

Arterial and venous cardiovascular conditions, such as coronary artery disease (CAD), peripheral artery disease (PAD), and venous thromboembolism (VTE), are genetically correlated. Interrogating distinct and overlapping mechanisms may shed new light on disease mechanisms.

Objective

In this study, we aimed to: identify and compare (1) epidemiologic and (2) causal, genetic relationships between metabolites and CAD, PAD, and VTE.

Methods

We used metabolomic data from 95,402 individuals in the UK Biobank, excluding individuals with prevalent cardiovascular disease. Logistic regression models adjusted for age, sex, genotyping array, first five principal components of ancestry, and statin use estimated the epidemiologic associations of 249 metabolites with incident CAD, PAD, or VTE. Bidirectional two-sample Mendelian randomization (MR) estimated the causal effects between metabolites and cardiovascular phenotypes using genome-wide association summary statistics for metabolites (N = 118,466 from UK Biobank), CAD (N = 184,305 from CARDIoGRAMplusC4D 2015), PAD (N = 243,060 from Million Veterans Project) and VTE (N = 650,119 from Million Veterans Project). Multivariable MR (MVMR) was performed in subsequent analyses.

Results

We found that 194, 111, and 69 metabolites were epidemiologically associated (P < 0.001) with CAD, PAD, and VTE, respectively. Metabolomic profiles exhibited variable similarity between disease pairs: CAD and PAD (N = 100 shared associations, R2 = 0.499), CAD and VTE (N = 68, R2 = 0.455), and PAD and VTE (N = 54, R2 = 0.752). MR revealed 28 metabolites that increased risk for both CAD and PAD and 2 metabolites that increased risk for CAD but decreased risk for VTE. Despite strong epidemiologic overlap, no metabolites had a shared genetic relationship between PAD and VTE. MVMR revealed several metabolites with shared causal effects on CAD and PAD related to cholesterol content within very-low-density lipoprotein particles.

Conclusions

While common arterial and venous conditions are associated with overlapping metabolomic profiles, MR prioritized the role of remnant cholesterol in arterial diseases but not venous thrombosis.

Association of Rare Protein-Truncating DNA Variants in APOB or PCSK9 With Low-density Lipoprotein Cholesterol Level and Risk of Coronary Heart Disease

Importance

Protein-truncating variants (PTVs) in apolipoprotein B (APOB) and proprotein convertase subtilisin/kexin type 9 (PCSK9) are associated with significantly lower low-density lipoprotein (LDL) cholesterol concentrations. The association of these PTVs with coronary heart disease (CHD) warrants further characterization in large, multiracial prospective cohort studies.

Objective

To evaluate the association of PTVs in APOB and PCSK9 with LDL cholesterol concentrations and CHD risk.

Design, Setting, and Participants

This studied included participants from 5 National Heart, Lung, and Blood Institute (NHLBI) studies and the UK Biobank. NHLBI study participants aged 5 to 84 years were recruited between 1971 and 2002 across the US and underwent whole-genome sequencing. UK Biobank participants aged 40 to 69 years were recruited between 2006 and 2010 in the UK and underwent whole-exome sequencing. Data were analyzed from June 2021 to October 2022.

Exposures

PTVs in APOB and PCSK9.

Main Outcomes and Measures

Estimated untreated LDL cholesterol levels and CHD.

Results

Among 19 073 NHLBI participants (10 598 [55.6%] female; mean [SD] age, 52 [17] years), 139 (0.7%) carried an APOB or PCSK9 PTV, which was associated with 49 mg/dL (95% CI, 43-56) lower estimated untreated LDL cholesterol level. Over a median (IQR) follow-up of 21.5 (13.9-29.4) years, incident CHD was observed in 12 of 139 carriers (8.6%) vs 3029 of 18 934 noncarriers (16.0%), corresponding to an adjusted hazard ratio of 0.51 (95% CI, 0.28-0.89; P = .02). Among 190 464 UK Biobank participants (104 831 [55.0%] female; mean [SD] age, 57 [8] years), 662 (0.4%) carried a PTV, which was associated with 45 mg/dL (95% CI, 42-47) lower estimated untreated LDL cholesterol level. Estimated CHD risk by age 75 years was 3.7% (95% CI, 2.0-5.3) in carriers vs 7.0% (95% CI, 6.9-7.2) in noncarriers, corresponding to an adjusted hazard ratio of 0.51 (95% CI, 0.32-0.81; P = .004).

Conclusions and Relevance

Among 209 537 individuals in this study, 0.4% carried an APOB or PCSK9 PTV that was associated with less exposure to LDL cholesterol and a 49% lower risk of CHD.

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

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

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Exploiting family history in aggregation unit-based genetic association tests

The development of sequencing technology calls for new powerful methods to detect disease associations and lower the cost of sequencing studies. Family history (FH) contains information on disease status of relatives, adding valuable information about the probands' health problems and risk of diseases. Incorporating data from FH is a cost-effective way to improve statistical evidence in genetic studies, and moreover, overcomes limitations in study designs with insufficient cases or missing genotype information for association analysis. We proposed family history aggregation unit-based test (FHAT) and optimal FHAT (FHAT-O) to exploit available FH for rare variant association analysis. Moreover, we extended liability threshold model of case-control status and FH (LT-FH) method in aggregated unit-based methods and compared that with FHAT and FHAT-O. The computational efficiency and flexibility of the FHAT and FHAT-O were demonstrated through both simulations and applications. We showed that FHAT, FHAT-O, and LT-FH methods offer reasonable control of the type I error unless case/control ratio is unbalanced, in which case they result in smaller inflation than that observed with conventional methods excluding FH. We also demonstrated that FHAT and FHAT-O are more powerful than LT-FH and conventional methods in many scenarios. By applying FHAT and FHAT-O to the analysis of all cause dementia and hypertension using the exome sequencing data from the UK Biobank, we showed that our methods can improve significance for known regions. Furthermore, we replicated the previous associations in all cause dementia and hypertension and detected novel regions through the exome-wide analysis.

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

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

Bayesian Genetic Colocalization Test of Two Traits Using coloc

Genetic colocalization is an approach for determining whether a genetic variant at a particular locus is shared across multiple phenotypes. Genome-wide association studies (GWAS) have successfully mapped genetic variants associated with thousands of complex traits and diseases. However, a large proportion of GWAS signals fall in non-coding regions of the genome, making functional interpretation a challenge. Colocalization relies on a Bayesian framework that can integrate summary statistics, for example those derived from GWAS and expression quantitative trait loci (eQTL) mapping, to assess whether two or more independent association signals at a region of interest are consistent with a shared causal variant. The results from a colocalization analysis may be used to evaluate putative causal relationships between omics-based molecular measurements and a complex disease, and can generate hypotheses that may be followed up by tailored experiments. In this article, we present an easy and straightforward protocol for conducting a Bayesian test for colocalization of two traits using the 'coloc' package in R with summary-level results derived from GWAS and eQTL studies. We also provide general guidelines that can assist in the interpretation of findings generated from colocalization analyses. © 2022 Wiley Periodicals LLC. Basic Protocol: Performing a genetic colocalization analysis using the 'coloc' package in R and summary-level data Support Protocol: Installing the 'coloc' R package.

Genome-wide discovery for diabetes-dependent triglycerides-associated loci

PURPOSE

We aimed to discover loci associated with triglyceride (TG) levels in the context of type 2 diabetes (T2D). We conducted a genome-wide association study (GWAS) in 424,120 genotyped participants of the UK Biobank (UKB) with T2D status and TG levels.

METHODS

We stratified the cohort based on T2D status and conducted association analyses of TG levels for genetic variants with minor allele count (MAC) at least 20 in each stratum. Effect differences of genetic variants by T2D status were determined by Cochran's Q-test and we validated the significantly associated variants in the Mass General Brigham Biobank (MGBB).

RESULTS

Among 21,176 T2D and 402,944 non-T2D samples from UKB, stratified GWAS identified 19 and 315 genomic risk loci significantly associated with TG levels, respectively. Only chr6p21.32 exhibited genome-wide significant heterogeneity (I2 = 98.4%; pheterogeneity = 2.1x10-15), with log(TG) effect estimates of -0.066 (95%CI: -0.082, -0.050) and 0.002 (95%CI: -0.002, 0.006) for T2D and non-T2D, respectively. The lead variant rs9274619:A (allele frequency 0.095) is located 2Kb upstream of the HLA-DQB1 gene, between HLA-DQB1 and HLA-DQA2 genes. We replicated this finding among 25,137 participants (6,951 T2D cases) of MGBB (pheterogeneity = 9.5x10-3). Phenome-wide interaction association analyses showed that the lead variant was strongly associated with a concomitant diagnosis of type 1 diabetes (T1D) as well as diabetes-associated complications.

CONCLUSION

In conclusion, we identified an intergenic variant near HLA-DQB1/DQA2 significantly associates with decreased triglycerides only among those with T2D and highlights an immune overlap with T1D.

Whole genome sequence analysis of blood lipid levels in >66,000 individuals

Blood lipids are heritable modifiable causal factors for coronary artery disease. Despite well-described monogenic and polygenic bases of dyslipidemia, limitations remain in discovery of lipid-associated alleles using whole genome sequencing (WGS), partly due to limited sample sizes, ancestral diversity, and interpretation of clinical significance. Among 66,329 ancestrally diverse (56% non-European) participants, we associate 428M variants from deep-coverage WGS with lipid levels; ~400M variants were not assessed in prior lipids genetic analyses. We find multiple lipid-related genes strongly associated with blood lipids through analysis of common and rare coding variants. We discover several associated rare non-coding variants, largely at Mendelian lipid genes. Notably, we observe rare LDLR intronic variants associated with markedly increased LDL-C, similar to rare LDLR exonic variants. In conclusion, we conducted a systematic whole genome scan for blood lipids expanding the alleles linked to lipids for multiple ancestries and characterize a clinically-relevant rare non-coding variant model for lipids.

Non-linear machine learning models incorporating SNPs and PRS improve polygenic prediction in diverse human populations

Polygenic risk scores (PRS) are commonly used to quantify the inherited susceptibility for a trait, yet they fail to account for non-linear and interaction effects between single nucleotide polymorphisms (SNPs). We address this via a machine learning approach, validated in nine complex phenotypes in a multi-ancestry population. We use an ensemble method of SNP selection followed by gradient boosted trees (XGBoost) to allow for non-linearities and interaction effects. We compare our results to the standard, linear PRS model developed using PRSice, LDpred2, and lassosum2. Combining a PRS as a feature in an XGBoost model results in a relative increase in the percentage variance explained compared to the standard linear PRS model by 22% for height, 27% for HDL cholesterol, 43% for body mass index, 50% for sleep duration, 58% for systolic blood pressure, 64% for total cholesterol, 66% for triglycerides, 77% for LDL cholesterol, and 100% for diastolic blood pressure. Multi-ancestry trained models perform similarly to specific racial/ethnic group trained models and are consistently superior to the standard linear PRS models. This work demonstrates an effective method to account for non-linearities and interaction effects in genetics-based prediction models.

Combining a standard polygenic risk score (PRS) as a feature in a machine learning model increases the percentage variance explained for those traits, helping to account for non-linearities or interaction effects in genetics-based prediction models.

Sugar-Sweetened Beverage Consumption and Plasma Lipoprotein Cholesterol, Apolipoprotein, and Lipoprotein Particle Size Concentrations in US Adults

BACKGROUND

Prospective cohort studies have found a relation between sugar-sweetened beverage (SSB) consumption (sodas and fruit drinks) and dyslipidemia. There is limited evidence linking SSB consumption to emerging features of dyslipidemia, which can be characterized by variation in lipoprotein particle size, remnant-like particle (RLP), and apolipoprotein concentrations.

OBJECTIVES

To examine the association between SSB consumption and plasma lipoprotein cholesterol, apolipoprotein, and lipoprotein particle size concentrations among US adults.

METHODS

We examined participants from the Framingham Offspring Study (FOS; 1987-1995, n = 3047) and the Women's Health Study (1992, n = 26,218). Concentrations of plasma LDL cholesterol, apolipoprotein B (apoB), HDL cholesterol, apolipoprotein A1 (apoA1), triglyceride (TG), and non-HDL cholesterol, as well as total cholesterol:HDL cholesterol ratio and apoB:apoA1 ratio, were quantified in both cohorts; concentrations of apolipoprotein E, apolipoprotein C3, RLP-TG, and RLP cholesterol (RLP-C) were measured in the FOS only. Lipoprotein particle sizes were calculated from nuclear magnetic resonance signals for lipoprotein particle subclass concentrations (TG-rich lipoprotein particles [TRL-Ps]: very large, large, medium, small, and very small; LDL particles [LDL-Ps]: large, medium, and small; HDL particles [HDL-Ps]: large, medium, and small). SSB consumption was estimated from food frequency questionnaire data. We examined the associations between SSB consumption and all lipoprotein and apoprotein measures in linear regression models, adjusting for confounding factors such as lifestyle, diet, and traditional lipoprotein risk factors.

RESULTS

SSB consumption was positively associated with LDL cholesterol, apoB, TG, RLP-TG, RLP-C, and non-HDL cholesterol concentrations and total cholesterol:HDL cholesterol and apoB:apoA1 ratios; and negatively associated with HDL cholesterol and apoA1 concentrations (P-trend range: <0.0001 to 0.008). After adjustment for traditional lipoprotein risk factors, SSB consumers had smaller LDL-P and HDL-P sizes; lower concentrations of large LDL-Ps and medium HDL-Ps; and higher concentrations of small LDL-Ps, small HDL-Ps, and large TRL-Ps (P-trend range: <0.0001 to 0.001).

CONCLUSIONS

Higher SSB consumption was associated with multiple emerging features of dyslipidemia that have been linked to higher cardiometabolic risk in US adults.

Association of Carotid Intima Media Thickening with Future Brain Region Specific Amyloid-β Burden

BACKGROUND

Carotid atherosclerosis is associated with cognitive impairment and dementia, though there is limited evidence of a direct link between carotid disease and amyloid-β (Aβ) burden.

OBJECTIVE

We studied the association of baseline and progressive carotid intima media thickness (CIMT) with Aβ on 11C-Pittsburgh Compound B (PiB) to determine if those with carotid atherosclerosis would have higher Aβ burden.

METHODS

We studied 47 participants from the Framingham Offspring cohort with carotid ultrasounds measuring CIMT at their 6th clinic examination (aged 49.5±5.7 years) and an average of 9.6 years later, and PiB imaging measuring Aβ on average 22.1 years post baseline. We used multivariate linear regression analyses to relate baseline, follow-up, mean, and progression of internal carotid artery (ICA) and common carotid artery (CCA) CIMT to Aβ in brain regions associated with Alzheimer's disease (AD) and related dementias (ADRD), adjusting for age, sex, and other vascular risk factors.

RESULTS

Participants with higher mean ICA IMT had more Aβ in the precuneus (beta±standard error [β±SE]: 0.466±0.171 mm, p = 0.01) and the frontal, lateral, and retrosplenial regions (β±SE: 0.392±0.164 mm, p = 0.022) after adjusting for age, sex, vascular risk factors, and medication use. We did not find an association between any CCA IMT measures and Aβ or progression of ICA or CCA IMT and Aβ.

CONCLUSION

Carotid atherosclerosis, as measured by ICA IMT, is associated with increased Aβ burden later in life. These findings support a link between vascular disease and AD/ADRD pathophysiology.

Whole Genome DNA and RNA Sequencing of Whole Blood Elucidates the Genetic Architecture of Gene Expression Underlying a Wide Range of Diseases

To create a scientific resource of expression quantitative trail loci (eQTL), we conducted a genome-wide association study (GWAS) using genotypes obtained from whole genome sequencing (WGS) of DNA and gene expression levels from RNA sequencing (RNA-seq) of whole blood in 2622 participants in Framingham Heart Study. We identified 6,778,286 cis -eQTL variant-gene transcript (eGene) pairs at p  < 5x10 - 8 (2,855,111 unique cis -eQTL variants and 15,982 unique eGenes) and 1,469,754 trans -eQTL variant-eGene pairs at p  < 1e-12 (526,056 unique trans -eQTL variants and 7,233 unique eGenes). In addition, 442,379 cis -eQTL variants were associated with expression of 1518 long non-protein coding RNAs (lncRNAs). Gene Ontology (GO) analyses revealed that the top GO terms for cis- eGenes are enriched for immune functions (FDR < 0.05). The cis -eQTL variants are enriched for SNPs reported to be associated with 815 traits in prior GWAS, including cardiovascular disease risk factors. As proof of concept, we used this eQTL resource in conjunction with genetic variants from public GWAS databases in causal inference testing (e.g., COVID-19 severity). After Bonferroni correction, Mendelian randomization analyses identified putative causal associations of 60 eGenes with systolic blood pressure, 13 genes with coronary artery disease, and seven genes with COVID-19 severity. This study created a comprehensive eQTL resource via BioData Catalyst that will be made available to the scientific community. This will advance understanding of the genetic architecture of gene expression underlying a wide range of diseases.

Genome-wide and phenome-wide analysis of ideal cardiovascular health in the VA Million Veteran Program

BACKGROUND

Genetic studies may help identify causal pathways; therefore, we sought to identify genetic determinants of ideal CVH and their association with CVD outcomes in the multi-population Veteran Administration Million Veteran Program.

METHODS

An ideal health score (IHS) was calculated from 3 clinical factors (blood pressure, total cholesterol, and blood glucose levels) and 3 behavioral factors (smoking status, physical activity, and BMI), ascertained at baseline. Multi-population genome-wide association study (GWAS) was performed on IHS and binary ideal health using linear and logistic regression, respectively. Using the genome-wide significant SNPs from the IHS GWAS, we created a weighted IHS polygenic risk score (PRSIHS) which was used (i) to conduct a phenome-wide association study (PheWAS) of associations between PRSIHS and ICD-9 phenotypes and (ii) to further test for associations with mortality and selected CVD outcomes using logistic and Cox regression and, as an instrumental variable, in Mendelian Randomization.

RESULTS

The discovery and replication cohorts consisted of 142,404 (119,129 European American (EUR); 16,495 African American (AFR)), and 45,766 (37,646 EUR; 5,366 AFR) participants, respectively. The mean age was 65.8 years (SD = 11.2) and 92.7% were male. Overall, 4.2% exhibited ideal CVH based on the clinical and behavioral factors. In the multi-population meta-analysis, variants at 17 loci were associated with IHS and each had known GWAS associations with multiple components of the IHS. PheWAS analysis in 456,026 participants showed that increased PRSIHS was associated with a lower odds ratio for many CVD outcomes and risk factors. Both IHS and PRSIHS measures of ideal CVH were associated with significantly less CVD outcomes and CVD mortality.

CONCLUSION

A set of high interest genetic variants contribute to the presence of ideal CVH in a multi-ethnic cohort of US Veterans. Genetically influenced ideal CVH is associated with lower odds of CVD outcomes and mortality.

Whole Genome DNA and RNA Sequencing of Whole Blood Elucidates the Genetic Architecture of Gene Expression Underlying a Wide Range of Diseases

To create a scientific resource of expression quantitative trail loci (eQTL), we conducted a genome-wide association study (GWAS) using genotypes obtained from whole genome sequencing (WGS) of DNA and gene expression levels from RNA sequencing (RNA-seq) of whole blood in 2622 participants in Framingham Heart Study. We identified 6,778,286 cis -eQTL variant-gene transcript (eGene) pairs at p <5×10 -8 (2,855,111 unique cis -eQTL variants and 15,982 unique eGenes) and 1,469,754 trans -eQTL variant-eGene pairs at p <1e-12 (526,056 unique trans -eQTL variants and 7,233 unique eGenes). In addition, 442,379 cis -eQTL variants were associated with expression of 1518 long non-protein coding RNAs (lncRNAs). Gene Ontology (GO) analyses revealed that the top GO terms for cis- eGenes are enriched for immune functions (FDR <0.05). The cis -eQTL variants are enriched for SNPs reported to be associated with 815 traits in prior GWAS, including cardiovascular disease risk factors. As proof of concept, we used this eQTL resource in conjunction with genetic variants from public GWAS databases in causal inference testing (e.g., COVID-19 severity). After Bonferroni correction, Mendelian randomization analyses identified putative causal associations of 60 eGenes with systolic blood pressure, 13 genes with coronary artery disease, and seven genes with COVID-19 severity. This study created a comprehensive eQTL resource via BioData Catalyst that will be made available to the scientific community. This will advance understanding of the genetic architecture of gene expression underlying a wide range of diseases.

Mendelian randomization supports bidirectional causality between telomere length and clonal hematopoiesis of indeterminate potential

Human genetic studies support an inverse causal relationship between leukocyte telomere length (LTL) and coronary artery disease (CAD), but directionally mixed effects for LTL and diverse malignancies. Clonal hematopoiesis of indeterminate potential (CHIP), characterized by expansion of hematopoietic cells bearing leukemogenic mutations, predisposes both hematologic malignancy and CAD. TERT (which encodes telomerase reverse transcriptase) is the most significantly associated germline locus for CHIP in genome-wide association studies. Here, we investigated the relationship between CHIP, LTL, and CAD in the Trans-Omics for Precision Medicine (TOPMed) program (n = 63,302) and UK Biobank (n = 47,080). Bidirectional Mendelian randomization studies were consistent with longer genetically imputed LTL increasing propensity to develop CHIP, but CHIP then, in turn, hastens to shorten measured LTL (mLTL). We also demonstrated evidence of modest mediation between CHIP and CAD by mLTL. Our data promote an understanding of potential causal relationships across CHIP and LTL toward prevention of CAD.

Family history aggregation unit-based tests to detect rare genetic variant associations with application to the Framingham Heart Study

A challenge in standard genetic studies is maintaining good power to detect associations, especially for low prevalent diseases and rare variants. The traditional methods are most powerful when evaluating the association between variants in balanced study designs. Without accounting for family correlation and unbalanced case-control ratio, these analyses could result in inflated type I error. One cost-effective solution to increase statistical power is exploitation of available family history (FH) that contains valuable information about disease heritability. Here, we develop methods to address the aforementioned type I error issues while providing optimal power to analyze aggregates of rare variants by incorporating additional information from FH. With enhanced power in these methods exploiting FH and accounting for relatedness and unbalanced designs, we successfully detect genes with suggestive associations with Alzheimer disease, dementia, and type 2 diabetes by using the exome chip data from the Framingham Heart Study.

Genetic Loci Associated With COVID-19 Positivity and Hospitalization in White, Black, and Hispanic Veterans of the VA Million Veteran Program

SARS-CoV-2 has caused symptomatic COVID-19 and widespread death across the globe. We sought to determine genetic variants contributing to COVID-19 susceptibility and hospitalization in a large biobank linked to a national United States health system. We identified 19,168 (3.7%) lab-confirmed COVID-19 cases among Million Veteran Program participants between March 1, 2020, and February 2, 2021, including 11,778 Whites, 4,893 Blacks, and 2,497 Hispanics. A multi-population genome-wide association study (GWAS) for COVID-19 outcomes identified four independent genetic variants (rs8176719, rs73062389, rs60870724, and rs73910904) contributing to COVID-19 positivity, including one novel locus found exclusively among Hispanics. We replicated eight of nine previously reported genetic associations at an alpha of 0.05 in at least one population-specific or the multi-population meta-analysis for one of the four MVP COVID-19 outcomes. We used rs8176719 and three additional variants to accurately infer ABO blood types. We found that A, AB, and B blood types were associated with testing positive for COVID-19 compared with O blood type with the highest risk for the A blood group. We did not observe any genome-wide significant associations for COVID-19 severity outcomes among those testing positive. Our study replicates prior GWAS findings associated with testing positive for COVID-19 among mostly White samples and extends findings at three loci to Black and Hispanic individuals. We also report a new locus among Hispanics requiring further investigation. These findings may aid in the identification of novel therapeutic agents to decrease the morbidity and mortality of COVID-19 across all major ancestral populations.

Multi-ancestry fine mapping implicates OAS1 splicing in risk of severe COVID-19

The OAS1/2/3 cluster has been identified as a risk locus for severe COVID-19 among individuals of European ancestry, with a protective haplotype of approximately 75 kilobases (kb) derived from Neanderthals in the chromosomal region 12q24.13. This haplotype contains a splice variant of OAS1, which occurs in people of African ancestry independently of gene flow from Neanderthals. Using trans-ancestry fine-mapping approaches in 20,779 hospitalized cases, we demonstrate that this splice variant is likely to be the SNP responsible for the association at this locus, thus strongly implicating OAS1 as an effector gene influencing COVID-19 severity.

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

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

Two-Sample Multivariable Mendelian Randomization Analysis Using R

Mendelian randomization is a framework that uses measured variation in genes for assessing and estimating the causal effect of an exposure on an outcome. Multivariable Mendelian randomization is an extension that can assess the causal effect of multiple exposures on an outcome, and can be advantageous when considering a set (>1) of potentially correlated candidate risk factors in evaluating the causal effect of each on a health outcome, accounting for measured pleiotropy. This can be seen, for example, in determining the causal effects of lipids and cholesterol on type 2 diabetes risk, where the correlated risk factors share genetic predictors. Similar to univariate Mendelian randomization, multivariable Mendelian randomization can be conducted using two-sample summary-level data where the gene-exposure and gene-outcome associations are derived from separate samples from the same underlying population. Here, we present a protocol for conducting a two-sample multivariable Mendelian randomization study using the 'MVMR' package in R and summary-level genetic data. We also provide a protocol for searching and obtaining instruments using available data sources in the 'MRInstruments' R package. Finally, we provide general guidelines and discuss the utility of performing a multivariable Mendelian randomization analysis for simultaneously assessing causality of multiple exposures. © 2021 Wiley Periodicals LLC. Basic Protocol: Performing a two-sample multivariable Mendelian randomization analysis using the 'MVMR' package in R and summarized genetic data Support Protocol 1: Installing the 'MVMR' R package Support Protocol 2: Obtaining instruments from the 'MRInstruments' R package.

History of Suicide Attempts and COVID-19 Infection in Veterans with Schizophrenia or Schizoaffective Disorder: Moderating Effects of Age and Body Mass Index

Introduction

Relative to the general population, patients with schizophrenia or schizoaffective disorder have higher rates of suicide attempts and mortality from COVID-19 infection. Therefore, determining whether a history of suicide attempt is associated with COVID-19 in patients with schizophrenia or schizoaffective disorder has implications for COVID-19 vulnerability stratification in this patient population.

Methods

We carried out cross-sectional analyses of electronic health records of veterans with a diagnosis of schizophrenia or schizoaffective disorder that received treatment at any United States Veterans Affairs Medical Center between January 1, 2020, and January 31, 2021. We used logistic regression to estimate unadjusted and adjusted (including age, sex, race, marital status, body mass index (BMI), and a medical comorbidity score) odds ratios (ORs) for COVID-19 positivity in suicide attempters relative to nonattempters.

Results

A total of 101,032 veterans (mean age 56.67 ± 13.13 years; males 91,715 [90.8%]) were included in the analyses. There were 2,703 (2.7%) suicide attempters and 719 (0.7%) patients were positive for COVID-19. The association between history of suicide attempt and COVID-19 positivity was modified by age and BMI, such that the relationship was only significant in patients younger than 59 years, and in obese (BMI ≥30) patients (adjusted OR 3.42, 95% CI 2.02–5.79 and OR 2.85, 95% CI 1.65–4.94, respectively).

Conclusions

Higher rates of COVID-19 in young or obese suicide attempters with a diagnosis of schizophrenia or schizoaffective disorder might be due to the elevated risk for the infection in this subgroup of patients.

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

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

A System for Phenotype Harmonization in the National Heart, Lung, and Blood Institute Trans-Omics for Precision Medicine (TOPMed) Program

Genotype-phenotype association studies often combine phenotype data from multiple studies to increase statistical power. Harmonization of the data usually requires substantial effort due to heterogeneity in phenotype definitions, study design, data collection procedures, and data-set organization. Here we describe a centralized system for phenotype harmonization that includes input from phenotype domain and study experts, quality control, documentation, reproducible results, and data-sharing mechanisms. This system was developed for the National Heart, Lung, and Blood Institute's Trans-Omics for Precision Medicine (TOPMed) program, which is generating genomic and other -omics data for more than 80 studies with extensive phenotype data. To date, 63 phenotypes have been harmonized across thousands of participants (recruited in 1948-2012) from up to 17 studies per phenotype. Here we discuss challenges in this undertaking and how they were addressed. The harmonized phenotype data and associated documentation have been submitted to National Institutes of Health data repositories for controlled access by the scientific community. We also provide materials to facilitate future harmonization efforts by the community, which include 1) the software code used to generate the 63 harmonized phenotypes, enabling others to reproduce, modify, or extend these harmonizations to additional studies, and 2) the results of labeling thousands of phenotype variables with controlled vocabulary terms.

Sugar-Sweetened Beverage Consumption May Modify Associations Between Genetic Variants in the CHREBP (Carbohydrate Responsive Element Binding Protein) Locus and HDL-C (High-Density Lipoprotein Cholesterol) and Triglyceride Concentrations

Supplemental Digital Content is available in the text.

Background:

ChREBP (carbohydrate responsive element binding protein) is a transcription factor that responds to sugar consumption. Sugar-sweetened beverage (SSB) consumption and genetic variants in the CHREBP locus have separately been linked to HDL-C (high-density lipoprotein cholesterol) and triglyceride concentrations. We hypothesized that SSB consumption would modify the association between genetic variants in the CHREBP locus and dyslipidemia.

Methods:

Data from 11 cohorts from the Cohorts for Heart and Aging Research in Genomic Epidemiology consortium (N=63 599) and the UK Biobank (N=59 220) were used to quantify associations of SSB consumption, genetic variants, and their interaction on HDL-C and triglyceride concentrations using linear regression models. A total of 1606 single nucleotide polymorphisms within or near CHREBP were considered. SSB consumption was estimated from validated questionnaires, and participants were grouped by their estimated intake.

Results:

In a meta-analysis, rs71556729 was significantly associated with higher HDL-C concentrations only among the highest SSB consumers (β, 2.12 [95% CI, 1.16–3.07] mg/dL per allele; P<0.0001), but not significantly among the lowest SSB consumers (P=0.81; P_Diff <0.0001). Similar results were observed for 2 additional variants (rs35709627 and rs71556736). For triglyceride, rs55673514 was positively associated with triglyceride concentrations only among the highest SSB consumers (β, 0.06 [95% CI, 0.02–0.09] ln-mg/dL per allele, P=0.001) but not the lowest SSB consumers (P=0.84; P_Diff=0.0005).

Conclusions:

Our results identified genetic variants in the CHREBP locus that may protect against SSB-associated reductions in HDL-C and other variants that may exacerbate SSB-associated increases in triglyceride concentrations.

Registration:

URL: https://www.clinicaltrials.gov; Unique identifier: NCT00005133, NCT00005121, NCT00005487, and NCT00000479.