Reiner, Alex P.

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 RATIONALE: Chronic kidney disease (CKD) is common in the general population and is associated with inflammation and compromised immune function. Since these mechanisms are hypothesized to promote severe COVID-19, we tested associations of pre-pandemic CKD with incident COVID-19 hospitalization or death in a large, multi-ethnic, US general population-based meta-cohort of adults. METHODS: Data were harmonized and pooled from 9 NIH-funded longitudinal cohorts participating in the Collaborative Cohort of Cohorts for COVID-19 Research (C4R). CKD was classified by KDIGO 2024 CKD guideline (Low Risk, Moderately Increased Risk, High Risk, Very High Risk) using pre-pandemic urinary albumin-to-creatinine ratio (uACR) and estimated glomerular filtration rate using creatinine (eGFRcr). Incident COVID-19 hospitalization or death was ascertained over 2020−2023 via self-report and active surveillance, with medical record adjudication in a subset. In our primary analysis, we assessed associations between CKD stages and incident COVID-19 hospitalization or death. A secondary analysis examined joint association of CKD stages and the number of comorbidities (obesity, hypertension, diabetes, and cardiovascular disease) with the severe COVID-19. Multivariable hazard ratios (HRs) were estimated by cause-specific Cox regression, adjusting for socio-demographics and smoking. RESULTS: Among 26,039 multi-ethnic adult participants (mean age 66.1±14.8, 60.6 % female, 28.6% Whites, 3.2% Asians, 20.6% Blacks, 41.3% Hispanics, 6.3% American Indians) with a median follow-up of 550 days since the onset of the pandemic (baseline: March 1, 2020), there were 734 incident COVID-19 hospitalizations or deaths. The mean time since CKD markers were measured before the baseline was 6.5 years. The cumulative incidence of severe COVID-19 was 2.3, 3.5, 6.5, and 7.3% across CKD risk stages, and similar patterns were observed within variant types (Figure 1). Compared to participants with Low Risk of CKD, there was a graded increased hazards of severe COVID-19 across CKD stages, with HRs (95% CIs) of 1.27 (1.03−1.55) for Moderately Increased Risk, 2.33 (1.77−3.06) for High Risk, and 2.72 (1.89−3.92) for Very High Risk. In a joint analysis, both increased CKD risk and a higher number of comorbidities were associated with a higher risk of severe COVID-19, showing the highest HR of 12.29 (5.75−26.27) in Very High Risk individuals with four comorbidities compared to Low Risk individuals with none. CONCLUSIONS: Higher CKD stages were associated with a greater risk of COVID-19 hospitalization or death and the risk is substantially heightened with additional comorbidities. These results suggest that adults with CKD and related comorbidities may merit enhanced COVID-19 risk mitigation.

Background: Thrombo-inflammation is a hallmark of acute Covid-19, with higher thrombo-inflammatory markers during infection associated with adverse outcomes. We hypothesized that higher pre-pandemic biomarkers, reflecting a higher thrombo-inflammatory setpoint, might identify people at greater risk of severe manifestations of Covid-19. Hypothesis: Pre-pandemic inflammation markers are associated with risk of hospitalized or fatal (hereafter, “severe”) Covid-19 in a multi-ethnic US sample. Methods: From 2020-2023, C4R harmonized pre-pandemic measures and ascertained Covid-19 across 14 longstanding cohort studies. In 12 cohorts, we assessed biomarker data and covariates from the study visit closest to the pandemic for 6 analytes shown in the Table. HRs of severe Covid-19 by pre-pandemic biomarkers were estimated using cause-specific Cox regression, treating non-severe Covid-19 as a competing risk, and adjusting for pre-pandemic sociodemographic and clinical characteristics. Results: Biomarkers were measured a mean of 20 years before the pandemic with different sample sizes for each (Table footnote). Mean age at the start of the pandemic was in the 60s or 70s depending on the biomarker. The most represented race or ethnic groups were White and Black. The largest sample was for C-reactive protein, and there were 1,001 severe Covid-19 cases among 38,784 participants. In the adjusted models, considered individually, higher CRP, leukocyte count, factor VIII and fibrinogen, but not other markers, were associated with greater risk of severe Covid-19 (Table).Associations were robust across subgroups and time since biomarker measurement (example for C-reactive protein in Figure). Conclusions: In this prospective multi-study analysis of a multi-ethnic population, higher basal levels of thrombo-inflammation markers were associated with greater risk of severe Covid-19. Since we considered non-severe Covid-19 as a competing risk, associations may reflect factors related to progression to more serious Covid-19 illness.

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.

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

Nononcogenic somatic mutations are thought to be uncommon and inconsequential. To test this, we analyzed 43,693 National Heart, Lung and Blood Institute Trans-Omics for Precision Medicine blood whole genomes from 37 cohorts and identified 7131 non-missense somatic mutations that are recurrently mutated in at least 50 individuals. These recurrent non-missense somatic mutations (RNMSMs) are not clearly explained by other clonal phenomena such as clonal hematopoiesis. RNMSM prevalence increased with age, with an average 50-year-old having 27 RNMSMs. Inherited germline variation associated with RNMSM acquisition. These variants were found in genes involved in adaptive immune function, proinflammatory cytokine production, and lymphoid lineage commitment. In addition, the presence of eight specific RNMSMs associated with blood cell traits at effect sizes comparable to Mendelian genetic mutations. Overall, we found that somatic mutations in blood are an unexpectedly common phenomenon with ancestry-specific determinants and human health consequences.

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

Background: Venous thromboembolism (VTE) is a life-threatening vascular event with environmental and genetic determinants. Recent VTE genome-wide association studies (GWAS) meta-analyses involved nearly 30 000 VTE cases and identified up to 40 genetic loci associated with VTE risk, including loci not previously suspected to play a role in hemostasis. The aim of our research was to expand discovery of new genetic loci associated with VTE by using cross-ancestry genomic resources. Methods: We present new cross-ancestry meta-analyzed GWAS results involving up to 81 669 VTE cases from 30 studies, with replication of novel loci in independent populations and loci characterization through in silico genomic interrogations. Results: In our genetic discovery effort that included 55 330 participants with VTE (47 822 European, 6320 African, and 1188 Hispanic ancestry), we identified 48 novel associations, of which 34 were replicated after correction for multiple testing. In our combined discovery-replication analysis (81 669 VTE participants) and ancestry-stratified meta-analyses (European, African, and Hispanic), we identified another 44 novel associations, which are new candidate VTE-associated loci requiring replication. In total, across all GWAS meta-analyses, we identified 135 independent genomic loci significantly associated with VTE risk. A genetic risk score of the significantly associated loci in Europeans identified a 6-fold increase in risk for those in the top 1% of scores compared with those with average scores. We also identified 31 novel transcript associations in transcriptome-wide association studies and 8 novel candidate genes with protein quantitative-trait locus Mendelian randomization analyses. In silico interrogations of hemostasis and hematology traits and a large phenome-wide association analysis of the 135 GWAS loci provided insights to biological pathways contributing to VTE, with some loci contributing to VTE through well-characterized coagulation pathways and others providing new data on the role of hematology traits, particularly platelet function. Many of the replicated loci are outside of known or currently hypothesized pathways to thrombosis. Conclusions: Our cross-ancestry GWAS meta-analyses identified new loci associated with VTE. These findings highlight new pathways to thrombosis and provide novel molecules that may be useful in the development of improved antithrombosis treatments.

Abstract With aging, somatic mutations in hematopoietic stem and progenitor cells (HSPC) can give rise to clonal hematopoiesis of indeterminate potential (CHIP), a premalignant state associated with diverse age-related diseases. Here we report the largest multi-ancestry genome-wide analysis of CHIP to date (N = 1,018,305), including individuals of African (N = 85,978), Admixed American (N = 191,371), East Asian (N = 13,532), European (N = 694,015), and South Asian (N = 13,193) ancestry. Multi-ancestry meta-analyses identified 72 genome-wide significant loci, including 44 novel associations implicating genes such as AFF1 , ATF7IP , ATP8B4 , BCL2 , CEBPA , CYRIA , DNM2 , ELF1 , NKX2-3 , PIK3CB , PRDM16 , RPN1 , TERC , and TRIM4 . Notably, variants at MECOM and PHF20L1 showed opposite allelic effects between DNMT3A - and non- DNMT3A -driven CHIP, highlighting driver-specific germline influences. Implicated loci converge on pathways regulating telomere maintenance, cell-cycle control, hematopoietic transcription, DNA damage response and immune signaling. These findings support a model in which germline variation both expands the HSPC pool and biases clonal selection in a driver-dependent manner, providing a mechanistic basis for inter-individual heterogeneity in CHIP. Phenome-wide association analyses further linked CHIP to hematologic, neoplastic, and circulatory traits, with enrichment across hematopoietic and non-hematopoietic cell types. Together, this work expands the genomic and phenomic landscape of CHIP and reveals germline–somatic interactions that shape clonal evolution during aging.

Human genetic studies support an inverse causal relationship between leukocyte telomere length (LTL) and coronary artery disease (CAD), but directionally mixed effects for LTL and diverse malignancies. Clonal hematopoiesis of indeterminate potential (CHIP), characterized by expansion of hematopoietic cells bearing leukemogenic mutations, predisposes both hematologic malignancy and CAD. TERT (which encodes telomerase reverse transcriptase) is the most significantly associated germline locus for CHIP in genome-wide association studies. Here, we investigated the relationship between CHIP, LTL, and CAD in the Trans-Omics for Precision Medicine (TOPMed) program ( n = 63,302) and UK Biobank ( n = 47,080). Bidirectional Mendelian randomization studies were consistent with longer genetically imputed LTL increasing propensity to develop CHIP, but CHIP then, in turn, hastens to shorten measured LTL (mLTL). We also demonstrated evidence of modest mediation between CHIP and CAD by mLTL. Our data promote an understanding of potential causal relationships across CHIP and LTL toward prevention of CAD.