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 Chronic obstructive pulmonary disease (COPD) is associated with musculoskeletal comorbidities, including cachexia. Weight loss (WL) is the major criterion for cachexia and increases risk for mortality in COPD. Risk factors for WL in COPD are incompletely understood. We performed this whole genome sequencing (WGS) analysis to identify genetic risk variants for WL in COPD. Methods We studied 16 972 participants from the Trans‐Omics for Precision Medicine (TOPMed) Initiative and All of Us Research Program. COPD was diagnosed using spirometry in TOPMed, while diagnosis codes were used in All of Us. WL was defined as WL ≥ 5% or a final body mass index (BMI) < 20 kg/m 2 . WGS data came from white blood cells in all cohorts. Single‐variant testing was conducted on both race‐ and study‐stratified cohorts and in a cosmopolitan, ancestry‐independent manner using GENESIS in TOPMed and SAIGE in All of Us. SAIGE‐GENE+ gene‐based analyses were performed on race‐stratified and cosmopolitan cohorts. Single variant meta‐analyses were conducted using METAL within (B/AA and NHW analyses of Black/African–American and non‐Hispanic white participants, respectively) and across racial groups (COSMO). Rare variant gene‐based results were combined using Fisher's method. Transcriptomic effects were predicted using MetaXcan. We used the GWAS Catalogue to analyse for colocalization with other related traits. Results Two single variants were associated with WL in COPD among All of Us participants: one intronic variant in HCN1 in Black/African–American participants (chr5:45271359:TACACAC:T, odds ratio with 95% confidence interval (OR (CI 95 )) = 2.43(1.78–3.31), p = 1.95 × 10 −8 ) and one intergenic variant between PPP4R2 and PDZRN3 in the cosmopolitan and NHW cohorts (chr3:73345901:A:G, OR (CI 95 ) = 0.21(0.12–0.35), p = 8.84 × 10 −9 in cosmopolitan and OR (CI 95 ) = 0.18(0.10–0.33), p = 9.44 × 10 −9 in NHW). Single‐variant meta‐analysis identified two loci associated with WL in COPD: five variants within DRAIC in the B/AA meta‐analysis (lead variant chr15:69571341:A:G, OR (CI 95 ) = 1.37(1.23–1.51), p = 1.29 × 10 −9 ) and two intronic variants within RFX3 in the cosmopolitan meta‐analysis (lead variant chr9:3390983:T:C, OR (CI 95 ) = 1.50(1.31–1.73), p = 1.06 × 10 −8 ). Rare variants within RNU6‐565P (NHW analysis in All of Us; p = 2.83 × 10 −7 ) and LOC339298 (B/AA analysis in All of Us, p = 1.85 × 10 −6 ) were associated with WL in COPD. The RNU‐565P signal remained significant after combination with TOPMed results (p combined = 1.23 × 10 −6 ). MetaXcan predicted differential expression of LNC00959 in visceral adipose tissue ( p = 1.16 × 10 −6 in COSMO analysis). Colocalization analyses identified genomic associations between BMI and variants in or near DRAIC , RFX3 , PDZRN3 , and LINC00959 . Conclusions In the first WGS analysis of WL in COPD, we have identified seven novel loci. Further characterization of these loci will validate our findings and improve our understanding of the molecular pathophysiology of this condition.

Objective To establish a framework for validating candidate biomarkers of cerebral small vessel diseases (SVD) associated with cognitive impairment and characterize individuals enrolled by the MarkVCID2 consortium under this framework. Methods Participants age 60 to 90 years were enrolled across 17 MarkVCID2 sites. Recruitment was targeted to enrich in cognitive symptoms (mild dementia, mild cognitive impairment, subjective cognitive decline), defined risk factors (diabetes mellitus, advanced hypertension), and Black/African American, White, and Hispanic/Latino subgroups. Enrolled participants underwent baseline visits that included cognitive testing, multimodal magnetic resonance imaging (MRI), and biofluid collection. Provisional risk for SVD‐related cognitive decline was estimated primarily by baseline cognitive symptoms plus SVD risk factors. Adjudicated risk status was estimated by cognitive symptoms plus presence of moderate‐to‐severe white matter hyperintensities, microbleeds, or lacunes on baseline MRI. Results MarkVCID2 enrolled 1883 individuals age 73.4 ± 7.5 years, 65.0% female, 24.2% Hispanic, and 27.1% non‐Hispanic Black. Among enrollees, 44.8% were provisionally designated high‐risk. After baseline MRI, 48.5% were categorized as adjudicated high‐risk status, with substantial recategorization both from low‐ to high‐risk (primarily because of MRI lesions without SVD risk factors) and high‐ to low‐risk (primarily suspected cognitive impairment at screening not confirmed by baseline testing). Interpretation MarkVCID2 baseline data indicate successful enrollment of diverse individuals enriched in factors associated with SVD‐related cognitive decline. Changes over 3 years of longitudinal follow‐up will be analyzed to validate the candidate biomarkers for 2 projected contexts of use: subject selection (identifying likelihood of future SVD progression) and study outcome (efficiently measuring SVD progression). ANN NEUROL 2026;99:449–458

BACKGROUND: Heart failure (HF) is a heterogeneous syndrome with high mortality. The need for a new taxonomy of HF is recognized; up to now, such phenomapping efforts have primarily used clinical data. Proteomics offers potential for more precise phenotypic identification and mechanistic insights. However, few phenomapping studies have used this approach, and all have focused on targeted cardiovascular proteomics panels and a restricted HF ejection fraction group. METHODS: We measured over 7000 plasma proteins in a population-based cohort of 1351 patients with HF, used k-means clustering to identify distinct phenogroups, and compared their clinical characteristics and all-cause mortality. RESULTS: Three proteomics-defined phenogroups were identified, with substantial differences in survival (phenogroup 1 5-year survival probability, 65% [95% CI, 61%–68%]; phenogroup 2, 45% [40%–51%]; phenogroup 3, 26% [22%–30%]), independent of clinical characteristics. Phenogroups also exhibited differences in several measures suggesting poorer health, including NT-proBNP (N-terminal pro-B-type natriuretic peptide), kidney function, and Meta-Analysis Global Group in Chronic Heart Failure scores, but did not differ by ejection fraction or New York Heart Association class. CONCLUSIONS: Our study demonstrates that molecular phenomapping can stratify patients with HF into distinct subgroups that go beyond predefined clinical classifications.

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 &lt; 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 &lt; 5.0 x 10-8. For rare variant set (MAF &lt;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 &lt; 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.

AbstractBackgroundCerebral small vessel disease (cSVD), as defined by neuroimaging characteristics such as white matter hyperintensities (WMHs), cerebral microhemorrhages (CMHs), and lacunar infarcts, is highly prevalent and has been associated with dementia risk and other clinical sequelae. Although risk factors for cSVD have been identified, little is known about the biological processes and molecular mediators that influence cSVD development and progression.MethodsWithin the Atherosclerosis Risk in Communities (ARIC) study, we used SomaScan Multiplexed Proteomic technology to relate 4,877 plasma proteins to concurrently measured MRI‐defined cSVD characteristics, including WMHs, CMHs, and lacunar infarcts, in late‐life (n=1508; mean age: 76). Candidate protein associations with WMH volume were replicated in the Cardiovascular Health Study (n=765; mean age: 73) and Baltimore Longitudinal Study of Aging (n=902; mean age: 66). Cohort‐specific effect estimates were pooled using a fixed‐effect, inverse variance weighted meta‐analysis. Within 4228 non‐demented older adults in ARIC, we examined the association of WMH‐associated proteins with eight‐year dementia risk (641 incident dementia cases) using multivariate Cox proportional hazard models. All multivariate regression models were adjusted for demographic, kidney function, and cardiovascular risk variables.ResultsThis proteome‐wide analysis of older adults identified 13 WMH‐associated plasma proteins involved in synaptic function, endothelial integrity, and angiogenesis (Figures 1a, 2a). Nine of these proteins replicated in one or more external cohort (Figures 1b, 2b) and two proteins, oligodendrocyte myelin glycoprotein (OMG) and neuronal pentraxin receptor (NPTXR), remained associated with late‐life WMH volume when measured nearly two decades earlier during midlife (Figures 1c, 2c). Eight of the 13 WMH‐associated proteins were associated with lacunar infarcts, 6 with subcortical CMHs, and 4 with lobar CMHs (Figures 1d, 2d). Although none of the genes coding for WMH‐associated proteins have been nominated by WMH GWAS, 7 of these genes have been previously linked to cardiometabolic traits known to influence cSVD, such as blood pressure and cholesterol (Figures 1e, 2e). Finally, 11 candidate proteins were associated with incident dementia, underscoring their clinical relevance (Figures 1f, 2f).ConclusionsThis large‐scale proteomic analysis within longitudinal community‐based samples identified several plasma proteins, including OMG and NPTXR, as top candidate biomarkers for elevated WMH volume and its clinical manifestations.

ImportancePersistent symptoms and disability following SARS-CoV-2 infection, known as post–COVID-19 condition or “long COVID,” are frequently reported and pose a substantial personal and societal burden.ObjectiveTo determine time to recovery following SARS-CoV-2 infection and identify factors associated with recovery by 90 days.Design, Setting, and ParticipantsFor this prospective cohort study, standardized ascertainment of SARS-CoV-2 infection was conducted starting in April 1, 2020, across 14 ongoing National Institutes of Health–funded cohorts that have enrolled and followed participants since 1971. This report includes data collected through February 28, 2023, on adults aged 18 years or older with self-reported SARS-CoV-2 infection.ExposurePreinfection health conditions and lifestyle factors assessed before and during the pandemic via prepandemic examinations and pandemic-era questionnaires.Main Outcomes and MeasuresProbability of nonrecovery by 90 days and restricted mean recovery times were estimated using Kaplan-Meier curves, and Cox proportional hazards regression was performed to assess multivariable-adjusted associations with recovery by 90 days.ResultsOf 4708 participants with self-reported SARS-CoV-2 infection (mean [SD] age, 61.3 [13.8] years; 2952 women [62.7%]), an estimated 22.5% (95% CI, 21.2%-23.7%) did not recover by 90 days post infection. Median (IQR) time to recovery was 20 (8-75) days. By 90 days post infection, there were significant differences in restricted mean recovery time according to sociodemographic, clinical, and lifestyle characteristics, particularly by acute infection severity (outpatient vs critical hospitalization, 32.9 days [95% CI, 31.9-33.9 days] vs 57.6 days [95% CI, 51.9-63.3 days]; log-rank P &amp;lt; .001). Recovery by 90 days post infection was associated with vaccination prior to infection (hazard ratio [HR], 1.30; 95% CI, 1.11-1.51) and infection during the sixth (Omicron variant) vs first wave (HR, 1.25; 95% CI, 1.06-1.49). These associations were mediated by reduced severity of acute infection (33.4% and 17.6%, respectively). Recovery was unfavorably associated with female sex (HR, 0.85; 95% CI, 0.79-0.92) and prepandemic clinical cardiovascular disease (HR, 0.84; 95% CI, 0.71-0.99). No significant multivariable-adjusted associations were observed for age, educational attainment, smoking history, obesity, diabetes, chronic kidney disease, asthma, chronic obstructive pulmonary disease, or elevated depressive symptoms. Results were similar for reinfections.Conclusions and RelevanceIn this cohort study, more than 1 in 5 adults did not recover within 3 months of SARS-CoV-2 infection. Recovery within 3 months was less likely in women and those with preexisting cardiovascular disease and more likely in those with COVID-19 vaccination or infection during the Omicron variant wave.

Background The relationship between mitochondrial DNA copy number (mtDNA CN) and cardiovascular disease remains elusive. Methods and Results We performed cross‐sectional and prospective association analyses of blood‐derived mtDNA CN and cardiovascular disease outcomes in 27 316 participants in 8 cohorts of multiple racial and ethnic groups with whole‐genome sequencing. We also performed Mendelian randomization to explore causal relationships of mtDNA CN with coronary heart disease (CHD) and cardiometabolic risk factors (obesity, diabetes, hypertension, and hyperlipidemia). P <0.01 was used for significance. We validated most of the previously reported associations between mtDNA CN and cardiovascular disease outcomes. For example, 1‐SD unit lower level of mtDNA CN was associated with 1.08 (95% CI, 1.04–1.12; P <0.001) times the hazard for developing incident CHD, adjusting for covariates. Mendelian randomization analyses showed no causal effect from a lower level of mtDNA CN to a higher CHD risk (β=0.091; P =0.11) or in the reverse direction (β=−0.012; P =0.076). Additional bidirectional Mendelian randomization analyses revealed that low‐density lipoprotein cholesterol had a causal effect on mtDNA CN (β=−0.084; P <0.001), but the reverse direction was not significant ( P =0.059). No causal associations were observed between mtDNA CN and obesity, diabetes, and hypertension, in either direction. Multivariable Mendelian randomization analyses showed no causal effect of CHD on mtDNA CN, controlling for low‐density lipoprotein cholesterol level ( P =0.52), whereas there was a strong direct causal effect of higher low‐density lipoprotein cholesterol on lower mtDNA CN, adjusting for CHD status (β=−0.092; P <0.001). Conclusions Our findings indicate that high low‐density lipoprotein cholesterol may underlie the complex relationships between mtDNA CN and vascular atherosclerosis.

Abstract Background Heterochronic parabiosis has identified growth differentiation factor (GDF)-11 as a potential means of cardiac rejuvenation, but findings have been inconsistent. A major barrier has been lack of assay specificity for GDF-11 and its homolog GDF-8. Methods We tested the hypothesis that GDF-11 and GDF-8, and their major antagonists follistatin and follistatin-like (FSTL)-3, are associated with incident heart failure (HF) and its subtypes in elders. Based on validation experiments, we used liquid chromatography–tandem mass spectrometry to measure total serum GDF-11 and GDF-8, along with follistatin and FSTL-3 by immunoassay, in 2 longitudinal cohorts of older adults. Results In 2 599 participants (age 75.2 ± 4.3) followed for 10.8 ± 5.6 years, 721 HF events occurred. After adjustment, neither GDF-11 (HR per doubling: 0.93 [0.67, 1.30]) nor GDF-8 (HR: 1.02 per doubling [0.83, 1.27]) was associated with incident HF or its subtypes. Positive associations with HF were detected for follistatin (HR: 1.15 [1.00, 1.32]) and FLST-3 (HR: 1.38 [1.03, 1.85]), and with HF with preserved ejection fraction for FSTL-3 (HR: 1.77 [1.03, 3.02]). (All HRs per doubling of biomarker.) FSTL-3 associations with HF appeared stronger at higher follistatin levels and vice versa, and also for men, Blacks, and lower kidney function. Conclusions Among older adults, serum follistatin and FSTL-3, but not GDF-11 or GDF-8, were associated with incident HF. These findings do not support the concept that low serum levels of total GDF-11 or GDF-8 contribute to HF late in life, but do implicate transforming growth factor-β superfamily pathways as potential therapeutic targets.