Citations

Most genetic variants associated with complex traits are hypothesized to regulate gene expression. To understand the genetics underlying gene expression variability, we characterized 14,324 RNA-sequencing samples from the Trans-Omics for Precision Medicine program and performed expression and splicing quantitative trait locus (e/sQTL) analyses in six tissues and cell types, including whole blood ( n = 6454) and lung ( n = 1291). We detected tens of thousands of secondary cis-e/sQTLs, showing that secondary cis-e/sQTL discovery remains unsaturated. We fine-mapped UK Biobank–derived genome-wide association study (GWAS) signals from 164 traits and identified e/sQTL colocalizations for 10,611 GWAS signals, including 7096 that colocalize with secondary e/sQTLs. Our results suggest that even larger e/sQTL analyses will uncover additional secondary e/sQTLs, further benefiting GWAS interpretation.

Abstract Background Heart failure with preserved ejection fraction (HFpEF) is a heterogeneous syndrome comprising multiple pathophysiological phenotypes. HFpEF trials have largely enrolled diverse populations and reported average treatment effects, consistently yielding neutral results that may obscure drug-specific benefits within distinct subgroups. To address this issue, we employ an interaction-based that incorporates treatment-by-variable interactions to uncover drug-specific responses. Methods We leveraged four HFpEF clinical trials (TOPCAT, RELAX, NEAT-HFpEF, INDIE-HFpEF) and developed a framework comprising two complementary approaches. The first employed a prognostic responder model to evaluate whether conventional responder definitions reflect treatment-specific benefit or instead capture favorable clinical trajectories common to both treatment and placebo groups. The second used an interaction-based individual treatment effect (ITE) modeling to identify baseline variables that modify therapy effect, distinguishing drug-specific response from prognostic phenotypes. Results Although the prognostic responder model demonstrated good discrimination, further analisys suggested it primarily captured a prognostic signal associated with favorable clinical trajectories common to both treatment and placebo arms. In contrast, the ITE model identified distinct, drug-specific effect modifiers across trials (cardiorenal-inflammatory for spironolactone (TOPCAT), NO-mediated anti-inflammatory for isosorbide mononitrate (NEAT-HFpEF), afterload-reducing for inorganic nitrite (INDIE-HFpEF), and anti–volume-overload for sildenafil (RELAX). Each ITE model demonstrated significance only within its own trial suggesting drug-specific signal. Conclusions The proposed method identifies mechanism-specific effect modifiers, and uncovers clinically meaningful heterogeneity in treatment response, which is not captured by conventional MCID-based approaches. Although exploratory, these findings support phenotype-guided therapy in HFpEF and argue for phenotype-informed trial design to enhance treatment-effect detection and therapy targeting. Graphical Abstract

ABSTRACT Obstructive sleep apnea (OSA) is a common, heritable disorder with diverse clinical presentations and etiologies. We conducted a genome-wide association meta-analysis of 492,107 individuals, including 46,028 OSA cases, identifying 14 genome-wide significant loci, eight of which are novel. Analyses adjusting for body mass index (BMI) revealed three loci independent of obesity, implicating alternative biological pathways contributing to disease. Integrative functional analyses, including chromatin interaction mapping, fine-mapping, and eQTL colocalization, prioritized candidate effector genes. Notably, implicated genes in chondrocytes were associated with craniofacial morphology in both mouse and human multivariate genotype-phenotype mapping, supporting a role for craniofacial structure in OSA risk. These findings highlight key genetic pathways (obesity-related, neurological, and craniofacial) underlying the development of OSA, offering new insights into its complex etiology.

Background: Susceptibility to cardiometabolic diseases (CDs) is shaped by genetic, socio-environmental, and lifestyle factors. We evaluated the associations of disease-specific polysocial risk scores (PSSs)— derived from social determinants of health (SDOH) and lifestyle factors—and genetic risk factors with six CDs, across self-identified race/ethnicity (SIRE) groups in the eMERGE IV cohort. Methods: The six CDs included: atrial fibrillation, coronary heart disease, chronic kidney disease, hyperlipidemia, obesity, and type 2 diabetes. Disease-specific PSSs were developed using machine-learning models based on 65 SDOH and lifestyle factors. Participants were grouped based on their SIRE as Asian, Black, Hispanic/Latino, and White. We evaluated PSS distributions across SIRE groups. Genetic risk was assessed through disease-specific polygenic risk scores (PRSs), monogenic etiologies, and family history (FamHx). High PSS and high PRS were defined as the top 5th percentile. Associations were tested using multivariable logistic regression; model performance was assessed using the C-statistic. Results: Among 20178 participants (50±15 years, 68% female, 40% non-White), the most frequently associated PSS components were physical activity, sleep, self-perceived health, and income. PSS distributions varied across SIRE groups, with Whites having the lowest and Black and Hispanic/Latinos the highest scores for most CDs (Figure 1). All PSSs were significantly associated with their respective CDs, independent of confounders (OR per 1 SD increase: 1.19–2.68; Figure 2-A). No significant interactions were observed between PSSs and PRSs or monogenic etiologies. However, significant interactions were noted between PSSs and FamHx for CKD and obesity. PSS and PRS effects on the odds of CDs were independent (Figure 2-B) and additive (Figure 2-C). Incorporating PSSs into models that included clinical and genetic risk factors improved predictive performance of risk prediction models (Figure 3). Conclusion: PSSs, reflecting socio-environmental and lifestyle disadvantage, varied significantly across SIRE groups and were associated with risk of six CDs, independent of and additive to PRS. Given the substantial variation in socio-environmental and lifestyle risks across SIRE groups, integrating these factors into risk prediction models for CDs may help address health disparities and enhance preventive efforts, especially in marginalized populations.

Abstract Background:Chronic obstructive pulmonary disease (COPD) is a major global health burden, with genome-wide association studies explaining only a fraction of its heritability. Recent advancements, including large-scale, multi-ancestry whole-genome sequencing and a focus on severe, early-onset cases, offer promising avenues to identify additional COPD-related variants. Furthermore, recent studies show that adjusting for polygenic risk scores (PRS) enhances statistical power to detect rare variants, advancing efforts to understand COPD's genetic basis. Methods:We performed a whole genome sequencing association study of severe, early-onset COPD (FEV1 < 50%, age ≤ 55), FEV1 and FEV1/FVC in eight population-based studies (ARIC, CARDIA, CFS, CHS, FHS, HCHS/SOL, JHS, MESA) and seven COPD-enriched studies (COPDGene, ECLIPSE, Boston Early-Onset COPD Study, LTRC, TCGS (Poland), ICGN, SPIROMICS). We performed single-variant analyses using regenie software, applying a genome-wide significance threshold of P < 5.0 x 10-8. For rare variant set (MAF <1%) analysis, we annotated variants with loss of function or deleterious missense using WGS Annotator and prediction tools LOFTEE, SIFT, and PolyPhen2, then analyzed these variants with SKATO-ACAT method in regenie, adjusting for covariates including sex, sequencing center, and principal components of ancestry. We additionally evaluated adjusting the analyses for PRSs for lung function traits. Results: We analyzed 53,028 participants including 34,751 European, 12,508 African, 1,088 Latino American, and 674 East Asian ancestry. Across cohorts, we analyzed 5,063 severe (FEV1 < 50%) COPD cases, among whom 1,103 were early-onset (age ≤ 55). Excluding variants within 1 Mbp of known variants, single-variant analysis identified 68 novel associations (MAC ≥ 30), with 8 common associations (MAF ≥ 1%) in 6 genomic regions. Our most significant novel association for severe early-onset COPD was rs74297390 (intergenic to POM121L12, MAF = 9.6%, Beta = 0.41, P = 3.5 x 10-8). For FEV1 and FEV1/FVC ratio, variants near DKK2 and AKR1E2 were significantly associated and were replicated in an independent dataset (Shrine et al., 2019). DKK2 is an antagonist in the Wnt pathway. In a preliminary variant set analysis, CHRM3-AS1 was associated with severe early-onset COPD, GPAT with FEV1, and COL8A1I and HMCN1 with FEV1/FVC. With PRS adjustment, ST6GALNAC5 was associated with FEV1/FVC. Conclusions: In a whole genome sequencing study of lung function and severe early-onset COPD, we identified new loci for COPD-related phenotypes, including two novel loci that replicated in an independent dataset. Ongoing work involves variant set-based analyses with additional covariate adjustments, as well as PRS adjustments for binary traits.

Introduction: Patients with heart failure with preserved ejection fraction (HFpEF) have a high burden of comorbidities. To recognize the complex interplay of obesity, diabetes, and chronic kidney disease (CKD) in patients with cardiovascular disease, the cardiovascular-kidney-metabolic (CKM) syndrome was recently defined, but data on the burden of CKM in HFpEF is limited. Aim: We sought to describe the prevalence of CKM components in HFpEF leveraging data from control participants in clinical trials and NHANES. Methods: We analyzed data presented in publications from four randomized clinical trials that enrolled participants with HFpEF (all of which are available via NHLBI BioLINCC and will be included in the extant HeartShare dataset): TOPCAT, RELAX, NEAT-HFpEF, and INDIE. We abstracted baseline demographic and clinical variables from study publications to define CKM components (e.g., body mass index [BMI], diabetes status, hypertension, CKD status). We compared prevalence of CKM in these trials with a representative sample of US adults with self-reported history of heart failure (HF) from NHANES 2011-2018. Results: We included 1937 patients with HFpEF from four trials with enrollment between 2006-2016 and 715 adults representing 5.6 million adults aged ≥40 years with a self-report of HF between 2011-2018 from NHANES ( Table ). Mean age was similar among enrolled trial participants with HFpEF and those with HF in NHANES (~68-69 years). Mean BMI was in the obesity range and was similar among trial participants with HFpEF and HF in NHANES (31-35 kg/m 2 ). A significant proportion of individuals had hypertension in HFpEF trials (>80%) and in the NHANES sample with HF (89%). The proportion of stage ≥ 3 CKD in HFpEF trials (≥ 37%) was similar to the NHANES sample with HF (37%). Results are shown in Table . Conclusions: Patients with HFpEF enrolled in clinical trials and those with HF in general population are significantly enriched with CKM components that may contribute to morbidity and mortality. Similar prevalence estimates were observed despite variability in trial inclusion criteria and heterogeneity in defining HFpEF and highlight CKM as a central target for HFpEF management.

Background: Coronary heart disease (CHD) is a leading cause of death globally. Although therapy with statins decreases circulating levels of low-density lipoprotein cholesterol and the incidence of CHD, additional events occur despite statin therapy in some individuals. The genetic determinants of this residual cardiovascular risk remain unknown. Methods: We performed a 2-stage genome-wide association study of CHD events during statin therapy. We first identified 3099 cases who experienced CHD events (defined as acute myocardial infarction or the need for coronary revascularization) during statin therapy and 7681 controls without CHD events during comparable intensity and duration of statin therapy from 4 sites in the Electronic Medical Records and Genomics Network. We then sought replication of candidate variants in another 160 cases and 1112 controls from a fifth Electronic Medical Records and Genomics site, which joined the network after the initial genome-wide association study. Finally, we performed a phenome-wide association study for other traits linked to the most significant locus. Results: The meta-analysis identified 7 single nucleotide polymorphisms at a genome-wide level of significance within the LPA/PLG locus associated with CHD events on statin treatment. The most significant association was for an intronic single nucleotide polymorphism within LPA/PLG (rs10455872; minor allele frequency, 0.069; odds ratio, 1.58; 95% confidence interval, 1.35–1.86; P =2.6×10 − 10 ). In the replication cohort, rs10455872 was also associated with CHD events (odds ratio, 1.71; 95% confidence interval, 1.14–2.57; P =0.009). The association of this single nucleotide polymorphism with CHD events was independent of statin-induced change in low-density lipoprotein cholesterol (odds ratio, 1.62; 95% confidence interval, 1.17–2.24; P =0.004) and persisted in individuals with low-density lipoprotein cholesterol ≤70 mg/dL (odds ratio, 2.43; 95% confidence interval, 1.18–4.75; P =0.015). A phenome-wide association study supported the effect of this region on coronary heart disease and did not identify noncardiovascular phenotypes. Conclusions: Genetic variations at the LPA locus are associated with CHD events during statin therapy independently of the extent of low-density lipoprotein cholesterol lowering. This finding provides support for exploring strategies targeting circulating concentrations of lipoprotein(a) to reduce CHD events in patients receiving statins.