Significant genotype by diet (G × D) interaction effects on cardiometabolic responses to a pedigree-wide, dietary challenge in vervet monkeys (Chlorocebus aethiops sabaeus)

Integrated omics analysis reveals sirtuin signaling is central to hepatic response to a high fructose diet

BACKGROUND

Dietary high fructose (HFr) is a known metabolic disruptor contributing to development of obesity and diabetes in Western societies. Initial molecular changes from exposure to HFr on liver metabolism may be essential to understand the perturbations leading to insulin resistance and abnormalities in lipid and carbohydrate metabolism. We studied vervet monkeys (Clorocebus aethiops sabaeus) fed a HFr (n=5) or chow diet (n=5) for 6 weeks, and obtained clinical measures of liver function, blood insulin, cholesterol and triglycerides. In addition, we performed untargeted global transcriptomics, proteomics, and metabolomics analyses on liver biopsies to determine the molecular impact of a HFr diet on coordinated pathways and networks that differed by diet.

RESULTS

We show that integration of omics data sets improved statistical significance for some pathways and networks, and decreased significance for others, suggesting that multiple omics datasets enhance confidence in relevant pathway and network identification. Specifically, we found that sirtuin signaling and a peroxisome proliferator activated receptor alpha (PPARA) regulatory network were significantly altered in hepatic response to HFr. Integration of metabolomics and miRNAs data further strengthened our findings.

CONCLUSIONS

Our integrated analysis of three types of omics data with pathway and regulatory network analysis demonstrates the usefulness of this approach for discovery of molecular networks central to a biological response. In addition, metabolites aspartic acid and docosahexaenoic acid (DHA), protein ATG3, and genes ATG7, and HMGCS2 link sirtuin signaling and the PPARA network suggesting molecular mechanisms for altered hepatic gluconeogenesis from consumption of a HFr diet.

Epigenetic clock and methylation studies in vervet monkeys

DNA methylation-based biomarkers of aging have been developed for many mammals but not yet for the vervet monkey (Chlorocebus sabaeus), which is a valuable non-human primate model for biomedical studies. We generated novel DNA methylation data from vervet cerebral cortex, blood, and liver using highly conserved mammalian CpGs represented on a custom array (HorvathMammalMethylChip40). We present six DNA methylation-based estimators of age: vervet multi-tissue epigenetic clock and tissue-specific clocks for brain cortex, blood, and liver. In addition, we developed two dual species clocks (human-vervet clocks) for measuring chronological age and relative age, respectively. Relative age was defined as ratio of chronological age to maximum lifespan to address the species differences in maximum lifespan. The high accuracy of the human-vervet clocks demonstrates that epigenetic aging processes are evolutionary conserved in primates. When applying these vervet clocks to tissue samples from another primate species, rhesus macaque, we observed high age correlations but strong offsets. We characterized CpGs that correlate significantly with age in the vervet. CpG probes that gain methylation with age across tissues were located near the targets of Polycomb proteins SUZ12 and EED and genes possessing the trimethylated H3K27 mark in their promoters. The epigenetic clocks are expected to be useful for anti-aging studies in vervets.

Resources for functional genomic studies of health and development in nonhuman primates

Primates display a wide range of phenotypic variation underlaid by complex genetically regulated mechanisms. The links among DNA sequence, gene function, and phenotype have been of interest from an evolutionary perspective, to understand functional genome evolution and its phenotypic consequences, and from a biomedical perspective to understand the shared and human-specific roots of health and disease. Progress in methods for characterizing genetic, transcriptomic, and DNA methylation (DNAm) variation is driving the rapid development of extensive omics resources, which are now increasingly available from humans as well as a growing number of nonhuman primates (NHPs). The fast growth of large-scale genomic data is driving the emergence of integrated tools and databases, thus facilitating studies of gene functionality across primates. This review describes NHP genomic resources that can aid in exploration of how genes shape primate phenotypes. It focuses on the gene expression trajectories across development in different tissues, the identification of functional genetic variation (including variants deleterious for protein function and regulatory variants modulating gene expression), and DNAm profiles as an emerging tool to understand the process of aging. These resources enable comparative functional genomics approaches to identify species-specific and primate-shared gene functionalities associated with health and development.

Neurodegenerative disease biomarkers Aβ1-40, Aβ1-42, tau, and p-tau181 in the vervet monkey cerebrospinal fluid: Relation to normal aging, genetic influences, and cerebral amyloid angiopathy

Background

The Caribbean vervet monkey (Chlorocebus aethiops sabaeus) is a potentially valuable animal model of neurodegenerative disease. However, the trajectory of aging in vervets and its relationship to human disease is incompletely understood.

Methods

To characterize biomarkers associated with neurodegeneration, we measured cerebrospinal fluid (CSF) concentrations of Aβ1-40, Aβ1-42, total tau, and p-tau181 in 329 members of a multigenerational pedigree. Linkage and genome-wide association were used to elucidate a genetic contribution to these traits.

Results

Aβ1-40 concentrations were significantly correlated with age, brain total surface area, and gray matter thickness. Levels of p-tau181 were associated with cerebral volume and brain total surface area. Among the measured analytes, only CSF Aβ1-40 was heritable. No significant linkage (LOD > 3.3) was found, though suggestive linkage was highlighted on chromosomes 4 and 12. Genome-wide association identified a suggestive locus near the chromosome 4 linkage peak.

Conclusions

Overall, these results support the vervet as a non-human primate model of amyloid-related neurodegeneration, such as Alzheimer's disease and cerebral amyloid angiopathy, and highlight Aβ1-40 and p-tau181 as potentially valuable biomarkers of these processes.

Obesity and obesogenic growth are both highly heritable and modified by diet in a nonhuman primate model, the African green monkey (Chlorocebus aethiops sabaeus)

Objective:

In humans, the ontogeny of obesity throughout the life course and the genetics underlying it has been historically difficult to study. We compared, in a non-human primate model, the lifelong growth trajectories of obese and non-obese adults to assess the heritability of and map potential genomic regions implicated in growth and obesity.

Study population:

A total of 905 African green monkeys, or vervets (Chlorocebus aethiops sabaeus) (472 females, 433 males) from a pedigreed captive colony.

Methods:

We measured fasted body weight (BW), crown-to-rump length (CRL), body-mass index (BMI) and waist circumference (WC) from 2000 to 2015. We used a longitudinal clustering algorithm to detect obesogenic growth, and logistic growth curves implemented in nonlinear mixed effects models to estimate three growth parameters. We used maximum likelihood variance decomposition methods to estimate the genetic contributions to obesity-related traits and growth parameters, including a test for the effects of a calorie-restricted dietary intervention. We used multipoint linkage analysis to map implicated genomic regions.

Results:

All measurements were significantly influenced by sex, and with the exception of WC, also influenced by maternal and post-natal diet. Chronic obesity outcomes were significantly associated with a pattern of extended growth duration with slow growth rates for BW. After accounting for environmental influences, all measurements were found to have a significant genetic component to variability. Linkage analysis revealed several regions suggested to be linked to obesity-related traits that are also implicated in human obesity and metabolic disorders.

Conclusions:

As in humans, growth patterns in vervets have a significant impact on adult obesity and are largely under genetic control with some evidence for maternal and dietary programming. These results largely mirror findings from human research, but reflect shorter developmental periods, suggesting that the vervet offers a strong genetic model for elucidating the ontogeny of human obesity.

Proteomics in non-human primates: utilizing RNA-Seq data to improve protein identification by mass spectrometry in vervet monkeys

Background

Shotgun proteomics utilizes a database search strategy to compare detected mass spectra to a library of theoretical spectra derived from reference genome information. As such, the robustness of proteomics results is contingent upon the completeness and accuracy of the gene annotation in the reference genome. For animal models of disease where genomic annotation is incomplete, such as non-human primates, proteogenomic methods can improve the detection of proteins by incorporating transcriptional data from RNA-Seq to improve proteomics search databases used for peptide spectral matching. Customized search databases derived from RNA-Seq data are capable of identifying unannotated genetic and splice variants while simultaneously reducing the number of comparisons to only those transcripts actively expressed in the tissue.

Results

We collected RNA-Seq and proteomic data from 10 vervet monkey liver samples and used the RNA-Seq data to curate sample-specific search databases which were analyzed in the program Morpheus. We compared these results against those from a search database generated from the reference vervet genome. A total of 284 previously unannotated splice junctions were predicted by the RNA-Seq data, 92 of which were confirmed by peptide spectral matches. More than half (53/92) of these unannotated splice variants had orthologs in other non-human primates, suggesting that failure to match these peptides in the reference analyses likely arose from incomplete gene model information. The sample-specific databases also identified 101 unique peptides containing single amino acid substitutions which were missed by the reference database. Because the sample-specific searches were restricted to actively expressed transcripts, the search databases were smaller, more computationally efficient, and identified more peptides at the empirically derived 1 % false discovery rate.

Conclusion

Proteogenomic approaches are ideally suited to facilitate the discovery and annotation of proteins in less widely studies animal models such as non-human primates. We expect that these approaches will help to improve existing genome annotations of non-human primate species such as vervet.

Electronic supplementary material

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

Regulators of mitochondrial quality control differ in subcutaneous fat of metabolically healthy and unhealthy obese monkeys

OBJECTIVE

Obesity exists with and without accompanying cardiometabolic disease, termed metabolically unhealthy obesity (MUO) and healthy obesity (MHO), respectively. Underlying differences in the ability of subcutaneous (SQ) fat to respond to nutrient excess are emerging as a key pathway. This study aimed to document the first spontaneous animal model of MHO and MUO and differences in SQ adipose tissue.

METHODS

Vervet monkeys (Chlorocebus aethiops; N = 171) were screened for metabolic syndrome. A subset of MHO and MUO monkeys (n = 6/group) had SQ fat biopsies collected for histological evaluations and examination of key mitochondrial proteins.

RESULTS

Obesity was seen in 20% of monkeys, and within this population, 31% were healthy, which mirrors human prevalence estimates. MUO monkeys had more than 60% lower adiponectin concentrations despite similar fat cell size, uncoupling protein 3, and activated macrophage abundance. However, alternatively activated/anti-inflammatory macrophages were 70% lower. Deficiencies of 50% or more in mitochondrial quality control regulators and selected mitochondrial fission and fusion markers were observed in the SQ fat of MUO monkeys despite comparable mitochondrial content.

CONCLUSIONS

A novel and translatable spontaneously obese animal model of MHO and MUO, occurring independently of dietary factors, was characterized. Differences in mitochondrial quality and inflammatory cell populations of subcutaneous fat may underpin divergent metabolic health.

Pair housing of Vervets/African Green Monkeys for biomedical research

Vervets, also known as African green monkeys, are a nonhuman primate species widely used in biomedical research. However, there are currently few references available describing techniques and rates of success for pair-housing this species. We present data from four cohorts of vervets from three different facilities: (i) the Wake Forest Vervet Research Colony (VRC; n = 72 female pairs, n= 52 male pairs), (ii) the University of Louisiana at Lafayette-New Iberia Research Center (UL-NIRC; n = 57 female pairs, n = 54 male pairs), (iii) the Tulane National Primate Research Center (TNRPC; n = 18 male pairs), and (iv) a cohort of imported males (n = 18 pairs) at Wake Forest. Compatibility was measured at 14, 30, and 60 days following introduction. Success rates for pair-housing at 14 days ranged from 96% to 98% for females and 96% to 100% for males at the VRC and UL-NIRC but were lower in the smaller imported male cohorts (TNPRC: 50%; WF: 28%). Among the UL-NIRC cohort and VRC male cohort, most of the pair separations after 14 days were due to reasons unrelated to social incompatibility. In contrast, a large proportion of TNPRC and imported male pairs successful at 14 days required separation within 60 days due to incompatibility. Multiple logistic regressions were performed using cohort, mean age of pair and weight difference between pair-mates as potential predictors of compatibility at 14 days. All three predicted the 14-day outcome in males but not females. A separate analysis in the VRC cohort found no evidence that prior familiarity in a group setting influenced outcomes. Variations in success rates across cohorts may have been influenced by introduction methodology. Behavioral differences between vervets and macaques, coupled with our findings, lead us to theorize that the gradual introduction techniques commonly implemented to pair house macaques may not be beneficial or suitable for this species. Am. J. Primatol. 79:e22501, 2017. © 2015 Wiley Periodicals, Inc.

Vitamin D heritability and effect of pregnancy status in Vervet monkeys (Chlorocebus aethiops sabaeus) under conditions of modest and high dietary supplementation

OBJECTIVES

The two objectives of the current study were to: 1) investigate the genetic contributions to variations in serum vitamin D concentrations under two dietary conditions (a standard monkey biscuit diet vs. a diet designed to model typical American consumption); and 2) explore the interaction of vitamin D with pregnancy status using a cohort of pedigreed female vervet/African green monkeys.

METHODS

This study includes 185 female (≥3.5 years) vervet/African green monkeys (Chlorocebus aethiops sabaeus) from a multi-generational, pedigreed breeding colony. The 25(OH)D3 concentrations were first measured seven to eight weeks after consuming a "typical American" diet (TAD), deriving 37, 18, and 45% of calories from fat, protein sources, and carbohydrates, and supplemented with vitamin D to a human equivalent of 1,000 IU/day. Vitamin D concentrations were assessed again when animals were switched to a low-fat, standard biscuit diet (LabDiet 5038) for 8 months, which provided a human equivalent of approximately 4,000 IU/day of vitamin D. All statistical analyses were implemented in SOLAR.

RESULTS

Pregnancy was associated with reduced 25(OH)D3 concentrations. Heritability analyses indicated a significant genetic contribution to 25(OH)D3 concentrations in the same monkeys consuming the biscuit diet (h(2) =0.66, P=0.0004) and TAD (h(2) =0.67, P=0.0078) diets, with higher 25(OH)D3 concentrations in animals consuming the biscuit diet. Additionally, there was a significant genotype-by-pregnancy status interaction on 25(OH)D3 concentrations (P<0.05) only among animals consuming the TAD diet.

DISCUSSION

These results support the existence of a genetic contribution to differences in serum 25(OH)D3 concentrations by pregnancy status and emphasize the role of diet (including vitamin D supplementation) in modifying genetic signals as well as vitamin D concentrations.

Variable responses of human and non-human primate gut microbiomes to a Western diet

BACKGROUND

The human gut microbiota interacts closely with human diet and physiology. To better understand the mechanisms behind this relationship, gut microbiome research relies on complementing human studies with manipulations of animal models, including non-human primates. However, due to unique aspects of human diet and physiology, it is likely that host-gut microbe interactions operate differently in humans and non-human primates.

RESULTS

Here, we show that the human microbiome reacts differently to a high-protein, high-fat Western diet than that of a model primate, the African green monkey, or vervet (Chlorocebus aethiops sabaeus). Specifically, humans exhibit increased relative abundance of Firmicutes and reduced relative abundance of Prevotella on a Western diet while vervets show the opposite pattern. Predictive metagenomics demonstrate an increased relative abundance of genes associated with carbohydrate metabolism in the microbiome of only humans consuming a Western diet.

CONCLUSIONS

These results suggest that the human gut microbiota has unique properties that are a result of changes in human diet and physiology across evolution or that may have contributed to the evolution of human physiology. Therefore, the role of animal models for understanding the relationship between the human gut microbiota and host metabolism must be re-focused.

Selection of models for the analysis of risk-factor trees: leveraging biological knowledge to mine large sets of risk factors with application to microbiome data

MOTIVATION

Establishment of a statistical association between microbiome features and clinical outcomes is of growing interest because of the potential for yielding insights into biological mechanisms and pathogenesis. Extracting microbiome features that are relevant for a disease is challenging and existing variable selection methods are limited due to large number of risk factor variables from microbiome sequence data and their complex biological structure.

RESULTS

We propose a tree-based scanning method, Selection of Models for the Analysis of Risk factor Trees (referred to as SMART-scan), for identifying taxonomic groups that are associated with a disease or trait. SMART-scan is a model selection technique that uses a predefined taxonomy to organize the large pool of possible predictors into optimized groups, and hierarchically searches and determines variable groups for association test. We investigate the statistical properties of SMART-scan through simulations, in comparison to a regular single-variable analysis and three commonly-used variable selection methods, stepwise regression, least absolute shrinkage and selection operator (LASSO) and classification and regression tree (CART). When there are taxonomic group effects in the data, SMART-scan can significantly increase power by using bacterial taxonomic information to split large numbers of variables into groups. Through an application to microbiome data from a vervet monkey diet experiment, we demonstrate that SMART-scan can identify important phenotype-associated taxonomic features missed by single-variable analysis, stepwise regression, LASSO and CART.

Systems biology of the vervet monkey

Nonhuman primates (NHP) provide crucial biomedical model systems intermediate between rodents and humans. The vervet monkey (also called the African green monkey) is a widely used NHP model that has unique value for genetic and genomic investigations of traits relevant to human diseases. This article describes the phylogeny and population history of the vervet monkey and summarizes the use of both captive and wild vervet monkeys in biomedical research. It also discusses the effort of an international collaboration to develop the vervet monkey as the most comprehensively phenotypically and genomically characterized NHP, a process that will enable the scientific community to employ this model for systems biology investigations.

Characterization of ovarian aging and reproductive senescence in vervet monkeys (Chlorocebus aethiops sabaeus)

Female vervet monkeys (Chlorocebus aethiops sabaeus) are used as an experimental model for chronic diseases relevant to women's health. However, reproductive senescence (menopause) has not yet been characterized for vervet monkeys. Here we describe the histologic, hormonal, and menstrual markers of reproductive senescence in vervet monkeys from the Wake Forest Vervet Research Colony. Ovaries from monkeys (age, 0 to 27 y) were serially sectioned (5 μm), stained, and photographed. In every 100th section, the numbers of primordial, primary, and secondary follicles were determined, and triplicate measurements were used to calculate mean numbers of follicles per ovary. Antimüllerian hormone (AMH), follicle stimulating hormone, and menstrual cycle length were measured in additional monkeys. Primordial follicles and AMH decreased significantly with age, and significant correlations between numbers of primordial and primary follicles and between numbers of primary and secondary follicles were noted. Histologic evaluation revealed that ovaries from 4 aged monkeys (older than 23 y) were senescent. One aged monkey transitioned to menopause, experiencing cycle irregularity over 4 y, eventual cessation of menses, and plasma AMH below the level of detection. Finally, with increasing age, the percentage of female vervets with offspring declined significantly. The present study provides insight into ovarian aging and reproductive senescence in vervet monkeys. Results highlight the importance of considering this nonhuman primate as a model to investigate the relationships between ovarian aging and chronic disease risk.