A loss-of-function CCR2 variant is associated with lower cardiovascular risk

Background and Aims

Ample evidence links CCL2, a key chemokine governing monocyte trafficking, with atherosclerosis. However, it remains unknown whether targeting the CCL2 receptor CCR2 could provide protection against cardiovascular disease.

Methods

Computationally predicted damaging (REVEL>0.5) variants within CCR2 were detected in whole-exome-sequencing data from 454,775 UK Biobank participants and tested for association with cardiovascular endpoints in gene-burden tests. Given the key role of CCR2 in monocyte mobilization, variants associated with lower monocyte count were prioritized for experimental validation. The response to CCL2 of human cells transfected with these variants was tested in migration and cAMP assays. Validated loss-of-function variants were tested for association with cardiovascular endpoints, atherosclerosis burden, and vascular risk factors. Significant associations were replicated in six independent datasets (n=1,062,595).

Results

Carriers of 45 predicted damaging CCR2 variants were at lower risk of myocardial infarction and coronary artery disease. One of these variants (M249K) was associated with lower monocyte count and decreased signaling and chemoattraction in response to CCL2. While M249K showed no association with conventional vascular risk factors, it was consistently associated with a lower risk of myocardial infarction (Odds Ratio: 0.66 95% Confidence Interval: 0.54-0.81,p=6.1×10-5) and coronary artery disease(Odds Ratio: 0.74 95% Confidence Interval: 0.62-0.87, p=2.9×10-4) in the UK Biobank and in six replication cohorts. In a phenome-wide association study, there was no evidence of higher infections risk among M249K carriers.

Conclusions

Carriers of an experimentally confirmed loss-of-function CCR2 variant are at a lower lifetime risk of myocardial infarction and coronary artery disease without carrying a higher infection risk. Our findings provide genetic support for the translational potential of CCR2-targeting as an atheroprotective approach.

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.

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.

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.

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.

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.

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.

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

Simple Summary

Over the last two decades of biological research, our understanding of how genes determine dental development and variation has expanded greatly. Here, we explore how this new knowledge can be applied to the fossil record of cercopithecid monkeys. We compare a traditional paleontological method for assessing dental size variation with measurement approaches derived from quantitative genetics and developmental biology. We find that these new methods for assessing dental variation provide novel insight to the evolution of the cercopithecid monkey dentition, different from the insight provided by traditional size measurements. When we explore the variation of these traits in the cercopithecid fossil record, we find that the variation is outside the range predicted based on extant variation alone. Our 21st century biological approach to paleontology reveals that we have even more to learn from fossils than previously recognized.

Abstract

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.

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.

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 (h²=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.

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.

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 (h² = 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 (V_P ) and additive genetic (V_A ) 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 V_P . EID V_A 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 V_A is shared with body mass variation, different methods of correcting for allometry by statistically controlling for size can alter residual V_P 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.

Genetic correlation of the plasma lipidome with type 2 diabetes, prediabetes and insulin resistance in Mexican American families

Background

Differential plasma concentrations of circulating lipid species are associated with pathogenesis of type 2 diabetes (T2D). Whether the wide inter-individual variability in the plasma lipidome contributes to the genetic basis of T2D is unknown. Here, we investigated the potential overlap in the genetic basis of the plasma lipidome and T2D-related traits.

Results

We used plasma lipidomic data (1202 pedigreed individuals, 319 lipid species representing 23 lipid classes) from San Antonio Family Heart Study in Mexican Americans. Bivariate trait analyses were used to estimate the genetic and environmental correlation of all lipid species with three T2D-related traits: risk of T2D, presence of prediabetes and homeostatic model of assessment – insulin resistance. We found that 44 lipid species were significantly genetically correlated with one or more of the three T2D-related traits. Majority of these lipid species belonged to the diacylglycerol (DAG, 17 species) and triacylglycerol (TAG, 17 species) classes. Six lipid species (all belonging to the triacylglycerol class and containing palmitate at the first position) were significantly genetically correlated with all the T2D-related traits.

Conclusions

Our results imply that: a) not all plasma lipid species are genetically informative for T2D pathogenesis; b) the DAG and TAG lipid classes partially share genetic basis of T2D; and c) 1-palmitate containing TAGs may provide additional insights into the genetic basis of T2D.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-017-0515-5) contains supplementary material, which is available to authorized users.

Lipidomic risk score independently and cost-effectively predicts risk of future type 2 diabetes: results from diverse cohorts

Background

Detection of type 2 diabetes (T2D) is routinely based on the presence of dysglycemia. Although disturbed lipid metabolism is a hallmark of T2D, the potential of plasma lipidomics as a biomarker of future T2D is unknown. Our objective was to develop and validate a plasma lipidomic risk score (LRS) as a biomarker of future type 2 diabetes and to evaluate its cost-effectiveness for T2D screening.

Methods

Plasma LRS, based on significantly associated lipid species from an array of 319 lipid species, was developed in a cohort of initially T2D-free individuals from the San Antonio Family Heart Study (SAFHS). The LRS derived from SAFHS as well as its recalibrated version were validated in an independent cohort from Australia – the AusDiab cohort. The participants were T2D-free at baseline and followed for 9197 person-years in the SAFHS cohort (n = 771) and 5930 person-years in the AusDiab cohort (n = 644). Statistically and clinically improved T2D prediction was evaluated with established statistical parameters in both cohorts. Modeling studies were conducted to determine whether the use of LRS would be cost-effective for T2D screening. The main outcome measures included accuracy and incremental value of the LRS over routinely used clinical predictors of T2D risk; validation of these results in an independent cohort and cost-effectiveness of including LRS in screening/intervention programs for T2D.

Results

The LRS was based on plasma concentration of dihydroceramide 18:0, lysoalkylphosphatidylcholine 22:1 and triacyglycerol 16:0/18:0/18:1. The score predicted future T2D independently of prediabetes with an accuracy of 76 %. Even in the subset of initially euglycemic individuals, the LRS improved T2D prediction. In the AusDiab cohort, the LRS continued to predict T2D significantly and independently. When combined with risk-stratification methods currently used in clinical practice, the LRS significantly improved the model fit (p < 0.001), information content (p < 0.001), discrimination (p < 0.001) and reclassification (p < 0.001) in both cohorts. Modeling studies demonstrated that LRS-based risk-stratification combined with metformin supplementation for high-risk individuals was the most cost-effective strategy for T2D prevention.

Conclusions

Considering the novelty, incremental value and cost-effectiveness of LRS it should be used for risk-stratification of future T2D.

Electronic supplementary material

The online version of this article (doi:10.1186/s12944-016-0234-3) contains supplementary material, which is available to authorized users.

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 (h² = 0.50 ± 0.05, p < 4 × 10⁻³⁵), and linkage analysis of serum uric acid concentrations confirmed our previous finding of a novel locus on 3p26 (LOD = 4.9, p < 1 × 10⁻⁵) in the extended sample. Additionally, we observed strong association of serum uric acid concentrations with variants in following candidate genes in the 3p26 region; inositol 1,4,5-trisphosphate receptor, type 1 (ITPR1), contactin 4 (CNTN4), decapping mRNA 1A (DCP1A); transglutaminase 4 (TGM4) and rho guanine nucleotide exchange factor (GEF) 26 (ARHGEF26) [p < 3 × 10⁻⁷; minor allele frequencies ranged between 0.003 and 0.42] and evidence of cis-regulation for ITPR1 transcripts.

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.

Soluble Forms of Intercellular and Vascular Cell Adhesion Molecules Independently Predict Progression to Type 2 Diabetes in Mexican American Families

Objective

While the role of type 2 diabetes (T2D) in inducing endothelial dysfunction is fairly well-established the etiological role of endothelial dysfunction in the onset of T2D is still a matter of debate. In the light of conflicting evidence in this regard, we conducted a prospective study to determine the association of circulating levels of soluble intercellular adhesion molecule 1 (sICAM-1) and soluble vessel cell adhesion molecule 1 (sVCAM-1) with incident T2D.

Methods

Data from this study came from 1,269 Mexican Americans of whom 821 initially T2D-free individuals were longitudinally followed up in the San Antonio Family Heart Study. These individuals were followed for 9752.95 person-years for development of T2D. Prospective association of sICAM-1 and sVCAM-1 with incident T2D was studied using Kaplan-Meier survival plots and mixed effects Cox proportional hazards modeling to account for relatedness among study participants. Incremental value of adhesion molecule biomarkers was studied using integrated discrimination improvement (IDI) and net reclassification improvement (NRI) indexes.

Results

Decreasing median values for serum concentrations of sICAM-1 and sVCAM-1 were observed in the following groups in this order: individuals with T2D at baseline, individuals who developed T2D during follow-up, individuals with prediabetes at baseline and normal glucose tolerant (NGT) individuals who remained T2D-free during follow-up. Top quartiles for sICAM-1 and sVCAM-1 were strongly and significantly associated with homeostatic model of assessment—insulin resistance (HOMA-IR). Mixed effects Cox proportional hazards modeling revealed that after correcting for important clinical confounders, high sICAM-1 and sVCAM-1 concentrations were associated with 2.52 and 1.99 times faster progression to T2D as compared to low concentrations, respectively. Individuals with high concentrations for both sICAM-1 and sVCAM-1 progressed to T2D 3.42 times faster than those with low values for both sICAM-1 and sVCAM-1. The results were similar in women in reproductive age group and the remainder of the cohort. Inclusion of sICAM-1 and sVCAM-1 in predictive models significantly improved reclassification and discrimination. The majority of these results were seen even when the analyses were restricted to NGT individuals.

Conclusion

Serum concentrations of sICAM-1 and sVCAM-1 independently and additively predict future T2D and represent important candidate biomarkers of T2D.

Genome- and epigenome-wide association study of hypertriglyceridemic waist in Mexican American families

BACKGROUND

There is growing interest in the hypertriglyceridemic waist (HTGW) phenotype, defined as high waist circumference (≥95 cm in males and ≥80 cm in females) combined with high serum triglyceride concentration (≥2.0 mmol/L in males and ≥1.5 mmol/L in females) as a marker of type 2 diabetes (T2D) and cardiovascular disease. However, the prevalence of this phenotype in high-risk populations, its association with T2D, and the genetic or epigenetic influences on HTGW are not well explored. Using data from large, extended families of Mexican Americans (a high-risk minority population in the USA) we aimed to: (1) estimate the prevalence of this phenotype, (2) test its association with T2D and related traits, and (3) dissect out the genetic and epigenetic associations with this phenotype using genome-wide and epigenome-wide studies, respectively.

RESULTS

Data for this study was from 850 Mexican American participants (representing 39 families) recruited under the ongoing San Antonio Family Heart Study, 26 % of these individuals had HTGW. This phenotype was significantly heritable (h (2) r = 0.52, p = 1.1 × 10(-5)) and independently associated with T2D as well as fasting glucose levels and insulin resistance. We conducted genome-wide association analyses using 759,809 single nucleotide polymorphisms (SNPs) and epigenome-wide association analyses using 457,331 CpG sites. There was no evidence of any SNP associated with HTGW at the genome-wide level but two CpG sites (cg00574958 and cg17058475) in CPT1A and one CpG site (cg06500161) in ABCG1 were significantly associated with HTGW and remained significant after adjusting for the closely related components of metabolic syndrome. CPT1A holds a cardinal position in the metabolism of long-chain fatty acids while ABCG1 plays a role in triglyceride metabolism.

CONCLUSIONS

Our results reemphasize the value of HTGW as a marker of T2D. This phenotype shows association with DNA methylation within CPT1A and ABCG1, genes involved in fatty acid and triglyceride metabolism. Our results underscore the importance of epigenetics in a clinically informative phenotype.

Lack of Association between SLC30A8 Variants and Type 2 Diabetes in Mexican American Families

SLC30A8 encodes zinc transporter 8 which is involved in packaging and release of insulin. Evidence for the association of SLC30A8 variants with type 2 diabetes (T2D) is inconclusive. We interrogated single nucleotide polymorphisms (SNPs) around SLC30A8 for association with T2D in high-risk, pedigreed individuals from extended Mexican American families. This study of 118 SNPs within 50 kb of the SLC30A8 locus tested the association with eight T2D-related traits at four levels: (i) each SNP using measured genotype approach (MGA); (ii) interaction of SNPs with age and sex; (iii) combinations of SNPs using Bayesian Quantitative Trait Nucleotide (BQTN) analyses; and (iv) entire gene locus using the gene burden test. Only one SNP (rs7817754) was significantly associated with incident T2D but a summary statistic based on all T2D-related traits identified 11 novel SNPs. Three SNPs and one SNP were weakly but interactively associated with age and sex, respectively. BQTN analyses could not demonstrate any informative combination of SNPs over MGA. Lastly, gene burden test results showed that at best the SLC30A8 locus could account for only 1-2% of the variability in T2D-related traits. Our results indicate a lack of association of the SLC30A8 SNPs with T2D in Mexican American families.

Axial quantitative ultrasound assessment of pediatric bone quality in eastern Nepal

This study presents quantitative ultrasonography (QUS) bone quality data for an underrepresented, south Asian pediatric population from Nepal. Data were collected as part of a longitudinal study of growth and development. This study offers normative data and documents the effect of stunting, wasting, and underweight on the bone properties measured by QUS.

INTRODUCTION

The purpose of this study was to (1) examine the bone quality of a rural, non-Western pediatric population using QUS, (2) explore variation in the trajectory of bone quality development between males and females, and (3) examine the impact of growth disruption(s) on bone quality.

METHODS

A cross-sectional study of 860 children and adolescents aged 5-18 years from the Jirel ethnic group in eastern Nepal was performed. The Sunlight Omnisense 7000P was used to assess bone quality of the distal 1/3 radius and midshaft tibia. WHO reference standards were used to assess growth disruptions of height, weight, and BMI.

RESULTS

QUS bone quality data for an underrepresented, non-Western pediatric population are presented for the radius and tibia. A sizable portion of the study participants were classified as stunted, wasted, and/or underweight. Despite this prevalence of growth disruption in the study sample, bone quality data conform to other documented populations with less growth disruption. Thus, this study offers normative data and documents the minimal effect of stunting, wasting, and underweight on the bone properties measured by QUS.

CONCLUSIONS

Non-Western pediatric populations are significantly underserved with regard to simple, non-invasive screening tools that may help identify developmental disorders and assess bone health. The children and adolescents examined here represent normal growth and development for an underrepresented south Asian population. While this work demonstrates that stunting, wasting, or underweight status at time of QUS assessment is not associated with poor bone quality, we do suggest that further study is needed to examine possible cumulative effects of persistent disruptions that may lead to compromised bone quality in later adolescence.

Novel epigenetic determinants of type 2 diabetes in Mexican-American families

Although DNA methylation is now recognized as an important mediator of complex diseases, the extent to which the genetic basis of such diseases is accounted for by DNA methylation is unknown. In the setting of large, extended families representing a minority, high-risk population of the USA, we aimed to characterize the role of epigenome-wide DNA methylation in type 2 diabetes (T2D). Using Illumina HumanMethylation450 BeadChip arrays, we tested for association of DNA methylation at 446 356 sites with age, sex and phenotypic traits related to T2D in 850 pedigreed Mexican-American individuals. Robust statistical analyses showed that (i) 15% of the methylome is significantly heritable, with a median heritability of 0.14; (ii) DNA methylation at 14% of CpG sites is associated with nearby sequence variants; (iii) 22% and 3% of the autosomal CpG sites are associated with age and sex, respectively; (iv) 53 CpG sites were significantly associated with liability to T2D, fasting blood glucose and insulin resistance; (v) DNA methylation levels at five CpG sites, mapping to three well-characterized genes (TXNIP, ABCG1 and SAMD12) independently explained 7.8% of the heritability of T2D (vi) methylation at these five sites was unlikely to be influenced by neighboring DNA sequence variation. Our study has identified novel epigenetic indicators of T2D risk in Mexican Americans who have increased risk for this disease. These results provide new insights into potential treatment targets of T2D.

Human plasma lipidome is pleiotropically associated with cardiovascular risk factors and death

BACKGROUND

Cardiovascular disease (CVD) is the most common cause of death in the United States and is associated with a high economic burden. Prevention of CVD focuses on controlling or improving the lipid profile of patients at risk. The human lipidome is made up of thousands of ubiquitous lipid species. By studying biologically simple canonical lipid species, we investigated whether the lipidome is genetically redundant and whether its genetic influences can provide clinically relevant clues of CVD risk.

METHODS AND RESULTS

We performed a genetic study of the human lipidome in 1212 individuals from 42 extended Mexican American families. High-throughput mass spectrometry enabled rapid capture of precise lipidomic profiles, providing 319 unique species. Using variance component-based heritability analyses and bivariate trait analyses, we detected significant genetic influences on each lipid assayed. Median heritability of the plasma lipid species was 0.37. Hierarchical clustering based on complex genetic correlation patterns identified 12 genetic clusters that characterized the plasma lipidome. These genetic clusters were differentially but consistently associated with risk factors of CVD, including central obesity, obesity, type 2 diabetes mellitus, raised serum triglycerides, and metabolic syndrome. Also, these clusters consistently predicted occurrence of cardiovascular deaths during follow-up.

CONCLUSIONS

The human plasma lipidome is heritable. Shared genetic influences reduce the dimensionality of the human lipidome into clusters that are associated with risk factors of CVD.

Baboons as a model to study genetics and epigenetics of human disease

A major challenge for understanding susceptibility to common human diseases is determining genetic and environmental factors that influence mechanisms underlying variation in disease-related traits. The most common diseases afflicting the US population are complex diseases that develop as a result of defects in multiple genetically controlled systems in response to environmental challenges. Unraveling the etiology of these diseases is exceedingly difficult because of the many genetic and environmental factors involved. Studies of complex disease genetics in humans are challenging because it is not possible to control pedigree structure and often not practical to control environmental conditions over an extended period of time. Furthermore, access to tissues relevant to many diseases from healthy individuals is quite limited. The baboon is a well-established research model for the study of a wide array of common complex diseases, including dyslipidemia, hypertension, obesity, and osteoporosis. It is possible to acquire tissues from healthy, genetically characterized baboons that have been exposed to defined environmental stimuli. In this review, we describe the genetic and physiologic similarity of baboons with humans, the ability and usefulness of controlling environment and breeding, and current genetic and genomic resources. We discuss studies on genetics of heart disease, obesity, diabetes, metabolic syndrome, hypertension, osteoporosis, osteoarthritis, and intrauterine growth restriction using the baboon as a model for human disease. We also summarize new studies and resources under development, providing examples of potential translational studies for targeted interventions and therapies for human disease.

Increased waist circumference is independently associated with hypothyroidism in Mexican Americans: replicative evidence from two large, population-based studies

BACKGROUND

Mexican Americans are at an increased risk of both thyroid dysfunction and metabolic syndrome (MS). Thus it is conceivable that some components of the MS may be associated with the risk of thyroid dysfunction in these individuals. Our objective was to investigate and replicate the potential association of MS traits with thyroid dysfunction in Mexican Americans.

METHODS

We conducted association testing for 18 MS traits in two large studies on Mexican Americans - the San Antonio Family Heart Study (SAFHS) and the National Health and Nutrition Examination Survey (NHANES) 2007-10. A total of 907 participants from 42 families in SAFHS and 1633 unrelated participants from NHANES 2007-10 were included in this study. The outcome measures were prevalence of clinical and subclinical hypothyroidism and thyroid function index (TFI) - a measure of thyroid function. For the SAFHS, we used polygenic regression analyses with multiple covariates to test associations in setting of family studies. For the NHANES 2007-10, we corrected for the survey design variables as needed for association analyses in survey data. In both datasets, we corrected for age, sex and their linear and quadratic interactions.

RESULTS

TFI was an accurate indicator of clinical thyroid status (area under the receiver-operating-characteristic curve to detect clinical hypothyroidism, 0.98) in both SAFHS and NHANES 2007-10. Of the 18 MS traits, waist circumference (WC) showed the most consistent association with TFI in both studies independently of age, sex and body mass index (BMI). In the SAFHS and NHANES 2007-10 datasets, each standard deviation increase in WC was associated with 0.13 (p < 0.001) and 0.11 (p < 0.001) unit increase in the TFI, respectively. In a series of polygenic and linear regression models, central obesity (defined as WC ≥ 102 cm in men and ≥88 cm in women) was associated with clinical and subclinical hypothyroidism independent of age, sex, BMI and type 2 diabetes in both datasets. Estimated prevalence of hypothyroidism was consistently high in those with central obesity, especially below 45y of age.

CONCLUSIONS

WC independently associates with increased risk of thyroid dysfunction. Use of WC to identify Mexican American subjects at high risk of thyroid dysfunction should be investigated in future studies.

Characterization of complex, co-adapted skeletal biomechanics phenotypes: a needed paradigm shift in the genetics of bone structure and function

The genetic architecture of skeletal biomechanical performance has tremendous potential to advance our knowledge of the biological mechanisms that drive variation in skeletal fragility and osteoporosis risk. Research using traditional approaches that focus on specific gene pathways is increasing our understanding of how and to what degree those pathways may affect population-level variation in fracture susceptibility, and shows that known pathways may affect bone fragility through unsuspected mechanisms. Non-traditional approaches that incorporate a new appreciation for the degree to which bone traits co-adapt to functional loading environments, using a wide variety of redundant compensatory mechanisms to meet both physiological and mechanical demands, represent a radical departure from the dominant reductionist paradigm and have the potential to rapidly advance our understanding of bone fragility and identification of new targets for therapeutic intervention.

Plasma lipidome is independently associated with variability in metabolic syndrome in Mexican American families

Plasma lipidome is now increasingly recognized as a potentially important marker of chronic diseases, but the exact extent of its contribution to the interindividual phenotypic variability in family studies is unknown. Here, we used the rich data from the ongoing San Antonio Family Heart Study (SAFHS) and developed a novel statistical approach to quantify the independent and additive value of the plasma lipidome in explaining metabolic syndrome (MS) variability in Mexican American families recruited in the SAFHS. Our analytical approach included two preprocessing steps: principal components analysis of the high-resolution plasma lipidomics data and construction of a subject-subject lipidomic similarity matrix. We then used the Sequential Oligogenic Linkage Analysis Routines software to model the complex family relationships, lipidomic similarities, and other important covariates in a variance components framework. Our results suggested that even after accounting for the shared genetic influences, indicators of lipemic status (total serum cholesterol, TGs, and HDL cholesterol), and obesity, the plasma lipidome independently explained 22% of variability in the homeostatic model of assessment-insulin resistance trait and 16% to 22% variability in glucose, insulin, and waist circumference. Our results demonstrate that plasma lipidomic studies can additively contribute to an understanding of the interindividual variability in MS.

Plasma dihydroceramide species associate with waist circumference in Mexican American families

OBJECTIVE

Waist circumference (WC), the clinical marker of central obesity, is gaining popularity as a screening tool for type 2 diabetes (T2D). While there is epidemiologic evidence favoring the WC-T2D association, its biological substantiation is generally weak. Our objective was to determine the independent association of plasma lipid repertoire with WC.

METHODS

Samples and data from the San Antonio Family Heart Study of 1208 Mexican Americans from 42 extended families were used. Association of plasma lipidomic profiles with the cross-sectionally assessed WC was determined. Plasma lipidomic profiling entailed liquid chromatography with mass spectrometry. Statistical analyses included multivariable polygenic regression models and bivariate trait analyses using the SOLAR software.

RESULTS

After adjusting for age and sex interactions, body mass index, homeostasis model of assessment-insulin resistance, total cholesterol, triglycerides, high density lipoproteins and use of lipid lowering drugs, dihydroceramides as a class were associated with WC. Dihydroceramide species 18:0, 20:0, 22:0, and 24:1 were significantly associated and genetically correlated with WC. Two sphingomyelin species (31:1 and 41:1) were also associated with WC.

CONCLUSIONS

Plasma dihydroceramide levels independently associate with WC. Thus, high resolution plasma lipidomic studies can provide further credence to the biological underpinnings of the association of WC with T2D.

Waist circumference is genetically correlated with incident Type 2 diabetes in Mexican-American families

AIMS

We aimed to determine the genetic and environmental correlation between various anthropometric indexes and incident Type 2 diabetes with a focus on waist circumference.

METHODS

We used the data on extended Mexican-American families (808 subjects, 7617.92 person-years follow-up) from the San Antonio Family Heart Study and estimated the genetic and environmental correlations of 16 anthropometric indexes with the genetic liability of incident Type 2 diabetes. We performed bivariate trait analyses using the solar software package.

RESULTS

All 16 anthropometric indexes were significantly heritable (range of heritabilities 0.24-0.99). Thirteen indexes were found to have significant environmental correlation with the liability of incident Type 2 diabetes. In contrast, only anthropometric indexes consisting of waist circumference (waist circumference, waist-hip ratio and waist-height ratio) were significantly genetically correlated (genetic correlation coefficients: 0.45, 0.55 and 0.44, respectively) with the liability of incident Type 2 diabetes. We did not observe such a correlation for BMI.

CONCLUSIONS

Waist circumference as a predictor of future Type 2 diabetes is supported by the finding that they share common genetic influences.

Genome-wide association analysis confirms and extends the association of SLC2A9 with serum uric acid levels to Mexican Americans

Increased serum uric acid (SUA) is a risk factor for gout and renal and cardiovascular disease (CVD). The purpose of this study was to identify genetic factors that affect the variation in SUA in 632 Mexican Americans participants of the San Antonio Family Heart Study (SAFHS). A genome-wide association (GWA) analysis was performed using the Illumina Human Hap 550K single nucleotide polymorphism (SNP) microarray. We used a linear regression-based association test under an additive model of allelic effect, while accounting for non-independence among family members via a kinship variance component. All analyses were performed in the software package SOLAR. SNPs rs6832439, rs13131257, and rs737267 in solute carrier protein 2 family, member 9 (SLC2A9) were associated with SUA at genome-wide significance (p < 1.3 × 10⁻⁷). The minor alleles of these SNPs had frequencies of 36.2, 36.2, and 38.2%, respectively, and were associated with decreasing SUA levels. All of these SNPs were located in introns 3–7 of SLC2A9, the location of the previously reported associations in European populations. When analyzed for association with cardiovascular-renal disease risk factors, conditional on SLC2A9 SNPs strongly associated with SUA, significant associations were found for SLC2A9 SNPs with BMI, body weight, and waist circumference (p < 1.4 × 10⁻³) and suggestive associations with albumin-creatinine ratio and total antioxidant status (TAS). The SLC2A9 gene encodes an urate transporter that has considerable influence on variation in SUA. In addition to the primary association locus, suggestive evidence (p < 1.9 × 10⁻⁶) for joint linkage/association (JLA) was found at a previously-reported urate quantitative trait locus (Logarithm of odds score = 3.6) on 3p26.3. In summary, our GWAS extends and confirms the association of SLC2A9 with SUA for the first time in a Mexican American cohort and also shows for the first time its association with cardiovascular-renal disease risk factors.

Intracortical bone remodeling variation shows strong genetic effects

Intracortical microstructure influences crack propagation and arrest within bone cortex. Genetic variation in intracortical remodeling may contribute to mechanical integrity and, therefore, fracture risk. Our aim was to determine the degree to which normal population-level variation in intracortical microstructure is due to genetic variation. We examined right femurs from 101 baboons (74 females, 27 males; aged 7-33 years) from a single, extended pedigree to determine osteon number, osteon area (On.Ar), haversian canal area, osteon population density, percent osteonal bone (%On.B), wall thickness (W.Th), and cortical porosity (Ct.Po). Through evaluation of the covariance in intracortical properties between pairs of relatives, we quantified the contribution of additive genetic effects (heritability [h (2)]) to variation in these traits using a variance decomposition approach. Significant age and sex effects account for 9 % (Ct.Po) to 21 % (W.Th) of intracortical microstructural variation. After accounting for age and sex, significant genetic effects are evident for On.Ar (h (2) = 0.79, p = 0.002), %On.B (h (2) = 0.82, p = 0.003), and W.Th (h (2) = 0.61, p = 0.013), indicating that 61-82 % of the residual variation (after accounting for age and sex effects) is due to additive genetic effects. This corresponds to 48-75 % of the total phenotypic variance. Our results demonstrate that normal, population-level variation in cortical microstructure is significantly influenced by genes. As a critical mediator of crack behavior in bone cortex, intracortical microstructural variation provides another mechanism through which genetic variation may affect fracture risk.

Linkage of type 2 diabetes on chromosome 9p24 in Mexican Americans: additional evidence from the Veterans Administration Genetic Epidemiology Study (VAGES)

OBJECTIVE

Type 2 diabetes (T2DM) is a complex metabolic disease and is more prevalent in certain ethnic groups such as the Mexican Americans. The goal of our study was to perform a genome-wide linkage (GWL) analysis to localize T2DM susceptibility loci in Mexican Americans.

METHODS

We used the phenotypic and genotypic data from 1,122 Mexican-American individuals (307 families) who participated in the Veterans Administration Genetic Epidemiology Study (VAGES). GWL analysis was performed using the variance components approach. Data from 2 additional Mexican-American family studies, the San Antonio Family Heart Study (SAFHS) and the San Antonio Family Diabetes/Gallbladder Study (SAFDGS), were combined with the VAGES data to test for improved linkage evidence.

RESULTS

After adjusting for covariate effects, T2DM was found to be under significant genetic influences (h2 = 0.62, p = 2.7 × 10(-6)). The strongest evidence for linkage of T2DM occurred between markers D9S1871 and D9S2169 on chromosome 9p24.2-p24.1 (LOD = 1.8). Given that we previously reported suggestive evidence for linkage of T2DM at this region also in SAFDGS, we found the significant and increased linkage evidence (LOD = 4.3, empirical p = 1.0 × 10(-5), genome-wide p = 1.6 × 10(-3)) for T2DM at the same chromosomal region, when we performed a GWL analysis of the VAGES data combined with the SAFHS and SAFDGS data.

CONCLUSION

Significant T2DM linkage evidence was found on chromosome 9p24 in Mexican Americans. Importantly, the chromosomal region of interest in this study overlaps with several recent genome-wide association studies involving T2DM-related traits. Given its overlap with such findings and our own initial T2DM association findings in the 9p24 chromosomal region, high throughput sequencing of the linked chromosomal region could identify the potential causal T2DM genes.

Plasma lipid profiling in a large population-based cohort

We have performed plasma lipid profiling using liquid chromatography electrospray ionization tandem mass spectrometry on a population cohort of more than 1,000 individuals. From 10 μl of plasma we were able to acquire comparative measures of 312 lipids across 23 lipid classes and subclasses including sphingolipids, phospholipids, glycerolipids, and cholesterol esters (CEs) in 20 min. Using linear and logistic regression, we identified statistically significant associations of lipid classes, subclasses, and individual lipid species with anthropometric and physiological measures. In addition to the expected associations of CEs and triacylglycerol with age, sex, and body mass index (BMI), ceramide was significantly higher in males and was independently associated with age and BMI. Associations were also observed for sphingomyelin with age but this lipid subclass was lower in males. Lysophospholipids were associated with age and higher in males, but showed a strong negative association with BMI. Many of these lipids have previously been associated with chronic diseases including cardiovascular disease and may mediate the interactions of age, sex, and obesity with disease risk.

Plasma lipidomic profile signature of hypertension in Mexican American families: specific role of diacylglycerols

Both as a component of metabolic syndrome and as an independent entity, hypertension poses a continued challenge with regard to its diagnosis, pathogenesis, and treatment. Previous studies have documented connections between hypertension and indicators of lipid metabolism. Novel technologies, such as plasma lipidomic profiling, promise a better understanding of disorders in which there is a derangement of the lipid metabolism. However, association of plasma lipidomic profiles with hypertension in a high-risk population, such as Mexican Americans, has not been evaluated before. Using the rich data and sample resource from the ongoing San Antonio Family Heart Study, we conducted plasma lipidomic profiling by combining high-performance liquid chromatography with tandem mass spectroscopy to characterize 319 lipid species in 1192 individuals from 42 large and extended Mexican American families. Robust statistical analyses using polygenic regression models, liability threshold models, and bivariate trait analyses implemented in the SOLAR software were conducted after accounting for obesity, insulin resistance, and relative abundance of various lipoprotein fractions. Diacylglycerols, in general, and the DG 16:0/22:5 and DG 16:0/22:6 lipid species, in particular, were significantly associated with systolic blood pressure (SBP), diastolic blood pressure (DBP), and mean arterial pressure (MAP), as well as liability of incident hypertension measured during 7140.17 person-years of follow-up. Four lipid species, including the DG 16:0/22:5 and DG 16:0/22:6 species, showed significant genetic correlations with the liability of hypertension in bivariate trait analyses. Our results demonstrate the value of plasma lipidomic profiling in the context of hypertension and identify disturbance of diacylglycerol metabolism as an independent biomarker of hypertension.

Identification of candidate genes encoding an LDL-C QTL in baboons

Cardiovascular disease (CVD) is the leading cause of death in developed countries, and dyslipidemia is a major risk factor for CVD. We previously identified a cluster of quantitative trait loci (QTL) on baboon chromosome 11 for multiple, related quantitative traits for serum LDL-cholesterol (LDL-C). Here we report differentially regulated hepatic genes encoding an LDL-C QTL that influences LDL-C levels in baboons. We performed hepatic whole-genome expression profiling for LDL-C-discordant baboons fed a high-cholesterol, high-fat (HCHF) diet for seven weeks. We detected expression of 117 genes within the QTL 2-LOD support interval. Three genes were differentially expressed in low LDL-C responders and 8 in high LDL-C responders in response to a HCHF diet. Seven genes (ACVR1B, CALCOCO1, DGKA, ERBB3, KRT73, MYL6B, TENC1) showed discordant expression between low and high LDL-C responders. To prioritize candidate genes, we integrated miRNA and mRNA expression profiles using network tools and found that four candidates (ACVR1B, DGKA, ERBB3, TENC1) were miRNA targets and that the miRNAs were inversely expressed to the target genes. Candidate gene expression was validated using QRT-PCR and Western blotting. This study reveals candidate genes that influence variation in LDL-C in baboons and potential genetic mechanisms for further investigation.

Variability in associations of phosphatidylcholine molecular species with metabolic syndrome in Mexican-American families

Plasma lipidomic studies using high performance liquid chromatography and mass spectroscopy offer detailed insights into metabolic processes. Taking the example of the most abundant plasma lipid class (phosphatidylcholines) we used the rich phenotypic and lipidomic data from the ongoing San Antonio Family Heart Study of large extended Mexican-American families to assess the variability of association of the plasma phosphatidylcholine species with metabolic syndrome. Using robust statistical analytical methods, our study made two important observations. First, there was a wide variability in the association of phosphatidylcholine species with risk measures of metabolic syndrome. Phosphatidylcholine 40:7 was associated with a low risk while phosphatidylcholines 32:1 and 38:3 were associated with a high risk of metabolic syndrome. Second, all the odd chain phosphatidylcholines were associated with a reduced risk of metabolic syndrome implying that phosphatidylcholines derived from dairy products might be beneficial against metabolic syndrome. Our results demonstrate the value of lipid species-specific information provided by the upcoming array of lipidomic studies and open potential avenues for prevention and control of metabolic syndrome in high prevalence settings.

Genetic architecture of carotid artery intima-media thickness in Mexican Americans

BACKGROUND- Intima-media thickness (IMT) of the common and internal carotid arteries is an established surrogate for atherosclerosis and predicts risk of stroke and myocardial infarction. Often IMT is measured as the average of these 2 arteries; yet, they are believed to result from separate biological mechanisms. The aim of this study was to conduct a family-based genome-wide association study (GWAS) for IMT to identify polymorphisms influencing IMT and to determine if distinct carotid artery segments are influenced by different genetic components. METHODS AND RESULTS- IMT for the common and internal carotid arteries was determined through B-mode ultrasound in 772 Mexican Americans from the San Antonio Family Heart Study. A GWAS using 931219 single-nucleotide polymorphisms was undertaken with 6 internal and common carotid artery IMT phenotypes using an additive measured genotype model. The most robust association detected was for 2 single-nucleotide polymorphisms (rs16983261, rs6113474; P=1.60e(-7)) in complete linkage disequilibrium on chromosome 20p11 for the internal carotid artery near wall, next to the gene PAX1. We also replicated previously reported GWAS regions on chromosomes 19q13 and 7q22. We found no overlapping associations between internal and common carotid artery phenotypes at P<5.0e(-6). The genetic correlation between the 2 carotid IMT arterial segments was 0.51. CONCLUSIONS- This study represents the first large-scale GWAS of carotid IMT in a non-European population and identified several novel loci. We do not detect any shared GWAS signals between common and internal carotid arterial segments, but the moderate genetic correlation implies both common and unique genetic components.

Genotype×age interaction in human transcriptional ageing

Individual differences in biological ageing (i.e., the rate of physiological response to the passage of time) may be due in part to genotype-specific variation in gene action. However, the sources of heritable variation in human age-related gene expression profiles are largely unknown. We have profiled genome-wide expression in peripheral blood mononuclear cells from 1240 individuals in large families and found 4472 human autosomal transcripts, representing ~4349 genes, significantly correlated with age. We identified 623 transcripts that show genotype by age interaction in addition to a main effect of age, defining a large set of novel candidates for characterization of the mechanisms of differential biological ageing. We applied a novel SNP genotype × age interaction test to one of these candidates, the ubiquilin-like gene UBQLNL, and found evidence of joint cis-association and genotype by age interaction as well as trans-genotype by age interaction for UBQLNL expression. Both UBQLNL expression levels at recruitment and cis genotype are associated with longitudinal cancer risk in our study cohort.

Modularity in the mammalian dentition: mice and monkeys share a common dental genetic architecture

The concept of modularity provides a useful tool for exploring the relationship between genotype and phenotype. Here, we use quantitative genetics to identify modularity within the mammalian dentition, connecting the genetics of organogenesis to the genetics of population-level variation for a phenotype well represented in the fossil record. We estimated the correlations between dental traits owing to the shared additive effects of genes (pleiotropy) and compared the pleiotropic relationships among homologous traits in two evolutionary distant taxa-mice and baboons. We find that in both mice and baboons, who shared a common ancestor >65 Ma, incisor size variation is genetically independent of molar size variation. Furthermore, baboon premolars show independent genetic variation from incisors, suggesting that a modular genetic architecture separates incisors from these posterior teeth as well. Such genetic independence between modules provides an explanation for the extensive diversity of incisor size variation seen throughout mammalian evolution-variation uncorrelated with equivalent levels of postcanine tooth size variation. The modularity identified here is supported by the odontogenic homeobox code proposed for the patterning of the rodent dentition. The baboon postcanine pattern of incomplete pleiotropy is also consistent with predictions from the morphogenetic field model.

A genome-wide linkage scan for quantitative trait loci influencing the craniofacial complex in humans (Homo sapiens sapiens)

The genetic architecture of the craniofacial complex has been the subject of intense scrutiny because of the high frequency of congenital malformations. Numerous animal models have been used to document the early development of the craniofacial complex, but few studies have focused directly on the genetic underpinnings of normal variation in the human craniofacial complex. This study examines 80 quantitative traits derived from lateral cephalographs of 981 participants in the Fels Longitudinal Study, Wright State University, Dayton, Ohio. Quantitative genetic analyses were conducted using the Sequential Oligogenic Linkage Analysis Routines analytic platform, a maximum-likelihood variance components method that incorporates all familial information for parameter estimation. Heritability estimates were significant and of moderate to high magnitude for all craniofacial traits. Additionally, significant quantitative trait loci (QTL) were identified for 10 traits from the three developmental components (basicranium, splanchnocranium, and neurocranium) of the craniofacial complex. These QTL were found on chromosomes 3, 6, 11, 12, and 14. This study of the genetic architecture of the craniofacial complex elucidates fundamental information of the genetic architecture of the craniofacial complex in humans.

Effects of diet on genetic regulation of lipoprotein metabolism in baboons

Several measures of lipoprotein phenotype are significant predictors of cardiovascular risk. Although such lipoprotein phenotypes are under strong genetic control, it is not clear to what extent they are controlled by the same - and by different - genes and whether these relationships may be altered in different dietary environments. Therefore, we measured six lipoprotein traits (three LDL traits - LDLC and apoB concentrations and LDL size - and three HDL traits - HDLC and apoA1 concentrations and HDL size) on each of three diets differing in level of fat and cholesterol. In bivariate analyses, all but two metabolically related trait pairs were genetically correlated, though none were completely correlated, implying additive genetic effects by both pleiotropic and unique genes. In comparing genetic correlations for the same pair of traits across diet, we detected evidence of diet effects on genetic control of these metabolically related traits; specifically, increasing level of dietary cholesterol was associated with a significant decrease in the genetic correlation of apoA1 with HDL size, and a significant increase in the genetic correlations of LDL size with LDLC and apoB. The results suggest a complex network of genes affecting lipoprotein metabolism: the genes may exert both unique and pleiotropic effects; the genes may exert detectable effects in many or only in specific dietary environments.

Genetic integration of molar cusp size variation in baboons

Many studies of primate diversity and evolution rely on dental morphology for insight into diet, behavior, and phylogenetic relationships. Consequently, variation in molar cusp size has increasingly become a phenotype of interest. In 2007 we published a quantitative genetic analysis of mandibular molar cusp size variation in baboons. Those results provided more questions than answers, as the pattern of genetic integration did not fit predictions from odontogenesis. To follow up, we expanded our study to include data from the maxillary molar cusps. Here we report on these later analyses, as well as inter-arch comparisons with the mandibular data. We analyzed variation in two-dimensional maxillary molar cusp size using data collected from a captive pedigreed breeding colony of baboons, Papio hamadryas, housed at the Southwest National Primate Research Center. These analyses show that variation in maxillary molar cusp size is heritable and sexually dimorphic. We also estimated additive genetic correlations between cusps on the same crown, homologous cusps along the tooth row, and maxillary and mandibular cusps. The pattern for maxillary molars yields genetic correlations of one between the paracone-metacone and protocone-hypocone. Bivariate analyses of cuspal homologues on adjacent teeth yield correlations that are high or not significantly different from one. Between dental arcades, the nonoccluding cusps consistently yield high genetic correlations, especially the metaconid-paracone and metaconid-metacone. This pattern of genetic correlation does not immediately accord with the pattern of development and/or calcification, however these results do follow predictions that can be made from the evolutionary history of the tribosphenic molar.

The etiology of eruption disorders - further evidence of a 'genetic paradigm'

The clinical spectrum of tooth eruption disorders includes both syndromic and non-syndromic problems ranging from delayed eruption to a complete failure of eruption. A defect in the differential apposition/resorption mechanism in alveolar bone can cause conditions such as tooth ankylosis, primary failure of eruption, failure of eruption due to inadequate arch length and canine impaction. As our knowledge of the molecular events underlying normal tooth eruption has increased, so too has our understanding of clinical eruption disorders. The recent finding that one gene, parathyroid hormone receptor 1 (PTH1R), is causative for familial cases of primary failure of eruption (PFE) suggests that other disturbances in tooth eruption may have a genetic etiology. In this report, we evaluated the current terminology (ankylosis, PFE, secondary retention, etc.) used to describe non-syndromic eruption disorders, in light of this genetic discovery. We observed that some individuals previously diagnosed with ankylosis were subsequently found to have alterations in the PTH1R gene, indicating the initial misdiagnosis of ankylosis and the necessary re-classification of PFE. We further investigated the relationship of the PTH1R gene, using a network pathway analysis, to determine its connectivity to previously identified genes that are critical to normal tooth eruption. We found that PTH1R acts in a pathway with genes such as PTHrP that have been shown to be important in bone remodeling, hence eruption, in a rat model. Thus, recent advances in our understanding of normal and abnormal tooth eruption should allow us in the future to develop a clinical nomenclature system based more on the molecular genetic cause of the eruption failures versus the clinical appearance of the various eruption disorders.