A genome-wide association study provides insights into the genetic etiology of 57 essential and non-essential trace elements in humans.

Trace elements are important for human health but may exert toxic or adverse effects. Mechanisms of uptake, distribution, metabolism, and excretion are partly under genetic control but have not yet been extensively mapped. Here we report a comprehensive multi-element genome-wide association study of 57 essential and non-essential trace elements. We perform genome-wide association meta-analyses of 14 trace elements in up to 6564 Scandinavian whole blood samples, and genome-wide association studies of 43 trace elements in up to 2819 samples measured only in the Trøndelag Health Study (HUNT). We identify 11 novel genetic loci associated with blood concentrations of arsenic, cadmium, manganese, selenium, and zinc in genome-wide association meta-analyses. In HUNT, several genome-wide significant loci are also indicated for other trace elements. Using two-sample Mendelian randomization, we find several indications of weak to moderate effects on health outcomes, the most precise being a weak harmful effect of increased zinc on prostate cancer. However, independent validation is needed. Our current understanding of trace element-associated genetic variants may help establish consequences of trace elements on human health.

Novel ancestry-specific primary open-angle glaucoma loci and shared biology with vascular mechanisms and cell proliferation.

Primary open-angle glaucoma (POAG), a leading cause of irreversible blindness globally, shows disparity in prevalence and manifestations across ancestries. We perform meta-analysis across 15 biobanks (of the Global Biobank Meta-analysis Initiative) (n = 1,487,441: cases = 26,848) and merge with previous multi-ancestry studies, with the combined dataset representing the largest and most diverse POAG study to date (n = 1,478,037: cases = 46,325) and identify 17 novel significant loci, 5 of which were ancestry specific. Gene-enrichment and transcriptome-wide association analyses implicate vascular and cancer genes, a fifth of which are primary ciliary related. We perform an extensive statistical analysis of SIX6 and CDKN2B-AS1 loci in human GTEx data and across large electronic health records showing interaction between SIX6 gene and causal variants in the chr9p21.3 locus, with expression effect on CDKN2A/B. Our results suggest that some POAG risk variants may be ancestry specific, sex specific, or both, and support the contribution of genes involved in programmed cell death in POAG pathogenesis.

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.

Multi-ancestry meta-analysis identifies 5 novel loci for ischemic stroke and reveals heterogeneity of effects between sexes and ancestries.

Stroke is the second leading cause of death and disability worldwide. Stroke prevalence varies by sex and ancestry, possibly due to genetic heterogeneity between subgroups. We performed a genome-wide meta-analysis of 16 biobanks across multiple ancestries to study the genetics of ischemic stroke (60,176 cases, 1,310,725 controls) as part of the Global Biobank Meta-analysis Initiative (GBMI) and further combined the results with previously published MegaStroke. Five novel loci for ischemic stroke (LAMC1, CALCRL, PLSCR1, CDKN1A, and SWAP70) were identified after replication in four additional datasets. One previously reported locus showed significant ancestry heterogeneity (ABO), and one showed significant sex heterogeneity (ALDH2). The ALDH2 association was male specific (males p = 1.67e-24, females p = 0.126) and was additionally observed only in the East Asian ancestry (male) samples. These findings emphasize the need for more diverse datasets with large sample sizes to further understand the genetic predisposition of stroke in different ancestry and sex groups.

Genome-wide association study of thoracic aortic aneurysm and dissection in the Million Veteran Program.

The current understanding of the genetic determinants of thoracic aortic aneurysms and dissections (TAAD) has largely been informed through studies of rare, Mendelian forms of disease. Here, we conducted a genome-wide association study (GWAS) of TAAD, testing ~25 million DNA sequence variants in 8,626 participants with and 453,043 participants without TAAD in the Million Veteran Program, with replication in an independent sample of 4,459 individuals with and 512,463 without TAAD from six cohorts. We identified 21 TAAD risk loci, 17 of which have not been previously reported. We leverage multiple downstream analytic methods to identify causal TAAD risk genes and cell types and provide human genetic evidence that TAAD is a non-atherosclerotic aortic disorder distinct from other forms of vascular disease. Our results demonstrate that the genetic architecture of TAAD mirrors that of other complex traits and that it is not solely inherited through protein-altering variants of large effect size.

A fast linkage method for population GWAS cohorts with related individuals.

Linkage analysis, a class of methods for detecting co-segregation of genomic segments and traits in families, was used to map disease-causing genes for decades before genotyping arrays and dense SNP genotyping enabled genome-wide association studies in population samples. Population samples often contain related individuals, but the segregation of alleles within families is rarely used because traditional linkage methods are computationally inefficient for larger datasets. Here, we describe Population Linkage, a novel application of Haseman-Elston regression as a method of moments estimator of variance components and their standard errors. We achieve additional computational efficiency by using modern methods for detection of IBD segments and variance component estimation, efficient preprocessing of input data, and minimizing redundant numerical calculations. We also refined variance component models to account for the biases in population-scale methods for IBD segment detection. We ran Population Linkage on four blood lipid traits in over 70,000 individuals from the HUNT and SardiNIA studies, successfully detecting 25 known genetic signals. One notable linkage signal that appeared in both was for low-density lipoprotein (LDL) cholesterol levels in the region near the gene APOE (LOD = 29.3, variance explained = 4.1%). This is the region where the missense variants rs7412 and rs429358, which together make up the ε2, ε3, and ε4 alleles each account for 2.4% and 0.8% of variation in circulating LDL cholesterol. Our results show the potential for linkage analysis and other large-scale applications of method of moments variance components estimation.

Global Biobank analyses provide lessons for developing polygenic risk scores across diverse cohorts.

Polygenic risk scores (PRSs) have been widely explored in precision medicine. However, few studies have thoroughly investigated their best practices in global populations across different diseases. We here utilized data from Global Biobank Meta-analysis Initiative (GBMI) to explore methodological considerations and PRS performance in 9 different biobanks for 14 disease endpoints. Specifically, we constructed PRSs using pruning and thresholding (P + T) and PRS-continuous shrinkage (CS). For both methods, using a European-based linkage disequilibrium (LD) reference panel resulted in comparable or higher prediction accuracy compared with several other non-European-based panels. PRS-CS overall outperformed the classic P + T method, especially for endpoints with higher SNP-based heritability. Notably, prediction accuracy is heterogeneous across endpoints, biobanks, and ancestries, especially for asthma, which has known variation in disease prevalence across populations. Overall, we provide lessons for PRS construction, evaluation, and interpretation using GBMI resources and highlight the importance of best practices for PRS in the biobank-scale genomics era.

Association of Rare Protein-Truncating DNA Variants in APOB or PCSK9 With Low-density Lipoprotein Cholesterol Level and Risk of Coronary Heart Disease.

Importance: Protein-truncating variants (PTVs) in apolipoprotein B (APOB) and proprotein convertase subtilisin/kexin type 9 (PCSK9) are associated with significantly lower low-density lipoprotein (LDL) cholesterol concentrations. The association of these PTVs with coronary heart disease (CHD) warrants further characterization in large, multiracial prospective cohort studies. Objective: To evaluate the association of PTVs in APOB and PCSK9 with LDL cholesterol concentrations and CHD risk. Design, Setting, and Participants: This studied included participants from 5 National Heart, Lung, and Blood Institute (NHLBI) studies and the UK Biobank. NHLBI study participants aged 5 to 84 years were recruited between 1971 and 2002 across the US and underwent whole-genome sequencing. UK Biobank participants aged 40 to 69 years were recruited between 2006 and 2010 in the UK and underwent whole-exome sequencing. Data were analyzed from June 2021 to October 2022. Exposures: PTVs in APOB and PCSK9. Main Outcomes and Measures: Estimated untreated LDL cholesterol levels and CHD. Results: Among 19 073 NHLBI participants (10 598 [55.6%] female; mean [SD] age, 52 [17] years), 139 (0.7%) carried an APOB or PCSK9 PTV, which was associated with 49 mg/dL (95% CI, 43-56) lower estimated untreated LDL cholesterol level. Over a median (IQR) follow-up of 21.5 (13.9-29.4) years, incident CHD was observed in 12 of 139 carriers (8.6%) vs 3029 of 18 934 noncarriers (16.0%), corresponding to an adjusted hazard ratio of 0.51 (95% CI, 0.28-0.89; P = .02). Among 190 464 UK Biobank participants (104 831 [55.0%] female; mean [SD] age, 57 [8] years), 662 (0.4%) carried a PTV, which was associated with 45 mg/dL (95% CI, 42-47) lower estimated untreated LDL cholesterol level. Estimated CHD risk by age 75 years was 3.7% (95% CI, 2.0-5.3) in carriers vs 7.0% (95% CI, 6.9-7.2) in noncarriers, corresponding to an adjusted hazard ratio of 0.51 (95% CI, 0.32-0.81; P = .004). Conclusions and Relevance: Among 209 537 individuals in this study, 0.4% carried an APOB or PCSK9 PTV that was associated with less exposure to LDL cholesterol and a 49% lower risk of CHD.

The Michigan Genomics Initiative: A biobank linking genotypes and electronic clinical records in Michigan Medicine patients.

Biobanks of linked clinical patient histories and biological samples are an efficient strategy to generate large cohorts for modern genetics research. Biobank recruitment varies by factors such as geographic catchment and sampling strategy, which affect biobank demographics and research utility. Here, we describe the Michigan Genomics Initiative (MGI), a single-health-system biobank currently consisting of >91,000 participants recruited primarily during surgical encounters at Michigan Medicine. The surgical enrollment results in a biobank enriched for many diseases and ideally suited for a disease genetics cohort. Compared with the much larger population-based UK Biobank, MGI has higher prevalence for nearly all diagnosis-code-based phenotypes and larger absolute case counts for many phenotypes. Genome-wide association study (GWAS) results replicate known findings, thereby validating the genetic and clinical data. Our results illustrate that opportunistic biobank sampling within single health systems provides a unique and complementary resource for exploring the genetics of complex diseases.

Sex-Specific Survival Bias and Interaction Modeling in Coronary Artery Disease Risk Prediction.

BACKGROUND: The 10-year Atherosclerotic Cardiovascular Disease risk score is the standard approach to predict risk of incident cardiovascular events, and recently, addition of coronary artery disease (CAD) polygenic scores has been evaluated. Although age and sex strongly predict the risk of CAD, their interaction with genetic risk prediction has not been systematically examined. This study performed an extensive evaluation of age and sex effects in genetic CAD risk prediction. METHODS: The population-based Norwegian HUNT2 (Trøndelag Health Study 2) cohort of 51 036 individuals was used as the primary dataset. Findings were replicated in the UK Biobank (372 410 individuals). Models for 10-year CAD risk were fitted using Cox proportional hazards, and Harrell concordance index, sensitivity, and specificity were compared. RESULTS: Inclusion of age and sex interactions of CAD polygenic score to the prediction models increased the C-index and sensitivity by accounting for nonadditive effects of CAD polygenic score and likely countering the observed survival bias in the baseline. The sensitivity for females was lower than males in all models including genetic information. We identified a total of 82.6% of incident CAD cases by using a 2-step approach: (1) Atherosclerotic Cardiovascular Disease risk score (74.1%) and (2) the CAD polygenic score interaction model for those in low clinical risk (additional 8.5%). CONCLUSIONS: These findings highlight the importance and complexity of genetic risk in predicting CAD. There is a need for modeling age- and sex-interaction terms with polygenic scores to optimize detection of individuals at high risk, those who warrant preventive interventions. Sex-specific studies are needed to understand and estimate CAD risk with genetic information.

Bicuspid Aortic Valve-Associated Regulatory Regions Reveal GATA4 Regulation and Function During Human-Induced Pluripotent Stem Cell-Based Endothelial-Mesenchymal Transition-Brief Report.

BACKGROUND: The endothelial-mesenchymal transition (EndoMT) is a fundamental process for heart valve formation and defects in EndoMT cause aortic valve abnormalities. Our previous genome-wide association study identified multiple variants in a large chromosome 8 segment as significantly associated with bicuspid aortic valve (BAV). The objective of this study is to determine the biological effects of this large noncoding segment in human induced pluripotent stem cell (hiPSC)-based EndoMT. METHODS: A large genomic segment enriched for BAV-associated variants was deleted in hiPSCs using 2-step CRISPR/Cas9 editing. To address the effects of the variants on GATA4 expression, we generated CRISPR repression hiPSC lines (CRISPRi) as well as hiPSCs from BAV patients. The resulting hiPSCs were differentiated to mesenchymal/myofibroblast-like cells through cardiovascular-lineage endothelial cells for molecular and cellular analysis. Single-cell RNA sequencing was also performed at different stages of EndoMT induction. RESULTS: The large deletion impaired hiPSC-based EndoMT in multiple biallelic clones compared with their isogenic control. It also reduced GATA4 transcript and protein levels during EndoMT, sparing the other genes nearby the deletion segment. Single-cell trajectory analysis revealed the molecular reprogramming during EndoMT. Putative GATA-binding protein targets during EndoMT were uncovered, including genes implicated in endocardial cushion formation and EndoMT process. Differentiation of cells derived from BAV patients carrying the rs117430032 variant as well as CRISPRi repression of the rs117430032 locus resulted in lower GATA4 expression in a stage-specific manner. TWIST1 was identified as a potential regulator of GATA4 expression, showing specificity to the locus tagged by rs117430032. CONCLUSIONS: BAV-associated distal regions regulate GATA4 expression during hiPSC-based EndoMT, which in turn promotes EndoMT progression, implicating its contribution to heart valve development.

Using Genomic Structural Equation Modeling to Partition the Genetic Covariance Between Birthweight and Cardiometabolic Risk Factors into Maternal and Offspring Components in the Norwegian HUNT Study.

The Barker Hypothesis posits that adverse intrauterine environments result in fetal growth restriction and increased risk of cardiometabolic disease through developmental compensations. Here we introduce a new statistical model using the genomic SEM software that is capable of simultaneously partitioning the genetic covariation between birthweight and cardiometabolic traits into maternally mediated and offspring mediated contributions. We model the covariance between birthweight and later life outcomes, such as blood pressure, non-fasting glucose, blood lipids and body mass index in the Norwegian HUNT study, consisting of 15,261 mother-eldest offspring pairs with genetic and phenotypic data. Application of this model showed some evidence for maternally mediated effects of systolic blood pressure on offspring birthweight, and pleiotropy between birthweight and non-fasting glucose mediated through the offspring genome. This underscores the importance of genetic links between birthweight and cardiometabolic phenotypes and offer alternative explanations to environmentally based hypotheses for the phenotypic correlation between these variables.

Increased soluble urokinase plasminogen activator levels modulate monocyte function to promote atherosclerosis.

People with kidney disease are disproportionately affected by atherosclerosis for unclear reasons. Soluble urokinase plasminogen activator receptor (suPAR) is an immune-derived mediator of kidney disease, levels of which are strongly associated with cardiovascular outcomes. We assessed suPAR's pathogenic involvement in atherosclerosis using epidemiologic, genetic, and experimental approaches. We found serum suPAR levels to be predictive of coronary artery calcification and cardiovascular events in 5,406 participants without known coronary disease. In a genome-wide association meta-analysis including over 25,000 individuals, we identified a missense variant in the plasminogen activator, urokinase receptor (PLAUR) gene (rs4760), confirmed experimentally to lead to higher suPAR levels. Mendelian randomization analysis in the UK Biobank using rs4760 indicated a causal association between genetically predicted suPAR levels and atherosclerotic phenotypes. In an experimental model of atherosclerosis, proprotein convertase subtilisin/kexin-9 (Pcsk9) transfection in mice overexpressing suPAR (suPARTg) led to substantially increased atherosclerotic plaques with necrotic cores and macrophage infiltration compared with those in WT mice, despite similar cholesterol levels. Prior to induction of atherosclerosis, aortas of suPARTg mice excreted higher levels of CCL2 and had higher monocyte counts compared with WT aortas. Aortic and circulating suPARTg monocytes exhibited a proinflammatory profile and enhanced chemotaxis. These findings characterize suPAR as a pathogenic factor for atherosclerosis acting at least partially through modulation of monocyte function.

ANGPTL7, a therapeutic target for increased intraocular pressure and glaucoma.

Glaucoma is a leading cause of blindness. Current glaucoma medications work by lowering intraocular pressure (IOP), a risk factor for glaucoma, but most treatments do not directly target the pathological changes leading to increased IOP, which can manifest as medication resistance as disease progresses. To identify physiological modulators of IOP, we performed genome- and exome-wide association analysis in >129,000 individuals with IOP measurements and extended these findings to an analysis of glaucoma risk. We report the identification and functional characterization of rare coding variants (including loss-of-function variants) in ANGPTL7 associated with reduction in IOP and glaucoma protection. We validated the human genetics findings in mice by establishing that Angptl7 knockout mice have lower (~2 mmHg) basal IOP compared to wild-type, with a trend towards lower IOP also in heterozygotes. Conversely, increasing murine Angptl7 levels via injection into mouse eyes increases the IOP. We also show that acute Angptl7 silencing in adult mice lowers the IOP (~2-4 mmHg), reproducing the observations in knockout mice. Collectively, our data suggest that ANGPTL7 is important for IOP homeostasis and is amenable to therapeutic modulation to help maintain a healthy IOP that can prevent onset or slow the progression of glaucoma.

Incorporating family disease history and controlling case-control imbalance for population-based genetic association studies.

MOTIVATION: In the genome-wide association analysis of population-based biobanks, most diseases have low prevalence, which results in low detection power. One approach to tackle the problem is using family disease history, yet existing methods are unable to address type I error inflation induced by increased correlation of phenotypes among closely related samples, as well as unbalanced phenotypic distribution. RESULTS: We propose a new method for genetic association test with family disease history, mixed-model-based Test with Adjusted Phenotype and Empirical saddlepoint approximation, which controls for increased phenotype correlation by adopting a two-variance-component mixed model, accounts for case-control imbalance by using empirical saddlepoint approximation, and is flexible to incorporate any existing adjusted phenotypes, such as phenotypes from the LT-FH method. We show through simulation studies and analysis of UK Biobank data of white British samples and the Korean Genome and Epidemiology Study of Korean samples that the proposed method is robust and yields better calibration compared to existing methods while gaining power for detection of variant-phenotype associations. AVAILABILITY AND IMPLEMENTATION: The summary statistics and code generated in this study are available at https://github.com/styvon/TAPE. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Polygenic Risk Scores for Cardiovascular Disease: A Scientific Statement From the American Heart Association.

Cardiovascular disease is the leading contributor to years lost due to disability or premature death among adults. Current efforts focus on risk prediction and risk factor mitigation' which have been recognized for the past half-century. However, despite advances, risk prediction remains imprecise with persistently high rates of incident cardiovascular disease. Genetic characterization has been proposed as an approach to enable earlier and potentially tailored prevention. Rare mendelian pathogenic variants predisposing to cardiometabolic conditions have long been known to contribute to disease risk in some families. However, twin and familial aggregation studies imply that diverse cardiovascular conditions are heritable in the general population. Significant technological and methodological advances since the Human Genome Project are facilitating population-based comprehensive genetic profiling at decreasing costs. Genome-wide association studies from such endeavors continue to elucidate causal mechanisms for cardiovascular diseases. Systematic cataloging for cardiovascular risk alleles also enabled the development of polygenic risk scores. Genetic profiling is becoming widespread in large-scale research, including in health care-associated biobanks, randomized controlled trials, and direct-to-consumer profiling in tens of millions of people. Thus, individuals and their physicians are increasingly presented with polygenic risk scores for cardiovascular conditions in clinical encounters. In this scientific statement, we review the contemporary science, clinical considerations, and future challenges for polygenic risk scores for cardiovascular diseases. We selected 5 cardiometabolic diseases (coronary artery disease, hypercholesterolemia, type 2 diabetes, atrial fibrillation, and venous thromboembolic disease) and response to drug therapy and offer provisional guidance to health care professionals, researchers, policymakers, and patients.

Genome-wide meta-analysis of iron status biomarkers and the effect of iron on all-cause mortality in HUNT.

Iron is essential for many biological processes, but iron levels must be tightly regulated to avoid harmful effects of both iron deficiency and overload. Here, we perform genome-wide association studies on four iron-related biomarkers (serum iron, serum ferritin, transferrin saturation, total iron-binding capacity) in the Trøndelag Health Study (HUNT), the Michigan Genomics Initiative (MGI), and the SardiNIA study, followed by their meta-analysis with publicly available summary statistics, analyzing up to 257,953 individuals. We identify 123 genetic loci associated with iron traits. Among 19 novel protein-altering variants, we observe a rare missense variant (rs367731784) in HUNT, which suggests a role for DNAJC13 in transferrin recycling. We further validate recently published results using genetic risk scores for each biomarker in HUNT (6% variance in serum iron explained) and present linear and non-linear Mendelian randomization analyses of the traits on all-cause mortality. We find evidence of a harmful effect of increased serum iron and transferrin saturation in linear analyses that estimate population-averaged effects. However, there was weak evidence of a protective effect of increasing serum iron at the very low end of its distribution. Our findings contribute to our understanding of the genes affecting iron status and its consequences on human health.