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.

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.

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.

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.