Integrating Genetic and Transcriptomic Data to Identify Genes Underlying Obesity Risk Loci.

Genome-wide association studies (GWAS) have identified numerous body mass index (BMI) loci. However, most underlying mechanisms from risk locus to BMI remain unknown. Leveraging omics data through integrative analyses could provide more comprehensive views of biological pathways on BMI. We analyzed genotype and blood gene expression data in up to 5,619 samples from the Framingham Heart Study (FHS). Using 3,992 single nucleotide polymorphisms (SNPs) at 97 BMI loci and 20,692 transcripts within 1 Mb, we performed separate association analyses of transcript with BMI and SNP with transcript (PBMI and PSNP, respectively) and then a correlated meta-analysis between the full summary data sets (PMETA). We identified transcripts that met Bonferroni-corrected significance for each omic, were more significant in the correlated meta-analysis than each omic, and were at least nominally associated with BMI in FHS data. Among 308 significant SNP-transcript-BMI associations, we identified seven genes (NT5C2, GSTM3, SNAPC3, SPNS1, TMEM245, YPEL3, and ZNF646) in five association regions. Using an independent sample of blood gene expression data, we validated results for SNAPC3 and YPEL3. We tested for generalization of these associations in hypothalamus, nucleus accumbens, and liver and observed significant (PMETA<0.05 & PMETA

Epigenetic and proteomic signatures associate with clonal hematopoiesis expansion rate.

Clonal hematopoiesis of indeterminate potential (CHIP), whereby somatic mutations in hematopoietic stem cells confer a selective advantage and drive clonal expansion, not only correlates with age but also confers increased risk of morbidity and mortality. Here, we leverage genetically predicted traits to identify factors that determine CHIP clonal expansion rate. We used the passenger-approximated clonal expansion rate method to quantify the clonal expansion rate for 4,370 individuals in the National Heart, Lung, and Blood Institute (NHLBI) Trans-Omics for Precision Medicine (TOPMed) cohort and calculated polygenic risk scores for DNA methylation aging, inflammation-related measures and circulating protein levels. Clonal expansion rate was significantly associated with both genetically predicted and measured epigenetic clocks. No associations were identified with inflammation-related lab values or diseases and CHIP expansion rate overall. A proteome-wide search identified predicted circulating levels of myeloid zinc finger 1 and anti-Müllerian hormone as associated with an increased CHIP clonal expansion rate and tissue inhibitor of metalloproteinase 1 and glycine N-methyltransferase as associated with decreased CHIP clonal expansion rate. Together, our findings identify epigenetic and proteomic patterns associated with the rate of hematopoietic clonal expansion.

Multi-ancestry polygenic risk scores for venous thromboembolism.

Venous thromboembolism (VTE) is a significant contributor to morbidity and mortality, with large disparities in incidence rates between Black and White Americans. Polygenic risk scores (PRSs) limited to variants discovered in genome-wide association studies in European-ancestry samples can identify European-ancestry individuals at high risk of VTE. However, there is limited evidence on whether high-dimensional PRS constructed using more sophisticated methods and more diverse training data can enhance the predictive ability and their utility across diverse populations. We developed PRSs for VTE using summary statistics from the International Network against Venous Thrombosis (INVENT) consortium genome-wide association studies meta-analyses of European- (71 771 cases and 1 059 740 controls) and African-ancestry samples (7482 cases and 129 975 controls). We used LDpred2 and PRS-CSx to construct ancestry-specific and multi-ancestry PRSs and evaluated their performance in an independent European- (6781 cases and 103 016 controls) and African-ancestry sample (1385 cases and 12 569 controls). Multi-ancestry PRSs with weights tuned in European-ancestry samples slightly outperformed ancestry-specific PRSs in European-ancestry test samples (e.g. the area under the receiver operating curve [AUC] was 0.609 for PRS-CSx_combinedEUR and 0.608 for PRS-CSxEUR [P = 0.00029]). Multi-ancestry PRSs with weights tuned in African-ancestry samples also outperformed ancestry-specific PRSs in African-ancestry test samples (PRS-CSxAFR: AUC = 0.58, PRS-CSx_combined AFR: AUC = 0.59), although this difference was not statistically significant (P = 0.34). The highest fifth percentile of the best-performing PRS was associated with 1.9-fold and 1.68-fold increased risk for VTE among European- and African-ancestry subjects, respectively, relative to those in the middle stratum. These findings suggest that the multi-ancestry PRS might be used to improve performance across diverse populations to identify individuals at highest risk for VTE.

Predicting Out-of-Hospital Cardiac Arrest in the General Population Using Electronic Health Records.

BACKGROUND: The majority of out-of-hospital cardiac arrests (OHCAs) occur among individuals in the general population, for whom there is no established strategy to identify risk. In this study, we assess the use of electronic health record (EHR) data to identify OHCA in the general population and define salient factors contributing to OHCA risk. METHODS: The analytical cohort included 2366 individuals with OHCA and 23 660 age- and sex-matched controls receiving health care at the University of Washington. Comorbidities, electrocardiographic measures, vital signs, and medication prescription were abstracted from the EHR. The primary outcome was OHCA. Secondary outcomes included shockable and nonshockable OHCA. Model performance including area under the receiver operating characteristic curve and positive predictive value were assessed and adjusted for observed rate of OHCA across the health system. RESULTS: There were significant differences in demographic characteristics, vital signs, electrocardiographic measures, comorbidities, and medication distribution between individuals with OHCA and controls. In external validation, discrimination in machine learning models (area under the receiver operating characteristic curve 0.80-0.85) was superior to a baseline model with conventional cardiovascular risk factors (area under the receiver operating characteristic curve 0.66). At a specificity threshold of 99%, correcting for baseline OHCA incidence across the health system, positive predictive value was 2.5% to 3.1% in machine learning models compared with 0.8% for the baseline model. Longer corrected QT interval, substance abuse disorder, fluid and electrolyte disorder, alcohol abuse, and higher heart rate were identified as salient predictors of OHCA risk across all machine learning models. Established cardiovascular risk factors retained predictive importance for shockable OHCA, but demographic characteristics (minority race, single marital status) and noncardiovascular comorbidities (substance abuse disorder) also contributed to risk prediction. For nonshockable OHCA, a range of salient predictors, including comorbidities, habits, vital signs, demographic characteristics, and electrocardiographic measures, were identified. CONCLUSIONS: In a population-based case-control study, machine learning models incorporating readily available EHR data showed reasonable discrimination and risk enrichment for OHCA in the general population. Salient factors associated with OCHA risk were myriad across the cardiovascular and noncardiovascular spectrum. Public health and tailored strategies for OHCA prediction and prevention will require incorporation of this complexity.

Genome-wide investigation of exogenous female hormones, genetic variation, and venous thromboembolism risk.

BACKGROUND: Increased risk of venous thromboembolism (VTE) is a life-threatening side effect for users of oral contraceptives (OCs) or hormone therapy (HT). OBJECTIVES: To investigate the potential for genetic predisposition to VTE in OC or HT users, we conducted a gene-by-environment case-only meta-analysis of genome-wide association studies (GWAS). METHODS: Use or nonuse of OCs (7 studies) or HT (8 studies) at the time of the VTE event was determined by pharmacy records or self-report. A synergy index (SI) was modeled for each variant in each study and submultiplicative/supramultiplicative gene-by-environment interactions were estimated. The SI parameters were first meta-analyzed across OC and HT studies and subsequently meta-analyzed to obtain an overall estimate. The primary analysis was agnostic GWAS and interrogated all imputed genotypes using a P value threshold of <5.0 × 10-8; secondary analyses were candidate-based. RESULTS: The VTE case-only OC meta-analysis included 2895 OC users and 6607 nonusers; the case-only HT meta-analysis included 2434 HT users and 12 793 nonusers. In primary GWAS meta-analyses, no variant reached genome-wide significance, but the smallest P value approached statistical significance: rs9386463 (P = 5.03 × 10-8). We tested associations for 138 candidate variants and identified 2 that exceeded statistical significance (0.05/138 = 3.62 × 10-4): F5 rs6025 (P = 1.87 × 10-5; SI, 1.29; previously observed) and F11 rs2036914 (P = 2.0 × 10-4; SI, 0.91; new observation). CONCLUSION: The candidate variant approach to identify submultiplictive/supramultiplicative associations between genetic variation and OC and HT use identified a new association with common genetic variation in F11, while the agnostic interrogations did not yield new discoveries.

Validation of human telomere length multi-ancestry meta-analysis association signals identifies POP5 and KBTBD6 as human telomere length regulation genes.

Genome-wide association studies (GWAS) have become well-powered to detect loci associated with telomere length. However, no prior work has validated genes nominated by GWAS to examine their role in telomere length regulation. We conducted a multi-ancestry meta-analysis of 211,369 individuals and identified five novel association signals. Enrichment analyses of chromatin state and cell-type heritability suggested that blood/immune cells are the most relevant cell type to examine telomere length association signals. We validated specific GWAS associations by overexpressing KBTBD6 or POP5 and demonstrated that both lengthened telomeres. CRISPR/Cas9 deletion of the predicted causal regions in K562 blood cells reduced expression of these genes, demonstrating that these loci are related to transcriptional regulation of KBTBD6 and POP5. Our results demonstrate the utility of telomere length GWAS in the identification of telomere length regulation mechanisms and validate KBTBD6 and POP5 as genes affecting telomere length regulation.

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.

Methylation patterns associated with C-reactive protein in racially and ethnically diverse populations.

Systemic low-grade inflammation is a feature of chronic disease. C-reactive protein (CRP) is a common biomarker of inflammation and used as an indicator of disease risk; however, the role of inflammation in disease is not completely understood. Methylation is an epigenetic modification in the DNA which plays a pivotal role in gene expression. In this study we evaluated differential DNA methylation patterns associated with blood CRP level to elucidate biological pathways and genetic regulatory mechanisms to improve the understanding of chronic inflammation. The racially and ethnically diverse participants in this study were included as 50% White, 41% Black or African American, 7% Hispanic or Latino/a, and 2% Native Hawaiian, Asian American, American Indian, or Alaska Native (total n = 13,433) individuals. We replicated 113 CpG sites from 87 unique loci, of which five were novel (CADM3, NALCN, NLRC5, ZNF792, and cg03282312), across a discovery set of 1,150 CpG sites associated with CRP level (p < 1.2E-7). The downstream pathways affected by DNA methylation included the identification of IFI16 and IRF7 CpG-gene transcript pairs which contributed to the innate immune response gene enrichment pathway along with NLRC5, NOD2, and AIM2. Gene enrichment analysis also identified the nuclear factor-kappaB transcription pathway. Using two-sample Mendelian randomization (MR) we inferred methylation at three CpG sites as causal for CRP levels using both White and Black or African American MR instrument variables. Overall, we identified novel CpG sites and gene transcripts that could be valuable in understanding the specific cellular processes and pathogenic mechanisms involved in inflammation.

Time-to-Event Genome-Wide Association Study for Incident Cardiovascular Disease in People With Type 2 Diabetes.

OBJECTIVE: To identify genetic risk factors for incident cardiovascular disease (CVD) among people with type 2 diabetes (T2D). RESEARCH DESIGN AND METHODS: We conducted a multiancestry time-to-event genome-wide association study for incident CVD among people with T2D. We also tested 204 known coronary artery disease (CAD) variants for association with incident CVD. RESULTS: Among 49,230 participants with T2D, 8,956 had incident CVD events (event rate 18.2%). We identified three novel genetic loci for incident CVD: rs147138607 (near CACNA1E/ZNF648, hazard ratio [HR] 1.23, P = 3.6 × 10-9), rs77142250 (near HS3ST1, HR 1.89, P = 9.9 × 10-9), and rs335407 (near TFB1M/NOX3, HR 1.25, P = 1.5 × 10-8). Among 204 known CAD loci, 5 were associated with incident CVD in T2D (multiple comparison-adjusted P < 0.00024, 0.05/204). A standardized polygenic score of these 204 variants was associated with incident CVD with HR 1.14 (P = 1.0 × 10-16). CONCLUSIONS: The data point to novel and known genomic regions associated with incident CVD among individuals with T2D.

A genetic association study of circulating coagulation factor VIII and von Willebrand factor levels.

ABSTRACT: Coagulation factor VIII (FVIII) and its carrier protein von Willebrand factor (VWF) are critical to coagulation and platelet aggregation. We leveraged whole-genome sequence data from the Trans-Omics for Precision Medicine (TOPMed) program along with TOPMed-based imputation of genotypes in additional samples to identify genetic associations with circulating FVIII and VWF levels in a single-variant meta-analysis, including up to 45 289 participants. Gene-based aggregate tests were implemented in TOPMed. We identified 3 candidate causal genes and tested their functional effect on FVIII release from human liver endothelial cells (HLECs) and VWF release from human umbilical vein endothelial cells. Mendelian randomization was also performed to provide evidence for causal associations of FVIII and VWF with thrombotic outcomes. We identified associations (P < 5 × 10-9) at 7 new loci for FVIII (ST3GAL4, CLEC4M, B3GNT2, ASGR1, F12, KNG1, and TREM1/NCR2) and 1 for VWF (B3GNT2). VWF, ABO, and STAB2 were associated with FVIII and VWF in gene-based analyses. Multiphenotype analysis of FVIII and VWF identified another 3 new loci, including PDIA3. Silencing of B3GNT2 and the previously reported CD36 gene decreased release of FVIII by HLECs, whereas silencing of B3GNT2, CD36, and PDIA3 decreased release of VWF by HVECs. Mendelian randomization supports causal association of higher FVIII and VWF with increased risk of thrombotic outcomes. Seven new loci were identified for FVIII and 1 for VWF, with evidence supporting causal associations of FVIII and VWF with thrombotic outcomes. B3GNT2, CD36, and PDIA3 modulate the release of FVIII and/or VWF in vitro.

Multi-ancestry polygenic risk scores for venous thromboembolism.

Venous thromboembolism (VTE) is a significant contributor to morbidity and mortality, with large disparities in incidence rates between Black and White Americans. Polygenic risk scores (PRSs) limited to variants discovered in genome-wide association studies in European-ancestry samples can identify European-ancestry individuals at high risk of VTE. However, there is limited evidence on whether high-dimensional PRS constructed using more sophisticated methods and more diverse training data can enhance the predictive ability and their utility across diverse populations. We developed PRSs for VTE using summary statistics from the International Network against Venous Thrombosis (INVENT) consortium GWAS meta-analyses of European- (71,771 cases and 1,059,740 controls) and African-ancestry samples (7,482 cases and 129,975 controls). We used LDpred2 and PRSCSx to construct ancestry-specific and multi-ancestry PRSs and evaluated their performance in an independent European- (6,261 cases and 88,238 controls) and African-ancestry sample (1,385 cases and 12,569 controls). Multi-ancestry PRSs with weights tuned in European- and African-ancestry samples, respectively, outperformed ancestry-specific PRSs in European- (PRSCSXEUR: AUC=0.61 (0.60, 0.61), PRSCSX_combinedEUR: AUC=0.61 (0.60, 0.62)) and African-ancestry test samples (PRSCSXAFR: AUC=0.58 (0.57, 0.6), PRSCSX_combined AFR: AUC=0.59 (0.57, 0.60)). The highest fifth percentile of the best-performing PRS was associated with 1.9-fold and 1.68-fold increased risk for VTE among European- and African-ancestry subjects, respectively, relative to those in the middle stratum. These findings suggest that the multi-ancestry PRS may be used to identify individuals at highest risk for VTE and provide guidance for the most effective treatment strategy across diverse populations.

A polygenic risk score of atrial fibrillation improves prediction of lifetime risk for heart failure.

AIMS: Heart failure (HF) has shared genetic architecture with its risk factors: atrial fibrillation (AF), body mass index (BMI), coronary heart disease (CHD), systolic blood pressure (SBP), and type 2 diabetes (T2D). We aim to assess the association and risk prediction performance of risk-factor polygenic risk scores (PRSs) for incident HF and its subtypes in bi-racial populations. METHODS AND RESULTS: Five PRSs were constructed for AF, BMI, CHD, SBP, and T2D in White participants of the Atherosclerosis Risk in Communities (ARIC) study. The associations between PRSs and incident HF and its subtypes were assessed using Cox models, and the risk prediction performance of PRSs was assessed using C statistics. Replication was performed in the ARIC study Black and Cardiovascular Health Study (CHS) White participants. In 8624 ARIC study Whites, 1922 (31% cumulative incidence) HF cases developed over 30 years of follow-up. PRSs of AF, BMI, and CHD were associated with incident HF (P < 0.001), where PRSAF showed the strongest association [hazard ratio (HR): 1.47, 95% confidence interval (CI): 1.41-1.53]. Only the addition of PRSAF to the ARIC study HF risk equation improved C statistics for 10 year risk prediction from 0.812 to 0.829 (∆C: 0.017, 95% CI: 0.009-0.026). The PRSAF was associated with both incident HF with reduced ejection fraction (HR: 1.43, 95% CI: 1.27-1.60) and incident HF with preserved ejection fraction (HR: 1.46, 95% CI: 1.33-1.62). The associations between PRSAF and incident HF and its subtypes, as well as the improved risk prediction, were replicated in the ARIC study Blacks and the CHS Whites (P < 0.050). Protein analyses revealed that N-terminal pro-brain natriuretic peptide and other 98 proteins were associated with PRSAF. CONCLUSIONS: The PRSAF was associated with incident HF and its subtypes and had significant incremental value over an established HF risk prediction equation.

Multi-trait analysis characterizes the genetics of thyroid function and identifies causal associations with clinical implications.

To date only a fraction of the genetic footprint of thyroid function has been clarified. We report a genome-wide association study meta-analysis of thyroid function in up to 271,040 individuals of European ancestry, including reference range thyrotropin (TSH), free thyroxine (FT4), free and total triiodothyronine (T3), proxies for metabolism (T3/FT4 ratio) as well as dichotomized high and low TSH levels. We revealed 259 independent significant associations for TSH (61% novel), 85 for FT4 (67% novel), and 62 novel signals for the T3 related traits. The loci explained 14.1%, 6.0%, 9.5% and 1.1% of the total variation in TSH, FT4, total T3 and free T3 concentrations, respectively. Genetic correlations indicate that TSH associated loci reflect the thyroid function determined by free T3, whereas the FT4 associations represent the thyroid hormone metabolism. Polygenic risk score and Mendelian randomization analyses showed the effects of genetically determined variation in thyroid function on various clinical outcomes, including cardiovascular risk factors and diseases, autoimmune diseases, and cancer. In conclusion, our results improve the understanding of thyroid hormone physiology and highlight the pleiotropic effects of thyroid function on various diseases.

Proteomic prediction of incident heart failure and its main subtypes.

AIM: To examine the ability of serum proteins in predicting future heart failure (HF) events, including HF with reduced or preserved ejection fraction (HFrEF or HFpEF), in relation to event time, and with or without considering established HF-associated clinical variables. METHODS AND RESULTS: In the prospective population-based Age, Gene/Environment Susceptibility Reykjavik Study (AGES-RS), 440 individuals developed HF after their first visit with a median follow-up of 5.45 years. Among them, 167 were diagnosed with HFrEF and 188 with HFpEF. A least absolute shrinkage and selection operator regression model with non-parametric bootstrap were used to select predictors from an analysis of 4782 serum proteins, and several pre-established clinical parameters linked to HF. A subset of 8-10 distinct or overlapping serum proteins predicted different future HF outcomes, and C-statistics were used to assess discrimination, revealing proteins combined with a C-index of 0.80 for all incident HF, 0.78 and 0.80 for incident HFpEF or HFrEF, respectively. In the AGES-RS, protein panels alone encompassed the risk contained in the clinical information and improved the performance characteristics of prediction models based on N-terminal pro-B-type natriuretic peptide and clinical risk factors. Finally, the protein predictors performed particularly well close to the time of an HF event, an outcome that was replicated in the Cardiovascular Health Study. CONCLUSION: A small number of circulating proteins accurately predicted future HF in the AGES-RS cohort of older adults, and they alone encompass the risk information found in a collection of clinical data. Incident HF events were predicted up to 8 years, with predictor performance significantly improving for events occurring less than 1 year ahead, a finding replicated in an external cohort study.

Polygenic Risk Score Associates With Atherosclerotic Plaque Characteristics at Autopsy.

BACKGROUND: Polygenic risk scores (PRSs) for coronary artery disease (CAD) potentially improve cardiovascular risk prediction. However, their relationship with histopathologic features of CAD has never been examined systematically. METHODS: From 4327 subjects referred to CVPath by the State of Maryland Office Chief Medical Examiner for sudden death between 1994 and 2015, 2455 cases were randomly selected for genotyping. We generated PRS from 291 known CAD risk loci. Detailed histopathologic examination of the coronary arteries was performed in all subjects. The primary study outcome measurements were histopathologic plaque features determining severity of atherosclerosis, including %stenosis, calcification, thin-cap fibroatheromas, and thrombotic CAD. RESULTS: After exclusion of cases with insufficient DNA sample quality or with missing data, 954 cases (mean age, 48.8±14.7 years; 75.7% men) remained in the final study cohort. Subjects in the highest PRS quintile exhibited more severe atherosclerosis compared with subjects in the lowest quintile, with greater %stenosis (80.3%±27.0% versus 50.4%±38.7%; adjusted P<0.001) and a higher frequency of calcification (69.6% versus 35.8%; adjusted P=0.004) and thin-cap fibroatheroma (26.7% versus 9.5%; adjusted P=0.007). Even after adjustment for traditional CAD risk factors, subjects within the highest PRS quintile had higher odds of severe atherosclerosis (ie, ≥75% stenosis; adjusted odds ratio, 3.77 [95% CI, 2.10-6.78]; P<0.001) and plaque rupture (adjusted odds ratio, 4.05 [95% CI, 2.26-7.24]; P<0.001). Moreover, subjects within the highest quintile had higher odds of CAD-associated cause of death, especially among those aged ≤50 years (adjusted odds ratio, 4.08 [95% CI, 2.01-8.30]; P<0.001). No statistically significant associations were observed with plaque erosion after adjusting for covariates. CONCLUSIONS: This is the first autopsy study investigating associations between PRS and atherosclerosis severity at the histopathologic level in subjects with sudden death. Our pathological analysis suggests PRS correlates with plaque burden and features of advanced atherosclerosis and may be useful as a method for CAD risk stratification, especially in younger subjects.

A statistical framework for powerful multi-trait rare variant analysis in large-scale whole-genome sequencing studies.

Large-scale whole-genome sequencing (WGS) studies have improved our understanding of the contributions of coding and noncoding rare variants to complex human traits. Leveraging association effect sizes across multiple traits in WGS rare variant association analysis can improve statistical power over single-trait analysis, and also detect pleiotropic genes and regions. Existing multi-trait methods have limited ability to perform rare variant analysis of large-scale WGS data. We propose MultiSTAAR, a statistical framework and computationally-scalable analytical pipeline for functionally-informed multi-trait rare variant analysis in large-scale WGS studies. MultiSTAAR accounts for relatedness, population structure and correlation among phenotypes by jointly analyzing multiple traits, and further empowers rare variant association analysis by incorporating multiple functional annotations. We applied MultiSTAAR to jointly analyze three lipid traits (low-density lipoprotein cholesterol, high-density lipoprotein cholesterol and triglycerides) in 61,861 multi-ethnic samples from the Trans-Omics for Precision Medicine (TOPMed) Program. We discovered new associations with lipid traits missed by single-trait analysis, including rare variants within an enhancer of NIPSNAP3A and an intergenic region on chromosome 1.

Identification of circulating proteins associated with general cognitive function among middle-aged and older adults.

Identifying circulating proteins associated with cognitive function may point to biomarkers and molecular process of cognitive impairment. Few studies have investigated the association between circulating proteins and cognitive function. We identify 246 protein measures quantified by the SomaScan assay as associated with cognitive function (p < 4.9E-5, n up to 7289). Of these, 45 were replicated using SomaScan data, and three were replicated using Olink data at Bonferroni-corrected significance. Enrichment analysis linked the proteins associated with general cognitive function to cell signaling pathways and synapse architecture. Mendelian randomization analysis implicated higher levels of NECTIN2, a protein mediating viral entry into neuronal cells, with higher Alzheimer's disease (AD) risk (p = 2.5E-26). Levels of 14 other protein measures were implicated as consequences of AD susceptibility (p < 2.0E-4). Proteins implicated as causes or consequences of AD susceptibility may provide new insight into the potential relationship between immunity and AD susceptibility as well as potential therapeutic targets.

Type 2 Diabetes Modifies the Association of CAD Genomic Risk Variants With Subclinical Atherosclerosis.

BACKGROUND: Individuals with type 2 diabetes (T2D) have an increased risk of coronary artery disease (CAD), but questions remain about the underlying pathology. Identifying which CAD loci are modified by T2D in the development of subclinical atherosclerosis (coronary artery calcification [CAC], carotid intima-media thickness, or carotid plaque) may improve our understanding of the mechanisms leading to the increased CAD in T2D. METHODS: We compared the common and rare variant associations of known CAD loci from the literature on CAC, carotid intima-media thickness, and carotid plaque in up to 29 670 participants, including up to 24 157 normoglycemic controls and 5513 T2D cases leveraging whole-genome sequencing data from the Trans-Omics for Precision Medicine program. We included first-order T2D interaction terms in each model to determine whether CAD loci were modified by T2D. The genetic main and interaction effects were assessed using a joint test to determine whether a CAD variant, or gene-based rare variant set, was associated with the respective subclinical atherosclerosis measures and then further determined whether these loci had a significant interaction test. RESULTS: Using a Bonferroni-corrected significance threshold of P<1.6×10-4, we identified 3 genes (ATP1B1, ARVCF, and LIPG) associated with CAC and 2 genes (ABCG8 and EIF2B2) associated with carotid intima-media thickness and carotid plaque, respectively, through gene-based rare variant set analysis. Both ATP1B1 and ARVCF also had significantly different associations for CAC in T2D cases versus controls. No significant interaction tests were identified through the candidate single-variant analysis. CONCLUSIONS: These results highlight T2D as an important modifier of rare variant associations in CAD loci with CAC.

The functional impact of rare variation across the regulatory cascade.

Each human genome has tens of thousands of rare genetic variants; however, identifying impactful rare variants remains a major challenge. We demonstrate how use of personal multi-omics can enable identification of impactful rare variants by using the Multi-Ethnic Study of Atherosclerosis, which included several hundred individuals, with whole-genome sequencing, transcriptomes, methylomes, and proteomes collected across two time points, 10 years apart. We evaluated each multi-omics phenotype's ability to separately and jointly inform functional rare variation. By combining expression and protein data, we observed rare stop variants 62 times and rare frameshift variants 216 times as frequently as controls, compared to 13-27 times as frequently for expression or protein effects alone. We extended a Bayesian hierarchical model, "Watershed," to prioritize specific rare variants underlying multi-omics signals across the regulatory cascade. With this approach, we identified rare variants that exhibited large effect sizes on multiple complex traits including height, schizophrenia, and Alzheimer's disease.

Mosaic chromosomal alterations in blood across ancestries using whole-genome sequencing.

Megabase-scale mosaic chromosomal alterations (mCAs) in blood are prognostic markers for a host of human diseases. Here, to gain a better understanding of mCA rates in genetically diverse populations, we analyzed whole-genome sequencing data from 67,390 individuals from the National Heart, Lung, and Blood Institute Trans-Omics for Precision Medicine program. We observed higher sensitivity with whole-genome sequencing data, compared with array-based data, in uncovering mCAs at low mutant cell fractions and found that individuals of European ancestry have the highest rates of autosomal mCAs and the lowest rates of chromosome X mCAs, compared with individuals of African or Hispanic ancestry. Although further studies in diverse populations will be needed to replicate our findings, we report three loci associated with loss of chromosome X, associations between autosomal mCAs and rare variants in DCPS, ADM17, PPP1R16B and TET2 and ancestry-specific variants in ATM and MPL with mCAs in cis.