Search

Blood pressure (BP) is influenced by genetic variation and sodium intake with sex-specific differences; however, studies to identify renal molecular mechanisms underlying the influence of sodium intake on BP in nonhuman primates (NHP) have focused on males. To address the gap in our understanding of molecular mechanisms regulating BP in female primates, we studied sodium-naïve female baboons ( n = 7) fed a high-sodium (HS) diet for 6 wk. We hypothesized that in female baboons variation in renal transcriptional networks correlates with variation in BP response to a high-sodium diet. BP was continuously measured for 64-h periods throughout the study by implantable telemetry devices. Sodium intake, blood samples for clinical chemistries, and ultrasound-guided kidney biopsies were collected before and after the HS diet for RNA-Seq and bioinformatic analyses. We found that on the LS diet but not the HS diet, sodium intake and serum 17 β-estradiol concentration correlated with BP. Furthermore, kidney transcriptomes differed by diet-unbiased weighted gene coexpression network analysis revealed modules of genes correlated with BP on the HS diet but not the LS diet. Our results showed variation in BP on the HS diet correlated with variation in novel kidney gene networks regulated by ESR1 and MYC; i.e., these regulators have not been associated with BP regulation in male humans or rodents. Validation of the mechanisms underlying regulation of BP-associated gene networks in female NHP will inform better therapies toward greater precision medicine for women.

Fetal liver tissue collected from a nonhuman primate (NHP) baboon model of maternal nutrient reduction (MNR) at four gestational time points (90, 120, 140, and 165 days gestation [dG], term in the baboon is ∼185 dG) was used to quantify MNR effects on the fetal liver transcriptome. 28 transcripts demonstrated different expression patterns between MNR and control livers during the second half of gestation, a developmental period when the fetus undergoes rapid weight gain and fat accumulation. Differentially expressed transcripts were enriched for fatty acid oxidation and RNA splicing-related pathways. Increased RNA splicing activity in MNR was reflected in greater abundances of transcript splice variant isoforms in the MNR group. It can be hypothesized that the increase in splice variants is deployed in an effort to adapt to the poor in utero environment and ensure near-normal development and energy metabolism. This study is the first to study developmental programming across four critical gestational stages during primate fetal liver development and reveals a potentially novel cellular response mechanism mediating fetal programming in response to MNR.

MicroRNAs (miRNAs) regulate gene expression post-transcriptionally and circulate in the blood, making them attractive biomarkers of disease state for tissues like bone that are challenging to interrogate directly. Here, we report on five miRNAs-miR-197-3p, miR-320a, miR-320b, miR-331-5p, and miR-423-5p-associated with bone mineral density (BMD) in 147 healthy adult baboons. These baboons ranged in age from 15 to 25 years (45-75 human equivalent years) and 65% were female with a broad range of BMD values including a minority of osteopenic animals. miRNAs were generated via RNA sequencing from buffy coats collected at necropsy and areal BMD (aBMD) measured postmortem via dual-energy X-ray absorptiometry (DXA) of the lumbar vertebrae. Differential expression analysis controlled for the underlying pedigree structure of these animals to account for genetic variation which may drive miRNA abundance and aBMD values. While many of these miRNAs have been associated with the risk of osteoporosis in humans, this finding is of interest because the cohort represents a model of normal aging and bone metabolism rather than a disease cohort. The replication of miRNA associations with osteoporosis or other bone metabolic disorders in animals with healthy aBMD suggests an overlap in normal variation and disease states. We suggest that these miRNAs are involved in the regulation of cellular proliferation, apoptosis, and protein composition in the extracellular matrix throughout life; and age-related dysregulation of these systems may lead to disease. These miRNAs may be early indicators of progression to disease in advance of clinically detectible osteoporosis.

Blood pressure (BP) is influenced by genetic variation and sodium intake. Ninety percent of Americans consume more than the AHA recommended amount of sodium. Studies on the impact of genetic variation and sodium intake on BP in nonhuman primates (NHP) to date have focused on males. We hypothesized that variation in renal transcriptional networks correlate with BP response to high dietary sodium in female baboons. Sodium-naïve female baboons (n=7) were fed a low-sodium (LS) diet for 6 weeks followed by a high sodium (HS) diet for 6 weeks. Sodium intake, serum 17 betaestradiol, and ultrasound-guided kidney biopsies for RNA-Seq BP were collected at the end of each diet. BP was continuously measured for 64-hour periods throughout the study by implantable telemetry devices. On the LS diet, Na+ intake and serum 17 beta-estradiol concentration correlated with BP.Kidney transcriptomes differed by diet; analysis by unbiased weighted gene co-expression network analysis revealed modules of genes correlated with BP on the HS diet. Cell type composition of renal biopsies was consistent among all animals for both diets. Network analysis of module genes showed causal networks linking hormone receptors, proliferation and differentiation, methylation, hypoxia, insulin and lipid regulation, and inflammation as regulators underlying variation in BP on the HS diet. Our results show variation in BP correlated with novel kidney gene networks with master regulators PPARG and MYC in female baboons on a HS diet. Identification of mechanisms underlying regulators that influence BP will inform better therapies towards greater precision medicine for women.

BackgroundDietary 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.ResultsWe 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.ConclusionsOur 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 genesATG7, HMGCS2link sirtuin signaling and the PPARA network suggesting molecular mechanisms for altered hepatic gluconeogenesis from consumption of a HFr diet.

Hypertension is a complex disease influenced by sex, and genetic and environmental factors. Blood pressure (BP) is a continuous trait that is heritable in primates, including humans and baboons. The kidneys play a role in systemically regulating BP. Sex differences in BP onset and control with antihypertensive drug therapies have been observed in humans and rodents. Hypertension studies in nonhuman primates (NHP) to date have focused on males. We hypothesized that there are differences in renal molecular networks associated with BP in female and male primates. Sodium-naïve female (n=8) and male (n=9) baboons were fed a low-sodium chow diet prior to and during the study. Implantable telemetry devices continuously monitored heart rate and blood pressure over 24-hours, and ultrasound-guided kidney biopsies were collected for RNA-Seq. Serum 17 beta-estradiol concentration correlated BP in females. BP in males correlated with Na+ intake, blood urea nitrogen, and glucose. Cell type composition of renal biopsies was consistent between females and males. Sex differences were observed in the kidney transcriptomes by principal components analysis and weighted gene co-expression network analysis. Network analysis revealed HNF4A, ESR1, ESR2, SMARCA4, TP53, and NR3C1 as BP regulators in males. Our results demonstrate sex differences in primate kidney molecular networks and provide evidence of a novel link between renal transcription factors and BP regulation in males. Understanding sex differences and transcriptome variation in primate kidneys correlated with BP and clinical measures associated with BP will inform better therapies towards the goal of precision medicine for women and men.

Key points Maternal obesity (MO) and exposure to a high-fat, high-simple-carbohydrate diet during pregnancy predisposes offspring to obesity, metabolic and cardiovascular disorders in later life. Underlying molecular pathways and potential epigenetic factors that are dysregulated in MO were identified using unbiased transcriptomic methods. There was increased lipid accumulation and severe steatosis in the MO baboon fetal liver suggesting that these offspring are on an early trajectory of non-alcoholic fatty liver disease and non-alcoholic steatohepatitis. Abstract Maternal obesity (MO) increases offspring cardiometabolic disease risk. Altered fetal liver development in response to the challenge of MO has metabolic consequences underlying adverse offspring life-course health outcomes. Little is known about the molecular pathways and potential epigenetic changes regulating primate fetal liver responses to MO. We hypothesized that MO would induce fetal baboon liver epigenetic changes resulting in dysregulation of key metabolic pathways that impact lipid metabolism. MO was induced prior to pregnancy by a high-fat, high-fructose diet. Unbiased gene and microRNA (small RNA Seq) abundance analyses were performed on fetal baboon livers at 0.9 gestation and subjected to pathway analyses to identify fetal liver molecular responses to MO. Fetal baboon liver lipid and glycogen content were quantified by the Computer Assisted Stereology Toolbox. In response to MO, fetal livers revealed dysregulation of TCA cycle, proteasome, oxidative phosphorylation, glycolysis and Wnt/β-catenin signalling pathways together with marked lipid accumulation supporting our hypothesis that multiple pathway dysregulation detrimentally impacts lipid management. This is the first study of MO programming of the non-human primate fetal liver using unbiased transcriptome analysis to detect changes in hepatic gene expression levels and identify potential microRNA epigenetic regulators of metabolic disruption.

BACKGROUND: Little is known about the repertoire of nonhuman primate kidney genes expressed throughout development. The present work establishes an understanding of the primate renal transcriptome during fetal development in the context of renal maturation. METHODS: The baboon kidney transcriptome was characterized at 60 days gestation (DG), 90DG, 125DG, 140DG, 160DG and adulthood (6–12 years) using gene arrays and validated by QRT-PCR. Pathway and cluster analyses were used to characterize gene expression in the context of biological pathways. RESULTS: Pathway analysis indicated activation of pathways not previously reported as relevant to kidney development. Cluster analysis also revealed gene splice variants with discordant expression profiles during development. CONCLUSIONS: This study provides the first detailed genetic analysis of the developing primate kidney, and our findings of discordant expression of gene splice variants suggest that gene arrays likely provide a simplified view and demonstrate the need to study the fetal renal proteome.

Schizophrenia is a debilitating disorder affecting just under 1% of the population. While the symptoms of this disorder do not appear until late adolescence, pathological alterations likely occur earlier, during development in utero. While there is an increasing literature examining transcriptome alterations in patients, it is not possible to examine the changes in gene expression that occur during development in humans that will develop schizophrenia. Here we utilize three distinct rodent developmental disruption models of schizophrenia to examine potential overlapping alterations in the transcriptome, with a specific focus on markers of interneuron development. Specifically, we administered either methylazoxymethanol acetate (MAM), Polyinosinic:polycytidylic acid (Poly I:C), or chronic protein malnutrition, on GD 17 and examined mRNA expression in the developing hippocampus of the offspring 18 hours later. Here, we report alterations in gene expression that may contribute to the pathophysiology of schizophrenia, including significant alterations in interneuron development and ribosome function.

RationalePlasma low-density lipoprotein cholesterol (plasma LDL-C), vascular endothelial cells and peripheral blood mononuclear cells (PBMCs), particularly monocytes, play key roles in initiating atherosclerosis, the primary cause of cardiovascular disease (CVD). Although the mechanisms underlying development of atherosclerosis are not well understood, LDL-C is known to influence expression of endothelial microRNAs (miRNAs) and gene-targets of miRNAs to promote cell senescence. However, the impact of LDL-C on expression of PBMC miRNAs and miRNA targeted genes in response to an atherogenic diet is not known. In this study, we used unbiased methods to identify coordinately responsive PBMC miRNA- gene networks that differ between low and high LDL-C baboons when fed a high-cholesterol, high-fat (HCHF) diet.Methods and resultsUsing RNA Seq, we quantified PBMC mRNAs and miRNAs from half-sib baboons discordant for LDL-C plasma concentrations (low LDL-C, n = 3; high LDL-C, n = 3) before and after a 7-week HCHF diet challenge. For low LDL-C baboons, 626 genes exhibited significant change in expression (255 down-regulated, 371 up-regulated) in response to the HCHF diet, and for high LDL-C baboons 379 genes exhibited significant change in expression (162 down-regulated, 217 up-regulated) in response to the HCHF diet. We identified 494 miRNAs identical to human miRNAs and 47 novel miRNAs. Fifty miRNAs were differentially expressed in low LDL-C baboons (21 up- and 29 down-regulated) and 20 in high LDL-C baboons (11 up- and 9 down-regulated) in response to the HCHF diet. Among the differentially expressed miRNAs were miR-221/222 and miR-34a-3p, which were down-regulated, and miR-148a/b-5p, which was up-regulated. In addition, gene-targets of these miRNAs, VEGFA, MAML3, SPARC, and DMGDH, were inversely expressed and are central hub genes in networks and signaling pathways that differ between low and high LDL-C baboon HCHF diet response.ConclusionsWe have identified coordinately regulated HCHF diet-responsive PBMC miRNA-gene networks that differ between baboons discordant for LDL-C concentrations. Our findings provide potential insights into molecular mechanisms underlying initiation of atherosclerosis where LDL-C concentrations influence expression of specific miRNAs, which in turn regulate expression of genes that play roles in initiation of lesions.

Baboons (genus Papio) are broadly studied in the wild and in captivity. They are widely used as a nonhuman primate model for biomedical studies, and the Southwest National Primate Research Center (SNPRC) at Texas Biomedical Research Institute has maintained a large captive baboon colony for more than 50 yr. Unlike other model organisms, however, the genomic resources for baboons are severely lacking. This has hindered the progress of studies using baboons as a model for basic biology or human disease. Here, we describe a data set of 100 high-coverage whole-genome sequences obtained from the mixed colony of olive (P. anubis) and yellow (P. cynocephalus) baboons housed at the SNPRC. These data provide a comprehensive catalog of common genetic variation in baboons, as well as a fine-scale genetic map. We show how the data can be used to learn about ancestry and admixture and to correct errors in the colony records. Finally, we investigated the consequences of inbreeding within the SNPRC colony and found clear evidence for increased rates of infant mortality and increased homozygosity of putatively deleterious alleles in inbred individuals.

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.

Maternal undernutrition increases the risk of perinatal complications and predisposes offspring to obesity, diabetes, and cardiovascular disease later in life. Emerging evidence suggests that changes in placental function play a role in linking altered maternal nutrition in pregnancy to the subsequent development of adult disease. The susceptibility for disease in response to an adverse intrauterine environment differs distinctly between boys and girls, with girls typically having better outcomes. Here, we tested the hypothesis that regulation of the placental transcriptome by maternal nutrient reduction (NR) is dependent on fetal sex. We used a nonhuman primate model of NR in which maternal global food intake was reduced by 30% in baboons starting at gestational day (GD) 30. At GD 165 (term = GD 183), placental genome expression profiling of 6 control (n = 3 females, 3 males) and 6 nutrient restricted (n = 3 females, 3 males) fetuses was carried out followed by bioinformatic analysis. Surprisingly, there was no coordinated placental molecular response to decreased nutrient availability when analyzing the data without accounting for fetal sex. In contrast, female placentas exhibited a highly coordinated response that included upregulation of genes in networks, pathways, and functional groups related to programmed cell death and downregulation of genes in networks, pathways, and functional groups associated with cell proliferation. These changes were not apparent in the male placentas. Our data support the concept that female placentas initiate complex adaptive responses to an adverse intrauterine environment, which may contribute to increased survival and better pregnancy outcomes in girls.

Currently available methods for interrogating DNA-protein interactions at individual genomic loci have significant limitations, and make it difficult to work with unmodified cells or examine single-copy regions without specific antibodies. In this study, we describe a physiological application of the Hybridization Capture of Chromatin-Associated Proteins for Proteomics (HyCCAPP) methodology we have developed. Both novel and known locus-specific DNA-protein interactions were identified at the ENO2 and GAL1 promoter regions of Saccharomyces cerevisiae, and revealed subgroups of proteins present in significantly different levels at the loci in cells grown on glucose versus galactose as the carbon source. Results were validated using chromatin immunoprecipitation. Overall, our analysis demonstrates that HyCCAPP is an effective and flexible technology that does not require specific antibodies nor prior knowledge of locally occurring DNA-protein interactions and can now be used to identify changes in protein interactions at target regions in the genome in response to physiological challenges.

The pregnant sheep has provided seminal insights into reproduction related to animal and human development (ovarian function, fertility, implantation, fetal growth, parturition and lactation). Fetal sheep physiology has been extensively studied since 1950, contributing significantly to the basis for our understanding of many aspects of fetal development and behaviour that remain in use in clinical practice today. Understanding mechanisms requires the combination of systems approaches uniquely available in fetal sheep with the power of genomic studies. Absence of the full range of sheep genomic resources has limited the full realization of the power of this model, impeding progress in emerging areas of pregnancy biology such as developmental programming. We have examined the expressed fetal sheep heart transcriptome using high-throughput sequencing technologies. In so doing we identified 36,737 novel transcripts and describe genes, gene variants and pathways relevant to fundamental developmental mechanisms. Genes with the highest expression levels and with novel exons in the fetal heart transcriptome are known to play central roles in muscle development. We show that high-throughput sequencing methods can generate extensive transcriptome information in the absence of an assembled and annotated genome for that species. The gene sequence data obtained provide a unique genomic resource for sheep specific genetic technology development and, combined with the polymorphism data, augment annotation and assembly of the sheep genome. In addition, identification and pathway analysis of novel fetal sheep heart transcriptome splice variants is a first step towards revealing mechanisms of genetic variation and gene environment interactions during fetal heart development.

Introduction: Cardiovascular disease is the leading cause of death in developed countries and commonly due to atherosclerosis. Atherosclerosis is a resultant of interplay of gene and environmental factors including diet. Liver is the primary target for lipid metabolism and circulatory microRNAs (miRNAs) are potential biomarkers for atherosclerosis. However, little information is available on miRNAs for early prediction and diagnosis of severity of atherosclerosis. Hypothesis: We hypothesized that peripheral blood mononuclear cell (PBMCs) miRNA(s) are potential biomarkers for atherosclerosis. We tested this by first identifying PBMCs and liver miRNAs with discordant expression in baboons with low and high serum low density-lipoprotein cholesterol (LDL-C) concentrations and then determined if the miRNAs correlated with atherosclerotic lesions in baboons fed a high fat diet for 2 years. Methods and Results: Using Next Gen methods, we sequenced PBMCs and liver small RNA libraries from baboons discordant for LDL-C concentrations (low LDL-C, n = 3; high LDL-C, n = 3); fed a high-cholesterol, high-fat (HCHF) diet for 7 weeks. We identified 47 and 27 novel baboon miRNAs, and 494 and 490 human homologues in PBMCs and in liver, respectively. Seventy-one and 30 miRNAs showed differential response to HCHF diet in PBMCs and in liver, respectively. Comparison of miRNA expression profiles from PBMCs and liver showed concordant expression for miR-221/222 family; it was down-regulated in high LDL-C baboons in response to HCHF diet. Using qRT-PCR, we quantified miR-221 abundance in PBMCs and in livers of high LDL-C baboons (n = 109) on 2-year HCHF challenge and revealed, for hepatic miR-221, a positive correlation with severity of lesions in descending aorta (Pearson correlation coefficient, r = 0.18, p = 0.05) and aortic arch (r = 0.22, p = 0.03). Conclusions: Liver miR-221 abundance is associated with extent of lesions in baboons fed a HCHF diet for 2 years and may be a biomarker for atherosclerosis.

The pregnant sheep has provided seminal insights into reproduction related to animal and human development (ovarian function, fertility, implantation, fetal growth, parturition and lactation). Fetal sheep physiology has been extensively studied since 1950, contributing significantly to the basis for our understanding of many aspects of fetal development and behaviour that remain in use in clinical practice today. Understanding mechanisms requires the combination of systems approaches uniquely available in fetal sheep with the power of genomic studies. Absence of the full range of sheep genomic resources has limited the full realization of the power of this model, impeding progress in emerging areas of pregnancy biology such as developmental programming. We have examined the expressed fetal sheep heart transcriptome using high-throughput sequencing technologies. In so doing we identified 36,737 novel transcripts and describe genes, gene variants and pathways relevant to fundamental developmental mechanisms. Genes with the highest expression levels and with novel exons in the fetal heart transcriptome are known to play central roles in muscle development. We show that high-throughput sequencing methods can generate extensive transcriptome information in the absence of an assembled and annotated genome for that species. The gene sequence data obtained provide a unique genomic resource for sheep specific genetic technology development and, combined with the polymorphism data, augment annotation and assembly of the sheep genome. In addition, identification and pathway analysis of novel fetal sheep heart transcriptome splice variants is a first step towards revealing mechanisms of genetic variation and gene environment interactions during fetal heart development.

1. Vascular inflammation plays a critical role in atherogenesis. Previously, we showed that baboon arterial endothelial cells (BAEC) were hyporesponsive to lipopolysaccharide (LPS) compared with human arterial endothelial cells (HAEC). 2. In the present study, we investigated mechanisms underlying differential responses between HAEC and BAEC to tumour necrosis factor (TNF)-alpha and LPS. 3. Both HAEC and BAEC responded similarly to TNF-alpha. However, BAEC showed retarded responses to LPS in expression of E-selectin, intercellular adhesion molecule-1, monocyte chemotactic protein-1 and interleukin-8 (P < 0.05). These changes were confirmed at the mRNA level. Tumour necrosis factor-alpha activated nuclear factor-kappaB members such as p50, p52, p65, c-rel and RelB in both HAEC and BAEC. In contrast, LPS activated p50 and p65 only in HAEC. Using microarray assays, we found that TNF receptor-associated factor 2 (TRAF-2), TNF receptor superfamily, member 1A-associated via death domain (TRADD) and nuclear factors such as nuclear factor of kappa in B-cells inhibitor, alpha (NFKBIA) and nuclear factor of kappa in B-cells inhibitor, beta (NFKBIB) were upregulated by LPS only in HAEC. Although the baseline expression of Toll-like receptor (TLR) 4 was low in both HAEC and BAEC, TNF-alpha activated TLR4 expression in both cell types. Although LPS increased TLR4 expression only in HAEC, human and baboon peripheral blood mononuclear cells exhibited similar TLR4 expression and response to LPS. Transfecting BAEC with TLR4/myeloid differentiation protein-2 overexpression vector conferred BAEC responsiveness to LPS. 4. The findings of the present study indicate that an altered TLR4 system may be responsible for the resistance of baboon endothelial cells to LPS. Given the importance of TLR4 in human immune responses and vascular diseases, the natural resistance of baboons to LPS/TLR4-initiated inflammation could make the baboon a valuable animal model in which to study how inflammation affects atherogenesis.

We have developed an integrated approach, using genetic and genomic methods, in conjunction with resources from the Southwest National Primate Research Center (SNPRC) baboon colony, for the identification of genes and their functional variants that encode quantitative trait loci (QTL). In addition, we use comparative genomic methods to overcome the paucity of baboon specific reagents and to augment translation of our findings in a nonhuman primate (NHP) to the human population. We are using the baboon as a model to study the genetics of cardiovascular disease (CVD). A key step for understanding gene-environment interactions in cardiovascular disease is the identification of genes and gene variants that influence CVD phenotypes. We have developed a sequential methodology that takes advantage of the SNPRC pedigreed baboon colony, the annotated human genome, and current genomic and bioinformatic tools. The process of functional polymorphism identification for genes encoding QTLs involves comparison of expression profiles for genes and predicted genes in the genomic region of the QTL for individuals discordant for the phenotypic trait mapping to the QTL. After comparison, genes of interest are prioritized, and functional polymorphisms are identified in candidate genes by genotyping and quantitative trait nucleotide analysis. This approach reduces the time and labor necessary to prioritize and identify genes and their polymorphisms influencing variation in a quantitative trait compared with traditional positional cloning methods.

Carboxylesterase (CES) is predominantly responsible for the detoxification of a wide range of drugs and narcotics, and catalyze several reactions in cholesterol and fatty acid metabolism. Studies of the genetic and biochemical properties of primate CES may contribute to an improved understanding of human disease, including atherosclerosis, obesity and drug addiction, for which non-human primates serve as useful animal models. We cloned and sequenced baboon CES1 and CES2 and used in vitro and in silico methods to predict protein secondary and tertiary structures, and examined evolutionary relationships for these enzymes with other primate and mouse CES orthologs. We found that baboon CES1 and CES2 proteins retained extensive similarity with human CES1 and CES2, shared key structural features reported for human CES1, and showed family specific sequences consistent with their multimeric and monomeric subunit structures respectively.

The baboon is an invaluable model for the study of human health and disease, including many complex diseases of the kidney. Although scientists have made great progress in developing this animal as a model for numerous areas of biomedical research, genomic resources for the baboon, such as a quality annotated genome, are still lacking. To this end, we characterized the baboon kidney transcriptome using high-throughput cDNA sequencing (RNA-Seq) to identify genes, gene variants, single nucleotide polymorphisms (SNPs), insertion-deletion polymorphisms (InDels), cellular functions, and key pathways in the baboon kidney to provide a genomic resource for the baboon. Analysis of our sequencing data revealed 45,499 high-confidence SNPs and 29,813 InDels comparing baboon cDNA sequences with the human hg18 reference assembly and identified 35,900 cDNAs in the baboon kidney, including 35,150 transcripts representing 15,369 genic genes that are novel for the baboon. Gene ontology analysis of our sequencing dataset also identified numerous biological functions and canonical pathways that were significant in the baboon kidney, including a large number of metabolic pathways that support known functions of the kidney. The results presented in this study catalogues the transcribed mRNAs, noncoding RNAs, and hypothetical proteins in the baboon kidney and establishes a genomic resource for scientists using the baboon as an experimental model.

Background Dysregulation of microRNA (miRNA) expression has been implicated in molecular genetic events leading to the progression and development of atherosclerosis. We hypothesized that miRNA expression profiles differ between baboons with low and high serum low-density lipoprotein cholesterol (LDL-C) concentrations in response to diet, and that a subset of these miRNAs regulate genes relevant to dyslipidemia and risk of atherosclerosis. Results Using Next Generation Illumina sequencing methods, we sequenced hepatic small RNA libraries from baboons differing in their LDL-C response to a high-cholesterol, high-fat (HCHF) challenge diet (low LDL-C, n = 3; high LDL-C, n = 3), resulting in 517 baboon miRNAs: 490 were identical to human miRNAs and 27 were novel. We compared miRNA expression profiles from liver biopsies collected before and after the challenge diet and observed that HCHF diet elicited expression of more miRNAs compared to baseline (chow) diet for both low and high LDL-C baboons. Eighteen miRNAs exhibited differential expression in response to HCHF diet in high LDL-C baboons compared to 10 miRNAs in low LDL-C baboons. We used TargetScan/Base tools to predict putative miRNA targets; miRNAs expressed in high LDL-C baboons had significantly more gene targets than miRNAs expressed in low LDL-C responders. Further, we identified miRNA isomers and other non-coding RNAs that were differentially expressed in response to the challenge diet in both high LDL-C and low LDL-C baboons. Conclusions We sequenced and annotated baboon liver miRNAs from low LDL-C and high LDL-C responders using high coverage Next Gen sequencing methods, determined expression changes in response to a HCHF diet challenge, and predicted target genes regulated by the differentially expressed miRNAs. The identified miRNAs will enrich the database for non-coding small RNAs including the extent of variation in these sequences. Further, we identified other small non-coding RNAs differentially expressed in response to diet. Our discovery of differentially expressed baboon miRNAs in response to a HCHF diet challenge that differ by LDL-C phenotype is a fundamental step in understating the role of non-coding RNAs in dyslipidemia.

Background MicroRNAs (miRNAs) are small noncoding RNAs (~22 nucleotides) that regulate gene expression by cleaving mRNAs or inhibiting translation. The baboon is a well-characterized cardiovascular disease model; however, no baboon miRNAs have been identified. Evidence indicates that the baboon and human genomes are highly conserved; based on this conservation, we hypothesized that comparative genomic methods could be used to identify baboon miRNAs. Methods We employed an in silico comparative genomics approach and human miRNA arrays to identify baboon expressed miRNAs in liver (n = 6) and lymphocytes (n = 6). Expression profiles for selected miRNAs in multiple tissues were validated by RT-PCR. Results We identified in silico 555 putative baboon pre-miRNAs, of which 41% exhibited 100% identity and an additional 58% shared more than 90% sequence identity with human pre-miRNAs. Some of these miRNAs are primate-specific and are clustered in the baboon genome like human miRNA clusters. We detected expression of 494 miRNAs on the microarray and validated expression of selected miRNAs in baboon liver and lymphocytes by RT-PCR. We also observed miRNA expression in additional tissues relevant to dyslipidemia and atherosclerosis. Approximately half of the miRNAs expressed on the array were not predicted in silico suggesting that we have identified novel baboon miRNAs, which could not be predicted using the current draft of the baboon genome. Conclusion We identified a subset of baboon miRNAs using a comparative genomic approach, identified additional baboon miRNAs using a human array and showed tissue-specific expression of baboon miRNAs. Our discovery of baboon miRNAs in liver and lymphocytes will provide resources for studies on the roles of miRNAs in dyslipidemia and atherosclerosis, and for translational studies.