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.

Interaction molecular QTL mapping discovers cellular and environmental modifiers of genetic regulatory effects.

Bulk-tissue molecular quantitative trait loci (QTLs) have been the starting point for interpreting disease-associated variants, and context-specific QTLs show particular relevance for disease. Here, we present the results of mapping interaction QTLs (iQTLs) for cell type, age, and other phenotypic variables in multi-omic, longitudinal data from the blood of individuals of diverse ancestries. By modeling the interaction between genotype and estimated cell-type proportions, we demonstrate that cell-type iQTLs could be considered as proxies for cell-type-specific QTL effects, particularly for the most abundant cell type in the tissue. The interpretation of age iQTLs, however, warrants caution because the moderation effect of age on the genotype and molecular phenotype association could be mediated by changes in cell-type composition. Finally, we show that cell-type iQTLs contribute to cell-type-specific enrichment of diseases that, in combination with additional functional data, could guide future functional studies. Overall, this study highlights the use of iQTLs to gain insights into the context specificity of regulatory effects.

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.

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.

Low MBOAT7 expression, a genetic risk for MASH, promotes a profibrotic pathway involving hepatocyte TAZ upregulation.

BACKGROUND AND AIMS: The common genetic variant rs641738 C>T is a risk factor for metabolic dysfunction-associated steatotic liver disease and metabolic dysfunction-associated steatohepatitis (MASH), including liver fibrosis, and is associated with decreased expression of the phospholipid-remodeling enzyme MBOAT7 (LPIAT1). However, whether restoring MBOAT7 expression in established metabolic dysfunction-associated steatotic liver disease dampens the progression to liver fibrosis and, importantly, the mechanism through which decreased MBOAT7 expression exacerbates MASH fibrosis remain unclear. APPROACH AND RESULTS: We first showed that hepatocyte MBOAT7 restoration in mice with diet-induced steatohepatitis slows the progression to liver fibrosis. Conversely, when hepatocyte-MBOAT7 was silenced in mice with established hepatosteatosis, liver fibrosis but not hepatosteatosis was exacerbated. Mechanistic studies revealed that hepatocyte-MBOAT7 restoration in MASH mice lowered hepatocyte-TAZ (WWTR1), which is known to promote MASH fibrosis. Conversely, hepatocyte-MBOAT7 silencing enhanced TAZ upregulation in MASH. Finally, we discovered that changes in hepatocyte phospholipids due to MBOAT7 loss-of-function promote a cholesterol trafficking pathway that upregulates TAZ and the TAZ-induced profibrotic factor Indian hedgehog (IHH). As evidence for relevance in humans, we found that the livers of individuals with MASH carrying the rs641738-T allele had higher hepatocyte nuclear TAZ, indicating higher TAZ activity and increased IHH mRNA. CONCLUSIONS: This study provides evidence for a novel mechanism linking MBOAT7-LoF to MASH fibrosis, adds new insight into an established genetic locus for MASH, and, given the druggability of hepatocyte TAZ for MASH fibrosis, suggests a personalized medicine approach for subjects at increased risk for MASH fibrosis due to inheritance of variants that lower MBOAT7.

Polygenic scores and Mendelian randomization identify plasma proteins causally implicated in Alzheimer's disease.

Background: An increasing body of evidence suggests that neuroinflammation is one of the key drivers of late-onset Alzheimer's disease (LOAD) pathology. Due to the increased permeability of the blood-brain barrier (BBB) in older adults, peripheral plasma proteins can infiltrate the central nervous system (CNS) and drive neuroinflammation through interactions with neurons and glial cells. Because these inflammatory factors are heritable, a greater understanding of their genetic relationship with LOAD could identify new biomarkers that contribute to LOAD pathology or offer protection against it. Methods: We used a genome-wide association study (GWAS) of 90 different plasma proteins (n = 17,747) to create polygenic scores (PGSs) in an independent discovery (cases = 1,852 and controls = 1,990) and replication (cases = 799 and controls = 778) cohort. Multivariate logistic regression was used to associate the plasma protein PGSs with LOAD diagnosis while controlling for age, sex, principal components 1-2, and the number of APOE-e4 alleles as covariates. After meta-analyzing the PGS-LOAD associations between the two cohorts, we then performed a two-sample Mendelian randomization (MR) analysis using the summary statistics of significant plasma protein level PGSs in the meta-analysis as an exposure, and a GWAS of clinically diagnosed LOAD (cases = 21,982, controls = 41,944) as an outcome to explore possible causal relationships between the two. Results: We identified four plasma protein level PGSs that were significantly associated (FDR-adjusted p < 0.05) with LOAD in a meta-analysis of the discovery and replication cohorts: CX3CL1, hepatocyte growth factor (HGF), TIE2, and matrix metalloproteinase-3 (MMP-3). When these four plasma proteins were used as exposures in MR with LOAD liability as the outcome, plasma levels of HGF were inferred to have a negative causal relationship with the disease when single-nucleotide polymorphisms (SNPs) used as instrumental variables were not restricted to cis-variants (OR/95%CI = 0.945/0.906-0.984, p = 0.005). Conclusion: Our results show that plasma HGF has a negative causal relationship with LOAD liability that is driven by pleiotropic SNPs possibly involved in other pathways. These findings suggest a low transferability between PGS and MR approaches, and future research should explore ways in which LOAD and the plasma proteome may interact.

Participant-derived cell line transcriptomic analyses and mouse studies reveal a role for ZNF335 in plasma cholesterol statin response.

BACKGROUND: Statins lower circulating low-density lipoprotein cholesterol (LDLC) levels and reduce cardiovascular disease risk. Though highly efficacious in general, there is considerable inter-individual variation in statin efficacy that remains largely unexplained. METHODS: To identify novel genes that may modulate statin-induced LDLC lowering, we used RNA-sequencing data from 426 control- and 2 µM simvastatin-treated lymphoblastoid cell lines (LCLs) derived from European and African American ancestry participants of the Cholesterol and Pharmacogenetics (CAP) 40 mg/day 6-week simvastatin clinical trial (ClinicalTrials.gov Identifier: NCT00451828). We correlated statin-induced changes in LCL gene expression with plasma LDLC statin response in the corresponding CAP participants. For the most correlated gene identified (ZNF335), we followed up in vivo by comparing plasma cholesterol levels, lipoprotein profiles, and lipid statin response between wild-type mice and carriers of a hypomorphic (partial loss of function) missense mutation in Zfp335 (the mouse homolog of ZNF335). RESULTS: The statin-induced expression changes of 147 human LCL genes were significantly correlated to the plasma LDLC statin responses of the corresponding CAP participants in vivo (FDR = 5%). The two genes with the strongest correlations were zinc finger protein 335 (ZNF335 aka NIF-1, rho = 0.237, FDR-adj p = 0.0085) and CCR4-NOT transcription complex subunit 3 (CNOT3, rho = 0.233, FDR-adj p = 0.0085). Chow-fed mice carrying a hypomorphic missense (R1092W; aka bloto) mutation in Zfp335 had significantly lower non-HDL cholesterol levels than wild-type C57BL/6J mice in a sex combined model (p = 0.04). Furthermore, male (but not female) mice carrying the Zfp335R1092W allele had significantly lower total and HDL cholesterol levels than wild-type mice. In a separate experiment, wild-type mice fed a control diet for 4 weeks and a matched simvastatin diet for an additional 4 weeks had significant statin-induced reductions in non-HDLC (-43 ± 18% and -23 ± 19% for males and females, respectively). Wild-type male (but not female) mice experienced significant reductions in plasma LDL particle concentrations, while male mice carrying Zfp335R1092W allele(s) exhibited a significantly blunted LDL statin response. CONCLUSIONS: Our in vitro and in vivo studies identified ZNF335 as a novel modulator of plasma cholesterol levels and statin response, suggesting that variation in ZNF335 activity could contribute to inter-individual differences in statin clinical efficacy.

Genetic Risk and Coronary Artery Calcium in Personalizing Antihypertensive Treatment: A Pooled Cohort Analysis.

OBJECTIVE: To assess the role of the systolic blood pressure polygenic risk score (SBP-PRS) in antihypertensive treatment initiation and its comparative efficacy with coronary artery calcium (CAC) scores. PATIENTS AND METHODS: This retrospective cohort study included participants with whole genome sequencing data who underwent CAC scanning between 1971 and 2008, were free of prevalent cardiovascular disease (CVD), and were not taking antihypertensive medications. The cohort was stratified by blood pressure (BP) treatment group and SBP-PRS (low/intermediate, first and second tertiles; high, third tertile) and CAC score (0 vs >0) subgroups. The primary outcome was the first occurence of adjudicated coronary heart disease, heart failure, or stroke during 10-year follow-up. The 10-year number needed to treat (NNT) to prevent 1 event of the primary outcome was estimated. A relative risk reduction of 25% for the primary outcome based on the treatment effect of intensive control (SBP <120 mm Hg) of hypertension in SPRINT (Systolic Blood Pressure Intervention Trial) was used for estimating the NNT. RESULTS: Among the 5267 study participants, the median age was 59 years (interquartile range, 51-68 years); 2817 (53.5%) were women and 2880 (54.7%) were non-White individuals. Among 1317 individuals with elevated BP/low-risk stage 1 hypertension not recommended treatment, the 10-year incidence rate of the primary outcome was 5.6% for low/intermediate SBP-PRS and 6.3% for high SBP-PRS with NNTs of 63 and 59, respectively. Similarly, the 10-year incidence rate of the primary outcome was 2.9% for CAC score 0 and 9.7% for CAC score greater than 0, with NNTs of 117 and 37, respectively. CONCLUSION: Including genetic information in risk estimation of individuals with elevated BP/low-risk stage 1 hypertension has modest value in the initiation of antihypertensive therapy. Genetic risk and CAC both have efficacy in personalizing antihypertensive therapy.

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 risks for development of type 2 diabetes (T2D). In this study, we used genetic data on 7,227 US Hispanic/Latino participants without diabetes from the Hispanic Community Health Study/Study of Latinos (HCHS/SOL) and 400,149 non-Hispanic White participants 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 an elevated risk of T2D during follow-up compared with cluster 1 (risk ratios [RRs] 1.29 [95% CI 1.08, 1.53] and 1.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 with cluster 1. Among individuals with a healthy lifestyle, those in cluster 5 had a similar risk of T2D compared with those in cluster 1 (RR 1.03 [0.91, 1.18]). This study identified genetic subtypes of prediabetes that differed in risk of progression to T2D and in benefits from a healthy lifestyle.

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 (PRSPT) and ancestry-based continuous shrinkage priors (PRSCSx) 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 PRSPT and PRSCSx outperformed ancestry-specific PRSPT and PRSCSx across a range of tuning values. Two best-performing multi-ancestry PRSCHD (ie, PRSPTmult and PRSCSxmult) and 1 ancestry-specific (PRSCSxEUR) were taken forward for validation. PRSPTmult demonstrated the strongest association with CHD in individuals of South Asian ancestry and European ancestry (odds ratio per 1 SD [95% CI, 2.75 [2.41-3.14], 1.65 [1.59-1.72]), followed by East Asian ancestry (1.56 [1.50-1.61]), Hispanic/Latino ancestry (1.38 [1.24-1.54]), and African ancestry (1.16 [1.11-1.21]). PRSCSxmult showed the strongest associations in South Asian ancestry (2.67 [2.38-3.00]) and European ancestry (1.65 [1.59-1.71]), lower in East Asian ancestry (1.59 [1.54-1.64]), Hispanic/Latino ancestry (1.51 [1.35-1.69]), and the lowest in African 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 ancestry. This highlights the need for larger genome-wide association study datasets of underrepresented populations to enhance the performance of PRSCHD.

Sexually dimorphic DNA methylation and gene expression patterns in human first trimester placenta.

BACKGROUND: Fetal sex and placental development impact pregnancy outcomes and fetal-maternal health, but the critical timepoint of placenta establishment in first trimester is understudied in human pregnancies. METHODS: Pregnant subjects were recruited in late first trimester (weeks 10-14) at time of chorionic villus sampling, a prenatal diagnostic test. Leftover placenta tissue was collected and stored until birth outcomes were known, then DNA and RNA were isolated from singleton, normal karyotype pregnancies resulting in live births. DNA methylation was measured with the Illumina Infinium MethylationEPIC BeadChip array (n = 56). Differential methylation analysis compared 25 females versus 31 males using a generalized linear model on 743,461 autosomal probes. Gene expression sex differences were analyzed with RNA-sequencing (n = 74). An integrated analysis was performed using linear regression to correlate gene expression and DNA methylation in 51 overlapping placentas. RESULTS: Methylation analysis identified 151 differentially methylated probes (DMPs) significant at false discovery rate < 0.05, including 89 (59%) hypermethylated in females. Probe cg17612569 (GABPA, ATP5J) was the most significant CpG site, hypermethylated in males. There were 11 differentially methylated regions affected by fetal sex, with transcription factors ZNF300 and ZNF311 most significantly hypermethylated in males and females, respectively. RNA-sequencing identified 152 genes significantly sexually dimorphic at false discovery rate < 0.05. The 151 DMPs were associated with 18 genes with gene downregulation (P < 0.05) in the direction of hypermethylation, including 2 genes significant at false discovery rate < 0.05 (ZNF300 and CUB and Sushi multiple domains 1, CSMD1). Both genes, as well as Family With Sequence Similarity 228 Member A (FAM228A), showed significant correlation between DNA methylation and sexually dimorphic gene expression, though FAM228A DNA methylation was less sexually dimorphic. Comparison with other sex differences studies found that cg17612569 is male-hypermethylated across gestation in placenta and in human blood up to adulthood. CONCLUSIONS: Overall, sex dimorphic differential methylation with associated differential gene expression in the first trimester placenta is small, but there remain significant genes that may be regulated through methylation leading to differences in the first trimester placenta.

Fetal sex and placenta development affect pregnancy outcomes for both the fetus and mother throughout pregnancy, including risk of miscarriages, preterm birth, preeclampsia, and other outcomes. Epigenetics, the "overlay" of regulatory signals on DNA which affects how DNA is read, is not well understood in early pregnancy when critical placenta developments are happening that affect the rest of pregnancy. Here, we use leftover placenta biopsy samples (n = 56) donated by Cedars-Sinai patients with informed consent to learn about first trimester human placenta DNA methylation differences due to fetal sex. Out of the total 743,461 sites analyzed, we identified 151 sites significantly affected by fetal sex after correcting p-values to reduce false positives (false discovery rate < 0.05). We also performed an analysis to look at multiple sites and identified 11 regions across the genome with significant DNA methylation changes due to fetal sex. Furthermore, because DNA methylation is a regulatory mark on DNA which typically dampens gene expression, we also compared the DNA methylation sex differences to placental RNA-sequencing gene expression analysis using the same tissue from a mostly overlapping patient group (n = 74 total sequenced, n = 51 overlap). We identify 18 genes which show both significant DNA methylation differences and gene expression changes. The most significant gene was transcription factor ZNF300 with higher DNA methylation in males and reduced gene expression in males (and thus higher gene expression in females). This study identifies some sex differences that continue until later pregnancy and others that are unique to first trimester.

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.

Genome-wide association identifies novel ROP risk loci in a multiethnic cohort.

We conducted a genome-wide association study (GWAS) in a multiethnic cohort of 920 at-risk infants for retinopathy of prematurity (ROP), a major cause of childhood blindness, identifying 1 locus at genome-wide significance level (p < 5×10-8) and 9 with significance of p < 5×10-6 for ROP ≥ stage 3. The most significant locus, rs2058019, reached genome-wide significance within the full multiethnic cohort (p = 4.96×10-9); Hispanic and European Ancestry infants driving the association. The lead single nucleotide polymorphism (SNP) falls in an intronic region within the Glioma-associated oncogene family zinc finger 3 (GLI3) gene. Relevance for GLI3 and other top-associated genes to human ocular disease was substantiated through in-silico extension analyses, genetic risk score analysis and expression profiling in human donor eye tissues. Thus, we identify a novel locus at GLI3 with relevance to retinal biology, supporting genetic susceptibilities for ROP risk with possible variability by race and ethnicity.

Identifying People at High Risk for Severe Aortic Stenosis: Aortic Valve Calcium Versus Lipoprotein(a) and Low-Density Lipoprotein Cholesterol.

BACKGROUND: Aortic valve calcification (AVC), Lp(a) [lipoprotein(a)], and low-density lipoprotein cholesterol (LDL-C) are associated with severe aortic stenosis (AS). We aimed to determine which of these risk factors were most strongly associated with the risk of incident severe AS. METHODS: A total of 6792 participants from the MESA study (Multi-Ethnic Study of Atherosclerosis) had computed tomography-quantified AVC, Lp(a), and LDL-C values at MESA visit 1 (2000-2002). We calculated the absolute event rate of incident adjudicated severe AS per 1000 person-years and performed multivariable adjusted Cox proportional hazards regression. RESULTS: The mean age was 62 years old, and 47% were women. Over a median 16.7-year follow-up, the rate of incident severe AS increased exponentially with higher AVC, regardless of Lp(a) or LDL-C values. Participants with AVC=0 had a very low rate of severe AS even with elevated Lp(a) ≥50 mg/dL (<0.1/1000 person-years) or LDL-C ≥130 mg/dL (0.1/1000 person-years). AVC >0 was strongly associated with severe AS when Lp(a) <50 mg/dL hazard ratio (HR) of 33.8 (95% CI, 16.4-70.0) or ≥50 mg/dL HR of 61.5 (95% CI, 7.7-494.2) and when LDL-C <130 mg/dL HR of 31.1 (95% CI, 14.4-67.1) or ≥130 mg/dL HR of 50.2 (95% CI, 13.2-191.9). CONCLUSIONS: AVC better identifies people at high risk for severe AS compared with Lp(a) or LDL-C, and people with AVC=0 have a very low long-term rate of severe AS regardless of Lp(a) or LDL-C level. These results suggest AVC should be the preferred prognostic risk marker to identify patients at high risk for severe AS, which may help inform participant selection for future trials testing novel strategies to prevent severe AS.

Association of aortic valve calcium with dementia and stroke: The Multi-Ethnic Study of Atherosclerosis.

BACKGROUND AND AIMS: Calcific aortic valve disease is associated with increased thrombin formation, platelet activation, decreased fibrinolysis, and subclinical brain infarcts. We examined the long-term association of aortic valve calcification (AVC) with newly diagnosed dementia and incident stroke in the Multi-Ethnic Study of Atherosclerosis (MESA). METHODS: AVC was measured using non-contrast cardiac CT at Visit 1. We examined AVC as a continuous (log-transformed) and categorical variable (0, 1-99, 100-299, ≥300). Newly diagnosed dementia was adjudicated using International Classification of Disease codes. Stroke was adjudicated from medical records. We calculated absolute event rates (per 1000 person-years) and multivariable adjusted Cox proportional hazards ratios (HR). RESULTS: Overall, 6812 participants had AVC quantified with a mean age of 62.1 years old, 52.9 % were women, and the median 10-year estimated atherosclerotic cardiovascular disease (ASCVD) risk was 13.5 %. Participants with AVC >0 were older and less likely to be women compared to those with AVC=0. Over a median 16-year follow-up, there were 535 cases of dementia and 376 cases of stroke. The absolute risk of newly diagnosed dementia increased in a stepwise pattern with higher AVC scores, and stroke increased in a logarithmic pattern. In multivariable analyses, AVC was significantly associated with newly diagnosed dementia as a log-transformed continuous variable (HR 1.09; 95 % CI 1.04-1.14) and persons with AVC ≥300 had nearly a two-fold higher risk (HR 1.77; 95 % CI 1.14-2.76) compared to those with AVC=0. AVC was associated with an increased risk of stroke after adjustment for age, sex, and race/ethnicity, but not after adjustment for ASCVD risk factors. CONCLUSIONS: After multivariable adjustment, AVC >0 was significantly associated with an increased risk of newly diagnosed dementia, but not incident stroke. This suggests that AVC may be an important risk factor for the long-term risk of dementia beyond traditional ASCVD risk factors.

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.

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.

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.