Search

To investigate additive and interactive associations of food allergies with three glutathione S-transferase (GST) genes in relation to ASD and ASD severity in Jamaican children. Using data from 344 1:1 age- and sex-matched ASD cases and typically developing controls, we assessed additive and interactive associations of food allergies with polymorphisms in GST genes (GSTM1, GSTP1 and GSTT1) in relation to ASD by applying conditional logistic regression models, and in relation to ASD severity in ASD cases as measured by the Autism Diagnostic Observation Schedule-2nd Edition (ADOS-2) total and domains specific comparison scores (CSs) by fitting general linear models. Although food allergies and GST genes were not associated with ASD, ASD cases allergic to non-dairy food had higher mean ADOS-2 Restricted and Repetitive Behaviors (RRB) CS (8.8 vs. 8.0, P = 0.04). In addition, allergy to dairy was associated with higher mean RRB CS only among ASD cases with GSTT1 DD genotype (9.9 vs. 7.8, P 0.01, interaction P = 0.01), and GSTP1 Val/Val genotype under a recessive genetic model (9.8 vs. 7.8, P = 0.02, interaction P = 0.06). Our findings are consistent with the role for GST genes in ASD and food allergies, though require replication in other populations.

We conducted cohort- and race-specific epigenome-wide association analyses of mitochondrial deoxyribonucleic acid (mtDNA) copy number (mtDNA CN) measured in whole blood from participants of African and European origins in five cohorts (n = 6182, mean age = 57–67 years, 65% women). In the meta-analysis of all the participants, we discovered 21 mtDNA CN-associated DNA methylation sites (CpG) (P

Genome-wide association studies (GWASs) have identified regions associated with chronic obstructive pulmonary disease (COPD). GWASs of other diseases have shown an approximately 10-fold overrepresentation of nonsynonymous variants, despite limited exonic coverage on genotyping arrays. We hypothesized that a large-scale analysis of coding variants could discover novel genetic associations with COPD, including rare variants with large effect sizes. We performed a meta-analysis of exome arrays from 218,399 controls and 33,851 moderate-to-severe COPD cases. All exome-wide significant associations were present in regions previously identified by GWAS. We did not identify any novel rare coding variants with large effect sizes. Within GWAS regions on chromosomes 5q, 6p, and 15q, four coding variants were conditionally significant ( P < 0.00015) when adjusting for lead GWAS single-nucleotide polymorphisms A common gasdermin B ( GSDMB) splice variant (rs11078928) previously associated with a decreased risk for asthma was nominally associated with a decreased risk for COPD [minor allele frequency (MAF) = 0.46, P = 1.8e-4]. Two stop variants in coiled-coil α-helical rod protein 1 ( CCHCR1), a gene involved in regulating cell proliferation, were associated with COPD (both P < 0.0001). The SERPINA1 Z allele was associated with a random-effects odds ratio of 1.43 for COPD (95% confidence interval = 1.17–1.74), though with marked heterogeneity across studies. Overall, COPD-associated exonic variants were identified in genes involved in DNA methylation, cell-matrix interactions, cell proliferation, and cell death. In conclusion, we performed the largest exome array meta-analysis of COPD to date and identified potential functional coding variants. Future studies are needed to identify rarer variants and further define the role of coding variants in COPD pathogenesis.

Aluminum (Al) is a metallic toxicant at high concentrations following natural or unnatural exposures. Dietary intake is considered as the main source of aluminum exposure in children. We used data from 366 typically developing (TD) children (ages 2–8 years) who participated as controls in an age- and sex-matched case–control study in Jamaica. We investigated additive and interactive associations among environmental factors and children’s genotypes for glutathione S-transferase (GST) genes (GSTT1, GSTM1, GSTP1), in relation to having a detectable blood aluminum concentration (BAlC) of >5.0 μg/L, using multivariable logistic regression models. Findings from interactive models revealed that the odds of having a detectable BAlC was significantly higher among children who ate string beans (p ≤ 0.01), whereas about 40% lower odds of having a detectable BAlC was observed in children with higher parental education level, (p = 0.02). A significant interaction between consumption of saltwater fish and GSTP1 in relation to having a detectable BAlC using either co-dominant or dominant genetic models (overall interaction p = 0.02 for both models) indicated that consumption of saltwater fish was associated with higher odds of having a detectable BAlC only among children with the GSTP1 Ile105Val Ile/Ile genotype using either co-dominant or dominant models [OR (95% CI) = 2.73 (1.07, 6.96), p = 0.04; and OR (95% CI) = 2.74 (1.08, 6.99), p = 0.03]. Since this is the first study from Jamaica that reports such findings, replication in other populations is warranted.

Mitochondrial DNA copy number (mtDNA-CN) measured from blood specimens is a minimally invasive marker of mitochondrial function that exhibits both inter-individual and intercellular variation. To identify genes involved in regulating mitochondrial function, we performed a genome-wide association study (GWAS) in 465,809 White individuals from the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) consortium and the UK Biobank (UKB). We identified 133 SNPs with statistically significant, independent effects associated with mtDNA-CN across 100 loci. A combination of fine-mapping, variant annotation, and co-localization analyses was used to prioritize genes within each of the 133 independent sites. Putative causal genes were enriched for known mitochondrial DNA depletion syndromes (p = 3.09 × 10–15) and the gene ontology (GO) terms for mtDNA metabolism (p = 1.43 × 10–8) and mtDNA replication (p = 1.2 × 10–7). A clustering approach leveraged pleiotropy between mtDNA-CN associated SNPs and 41 mtDNA-CN associated phenotypes to identify functional domains, revealing three distinct groups, including platelet activation, megakaryocyte proliferation, and mtDNA metabolism. Finally, using mitochondrial SNPs, we establish causal relationships between mitochondrial function and a variety of blood cell-related traits, kidney function, liver function and overall (p = 0.044) and non-cancer mortality (p = 6.56 × 10–4).

Arsenic (As) is a metalloid that has been classified as a xenobiotic with toxic effects on human beings, especially on children. Since the soil in Jamaica contains As, dietary intake is considered the main source of As exposure in Jamaicans. In addition, glutathione S-transferase (GST) genes, including GSTT1, GSTP1, and GSTM1, play an important role in the metabolism of xenobiotics including As in humans. Using data from 375 typically developing children (2–8 years) in Jamaica, we investigated the environmental and sociodemographic factors, as well as their possible interactions with the children’s genotype for GST genes in relation to having a detectable level of blood As concentration (i.e., >1.3 μg/L). Using multivariable logistic regression, we have identified environmental factors significantly associated with blood As concentrations that include a child’s age, parental education levels, and the consumption of saltwater fish, cabbage, broad beans, and avocado (all p < 0.01). Based on the multivariable analysis including gene x environment interactions, we found that among children with the Ile/Ile genotype for GSTP1 Ile105Val, children who consumed avocado had higher odds of having a detectable blood As concentration compared to children who did not eat avocado.

Cigarette smoking is associated with epigenetic changes that may be reversible following smoking cessation. Whole blood DNA methylation was evaluated in Framingham Heart Study Offspring (n = 169) and Third Generation (n = 30) cohort participants at two study visits 6 years apart and in Atherosclerosis Risk in Communities (ARIC) study (n = 222) participants at two study visits 20 years apart. Changes in DNA methylation (delta β values) at 483,565 cytosine-phosphate-guanine (CpG) sites and differentially methylated regions (DMRs) were compared between participants who were current, former, or never smokers at both visits (current-current, former-former, never-never, respectively), versus those who quit in the interim (current-former). Interim quitters had more hypermethylation at four CpGs annotated to AHRR, one CpG annotated to F2RL3, and one intergenic CpG (cg21566642) compared with current-current smokers (FDR < 0.02 for all), and two significant DMRs were identified. While there were no significant differentially methylated CpGs in the comparison of interim quitters and former-former smokers, 106 DMRs overlapping with small nucleolar RNA were identified. As compared with all non-smokers, current-current smokers additionally had more hypermethylation at two CpG sites annotated to HIVEP3 and TMEM126A, respectively, and another intergenic CpG (cg14339116). Gene transcripts associated with smoking cessation were implicated in immune responses, cell homoeostasis, and apoptosis. Smoking cessation is associated with early reversion of blood DNA methylation changes at CpG sites annotated to AHRR and F2RL3 towards those of never smokers. Associated gene expression suggests a role of longitudinal smoking-related DNA methylation changes in immune response processes.

This study investigated whether the concentrations of four metals [lead (Pb), mercury (Hg), manganese (Mn), and aluminum (Al)] are correlated in cord blood and childhood blood samples from Jamaican children. Cord blood samples were obtained from 21 pregnant women enrolled in the second Jamaican Birth Cohort Study from July 1, 2011 to September 30, 2011, and blood samples were drawn from their children who participated in a follow up study when the children were 4–8 years old. Correlations were assessed by the Pearson or the Spearman’s rank correlation coefficient. The mean ages of children at the childhood visit and their mother at the child’s birth were 5.5 years and 29.8 years, respectively. About 47.6% of children were male. Statistically significant correlations between cord blood and childhood blood concentrations of Pb (r(Spearman)= 0.45; P = 0.04) and Mn (r(Pearson) = 0.48; P = 0.03) were found, and these remained significant when adjusted for the child’s sex, age, or both. For Al and Hg, r(Spearman)= 0.29 and 0.08, respectively, but the correlations were not statistically significant (both P ≥ 0.20). A significant correlation between cord blood and childhood blood Pb concentrations for children 4–8 years old has not been previously reported.

Here we examine the association between DNA methylation in circulating leukocytes and blood lipids in a multi-ethnic sample of 16,265 subjects. We identify 148, 35, and 4 novel associations among Europeans, African Americans, and Hispanics, respectively, and an additional 186 novel associations through a trans-ethnic meta-analysis. We observe a high concordance in the direction of effects across racial/ethnic groups, a high correlation of effect sizes between high-density lipoprotein and triglycerides, a modest overlap of associations with epigenome-wide association studies of other cardio-metabolic traits, and a largely non-overlap with lipid loci identified to date through genome-wide association studies. Thirty CpGs reached significance in at least 2 racial/ethnic groups including 7 that showed association with the expression of an annotated gene. CpGs annotated to CPT1A showed evidence of being influenced by triglycerides levels. DNA methylation levels of circulating leukocytes show robust and consistent association with blood lipid levels across multiple racial/ethnic groups., Abnormal blood lipid levels are important risk factors for cardiovascular and other various diseases. Here the authors conduct a large-scale multi-ethnic epigenome-wide association study combined with epigenetic (cis-QTL and eQTM) data, and identify CpG-lipid traits associations that are specific to or common across racial/ethnic groups.

ObjectiveTryptophan can be catabolised to various metabolites through host kynurenine and microbial indole pathways. We aimed to examine relationships of host and microbial tryptophan metabolites with incident type 2 diabetes (T2D), host genetics, diet and gut microbiota.MethodWe analysed associations between circulating levels of 11 tryptophan metabolites and incident T2D in 9180 participants of diverse racial/ethnic backgrounds from five cohorts. We examined host genome-wide variants, dietary intake and gut microbiome associated with these metabolites.ResultsTryptophan, four kynurenine-pathway metabolites (kynurenine, kynurenate, xanthurenate and quinolinate) and indolelactate were positively associated with T2D risk, while indolepropionate was inversely associated with T2D risk. We identified multiple host genetic variants, dietary factors, gut bacteria and their potential interplay associated with these T2D-relaetd metabolites. Intakes of fibre-rich foods, but not protein/tryptophan-rich foods, were the dietary factors most strongly associated with tryptophan metabolites. The fibre-indolepropionate association was partially explained by indolepropionate-associated gut bacteria, mostly fibre-usingFirmicutes. We identified a novel association between a host functionalLCTvariant (determining lactase persistence) and serum indolepropionate, which might be related to a host gene-diet interaction on gutBifidobacterium, a probiotic bacterium significantly associated with indolepropionate independent of other fibre-related bacteria. Higher milk intake was associated with higher levels of gutBifidobacteriumand serum indolepropionate only among genetically lactase non-persistent individuals.ConclusionHigher milk intake among lactase non-persistent individuals, and higher fibre intake were associated with a favourable profile of circulating tryptophan metabolites for T2D, potentially through the host–microbial cross-talk shifting tryptophan metabolism toward gut microbial indolepropionate production.

Using data from 266 age- and sex-matched pairs of Jamaican children with autism spectrum disorder (ASD) and typically developing (TD) controls (2–8 years), we investigated whether glutathione S-transferase theta 1 (GSTT1) modifies the association between blood manganese concentrations (BMC) and ASD. After adjusting conditional logistic regression models for socioeconomic status and the interaction between GSTT1 and GSTP1 (glutathione S-transferase pi 1), using a recessive genetic model for GSTT1 and either a co-dominant or dominant model for GSTP1, the interaction between GSTT1 and BMC was significant (P = 0.02, P = 0.01, respectively). Compared to controls, ASD cases with GSTT1-DD genotype had 4.33 and 4.34 times higher odds of BMC > 12 vs. ≤ 8.3 μg/L, respectively. Replication in other populations is warranted.

Background GSTM1 encodes glutathione S-transferase μ-1 (GSTM1), which belongs to a superfamily of phase 2 antioxidant enzymes. The highly prevalent GSTM1 deletion variant is associated with kidney disease progression in human cohorts: the African American Study of Kidney Disease and Hypertension and the Atherosclerosis Risk in Communities (ARIC) Study. Methods We generated a Gstm1 knockout mouse line to study its role in a CKD model (involving subtotal nephrectomy) and a hypertension model (induced by angiotensin II). We examined the effect of intake of cruciferous vegetables and GSTM1 genotypes on kidney disease in mice as well as in human ARIC study participants. We also examined the importance of superoxide in the mediating pathways and of hematopoietic GSTM1 on renal inflammation. Results Gstm1 knockout mice displayed increased oxidative stress, kidney injury, and inflammation in both models. The central mechanism for kidney injury is likely mediated by oxidative stress, because treatment with Tempol, an superoxide dismutase mimetic, rescued kidney injury in knockout mice without lowering BP. Bone marrow crosstransplantation revealed that Gstm1 deletion in the parenchyma, and not in bone marrow-derived cells, drives renal inflammation. Furthermore, supplementation with cruciferous broccoli powder rich in the precursor to antioxidant-activating sulforaphane significantly ameliorated kidney injury in Gstm1 knockout, but not wild-type mice. Similarly, among humans (ARIC study participants), high consumption of cruciferous vegetables was associated with fewer kidney failure events compared with low consumption, but this association was observed primarily in participants homozygous for the GSTM1 deletion variant. Conclusions Our data support a role for the GSTM1 enzyme in the modulation of oxidative stress, inflammation, and protective metabolites in CKD.

Mode of delivery, preterm birth, and low birth weight (LBW) are hypothesized to be associated with Autism Spectrum Disorder (ASD) in the offspring. Using data from 343 ASD cases (2–8 years) and their age- and sex-matched typically developing controls in Jamaica we investigated these hypotheses. Our statistical analyses revealed that the parish of residence could modify the association between cesarean delivery and ASD, with a difference found in this relationship in Kingston parish [matched odds ratio (MOR) (95% confidence interval (CI)): 2.30 (1.17–4.53)] and other parishes [MOR (95% CI): 0.87 (0.48–1.59)]. Although the associations of LBW and preterm birth with ASD were not significant, we observed a significant interaction between LBW and the household socioeconomic status. These findings require replication.

Interpretation of genetic association results is difficult because signals often lack biological context. To generate hypotheses of the functional genetic etiology of complex cardiometabolic traits, we estimated the genetically determined component of gene expression from common variants using PrediXcan (1) and determined genes with differential predicted expression by trait. PrediXcan imputes tissue-specific expression levels from genetic variation using variant-level effect on gene expression in transcriptome data. To explore the value of imputed genetically regulated gene expression (GReX) models across different ancestral populations, we evaluated imputed expression levels for predictive accuracy genome-wide in RNA sequence data in samples drawn from European-ancestry and African-ancestry populations and identified substantial predictive power using European-derived models in a non-European target population. We then tested the association of GReX on 15 cardiometabolic traits including blood lipid levels, body mass index, height, blood pressure, fasting glucose and insulin, RR interval, fibrinogen level, factor VII level and white blood cell and platelet counts in 15 755 individuals across three ancestry groups, resulting in 20 novel gene-phenotype associations reaching experiment-wide significance across ancestries. In addition, we identified 18 significant novel gene-phenotype associations in our ancestry-specific analyses. Top associations were assessed for additional support via query of S-PrediXcan (2) results derived from publicly available genome-wide association studies summary data. Collectively, these findings illustrate the utility of transcriptome-based imputation models for discovery of cardiometabolic effect genes in a diverse dataset.

© The Author(s) 2019. A person's lipid profile is influenced by genetic variants and alcohol consumption, but the contribution of interactions between these exposures has not been studied. We therefore incorporated gene-Alcohol interactions into a multiancestry genome-wide association study of levels of high-density lipoprotein cholesterol, low-density lipoprotein cholesterol, and triglycerides. We included 45 studies in stage 1 (genome-wide discovery) and 66 studies in stage 2 (focused follow-up), for a total of 394,584 individuals from 5 ancestry groups. Analyses covered the period July 2014-November 2017. Genetic main effects and interaction effects were jointly assessed by means of a 2- degrees-of-freedom (df) test, and a 1-df test was used to assess the interaction effects alone. Variants at 495 loci were at least suggestively associated (P < 1 × 10?6) with lipid levels in stage 1 and were evaluated in stage 2, followed by combined analyses of stage 1 and stage 2. In the combined analysis of stages 1 and 2, a total of 147 independent loci were associated with lipid levels at P < 5 × 10?8 using 2-df tests, of which 18 were novel. No genome-wide-significant associations were found testing the interaction effect alone. The novel loci included several genes (proprotein convertase subtilisin/kexin type 5 (PCSK5), vascular endothelial growth factor B (VEGFB), and apolipoprotein B mRNA editing enzyme, catalytic polypeptide 1 (APOBEC1) complementation factor (A1CF)) that have a putative role in lipid metabolism on the basis of existing evidence from cellular and experimental models. National Heart, Lung, and Blood Institute; American Heart Association Grant

H63D has been identified as a novel locus associated with the development of hypertension. The quantitative risks for hypertension, cardiac remodeling, and adverse events are not well studied. We analyzed white participants from the ARIC study (Atherosclerosis Risk in Communities) with H63D genotyping (N=10 902). We related genotype status to prevalence of hypertension at each of 5 study visits and risk for adverse cardiovascular events. Among visit 5 participants (N=4507), we related genotype status to echocardiographic features. Frequencies of wild type (WT)/WT, H63D/WT, and H63D/H63D were 73%, 24.6%, and 2.4%. The average age at baseline was 54.9±5.7 years and 47% were men. Participants carrying the H63D variant had higher systolic blood pressure ( P =0.004), diastolic blood pressure (0.012), and more frequently had hypertension ( P 65 years. After 25 years of follow-up, there was no relationship between genotype status and any outcome ( P >0.05). H63D/WT and H63D/H63D genotypes were associated with small differences in cardiac remodeling. In conclusion, the HFE H63D variant confers an increased risk for hypertension per allele and, given its frequency, accounts for a significant number of cases of hypertension. However, there was no increased risk for adverse cardiovascular events or substantial left ventricular remodeling.

Background: DNA methylation (DNAm)-based measures of aging, termed epigenetic clocks (EC), are associated with aging-related outcomes including cardiovascular disease (CVD) and all-cause mortality. Associations of ECs with heart failure (HF) are unclear. We tested whether ECs were positively associated with risk of incident HF in the Atherosclerosis Risk in Communities (ARIC) study and evaluated whether adding ECs to Pooled Cohort Equation (PCE) variables improved risk prediction. Methods: We measured DNAm in peripheral blood leukocytes in 2,263 African American (mean age=56.3 years) and 925 European American (mean age=59.5 years) participants using the Illumina HM450K and calculated 7 ECs: Horvath, Hannum, extrinsic (EEAA) and intrinsic (IEAA) epigenetic age acceleration, Hannum, PhenoAge, and GrimAge. HF was ascertained by ICD- 9 code 428 and adjudication by an expert panel. We carried out race stratified proportional hazards regression to test associations ECs with incident HF, adjusting for PCE variables: age, sex, smoking, total cholesterol, HDL, systolic blood pressure (SBP), antihypertensive medication use, and diabetes. We calculated area under the curve (AUC) and integrated discrimination index (IDI) to evaluate improvement in risk prediction when adding the ECs to PCE variables. Results: The number of incident HF events and mean follow-up time in African Americans and European Americans were 640 (189 in the first 10 years) and 19.3 years, and 191 and 21.7 years, respectively. All 7 ECs were positively associated with HF in both African Americans and European Americans. In African Americans with follow-up restricted to the first 10 years, the HR for a one SD increment in GrimAge (5.64 years) was 1.57 (95% CI=1.31, 1.88), comparable to that for a one-SD (5.82 years) increment in age (HR=1.58, 95% CI=1.36, 1.83) and greater than that for a one-SD (20.2mmHg) increment in SBP (HR=1.33, 95% CI=1.18, 1.51). In European Americans across the entire follow-up period, the HR for a one-SD increment in GrimAge (6.13 years) was 1.22 (95% CI=1.06, 1.41), smaller than that for a one-SD (5.50 years) increment in age (HR=1.93, 95% CI=1.63, 2.29) and larger than that for a one-SD (17.9 mmHG) increment in SBP (HR=1.13, 95% CI=0.98, 1.30). In African Americans with follow-up restricted to the first 10 years, adding GrimAge to PCE variables increased AUC by 0.019 (95% CI=0.003, 0.035) and the IDI was 0.010 (95% CI=0.002, 0.019). In European Americans, adding GrimAge did not change AUC appreciably (0.004, 95% CI=-0.006, 0.014) and the IDI was 0.002 (95% CI=0.000, 0.005). Conclusion: ECs are positively associated with HF in African American and European American participants independent of traditional CVD risk factors. GrimAge modestly improved heart failure risk prediction in African Americans. HF-specific DNAm-based measures should be developed and evaluated for improvement in risk prediction.

Aims Inflammation plays an important role in cardiovascular disease (CVD) development. The NOD-like receptor protein-3 (NLRP3) inflammasome contributes to the development of atherosclerosis in animal models. Components of the NLRP3 inflammasome pathway such as interleukin-1β can therapeutically be targeted. Associations of genetically determined inflammasome-mediated systemic inflammation with CVD and mortality in humans are unknown. Methods and results We explored the association of genetic NLRP3 variants with prevalent CVD and cardiovascular mortality in 538 167 subjects on the individual participant level in an explorative gene-centric approach without performing multiple testing. Functional relevance of single-nucleotide polymorphisms on NLRP3 inflammasome activation has been evaluated in monocyte-enriched peripheral blood mononuclear cells (PBMCs). Genetic analyses identified the highly prevalent (minor allele frequency 39.9%) intronic NLRP3 variant rs10754555 to affect NLRP3 gene expression. rs10754555 carriers showed significantly higher C-reactive protein and serum amyloid A plasma levels. Carriers of the G allele showed higher NLRP3 inflammasome activation in isolated human PBMCs. In carriers of the rs10754555 variant, the prevalence of coronary artery disease was significantly higher as compared to non-carriers with a significant interaction between rs10754555 and age. Importantly, rs10754555 carriers had significantly higher risk for cardiovascular mortality during follow-up. Inflammasome inducers (e.g. urate, triglycerides, apolipoprotein C3) modulated the association between rs10754555 and mortality. Conclusion The NLRP3 intronic variant rs10754555 is associated with increased systemic inflammation, inflammasome activation, prevalent coronary artery disease, and mortality. This study provides evidence for a substantial role of genetically driven systemic inflammation in CVD and highlights the NLRP3 inflammasome as a therapeutic target.

We investigated interactive roles of three metabolic glutathione S-transferase (GST) genes (GSTP1, GSTT1, and GSTM1) and autism spectrum disorder (ASD) status in relation to blood Hg concentrations (BHC) of Jamaican children. We used data from 266 children (2-8 years) with ASD and their 1:1 age- and sex-matched typically developing (TD) controls. After adjusting General Linear Models for child’s age, socioeconomic status, consumption of leafy vegetables, fried plantain, canned fish, and the interaction between GSTP1 and GSTT1, we found significant interactions between GSTP1 and ASD status in relation to BHC either in a co-dominant or dominant genetic model for GSTP1(P &lt, 0.001, P = 0.007, respectively). In the co-dominant model for the Ile105Val GSTP1 polymorphism, geometric mean (GM) BHC in ASD cases with genotype Ile/Ile were significantly higher than in cases with the Ile/Val genotype (0.73 vs. 0.48 µg/L, P = 0.01). In contrast, in TD controls with the Ile/Val genotype GM BHC were significantly higher than in those with the Ile/Ile genotype (0.72 vs. 0.49 µg/L, P = 0.03) or the Val/Val genotype (0.72 vs. 0.51 µg/L, P = 0.04). Although our findings are consistent with the role of GSTP1 in detoxification of Hg, replication in other populations is warranted.

Mitochondrial DNA copy number (mtDNA-CN) measured from blood specimens is a minimally invasive marker of mitochondrial function that exhibits both inter-individual and intercellular variation. To identify genes involved in regulating mitochondrial function, we performed a genome-wide association study (GWAS) in 465,809 White individuals from the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) consortium and the UK Biobank (UKB). We identified 133 SNPs with statistically significant, independent effects associated with mtDNA-CN across 100 loci. A combination of fine-mapping, variant annotation, and co-localization analyses were used to prioritize genes within each of the 133 independent sites. Putative causal genes were enriched for known mitochondrial DNA depletion syndromes (p= 3.09 x 10−15) and the gene ontology (GO) terms for mtDNA metabolism (p= 1.43 x 10−8) and mtDNA replication (p= 1.2 x 10−7). A clustering approach leveraged pleiotropy between mtDNA-CN associated SNPs and 41 mtDNA-CN associated phenotypes to identify functional domains, revealing three distinct groups, including platelet activation, megakaryocyte proliferation, and mtDNA metabolism. Finally, using mitochondrial SNPs, we establish causal relationships between mitochondrial function and a variety of blood cell related traits, kidney function, liver function and overall (p= 0.044) and non-cancer mortality (p= 6.56 x 10−4).

Low muscle strength is an important heritable indicator of poor health linked to morbidity and mortality in older people. In a genome-wide association study meta-analysis of 256,523 Europeans aged 60 years and over from 22 cohorts we identify 15 loci associated with muscle weakness (European Working Group on Sarcopenia in Older People definition: n = 48,596 cases, 18.9% of total), including 12 loci not implicated in previous analyses of continuous measures of grip strength. Loci include genes reportedly involved in autoimmune disease (HLA-DQA1 p = 4 × 10−17), arthritis (GDF5 p = 4 × 10−13), cell cycle control and cancer protection, regulation of transcription, and others involved in the development and maintenance of the musculoskeletal system. Using Mendelian randomization we report possible overlapping causal pathways, including diabetes susceptibility, haematological parameters, and the immune system. We conclude that muscle weakness in older adults has distinct mechanisms from continuous strength, including several pathways considered to be hallmarks of ageing., Muscle weakness has been associated with morbidity and mortality in older people. Here, the authors have investigated this trait further by performing a genome-wide meta-analysis of grip strength and Mendelian randomization to discover causal relationships between muscle weakness and other diseases.

Genotype-phenotype association studies often combine phenotype data from multiple studies to increase statistical power. Harmonization of the data usually requires substantial effort due to heterogeneity in phenotype definitions, study design, data collection procedures, and data-set organization. Here we describe a centralized system for phenotype harmonization that includes input from phenotype domain and study experts, quality control, documentation, reproducible results, and data-sharing mechanisms. This system was developed for the National Heart, Lung, and Blood Institute’s Trans-Omics for Precision Medicine (TOPMed) program, which is generating genomic and other -omics data for more than 80 studies with extensive phenotype data. To date, 63 phenotypes have been harmonized across thousands of participants (recruited in 1948–2012) from up to 17 studies per phenotype. Here we discuss challenges in this undertaking and how they were addressed. The harmonized phenotype data and associated documentation have been submitted to National Institutes of Health data repositories for controlled access by the scientific community. We also provide materials to facilitate future harmonization efforts by the community, which include 1) the software code used to generate the 63 harmonized phenotypes, enabling others to reproduce, modify, or extend these harmonizations to additional studies, and 2) the results of labeling thousands of phenotype variables with controlled vocabulary terms.

Glutathione S-transferases (GST) are involved in the detoxification of exogenous chemicals including lead (Pb). Using data from 344 pairs of autism spectrum disorder (ASD) cases and age- and sex-matched typically developing (TD) controls (2–8 years old) from Jamaica, we investigated the interaction between three GST genes and ASD status as determinants of blood Pb concentrations (BPbCs). We found that ASD cases had lower geometric mean BPbCs than TD children (1.74 vs. 2.27 µg/dL, p < 0.01). Using a co-dominant genetic model, ASD cases with the Ile/Val genotype for the GSTP1 Ile105Val polymorphism had lower GM BPbCs than TD controls, after adjusting for a known interaction between GSTP1 and GSTT1, child’s parish, socioeconomic status, consumption of lettuce, fried plantains, and canned fish (Ile/Val: 1.78 vs. 2.13 µg/dL, p = 0.03). Similarly, among carriers of the I/I or I/D (I*) genotype for GSTT1 and GSTM1, ASD cases had lower adjusted GM BPbCs than TD controls (GSTT1 I*: 1.61 vs. 1.91 µg/dL, p = 0.01; GSTM1 I*: 1.71 vs. 2.04 µg/dL, p = 0.01). Our findings suggest that genetic polymorphisms in GST genes may influence detoxification of Pb by the enzymes they encode in Jamaican children with and without ASD.

BACKGROUND: Exposure to many environmental chemicals, including metals, often does not occur in isolation, hence requires assessment of the associations between exposure to mixtures of chemicals and human health. OBJECTIVES: To investigate associations of a metal mixture of lead (Pb), mercury (Hg), arsenic (As), cadmium (Cd), manganese (Mn), and aluminum (Al) in children with autism spectrum disorder (ASD), additively or interactively with each of three glutathione S-transferase (GST) genes (GSTP1, GSTT1, and GSTM1). METHOD: Using data from 266 case-control pairs of Jamaican children (2–8 years old), we fitted negative and positive generalized weighted quantile sum (gWQS) regression models to assess the aforementioned associations. RESULTS: Based on additive and interactive negative gWQS models adjusted for maternal age, parental education, child’s parish, and seafood consumption, we found inverse associations of the overall mixture score with ASD [MOR (95% CI): 0.70 (0.49, 0.99); P < 0.05) and [MOR (95%CI): 0.46 (0.25, 0.84); P = 0.01], respectively. In an unadjusted negative gWQS model, we found a marginally significant interaction between GSTP1 and a mixture of three metals (Pb, Hg, and Mn) (P = 0.07) while the association was no longer significant after adjustment for the same covariates (P = 0.24). CONCLUSIONS: Differences in diet between ASD and control groups may play a role in the inverse associations we found. The possible interactive association between Mn and GSTP1 in ASD based on gWQS is consistent with our previous reports. However, possible interaction of GSTP1 with Pb and Hg in ASD requires further investigation and replication.

Variation in levels of the human metabolome reflect changes in homeostasis, providing a window into health and disease. The genetic impact on circulating metabolites in Hispanics, a population with high cardiometabolic disease burden, is largely unknown. We conducted genome-wide association analyses on 640 circulating metabolites in 3,926 Hispanic Community Health Study/Study of Latinos participants. The estimated heritability for 640 metabolites ranged between 0%–54% with a median at 2.5%. We discovered 46 variant-metabolite pairs (p value < 1.2 × 10(−10), minor allele frequency ≥ 1%, proportion of variance explained [PEV] mean = 3.4%, PEV(range) = 1%–22%) with generalized effects in two population-based studies and confirmed 301 known locus-metabolite associations. Half of the identified variants with generalized effect were located in genes, including five nonsynonymous variants. We identified co-localization with the expression quantitative trait loci at 105 discovered and 151 known loci-metabolites sets. rs5855544, upstream of SLC51A, was associated with higher levels of three steroid sulfates and co-localized with expression levels of SLC51A in several tissues. Mendelian randomization (MR) analysis identified several metabolites associated with coronary heart disease (CHD) and type 2 diabetes. For example, two variants located in or near CYP4F2 (rs2108622 and rs79400241, respectively), involved in vitamin E metabolism, were associated with the levels of octadecanedioate and vitamin E metabolites (gamma-CEHC and gamma-CEHC glucuronide); MR analysis showed that genetically high levels of these metabolites were associated with lower odds of CHD. Our findings document the genetic architecture of circulating metabolites in an underrepresented Hispanic/Latino community, shedding light on disease etiology.

Background Mitochondrial DNA copy number (mtDNA-CN) has been associated with a variety of aging-related diseases, including all-cause mortality. However, the mechanism by which mtDNA-CN influences disease is not currently understood. One such mechanism may be through regulation of nuclear gene expression via the modification of nuclear DNA (nDNA) methylation. Methods To investigate this hypothesis, we assessed the relationship between mtDNA-CN and nDNA methylation in 2507 African American (AA) and European American (EA) participants from the Atherosclerosis Risk in Communities (ARIC) study. To validate our findings, we assayed an additional 2528 participants from the Cardiovascular Health Study (CHS) (N = 533) and Framingham Heart Study (FHS) (N = 1995). We further assessed the effect of experimental modification of mtDNA-CN through knockout of TFAM, a regulator of mtDNA replication, via CRISPR-Cas9. Results Thirty-four independent CpGs were associated with mtDNA-CN at genome-wide significance (P − 8). Meta-analysis across all cohorts identified six mtDNA-CN-associated CpGs at genome-wide significance (P − 8). Additionally, over half of these CpGs were associated with phenotypes known to be associated with mtDNA-CN, including coronary heart disease, cardiovascular disease, and mortality. Experimental modification of mtDNA-CN demonstrated that modulation of mtDNA-CN results in changes in nDNA methylation and gene expression of specific CpGs and nearby transcripts. Strikingly, the “neuroactive ligand receptor interaction” KEGG pathway was found to be highly overrepresented in the ARIC cohort (P = 5.24 × 10− 12), as well as the TFAM knockout methylation (P = 4.41 × 10− 4) and expression (P = 4.30 × 10− 4) studies. Conclusions These results demonstrate that changes in mtDNA-CN influence nDNA methylation at specific loci and result in differential expression of specific genes that may impact human health and disease via altered cell signaling.

Genotype-phenotype association studies often combine phenotype data from multiple studies to increase power. Harmonization of the data usually requires substantial effort due to heterogeneity in phenotype definitions, study design, data collection procedures, and data set organization. Here we describe a centralized system for phenotype harmonization that includes input from phenotype domain and study experts, quality control, documentation, reproducible results, and data sharing mechanisms. This system was developed for the National Heart, Lung and Blood Institute’s Trans-Omics for Precision Medicine (TOPMed) program, which is generating genomic and other omics data for >80 studies with extensive phenotype data. To date, 63 phenotypes have been harmonized across thousands of participants from up to 17 TOPMed studies per phenotype. We discuss the challenges faced in this undertaking and how they were addressed. The harmonized phenotype data and associated documentation have been submitted to National Institutes of Health data repositories for controlled-access by the scientific community. We also provide materials to facilitate future harmonization efforts by the community, which include (1) the code used to generate the 63 harmonized phenotypes, enabling others to reproduce, modify or extend these harmonizations to additional studies; and (2) results of labeling thousands of phenotype variables with controlled vocabulary terms.

BackgroundThe association between mitochondrial DNA-copy number (mtDNA-CN) and incident heart failure (HF) in the general population is unclear.MethodsWe examined the association between mtDNA-CN and the risk of incident HF among 10,802 participants free of HF at baseline from the Atherosclerosis Risk in Communities (ARIC) study, a large bi-racial population-based cohort. mtDNA-CN was estimated using probe intensities on the Affymetrix Genome-Wide Human single nucleotide polymorphisms Array 6.0. Incident HF events were identified through hospital discharge codes from 1987 until 2005 and through adjudication by the ARIC HF Classification Committee since 2005.ResultsDuring a median follow-up of 23.1 years, there were 2,227 incident HF events (incidence rate 10.3 per 1000 person-years). In fully adjusted models, the hazard ratios (95% confidence intervals) for HF comparing the 2nd through 5th quintiles of mtDNA-CN to the 1st quintile were 0.91 (0.80–1.04), 0.82 (0.72–0.93), 0.81 (0.71–0.92), and 0.74 (0.65–0.85), respectively (P for trend < 0.001). In stratified analyses, the associations between mtDNA-CN and HF were similar across examined subgroups. The inverse association between mtDNA-CN and incident HF was stronger in HF with reduced ejection fraction (HFrEF) than in HF with preserved ejection fraction (HFpEF).ConclusionsIn this prospective cohort, mtDNA-CN was inversely associated with the risk of incident HF suggesting that reduced levels of mtDNA-CN, a biomarker of mitochondrial dysfunction, could reflect early susceptibility to HF.

BACKGROUND: Genetic or pharmacological inhibition of de novo sphingolipid synthases prevented diabetes in animal studies. OBJECTIVES: We sought to evaluate prospective associations of serum sphingolipids with incident diabetes in a population-based cohort. METHODS: We included 2010 participants of the Hispanic Community Health Study/Study of Latinos (HCHS/SOL) aged 18–74 y who were free of diabetes and other major chronic diseases at baseline (2008–2011). Metabolomic profiling of fasting serum was performed using a global, untargeted approach. A total of 43 sphingolipids were quantified and, considering subclasses and chemical structures of individual species, 6 sphingolipid scores were constructed. Diabetes status was assessed using standard procedures including blood tests. Multivariable survey Poisson regressions were applied to estimate RR and 95% CI of incident diabetes associated with individual sphingolipids or sphingolipid scores. RESULTS: There were 224 incident cases of diabetes identified during, on average, 6 y of follow-up. After adjustment for socioeconomic and lifestyle factors, a ceramide score (RR (Q4 versus Q1) = 2.40; 95% CI: 1.24, 4.65; P-trend = 0.003) and a score of sphingomyelins with fully saturated sphingoid-fatty acid pairs (RR (Q4 versus Q1) = 3.15; 95% CI: 1.75, 5.67; P-trend

AimsWe tested the hypothesis that mitochondrial DNA copy number (CN) is associated with cardiometabolic disease (CMD) traits.Methods and resultsWe determined the cross-sectional association of mtDNA CN measured in whole blood with several CMD traits in 65,996 individuals (mean age 60, 54% women, and 79% European descent). Cohort- and ancestry/ethnicity-specific association analysis was performed adjusting for trait- and cohort-specific covariates. Age was slightly positively associated with age (0.03 s.d. / 10 years (95% CI=0.01, 0.05)) before 65 years, while every 10 years older age was associated with 0.14 s.d. lower level of mtDNA CN after 65 years (95% CI= -0.18, -0.10). In meta-analysis without adjustment for white blood cell (WBC) and differential count in participants of European descent (N=52,491), low mtDNA CN was associated with increased odds of obesity (OR with 95% CI=1.13 (1.11, 1.16), P=3.3e-30) and hypertension (OR=1.05 (1.03, 1.08), P=4.0e-07). Further adjusting for WBC and differential count in the same participants of European descent (N=44,035), associations became non-significant (P>0.05) for hypertension, attenuated for obesity (ORwithout cell count=1.15 (1.12, 1.18) versus ORcell count=1.06 (1.03, 1.08)) but strengthened for hyperlipidemia (ORwithout cell counts=1.03 (1.00, 1.06) versus ORcell counts=1.06 (1.03, 1.09)). The magnitude and directionality of most associations were consistent between participants of European descent and other ethnicity/ancestry origins.ConclusionLow levels of mtDNA CN in peripheral blood were associated with an increased risk of CMD diseases.