Analysis of serum changes in response to a high fat high cholesterol diet challenge reveals metabolic biomarkers of atherosclerosis

Atherosclerotic plaques are characterized by an accumulation of macrophages, lipids, smooth muscle cells, and fibroblasts, and, in advanced stages, necrotic debris within the arterial walls. Dietary habits such as high fat and high cholesterol (HFHC) consumption are known risk factors for atherosclerosis. However, the key metabolic contributors to diet-induced atherosclerosis are far from established. Herein, we investigate the role of a 2-year HFHC diet challenge in the metabolic changes of development and progression of atherosclerosis. We used a non-human primate (NHP) model (baboons, n = 60) fed a HFHC diet for two years and compared metabolomic profiles in serum from animals on baseline chow with serum collected after the challenge diet using two-dimensional gas chromatography time-of-flight mass-spectrometry (2D GC-ToF-MS) for untargeted metabolomic analysis, to quantify metabolites that contribute to atherosclerotic lesion formation. Further, clinical biomarkers associated with atherosclerosis, lipoprotein measures, fat indices, and arterial plaque formation (lesions) were quantified. Using two chemical derivatization (i.e., silylation) approaches, we quantified 321 metabolites belonging to 66 different metabolic pathways, which revealed significantly different metabolic profiles of HFHC diet and chow diet fed baboon sera. We found heritability of two important metabolites, lactic acid and asparagine, in the context of diet-induced metabolic changes. In addition, abundance of cholesterol, lactic acid, and asparagine were sex-dependent. Finally, 35 metabolites correlated (R2, 0.068-0.271, P \u3c 0.05) with total lesion burden assessed in three arteries (aortic arch, common iliac artery, and descending aorta) which could serve as potential biomarkers pending further validation. This study demonstrates the feasibility of detecting sex-specific and heritable metabolites in NHPs with diet-induced atherosclerosis using untargeted metabolomics allowing understanding of atherosclerotic disease progression in humans

Blood-Based miRNA Biomarkers as Correlates of Brain-Based miRNA Expression

The use of easily accessible peripheral samples, such as blood or saliva, to investigate neurological and neuropsychiatric disorders is well-established in genetic and epigenetic research, but the pathological implications of such biomarkers are not easily discerned. To better understand the relationship between peripheral blood- and brain-based epigenetic activity, we conducted a pilot study on captive baboons (Papio hamadryas) to investigate correlations between miRNA expression in peripheral blood mononuclear cells (PBMCs) and 14 different cortical and subcortical brain regions, represented by two study groups comprised of 4 and 6 animals. Using next-generation sequencing, we identified 362 miRNAs expressed at ≥ 10 read counts in 80% or more of the brain samples analyzed. Nominally significant pairwise correlations (one-sided P \u3c 0.05) between peripheral blood and mean brain expression levels of individual miRNAs were observed for 39 and 44 miRNAs in each group. When miRNA expression levels were averaged for tissue type across animals within the groups, Spearman\u27s rank correlations between PBMCs and the brain regions are all highly significant (r s = 0.47-0.57; P \u3c 2.2 × 10-16), although pairwise correlations among the brain regions are markedly stronger (r s = 0.86-0.99). Principal component analysis revealed differentiation in miRNA expression between peripheral blood and the brain regions for the first component (accounting for ∼75% of variance). Linear mixed effects modeling attributed most of the variance in expression to differences between miRNAs (\u3e70%), with non-significant 7.5% and 13.1% assigned to differences between blood and brain-based samples in the two study groups. Hierarchical UPGMA clustering revealed a major co-expression branch in both study groups, comprised of miRNAs globally upregulated in blood relative to the brain samples, exhibiting an enrichment of miRNAs expressed in immune cells (CD14+, CD15+, CD19+, CD3+, and CD56 + leukocytes) among the top blood-brain correlates, with the gene MYC, encoding a master transcription factor that regulates angiogenesis and neural stem cell activation, representing the most prevalent miRNA target. Although some differentiation was observed between tissue types, these preliminary findings reveal wider correlated patterns between blood- and brain-expressed miRNAs, suggesting the potential utility of blood-based miRNA profiling for investigating by proxy certain miRNA activity in the brain, with implications for neuroinflammatory and c-Myc-mediated processes

Physical activity and FTO genotype by physical activity interactive influences on obesity

Background: Although the effect of the fat mass and obesity-associated (FTO) gene on adiposity is well established, there is a lack of evidence whether physical activity (PA) modifies the effect of FTO variants on obesity in Latino populations. Therefore, the purpose of this study was to examine PA influences and interactive effects between FTO variants and PA on measures of adiposity in Latinos. Results: After controlling for age and sex, participants who did not engage in regular PA exhibited higher BMI, fat mass, HC, and WC with statistical significance (P \u3c 0.001). Although significant associations between the three FTO genotypes and adiposity measures were found, none of the FTO genotype by PA interaction assessments revealed nominally significant associations. However, several of such interactive influences exhibited considerable trend towards association. Conclusions: These data suggest that adiposity measures are associated with PA and FTO variants in Latinos, but the impact of their interactive influences on these obesity measures appear to be minimal. Future studies with large sample sizes may help to determine whether individuals with specific FTO variants exhibit differential responses to PA interventions

Hepatic transcript signatures predict atherosclerotic lesion burden prior to a 2-year high cholesterol, high fat diet challenge

The purpose of this study was to identify molecular mechanisms by which the liver influences total lesion burden in a nonhuman primate model (NHP) of cardiovascular disease with acute and chronic feeding of a high cholesterol, high fat (HCHF) diet. Baboons (47 females, 64 males) were fed a HCHF diet for 2 years (y); liver biopsies were collected at baseline, 7 weeks (w) and 2y, and lesions were quantified in aortic arch, descending aorta, and common iliac at 2y. Unbiased weighted gene co-expression network analysis (WGCNA) revealed several modules of hepatic genes correlated with lesions at different time points of dietary challenge. Pathway and network analyses were performed to study the roles of hepatic module genes. More significant pathways were observed in males than females. In males, we found modules enriched for genes in oxidative phosphorylation at baseline, opioid signaling at 7w, and EIF2 signaling and HNF1A and HNF4A networks at baseline and 2y. One module enriched for fatty acid β oxidation pathway genes was found in males and females at 2y. To our knowledge, this is the first study of a large NHP cohort to identify hepatic genes that correlate with lesion burden. Correlations of baseline and 7w module genes with lesions at 2y were observed in males but not in females. Pathway analyses of baseline and 7w module genes indicate EIF2 signaling, oxidative phosphorylation, and μ-opioid signaling are possible mechanisms that predict lesion formation induced by HCHF diet consumption in males. Our findings of coordinated hepatic transcriptional response in male baboons but not female baboons indicate underlying molecular mechanisms differ between female and male primate atherosclerosis

The Arg59Trp variant in ANGPTL8 (betatrophin) is associated with total and HDL-cholesterol in American Indians and Mexican Americans and differentially affects cleavage of ANGPTL3

We previously identified a locus linked to total cholesterol (TC) concentration in Pima Indians on chromosome 19p. To characterize this locus, we genotyped \u3e2000 SNPs in 1838 Pimas and assessed association with log(TC). We observed evidence for association with log(TC) with rs2278426 (3.5% decrease/copy of the T allele; P=5.045×10(-6)) in the ANGPTL8 (angiopoietin-like 8) gene. We replicated this association in 2413 participants of the San Antonio Mexican American Family Study (SAMAFS: 2.0% decrease per copy of the T allele; P=0.005842). In a meta-analysis of the combined data, we found the strongest estimated effect with rs2278426 (P=2.563×10(-7)). The variant T allele at rs2278426 predicts an Arg59Trp substitution and has previously been associated with LDL-C and HDL-C. In Pimas and SAMAFS participants, the T allele of rs2278426 was associated with reduced HDL-C levels (P=0.000741 and 0.00002, respectively), and the combined estimated effect for the two cohorts was -3.8% (P=8.526×10(-8)). ANGPTL8 transcript and protein levels increased in response to both glucose and insulin. The variant allele was associated with increased levels of cleaved ANGPTL3. We conclude that individuals with the variant allele may have lower TC and HDL-C levels due to increased activation of ANGPTL3 by ANGPTL8

Physical activity and FTO genotype by physical activity interactive influences on obesity

Abstract Background Although the effect of the fat mass and obesity-associated (FTO) gene on adiposity is well established, there is a lack of evidence whether physical activity (PA) modifies the effect of FTO variants on obesity in Latino populations. Therefore, the purpose of this study was to examine PA influences and interactive effects between FTO variants and PA on measures of adiposity in Latinos. Results After controlling for age and sex, participants who did not engage in regular PA exhibited higher BMI, fat mass, HC, and WC with statistical significance (P < 0.001). Although significant associations between the three FTO genotypes and adiposity measures were found, none of the FTO genotype by PA interaction assessments revealed nominally significant associations. However, several of such interactive influences exhibited considerable trend towards association. Conclusions These data suggest that adiposity measures are associated with PA and FTO variants in Latinos, but the impact of their interactive influences on these obesity measures appear to be minimal. Future studies with large sample sizes may help to determine whether individuals with specific FTO variants exhibit differential responses to PA interventions

Data for Genetic Analysis Workshop 18: human whole genome sequence, blood pressure, and simulated phenotypes in extended pedigrees

Genetic Analysis Workshop 18 (GAW18) focused on identification of genes and functional variants that influence complex phenotypes in human sequence data. Data for the workshop were donated by the T2D-GENES Consortium and included whole genome sequences for odd-numbered autosomes in 464 key individuals selected from 20 Mexican American families, a dense set of single-nucleotide polymorphisms in 959 individuals in these families, and longitudinal data on systolic and diastolic blood pressure measured at 1-4 examinations over a period of 20 years. Simulated phenotypes were generated based on the real sequence data and pedigree structures. In the design of the simulation model, gene expression measures from the San Antonio Family Heart Study (not distributed as part of the GAW18 data) were used to identify genes whose mRNA levels were correlated with blood pressure. Observed variants within these genes were designated as functional in the GAW18 simulation if they were nonsynonymous and predicted to have deleterious effects on protein function or if they were noncoding and associated with mRNA levels. Two simulated longitudinal phenotypes were modeled to have the same trait distributions as the real systolic and diastolic blood pressure data, with effects of age, sex, and medication use, including a genotype-medication interaction. For each phenotype, more than 1000 sequence variants in more than 200 genes present on the odd-numbered autosomes individually explained less than 0.01-2.78% of phenotypic variance. Cumulatively, variants in the most influential gene explained 7.79% of trait variance. An additional simulated phenotype, Q1, was designed to be correlated among family members but to not be associated with any sequence variants. Two hundred replicates of the phenotypes were simulated, with each including data for 849 individuals

Genetic and environmental (physical fitness and sedentary activity) interaction effects on cardiometabolic risk factors in Mexican American children and adolescents

Knowledge on genetic and environmental (G × E) interaction effects on cardiometabolic risk factors (CMRFs) in children is limited. The purpose of this study was to examine the impact of G × E interaction effects on CMRFs in Mexican American (MA) children (n = 617, ages 6–17 years). The environments examined were sedentary activity (SA), assessed by recalls from “yesterday” (SAy) and “usually” (SAu) and physical fitness (PF) assessed by Harvard PF scores (HPFS). CMRF data included body mass index (BMI), waist circumference (WC), fat mass (FM), fasting insulin (FI), homeostasis model of assessment—insulin resistance (HOMA‐IR), high‐density lipoprotein cholesterol (HDL‐C), triglycerides (TG), systolic (SBP) and diastolic (DBP) blood pressure, and number of metabolic syndrome components (MSC). We examined potential G × E interaction in the phenotypic expression of CMRFs using variance component models and likelihood‐based statistical inference. Significant G × SA interactions were identified for six CMRFs: BMI, WC, FI, HOMA‐IR, MSC, and HDL, and significant G × HPFS interactions were observed for four CMRFs: BMI, WC, FM, and HOMA‐IR. However, after correcting for multiple hypothesis testing, only WC × SAy, FM × SAy, and FI × SAu interactions became marginally significant. After correcting for multiple testing, most of CMRFs exhibited significant G × E interactions (Reduced G × E model vs. Constrained model). These findings provide evidence that genetic factors interact with SA and PF to influence variation in CMRFs, and underscore the need for better understanding of these relationships to develop strategies and interventions to effectively reduce or prevent cardiometabolic risk in children

GWAS and transcriptional analysis prioritize ITPR1 and CNTN4 for a serum uric acid 3p26 QTL in Mexican Americans

Background: The variation in serum uric acid concentrations is under significant genetic influence. Elevated SUA concentrations have been linked to increased risk for gout, kidney stones, chronic kidney disease, and cardiovascular disease whereas reduced serum uric acid concentrations have been linked to multiple sclerosis, Parkinson’s disease and Alzheimer’s disease. Previously, we identified a novel locus on chromosome 3p26 affecting serum uric acid concentrations in Mexican Americans from San Antonio Family Heart Study. As a follow up, we examined genome-wide single nucleotide polymorphism data in an extended cohort of 1281 Mexican Americans from multigenerational families of the San Antonio Family Heart Study and the San Antonio Family Diabetes/ Gallbladder Study. We used a linear regression-based joint linkage/association test under an additive model of allelic effect, while accounting for non-independence among family members via a kinship variance component. Results:Univariate genetic analysis indicated serum uric acid concentrations to be significant heritable (h2 = 0.50 ± 0.05, p Conclusion: Our results confirm the importance of the chromosome 3p26 locus and genetic variants in this region in the regulation of serum uric acid concentrations