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.

Admixture mapping of cognitive function in diverse Hispanic and Latino adults: Results from the Hispanic Community Health Study/Study of Latinos.

INTRODUCTION: We conducted admixture mapping and fine-mapping analyses to identify ancestry-of-origin loci influencing cognitive abilities. METHODS: We estimated the association of local ancestry intervals across the genome with five neurocognitive measures in 7140 diverse Hispanic and Latino adults (mean age 55 years). We prioritized genetic variants in associated loci and tested them for replication in four independent cohorts. RESULTS: We identified nine local ancestry-associated regions for the five neurocognitive measures. There was strong biological support for the observed associations to cognitive function at all loci and there was statistical evidence of independent replication at 4q12, 9p22.1, and 13q12.13. DISCUSSION: Our study identified multiple novel loci harboring genes implicated in cognitive functioning and dementia, and uncovered ancestry-relevant genetic variants. It adds to our understanding of the genetic architecture of cognitive function in Hispanic and Latino adults and demonstrates the power of admixture mapping to discover unique haplotypes influencing cognitive function, complementing genome-wide association studies. HIGHLIGHTS: We identified nine ancestry-of-origin chromosomal regions associated with five neurocognitive traits. In each associated region, we identified single nucleotide polymorphisms (SNPs) that explained, at least in part, the admixture signal and were tested for replication in independent samples of Black, non-Hispanic White, and Hispanic/Latino adults with the same or similar neurocognitive tests. Statistical evidence of independent replication of the prioritized SNPs was observed for three of the nine associations, at chr4q12, chr9p22.1, and chr13q12.13. At all loci, there was strong biological support for the observed associations to cognitive function and dementia, prioritizing genes such as KIT, implicated in autophagic clearance of neurotoxic proteins and on mast cell and microglial-mediated inflammation; SLC24A2, implicated in synaptic plasticity associated with learning and memory; and MTMR6, implicated in phosphoinositide lipids metabolism.

Key variants via the Alzheimer's Disease Sequencing Project whole genome sequence data.

INTRODUCTION: Genome-wide association studies (GWAS) have identified loci associated with Alzheimer's disease (AD) but did not identify specific causal genes or variants within those loci. Analysis of whole genome sequence (WGS) data, which interrogates the entire genome and captures rare variations, may identify causal variants within GWAS loci. METHODS: We performed single common variant association analysis and rare variant aggregate analyses in the pooled population (N cases = 2184, N controls = 2383) and targeted analyses in subpopulations using WGS data from the Alzheimer's Disease Sequencing Project (ADSP). The analyses were restricted to variants within 100 kb of 83 previously identified GWAS lead variants. RESULTS: Seventeen variants were significantly associated with AD within five genomic regions implicating the genes OARD1/NFYA/TREML1, JAZF1, FERMT2, and SLC24A4. KAT8 was implicated by both single variant and rare variant aggregate analyses. DISCUSSION: This study demonstrates the utility of leveraging WGS to gain insights into AD loci identified via GWAS.

Unsupervised deep representation learning enables phenotype discovery for genetic association studies of brain imaging.

Understanding the genetic architecture of brain structure is challenging, partly due to difficulties in designing robust, non-biased descriptors of brain morphology. Until recently, brain measures for genome-wide association studies (GWAS) consisted of traditionally expert-defined or software-derived image-derived phenotypes (IDPs) that are often based on theoretical preconceptions or computed from limited amounts of data. Here, we present an approach to derive brain imaging phenotypes using unsupervised deep representation learning. We train a 3-D convolutional autoencoder model with reconstruction loss on 6130 UK Biobank (UKBB) participants' T1 or T2-FLAIR (T2) brain MRIs to create a 128-dimensional representation known as Unsupervised Deep learning derived Imaging Phenotypes (UDIPs). GWAS of these UDIPs in held-out UKBB subjects (n = 22,880 discovery and n = 12,359/11,265 replication cohorts for T1/T2) identified 9457 significant SNPs organized into 97 independent genetic loci of which 60 loci were replicated. Twenty-six loci were not reported in earlier T1 and T2 IDP-based UK Biobank GWAS. We developed a perturbation-based decoder interpretation approach to show that these loci are associated with UDIPs mapped to multiple relevant brain regions. Our results established unsupervised deep learning can derive robust, unbiased, heritable, and interpretable brain imaging phenotypes.

Intracranial Volume Is Driven by Both Genetics and Early Life Exposures: The SOL-INCA-MRI Study.

Intracranial volume (ICV) reflects maximal brain development and is associated with later-life cognitive abilities. We quantified ICV among first- and second-generation Hispanic and Latino adults from the Study of Latinos-Investigation of Cognitive Aging - MRI (SOL-INCA-MRI), estimated ICV heritability, and tested its associations with previously reported genetic variants, both individually and as a genetic risk score (GRS). We also estimated the association of ICV with early life environmental measures: nativity or age of immigration and parental education. The estimated heritability of ICV was 19% (95% CI, 0.1%-56%) in n=1781 unrelated SOL-INCA-MRI individuals. Four of 10 tested genetic variants were associated with ICV and an increase of 1 SD of the ICV-GRS was associated with an increase of 10.37 cm3 in the ICV (95% CI, 5.29-15.45). Compared to being born in the continental United States, immigrating to the United States at age 11 years or older was associated with 24 cm3 smaller ICV (95% CI, -39.97 to -8.06). Compared to both parents having less than high-school education, at least 1 parent completing high-school education was associated with 15.4 cm3 greater ICV (95% CI, 4.46-26.39). These data confirm the importance of early life health on brain development.

White matter hyperintensity genetic risk factor TRIM47 regulates autophagy in brain endothelial cells.

White matter hyperintensity (WMH) is strongly correlated with age-related dementia and hypertension, but its pathogenesis remains obscure. GWAS identified TRIM47 at 17q25 locus as a top genetic risk factor for WMH formation. TRIM family is a class of E3 ubiquitin ligase with pivotal functions in autophagy, which is critical for brain endothelial cell (ECs) remodeling during hypertension. We hypothesize that TRIM47 regulates autophagy and its loss-of-function disturbs cerebrovasculature. Based on transcriptomics and immunohistochemistry, TRIM47 is found selectively expressed by brain ECs in human and mouse, and its transcription is upregulated by artificially-induced autophagy while downregulated in hypertension-like conditions. Using in silico simulation, immunocytochemistry and super-resolution microscopy, we identified the highly conserved binding site between TRIM47 and the LIR (LC3-interacting region) motif of LC3B. Importantly, pharmacological autophagy induction increased Trim47 expression on mouse ECs (b.End3) culture, while silencing Trim47 significantly increased autophagy with ULK1 phosphorylation induction, transcription and vacuole formation. Together, we confirm that TRIM47 is an endogenous inhibitor of autophagy in brain ECs, and such TRIM47-mediated regulation connects genetic and physiological risk factors for WMH formation but warrants further investigation. SUMMARY STATEMENT: TRIM47, top genetic risk factor for white matter hyperintensity formation, is a negative regulator of autophagy in brain endothelial cells and implicates a novel cellular mechanism for age-related cerebrovascular changes.

Genome-Wide Association Study Meta-Analysis Uncovers Novel Genetic Variants Associated with Olfactory Dysfunction.

IMPORTANCE: Olfactory dysfunction is among the earliest signs of many age-related neurodegenerative diseases and has been associated with increased mortality in older adults; however, its genetic basis remains largely unknown. OBJECTIVE: To identify the genetic loci associated with olfactory dysfunction in the general population. DESIGN SETTING AND PARTIICIPANTS: This genome-wide association study meta-analysis (GWMA) included participants of European ancestry (N = 22,730) enrolled in four different large population-based studies, followed by a multi-ancestry GWMA including participants of African ancestry (N = 1,030). The data analysis was performed from March 2023 through June 2024. EXPOSURES: Genome-wide single nucleotide polymorphisms. MAIN OUTCOMES AND MEASURES: Olfactory dysfunction was the outcome and assessed using a 12-item smell identification test. RESULTS: GWMA revealed a novel genome-wide significant locus (tagged by rs11228623 at 11q12) associated with olfactory dysfunction. Gene-based analysis revealed a high enrichment for olfactory receptor genes in this region. Phenome-wide association studies demonstrated associations between genetic variants related to olfactory dysfunction and blood cell counts, kidney function, skeletal muscle mass, cholesterol levels and cardiovascular disease. Using individual-level data, we also confirmed and quantified the strength of these associations on a phenotypic level. Moreover, employing two-sample Mendelian Randomization analyses, we found evidence for causal associations between olfactory dysfunction and these phenotypes. CONCLUSIONS: These findings provide novel insights into the genetic architecture of the sense of smell and highlight its importance for many aspects of human health.

Association of Plasma YKL-40 With MRI, CSF, and Cognitive Markers of Brain Health and Dementia.

BACKGROUND AND OBJECTIVES: Higher YKL-40 levels in the CSF are a known biomarker of brain inflammation. We explored the utility of plasma YKL-40 as a biomarker for accelerated brain aging and dementia risk. METHODS: We performed cross-sectional and prospective analyses of 4 community-based cohorts in the United States or Europe: the Age, Gene/Environment Susceptibility-Reykjavik Study, Atherosclerosis Risk in the Communities study, Coronary Artery Risk Development in Young Adults study, and Framingham Heart Study (FHS). YKL-40 was measured from stored plasma by a single laboratory using Mesoscale Discovery with levels log transformed and standardized within each cohort. Outcomes included MRI total brain volume, hippocampal volume, and white matter hyperintensity volume (WMHV) as a percentage of intracranial volume, a general cognitive composite derived from neuropsychological testing (SD units [SDU]), and the risk of incident dementia. We sought to replicate associations with dementia in the clinic-based ACE csf cohort, which also had YKL-40 measured from the CSF. RESULTS: Meta-analyses of MRI outcomes included 6,558 dementia-free participants, and for analysis of cognition, 6,670. The blood draw preceded MRI/cognitive assessment by up to 10.6 years across cohorts. The mean ages ranged from 50 to 76 years, with 39%-48% male individuals. In random-effects meta-analysis of study estimates, each SDU increase in log-transformed YKL-40 levels was associated with smaller total brain volume (ß = -0.33; 95% CI -0.45 to -0.22; p < 0.0001) and poorer cognition (ß = -0.04; 95% CI -0.07 to -0.02; p < 0.01), following adjustments for demographic variables. YKL-40 levels did not associate with hippocampal volume or WMHV. In the FHS, each SDU increase in log YKL-40 levels was associated with a 64% increase in incident dementia risk over a median of 5.8 years of follow-up, following adjustments for demographic variables (hazard ratio 1.64; 95% CI 1.25-2.16; p < 0.001). In the ACE csf cohort, plasma and CSF YKL-40 were correlated (r = 0.31), and both were associated with conversion from mild cognitive impairment to dementia, independent of amyloid, tau, and neurodegeneration status. DISCUSSION: Higher plasma YKL-40 levels were associated with lower brain volume, poorer cognition, and incident dementia. Plasma YKL-40 may be useful for studying the association of inflammation and its treatment on dementia risk.

Epigenome-wide DNA Methylation Association Study of CHIP Provides Insight into Perturbed Gene Regulation.

With age, hematopoietic stem cells can acquire somatic mutations in leukemogenic genes that confer a proliferative advantage in a phenomenon termed "clonal hematopoiesis of indeterminate potential" (CHIP). How these mutations confer a proliferative advantage and result in increased risk for numerous age-related diseases remains poorly understood. We conducted a multiracial meta-analysis of epigenome-wide association studies (EWAS) of CHIP and its subtypes in four cohorts (N=8196) to elucidate the molecular mechanisms underlying CHIP and illuminate how these changes influence cardiovascular disease risk. The EWAS findings were functionally validated using human hematopoietic stem cell (HSC) models of CHIP. A total of 9615 CpGs were associated with any CHIP, 5990 with DNMT3A CHIP, 5633 with TET2 CHIP, and 6078 with ASXL1 CHIP (P <1×10-7). CpGs associated with CHIP subtypes overlapped moderately, and the genome-wide DNA methylation directions of effect were opposite for TET2 and DNMT3A CHIP, consistent with their opposing effects on global DNA methylation. There was high directional concordance between the CpGs shared from the meta-EWAS and human edited CHIP HSCs. Expression quantitative trait methylation analysis further identified transcriptomic changes associated with CHIP-associated CpGs. Causal inference analyses revealed 261 CHIP-associated CpGs associated with cardiovascular traits and all-cause mortality (FDR adjusted p-value <0.05). Taken together, our study sheds light on the epigenetic changes impacted by CHIP and their associations with age-related disease outcomes. The novel genes and pathways linked to the epigenetic features of CHIP may serve as therapeutic targets for preventing or treating CHIP-mediated diseases.

Familial Hypercholesterolemia Variant and Cardiovascular Risk in Individuals With Elevated Cholesterol.

Importance: Familial hypercholesterolemia (FH) is a genetic disorder that often results in severely high low-density lipoprotein cholesterol (LDL-C) and high risk of premature coronary heart disease (CHD). However, the impact of FH variants on CHD risk among individuals with moderately elevated LDL-C is not well quantified. Objective: To assess CHD risk associated with FH variants among individuals with moderately (130-189 mg/dL) and severely (≥190 mg/dL) elevated LDL-C and to quantify excess CHD deaths attributable to FH variants in US adults. Design, Setting, and Participants: A total of 21 426 individuals without preexisting CHD from 6 US cohort studies (Atherosclerosis Risk in Communities study, Coronary Artery Risk Development in Young Adults study, Cardiovascular Health Study, Framingham Heart Study Offspring cohort, Jackson Heart Study, and Multi-Ethnic Study of Atherosclerosis) were included, 63 of whom had an FH variant. Data were collected from 1971 to 2018, and the median (IQR) follow-up was 18 (13-28) years. Data were analyzed from March to May 2023. Exposures: LDL-C, cumulative past LDL-C, FH variant status. Main Outcomes and Measures: Cox proportional hazards models estimated associations between FH variants and incident CHD. The Cardiovascular Disease Policy Model projected excess CHD deaths associated with FH variants in US adults. Results: Of the 21 426 individuals without preexisting CHD (mean [SD] age 52.1 [15.5] years; 12 041 [56.2%] female), an FH variant was found in 22 individuals with moderately elevated LDL-C (0.3%) and in 33 individuals with severely elevated LDL-C (2.5%). The adjusted hazard ratios for incident CHD comparing those with and without FH variants were 2.9 (95% CI, 1.4-6.0) and 2.6 (95% CI, 1.4-4.9) among individuals with moderately and severely elevated LDL-C, respectively. The association between FH variants and CHD was slightly attenuated when further adjusting for baseline LDL-C level, whereas the association was no longer statistically significant after adjusting for cumulative past LDL-C exposure. Among US adults 20 years and older with no history of CHD and LDL-C 130 mg/dL or higher, more than 417 000 carry an FH variant and were projected to experience more than 12 000 excess CHD deaths in those with moderately elevated LDL-C and 15 000 in those with severely elevated LDL-C compared with individuals without an FH variant. Conclusions and Relevance: In this pooled cohort study, the presence of FH variants was associated with a 2-fold higher CHD risk, even when LDL-C was only moderately elevated. The increased CHD risk appeared to be largely explained by the higher cumulative LDL-C exposure in individuals with an FH variant compared to those without. Further research is needed to assess the value of adding genetic testing to traditional phenotypic FH screening.

Machine learning models for predicting blood pressure phenotypes by combining multiple polygenic risk scores.

We construct non-linear machine learning (ML) prediction models for systolic and diastolic blood pressure (SBP, DBP) using demographic and clinical variables and polygenic risk scores (PRSs). We developed a two-model ensemble, consisting of a baseline model, where prediction is based on demographic and clinical variables only, and a genetic model, where we also include PRSs. We evaluate the use of a linear versus a non-linear model at both the baseline and the genetic model levels and assess the improvement in performance when incorporating multiple PRSs. We report the ensemble model's performance as percentage variance explained (PVE) on a held-out test dataset. A non-linear baseline model improved the PVEs from 28.1 to 30.1% (SBP) and 14.3% to 17.4% (DBP) compared with a linear baseline model. Including seven PRSs in the genetic model computed based on the largest available GWAS of SBP/DBP improved the genetic model PVE from 4.8 to 5.1% (SBP) and 4.7 to 5% (DBP) compared to using a single PRS. Adding additional 14 PRSs computed based on two independent GWASs further increased the genetic model PVE to 6.3% (SBP) and 5.7% (DBP). PVE differed across self-reported race/ethnicity groups, with primarily all non-White groups benefitting from the inclusion of additional PRSs. In summary, non-linear ML models improves BP prediction in models incorporating diverse populations.

Multi-omics and pathway analyses of genome-wide associations implicate regulation and immunity in verbal declarative memory performance.

BACKGROUND: Uncovering the functional relevance underlying verbal declarative memory (VDM) genome-wide association study (GWAS) results may facilitate the development of interventions to reduce age-related memory decline and dementia. METHODS: We performed multi-omics and pathway enrichment analyses of paragraph (PAR-dr) and word list (WL-dr) delayed recall GWAS from 29,076 older non-demented individuals of European descent. We assessed the relationship between single-variant associations and expression quantitative trait loci (eQTLs) in 44 tissues and methylation quantitative trait loci (meQTLs) in the hippocampus. We determined the relationship between gene associations and transcript levels in 53 tissues, annotation as immune genes, and regulation by transcription factors (TFs) and microRNAs. To identify significant pathways, gene set enrichment was tested in each cohort and meta-analyzed across cohorts. Analyses of differential expression in brain tissues were conducted for pathway component genes. RESULTS: The single-variant associations of VDM showed significant linkage disequilibrium (LD) with eQTLs across all tissues and meQTLs within the hippocampus. Stronger WL-dr gene associations correlated with reduced expression in four brain tissues, including the hippocampus. More robust PAR-dr and/or WL-dr gene associations were intricately linked with immunity and were influenced by 31 TFs and 2 microRNAs. Six pathways, including type I diabetes, exhibited significant associations with both PAR-dr and WL-dr. These pathways included fifteen MHC genes intricately linked to VDM performance, showing diverse expression patterns based on cognitive status in brain tissues. CONCLUSIONS: VDM genetic associations influence expression regulation via eQTLs and meQTLs. The involvement of TFs, microRNAs, MHC genes, and immune-related pathways contributes to VDM performance in older individuals.

Determinants of mosaic chromosomal alteration fitness.

Clonal hematopoiesis (CH) is characterized by the acquisition of a somatic mutation in a hematopoietic stem cell that results in a clonal expansion. These driver mutations can be single nucleotide variants in cancer driver genes or larger structural rearrangements called mosaic chromosomal alterations (mCAs). The factors that influence the variations in mCA fitness and ultimately result in different clonal expansion rates are not well understood. We used the Passenger-Approximated Clonal Expansion Rate (PACER) method to estimate clonal expansion rate as PACER scores for 6,381 individuals in the NHLBI TOPMed cohort with gain, loss, and copy-neutral loss of heterozygosity mCAs. Our mCA fitness estimates, derived by aggregating per-individual PACER scores, were correlated (R2 = 0.49) with an alternative approach that estimated fitness of mCAs in the UK Biobank using population-level distributions of clonal fraction. Among individuals with JAK2 V617F clonal hematopoiesis of indeterminate potential or mCAs affecting the JAK2 gene on chromosome 9, PACER score was strongly correlated with erythrocyte count. In a cross-sectional analysis, genome-wide association study of estimates of mCA expansion rate identified a TCL1A locus variant associated with mCA clonal expansion rate, with suggestive variants in NRIP1 and TERT.

Genomic and phenotypic correlates of mosaic loss of chromosome Y in blood.

Mosaic loss of Y (mLOY) is the most common somatic chromosomal alteration detected in human blood. The presence of mLOY is associated with altered blood cell counts and increased risk of Alzheimer's disease, solid tumors, and other age-related diseases. We sought to gain a better understanding of genetic drivers and associated phenotypes of mLOY through analyses of whole genome sequencing of a large set of genetically diverse males from the Trans-Omics for Precision Medicine (TOPMed) program. This approach enabled us to identify differences in mLOY frequencies across populations defined by genetic similarity, revealing a higher frequency of mLOY in the European American (EA) ancestry group compared to those of Hispanic American (HA), African American (AA), and East Asian (EAS) ancestry. Further, we identified two genes ( CFHR1 and LRP6 ) that harbor multiple rare, putatively deleterious variants associated with mLOY susceptibility, show that subsets of human hematopoietic stem cells are enriched for activity of mLOY susceptibility variants, and that certain alleles on chromosome Y are more likely to be lost than others.

Proteogenomics in cerebrospinal fluid and plasma reveals new biological fingerprint of cerebral small vessel disease.

Cerebral small vessel disease (cSVD) is a leading cause of stroke and dementia with no specific mechanism-based treatment. We used Mendelian randomization to combine a unique cerebrospinal fluid (CSF) and plasma pQTL resource with the latest European-ancestry GWAS of MRI-markers of cSVD (white matter hyperintensities, perivascular spaces). We describe a new biological fingerprint of 49 protein-cSVD associations, predominantly in the CSF. We implemented a multipronged follow-up, across fluids, platforms, and ancestries (Europeans and East-Asian), including testing associations of direct plasma protein measurements with MRI-cSVD. We highlight 16 proteins robustly associated in both CSF and plasma, with 24/4 proteins identified in CSF/plasma only. cSVD-proteins were enriched in extracellular matrix and immune response pathways, and in genes enriched in microglia and specific microglial states (integration with single-nucleus RNA sequencing). Immune-related proteins were associated with MRI-cSVD already at age twenty. Half of cSVD-proteins were associated with stroke, dementia, or both, and seven cSVD-proteins are targets for known drugs (used for other indications in directions compatible with beneficial therapeutic effects. This first cSVD proteogenomic signature opens new avenues for biomarker and therapeutic developments.

Machine learning models for blood pressure phenotypes combining multiple polygenic risk scores.

We construct non-linear machine learning (ML) prediction models for systolic and diastolic blood pressure (SBP, DBP) using demographic and clinical variables and polygenic risk scores (PRSs). We developed a two-model ensemble, consisting of a baseline model, where prediction is based on demographic and clinical variables only, and a genetic model, where we also include PRSs. We evaluate the use of a linear versus a non-linear model at both the baseline and the genetic model levels and assess the improvement in performance when incorporating multiple PRSs. We report the ensemble model's performance as percentage variance explained (PVE) on a held-out test dataset. A non-linear baseline model improved the PVEs from 28.1% to 30.1% (SBP) and 14.3% to 17.4% (DBP) compared with a linear baseline model. Including seven PRSs in the genetic model computed based on the largest available GWAS of SBP/DBP improved the genetic model PVE from 4.8% to 5.1% (SBP) and 4.7% to 5% (DBP) compared to using a single PRS. Adding additional 14 PRSs computed based on two independent GWASs further increased the genetic model PVE to 6.3% (SBP) and 5.7% (DBP). PVE differed across self-reported race/ethnicity groups, with primarily all non-White groups benefitting from the inclusion of additional PRSs.