The power of genetic diversity in genome-wide association studies of lipids.

Increased blood lipid levels are heritable risk factors of cardiovascular disease with varied prevalence worldwide owing to different dietary patterns and medication use1. Despite advances in prevention and treatment, in particular through reducing low-density lipoprotein cholesterol levels2, heart disease remains the leading cause of death worldwide3. Genome-wideassociation studies (GWAS) of blood lipid levels have led to important biological and clinical insights, as well as new drug targets, for cardiovascular disease. However, most previous GWAS4-23 have been conducted in European ancestry populations and may have missed genetic variants that contribute to lipid-level variation in other ancestry groups. These include differences in allele frequencies, effect sizes and linkage-disequilibrium patterns24. Here we conduct a multi-ancestry, genome-wide genetic discovery meta-analysis of lipid levels in approximately 1.65 million individuals, including 350,000 of non-European ancestries. We quantify the gain in studying non-European ancestries and provide evidence to support the expansion of recruitment of additional ancestries, even with relatively small sample sizes. We find that increasing diversity rather than studying additional individuals of European ancestry results in substantial improvements in fine-mapping functional variants and portability of polygenic prediction (evaluated in approximately 295,000 individuals from 7 ancestry groupings). Modest gains in the number of discovered loci and ancestry-specific variants were also achieved. As GWAS expand emphasis beyond the identification of genes and fundamental biology towards the use of genetic variants for preventive and precision medicine25, we anticipate that increased diversity of participants will lead to more accurate and equitable26 application of polygenic scores in clinical practice.

Association of Physical Activity With Bioactive Lipids and Cardiovascular Events.

BACKGROUND: To clarify the mechanisms underlying physical activity (PA)-related cardioprotection, we examined the association of PA with plasma bioactive lipids (BALs) and cardiovascular disease (CVD) events. We additionally performed genome-wide associations. METHODS: PA-bioactive lipid associations were examined in VITAL (VITamin D and OmegA-3 TriaL)-clinical translational science center (REGISTRATION: URL: https://www. CLINICALTRIALS: gov; Unique identifier: NCT01169259; N=1032) and validated in JUPITER (Justification for the Use of statins in Prevention: an Intervention Trial Evaluating Rosuvastatin)-NC (REGISTRATION: URL: https://www. CLINICALTRIALS: gov; Unique identifier: NCT00239681; N=589), using linear models adjusted for age, sex, race, low-density lipoprotein-cholesterol, total-C, and smoking. Significant BALs were carried over to examine associations with incident CVD in 2 nested CVD case-control studies: VITAL-CVD (741 case-control pairs) and JUPITER-CVD (415 case-control pairs; validation). RESULTS: We detected 145 PA-bioactive lipid validated associations (false discovery rate <0.1). Annotations were found for 6 of these BALs: 12,13-diHOME, 9,10-diHOME, lysoPC(15:0), oxymorphone-3b-D-glucuronide, cortisone, and oleoyl-glycerol. Genetic analysis within JUPITER-NC showed associations of 32 PA-related BALs with 22 single-nucleotide polymorphisms. From PA-related BALs, 12 are associated with CVD. CONCLUSIONS: We identified a PA-related bioactive lipidome profile out of which 12 BALs also had opposite associations with incident CVD events.

Cross-Ancestry Investigation of Venous Thromboembolism Genomic Predictors.

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

Pleiotropy-Based Decomposition of Genetic Risk Scores: Association and Interaction Analysis for Type 2 Diabetes and CAD.

Genetic risk for a disease in the population may be represented as a genetic risk score (GRS) constructed as the sum of inherited risk alleles, weighted by allelic effects established in an independent population. While this formulation captures overall genetic risk, it typically does not address risk due to specific biological mechanisms or pathways that may nevertheless be important for interpretation or treatment response. Here, a GRS for disease is resolved into independent or nearly independent components pertaining to biological mechanisms inferred from pleiotropic relationships. The component GRSs' weights are derived from the singular value decomposition (SVD) of the matrix of appropriately scaled genetic effects, i.e., beta coefficients, of the disease variants across a panel of the disease-related phenotypes. The SVD-based formalism also associates combinations of disease-related phenotypes with inferred disease pathways. Applied to incident type 2 diabetes (T2D) in the Women's Genome Health Study (N = 23,294), component GRSs discriminate glycemic control and lipid-based genetic risk, while revealing significant interactions between specific components and BMI or physical activity, the latter not observed with a GRS for overall T2D genetic liability. Applied to coronary artery disease (CAD) in both the WGHS and in JUPITER (N = 8,749), a randomized trial of rosuvastatin for primary prevention of CVD, component GRSs discriminate genetic risk associated with LDL-C from risk associated with reciprocal genetic effects on triglycerides and HDL-C. They also inform the pharmacogenetics of statin treatment by demonstrating that benefit from rosuvastatin is as strongly related to genetic risk from triglycerides and HDL-C as from LDL-C.

Associations of Mitochondrial and Nuclear Mitochondrial Variants and Genes with Seven Metabolic Traits.

Mitochondria (MT), the major site of cellular energy production, are under dual genetic control by 37 mitochondrial DNA (mtDNA) genes and numerous nuclear genes (MT-nDNA). In the CHARGEmtDNA+ Consortium, we studied genetic associations of mtDNA and MT-nDNA associations with body mass index (BMI), waist-hip-ratio (WHR), glucose, insulin, HOMA-B, HOMA-IR, and HbA1c. This 45-cohort collaboration comprised 70,775 (insulin) to 170,202 (BMI) pan-ancestry individuals. Validation and imputation of mtDNA variants was followed by single-variant and gene-based association testing. We report two significant common variants, one in MT-ATP6 associated (p ≤ 5E-04) with WHR and one in the D-loop with glucose. Five rare variants in MT-ATP6, MT-ND5, and MT-ND6 associated with BMI, WHR, or insulin. Gene-based meta-analysis identified MT-ND3 associated with BMI (p ≤ 1E-03). We considered 2,282 MT-nDNA candidate gene associations compiled from online summary results for our traits (20 unique studies with 31 dataset consortia's genome-wide associations [GWASs]). Of these, 109 genes associated (p ≤ 1E-06) with at least 1 of our 7 traits. We assessed regulatory features of variants in the 109 genes, cis- and trans-gene expression regulation, and performed enrichment and protein-protein interactions analyses. Of the identified mtDNA and MT-nDNA genes, 79 associated with adipose measures, 49 with glucose/insulin, 13 with risk for type 2 diabetes, and 18 with cardiovascular disease, indicating for pleiotropic effects with health implications. Additionally, 21 genes related to cholesterol, suggesting additional important roles for the genes identified. Our results suggest that mtDNA and MT-nDNA genes and variants reported make important contributions to glucose and insulin metabolism, adipocyte regulation, diabetes, and cardiovascular disease.

A multi-ancestry genome-wide study incorporating gene-smoking interactions identifies multiple new loci for pulse pressure and mean arterial pressure.

Elevated blood pressure (BP), a leading cause of global morbidity and mortality, is influenced by both genetic and lifestyle factors. Cigarette smoking is one such lifestyle factor. Across five ancestries, we performed a genome-wide gene-smoking interaction study of mean arterial pressure (MAP) and pulse pressure (PP) in 129 913 individuals in stage 1 and follow-up analysis in 480 178 additional individuals in stage 2. We report here 136 loci significantly associated with MAP and/or PP. Of these, 61 were previously published through main-effect analysis of BP traits, 37 were recently reported by us for systolic BP and/or diastolic BP through gene-smoking interaction analysis and 38 were newly identified (P < 5 × 10-8, false discovery rate < 0.05). We also identified nine new signals near known loci. Of the 136 loci, 8 showed significant interaction with smoking status. They include CSMD1 previously reported for insulin resistance and BP in the spontaneously hypertensive rats. Many of the 38 new loci show biologic plausibility for a role in BP regulation. SLC26A7 encodes a chloride/bicarbonate exchanger expressed in the renal outer medullary collecting duct. AVPR1A is widely expressed, including in vascular smooth muscle cells, kidney, myocardium and brain. FHAD1 is a long non-coding RNA overexpressed in heart failure. TMEM51 was associated with contractile function in cardiomyocytes. CASP9 plays a central role in cardiomyocyte apoptosis. Identified only in African ancestry were 30 novel loci. Our findings highlight the value of multi-ancestry investigations, particularly in studies of interaction with lifestyle factors, where genomic and lifestyle differences may contribute to novel findings.

A Large-Scale Multi-ancestry Genome-wide Study Accounting for Smoking Behavior Identifies Multiple Significant Loci for Blood Pressure.

Genome-wide association analysis advanced understanding of blood pressure (BP), a major risk factor for vascular conditions such as coronary heart disease and stroke. Accounting for smoking behavior may help identify BP loci and extend our knowledge of its genetic architecture. We performed genome-wide association meta-analyses of systolic and diastolic BP incorporating gene-smoking interactions in 610,091 individuals. Stage 1 analysis examined ∼18.8 million SNPs and small insertion/deletion variants in 129,913 individuals from four ancestries (European, African, Asian, and Hispanic) with follow-up analysis of promising variants in 480,178 additional individuals from five ancestries. We identified 15 loci that were genome-wide significant (p < 5 × 10-8) in stage 1 and formally replicated in stage 2. A combined stage 1 and 2 meta-analysis identified 66 additional genome-wide significant loci (13, 35, and 18 loci in European, African, and trans-ancestry, respectively). A total of 56 known BP loci were also identified by our results (p < 5 × 10-8). Of the newly identified loci, ten showed significant interaction with smoking status, but none of them were replicated in stage 2. Several loci were identified in African ancestry, highlighting the importance of genetic studies in diverse populations. The identified loci show strong evidence for regulatory features and support shared pathophysiology with cardiometabolic and addiction traits. They also highlight a role in BP regulation for biological candidates such as modulators of vascular structure and function (CDKN1B, BCAR1-CFDP1, PXDN, EEA1), ciliopathies (SDCCAG8, RPGRIP1L), telomere maintenance (TNKS, PINX1, AKTIP), and central dopaminergic signaling (MSRA, EBF2).

Genomic and transcriptomic association studies identify 16 novel susceptibility loci for venous thromboembolism.

Venous thromboembolism (VTE) is a significant contributor to morbidity and mortality. To advance our understanding of the biology contributing to VTE, we conducted a genome-wide association study (GWAS) of VTE and a transcriptome-wide association study (TWAS) based on imputed gene expression from whole blood and liver. We meta-analyzed GWAS data from 18 studies for 30 234 VTE cases and 172 122 controls and assessed the association between 12 923 718 genetic variants and VTE. We generated variant prediction scores of gene expression from whole blood and liver tissue and assessed them for association with VTE. Mendelian randomization analyses were conducted for traits genetically associated with novel VTE loci. We identified 34 independent genetic signals for VTE risk from GWAS meta-analysis, of which 14 are newly reported associations. This included 11 newly associated genetic loci (C1orf198, PLEK, OSMR-AS1, NUGGC/SCARA5, GRK5, MPHOSPH9, ARID4A, PLCG2, SMG6, EIF5A, and STX10) of which 6 replicated, and 3 new independent signals in 3 known genes. Further, TWAS identified 5 additional genetic loci with imputed gene expression levels differing between cases and controls in whole blood (SH2B3, SPSB1, RP11-747H7.3, RP4-737E23.2) and in liver (ERAP1). At some GWAS loci, we found suggestive evidence that the VTE association signal for novel and previously known regions colocalized with expression quantitative trait locus signals. Mendelian randomization analyses suggested that blood traits may contribute to the underlying risk of VTE. To conclude, we identified 16 novel susceptibility loci for VTE; for some loci, the association signals are likely mediated through gene expression of nearby genes.

Metabolic Predictors of Incident Coronary Heart Disease in Women.

BACKGROUND: Although metabolomic profiling offers promise for the prediction of coronary heart disease (CHD), and metabolic risk factors are more strongly associated with CHD in women than men, limited data are available for women. METHODS: We applied a liquid chromatography-tandem mass spectrometry metabolomics platform to measure 371 metabolites in a discovery set of postmenopausal women (472 incident CHD cases, 472 controls) with validation in an independent set of postmenopausal women (312 incident CHD cases, 315 controls). RESULTS: Eight metabolites, primarily oxidized lipids, were significantly dysregulated in cases after the adjustment for matching and CHD risk factors in both the discovery and validation data sets. One oxidized phospholipid, C34:2 hydroxy-phosphatidylcholine, remained associated with CHD after further adjustment for other validated metabolites. Subjects with C34:2 hydroxy-phosphatidylcholine levels in the highest quartile had a 4.7-fold increase in CHD odds in comparison with the lowest quartile; C34:2 hydroxy-phosphatidylcholine also significantly improved the area under the curve (P<0.01) for CHD. The C34:2 hydroxy-phosphatidylcholine findings were replicated in a third replication data set of 980 men and women (230 cardiovascular events) with a stronger association observed in women. CONCLUSIONS: These data replicate known metabolite predictors, identify novel markers, and support the relationship between lipid oxidation and subsequent CHD.

Genetic Obesity and the Risk of Atrial Fibrillation: Causal Estimates from Mendelian Randomization.

BACKGROUND: Observational studies have identified an association between body mass index (BMI) and incident atrial fibrillation (AF). Inferring causality from observational studies, however, is subject to residual confounding, reverse causation, and bias. The primary objective of this study was to evaluate the causal association between BMI and AF by using genetic predictors of BMI. METHODS: We identified 51 646 individuals of European ancestry without AF at baseline from 7 prospective population-based cohorts initiated between 1987 and 2002 in the United States, Iceland, and the Netherlands with incident AF ascertained between 1987 and 2012. Cohort-specific mean follow-up ranged from 7.4 to 19.2 years, over which period there was a total of 4178 cases of incident AF. We performed a Mendelian randomization with instrumental variable analysis to estimate a cohort-specific causal hazard ratio for the association between BMI and AF. Two genetic instruments for BMI were used: FTO genotype (rs1558902) and a BMI gene score comprising 39 single-nucleotide polymorphisms identified by genome-wide association studies to be associated with BMI. Cohort-specific estimates were combined by random-effects, inverse variance-weighted meta-analysis. RESULTS: In age- and sex-adjusted meta-analysis, both genetic instruments were significantly associated with BMI (FTO: 0.43 [95% confidence interval, 0.32-0.54] kg/m2 per A-allele, P<0.001; BMI gene score: 1.05 [95% confidence interval, 0.90-1.20] kg/m2 per 1-U increase, P<0.001) and incident AF (FTO, hazard ratio, 1.07 [1.02-1.11] per A-allele, P=0.004; BMI gene score, hazard ratio, 1.11 [1.05-1.18] per 1-U increase, P<0.001). Age- and sex-adjusted instrumental variable estimates for the causal association between BMI and incident AF were hazard ratio, 1.15 (1.04-1.26) per kg/m2, P=0.005 (FTO) and 1.11 (1.05-1.17) per kg/m2, P<0.001 (BMI gene score). Both of these estimates were consistent with the meta-analyzed estimate between observed BMI and AF (age- and sex-adjusted hazard ratio 1.05 [1.04-1.06] per kg/m2, P<0.001). Multivariable adjustment did not significantly change findings. CONCLUSIONS: Our data are consistent with a causal relationship between BMI and incident AF. These data support the possibility that public health initiatives targeting primordial prevention of obesity may reduce the incidence of AF.

An Empirical Dietary Inflammatory Pattern Score Is Associated with Circulating Inflammatory Biomarkers in a Multi-Ethnic Population of Postmenopausal Women in the United States.

BACKGROUND: The empirical dietary inflammatory pattern (EDIP) score has been associated with concentrations of circulating inflammatory biomarkers in European Americans. OBJECTIVE: We used the EDIP score, a weighted sum of 18 food groups that characterizes dietary inflammatory potential based on circulating concentrations of inflammatory biomarkers, to test the hypothesis that a pro-inflammatory dietary pattern is associated with inflammatory biomarker concentrations in a US multi-ethnic population. METHODS: In this cross-sectional study, we calculated EDIP scores using baseline food frequency questionnaire data from 31,472 women, aged 50-79 y, in the Women's Health Initiative observational study and clinical trials. Circulating biomarkers outcomes at baseline were: C-reactive protein (CRP), interleukin-6 (IL-6), tumor necrosis factor (TNF)-α, TNF receptor (TNFR) 1 and 2, and adiponectin. We used multivariable-adjusted linear regression analyses to estimate absolute concentrations and relative differences in biomarker concentrations, overall and in subgroups of race/ethnicity and BMI (body mass index) categories. RESULTS: Independent of energy intake, BMI, physical activity, and other potential confounding variables, higher EDIP scores were significantly associated with higher (lower for adiponectin) absolute concentrations of all 6 biomarkers. On the relative scale, the percentage of difference in the concentration of biomarkers, among women in the highest compared to the lowest EDIP quintile, was: CRP, +13% (P-trend < 0.0001); IL-6, +15% (P-trend < 0.0001); TNF-α, +7% (P-trend = 0.0007); TNFR1, +4% (P-trend = 0.0009); TNFR2, +5% (P-trend < 0.0001); and adiponectin, -13% (P-trend <0.0001). These associations differed by racial/ethnic groups and by BMI categories. Whereas the absolute biomarker concentrations were lower among European-American women and among normal-weight women, the associations with diet were stronger than among women of African-American or Hispanic/Latino origin and among overweight and obese women. CONCLUSIONS: Findings demonstrate the successful replication of an empirical hypothesis-oriented a posteriori dietary pattern score in a multi-ethnic population of postmenopausal women, with subgroup differences by race/ethnicity and body weight. Future research needs to apply the score in non-US populations.

SOS2 and ACP1 Loci Identified through Large-Scale Exome Chip Analysis Regulate Kidney Development and Function.

Genome-wide association studies have identified >50 common variants associated with kidney function, but these variants do not fully explain the variation in eGFR. We performed a two-stage meta-analysis of associations between genotypes from the Illumina exome array and eGFR on the basis of serum creatinine (eGFRcrea) among participants of European ancestry from the CKDGen Consortium (nStage1: 111,666; nStage2: 48,343). In single-variant analyses, we identified single nucleotide polymorphisms at seven new loci associated with eGFRcrea (PPM1J, EDEM3, ACP1, SPEG, EYA4, CYP1A1, and ATXN2L; PStage1<3.7×10-7), of which most were common and annotated as nonsynonymous variants. Gene-based analysis identified associations of functional rare variants in three genes with eGFRcrea, including a novel association with the SOS Ras/Rho guanine nucleotide exchange factor 2 gene, SOS2 (P=5.4×10-8 by sequence kernel association test). Experimental follow-up in zebrafish embryos revealed changes in glomerular gene expression and renal tubule morphology in the embryonic kidney of acp1- and sos2-knockdowns. These developmental abnormalities associated with altered blood clearance rate and heightened prevalence of edema. This study expands the number of loci associated with kidney function and identifies novel genes with potential roles in kidney formation.

Impact of common genetic determinants of Hemoglobin A1c on type 2 diabetes risk and diagnosis in ancestrally diverse populations: A transethnic genome-wide meta-analysis.

BACKGROUND: Glycated hemoglobin (HbA1c) is used to diagnose type 2 diabetes (T2D) and assess glycemic control in patients with diabetes. Previous genome-wide association studies (GWAS) have identified 18 HbA1c-associated genetic variants. These variants proved to be classifiable by their likely biological action as erythrocytic (also associated with erythrocyte traits) or glycemic (associated with other glucose-related traits). In this study, we tested the hypotheses that, in a very large scale GWAS, we would identify more genetic variants associated with HbA1c and that HbA1c variants implicated in erythrocytic biology would affect the diagnostic accuracy of HbA1c. We therefore expanded the number of HbA1c-associated loci and tested the effect of genetic risk-scores comprised of erythrocytic or glycemic variants on incident diabetes prediction and on prevalent diabetes screening performance. Throughout this multiancestry study, we kept a focus on interancestry differences in HbA1c genetics performance that might influence race-ancestry differences in health outcomes. METHODS & FINDINGS: Using genome-wide association meta-analyses in up to 159,940 individuals from 82 cohorts of European, African, East Asian, and South Asian ancestry, we identified 60 common genetic variants associated with HbA1c. We classified variants as implicated in glycemic, erythrocytic, or unclassified biology and tested whether additive genetic scores of erythrocytic variants (GS-E) or glycemic variants (GS-G) were associated with higher T2D incidence in multiethnic longitudinal cohorts (N = 33,241). Nineteen glycemic and 22 erythrocytic variants were associated with HbA1c at genome-wide significance. GS-G was associated with higher T2D risk (incidence OR = 1.05, 95% CI 1.04-1.06, per HbA1c-raising allele, p = 3 × 10-29); whereas GS-E was not (OR = 1.00, 95% CI 0.99-1.01, p = 0.60). In Europeans and Asians, erythrocytic variants in aggregate had only modest effects on the diagnostic accuracy of HbA1c. Yet, in African Americans, the X-linked G6PD G202A variant (T-allele frequency 11%) was associated with an absolute decrease in HbA1c of 0.81%-units (95% CI 0.66-0.96) per allele in hemizygous men, and 0.68%-units (95% CI 0.38-0.97) in homozygous women. The G6PD variant may cause approximately 2% (N = 0.65 million, 95% CI 0.55-0.74) of African American adults with T2D to remain undiagnosed when screened with HbA1c. Limitations include the smaller sample sizes for non-European ancestries and the inability to classify approximately one-third of the variants. Further studies in large multiethnic cohorts with HbA1c, glycemic, and erythrocytic traits are required to better determine the biological action of the unclassified variants. CONCLUSIONS: As G6PD deficiency can be clinically silent until illness strikes, we recommend investigation of the possible benefits of screening for the G6PD genotype along with using HbA1c to diagnose T2D in populations of African ancestry or groups where G6PD deficiency is common. Screening with direct glucose measurements, or genetically-informed HbA1c diagnostic thresholds in people with G6PD deficiency, may be required to avoid missed or delayed diagnoses.

Rare and low-frequency variants and their association with plasma levels of fibrinogen, FVII, FVIII, and vWF.

Fibrinogen, coagulation factor VII (FVII), and factor VIII (FVIII) and its carrier von Willebrand factor (vWF) play key roles in hemostasis. Previously identified common variants explain only a small fraction of the trait heritabilities, and additional variations may be explained by associations with rarer variants with larger effects. The aim of this study was to identify low-frequency (minor allele frequency [MAF] ≥0.01 and <0.05) and rare (MAF <0.01) variants that influence plasma concentrations of these 4 hemostatic factors by meta-analyzing exome chip data from up to 76,000 participants of 4 ancestries. We identified 12 novel associations of low-frequency (n = 2) and rare (n = 10) variants across the fibrinogen, FVII, FVIII, and vWF traits that were independent of previously identified associations. Novel loci were found within previously reported genes and had effect sizes much larger than and independent of previously identified common variants. In addition, associations at KCNT1, HID1, and KATNB1 identified new candidate genes related to hemostasis for follow-up replication and functional genomic analysis. Newly identified low-frequency and rare-variant associations accounted for modest amounts of trait variance and therefore are unlikely to increase predicted trait heritability but provide new information for understanding individual variation in hemostasis pathways.

Novel locus including FGF21 is associated with dietary macronutrient intake.

Dietary intake of macronutrients (carbohydrate, protein, and fat) has been associated with risk of chronic conditions such as obesity and diabetes. Family studies have reported a moderate contribution of genetics to variation in macronutrient intake. In a genome-wide meta-analysis of a population-based discovery cohort (n = 33 533), rs838133 in FGF21 (19q13.33), rs197273 near TRAF family member-associated NF-kappa-B activator (TANK) (2p24.2), and rs10163409 in FTO (16q12.2) were among the top associations (P < 10(-5)) for percentage of total caloric intake from protein and carbohydrate. rs838133 was replicated in silico in an independent sample from the Cohorts for Heart and Aging Research in Genomic Epidemiology Consortium (CHARGE) Nutrition Working Group (n = 38 360) and attained genome-wide significance in combined analysis (Pjoint = 7.9 × 10(-9)). A cytokine involved in cellular metabolism, FGF21 is a potential susceptibility gene for obesity and type 2 diabetes. Our results highlight the potential of genetic variation for determining dietary macronutrient intake.

Tailoring Risk Prediction Models to Local Populations.

Importance: Risk estimation is an integral part of cardiovascular care. Local recalibration of guideline-recommended models could address the limitations of existing tools. Objective: To provide a machine learning (ML) approach to augment the performance of the American Heart Association's Predicting Risk of Cardiovascular Disease Events (AHA-PREVENT) equations when applied to a local population while preserving clinical interpretability. Design, Setting, and Participants: This cohort study used a New England-based electronic health record cohort of patients without prior atherosclerotic cardiovascular disease (ASCVD) who had the data necessary to calculate the AHA-PREVENT 10-year risk of developing ASCVD in the event period (2007-2016). Patients with prior ASCVD events, death prior to 2007, or age 79 years or older in 2007 were subsequently excluded. The final study population of 95 326 patients was split into 3 nonoverlapping subsets for training, testing, and validation. The AHA-PREVENT model was adapted to this local population using the open-source ML model (MLM) Extreme Gradient Boosting model (XGBoost) with minimal predictor variables, including age, sex, and AHA-PREVENT. The MLM was monotonically constrained to preserve known associations between risk factors and ASCVD risk. Along with sex, race and ethnicity data from the electronic health record were collected to validate the performance of ASCVD risk prediction in subgroups. Data were analyzed from August 2021 to February 2024. Main Outcomes and Measures: Consistent with the AHA-PREVENT model, ASCVD events were defined as the first occurrence of either nonfatal myocardial infarction, coronary artery disease, ischemic stroke, or cardiovascular death. Cardiovascular death was coded via government registries. Discrimination, calibration, and risk reclassification were assessed using the Harrell C index, a modified Hosmer-Lemeshow goodness-of-fit test and calibration curves, and reclassification tables, respectively. Results: In the test set of 38 137 patients (mean [SD] age, 64.8 [6.9] years, 22 708 [59.5]% women and 15 429 [40.5%] men; 935 [2.5%] Asian, 2153 [5.6%] Black, 1414 [3.7%] Hispanic, 31 400 [82.3%] White, and 2235 [5.9%] other, including American Indian, multiple races, unspecified, and unrecorded, consolidated owing to small numbers), MLM-PREVENT had improved calibration (modified Hosmer-Lemeshow P > .05) compared to the AHA-PREVENT model across risk categories in the overall cohort (χ23 = 2.2; P = .53 vs χ23 > 16.3; P < .001) and sex subgroups (men: χ23 = 2.1; P = .55 vs χ23 > 16.3; P < .001; women: χ23 = 6.5; P = .09 vs. χ23 > 16.3; P < .001), while also surpassing a traditional recalibration approach. MLM-PREVENT maintained or improved AHA-PREVENT's calibration in Asian, Black, and White individuals. Both MLM-PREVENT and AHA-PREVENT performed equally well in discriminating risk (approximate ΔC index, ±0.01). Using a clinically significant 7.5% risk threshold, MLM-PREVENT reclassified a total of 11.5% of patients. We visualize the recalibration through MLM-PREVENT ASCVD risk charts that highlight preserved risk associations of the original AHA-PREVENT model. Conclusions and Relevance: The interpretable ML approach presented in this article enhanced the accuracy of the AHA-PREVENT model when applied to a local population while still preserving the risk associations found by the original model. This method has the potential to recalibrate other established risk tools and is implementable in electronic health record systems for improved cardiovascular risk assessment.

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.

X-chromosome and kidney function: evidence from a multi-trait genetic analysis of 908,697 individuals reveals sex-specific and sex-differential findings in genes regulated by androgen response elements.

X-chromosomal genetic variants are understudied but can yield valuable insights into sexually dimorphic human traits and diseases. We performed a sex-stratified cross-ancestry X-chromosome-wide association meta-analysis of seven kidney-related traits (n = 908,697), identifying 23 loci genome-wide significantly associated with two of the traits: 7 for uric acid and 16 for estimated glomerular filtration rate (eGFR), including four novel eGFR loci containing the functionally plausible prioritized genes ACSL4, CLDN2, TSPAN6 and the female-specific DRP2. Further, we identified five novel sex-interactions, comprising male-specific effects at FAM9B and AR/EDA2R, and three sex-differential findings with larger genetic effect sizes in males at DCAF12L1 and MST4 and larger effect sizes in females at HPRT1. All prioritized genes in loci showing significant sex-interactions were located next to androgen response elements (ARE). Five ARE genes showed sex-differential expressions. This study contributes new insights into sex-dimorphisms of kidney traits along with new prioritized gene targets for further molecular research.

A comparison of the theoretical relationship between HDL size and the ratio of HDL cholesterol to apolipoprotein A-I with experimental results from the Women's Health Study.

BACKGROUND: HDL size and composition vary among individuals and may be associated with cardiovascular disease and diabetes. We investigated the theoretical relationship between HDL size and composition using an updated version of the spherical model of lipoprotein structure proposed by Shen et al. (Proc Natl Acad Sci U S A 1977;74:837-41.) and compared its predictions with experimental data from the Women's Health Study (WHS). METHODS: The Shen model was updated to predict the relationship between HDL diameter and the ratio of HDL-cholesterol (HDL-C) to apolipoprotein A-I (ApoA-I) plasma concentrations (HDL-C/ApoA-I ratio). In the WHS (n = 26 772), nuclear magnetic resonance spectroscopy (NMR) was used to measure the mean HDL diameter (d(mean,NMR)) and particle concentration (HDL-P); HDL-C and ApoA-I (mg/dL) were measured by standardized assays. RESULTS: The updated Shen model predicts a quasilinear increase of HDL diameter with the HDL-C/ApoA-I ratio, consistent with the d(mean,NMR) values from WHS, which ranged between 8.0 and 10.8 nm and correlated positively with the HDL-C/ApoA-I ratio (r = 0.608, P < 2.2 × 10(-16)). The WHS data were further described by a linear regression equation: d(WHS) = 4.66 nm + 12.31(HDL-C/Apo-I), where d(WHS) is expressed in nanometers. The validity of this equation for estimating HDL size was assessed with data from cholesteryl ester transfer protein deficiency and pharmacologic inhibition. We also illustrate how HDL-P can be estimated from the HDL size and ApoA-I concentration. CONCLUSIONS: This study provides a large-scale experimental examination of the updated Shen model. The results offer new insights into HDL structure, composition and remodeling and suggest that the HDL-C/ApoA-I ratio might be a readily available biomarker for estimating HDL size and HDL-P.