Search

Background: Genetic factors that influence kidney traits have been understudied for low frequency and ancestry-specific variants. Methods: We combined whole genome sequencing (WGS) data from 23,732 participants from 10 NHLBI Trans-Omics for Precision Medicine (TOPMed) Program multi-ethnic studies to identify novel loci for estimated glomerular filtration rate (eGFR). Participants included European, African, East Asian, and Hispanic ancestries. We applied linear mixed models using a genetic relationship matrix estimated from the WGS data and adjusted for age, sex, study, and ethnicity. Findings: When testing single variants, we identified three novel loci driven by low frequency variants more commonly observed in non-European ancestry (PRKAA2, rs180996919, minor allele frequency [MAF] 0.04%, P = 6.1 × 10−11; METTL8, rs116951054, MAF 0.09%, P = 4.5 × 10−9; and MATK, rs539182790, MAF 0.05%, P = 3.4 × 10−9). We also replicated two known loci for common variants (rs2461702, MAF=0.49, P = 1.2 × 10−9, nearest gene GATM, and rs71147340, MAF=0.34, P = 3.3 × 10−9, CDK12). Testing aggregated variants within a gene identified the MAF gene. A statistical approach based on local ancestry helped to identify replication samples for ancestry-specific variants. Interpretation: This study highlights challenges in studying variants influencing kidney traits that are low frequency in populations and more common in non-European ancestry.

There has been great interest in genetic risk prediction using risk scores in recent years, however, the utility of scores developed in European populations and later applied to non-European populations has not been successful. The goal of this study was to create a methylation risk score (MRS) for metabolic syndrome (MetS), demonstrating the utility of MRS across race groups using cross-sectional data from the Hypertension Genetic Epidemiology Network (HyperGEN, N = 614 African Americans (AA)) and the Genetics of Lipid Lowering Drugs and Diet Network (GOLDN, N = 995 European Americans (EA)). To demonstrate this, we first selected cytosine-guanine dinucleotides (CpG) sites measured on Illumina Methyl450 arrays previously reported to be significantly associated with MetS and/or component conditions in more than one race/ethnic group (CPT1A cg00574958, PHOSPHO1 cg02650017, ABCG1 cg06500161, SREBF1 cg11024682, SOCS3 cg18181703, TXNIP cg19693031). Second, we calculated the parameter estimates for the 6 CpGs in the HyperGEN data (AA) and used the beta estimates as weights to construct a MRS in HyperGEN (AA), which was validated in GOLDN (EA). We performed association analyses using logistic mixed models to test the association between the MRS and MetS, adjusting for covariates. Results showed the MRS was significantly associated with MetS in both populations. In summary, a MRS for MetS was a strong predictor for the condition across two race groups, suggesting MRS may be useful to examine metabolic disease risk or related complications across race/ethnic groups.

Handgrip strength is a widely used measure of muscle strength and a predictor of a range of morbidities including cardiovascular diseases and all-cause mortality. Previous genome-wide association studies of handgrip strength have focused on common variants primarily in persons of European descent. We aimed to identify rare and ancestry-specific genetic variants associated with handgrip strength by conducting whole-genome sequence association analyses using 13,552 participants from six studies representing diverse population groups from the Trans-Omics in Precision Medicine (TOPMed) Program. By leveraging multiple handgrip strength measures performed in study participants over time, we increased our effective sample size by 7-12%. Single-variant analyses identified ten handgrip strength loci among African-Americans: four rare variants, five low-frequency variants, and one common variant. One significant and four suggestive genes were identified associated with handgrip strength when aggregating rare and functional variants; all associations were ancestry-specific. We additionally leveraged the different ancestries available in the UK Biobank to further explore the ancestry-specific association signals from the single-variant association analyses. In conclusion, our study identified 11 new loci associated with handgrip strength with rare and/or ancestry-specific genetic variations, highlighting the added value of whole-genome sequencing in diverse samples. Several of the associations identified using single-variant or aggregate analyses lie in genes with a function relevant to the brain or muscle or were reported to be associated with muscle or age-related traits. Further studies in samples with sequence data and diverse ancestries are needed to confirm these findings.

Smoking is a potentially causal behavioral risk factor for type 2 diabetes (T2D), but not all smokers develop T2D. It is unknown whether genetic factors partially explain this variation. We performed genome-environment-wide interaction studies to identify loci exhibiting potential interaction with baseline smoking status (ever vs. never) on incident T2D and fasting glucose (FG). Analyses were performed in participants of European (EA) and African ancestry (AA) separately. Discovery analyses were conducted using genotype data from the 50,000-single-nucleotide polymorphism (SNP) ITMAT-Broad-CARe (IBC) array in 5 cohorts from from the Candidate Gene Association Resource Consortium (n = 23,189). Replication was performed in up to 16 studies from the Cohorts for Heart Aging Research in Genomic Epidemiology Consortium (n = 74,584). In meta-analysis of discovery and replication estimates, 5 SNPs met at least one criterion for potential interaction with smoking on incident T2D at p

Most genome-wide association and fine-mapping studies to date have been conducted in individuals of European descent, and genetic studies of populations of Hispanic/Latino and African ancestry are limited. In addition, these populations have more complex linkage disequilibrium structure. In order to better define the genetic architecture of these understudied populations, we leveraged >100,000 phased sequences available from deep-coverage whole genome sequencing through the multi-ethnic NHLBI Trans-Omics for Precision Medicine (TOPMed) program to impute genotypes into admixed African and Hispanic/Latino samples with genome-wide genotyping array data. We demonstrated that using TOPMed sequencing data as the imputation reference panel improves genotype imputation quality in these populations, which subsequently enhanced gene-mapping power for complex traits. For rare variants with minor allele frequency (MAF) < 0.5%, we observed a 2.3- to 6.1-fold increase in the number of well-imputed variants, with 11-34% improvement in average imputation quality, compared to the state-of-the-art 1000 Genomes Project Phase 3 and Haplotype Reference Consortium reference panels. Impressively, even for extremely rare variants with minor allele count 86%. Subsequent association analyses of TOPMed reference panel-imputed genotype data with hematological traits (hemoglobin (HGB), hematocrit (HCT), and white blood cell count (WBC)) in ~21,600 African-ancestry and ~21,700 Hispanic/Latino individuals identified associations with two rare variants in the HBB gene (rs33930165 with higher WBC [p = 8.8x10-15] in African populations, rs11549407 with lower HGB [p = 1.5x10-12] and HCT [p = 8.8x10-10] in Hispanics/Latinos). By comparison, neither variant would have been genome-wide significant if either 1000 Genomes Project Phase 3 or Haplotype Reference Consortium reference panels had been used for imputation. Our findings highlight the utility of the TOPMed imputation reference panel for identification of novel rare variant associations not previously detected in similarly sized genome-wide studies of under-represented African and Hispanic/Latino populations.

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

Abstract The role of rare non-coding variation in complex human phenotypes is still largely unknown. To elucidate the impact of rare variants in regulatory elements, we performed a whole-genome sequencing association analysis for height using 333,100 individuals from three datasets: UK Biobank (N = 200,003), TOPMed (N = 87,652) and All of Us (N = 45,445). We performed rare (

[This corrects the article DOI: 10.1371/journal.pgen.1006728.].

BACKGROUND:Odd-numbered chain saturated fatty acids (OCSFA) have been associated with potential health benefits. Although some OCSFA (e.g., C15:0 and C17:0) are found in meats and dairy products, sources and metabolism of C19:0 and C23:0 are relatively unknown, and the influence of non-dietary determinants, including genetic factors, on circulating levels of OCSFA is not established. OBJECTIVE:To elucidate the biological processes that influence circulating levels of OCSFA by investigating associations between genetic variation and OCSFA. DESIGN:We performed a meta-analysis of genome-wide association studies (GWAS) of plasma phospholipid/erythrocyte levels of C15:0, C17:0, C19:0, and C23:0 among 11,494 individuals of European descent. We also investigated relationships between specific single nucleotide polymorphisms (SNPs) in the lactase (LCT) gene, associated with adult-onset lactase intolerance, with circulating levels of dairy-derived OCSFA, and evaluated associations of candidate sphingolipid genes with C23:0 levels. RESULTS:We found no genome-wide significant evidence that common genetic variation is associated with circulating levels of C15:0 or C23:0. In two cohorts with available data, we identified one intronic SNP (rs13361131) in myosin X gene (MYO10) associated with C17:0 level (P = 1.37×10-8), and two intronic SNP (rs12874278 and rs17363566) in deleted in lymphocytic leukemia 1 (DLEU1) region associated with C19:0 level (P = 7.07×10-9). In contrast, when using a candidate-gene approach, we found evidence that three SNPs in LCT (rs11884924, rs16832067, and rs3816088) are associated with circulating C17:0 level (adjusted P = 4×10-2). In addition, nine SNPs in the ceramide synthase 4 (CERS4) region were associated with circulating C23:0 levels (adjusted P

Hypertension is a leading cause of global disease, mortality, and disability. While individuals of African descent suffer a disproportionate burden of hypertension and its complications, they have been underrepresented in genetic studies. To identify novel susceptibility loci for blood pressure and hypertension in people of African ancestry, we performed both single and multiple-trait genome-wide association analyses. We analyzed 21 genome-wide association studies comprised of 31,968 individuals of African ancestry, and validated our results with additional 54,395 individuals from multi-ethnic studies. These analyses identified nine loci with eleven independent variants which reached genome-wide significance (P < 1.25×10-8) for either systolic and diastolic blood pressure, hypertension, or for combined traits. Single-trait analyses identified two loci (TARID/TCF21 and LLPH/TMBIM4) and multiple-trait analyses identified one novel locus (FRMD3) for blood pressure. At these three loci, as well as at GRP20/CDH17, associated variants had alleles common only in African-ancestry populations. Functional annotation showed enrichment for genes expressed in immune and kidney cells, as well as in heart and vascular cells/tissues. Experiments driven by these findings and using angiotensin-II induced hypertension in mice showed altered kidney mRNA expression of six genes, suggesting their potential role in hypertension. Our study provides new evidence for genes related to hypertension susceptibility, and the need to study African-ancestry populations in order to identify biologic factors contributing to hypertension.

Genome-wide association studies (GWAS) have identified >300 loci associated with measures of adiposity including body mass index (BMI) and waist-to-hip ratio (adjusted for BMI, WHRadjBMI), but few have been identified through screening of the African ancestry genomes. We performed large scale meta-analyses and replications in up to 52,895 individuals for BMI and up to 23,095 individuals for WHRadjBMI from the African Ancestry Anthropometry Genetics Consortium (AAAGC) using 1000 Genomes phase 1 imputed GWAS to improve coverage of both common and low frequency variants in the low linkage disequilibrium African ancestry genomes. In the sex-combined analyses, we identified one novel locus (TCF7L2/HABP2) for WHRadjBMI and eight previously established loci at P < 5×10-8: seven for BMI, and one for WHRadjBMI in African ancestry individuals. An additional novel locus (SPRYD7/DLEU2) was identified for WHRadjBMI when combined with European GWAS. In the sex-stratified analyses, we identified three novel loci for BMI (INTS10/LPL and MLC1 in men, IRX4/IRX2 in women) and four for WHRadjBMI (SSX2IP, CASC8, PDE3B and ZDHHC1/HSD11B2 in women) in individuals of African ancestry or both African and European ancestry. For four of the novel variants, the minor allele frequency was low (

BackgroundThe link between DNA methylation, obesity, and adiposity-related diseases in the general population remains uncertain.Methods and findingsWe conducted an association study of body mass index (BMI) and differential methylation for over 400,000 CpGs assayed by microarray in whole-blood-derived DNA from 3,743 participants in the Framingham Heart Study and the Lothian Birth Cohorts, with independent replication in three external cohorts of 4,055 participants. We examined variations in whole blood gene expression and conducted Mendelian randomization analyses to investigate the functional and clinical relevance of the findings. We identified novel and previously reported BMI-related differential methylation at 83 CpGs that replicated across cohorts; BMI-related differential methylation was associated with concurrent changes in the expression of genes in lipid metabolism pathways. Genetic instrumental variable analysis of alterations in methylation at one of the 83 replicated CpGs, cg11024682 (intronic to sterol regulatory element binding transcription factor 1 [SREBF1]), demonstrated links to BMI, adiposity-related traits, and coronary artery disease. Independent genetic instruments for expression of SREBF1 supported the findings linking methylation to adiposity and cardiometabolic disease. Methylation at a substantial proportion (16 of 83) of the identified loci was found to be secondary to differences in BMI. However, the cross-sectional nature of the data limits definitive causal determination.ConclusionsWe present robust associations of BMI with differential DNA methylation at numerous loci in blood cells. BMI-related DNA methylation and gene expression provide mechanistic insights into the relationship between DNA methylation, obesity, and adiposity-related diseases.

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

In this study, we conducted an epigenome-wide association study of metabolic syndrome (MetS) among 846 participants of European descent in the Genetics of Lipid Lowering Drugs and Diet Network (GOLDN). DNA was isolated from CD4+ T cells and methylation at ~470,000 cytosine-phosphate-guanine dinucleotide (CpG) pairs was assayed using the Illumina Infinium HumanMethylation450 BeadChip. We modeled the percentage methylation at individual CpGs as a function of MetS using linear mixed models. A Bonferroni-corrected P-value of 1.1 x 10(-7) was considered significant. Methylation at two CpG sites in CPT1A on chromosome 11 was significantly associated with MetS (P for cg00574958 = 2.6x10(-14) and P for cg17058475 = 1.2x10(-9)). Significant associations were replicated in both European and African ancestry participants of the Bogalusa Heart Study. Our findings suggest that methylation in CPT1A is a promising epigenetic marker for MetS risk which could become useful as a treatment target in the future.

DNA methylation at CpG sites is both heritable and influenced by environment, but the relative contributions of each to DNA methylation levels are unclear. We conducted a heritability analysis of CpG methylation in human CD4+ cells across 975 individuals from 163 families in the Genetics of Lipid-lowering Drugs and Diet Network (GOLDN). Based on a broad-sense heritability (H2) value threshold of 0.4, we identified 20,575 highly heritable CpGs among the 174,445 most variable autosomal CpGs (SD > 0.02). Tests for associations of heritable CpGs with genotype at 2,145,360 SNPs using 717 of 975 individuals showed that ~74% were cis-meQTLs (< 1 Mb away from the CpG), 6% of CpGs exhibited trans-meQTL associations (>1 Mb away from the CpG or located on a different chromosome), and 20% of CpGs showed no strong significant associations with genotype (based on a p-value threshold of 1e-7). Genes proximal to the genotype independent heritable CpGs were enriched for functional terms related to regulation of T cell activation. These CpGs were also among those that distinguished T cells from other blood cell lineages. Compared to genes proximal to meQTL-associated heritable CpGs, genotype independent heritable CpGs were moderately enriched in the same genomic regions that escape erasure during primordial germ cell development and could carry potential for generational transmission.

Fenofibrate lowers triglycerides (TG) and raises high density lipoprotein cholesterol (HDLc) in dyslipidemic individuals. Several studies have shown genetic variability in lipid responses to fenofibrate treatment. It is, however, not known whether epigenetic patterns are also correlated with the changes in lipids due to fenofibrate treatment. The present study was therefore undertaken to examine the changes in DNA methylation among the participants of Genetics of Lipid Lowering Drugs and Diet Network (GOLDN) study. A total of 443 individuals were studied for epigenome-wide changes in DNA methylation, assessed using the Illumina Infinium HumanMethylation450 array, before and after a 3-week daily treatment with 160 mg of fenofibrate. The association between the change in DNA methylation and changes in TG, HDLc, and low-density lipoprotein cholesterol (LDLc) were assessed using linear mixed models adjusted for age, sex, baseline lipids, and study center as fixed effects and family as a random effect. Changes in DNA methylation were not significantly associated with changes in TG, HDLc, or LDLc after 3 weeks of fenofibrate for any CpG. CpG changes in genes known to be involved in fenofibrate response, e.g., PPAR-α, APOA1, LPL, APOA5, APOC3, CETP, and APOB, also did not show evidence of association. In conclusion, changes in lipids in response to 3-week treatment with fenofibrate were not associated with changes in DNA methylation. Studies of longer duration may be required to detect treatment-induced changes in methylation.

Abstract 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.

Proprotein convertase subtilisin/kexin type 9 (encoded by PCSK9) plays a well-known role in the regulation of low-density lipoprotein (LDL) receptors, and an inhibitor of this enzyme is a promising new therapeutic for hyperlipidemia. Recently, animal and human studies also implicate PCSK9 genetic variation in the regulation of blood pressure. The goal of this study was to examine if common and rare polymorphisms in PCSK9 are associated with blood pressure in an African-American population at high risk for cardiovascular disease. Using genomic data assayed on the Affymetrix 6.0 array (n = 1199) and the Illumina HumanExome Beadchip (n = 1966) from the Hypertension Genetic Epidemiology Network (HyperGEN), we tested the association of PCSK9 polymorphisms with blood pressure. We used linear mixed models and the sequence kernel association test (SKAT) to assess the association of 31 common and 19 rare variants with blood pressure. The models were adjusted for age, sex, center, smoking status, principal components for ancestry and diabetes as fixed effects and family as a random effect. The results showed a marginally significant effect of two genome-wide association study (GWAS) single-nucleotide polymorphisms (SNPs) (rs12048828: β = 1.8, P = 0.05 and rs9730100: β = 1.0, P = 0.05) with diastolic blood pressure (DBP); however these results were not significant after correction for multiple testing. Rare variants were cumulatively associated with DBP (P = 0.04), an effect that was strengthened by restriction to nonsynonymous or stop-gain SNPs (P = 0.02). While gene-based results for DBP did not replicate (P = 0.36), we found an association with SBP (P = 0.04) in the Reasons for Geographic And Racial Differences in Stroke study (REGARDS). The findings here suggest rare variants in PCSK9 may influence blood pressure among African Americans, laying the ground work for further validation studies.

Previous studies have reported that risk of cardiovascular morbidity and mortality substantially increases in hypertensive patients, especially among those with inadequate blood pressure control. Two common antihypertensive drug classes including thiazide diuretics and angiotensinogen converting enzyme (ACE) inhibitors affect different enzymes in the renin angiotensinogen aldosterone system (RAAS). In the RAAS, angiotensinogen is converted into angiotensin II which increases blood pressure through vasoconstriction. Using a case-only design with 3,448 high-risk hypertensive individuals from the Genetics of Hypertension Associated Treatment (GenHAT) study, we examined whether 7 single nucleotide polymorphisms (SNPs) in the angiotensinogen gene (AGT) interact with three classes of antihypertensive drugs including chlorthalidone (a thiazide diuretic), lisinopril (an ACE inhibitor), and amlodipine (a calcium channel blocker) to modify the risk of incident coronary heart disease (CHD) and heart failure (HF) among Caucasian and African American participants, separately. We found no gene by treatment interactions to be statistically significant after correction for multiple testing. However, some suggestive results were found. African American participants with the minor allele of rs11122576 had over two-fold higher risk of CHD when using chlorthalidone compared to using amlodipine, or lisinopril compared to amlodipine (p=0.006, and p=0.01, respectively). Other marginal associations are also reported among both race groups. The findings reported here suggest that rs11122576 could contribute to future personalization of antihypertensive treatment among African Americans though more studies are needed.

Cardiac hypertrophy is an independent risk factor for cardiovascular disease and heart failure. There is increasing evidence that microRNAs (miRNAs) play an important role in the regulation of messenger RNA (mRNA) and the pathogenesis of various cardiovascular diseases. However, the ability to comprehensively study cardiac hypertrophy on a gene regulatory level is impacted by the limited availability of human cardiomyocytes. Human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) offer the opportunity for disease modeling. Here we utilize a previously established in vitro model of cardiac hypertrophy to interrogate the regulatory mechanism associated with the cardiac disease process. We perform miRNA sequencing and mRNA expression analysis on endothelin 1 (ET-1) stimulated hiPSC-CMs to describe associated RNA expression profiles. MicroRNA sequencing revealed over 250 known and 34 predicted novel miRNAs to be differentially expressed between ET-1 stimulated and unstimulated control hiPSC-CMs. Messenger RNA expression analysis identified 731 probe sets with significant differential expression. Computational target prediction on significant differentially expressed miRNAs and mRNAs identified nearly 2000 target pairs. A principal component analysis approach comparing the in vitro data with human myocardial biopsies detected overlapping expression changes between the in vitro samples and myocardial biopsies with Left Ventricular Hypertrophy. These results provide further insights into the complex RNA regulatory mechanism associated with cardiac hypertrophy.

Increased postprandial lipid (PPL) response to dietary fat intake is a heritable risk factor for cardiovascular disease (CVD). Variability in postprandial lipids results from the complex interplay of dietary and genetic factors. We hypothesized that detailed lipid profiles (eg, sterols and fatty acids) may help elucidate specific genetic and dietary pathways contributing to the PPL response.We used gas chromatography mass spectrometry to quantify the change in plasma concentration of 35 fatty acids and 11 sterols between fasting and 3.5 hours after the consumption of a high-fat meal (PPL challenge) among 40 participants from the GOLDN study. Correlations between sterols, fatty acids and clinical measures were calculated. Mixed linear regression was used to evaluate associations between lipidomic profiles and genomic markers including single nucleotide polymorphisms (SNPs) and methylation markers derived from the Affymetrix 6.0 array and the Illumina Methyl450 array, respectively. After the PPL challenge, fatty acids increased as well as sterols associated with cholesterol absorption, while sterols associated with cholesterol synthesis decreased. PPL saturated fatty acids strongly correlated with triglycerides, very low-density lipoprotein, and chylomicrons. Two SNPs (rs12247017 and rs12240292) in the sorbin and SH3 domain containing 1 (SORBS1) gene were associated with b-Sitosterol after correction for multiple testing (P≤4.5*10(-10)). SORBS1 has been linked to obesity and insulin signaling. No other markers reached the genome-wide significance threshold, yet several other biologically relevant loci are highlighted (eg, PRIC285, a co-activator of PPARa).Integration of lipidomic and genomic data has the potential to identify new biomarkers of CVD risk.

Genome-wide association studies (GWAS) have identified ~100 loci associated with blood lipid levels, but much of the trait heritability remains unexplained, and at most loci the identities of the trait-influencing variants remain unknown. We conducted a trans-ethnic fine-mapping study at 18, 22, and 18 GWAS loci on the Metabochip for their association with triglycerides (TG), high-density lipoprotein cholesterol (HDL-C), and low-density lipoprotein cholesterol (LDL-C), respectively, in individuals of African American (n = 6,832), East Asian (n = 9,449), and European (n = 10,829) ancestry. We aimed to identify the variants with strongest association at each locus, identify additional and population-specific signals, refine association signals, and assess the relative significance of previously described functional variants. Among the 58 loci, 33 exhibited evidence of association at P

While genome-wide association studies (GWAS) and candidate gene approaches have identified many genetic variants that contribute to disease risk as main effects, the impact of genotype by environment (GxE) interactions remains rather under-surveyed. To explore the importance of GxE interactions for diabetes-related traits, a tool for Genome-wide Complex Trait Analysis (GCTA) was used to examine GxE variance contribution of 15 macronutrients and lifestyle to the total phenotypic variance of diabetes-related traits at the genome-wide level in a European American population. GCTA identified two key environmental factors making significant contributions to the GxE variance for diabetes-related traits: carbohydrate for fasting insulin (25.1% of total variance, P-nominal = 0.032) and homeostasis model assessment of insulin resistance (HOMA-IR) (24.2% of total variance, P-nominal = 0.035), n-6 polyunsaturated fatty acid (PUFA) for HOMA-β-cell-function (39.0% of total variance, P-nominal = 0.005). To demonstrate and support the results from GCTA, a GxE GWAS was conducted with each of the significant dietary factors and a control E factor (dietary protein), which contributed a non-significant GxE variance. We observed that GxE GWAS for the environmental factor contributing a significant GxE variance yielded more significant SNPs than the control factor. For each trait, we selected all significant SNPs produced from GxE GWAS, and conducted anew the GCTA to estimate the variance they contributed. We noted the variance contributed by these SNPs is higher than that of the control. In conclusion, we utilized a novel method that demonstrates the importance of genome-wide GxE interactions in explaining the variance of diabetes-related traits.

BACKGROUND:Anthrax and its etiologic agent remain a biological threat. Anthrax vaccine is highly effective, but vaccine-induced IgG antibody responses vary widely following required doses of vaccinations. Such variation can be related to genetic factors, especially genomic copy number variants (CNVs) that are known to be enriched among genes with immunologic function. We have tested this hypothesis in two study populations from a clinical trial of anthrax vaccination. METHODS:We performed CNV-based genome-wide association analyses separately on 794 European Americans and 200 African-Americans. Antibodies to protective antigen were measured at week 8 (early response) and week 30 (peak response) using an enzyme-linked immunosorbent assay. We used DNA microarray data (Affymetrix 6.0) and two CNV detection algorithms, hidden markov model (PennCNV) and circular binary segmentation (GeneSpring) to determine CNVs in all individuals. Multivariable regression analyses were used to identify CNV-specific associations after adjusting for relevant non-genetic covariates. RESULTS:Within the 22 autosomal chromosomes, 2,943 non-overlapping CNV regions were detected by both algorithms. Genomic insertions containing HLA-DRB5, DRB1 and DQA1/DRA genes in the major histocompatibility complex (MHC) region (chromosome 6p21.3) were moderately associated with elevated early antibody response (β = 0.14, p = 1.78×10(-3)) among European Americans, and the strongest association was observed between peak antibody response and a segmental insertion on chromosome 1, containing NBPF4, NBPF5, STXMP3, CLCC1, and GPSM2 genes (β = 1.66, p = 6.06×10(-5)). For African-Americans, segmental deletions spanning PRR20, PCDH17 and PCH68 genes on chromosome 13 were associated with elevated early antibody production (β = 0.18, p = 4.47×10(-5)). Population-specific findings aside, one genomic insertion on chromosome 17 (containing NSF, ARL17 and LRRC37A genes) was associated with elevated peak antibody response in both populations. CONCLUSION:Multiple CNV regions, including the one consisting of MHC genes that is consistent with earlier research, can be important to humoral immune responses to anthrax vaccine adsorbed.

BACKGROUND:Neighboring genes PIK3CA and KCNMB3 are both important for insulin signaling and β-cell function, but their associations with glucose-related traits are unclear. OBJECTIVE:The objective was to examine associations of PIK3CA-KCNMB3 variants with glucose-related traits and potential interaction with dietary fat. DESIGN:We first investigated genetic associations and their modulation by dietary fat in the Genetics of Lipid Lowering Drugs and Diet Network (GOLDN) study (n = 820). Nine single-nucleotide polymorphisms (SNPs) were selected for analysis, covering more than 80% of the SNPs in the region. We then sought to replicate the findings in the Boston Puerto Rican Health Study (BPRHS) (n = 844). RESULTS:For KCNMB3 missense mutation rs7645550, meta-analysis indicated that homeostasis model assessment of insulin resistance (HOMA-IR) was significantly lower in minor allele T homozygotes compared with major allele C carriers (pooled P-value = 0.004); for another SNP rs1183319, which is in moderate LD with rs7645550, minor allele G carriers had higher HOMA-IR compared with non-carriers in both populations (pooled P-value = 0.028). In GOLDN, rs7645550 T allele homozygotes had lower HOMA-IR only when dietary n-3: n-6 PUFA ratio was low (≤0.11, P = 0.001), but not when it was high (>0.11, P-interaction = 0.033). Similar interaction was observed between rs1183319 and n-3: n-6 PUFA ratio on HOMA-IR (P-interaction = 0.001) in GOLDN. Variance contribution analyses in GOLDN confirmed the genetic association and gene-diet interaction. In BPRHS, dietary n-3: n-6 PUFA ratio significantly modulated the association between rs1183319 and HbA1c (P-interaction = 0.034). CONCLUSION:PIK3CA-KCNMB3 variants are associated with insulin resistance in populations of different ancestries, and are modified by dietary PUFA.

Rationale: Left ventricular hypertrophy (LVH) is a heritable predictor of cardiovascular disease, particularly in blacks. Objective: Determine the feasibility of combining evidence from two distinct but complimentary experimental approaches to identify novel genetic predictors of increased LV mass . Methods: Whole exome sequencing (WES) was conducted in 7 African American sibling trios ascertained on high average familial LV mass indexed to height (LVMHT). WES identified 31,426 missense or nonsense mutations (MS/NS) which were examined for association with LVMHT using linear mixed models adjusted for age, sex, body weight, and family relationship. To functionally assess WES findings, human induced pluripotent stem cell-derived cardiomyocytes (iPSC-CM) were stimulated to induce hypertrophy; mRNA sequencing was used to determine expression differences associated with hypertrophy onset. Results: After correction for multiple testing, 295 MS/NS variants in 265 genes were associated with LVMHT. We identified 44 of 265 WES genes differentially expressed (P

A recent large-scale meta-analysis of genome-wide studies has identified 95 loci, 59 of them novel, as statistically significant predictors of blood lipid traits; we tested whether the same loci explain the observed heterogeneity in response to lipid-lowering therapy with fenofibrate. Using data from the Genetics of Lipid Lowering Drugs and Diet Network (GOLDN, n = 861) we fit linear mixed models with the genetic markers as predictors and high-density lipoprotein (HDL) cholesterol, low-density lipoprotein (LDL) cholesterol, total cholesterol, and triglyceride concentrations as outcomes. For all four traits, we analyzed both baseline levels and changes in response to treatment with fenofibrate. For the markers that were significantly associated with fenofibrate response, we fit additional models evaluating potential epistatic interactions. All models were adjusted for age, sex, and study center as fixed effects, and pedigree as a random effect. Statistically significant associations were observed between the rs964184 polymorphism near APOA1 (P-value≤0.0001) and fenofibrate response for HDL and triglycerides. The association was replicated in the Pharmacogenetics of Hypertriglyceridemia in Hispanics study (HyperTG, n = 267). Suggestive associations with fenofibrate response were observed for markers in or near PDE3A, MOSC1, FLJ36070, CETP, the APOE-APOC1-APOC4-APOC2, and CILP2. Finally, we present strong evidence for epistasis (P-value for interaction = 0.0006 in GOLDN, 0.05 in HyperTG) between rs10401969 near CILP2 and rs4420638 in the APOE-APOC1-APOC4-APOC2 cluster with total cholesterol response to fenofibrate. In conclusion, we present evidence linking several novel and biologically relevant genetic polymorphisms to lipid lowering drug response, as well as suggesting novel gene-gene interactions in fenofibrate pharmacogenetics.

Nitric oxide synthase 3 (NOS3) catalyzes production of NO in the endothelium and may play a role in cardiovascular disease (CVD). We assessed the pharmacogenetic associations of three NOS3 polymorphisms and three antihypertensive drugs with CVD outcomes. Hypertensive subjects (n = 30,280) from a multi-center, double-blind clinical trial were randomized to chlorthalidone, amlodipine, or lisinopril treatment (mean follow up, 4.9 years). Outcomes included coronary heart disease (CHD: fatal CHD and nonfatal myocardial infarction); stroke; heart failure (fatal, requiring hospitalization, or outpatient treatment); all-cause mortality; and end-stage renal disease (ESRD). Main effects of NOS3 variants on outcome and genotype-treatment interactions were tested. For NOS3 -690 C>T (rs3918226), a higher hazard ratio (HR) was found in minor allele carriers for CHD (CC = 1.00, CT+TT = 1.12 (95% confidence interval (CI) = 1.00-1.26), P = 0.048). For NOS3 -922 A>G (rs1800779), a higher HR was found in minor allele carriers for heart failure (AA = 1.00, AG+GG = 1.10 (CI = 1.00-1.21), P = 0.046). Significant pharmacogenetic findings were observed for stroke and all-cause mortality. For -690 C>T, a lower HR was observed for stroke in minor allele carriers when treated with amlodipine versus lisinopril (CC = 0.85 (CI = 0.73-0.99), CT+TT = 0.49 (CI = 0.31-0.80), P = 0.04). For glu298asp G>T (rs1799983), a lower HR was observed for all-cause mortality in minor allele carriers when treated with amlodipine versus lisinopril (GG = 1.01 (CI = 0.91-1.13), GT+TT = 0.85 (CI = 0.75-0.97), P = 0.04). We observed significant associations with NOS3 variants and CHD and heart failure and significant pharmacogenetic effects for stroke and all cause mortality. This suggests that NOS3 variants may potentially provide useful clinical information with respect to treatment decisions in the future.

Trans fatty acids (TFA) lower HDL and increase triglyceride concentrations while polyunsaturated fatty acids (PUFA) lower triglycerides and may decrease HDL concentrations. The effect of the interaction between trans fat and PUFA on lipids is uncertain.Men and women (n = 1032) in the Genetics of Lipid-Lowering Drugs and Diet Network (GOLDN) study were included. Fatty acids in erythrocyte membranes were measured with gas chromatography while data on potential confounders were obtained from questionnaires. To test the interaction between total erythrocyte PUFA (ePUFA) and TFA (eTFA) on lipid concentrations we distributed eTFA into tertiles and dichotomized ePUFA at the median concentration.For the 1(st), 2(nd) and 3(rd) tertiles of eTFA, multivariate-adjusted means±s.e.m for HDL were 46.2±1.1, 46.3±1.1 and 45.5±1.0 mg/dL among those with low ePUFA, respectively, while they were 50.0±1.1, 46.9±1.1 and 44.7±1.1 mg/dL among those with high ePUFA, respectively (P for interaction = 0.01). For the 1(st), 2(nd) and 3(rd) tertiles of eTFA, multivariate-adjusted means±s.e.m for triglycerides were 178.6±11.3, 144.7±10.9 and 140.8±10.6, respectively, among those with low ePUFA, while they were 133.8±11.3, 145.7±10.9 and 149.3±11.5, respectively, among those with high ePUFA (P for interaction = 0.005). Results for VLDL were similar to those for triglycerides. No significant interactions were observed for LDL or total cholesterol.The relation between trans fat and HDL, VLDL and triglycerides may depend on PUFA. The benefit of avoiding trans fat may be greater among individuals with higher PUFA intake. Supplementation with PUFA among individuals with relatively high trans fat intake may have limited benefits on lipid profiles.

The LDL receptor-related protein-1 gene (LRP-1) has been associated with obesity in animal models, but no such association has yet been reported in humans. As data suggest this increase in fat mass may be mediated through a mechanism involving the clearance of plasma triglyceride-rich lipoproteins (TGRL), where the LRP interacts with apolipoprotein E (ApoE) on chylomicron remnants, we aimed to examine (1) whether there was an association between 3 single nucleotide polymorphisms (SNPs) on LRP-1 with body mass index (BMI) and (2) whether any association between LRP-1 SNPs and BMI could be modified by polymorphisms on the ApoE gene when comparing the wild type ε3/ε3 genotype against mutant ApoE allele (ε2/ε4) carriers.We used data from 1,036 men and women (mean age ± SD = 49 ± 16 y) participating in the Genetics of Lipid Lowering Drugs and Diet Network (GOLDN) Study. Mixed linear models, which controlled for age, sex, alcohol intake and smoking, as well as family pedigree and center of data collection were calculated. Models that used LRP-1 genotype as a predictor of BMI revealed that individuals who were homozygous for the minor allele at the LRP-1 I10701 locus had BMIs, on average, 1.03 kg/m(2) higher than major allele carriers (P = 0.03). In subsequent mixed linear models that included main effects of LRP-1 I10701 SNP and ApoE alleles, and an interaction term the two genotypes, there was no interaction detected between the LRP-1 I70701 genotype with either the ApoE ε2 or ε4 allele carriers (P>0.05).This has implications for starting to understand pathways from genotype to human BMI, which may operate through TGRL uptake at the LRP-1 receptor. This may pave the way for future research into individualized dietary interventions.

Fibrates lower triglycerides and raise HDL cholesterol in dyslipidemic patients, but show heterogeneous treatment response. We used k-means clustering to identify three representative NMR lipoprotein profiles for 775 subjects from the GOLDN population, and study the response to fenofibrate in corresponding subgroups. The subjects in each subgroup showed differences in conventional lipid characteristics and in presence/absence of cardiovascular risk factors at baseline; there were subgroups with a low, medium and high degree of dyslipidemia. Modeling analysis suggests that the difference between the subgroups with low and medium dyslipidemia is influenced mainly by hepatic uptake dysfunction, while the difference between subgroups with medium and high dyslipidemia is influenced mainly by extrahepatic lipolysis disfunction. The medium and high dyslipidemia subgroups showed a positive, yet distinct lipid response to fenofibrate treatment. When comparing our subgroups to known subgrouping methods, we identified an additional 33% of the population with favorable lipid response to fenofibrate compared to a standard baseline triglyceride cutoff method. Compared to a standard HDL cholesterol cutoff method, the addition was 18%. In conclusion, by using constructing subgroups based on representative lipoprotein profiles, we have identified two subgroups of subjects with positive lipid response to fenofibrate therapy and with different underlying disturbances in lipoprotein metabolism. The total subgroup with positive lipid response to fenofibrate is larger than subgroups identified with baseline triglyceride and HDL cholesterol cutoffs.

Long-chain n-3 polyunsaturated fatty acids (PUFAs) can derive from diet or from α-linolenic acid (ALA) by elongation and desaturation. We investigated the association of common genetic variation with plasma phospholipid levels of the four major n-3 PUFAs by performing genome-wide association studies in five population-based cohorts comprising 8,866 subjects of European ancestry. Minor alleles of SNPs in FADS1 and FADS2 (desaturases) were associated with higher levels of ALA (p = 3 x 10⁻⁶⁴) and lower levels of eicosapentaenoic acid (EPA, p = 5 x 10⁻⁵⁸) and docosapentaenoic acid (DPA, p = 4 x 10⁻¹⁵⁴). Minor alleles of SNPs in ELOVL2 (elongase) were associated with higher EPA (p = 2 x 10⁻¹²) and DPA (p = 1 x 10⁻⁴³) and lower docosahexaenoic acid (DHA, p = 1 x 10⁻¹⁵). In addition to genes in the n-3 pathway, we identified a novel association of DPA with several SNPs in GCKR (glucokinase regulator, p = 1 x 10⁻⁸). We observed a weaker association between ALA and EPA among carriers of the minor allele of a representative SNP in FADS2 (rs1535), suggesting a lower rate of ALA-to-EPA conversion in these subjects. In samples of African, Chinese, and Hispanic ancestry, associations of n-3 PUFAs were similar with a representative SNP in FADS1 but less consistent with a representative SNP in ELOVL2. Our findings show that common variation in n-3 metabolic pathway genes and in GCKR influences plasma phospholipid levels of n-3 PUFAs in populations of European ancestry and, for FADS1, in other ancestries.

PLIN4 is a member of the PAT family of lipid storage droplet (LSD) proteins. Associations between seven single nucleotide polymorphisms (SNPs) at human PLIN4 with obesity related phenotypes were investigated using meta-analysis followed by a determination if these phenotypes are modulated by interactions between PLIN4 SNPs and dietary PUFA. Samples consisted of subjects from two populations of European ancestry. We demonstrated association of rs8887 with anthropometrics. Meta-analysis demonstrated significant interactions between the rs8887 minor allele with PUFA n3 modulating anthropometrics. rs884164 showed interaction with both n3 and n6 PUFA modulating anthropometric and lipid phenotypes. In silico analysis of the PLIN4 3'UTR sequence surrounding the rs8887 minor A allele predicted a seed site for the human microRNA-522 (miR-522), suggesting a functional mechanism. Our data showed that a PLIN4 3'UTR luciferase reporter carrying the A allele of rs8887 was reduced in response to miR-522 mimics compared to the G allele. These results suggest variation at the PLIN4 locus, and its interaction with PUFA as a modulator of obesity related phenotypes, acts in part through creation of a miR-522 regulatory site.

MMP-9 and -12 function in tissue remodeling and may play roles in cardiovascular disease (CVD). We assessed associations of four MMP polymorphisms and three antihypertensive drugs with cardiovascular outcomes.Hypertensives (n = 42,418) from a double-blind, randomized, clinical trial were randomized to chlorthalidone, amlodipine, lisinopril, or doxazosin treatment (mean follow up, 4.9 years). The primary outcome was coronary heart disease (CHD). Secondary outcomes included combined CHD, all CVD outcomes combined, stroke, heart failure (HF), and mortality. Genotype-treatment interactions were tested.There were 38,698 participants genotyped for at least one of the polymorphisms included here. For MMP9 R668Q (rs2274756), lower hazard ratios (HRs) were found for AA subjects for most outcomes when treated with chlorthalidone versus amlodipine (eg., CCHD: GG = 1.00, GA = 1.01, AA = 0.64; P = 0.038). For MMP9 R279Q (rs17576), modest pharmacogenetic findings were observed for combined CHD and the composite CVD outcome. For MMP12 N122S (rs652438), lower HRs were observed for CHD in subjects carrying at least one G allele and being treated with chlorthalidone versus lisinopril (CHD: AA = 1.07, AG = 0.80, GG = 0.49; P = 0.005). In the lisinopril-amlodipine comparison, higher HRs were observed for participants having at least one G allele at the MMP12 N122S locus (CHD: AA = 0.94, AG = 1.19, GG = 1.93; P = 0.041). For MMP12 -82A>G (rs2276109), no pharmacogenetic effect was found for the primary outcome, although lower HRs were observed for AA homozygotes in the chlorthalidone-amlodipine comparison for HF (P = 0.015).We observed interactions between antihypertensive drugs and MMP9 and MMP12 for CHD and composite CVD. The data suggest that these genes may provide useful clinical information with respect to treatment decisions.

African Americans have been understudied in genome wide association studies of diabetes and related traits. In the current study, we examined the joint association of single nucleotide polymorphisms (SNPs) and copy number variants (CNVs) with fasting insulin and an index of insulin resistance (HOMA-IR) in the HyperGEN study, a family based study with proband ascertainment for hypertension. This analysis is restricted to 1,040 African Americans without diabetes. We generated allele specific CNV genotypes at 872,243 autosomal loci using Birdsuite, a freely available multi-stage program. Joint tests of association for SNPs and CNVs were performed using linear mixed models adjusting for covariates and familial relationships. Our results highlight SNPs associated with fasting insulin and HOMA-IR (rs6576507 and rs8026527, 3.7*10(-7)≤P≤1.1*10(-5)) near ATPase, class V, type 10A (ATP10A), and the L Type voltage dependent calcium channel (CACNA1D, rs1401492, P≤5.2*10(-6)). ATP10A belongs to a family of aminophospholipid-transporting ATPases and has been associated with type 2 diabetes in mice. CACNA1D has been linked to pancreatic beta cell generation in mice. The two most significant copy variable markers (rs10277702 and rs361367; P

Abstract We assess performance and limitations of polygenic risk scores (PRSs) for multiple blood pressure (BP) phenotypes in diverse population groups. We compare “clumping-and-thresholding” (PRSice2) and LD-based (LDPred2) methods to construct PRSs from each of multiple GWAS, as well as multi-PRS approaches that sum PRSs with and without weights, including PRS-CSx. We use datasets from the MGB Biobank, TOPMed study, UK biobank, and from All of Us to train, assess, and validate PRSs in groups defined by self-reported race/ethnic background (Asian, Black, Hispanic/Latino, and White). For both SBP and DBP, the PRS-CSx based PRS, constructed as a weighted sum of PRSs developed from multiple independent GWAS, perform best across all race/ethnic backgrounds. Stratified analysis in All of Us shows that PRSs are better predictive of BP in females compared to males, individuals without obesity, and middle-aged (40-60 years) compared to older and younger individuals.

Introduction: Apparent treatment-resistant hypertension (aTRH) is characterized by the use of four or more antihypertensive (AHT) classes to achieve blood pressure (BP) control. In the current study, we conducted single-variant and gene-based analyses of aTRH among individuals from 12 Trans-Omics for Precision Medicine cohorts with whole-genome sequencing data.Methods: Cases were defined as individuals treated for hypertension (HTN) taking three different AHT classes, with average systolic BP ≥ 140 or diastolic BP ≥ 90 mmHg, or four or more medications regardless of BP (n = 1,705). A normotensive control group was defined as individuals with BP < 140/90 mmHg (n = 22,079), not on AHT medication. A second control group comprised individuals who were treatment responsive on one AHT medication with BP < 140/ 90 mmHg (n = 5,424). Logistic regression with kinship adjustment using the Scalable and Accurate Implementation of Generalized mixed models (SAIGE) was performed, adjusting for age, sex, and genetic ancestry. We assessed variants using SKAT-O in rare-variant analyses. Single-variant and gene-based tests were conducted in a pooled multi-ethnicity stratum, as well as self-reported ethnic/racial strata (European and African American).Results: One variant in the known HTN locus, KCNK3, was a top finding in the multi-ethnic analysis (p = 8.23E-07) for the normotensive control group [rs12476527, odds ratio (95% confidence interval) = 0.80 (0.74–0.88)]. This variant was replicated in the Vanderbilt University Medical Center’s DNA repository data. Aggregate gene-based signals included the genes AGTPBP, MYL4, PDCD4, BBS9, ERG, and IER3.Discussion: Additional work validating these loci in larger, more diverse populations, is warranted to determine whether these regions influence the pathobiology of aTRH.

Introduction: Educational attainment, widely used in epidemiologic studies as a surrogate for socioeconomic status, is a predictor of cardiovascular health outcomes.Methods: A two-stage genome-wide meta-analysis of low-density lipoprotein cholesterol (LDL), high-density lipoprotein cholesterol (HDL), and triglyceride (TG) levels was performed while accounting for gene-educational attainment interactions in up to 226,315 individuals from five population groups. We considered two educational attainment variables: “Some College” (yes/no, for any education beyond high school) and “Graduated College” (yes/no, for completing a 4-year college degree). Genome-wide significant (p < 5 × 10−8) and suggestive (p < 1 × 10−6) variants were identified in Stage 1 (in up to 108,784 individuals) through genome-wide analysis, and those variants were followed up in Stage 2 studies (in up to 117,531 individuals).Results: In combined analysis of Stages 1 and 2, we identified 18 novel lipid loci (nine for LDL, seven for HDL, and two for TG) by two degree-of-freedom (2 DF) joint tests of main and interaction effects. Four loci showed significant interaction with educational attainment. Two loci were significant only in cross-population analyses. Several loci include genes with known or suggested roles in adipose (FOXP1, MBOAT4, SKP2, STIM1, STX4), brain (BRI3, FILIP1, FOXP1, LINC00290, LMTK2, MBOAT4, MYO6, SENP6, SRGAP3, STIM1, TMEM167A, TMEM30A), and liver (BRI3, FOXP1) biology, highlighting the potential importance of brain-adipose-liver communication in the regulation of lipid metabolism. An investigation of the potential druggability of genes in identified loci resulted in five gene targets shown to interact with drugs approved by the Food and Drug Administration, including genes with roles in adipose and brain tissue.Discussion: Genome-wide interaction analysis of educational attainment identified novel lipid loci not previously detected by analyses limited to main genetic effects.

Introduction: Metabolic syndrome (MetS) increases the risk of cardiovascular disease and death. Previous ‘-omics’ studies have identified dysregulated serum metabolites and aberrant DNA methylation in the setting of MetS. However, the relationship between the metabolome and epigenome have not been elucidated. In this study, we identified serum metabolites associated with MetS and DNA methylation, and we conducted bidirectional Mendelian randomization (MR) to assess causal relationships between metabolites and methylation.Methods: We leveraged metabolomic and genomic data from a national United States cohort of older adults (REGARDS), as well as metabolomic, epigenomic, and genomic data from a family-based study of hypertension (HyperGEN). We conducted metabolite profiling for MetS in REGARDS using weighted logistic regression models and validated them in HyperGEN. Validated metabolites were selected for methylation studies which fit linear mixed models between metabolites and six CpG sites previously linked to MetS. Statistically significant metabolite-CpG pairs were selected for two-sample, bidirectional MR.Results: Forward MR indicated that glucose and serine metabolites were causal on CpG methylation near CPT1A [B(SE): −0.003 (0.002), p = 0.028 and B(SE): 0.029 (0.011), p = 0.030, respectively] and that serine metabolites were causal on ABCG1 [B(SE): −0.008(0.003), p = 0.006] and SREBF1 [B(SE): −0.009(0.004), p = 0.018] methylation, which suggested a protective effect of serine. Reverse MR showed a bidirectional relationship between cg06500161 (ABCG1) and serine [B(SE): −1.534 (0.668), p = 0.023].Discussion: The metabolome may contribute to the relationship between MetS and epigenetic modifications.

Summary: Background: Identifying novel epigenetic signatures associated with serum immunoglobulin E (IgE) may improve our understanding of molecular mechanisms underlying asthma and IgE-mediated diseases. Methods: We performed an epigenome-wide association study using whole blood from Framingham Heart Study (FHS; n = 3,471, 46% females) participants and validated results using the Childhood Asthma Management Program (CAMP; n = 674, 39% females) and the Genetic Epidemiology of Asthma in Costa Rica Study (CRA; n = 787, 41% females). Using the closest gene to each IgE-associated CpG, we highlighted biologically plausible pathways underlying IgE regulation and analyzed the transcription patterns linked to IgE-associated CpGs (expression quantitative trait methylation loci; eQTMs). Using prior UK Biobank summary data from genome-wide association studies of asthma and allergy, we performed Mendelian randomization (MR) for causal inference testing using the IgE-associated CpGs from FHS with methylation quantitative trait loci (mQTLs) as instrumental variables. Findings: We identified 490 statistically significant differentially methylated CpGs associated with IgE in FHS, of which 193 (39.3%) replicated in CAMP and CRA (FDR < 0.05). Gene ontology analysis revealed enrichment in pathways related to transcription factor binding, asthma, and other immunological processes. eQTM analysis identified 124 cis-eQTMs for 106 expressed genes (FDR < 0.05). MR in combination with drug-target analysis revealed CTSB and USP20 as putatively causal regulators of IgE levels (Bonferroni adjusted P

Background: Many epigenetic loci have been associated with plasma triglyceride (TG) levels, but epigenetic connections between those loci and dietary exposures are largely unknown. This study aimed to characterize the epigenetic links between diet, lifestyle, and TG.Methods: We first conducted an epigenome-wide association study (EWAS) for TG in the Framingham Heart Study Offspring population (FHS, n = 2,264). We then examined relationships between dietary and lifestyle-related variables, collected four times in 13 years, and differential DNA methylation sites (DMSs) associated with the last TG measures. Third, we conducted a mediation analysis to evaluate the causal relationships between diet-related variables and TG. Finally, we replicated three steps to validate identified DMSs associated with alcohol and carbohydrate intake in the Genetics of Lipid-Lowering Drugs and Diet Network (GOLDN) study (n = 993).Results: In the FHS, the EWAS revealed 28 TG-associated DMSs at 19 gene regions. We identified 102 unique associations between these DMSs and one or more dietary and lifestyle-related variables. Alcohol and carbohydrate intake showed the most significant and consistent associations with 11 TG-associated DMSs. Mediation analyses demonstrated that alcohol and carbohydrate intake independently affect TG via DMSs as mediators. Higher alcohol intake was associated with lower methylation at seven DMSs and higher TG. In contrast, increased carbohydrate intake was associated with higher DNA methylation at two DMSs (CPT1A and SLC7A11) and lower TG. Validation in the GOLDN further supports the findings.Conclusion: Our findings imply that TG-associated DMSs reflect dietary intakes, particularly alcoholic drinks, which could affect the current cardiometabolic risk via epigenetic changes. This study illustrates a new method to map epigenetic signatures of environmental factors for disease risk. Identification of epigenetic markers of dietary intake can provide insight into an individual’s risk of cardiovascular disease and support the application of precision nutrition.Clinical Trial Registration:www.ClinicalTrials.gov, the Framingham Heart Study (FHS), NCT00005121; the Genetics of Lipid Lowering Drugs and Diet Network (GOLDN), NCT01023750.

Target organ damage (TOD) manifests as vascular injuries in the body organ systems associated with long-standing hypertension. DNA methylation in peripheral blood leukocytes can capture inflammatory processes and gene expression changes underlying TOD. We investigated the association between epigenome-wide DNA methylation and five measures of TOD (estimated glomerular filtration rate (eGFR), urinary albumin-creatinine ratio (UACR), left ventricular mass index (LVMI), relative wall thickness (RWT), and white matter hyperintensity (WMH)) in 961 African Americans from hypertensive sibships. A multivariate (multi-trait) model of eGFR, UACR, LVMI, and RWT identified seven CpGs associated with at least one of the traits (cg21134922, cg04816311 near C7orf50, cg09155024, cg10254690 near OAT, cg07660512, cg12661888 near IFT43, and cg02264946 near CATSPERD) at FDR q

DNA methylation (DNAm) is an important epigenetic process involved in the regulation of gene expression. While many studies have identified thousands of loci associated with age, few have differentiated between linear and non-linear DNAm trends with age. Non-linear trends could indicate early- or late-life gene regulatory processes. Using data from the Illumina 450K array on 336 human peripheral blood samples, we identified 21 CpG sites that associated with age (P