Metabolic syndrome traits exhibit genotype-by-environment interaction in relation to socioeconomic status in the Mexican American family heart study.

Background: Socioeconomic Status (SES) is a potent environmental determinant of health. To our knowledge, no assessment of genotype-environment interaction has been conducted to consider the joint effects of socioeconomic status and genetics on risk for metabolic disease. We analyzed data from the Mexican American Family Studies (MAFS) to evaluate the hypothesis that genotype-by-environment interaction (GxE) is an essential determinant of variation in risk factors for metabolic syndrome (MS). Methods: We employed a maximum likelihood estimation of the decomposition of variance components to detect GxE interaction. After excluding individuals with diabetes and individuals on medication for diabetes, hypertension, or dyslipidemia, we analyzed 12 MS risk factors: fasting glucose (FG), fasting insulin (FI), 2-h glucose (2G), 2-h insulin (2I), body mass index (BMI), waist circumference (WC), leptin (LP), high-density lipoprotein-cholesterol (HDL-C), triglycerides (TG), total serum cholesterol (TSC), systolic blood pressure (SBP), and diastolic blood pressure (DBP). Our SES variable used a combined score of Duncan's socioeconomic index and education years. Heterogeneity in the additive genetic variance across the SES continuum and a departure from unity in the genetic correlation coefficient were taken as evidence of GxE interaction. Hypothesis tests were conducted using standard likelihood ratio tests. Results: We found evidence of GxE for fasting glucose, 2-h glucose, 2-h insulin, BMI, and triglycerides. The genetic effects underlying the insulin/glucose metabolism component of MS are upregulated at the lower end of the SES spectrum. We also determined that the household variance for systolic blood pressure decreased with increasing SES. Conclusion: These results show a significant change in the GxE interaction underlying the major components of MS in response to changes in socioeconomic status. Further mRNA sequencing studies will identify genes and canonical gene pathways to support our molecular-level hypotheses.

A statistical framework for powerful multi-trait rare variant analysis in large-scale whole-genome sequencing studies.

Large-scale whole-genome sequencing (WGS) studies have improved our understanding of the contributions of coding and noncoding rare variants to complex human traits. Leveraging association effect sizes across multiple traits in WGS rare variant association analysis can improve statistical power over single-trait analysis, and also detect pleiotropic genes and regions. Existing multi-trait methods have limited ability to perform rare variant analysis of large-scale WGS data. We propose MultiSTAAR, a statistical framework and computationally-scalable analytical pipeline for functionally-informed multi-trait rare variant analysis in large-scale WGS studies. MultiSTAAR accounts for relatedness, population structure and correlation among phenotypes by jointly analyzing multiple traits, and further empowers rare variant association analysis by incorporating multiple functional annotations. We applied MultiSTAAR to jointly analyze three lipid traits (low-density lipoprotein cholesterol, high-density lipoprotein cholesterol and triglycerides) in 61,861 multi-ethnic samples from the Trans-Omics for Precision Medicine (TOPMed) Program. We discovered new associations with lipid traits missed by single-trait analysis, including rare variants within an enhancer of NIPSNAP3A and an intergenic region on chromosome 1.

Type 2 Diabetes Modifies the Association of CAD Genomic Risk Variants With Subclinical Atherosclerosis.

BACKGROUND: Individuals with type 2 diabetes (T2D) have an increased risk of coronary artery disease (CAD), but questions remain about the underlying pathology. Identifying which CAD loci are modified by T2D in the development of subclinical atherosclerosis (coronary artery calcification [CAC], carotid intima-media thickness, or carotid plaque) may improve our understanding of the mechanisms leading to the increased CAD in T2D. METHODS: We compared the common and rare variant associations of known CAD loci from the literature on CAC, carotid intima-media thickness, and carotid plaque in up to 29 670 participants, including up to 24 157 normoglycemic controls and 5513 T2D cases leveraging whole-genome sequencing data from the Trans-Omics for Precision Medicine program. We included first-order T2D interaction terms in each model to determine whether CAD loci were modified by T2D. The genetic main and interaction effects were assessed using a joint test to determine whether a CAD variant, or gene-based rare variant set, was associated with the respective subclinical atherosclerosis measures and then further determined whether these loci had a significant interaction test. RESULTS: Using a Bonferroni-corrected significance threshold of P<1.6×10-4, we identified 3 genes (ATP1B1, ARVCF, and LIPG) associated with CAC and 2 genes (ABCG8 and EIF2B2) associated with carotid intima-media thickness and carotid plaque, respectively, through gene-based rare variant set analysis. Both ATP1B1 and ARVCF also had significantly different associations for CAC in T2D cases versus controls. No significant interaction tests were identified through the candidate single-variant analysis. CONCLUSIONS: These results highlight T2D as an important modifier of rare variant associations in CAD loci with CAC.

Rare variants in long non-coding RNAs are associated with blood lipid levels in the TOPMed whole-genome sequencing study.

Long non-coding RNAs (lncRNAs) are known to perform important regulatory functions in lipid metabolism. Large-scale whole-genome sequencing (WGS) studies and new statistical methods for variant set tests now provide an opportunity to assess more associations between rare variants in lncRNA genes and complex traits across the genome. In this study, we used high-coverage WGS from 66,329 participants of diverse ancestries with measurement of blood lipids and lipoproteins (LDL-C, HDL-C, TC, and TG) in the National Heart, Lung, and Blood Institute (NHLBI) Trans-Omics for Precision Medicine (TOPMed) program to investigate the role of lncRNAs in lipid variability. We aggregated rare variants for 165,375 lncRNA genes based on their genomic locations and conducted rare-variant aggregate association tests using the STAAR (variant-set test for association using annotation information) framework. We performed STAAR conditional analysis adjusting for common variants in known lipid GWAS loci and rare-coding variants in nearby protein-coding genes. Our analyses revealed 83 rare lncRNA variant sets significantly associated with blood lipid levels, all of which were located in known lipid GWAS loci (in a ±500-kb window of a Global Lipids Genetics Consortium index variant). Notably, 61 out of 83 signals (73%) were conditionally independent of common regulatory variation and rare protein-coding variation at the same loci. We replicated 34 out of 61 (56%) conditionally independent associations using the independent UK Biobank WGS data. Our results expand the genetic architecture of blood lipids to rare variants in lncRNAs.

Genotype-by-socioeconomic status interaction influences heart disease risk scores and carotid artery thickness in Mexican Americans: the predominant role of education in comparison to household income and socioeconomic index.

Background: Socioeconomic status (SES) is a potent environmental determinant of health. To our knowledge, no assessment of genotype-environment interaction has been conducted to consider the joint effects of socioeconomic status and genetics on risk for cardiovascular disease (CVD). We analyzed Mexican American Family Studies (MAFS) data to evaluate the hypothesis that genotype-by-environment interaction (GxE) is an important determinant of variation in CVD risk factors. Methods: We employed a linear mixed model to investigate GxE in Mexican American extended families. We studied two proxies for CVD [Pooled Cohort Equation Risk Scores/Framingham Risk Scores (FRS/PCRS) and carotid artery intima-media thickness (CA-IMT)] in relation to socioeconomic status as determined by Duncan's Socioeconomic Index (SEI), years of education, and household income. Results: We calculated heritability for FRS/PCRS and carotid artery intima-media thickness. There was evidence of GxE due to additive genetic variance heterogeneity and genetic correlation for FRS, PCRS, and CA-IMT measures for education (environment) but not for household income or SEI. Conclusion: The genetic effects underlying CVD are dynamically modulated at the lower end of the SES spectrum. There is a significant change in the genetic architecture underlying the major components of CVD in response to changes in education.

Carriers of rare damaging CCR2 genetic variants are at lower risk of atherosclerotic disease.

Background: The CCL2/CCR2 axis governs monocyte trafficking and recruitment to atherosclerotic lesions. Human genetic analyses and population-based studies support an association between circulating CCL2 levels and atherosclerosis. Still, it remains unknown whether pharmacological targeting of CCR2, the main CCL2 receptor, would provide protection against human atherosclerotic disease. Methods: In whole-exome sequencing data from 454,775 UK Biobank participants (40-69 years), we identified predicted loss-of-function (LoF) or damaging missense (REVEL score >0.5) variants within the CCR2 gene. We prioritized variants associated with lower monocyte count (p<0.05) and tested associations with vascular risk factors and risk of atherosclerotic disease over a mean follow-up of 14 years. The results were replicated in a pooled cohort of three independent datasets (TOPMed, deCODE and Penn Medicine BioBank; total n=441,445) and the effect of the most frequent damaging variant was experimentally validated. Results: A total of 45 predicted LoF or damaging missense variants were identified in the CCR2 gene, 4 of which were also significantly associated with lower monocyte count, but not with other white blood cell counts. Heterozygous carriers of these variants were at a lower risk of a combined atherosclerosis outcome, showed a lower burden of atherosclerosis across four vascular beds, and were at a lower lifetime risk of coronary artery disease and myocardial infarction. There was no evidence of association with vascular risk factors including LDL-cholesterol, blood pressure, glycemic status, or C-reactive protein. Using a cAMP assay, we found that cells transfected with the most frequent CCR2 damaging variant (3:46358273:T:A, M249K, 547 carriers, frequency: 0.14%) show a decrease in signaling in response to CCL2. The associations of the M249K variant with myocardial infarction were consistent across cohorts (ORUKB: 0.62 95%CI: 0.39-0.96; ORexternal: 0.64 95%CI: 0.34-1.19; ORpooled: 0.64 95%CI: 0.450.90). In a phenome-wide association study, we found no evidence for higher risk of common infections or mortality among carriers of damaging CCR2 variants. Conclusions: Heterozygous carriers of damaging CCR2 variants have a lower burden of atherosclerosis and lower lifetime risk of myocardial infarction. In conjunction with previous evidence from experimental and epidemiological studies, our findings highlight the translational potential of CCR2-targeting as an atheroprotective approach.

Cerebral small vessel disease burden is associated with decreased abundance of gut Barnesiella intestinihominis bacterium in the Framingham Heart Study.

A bidirectional communication exists between the brain and the gut, in which the gut microbiota influences cognitive function and vice-versa. Gut dysbiosis has been linked to several diseases, including Alzheimer's disease and related dementias (ADRD). However, the relationship between gut dysbiosis and markers of cerebral small vessel disease (cSVD), a major contributor to ADRD, is unknown. In this cross-sectional study, we examined the connection between the gut microbiome, cognitive, and neuroimaging markers of cSVD in the Framingham Heart Study (FHS). Markers of cSVD included white matter hyperintensities (WMH), peak width of skeletonized mean diffusivity (PSMD), and executive function (EF), estimated as the difference between the trail-making tests B and A. We included 972 FHS participants with MRI scans, neurocognitive measures, and stool samples and quantified the gut microbiota composition using 16S rRNA sequencing. We used multivariable association and differential abundance analyses adjusting for age, sex, BMI, and education level to estimate the association between gut microbiota and WMH, PSMD, and EF measures. Our results suggest an increased abundance of Pseudobutyrivibrio and Ruminococcus genera was associated with lower WMH and PSMD (p values < 0.001), as well as better executive function (p values < 0.01). In addition, in both differential and multivariable analyses, we found that the gram-negative bacterium Barnesiella intestinihominis was strongly associated with markers indicating a higher cSVD burden. Finally, functional analyses using PICRUSt implicated various KEGG pathways, including microbial quorum sensing, AMP/GMP-activated protein kinase, phenylpyruvate, and ß-hydroxybutyrate production previously associated with cognitive performance and dementia. Our study provides important insights into the association between the gut microbiome and cSVD, but further studies are needed to replicate the findings.

Rare variants in long non-coding RNAs are associated with blood lipid levels in the TOPMed Whole Genome Sequencing Study.

Long non-coding RNAs (lncRNAs) are known to perform important regulatory functions. Large-scale whole genome sequencing (WGS) studies and new statistical methods for variant set tests now provide an opportunity to assess the associations between rare variants in lncRNA genes and complex traits across the genome. In this study, we used high-coverage WGS from 66,329 participants of diverse ancestries with blood lipid levels (LDL-C, HDL-C, TC, and TG) in the National Heart, Lung, and Blood Institute (NHLBI) Trans-Omics for Precision Medicine (TOPMed) program to investigate the role of lncRNAs in lipid variability. We aggregated rare variants for 165,375 lncRNA genes based on their genomic locations and conducted rare variant aggregate association tests using the STAAR (variant-Set Test for Association using Annotation infoRmation) framework. We performed STAAR conditional analysis adjusting for common variants in known lipid GWAS loci and rare coding variants in nearby protein coding genes. Our analyses revealed 83 rare lncRNA variant sets significantly associated with blood lipid levels, all of which were located in known lipid GWAS loci (in a ±500 kb window of a Global Lipids Genetics Consortium index variant). Notably, 61 out of 83 signals (73%) were conditionally independent of common regulatory variations and rare protein coding variations at the same loci. We replicated 34 out of 61 (56%) conditionally independent associations using the independent UK Biobank WGS data. Our results expand the genetic architecture of blood lipids to rare variants in lncRNA, implicating new therapeutic opportunities.

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.

Keeping 21st Century Paleontology Grounded: Quantitative Genetic Analyses and Ancestral State Reconstruction Re-Emphasize the Essentiality of Fossils.

Advances in genetics and developmental biology are revealing the relationship between genotype and dental phenotype (G:P), providing new approaches for how paleontologists assess dental variation in the fossil record. Our aim was to understand how the method of trait definition influences the ability to reconstruct phylogenetic relationships and evolutionary history in the Cercopithecidae, the Linnaean Family of monkeys currently living in Africa and Asia. We compared the two-dimensional assessment of molar size (calculated as the mesiodistal length of the crown multiplied by the buccolingual breadth) to a trait that reflects developmental influences on molar development (the inhibitory cascade, IC) and two traits that reflect the genetic architecture of postcanine tooth size variation (defined through quantitative genetic analyses: MMC and PMM). All traits were significantly influenced by the additive effects of genes and had similarly high heritability estimates. The proportion of covariate effects was greater for two-dimensional size compared to the G:P-defined traits. IC and MMC both showed evidence of selection, suggesting that they result from the same genetic architecture. When compared to the fossil record, Ancestral State Reconstruction using extant taxa consistently underestimated MMC and PMM values, highlighting the necessity of fossil data for understanding evolutionary patterns in these traits. Given that G:P-defined dental traits may provide insight to biological mechanisms that reach far beyond the dentition, this new approach to fossil morphology has the potential to open an entirely new window onto extinct paleobiologies. Without the fossil record, we would not be able to grasp the full range of variation in those biological mechanisms that have existed throughout evolution.

Genetic influences on dentognathic morphology in the Jirel population of Nepal.

Patterns of genetic variation and covariation impact the evolution of the craniofacial complex and contribute to clinically significant malocclusions in modern human populations. Previous quantitative genetic studies have estimated the heritabilities and genetic correlations of skeletal and dental traits in humans and nonhuman primates, but none have estimated these quantitative genetic parameters across the dentognathic complex. A large and powerful pedigree from the Jirel population of Nepal was leveraged to estimate heritabilities and genetic correlations in 62 maxillary and mandibular arch dimensions, incisor and canine lengths, and post-canine tooth crown areas (N ≥ 739). Quantitative genetic parameter estimation was performed using maximum likelihood-based variance decomposition. Residual heritability estimates were significant for all traits, ranging from 0.269 to 0.898. Genetic correlations were positive for all trait pairs. Principal components analyses of the phenotypic and genetic correlation matrices indicate an overall size effect across all measurements on the first principal component. Additional principal components demonstrate positive relationships between post-canine tooth crown areas and arch lengths and negative relationships between post-canine tooth crown areas and arch widths, and between arch lengths and arch widths. Based on these findings, morphological variation in the human dentognathic complex may be constrained by genetic relationships between dental dimensions and arch lengths, with weaker genetic correlations between these traits and arch widths allowing for variation in arch shape. The patterns identified are expected to have impacted the evolution of the dentognathic complex and its genetic architecture as well as the prevalence of dental crowding in modern human populations.

Identifying the Lipidomic Effects of a Rare Loss-of-Function Deletion in ANGPTL3.

BACKGROUND: The identification and understanding of therapeutic targets for atherosclerotic cardiovascular disease is of fundamental importance given its global health and economic burden. Inhibition of ANGPTL3 (angiopoietin-like 3) has demonstrated a cardioprotective effect, showing promise for atherosclerotic cardiovascular disease treatment, and is currently the focus of ongoing clinical trials. Here, we assessed the genetic basis of variation in ANGPTL3 levels in the San Antonio Family Heart Study. METHODS: We assayed ANGPTL3 protein levels in ≈1000 Mexican Americans from extended pedigrees. By drawing upon existing plasma lipidome profiles and genomic data we conducted analyses to understand the genetic basis to variation in ANGPTL3 protein levels, and accordingly the correlation with the plasma lipidome. RESULTS: In a variance components framework, we identified that variation in ANGPTL3 was significantly heritable (h2=0.33, P=1.31×10-16). To explore the genetic basis of this heritability, we conducted a genome-wide linkage scan and identified significant linkage (logarithm of odds =6.18) to a locus on chromosome 1 at 90 centimorgans, corresponding to the ANGPTL3 gene location. In the genomes of 23 individuals from a single pedigree, we identified a loss-of-function variant, rs398122988 (N121Kfs*2), in ANGPTL3, that was significantly associated with lower ANGPTL3 levels (ß=-1.69 SD units, P=3.367×10-13), and accounted for the linkage signal at this locus. Given the known role of ANGPTL3 as an inhibitor of endothelial and lipoprotein lipase, we explored the association of ANGPTL3 protein levels and rs398122988 with the plasma lipidome and related phenotypes, identifying novel associations with phosphatidylinositols. CONCLUSIONS: Variation in ANGPTL3 protein levels is heritable and under significant genetic control. Both ANGPTL3 levels and loss-of-function variants in ANGPTL3 have significant associations with the plasma lipidome. These findings further our understanding of ANGPTL3 as a therapeutic target for atherosclerotic cardiovascular disease.

Dynamic incorporation of multiple in silico functional annotations empowers rare variant association analysis of large whole-genome sequencing studies at scale.

Large-scale whole-genome sequencing studies have enabled the analysis of rare variants (RVs) associated with complex phenotypes. Commonly used RV association tests have limited scope to leverage variant functions. We propose STAAR (variant-set test for association using annotation information), a scalable and powerful RV association test method that effectively incorporates both variant categories and multiple complementary annotations using a dynamic weighting scheme. For the latter, we introduce 'annotation principal components', multidimensional summaries of in silico variant annotations. STAAR accounts for population structure and relatedness and is scalable for analyzing very large cohort and biobank whole-genome sequencing studies of continuous and dichotomous traits. We applied STAAR to identify RVs associated with four lipid traits in 12,316 discovery and 17,822 replication samples from the Trans-Omics for Precision Medicine Program. We discovered and replicated new RV associations, including disruptive missense RVs of NPC1L1 and an intergenic region near APOC1P1 associated with low-density lipoprotein cholesterol.

Diet-induced leukocyte telomere shortening in a baboon model for early stage atherosclerosis.

Reported associations between leukocyte telomere length (LTL) attrition, diet and cardiovascular disease (CVD) are inconsistent. This study explores effects of prolonged exposure to a high cholesterol high fat (HCHF) diet on LTL in a baboon model of atherosclerosis. We measured LTL by qPCR in pedigreed baboons fed a chow (n = 105) or HCHF (n = 106) diet for 2 years, tested for effects of diet on LTL, and association between CVD risk factors and atherosclerotic lesions with LTL. Though not different at baseline, after 2 years median LTL is shorter in HCHF fed baboons (P < 0.0001). Diet predicts sex- and age-adjusted LTL and LTL attrition (P = 0.0009 and 0.0156, respectively). Serum concentrations of CVD biomarkers are associated with LTL at the 2-year endpoint and LTL accounts approximately 6% of the variance in aortic lesions (P = 0.04). Although heritable at baseline (h2 = 0.27, P = 0.027) and after 2 years (h2 = 0.46, P = 0.0038), baseline LTL does not predict lesion extent after 2 years. Atherogenic diet influences LTL, and LTL is a potential biomarker for early atherosclerosis. Prolonged exposure to an atherogenic diet decreases LTL and increases LTL attrition, and shortened LTL is associated with early-stage atherosclerosis in pedigreed baboons.

Rare DEGS1 variant significantly alters de novo ceramide synthesis pathway.

The de novo ceramide synthesis pathway is essential to human biology and health, but genetic influences remain unexplored. The core function of this pathway is the generation of biologically active ceramide from its precursor, dihydroceramide. Dihydroceramides have diverse, often protective, biological roles; conversely, increased ceramide levels are biomarkers of complex disease. To explore the genetics of the ceramide synthesis pathway, we searched for deleterious nonsynonymous variants in the genomes of 1,020 Mexican Americans from extended pedigrees. We identified a Hispanic ancestry-specific rare functional variant, L175Q, in delta 4-desaturase, sphingolipid 1 (DEGS1), a key enzyme in the pathway that converts dihydroceramide to ceramide. This amino acid change was significantly associated with large increases in plasma dihydroceramides. Indexes of DEGS1 enzymatic activity were dramatically reduced in heterozygotes. CRISPR/Cas9 genome editing of HepG2 cells confirmed that the L175Q variant results in a partial loss of function for the DEGS1 enzyme. Understanding the biological role of DEGS1 variants, such as L175Q, in ceramide synthesis may improve the understanding of metabolic-related disorders and spur ongoing research of drug targets along this pathway.

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

Disentangling the genetic overlap between cholesterol and suicide risk.

Suicide is major public health concern; one million individuals worldwide die by suicide each year of which there are many more attempts. Thus, it is imperative that robust and reliable indicators, or biomarkers, of suicide risk be identified so that individuals at risk can be identified and provided appropriate interventions as quickly as possible. Previous work has revealed a relationship between low levels of circulating cholesterol and suicide risk, implicating cholesterol level as one such potential biomarker, but the factors underlying this relationship remain unknown. In the present study, we applied a combination of bivariate polygenic and coefficient-of-relatedness analysis, followed by mediation analysis, in a large sample of Mexican-American individuals from extended pedigrees [N = 1897; 96 pedigrees (average size = 19.17 individuals, range = 2-189) 60% female; mean age = 42.58 years, range = 18-97 years, sd = 15.75 years] with no exclusion criteria for any given psychiatric disorder. We observed that total esterified cholesterol measured at the time of psychiatric assessment shared a significant genetic overlap with risk for suicide attempt (ρg = -0.64, p = 1.24 × 10-04). We also found that total unesterified cholesterol measured around 20 years prior to assessment varied as a function of genetic proximity to an affected individual (h2 = 0.21, se = 0.10, p = 8.73 × 10-04; ßsuicide = -0.70, se = 0.25, p = 8.90 × 10-03). Finally, we found that the relationship between total unesterified cholesterol and suicide risk was significantly mediated by ABCA-1-specific cholesterol efflux capacity (ßsuicide-efflux = -0.45, p = 0.039; ßefflux-cholexterol = -0.34, p < 0.0001; ßindirect = -0.15, p = 0.044). These findings suggest that the relatively well-delineated process of cholesterol metabolism and associated molecular pathways will be informative for understanding the neurobiological underpinnings of risk for suicide attempt.

Additive genetic variation in the craniofacial skeleton of baboons (genus Papio) and its relationship to body and cranial size.

OBJECTIVES: Determining the genetic architecture of quantitative traits and genetic correlations among them is important for understanding morphological evolution patterns. We address two questions regarding papionin evolution: (1) what effect do body and cranial size, age, and sex have on phenotypic (VP ) and additive genetic (VA ) variation in baboon crania, and (2) how might additive genetic correlations between craniofacial traits and body mass affect morphological evolution? MATERIALS AND METHODS: We use a large captive pedigreed baboon sample to estimate quantitative genetic parameters for craniofacial dimensions (EIDs). Our models include nested combinations of the covariates listed above. We also simulate the correlated response of a given EID due to selection on body mass alone. RESULTS: Covariates account for 1.2-91% of craniofacial VP . EID VA decreases across models as more covariates are included. The median genetic correlation estimate between each EID and body mass is 0.33. Analysis of the multivariate response to selection reveals that observed patterns of craniofacial variation in extant baboons cannot be attributed solely to correlated response to selection on body mass, particularly in males. DISCUSSION: Because a relatively large proportion of EID VA is shared with body mass variation, different methods of correcting for allometry by statistically controlling for size can alter residual VP patterns. This may conflate direct selection effects on craniofacial variation with those resulting from a correlated response to body mass selection. This shared genetic variation may partially explain how selection for increased body mass in two different papionin lineages produced remarkably similar craniofacial phenotypes.

Diet-induced early-stage atherosclerosis in baboons: Lipoproteins, atherogenesis, and arterial compliance.

BACKGROUND: The purpose of this study was to determine whether dietary manipulation can reliably induce early-stage atherosclerosis and clinically relevant changes in vascular function in an established, well-characterized non-human primate model. METHODS: We fed 112 baboons a high-cholesterol, high-fat challenge diet for two years. We assayed circulating biomarkers of cardiovascular disease (CVD) risk, at 0, 7, and 104 weeks into the challenge; assessed arterial compliance noninvasively at 104 weeks; and measured atherosclerotic lesions in three major arteries at necropsy. RESULTS: We observed evidence of atherosclerosis in all but one baboon fed the two-year challenge diet. CVD risk biomarkers, the prevalence, size, and complexity of arterial lesions, plus consequent arterial stiffness, were increased in comparison with dietary control animals. CONCLUSIONS: Feeding baboons a high-cholesterol, high-fat diet for two years reliably induces atherosclerosis, with risk factor profiles, arterial lesions, and changes in vascular function also seen in humans.

TRAK2, a novel regulator of ABCA1 expression, cholesterol efflux and HDL biogenesis.

Aims: The recent failures of HDL-raising therapies have underscored our incomplete understanding of HDL biology. Therefore there is an urgent need to comprehensively investigate HDL metabolism to enable the development of effective HDL-centric therapies. To identify novel regulators of HDL metabolism, we performed a joint analysis of human genetic, transcriptomic, and plasma HDL-cholesterol (HDL-C) concentration data and identified a novel association between trafficking protein, kinesin binding 2 (TRAK2) and HDL-C concentration. Here we characterize the molecular basis of the novel association between TRAK2 and HDL-cholesterol concentration. Methods and results: Analysis of lymphocyte transcriptomic data together with plasma HDL from the San Antonio Family Heart Study (n = 1240) revealed a significant negative correlation between TRAK2 mRNA levels and HDL-C concentration, HDL particle diameter and HDL subspecies heterogeneity. TRAK2 siRNA-mediated knockdown significantly increased cholesterol efflux to apolipoprotein A-I and isolated HDL from human macrophage (THP-1) and liver (HepG2) cells by increasing the mRNA and protein expression of the cholesterol transporter ATP-binding cassette, sub-family A member 1 (ABCA1). The effect of TRAK2 knockdown on cholesterol efflux was abolished in the absence of ABCA1, indicating that TRAK2 functions in an ABCA1-dependent efflux pathway. TRAK2 knockdown significantly increased liver X receptor (LXR) binding at the ABCA1 promoter, establishing TRAK2 as a regulator of LXR-mediated transcription of ABCA1. Conclusion: We show, for the first time, that TRAK2 is a novel regulator of LXR-mediated ABCA1 expression, cholesterol efflux, and HDL biogenesis. TRAK2 may therefore be an important target in the development of anti-atherosclerotic therapies.