Determinants of mosaic chromosomal alteration fitness.

Clonal hematopoiesis (CH) is characterized by the acquisition of a somatic mutation in a hematopoietic stem cell that results in a clonal expansion. These driver mutations can be single nucleotide variants in cancer driver genes or larger structural rearrangements called mosaic chromosomal alterations (mCAs). The factors that influence the variations in mCA fitness and ultimately result in different clonal expansion rates are not well understood. We used the Passenger-Approximated Clonal Expansion Rate (PACER) method to estimate clonal expansion rate as PACER scores for 6,381 individuals in the NHLBI TOPMed cohort with gain, loss, and copy-neutral loss of heterozygosity mCAs. Our mCA fitness estimates, derived by aggregating per-individual PACER scores, were correlated (R2 = 0.49) with an alternative approach that estimated fitness of mCAs in the UK Biobank using population-level distributions of clonal fraction. Among individuals with JAK2 V617F clonal hematopoiesis of indeterminate potential or mCAs affecting the JAK2 gene on chromosome 9, PACER score was strongly correlated with erythrocyte count. In a cross-sectional analysis, genome-wide association study of estimates of mCA expansion rate identified a TCL1A locus variant associated with mCA clonal expansion rate, with suggestive variants in NRIP1 and TERT.

LXR signaling pathways link cholesterol metabolism with risk for prediabetes and diabetes.

BACKGROUNDPreclinical studies suggest that cholesterol accumulation leads to insulin resistance. We previously reported that alterations in a monocyte cholesterol metabolism transcriptional network (CMTN) - suggestive of cellular cholesterol accumulation - were cross-sectionally associated with obesity and type 2 diabetes (T2D). Here, we sought to determine whether the CMTN alterations independently predict incident prediabetes/T2D risk, and correlate with cellular cholesterol accumulation.METHODSMonocyte mRNA expression of 11 CMTN genes was quantified among 934 Multi-Ethnic Study of Atherosclerosis (MESA) participants free of prediabetes/T2D; cellular cholesterol was measured in a subset of 24 monocyte samples.RESULTSDuring a median 6-year follow-up, lower expression of 3 highly correlated LXR target genes - ABCG1 and ABCA1 (cholesterol efflux) and MYLIP (cholesterol uptake suppression) - and not other CMTN genes, was significantly associated with higher risk of incident prediabetes/T2D. Lower expression of the LXR target genes correlated with higher cellular cholesterol levels (e.g., 47% of variance in cellular total cholesterol explained by ABCG1 expression). Further, adding the LXR target genes to overweight/obesity and other known predictors significantly improved prediction of incident prediabetes/T2D.CONCLUSIONThese data suggest that the aberrant LXR/ABCG1-ABCA1-MYLIP pathway (LAAMP) is a major T2D risk factor and support a potential role for aberrant LAAMP and cellular cholesterol accumulation in diabetogenesis.FUNDINGThe MESA Epigenomics and Transcriptomics Studies were funded by NIH grants 1R01HL101250, 1RF1AG054474, R01HL126477, R01DK101921, and R01HL135009. This work was supported by funding from NIDDK R01DK103531 and NHLBI R01HL119962.

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.

Association analysis between an epigenetic risk score and blood pressure.

Background: Epigenome-wide association studies have revealed multiple DNA methylation sites (CpGs) associated with alcohol consumption, an important lifestyle risk factor for cardiovascular diseases. Results: We generated an alcohol consumption epigenetic risk score (ERS) based on previously reported 144 alcohol-associated CpGs and examined the association of the ERS with systolic blood pressure (SBP), diastolic blood pressure (DBP), and hypertension (HTN) in 3,898 Framingham Heart Study (FHS) participants. We found an association of alcohol intake with the ERS in the meta-analysis with 0.09 units higher ERS per drink consumed per day (p < 0.0001). Cross-sectional analyses in FHS revealed that a one-unit increment of the ERS was associated with 1.93 mm Hg higher SBP (p = 4.64E-07), 0.68 mm Hg higher DBP (p = 0.006), and an odds ratio of 1.78 for HTN (p < 2E-16). Meta-analysis of the cross-sectional association of the ERS with BP traits in eight independent external cohorts (n = 11,544) showed similar relationships with blood pressure levels, i.e., a one-unit increase in ERS was associated with 0.74 (p = 0.002) and 0.50 (p = 0.0006) mm Hg higher SBP and DBP, but could not confirm the association with hypertension. Longitudinal analyses in FHS (n = 3,260) and five independent external cohorts (n = 4,021) showed that the baseline ERS was not associated with a change in blood pressure over time or with incident HTN. Conclusions: Our findings provide proof-of-concept that utilizing an ERS is a useful approach to capture the recent health consequences of lifestyle behaviors such as alcohol consumption.

Longitudinal change in mitochondrial heteroplasmy exhibits positive selection for deleterious variants.

A common feature of human aging is the acquisition of somatic mutations, and mitochondria are particularly prone to mutation due to their inefficient DNA repair and close proximity to reactive oxygen species, leading to a state of mitochondrial DNA heteroplasmy1,2. Cross-sectional studies have demonstrated that detection of heteroplasmy increases with participant age3, a phenomenon that has been attributed to genetic drift4-7. In this first large-scale longitudinal study, we measured heteroplasmy in two prospective cohorts (combined n=1405) at two timepoints (mean time between visits, 8.6 years), demonstrating that deleterious heteroplasmies were more likely to increase in variant allele fraction (VAF). We further demonstrated that increase in VAF was associated with increased risk of overall mortality. These results challenge the claim that somatic mtDNA mutations arise mainly due to genetic drift, instead demonstrating positive selection for predicted deleterious mutations at the cellular level, despite an negative impact on overall mortality.

Genetic subtypes of prediabetes, healthy lifestyle, and risk of type 2 diabetes.

Prediabetes is a heterogenous metabolic state with various risk for development of type 2 diabetes (T2D). In this study, we used genetic data on 7,227 US Hispanic/Latinos without diabetes from the Hispanic Community Health Study/Study of Latinos (HCHS/SOL) and 400,149 non-Hispanic whites without diabetes from the UK Biobank (UKBB) to calculate five partitioned polygenetic risk scores (pPRSs) representing various pathways related to T2D. Consensus clustering was performed in participants with prediabetes in HCHS/SOL (n=3,677) and UKBB (n=16,284) separately, based on these pPRSs. Six clusters of individuals with prediabetes with distinctive patterns of pPRSs and corresponding metabolic traits were identified in the HCHS/SOL, five of which were confirmed in the UKBB. Although baseline glycemic traits were similar across clusters, individuals in Cluster 5 and Cluster 6 showed elevated risk of T2D during follow-up compared to Cluster 1 (RR=1.29 [95% CI 1.08-1.53] and1.34 [1.13-1.60], respectively). Inverse associations between a healthy lifestyle score and risk of T2D were observed across different clusters, with a suggestively stronger association observed in Cluster 5 compared to Cluster 1. Among individuals with healthy lifestyle, those in Cluster 5 had a similar risk of T2D compared to those in Cluster 1 (RR=1.03 [0.91-1.18]). This study identified genetic subtypes of prediabetes which differed in risk of progression to T2D and in benefits from healthy lifestyle.

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.

Genome-wide analysis in over 1 million individuals of European ancestry yields improved polygenic risk scores for blood pressure traits.

Hypertension affects more than one billion people worldwide. Here we identify 113 novel loci, reporting a total of 2,103 independent genetic signals (P < 5 × 10-8) from the largest single-stage blood pressure (BP) genome-wide association study to date (n = 1,028,980 European individuals). These associations explain more than 60% of single nucleotide polymorphism-based BP heritability. Comparing top versus bottom deciles of polygenic risk scores (PRSs) reveals clinically meaningful differences in BP (16.9 mmHg systolic BP, 95% CI, 15.5-18.2 mmHg, P = 2.22 × 10-126) and more than a sevenfold higher odds of hypertension risk (odds ratio, 7.33; 95% CI, 5.54-9.70; P = 4.13 × 10-44) in an independent dataset. Adding PRS into hypertension-prediction models increased the area under the receiver operating characteristic curve (AUROC) from 0.791 (95% CI, 0.781-0.801) to 0.826 (95% CI, 0.817-0.836, ∆AUROC, 0.035, P = 1.98 × 10-34). We compare the 2,103 loci results in non-European ancestries and show significant PRS associations in a large African-American sample. Secondary analyses implicate 500 genes previously unreported for BP. Our study highlights the role of increasingly large genomic studies for precision health research.

Methylation patterns associated with C-reactive protein in racially and ethnically diverse populations.

Systemic low-grade inflammation is a feature of chronic disease. C-reactive protein (CRP) is a common biomarker of inflammation and used as an indicator of disease risk; however, the role of inflammation in disease is not completely understood. Methylation is an epigenetic modification in the DNA which plays a pivotal role in gene expression. In this study we evaluated differential DNA methylation patterns associated with blood CRP level to elucidate biological pathways and genetic regulatory mechanisms to improve the understanding of chronic inflammation. The racially and ethnically diverse participants in this study were included as 50% White, 41% Black or African American, 7% Hispanic or Latino/a, and 2% Native Hawaiian, Asian American, American Indian, or Alaska Native (total n = 13,433) individuals. We replicated 113 CpG sites from 87 unique loci, of which five were novel (CADM3, NALCN, NLRC5, ZNF792, and cg03282312), across a discovery set of 1,150 CpG sites associated with CRP level (p < 1.2E-7). The downstream pathways affected by DNA methylation included the identification of IFI16 and IRF7 CpG-gene transcript pairs which contributed to the innate immune response gene enrichment pathway along with NLRC5, NOD2, and AIM2. Gene enrichment analysis also identified the nuclear factor-kappaB transcription pathway. Using two-sample Mendelian randomization (MR) we inferred methylation at three CpG sites as causal for CRP levels using both White and Black or African American MR instrument variables. Overall, we identified novel CpG sites and gene transcripts that could be valuable in understanding the specific cellular processes and pathogenic mechanisms involved in inflammation.

Lac-Phe mediates the effects of metformin on food intake and body weight.

Metformin is a widely prescribed anti-diabetic medicine that also reduces body weight. There is ongoing debate about the mechanisms that mediate metformin's effects on energy balance. Here, we show that metformin is a powerful pharmacological inducer of the anorexigenic metabolite N-lactoyl-phenylalanine (Lac-Phe) in cells, in mice and two independent human cohorts. Metformin drives Lac-Phe biosynthesis through the inhibition of complex I, increased glycolytic flux and intracellular lactate mass action. Intestinal epithelial CNDP2+ cells, not macrophages, are the principal in vivo source of basal and metformin-inducible Lac-Phe. Genetic ablation of Lac-Phe biosynthesis in male mice renders animals resistant to the effects of metformin on food intake and body weight. Lastly, mediation analyses support a role for Lac-Phe as a downstream effector of metformin's effects on body mass index in participants of a large population-based observational cohort, the Multi-Ethnic Study of Atherosclerosis. Together, these data establish Lac-Phe as a critical mediator of the body weight-lowering effects of metformin.

Association analysis between an epigenetic alcohol risk score and blood pressure.

Background: Epigenome-wide association studies have revealed multiple DNA methylation sites (CpGs) associated with alcohol consumption, an important lifestyle risk factor for cardiovascular diseases. Results: We generated an alcohol consumption epigenetic risk score (ERS) based on previously reported 144 alcohol-associated CpGs and examined the association of the ERS with systolic blood pressure (SBP), diastolic blood pressure (DBP), and hypertension (HTN) in 3,898 Framingham Heart Study (FHS) participants. We found an association of alcohol intake with the ERS in the meta-analysis with 0.09 units higher ERS per drink consumed per day (p < 0.0001). Cross-sectional analyses in FHS revealed that a one-unit increment of the ERS was associated with 1.93 mm Hg higher SBP (p = 4.64E-07), 0.68 mm Hg higher DBP (p = 0.006), and an odds ratio of 1.78 for HTN (p < 2E-16). Meta-analysis of the cross-sectional association of the ERS with BP traits in eight independent external cohorts (n = 11,544) showed similar relationships with blood pressure levels, i.e., a one-unit increase in ERS was associated with 0.74 (p = 0.002) and 0.50 (p = 0.0006) mm Hg higher SBP and DBP, but could not confirm the association with hypertension. Longitudinal analyses in FHS (n = 3,260) and five independent external cohorts (n = 4,021) showed that the baseline ERS was not associated with a change in blood pressure over time or with incident HTN. Conclusions: Our findings provide proof-of-concept that utilizing an ERS is a useful approach to capture the recent health consequences of lifestyle behaviors such as alcohol consumption.

Peripheral Blood Mononuclear Cell Gene Expression Associated with Pulmonary Microvascular Perfusion: The Multi-Ethnic Study of Atherosclerosis Chronic Obstructive Pulmonary Disease Study.

Rationale Chronic obstructive pulmonary disease (COPD) and emphysema are associated with endothelial damage and altered pulmonary microvascular perfusion. Molecular mechanisms underlying these changes are poorly understood in patients due, in part, to the inaccessibility of the pulmonary vasculature. Peripheral blood mononuclear cells (PBMC) interact with the pulmonary endothelium. Objective To test the association between gene expression in PBMCs and pulmonary microvascular perfusion in COPD. Methods The Multi-Ethnic Study of Atherosclerosis (MESA) COPD Study recruited two independent samples of COPD cases and controls with 10 or more pack-years. In both samples, pulmonary microvascular blood flow, pulmonary microvascular blood volume (PMBV), and mean transit time were assessed on contrast-enhanced MRI, and PBMC gene expression was assessed by microarray. Additional replication was performed in a third sample with PMBV measures on contrast-enhanced, dual-energy CT. Differential expression analyses were adjusted for age, gender, race-ethnicity, educational attainment, height, weight, smoking status, and pack-years. Results The 79 participants in the discovery sample had mean age of 69±6 years, 44% were female, 25% were non-white, 34% were current smokers and 66% had COPD. There were large PBMC gene expression signatures associated with pulmonary microvascular perfusion traits, with several replicated in the replication sets with MRI (n=47) or dual-energy CT scan (n=157) measures. Many of the identified genes are involved in inflammatory processes, including NF-κB and chemokine signaling pathways. Conclusions PBMC gene expression in NF-κB, inflammatory and chemokine signaling pathways was associated pulmonary microvascular perfusion in COPD, potentially offering new targetable candidates for novel therapies.

Transcriptome-wide association study of the plasma proteome reveals cis and trans regulatory mechanisms underlying complex traits.

Regulation of transcription and translation are mechanisms through which genetic variants affect complex traits. Expression quantitative trait locus (eQTL) studies have been more successful at identifying cis-eQTL (within 1 Mb of the transcription start site) than trans-eQTL. Here, we tested the cis component of gene expression for association with observed plasma protein levels to identify cis- and trans-acting genes that regulate protein levels. We used transcriptome prediction models from 49 Genotype-Tissue Expression (GTEx) Project tissues to predict the cis component of gene expression and tested the predicted expression of every gene in every tissue for association with the observed abundance of 3,622 plasma proteins measured in 3,301 individuals from the INTERVAL study. We tested significant results for replication in 971 individuals from the Trans-omics for Precision Medicine (TOPMed) Multi-Ethnic Study of Atherosclerosis (MESA). We found 1,168 and 1,210 cis- and trans-acting associations that replicated in TOPMed (FDR < 0.05) with a median expected true positive rate (π1) across tissues of 0.806 and 0.390, respectively. The target proteins of trans-acting genes were enriched for transcription factor binding sites and autoimmune diseases in the GWAS catalog. Furthermore, we found a higher correlation between predicted expression and protein levels of the same underlying gene (R = 0.17) than observed expression (R = 0.10, p = 7.50 × 10-11). This indicates the cis-acting genetically regulated (heritable) component of gene expression is more consistent across tissues than total observed expression (genetics + environment) and is useful in uncovering the function of SNPs associated with complex traits.

Multi-Ancestry Polygenic Risk Score for Coronary Heart Disease Based on an Ancestrally Diverse Genome-Wide Association Study and Population-Specific Optimization.

BACKGROUND: Predictive performance of polygenic risk scores (PRS) varies across populations. To facilitate equitable clinical use, we developed PRS for coronary heart disease (CHD; PRSCHD) for 5 genetic ancestry groups. METHODS: We derived ancestry-specific and multi-ancestry PRSCHD based on pruning and thresholding and continuous shrinkage priors (polygenic risk score for CHD developed using ancestry-based continuous shrinkage methods) applied to summary statistics from the largest multi-ancestry genome-wide association study meta-analysis for CHD to date, including 1.1 million participants from 5 major genetic ancestry groups. Following training and optimization in the Million Veteran Program, we evaluated the best-performing PRSCHD in 176 988 individuals across 9 diverse cohorts. RESULTS: Multi-ancestry polygenic risk score for CHD developed using pruning and thresholding methods and polygenic risk score for CHD developed using ancestry-based continuous shrinkage methods outperformed ancestry-specific Polygenic risk score for CHD developed using pruning and thresholding methods and polygenic risk score for CHD developed using ancestry-based continuous shrinkage methods across a range of tuning values. Two best-performing multi-ancestry PRSCHD (ie, polygenic risk score for CHD developed using pruning and thresholding methods optimized using a multi-ancestry population and polygenic risk score for CHD developed using ancestry-based continuous shrinkage methods optimized using a multi-ancestry population) and 1 ancestry-specific (PRSCSxEUR) were taken forward for validation. Polygenic risk score for CHD developed using pruning and thresholding methods (PT) optimized using a multi-ancestry population demonstrated the strongest association with CHD in individuals of South Asian genetic ancestry and European genetic ancestry (odds ratio per 1 SD [95% CI, 2.75 [2.41-3.14], 1.65 [1.59-1.72]), followed by East Asian genetic ancestry (1.56 [1.50-1.61]), Hispanic/Latino genetic ancestry (1.38 [1.24-1.54]), and African genetic ancestry (1.16 [1.11-1.21]). Polygenic risk score for CHD developed using ancestry-based continuous shrinkage methods optimized using a multi-ancestry population showed the strongest associations in South Asian genetic ancestry (2.67 [2.38-3.00]) and European genetic ancestry (1.65 [1.59-1.71]), lower in East Asian genetic ancestry (1.59 [1.54-1.64]), Hispanic/Latino genetic ancestry (1.51 [1.35-1.69]), and the lowest in African genetic ancestry (1.20 [1.15-1.26]). CONCLUSIONS: The use of summary statistics from a large multi-ancestry genome-wide meta-analysis improved the performance of PRSCHD in most ancestry groups compared with single-ancestry methods. Despite the use of one of the largest and most diverse sets of training and validation cohorts to date, improvement of predictive performance was limited in African genetic ancestry. This highlights the need for larger Genome-wide association study datasets of underrepresented populations to enhance the performance of PRSCHD.

Refinement of a published gene-physical activity interaction impacting HDL-cholesterol: role of sex and lipoprotein subfractions.

Large-scale gene-environment interaction (GxE) discovery efforts often involve compromises in the definition of outcomes and choice of covariates for the sake of data harmonization and statistical power. Consequently, refinement of exposures, covariates, outcomes, and population subsets may be helpful to establish often-elusive replication and evaluate potential clinical utility. Here, we used additional datasets, an expanded set of statistical models, and interrogation of lipoprotein metabolism via nuclear magnetic resonance (NMR)-based lipoprotein subfractions to refine a previously discovered GxE modifying the relationship between physical activity (PA) and HDL-cholesterol (HDL-C). This GxE was originally identified by Kilpeläinen et al., with the strongest cohort-specific signal coming from the Women's Genome Health Study (WGHS). We thus explored this GxE further in the WGHS (N = 23,294), with follow-up in the UK Biobank (UKB; N = 281,380), and the Multi-Ethnic Study of Atherosclerosis (MESA; N = 4,587). Self-reported PA (MET-hrs/wk), genotypes at rs295849 (nearest gene: LHX1), and NMR metabolomics data were available in all three cohorts. As originally reported, minor allele carriers of rs295849 in WGHS had a stronger positive association between PA and HDL-C (pint = 0.002). When testing a range of NMR metabolites (primarily lipoprotein and lipid subfractions) to refine the HDL-C outcome, we found a stronger interaction effect on medium-sized HDL particle concentrations (M-HDL-P; pint = 1.0×10-4) than HDL-C. Meta-regression revealed a systematically larger interaction effect in cohorts from the original meta-analysis with a greater fraction of women (p = 0.018). In the UKB, GxE effects were stronger both in women and using M-HDL-P as the outcome. In MESA, the primary interaction for HDL-C showed nominal significance (pint = 0.013), but without clear differences by sex and with a greater magnitude using large, rather than medium, HDL-P as an outcome. Towards reconciling these observations, further exploration leveraging NMR platform-specific HDL subfraction diameter annotations revealed modest agreement across all cohorts in the interaction affecting medium-to-large particles. Taken together, our work provides additional insights into a specific known gene-PA interaction while illustrating the importance of phenotype and model refinement towards understanding and replicating GxEs.

Proteomic Signature of HIV-Associated Subclinical Left Atrial Remodeling and Incident Heart Failure.

Background: People living with HIV (PLWH) are at higher risk of heart failure (HF) and preceding subclinical cardiac abnormalities, including left atrial dilation, compared to people without HIV (PWOH). Hypothesized mechanisms include premature aging linked to chronic immune activation. We leveraged plasma proteomics to identify potential novel contributors to HIV-associated differences in indexed left atrial volume (LAVi) among PLWH and PWOH and externally validated identified proteomic signatures with incident HF among a cohort of older PWOH. Methods: We performed proteomics (Olink Explore 3072) on plasma obtained concurrently with cardiac magnetic resonance imaging among PLWH and PWOH in the United States. Proteins were analyzed individually and as agnostically defined clusters. Cross-sectional associations with HIV and LAVi were estimated using multivariable regression with robust variance. Among an independent general population cohort, we estimated associations between identified signatures and LAVi using linear regression and incident HF using Cox regression. Results: Among 352 participants (age 55±6 years; 25% female), 61% were PLWH (88% on ART; 73% with undetectable HIV RNA) and mean LAVi was 29±9 mL/m 2 . Of 2594 analyzed proteins, 439 were associated with HIV serostatus, independent of demographics, hepatitis C virus infection, renal function, and substance use (FDR<0.05). We identified 73 of these proteins as candidate contributors to the independent association between positive HIV serostatus and higher LAVi, enriched in tumor necrosis factor (TNF) signaling and immune checkpoint proteins regulating T cell, B cell, and NK cell activation. We identified one protein cluster associated with LAVi and HIV regardless of HIV viral suppression status, which comprised 42 proteins enriched in TNF signaling, ephrin signaling, and extracellular matrix (ECM) organization. This protein cluster and 30 of 73 individual proteins were associated with incident HF among 2273 older PWOH (age 68±9 years; 52% female; 8.5±1.4 years of follow-up). Conclusion: Proteomic signatures that may contribute to HIV-associated LA remodeling were enriched in immune checkpoint proteins, cytokine signaling, and ECM organization. These signatures were also associated with incident HF among older PWOH, suggesting specific markers of chronic immune activation, systemic inflammation, and fibrosis may identify shared pathways in HIV and aging that contribute to risk of HF.

Genome-wide association analysis of composite sleep health scores in 413,904 individuals.

Recent genome-wide association studies (GWASs) of several individual sleep traits have identified hundreds of genetic loci, suggesting diverse mechanisms. Moreover, sleep traits are moderately correlated, and together may provide a more complete picture of sleep health, while also illuminating distinct domains. Here we construct novel sleep health scores (SHSs) incorporating five core self-report measures: sleep duration, insomnia symptoms, chronotype, snoring, and daytime sleepiness, using additive (SHS-ADD) and five principal components-based (SHS-PCs) approaches. GWASs of these six SHSs identify 28 significant novel loci adjusting for multiple testing on six traits (p<8.3e-9), along with 341 previously reported loci (p<5e-08). The heritability of the first three SHS-PCs equals or exceeds that of SHS-ADD (SNP-h2=0.094), while revealing sleep-domain-specific genetic discoveries. Significant loci enrich in multiple brain tissues and in metabolic and neuronal pathways. Post GWAS analyses uncover novel genetic mechanisms underlying sleep health and reveal connections to behavioral, psychological, and cardiometabolic traits.

Association analysis of mitochondrial DNA heteroplasmic variants: methods and application.

We rigorously assessed a comprehensive association testing framework for heteroplasmy, employing both simulated and real-world data. This framework employed a variant allele fraction (VAF) threshold and harnessed multiple gene-based tests for robust identification and association testing of heteroplasmy. Our simulation studies demonstrated that gene-based tests maintained an appropriate type I error rate at α=0.001. Notably, when 5% or more heteroplasmic variants within a target region were linked to an outcome, burden-extension tests (including the adaptive burden test, variable threshold burden test, and z-score weighting burden test) outperformed the sequence kernel association test (SKAT) and the original burden test. Applying this framework, we conducted association analyses on whole-blood derived heteroplasmy in 17,507 individuals of African and European ancestries (31% of African Ancestry, mean age of 62, with 58% women) with whole genome sequencing data. We performed both cohort- and ancestry-specific association analyses, followed by meta-analysis on both pooled samples and within each ancestry group. Our results suggest that mtDNA-encoded genes/regions are likely to exhibit varying rates in somatic aging, with the notably strong associations observed between heteroplasmy in the RNR1 and RNR2 genes (p<0.001) and advance aging by the Original Burden test. In contrast, SKAT identified significant associations (p<0.001) between diabetes and the aggregated effects of heteroplasmy in several protein-coding genes. Further research is warranted to validate these findings. In summary, our proposed statistical framework represents a valuable tool for facilitating association testing of heteroplasmy with disease traits in large human populations.

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.