Whole genome sequencing based analysis of inflammation biomarkers in the Trans-Omics for Precision Medicine (TOPMed) consortium.

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 (with varying sample size by trait, where the minimum sample size was n = 737 for MMP-1). 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.

MagicalRsq-X: A cross-cohort transferable genotype imputation quality metric.

Since genotype imputation was introduced, researchers have been relying on the estimated imputation quality from imputation software to perform post-imputation quality control (QC). However, this quality estimate (denoted as Rsq) performs less well for lower-frequency variants. We recently published MagicalRsq, a machine-learning-based imputation quality calibration, which leverages additional typed markers from the same cohort and outperforms Rsq as a QC metric. In this work, we extended the original MagicalRsq to allow cross-cohort model training and named the new model MagicalRsq-X. We removed the cohort-specific estimated minor allele frequency and included linkage disequilibrium scores and recombination rates as additional features. Leveraging whole-genome sequencing data from TOPMed, specifically participants in the BioMe, JHS, WHI, and MESA studies, we performed comprehensive cross-cohort evaluations for predominantly European and African ancestral individuals based on their inferred global ancestry with the 1000 Genomes and Human Genome Diversity Project data as reference. Our results suggest MagicalRsq-X outperforms Rsq in almost every setting, with 7.3%-14.4% improvement in squared Pearson correlation with true R2, corresponding to 85-218 K variant gains. We further developed a metric to quantify the genetic distances of a target cohort relative to a reference cohort and showed that such metric largely explained the performance of MagicalRsq-X models. Finally, we found MagicalRsq-X saved up to 53 known genome-wide significant variants in one of the largest blood cell trait GWASs that would be missed using the original Rsq for QC. In conclusion, MagicalRsq-X shows superiority for post-imputation QC and benefits genetic studies by distinguishing well and poorly imputed lower-frequency variants.

Time-to-Event Genome-Wide Association Study for Incident Cardiovascular Disease in People With Type 2 Diabetes.

OBJECTIVE: To identify genetic risk factors for incident cardiovascular disease (CVD) among people with type 2 diabetes (T2D). RESEARCH DESIGN AND METHODS: We conducted a multiancestry time-to-event genome-wide association study for incident CVD among people with T2D. We also tested 204 known coronary artery disease (CAD) variants for association with incident CVD. RESULTS: Among 49,230 participants with T2D, 8,956 had incident CVD events (event rate 18.2%). We identified three novel genetic loci for incident CVD: rs147138607 (near CACNA1E/ZNF648, hazard ratio [HR] 1.23, P = 3.6 × 10-9), rs77142250 (near HS3ST1, HR 1.89, P = 9.9 × 10-9), and rs335407 (near TFB1M/NOX3, HR 1.25, P = 1.5 × 10-8). Among 204 known CAD loci, 5 were associated with incident CVD in T2D (multiple comparison-adjusted P < 0.00024, 0.05/204). A standardized polygenic score of these 204 variants was associated with incident CVD with HR 1.14 (P = 1.0 × 10-16). CONCLUSIONS: The data point to novel and known genomic regions associated with incident CVD among individuals with T2D.

Adjusting for principal components can induce spurious associations in genome-wide association studies in admixed populations.

Principal component analysis (PCA) is widely used to control for population structure in genome-wide association studies (GWAS). Top principal components (PCs) typically reflect population structure, but challenges arise in deciding how many PCs are needed and ensuring that PCs do not capture other artifacts such as regions with atypical linkage disequilibrium (LD). In response to the latter, many groups suggest performing LD pruning or excluding known high LD regions prior to PCA. However, these suggestions are not universally implemented and the implications for GWAS are not fully understood, especially in the context of admixed populations. In this paper, we investigate the impact of pre-processing and the number of PCs included in GWAS models in African American samples from the Women's Women's Health Initiative SNP Health Association Resource and two Trans-Omics for Precision Medicine Whole Genome Sequencing Project contributing studies (Jackson Heart Study and Genetic Epidemiology of Chronic Obstructive Pulmonary Disease Study). In all three samples, we find the first PC is highly correlated with genome-wide ancestry whereas later PCs often capture local genomic features. The pattern of which, and how many, genetic variants are highly correlated with individual PCs differs from what has been observed in prior studies focused on European populations and leads to distinct downstream consequences: adjusting for such PCs yields biased effect size estimates and elevated rates of spurious associations due to the phenomenon of collider bias. Excluding high LD regions identified in previous studies does not resolve these issues. LD pruning proves more effective, but the optimal choice of thresholds varies across datasets. Altogether, our work highlights unique issues that arise when using PCA to control for ancestral heterogeneity in admixed populations and demonstrates the importance of careful pre-processing and diagnostics to ensure that PCs capturing multiple local genomic features are not included in GWAS models.

Radon Exposure and Incident Stroke Risk in the Women's Health Initiative.

BACKGROUND AND OBJECTIVES: Little is known about the role of radon in the epidemiology of stroke among women. We therefore examined the association between home radon exposure and risk of stroke among middle-aged and older women in the United States. METHODS: We conducted a prospective cohort study of postmenopausal women aged 50-79 years at baseline (1993-1998) in the Women's Health Initiative. We measured exposures as 2-day, indoor, lowest living-level average radon concentrations in picocuries per liter (pCi/L) as estimated in 1993 by the US Geological Survey and reviewed by the Association of American State Geologists under the Indoor Radon Abatement Act. We used Cox proportional hazards models to estimate risk of incident, neurologist-adjudicated stroke during follow-up through 2020 as a hazard ratio and 95% CI, adjusting for study design and participant demographic, social, behavioral, and clinical characteristics. RESULTS: Among 158,910 women without stroke at baseline (mean age 63.2 years; 83% white), 6,979 incident strokes were identified over follow-up (mean 13.4 years). Incidence rates were 333, 343, and 349 strokes per 100,000 woman-years at radon concentrations of <2, 2-4, and >4 pCi/L, respectively. Compared with women living at concentrations <2 pCi/L, those at 2-4 and >4 pCi/L had higher covariate-adjusted risks of incident stroke: hazard ratio (95% CI) 1.06 (0.99-1.13) and 1.14 (1.05-1.22). Using nonlinear spline functions to model radon, stroke risk was significantly elevated at concentrations ranging from 2 to 4 pCi/L (p = 0.0004), that is, below the United States Environmental Protection Agency Radon Action Level for mitigation (4 pCi/L). Associations were slightly stronger for ischemic (especially cardioembolic, small vessel occlusive, and large artery atherosclerotic) than hemorrhagic stroke, but otherwise robust in sensitivity analyses. DISCUSSION: Radon exposure is associated with moderately increased stroke risk among middle-aged and older women in the United States, suggesting that promulgation of a lower Radon Action Level may help reduce the domestic impact of cerebrovascular disease on public health.

Improving polygenic risk prediction in admixed populations by explicitly modeling ancestral-differential effects via GAUDI.

Polygenic risk scores (PRS) have shown successes in clinics, but most PRS methods focus only on participants with distinct primary continental ancestry without accommodating recently-admixed individuals with mosaic continental ancestry backgrounds for different segments of their genomes. Here, we develop GAUDI, a novel penalized-regression-based method specifically designed for admixed individuals. GAUDI explicitly models ancestry-differential effects while borrowing information across segments with shared ancestry in admixed genomes. We demonstrate marked advantages of GAUDI over other methods through comprehensive simulation and real data analyses for traits with associated variants exhibiting ancestral-differential effects. Leveraging data from the Women's Health Initiative study, we show that GAUDI improves PRS prediction of white blood cell count and C-reactive protein in African Americans by > 64% compared to alternative methods, and even outperforms PRS-CSx with large European GWAS for some scenarios. We believe GAUDI will be a valuable tool to mitigate disparities in PRS performance in admixed individuals.

Mortality After Major Cardiovascular Events in Survivors of Childhood Cancer.

BACKGROUND: Adult survivors of childhood cancer are at risk for cardiovascular events. OBJECTIVES: In this study, we sought to determine the risk for mortality after a major cardiovascular event among childhood cancer survivors compared with noncancer populations. METHODS: All-cause and cardiovascular cause-specific mortality risks after heart failure (HF), coronary artery disease (CAD), or stroke were compared among survivors and siblings in the Childhood Cancer Survivor Study (CCSS) and participants in the Coronary Artery Risk Development in Young Adults (CARDIA) study. Cox proportional hazard regression models were used to estimate HRs and 95% CIs between groups, adjusted for demographic and clinical factors. RESULTS: Among 25,658 childhood cancer survivors (median age at diagnosis 7 years, median age at follow-up or death 38 years) and 5,051 siblings, 1,780 survivors and 91 siblings had a cardiovascular event. After HF, CAD, and stroke, 10-year all-cause mortalities were 30% (95% CI: 26%-33%), 36% (95% CI: 31%-40%), and 29% (95% CI: 24%-33%), respectively, among survivors vs 14% (95% CI: 0%-25%), 14% (95% CI: 2%-25%), and 4% (95% CI: 0%-11%) among siblings. All-cause mortality risks among childhood cancer survivors were increased after HF (HR: 7.32; 95% CI: 2.56-20.89), CAD (HR: 5.54; 95% CI: 2.37-12.93), and stroke (HR: 3.57; 95% CI: 1.12-11.37). CAD-specific mortality risk was increased (HR: 3.70; 95% CI: 1.05-13.02). Among 5,114 CARDIA participants, 345 had a major event. Although CARDIA participants were on average decades older at events (median age 57 years vs 31 years), mortality risks were similar, except that all-cause mortality after CAD was significantly increased among childhood cancer survivors (HR: 1.85; 95% CI: 1.16-2.95). CONCLUSIONS: Survivors of childhood cancer represent a population at high risk for mortality after major cardiovascular events.

Clonal Hematopoiesis and Incident Heart Failure With Preserved Ejection Fraction.

Importance: Clonal hematopoiesis of indeterminate potential (CHIP), the age-related clonal expansion of hematopoietic stem cells with leukemogenic acquired genetic variants, is associated with incident heart failure (HF). Objective: To evaluate the associations of CHIP and key gene-specific CHIP subtypes with incident HF with preserved ejection fraction (HFpEF) and reduced ejection fraction (HFrEF). Design, Setting, and Participants: This population-based cohort study included participants from 2 racially diverse prospective cohort studies with uniform HF subtype adjudication: the Jackson Heart Study (JHS) and Women's Health Initiative (WHI). JHS participants were enrolled during 2000 to 2004 and followed up through 2016. WHI participants were enrolled during 1993 to 1998 and followed up through 2022. Participants who underwent whole-genome sequencing, lacked prevalent HF at baseline, and were followed up for HF adjudication were included. Follow-up occurred over a median (IQR) of 12.0 (11.0-12.0) years in the JHS and 15.3 (9.0-22.0) years in the WHI. Statistical analysis was performed from June to December 2023. Exposures: Any CHIP and the most common gene-specific CHIP subtypes (DNMT3A and TET2 CHIP). Main Outcomes and Measures: First incident hospitalized HF events were adjudicated from hospital records and classified as HFpEF (left ventricular ejection fraction ≥50%) or HFrEF (ejection fraction <50%). Results: A total of 8090 participants were included; 2927 from the JHS (median [IQR] age, 56 [46-65] years; 1846 [63.1%] female; 2927 [100.0%] Black or African American) and 5163 from the WHI (median [IQR] age, 67 [62-72] years; 5163 [100.0%] female; 29 [0.6%] American Indian or Alaska Native, 37 [0.7%] Asian or Pacific Islander, 1383 [26.8%] Black or African American, 293 [5.7%] Hispanic or Latinx, 3407 [66.0%] non-Hispanic White, and 14 [0.3%] with other race and ethnicity). The multivariable-adjusted hazard ratio (HR) for composite CHIP and HFpEF was 1.28 (95% CI, 0.93-1.76; P = .13), and for CHIP and HFrEF it was 0.79 (95% CI, 0.49-1.25; P = .31). TET2 CHIP was associated with HFpEF in both cohorts (meta-analyzed HR, 2.35 [95% CI, 1.34 to 4.11]; P = .003) independent of cardiovascular risk factors and coronary artery disease. Analyses stratified by C-reactive protein (CRP) in the WHI found an increased risk of incident HFpEF in individuals with CHIP and CRP greater than or equal to 2 mg/L (HR, 1.94 [95% CI, 1.20-3.15]; P = .007), but not in those with CHIP and CRP less than 2 mg/L or those with CRP greater than or equal to 2 mg/L without CHIP, when compared with participants without CHIP and CRP less than 2 mg/L. Conclusions and Relevance: In this cohort study, TET2 CHIP was an independent risk factor associated with incident HFpEF. This finding may have implications for the prevention and management of HFpEF, including development of targeted therapies.

Multi-ancestry polygenic risk scores for venous thromboembolism.

Venous thromboembolism (VTE) is a significant contributor to morbidity and mortality, with large disparities in incidence rates between Black and White Americans. Polygenic risk scores (PRSs) limited to variants discovered in genome-wide association studies in European-ancestry samples can identify European-ancestry individuals at high risk of VTE. However, there is limited evidence on whether high-dimensional PRS constructed using more sophisticated methods and more diverse training data can enhance the predictive ability and their utility across diverse populations. We developed PRSs for VTE using summary statistics from the International Network against Venous Thrombosis (INVENT) consortium GWAS meta-analyses of European- (71,771 cases and 1,059,740 controls) and African-ancestry samples (7,482 cases and 129,975 controls). We used LDpred2 and PRSCSx to construct ancestry-specific and multi-ancestry PRSs and evaluated their performance in an independent European- (6,261 cases and 88,238 controls) and African-ancestry sample (1,385 cases and 12,569 controls). Multi-ancestry PRSs with weights tuned in European- and African-ancestry samples, respectively, outperformed ancestry-specific PRSs in European- (PRSCSXEUR: AUC=0.61 (0.60, 0.61), PRSCSX_combinedEUR: AUC=0.61 (0.60, 0.62)) and African-ancestry test samples (PRSCSXAFR: AUC=0.58 (0.57, 0.6), PRSCSX_combined AFR: AUC=0.59 (0.57, 0.60)). The highest fifth percentile of the best-performing PRS was associated with 1.9-fold and 1.68-fold increased risk for VTE among European- and African-ancestry subjects, respectively, relative to those in the middle stratum. These findings suggest that the multi-ancestry PRS may be used to identify individuals at highest risk for VTE and provide guidance for the most effective treatment strategy across diverse populations.

Radon Exposure, Clonal Hematopoiesis, and Stroke Susceptibility in the Women's Health Initiative.

BACKGROUND AND OBJECTIVES: Studies suggest that clonal hematopoiesis of indeterminate potential (CHIP) may increase risk of hematologic malignancy and cardiovascular disease, including stroke. However, few studies have investigated plausible environmental risk factors for CHIP such as radon, despite the climate-related increases in and documented infrequency of testing for this common indoor air pollutant.The purpose of this study was to estimate the risk of CHIP related to radon, an established environmental mutagen. METHODS: We linked geocoded addresses of 10,799 Women's Health Initiative Trans-Omics for Precision Medicine (WHI TOPMed) participants to US Environmental Protection Agency-predicted, county-level, indoor average screening radon concentrations, categorized as follows: Zone 1 (>4 pCi/L), Zone 2 (2-4 pCi/L), and Zone 3 (<2 pCi/L). We defined CHIP as the presence of one or more leukemogenic driver mutations with variant allele frequency >0.02. We identified prevalent and incident ischemic and hemorrhagic strokes; subtyped ischemic stroke using Trial of ORG 10172 in Acute Stroke Treatment (TOAST) criteria; and then estimated radon-related risk of CHIP as an odds ratio (OR) and 95% CI using multivariable-adjusted, design-weighted logistic regression stratified by age, race/ethnicity, smoking status, and stroke type/subtype. RESULTS: The percentages of participants with CHIP in Zones 1, 2, and 3 were 9.0%, 8.4%, and 7.7%, respectively (ptrend = 0.06). Among participants with ischemic stroke, Zones 2 and 1 were associated with higher estimated risks of CHIP relative to Zone 3: 1.39 (1.15-1.68) and 1.46 (1.15-1.87), but not among participants with hemorrhagic stroke: 0.98 (0.68-1.40) and 1.03 (0.70-1.52), or without stroke: 1.04 (0.74-1.46) and 0.95 (0.63-1.42), respectively (pinteraction = 0.03). Corresponding estimates were particularly high among TOAST-subtyped cardioembolism: 1.78 (1.30-2.47) and 1.88 (1.31-2.72), or other ischemic etiologies: 1.37 (1.06-1.78) and 1.50 (1.11-2.04), but not small vessel occlusion: 1.05 (0.74-1.49) and 1.00 (0.68-1.47), respectively (pinteraction = 0.10). Observed patterns of association among strata were insensitive to attrition weighting, ancestry adjustment, prevalent stroke exclusion, separate analysis of DNMT3A driver mutations, and substitution with 3 alternative estimates of radon exposure. DISCUSSION: The robust elevation of radon-related risk of CHIP among postmenopausal women who develop incident cardioembolic stroke is consistent with a potential role of somatic genomic mutation in this societally burdensome form of cerebrovascular disease, although the mechanism has yet to be confirmed.

Mosaic chromosomal alterations in blood across ancestries using whole-genome sequencing.

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.

Whole Genome Sequencing Based Analysis of Inflammation Biomarkers in the Trans-Omics for Precision Medicine (TOPMed) Consortium.

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.

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.

Time-to-Event Genome-Wide Association Study for Incident Cardiovascular Disease in People with Type 2 Diabetes Mellitus.

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 (P<5.0×10-8): rs147138607 (intergenic variant between CACNA1E and ZNF648) with a hazard ratio (HR) 1.23, 95% confidence interval (CI) 1.15 - 1.32, P=3.6×10-9, rs11444867 (intergenic variant near HS3ST1) with HR 1.89, 95% CI 1.52 - 2.35, P=9.9×10-9, and rs335407 (intergenic variant between TFB1M and NOX3) HR 1.25, 95% CI 1.16 - 1.35, P=1.5×10-8. Among 204 known CAD loci, 32 were associated with incident CVD in people with T2D with P<0.05, and 5 were significant after Bonferroni correction (P<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 (P=1.0×10-16). CONCLUSIONS: The data point to novel and known genomic regions associated with incident CVD among individuals with T2D.

Genetic examination of hematological parameters in SARS-CoV-2 infection and COVID-19.

People hospitalized with COVID-19 often exhibit altered hematological traits associated with disease prognosis (e.g., lower lymphocyte and platelet counts). We investigated whether inter-individual variability in baseline hematological traits influences risk of acute SARS-CoV-2 infection or progression to severe COVID-19. We report inconsistent associations between blood cell traits with incident SARS-CoV-2 infection and severe COVID-19 in UK Biobank and the Vanderbilt University Medical Center Synthetic Derivative (VUMC SD). Since genetically determined blood cell measures better represent cell abundance across the lifecourse, we also assessed the shared genetic architecture of baseline blood cell traits on COVID-19 related outcomes by Mendelian randomization (MR) analyses. We found significant relationships between COVID-19 severity and mean sphered cell volume after adjusting for multiple testing. However, MR results differed significantly across different freezes of COVID-19 summary statistics and genetic correlation between these traits was modest (0.1), decreasing our confidence in these results. We observed overlapping genetic association signals between other hematological and COVID-19 traits at specific loci such as MAPT and TYK2. In conclusion, we did not find convincing evidence of relationships between the genetic architecture of blood cell traits and either SARS-CoV-2 infection or COVID-19 hospitalization, though we do see evidence of shared signals at specific loci.

Age at Menopause, Leukocyte Telomere Length, and Coronary Artery Disease in Postmenopausal Women.

BACKGROUND: Premature menopause is a risk factor for accelerated cardiovascular aging, but underlying mechanisms remain incompletely understood. This study investigated the role of leukocyte telomere length (LTL), a marker of cellular aging and genomic instability, in the association of premature menopause with cardiovascular disease. METHODS: Participants from the UK Biobank and Women's Health Initiative with complete reproductive history and LTL measurements were included. Primary analyses tested the association between age at menopause and LTL using multivariable-adjusted linear regression. Secondary analyses stratified women by history of gynecologic surgery. Mendelian randomization was used to infer causal relationships between LTL and age at natural menopause. Multivariable-adjusted Cox regression and mediation analyses tested the joint associations of premature menopause and LTL with incident coronary artery disease. RESULTS: This study included 130 254 postmenopausal women (UK Biobank: n=122 224; Women's Health Initiative: n=8030), of whom 4809 (3.7%) had experienced menopause before age 40. Earlier menopause was associated with shorter LTL (meta-analyzed ß=-0.02 SD/5 years of earlier menopause [95% CI, -0.02 to -0.01]; P=7.2×10-12). This association was stronger and significant in both cohorts for women with natural/spontaneous menopause (meta-analyzed ß=-0.04 SD/5 years of earlier menopause [95% CI, -0.04 to -0.03]; P<2.2×10-16) and was independent of hormone therapy use. Mendelian randomization supported a causal association of shorter genetically predicted LTL with earlier age at natural menopause. LTL and age at menopause were independently associated with incident coronary artery disease, and mediation analyses indicated small but significant mediation effects of LTL in the association of menopausal age with coronary artery disease. CONCLUSIONS: Earlier age at menopause is associated with shorter LTL, especially among women with natural menopause. Accelerated telomere shortening may contribute to the heightened cardiovascular risk associated with premature menopause.

Traditional and novel risk factors for incident aortic stenosis in community-dwelling older adults.

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

Gene Expression and Splicing QTL Analysis of Blood Cells in African American Participants from the Jackson Heart Study.

Most gene expression and alternative splicing quantitative trait loci (eQTL/sQTL) studies have been biased toward European ancestry individuals. Here, we performed eQTL and sQTL analysis using TOPMed whole genome sequencing-derived genotype data and RNA sequencing data from stored peripheral blood mononuclear cells in 1,012 African American participants from the Jackson Heart Study (JHS). At a false discovery rate (FDR) of 5%, we identified 4,798,604 significant eQTL-gene pairs, covering 16,538 unique genes; and 5,921,368 sQTL-gene-cluster pairs, covering 9,605 unique genes. About 31% of detected eQTL and sQTL variants with a minor allele frequency (MAF) > 1% in JHS were rare (MAF < 0.1%), and therefore unlikely to be detected, in European ancestry individuals. We also generated 17,630 eQTL credible sets and 24,525 sQTL credible sets for genes (gene-clusters) with lead QTL p < 5e-8. Finally, we created an open database, which is freely available online, allowing fast query and bulk download of our QTL results.

Unraveling Neuro-Proteogenomic Landscape and Therapeutic Implications for Human Behaviors and Psychiatric Disorders.

Understanding the genetic basis of neuro-related proteins is essential for dissecting the molecular basis of human behavioral traits and the disease etiology of neuropsychiatric disorders. Here, the SCALLOP Consortium conducted a genome-wide association meta-analysis of over 12,500 individuals for 184 neuro-related proteins in human plasma. The analysis identified 117 cis-regulatory protein quantitative trait loci (cis-pQTL) and 166 trans-pQTL. The mapped pQTL capture on average 50% of each protein's heritability. Mendelian randomization analyses revealed multiple proteins showing potential causal effects on neuro-related traits such as sleeping, smoking, feelings, alcohol intake, mental health, and psychiatric disorders. Integrating with established drug information, we validated 13 out of 13 matched combinations of protein targets and diseases or side effects with available drugs, while suggesting hundreds of re-purposing and new therapeutic targets. This consortium effort provides a large-scale proteogenomic resource for biomedical research on human behaviors and other neuro-related phenotypes.

Whole Genome Analysis of Venous Thromboembolism: the Trans-Omics for Precision Medicine Program.

BACKGROUND: Risk for venous thromboembolism has a strong genetic component. Whole genome sequencing from the TOPMed program (Trans-Omics for Precision Medicine) allowed us to look for new associations, particularly rare variants missed by standard genome-wide association studies. METHODS: The 3793 cases and 7834 controls (11.6% of cases were individuals of African, Hispanic/Latino, or Asian ancestry) were analyzed using a single variant approach and an aggregate gene-based approach using our primary filter (included only loss-of-function and missense variants predicted to be deleterious) and our secondary filter (included all missense variants). RESULTS: Single variant analyses identified associations at 5 known loci. Aggregate gene-based analyses identified only PROC (odds ratio, 6.2 for carriers of rare variants; P=7.4×10-14) when using our primary filter. Employing our secondary variant filter led to a smaller effect size at PROC (odds ratio, 3.8; P=1.6×10-14), while excluding variants found only in rare isoforms led to a larger one (odds ratio, 7.5). Different filtering strategies improved the signal for 2 other known genes: PROS1 became significant (minimum P=1.8×10-6 with the secondary filter), while SERPINC1 did not (minimum P=4.4×10-5 with minor allele frequency <0.0005). Results were largely the same when restricting the analyses to include only unprovoked cases; however, one novel gene, MS4A1, became significant (P=4.4×10-7 using all missense variants with minor allele frequency <0.0005). CONCLUSIONS: Here, we have demonstrated the importance of using multiple variant filtering strategies, as we detected additional genes when filtering variants based on their predicted deleteriousness, frequency, and presence on the most expressed isoforms. Our primary analyses did not identify new candidate loci; thus larger follow-up studies are needed to replicate the novel MS4A1 locus and to identify additional rare variation associated with venous thromboembolism.