Biomedical subjects

Stephen S Rich

Publications and source records attributed to Stephen S Rich.

34 records · Page 2Linked to original sources

PubMed2025

Alterations in DNA Methylation, Proteomic, and Metabolomic Profiles in African Ancestry Populations with APOL1 Risk Alleles.

KEY POINTS: We aimed to elucidate potential methylation, proteomic, and metabolomic mechanisms by which APOL1 variants may be linked to kidney disease. We report distinct methylation profiling between APOL1 risk allele carriers and noncarriers, many near APOL gene family. We report higher APOL1 protein and lower C18:1 cholesteryl ester in two risk allele carriers. BACKGROUND: The APOL1 high-risk haplotype has been associated with CKD and the deterioration of kidney function, particularly in populations with West African ancestry. However, the mechanisms by which APOL1 risk variants increase the risk for kidney disease and its progression have not been fully elucidated. METHODS: We compared methylation (N=3191; 715 [22%] carriers), proteomic (N=1240; 169 [14%] carriers), and metabolomic (N=6309; 674 [11%] carriers) profiles in African and Hispanic/Latino carriers of two APOL1 high-risk alleles (G1/G1, G2/G2, G1/G2) and noncarriers (G0/G0), excluding heterozygotes (G0/G1, G0/G2), from the Population Architecture using Genomics and Epidemiology Consortium and UK Biobank. In each study, the associations between the APOL1 high-risk haplotype and up to 722,719 cytosine-phosphate-guanine (CpG) sites, 2923 proteins, or 836 metabolites were estimated using covariate-adjusted linear regression models, followed by fixed-effects sample size–weighted meta-analyses. RESULTS: Significant associations were observed between APOL1 high-risk haplotype and methylation at 52 CpG sites, with 48 located on chromosome 22 and 18 in the vicinity of APOL1–4 and MYH9. All significant CpG sites near APOL2 were hypomethylated, whereas those near APOL3 and APOL4 were hypermethylated. APOL1-associated CpG sites were also identified in genes involved in ion transport and mitochondrial stress pathways. Sensitivity analyses indicated consistent yet attenuated effects among heterozygotes, supporting an additive effect of APOL1 risk alleles. Further analyses of the 52 CpG sites identified two near APOL4 exhibiting G1-specific effects, eight associated with CKD but none with eGFR, and three showing heterogeneity by CKD status. In addition, carrying two APOL1 risk alleles was associated with higher plasma APOL1 protein (β=1.12, PFDR = 2.26e-70) and lower C18:1 cholesteryl ester metabolite (Z=−4.50, PFDR = 4.83e-3). CONCLUSIONS: Our results demonstrate differential methylation, proteomic, and metabolomic profiles associated with APOL1 high-risk haplotypes.

APOL1↗

PubMed2025

Rare damaging CCR2 variants are associated with lower lifetime cardiovascular risk.

BACKGROUND: Previous work has shown a role of CCL2, a key chemokine governing monocyte trafficking, in atherosclerosis. However, it remains unknown whether targeting CCR2, the cognate receptor of CCL2, provides protection against human atherosclerotic cardiovascular disease. METHODS: Computationally predicted damaging or loss-of-function (REVEL > 0.5) variants within CCR2 were detected in whole-exome-sequencing data from 454,775 UK Biobank participants and tested for association with cardiovascular endpoints in gene-burden tests. Given the key role of CCR2 in monocyte mobilization, variants associated with lower monocyte count were prioritized for experimental validation. The response to CCL2 of human cells transfected with these variants was tested in migration and cAMP assays. Validated damaging variants were tested for association with cardiovascular endpoints, atherosclerosis burden, and vascular risk factors. Significant associations were replicated in six independent datasets (n = 1,062,595). RESULTS: Carriers of 45 predicted damaging or loss-of-function CCR2 variants (n = 787 individuals) were at lower risk of myocardial infarction and coronary artery disease. One of these variants (M249K, n = 585, 0.15% of European ancestry individuals) was associated with lower monocyte count and with both decreased downstream signaling and chemoattraction in response to CCL2. While M249K showed no association with conventional vascular risk factors, it was consistently associated with a lower risk of myocardial infarction (odds ratio [OR]: 0.66, 95% confidence interval [CI]: 0.54-0.81, p = 6.1 × 10-5) and coronary artery disease (OR: 0.74, 95%CI: 0.63-0.87, p = 2.9 × 10-4) in the UK Biobank and in six replication cohorts. In a phenome-wide association study, there was no evidence of a higher risk of infections among M249K carriers. CONCLUSIONS: Carriers of an experimentally confirmed damaging CCR2 variant are at a lower lifetime risk of myocardial infarction and coronary artery disease without carrying a higher risk of infections. Our findings provide genetic support for the translational potential of CCR2-targeting as an atheroprotective approach.

Humans↗

PubMed2025

Whole genome sequence-based association analysis of African American individuals with bipolar disorder and schizophrenia.

In studies of individuals of primarily European genetic ancestry, common and low-frequency variants and rare coding variants have been found to be associated with the risk of bipolar disorder (BD) and schizophrenia (SZ). However, less is known for individuals of other genetic ancestries or the role of rare non-coding variants in BD and SZ risk. We performed whole genome sequencing of African American individuals: 1,598 with BD, 3,295 with SZ, and 2,651 unaffected controls (InPSYght study). We increased power by incorporating 14,812 jointly called psychiatrically unscreened ancestry-matched controls from the Trans-Omics for Precision Medicine (TOPMed) Program for a total of 17,463 controls. To identify variants and sets of variants associated with BD and/or SZ, we performed single-variant tests, gene-based tests for singleton protein truncating variants, and rare and low-frequency variant annotation-based tests with conservation and universal chromatin states and sliding windows. We found suggestive evidence of BD association with single-variants on chromosome 18 and of lower BD risk associated with rare and low-frequency variants on chromosome 11 in a region with multiple BD GWAS loci, using a sliding window approach. We also found that chromatin and conservation state tests can be used to detect differential calling of variants in controls sequenced at different centers and to assess the effectiveness of sequencing metric covariate adjustments. Our findings reinforce the need for continued whole genome sequencing in additional samples of African American individuals and more comprehensive functional annotation of non-coding variants.

Journal Article↗

PubMed2025

Unveiling the Genetic Landscape of Coronary Artery Disease Through Common and Rare Structural Variants.

BACKGROUND: Genome-wide association studies have identified several hundred susceptibility single nucleotide variants for coronary artery disease (CAD). Despite single nucleotide variant-based genome-wide association studies improving our understanding of the genetics of CAD, the contribution of structural variants (SVs) to the risk of CAD remains largely unclear. METHOD AND RESULTS: We leveraged SVs detected from high-coverage whole genome sequencing data in a diverse group of participants from the National Heart Lung and Blood Institute's Trans-Omics for Precision Medicine program. Single variant tests were performed on 58 706 SVs in a study sample of 11 556 CAD cases and 42 907 controls. Additionally, aggregate tests using sliding windows were performed to examine rare SVs. One genome-wide significant association was identified for a common biallelic intergenic duplication on chromosome 6q21 (P=1.54E-09, odds ratio=1.34). The sliding window-based aggregate tests found 1 region on chromosome 17q25.3, overlapping USP36, to be significantly associated with coronary artery disease (P=1.03E-10). USP36 is highly expressed in arterial and adipose tissues while broadly affecting several cardiometabolic traits. CONCLUSIONS: Our results suggest that SVs, both common and rare, may influence the risk of coronary artery disease.

Humans↗

PubMed2025

Soluble Immune Checkpoint Protein and Lipid Network Associations with All-Cause Mortality Risk: Trans-Omics for Precision Medicine (TOPMed) Program.

Adverse cardiovascular events are emerging with the use of immune checkpoint therapies in oncology. Using datasets in the Trans-Omics for Precision Medicine program (Multi-Ethnic Study of Atherosclerosis, Jackson Heart Study [JHS], and Framingham Heart Study), we examined the association of immune checkpoint plasma proteins with each other, their associated protein network with high-density lipoprotein cholesterol (HDL-C) and low-density lipoprotein cholesterol (LDL-C), and the association of HDL-C- and LDL-C-associated protein networks with all-cause mortality risk. Plasma levels of LAG3 and HAVCR2 showed statistically significant associations with mortality risk. Colocalization analysis using genome wide-association studies of HDL-C or LDL-C and protein quantitative trait loci from JHS and the Atherosclerosis Risk in Communities identified TFF3 rs60467699 and CD36 rs3211938 variants as significantly colocalized with HDL-C; in contrast, none colocalized with LDL-C. The measurement of plasma LAG3, HAVCR2, and associated proteins plus targeted genotyping may identify patients at increased mortality risk.

Journal Article↗

PubMed2025

Associations of High Attenuation Area-Related Proteomic Biomarkers with Fibrotic or Subpleural Interstitial Lung Abnormalities.

Rationale: High-attenuation area (HAA) is a computed tomography (CT) tool that correlates with lung inflammation and fibrosis. Systemic molecular correlates of HAA (e.g., plasma proteins) may inform biological processes involved in interstitial lung disease. Objectives: To identify plasma proteins that associate with HAA and correlate with a higher probability of developing new-onset fibrotic or subpleural interstitial lung abnormalities (ILAs). Methods: Plasma protein levels were measured using a semiquantitative aptamer-based platform in MESA (the Multi-Ethnic Study of Atherosclerosis; N = 5,486) and SPIROMICS (Subpopulations and Intermediate Outcome Measures in COPD Study; N = 1,781). Linear regression models identified HAA-associated proteins after adjustment for demographic and socioeconomic factors, CT scanner parameters, study center, and batch. Associations of HAA-related proteins with new-onset fibrotic or subpleural ILAs were examined in MESA participants with ILA assessments on full-lung CT 10 years later. Immunohistochemical staining of select proteins was performed in lung tissue from pulmonary fibrosis cases. Measurements and Main Results: There were 75 proteins detected that were significantly associated with HAA in MESA and SPIROMICS. Gene Ontology analysis of these proteins identified processes involved in immune cell chemotaxis and cellular growth and apoptosis. Seven proteins were associated with a higher probability of new-onset fibrotic or subpleural ILAs in MESA, and two of these, junctional adhesion molecule-like protein and GTP cyclohydrolase 1 feedback regulatory protein, stained in areas of fibrosis in lung tissue from patients with interstitial lung disease. Conclusions: Plasma proteins associated with more HAA are involved in immune and cellular processes and associate with new-onset fibrotic-subpleural ILA.

Humans↗

PubMed2024

The expected polygenic risk score (ePRS) framework: an equitable metric for quantifying polygenetic risk via modeling of ancestral makeup.

Polygenic risk scores (PRSs) depend on genetic ancestry due to differences in allele frequencies between ancestral populations. This leads to implementation challenges in diverse populations. We propose a framework to calibrate PRS based on ancestral makeup. We define a metric called "expected PRS" (ePRS), the expected value of a PRS based on one's global or local admixture patterns. We further define the "residual PRS" (rPRS), measuring the deviation of the PRS from the ePRS. Simulation studies confirm that it suffices to adjust for ePRS to obtain nearly unbiased estimates of the PRS-outcome association without further adjusting for PCs. Using the TOPMed dataset, the estimated effect size of the rPRS adjusting for the ePRS is similar to the estimated effect of the PRS adjusting for genetic PCs. Similarly, we applied the ePRS framework to six cardiovascular-related traits in the All of Us dataset, and the results are consistent with those from the TOPMed analysis. The ePRS framework can protect from population stratification in association analysis and provide an equitable strategy to quantify genetic risk across diverse populations.

Journal Article↗

PubMed2024

Cardiovascular Risk Factors and Genetic Risk in Transthyretin V142I Carriers.

BACKGROUND: Nearly 3% to 4% of Black individuals in the United States carry the transthyretin V142I variant, which increases their risk of heart failure. However, the role of cardiovascular (CV) risk factors (RFs) in influencing the risk of clinical outcomes among V142I variant carriers is unknown. OBJECTIVES: This study aimed to assess the impact of CV RFs on the risk of heart failure in V142I carriers. METHODS: This study included self-identified Black individuals without prevalent heart failure from 6 TOPMed (Trans-Omics for Precision Medicine) cohorts, the REGARDS (Reasons for Geographic And Racial Differences in Stroke) study, and the All of Us Research Program. The cohort was stratified based on the V142I genotype and the number of CV RFs (hypertension, diabetes, obesity, and hypercholesterolemia). Adjusted Cox models were used to assess the association of heart failure with the V142I genotype and CV RF profile, taking noncarriers with a favorable CV RF profile as reference. RESULTS: The cross-sectional analysis, including 1,625 V142I carriers among 48,365 Black individuals, found that the prevalence of CV RFs did not vary by V142I carrier status. In the longitudinal analysis, there were 587 (3.2%) V142I carriers among 18,407 Black individuals (median age: 60 years [Q1-Q3: 52-68 years], 63.0% female). Among carriers, the heart failure risk was attenuated with a favorable (0 or 1 RF) CV RF profile (adjusted HR: 2.26; 95% CI: 1.58-3.23) compared with an unfavorable (3 or 4 RFs) CV RF profile (adjusted HR: 4.14; 95% CI: 2.79-6.14). CONCLUSIONS: A favorable CV RF profile lowers but does not abrogate V142I variant-associated heart failure risk. This study highlights the importance of having a favorable CV RF profile among V142I carriers for risk reduction of heart failure.

Aged↗

PubMed2024

Genetics of Latin American Diversity Project: Insights into population genetics and association studies in admixed groups in the Americas.

Latin Americans are underrepresented in genetic studies, increasing disparities in personalized genomic medicine. Despite available genetic data from thousands of Latin Americans, accessing and navigating the bureaucratic hurdles for consent or access remains challenging. To address this, we introduce the Genetics of Latin American Diversity (GLAD) Project, compiling genome-wide information from 53,738 Latin Americans across 39 studies representing 46 geographical regions. Through GLAD, we identified heterogeneous ancestry composition and recent gene flow across the Americas. Additionally, we developed GLAD-match, a simulated annealing-based algorithm, to match the genetic background of external samples to our database, sharing summary statistics (i.e., allele and haplotype frequencies) without transferring individual-level genotypes. Finally, we demonstrate the potential of GLAD as a critical resource for evaluating statistical genetic software in the presence of admixture. By providing this resource, we promote genomic research in Latin Americans and contribute to the promises of personalized medicine to more people.

Humans↗

PubMed2024

Rare variant contribution to the heritability of coronary artery disease.

Whole genome sequences (WGS) enable discovery of rare variants which may contribute to missing heritability of coronary artery disease (CAD). To measure their contribution, we apply the GREML-LDMS-I approach to WGS of 4949 cases and 17,494 controls of European ancestry from the NHLBI TOPMed program. We estimate CAD heritability at 34.3% assuming a prevalence of 8.2%. Ultra-rare (minor allele frequency ≤ 0.1%) variants with low linkage disequilibrium (LD) score contribute ~50% of the heritability. We also investigate CAD heritability enrichment using a diverse set of functional annotations: i) constraint; ii) predicted protein-altering impact; iii) cis-regulatory elements from a cell-specific chromatin atlas of the human coronary; and iv) annotation principal components representing a wide range of functional processes. We observe marked enrichment of CAD heritability for most functional annotations. These results reveal the predominant role of ultra-rare variants in low LD on the heritability of CAD. Moreover, they highlight several functional processes including cell type-specific regulatory mechanisms as key drivers of CAD genetic risk.

Humans↗

PubMed2024

Metagenomic Study of the MESA: Detection of Gemella Morbillorum and Association With Coronary Heart Disease.

BACKGROUND: Inflammation is a feature of coronary heart disease (CHD), but the role of proinflammatory microbial infection in CHD remains understudied. METHODS AND RESULTS: CHD was defined in the MESA (Multi-Ethnic Study of Atherosclerosis) as myocardial infarction (251 participants), resuscitated arrest (2 participants), and CHD death (80 participants). We analyzed sequencing reads from 4421 MESA participants in the National Heart, Lung, and Blood Institute Trans-Omics for Precision Medicine program using the PathSeq workflow of the Genome Analysis Tool Kit and a 65-gigabase microbial reference. Paired reads aligning to 840 microbes were detected in >1% of participants. The association of the presence of microbe reads with incident CHD (follow-up, ~18 years) was examined. First, important variables were ascertained using a single regularized Cox proportional hazard model, examining change of risk as a function of presence of microbe with age, sex, education level, Life's Simple 7, and inflammation. For variables of importance, the hazard ratio (HR) was estimated in separate (unregularized) Cox proportional hazard models including the same covariates (significance threshold Bonferroni corrected P<6×10-5, 0.05/840). Reads from 2 microbes were significantly associated with CHD: Gemella morbillorum (HR, 3.14 [95% CI, 1.92-5.12]; P=4.86×10-6) and Pseudomonas species NFACC19-2 (HR, 3.22 [95% CI, 2.03-5.41]; P=1.58×10-6). CONCLUSIONS: Metagenomics of whole-genome sequence reads opens a possible frontier for detection of pathogens for chronic diseases. The association of G morbillorum and Pseudomonas species reads with CHD raises the possibilities that microbes may drive atherosclerotic inflammation and that treatments for specific pathogens may provide clinical utility for CHD reduction.

Humans↗

PubMed2024

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.

Humans↗

PubMed2024

The Genetic Determinants and Genomic Consequences of Non-Leukemogenic Somatic Point Mutations.

Clonal hematopoiesis (CH) is defined by the expansion of a lineage of genetically identical cells in blood. Genetic lesions that confer a fitness advantage, such as point mutations or mosaic chromosomal alterations (mCAs) in genes associated with hematologic malignancy, are frequent mediators of CH. However, recent analyses of both single cell-derived colonies of hematopoietic cells and population sequencing cohorts have revealed CH frequently occurs in the absence of known driver genetic lesions. To characterize CH without known driver genetic lesions, we used 51,399 deeply sequenced whole genomes from the NHLBI TOPMed sequencing initiative to perform simultaneous germline and somatic mutation analyses among individuals without leukemogenic point mutations (LPM), which we term CH-LPMneg. We quantified CH by estimating the total mutation burden. Because estimating somatic mutation burden without a paired-tissue sample is challenging, we developed a novel statistical method, the Genomic and Epigenomic informed Mutation (GEM) rate, that uses external genomic and epigenomic data sources to distinguish artifactual signals from true somatic mutations. We performed a genome-wide association study of GEM to discover the germline determinants of CH-LPMneg. After fine-mapping and variant-to-gene analyses, we identified seven genes associated with CH-LPMneg (TCL1A, TERT, SMC4, NRIP1, PRDM16, MSRA, SCARB1), and one locus associated with a sex-associated mutation pathway (SRGAP2C). We performed a secondary analysis excluding individuals with mCAs, finding that the genetic architecture was largely unaffected by their inclusion. Functional analyses of SMC4 and NRIP1 implicated altered HSC self-renewal and proliferation as the primary mediator of mutation burden in blood. We then performed comprehensive multi-tissue transcriptomic analyses, finding that the expression levels of 404 genes are associated with GEM. Finally, we performed phenotypic association meta-analyses across four cohorts, finding that GEM is associated with increased white blood cell count and increased risk for incident peripheral artery disease, but is not significantly associated with incident stroke or coronary disease events. Overall, we develop GEM for quantifying mutation burden from WGS without a paired-tissue sample and use GEM to discover the genetic, genomic, and phenotypic correlates of CH-LPMneg.

Journal Article↗

PubMed2024

A methylation risk score for chronic kidney disease: a HyperGEN study.

Chronic kidney disease (CKD) impacts about 1 in 7 adults in the United States, but African Americans (AAs) carry a disproportionately higher burden of disease. Epigenetic modifications, such as DNA methylation at cytosine-phosphate-guanine (CpG) sites, have been linked to kidney function and may have clinical utility in predicting the risk of CKD. Given the dynamic relationship between the epigenome, environment, and disease, AAs may be especially sensitive to environment-driven methylation alterations. Moreover, risk models incorporating CpG methylation have been shown to predict disease across multiple racial groups. In this study, we developed a methylation risk score (MRS) for CKD in cohorts of AAs. We selected nine CpG sites that were previously reported to be associated with estimated glomerular filtration rate (eGFR) in epigenome-wide association studies to construct a MRS in the Hypertension Genetic Epidemiology Network (HyperGEN). In logistic mixed models, the MRS was significantly associated with prevalent CKD and was robust to multiple sensitivity analyses, including CKD risk factors. There was modest replication in validation cohorts. In summary, we demonstrated that an eGFR-based CpG score is an independent predictor of prevalent CKD, suggesting that MRS should be further investigated for clinical utility in evaluating CKD risk and progression.

Humans↗

PubMed2023

Genetic variants in the complement system and their potential link in the aetiology of type 1 diabetes.

Type 1 diabetes is an autoimmune disease in which one's own immune system destroys insulin-secreting beta cells in the pancreas. This process results in life-long dependence on exogenous insulin for survival. Both genetic and environmental factors play a role in disease initiation, progression, and ultimate clinical diagnosis of type 1 diabetes. This review will provide background on the natural history of type 1 diabetes and the role of genetic factors involved in the complement system, as several recent studies have identified changes in levels of these proteins as the disease evolves from pre-clinical through to clinically apparent disease.

Humans↗

PubMed2022

Integration of Infant Metabolite, Genetic, and Islet Autoimmunity Signatures to Predict Type 1 Diabetes by Age 6 Years.

CONTEXT: Biomarkers that can accurately predict risk of type 1 diabetes (T1D) in genetically predisposed children can facilitate interventions to delay or prevent the disease. OBJECTIVE: This work aimed to determine if a combination of genetic, immunologic, and metabolic features, measured at infancy, can be used to predict the likelihood that a child will develop T1D by age 6 years. METHODS: Newborns with human leukocyte antigen (HLA) typing were enrolled in the prospective birth cohort of The Environmental Determinants of Diabetes in the Young (TEDDY). TEDDY ascertained children in Finland, Germany, Sweden, and the United States. TEDDY children were either from the general population or from families with T1D with an HLA genotype associated with T1D specific to TEDDY eligibility criteria. From the TEDDY cohort there were 702 children will all data sources measured at ages 3, 6, and 9 months, 11.4% of whom progressed to T1D by age 6 years. The main outcome measure was a diagnosis of T1D as diagnosed by American Diabetes Association criteria. RESULTS: Machine learning-based feature selection yielded classifiers based on disparate demographic, immunologic, genetic, and metabolite features. The accuracy of the model using all available data evaluated by the area under a receiver operating characteristic curve is 0.84. Reducing to only 3- and 9-month measurements did not reduce the area under the curve significantly. Metabolomics had the largest value when evaluating the accuracy at a low false-positive rate. CONCLUSION: The metabolite features identified as important for progression to T1D by age 6 years point to altered sugar metabolism in infancy. Integrating this information with classic risk factors improves prediction of the progression to T1D in early childhood.

Autoantibodies↗