A genetic linkage map of the baboon (Papio hamadryas) genome based on human microsatellite polymorphisms

Global and local ancestry estimation in a captive baboon colony

The last couple of decades have highlighted the importance of studying hybridization, particularly among primate species, as it allows us to better understand our own evolutionary trajectory. Here, we report on genetic ancestry estimates using dense, full genome data from 881 olive (Papio anubus), yellow (Papio cynocephalus), or olive-yellow crossed captive baboons from the Southwest National Primate Research Center. We calculated global and local ancestry information, imputed low coverage genomes (n = 830) to improve marker quality, and updated the genetic resources of baboons available to assist future studies. We found evidence of historical admixture in some putatively purebred animals and identified errors within the Southwest National Primate Research Center pedigree. We also compared the outputs between two different phasing and imputation pipelines along with two different global ancestry estimation software. There was good agreement between the global ancestry estimation software, with R2 > 0.88, while evidence of phase switch errors increased depending on what phasing and imputation pipeline was used. We also generated updated genetic maps and created a concise set of ancestry informative markers (n = 1,747) to accurately obtain global ancestry estimates.

The temporal and genomic scale of selection following hybridization

Genomic evidence supports an important role for selection in shaping patterns of introgression along the genome, but frameworks for understanding the evolutionary dynamics within hybrid populations that underlie these patterns have been lacking. Due to the clock-like effect of recombination in hybrids breaking up parental haplotypes, drift and selection produce predictable patterns of ancestry variation at varying spatial genomic scales through time. Here, we develop methods based on the Discrete Wavelet Transform to study the genomic scale of local ancestry variation and its association with recombination rates and show that these methods capture temporal dynamics of drift and genome-wide selection after hybridization. We apply these methods to published datasets from hybrid populations of swordtail fish (Xiphophorus) and baboons (Papio) and to inferred Neanderthal introgression in modern humans. Across systems, upward of 20% of variation in local ancestry at the broadest genomic scales can be attributed to systematic selection against introgressed alleles, consistent with strong selection acting on early-generation hybrids. Signatures of selection at fine genomic scales suggest selection over longer time scales; however, we suggest that our ability to confidently infer selection at fine scales is likely limited by inherent biases in current methods for estimating local ancestry from contiguous segments of genomic similarity. Wavelet approaches will become widely applicable as genomic data from systems with introgression become increasingly available and can help shed light on generalities of the genomic consequences of interspecific hybridization.

Novel Insights into the Landscape of Crossover and Noncrossover Events in Rhesus Macaques (Macaca mulatta)

Meiotic recombination landscapes differ greatly between distantly and closely related taxa, populations, individuals, sexes, and even within genomes; however, the factors driving this variation are yet to be well elucidated. Here, we directly estimate contemporary crossover rates and, for the first time, noncrossover rates in rhesus macaques (Macaca mulatta) from four three-generation pedigrees comprising 32 individuals. We further compare these results with historical, demography-aware, linkage disequilibrium-based recombination rate estimates. From paternal meioses in the pedigrees, 165 crossover events with a median resolution of 22.3 kb were observed, corresponding to a male autosomal map length of 2,357 cM-approximately 15% longer than an existing linkage map based on human microsatellite loci. In addition, 85 noncrossover events with a mean tract length of 155 bp were identified-similar to the tract lengths observed in the only other two primates in which noncrossovers have been studied to date, humans and baboons. Consistent with observations in other placental mammals with PRDM9-directed recombination, crossover (and to a lesser extent noncrossover) events in rhesus macaques clustered in intergenic regions and toward the chromosomal ends in males-a pattern in broad agreement with the historical, sex-averaged recombination rate estimates-and evidence of GC-biased gene conversion was observed at noncrossover sites.

The temporal and genomic scale of selection following hybridization

Genomic evidence supports an important role for selection in shaping patterns of introgression along the genome, but frameworks for understanding the dynamics underlying these patterns within hybrid populations have been lacking. Here, we develop methods based on the Wavelet Transform to understand the spatial genomic scale of local ancestry variation and its association with recombination rates. We present theory and use simulations to show how wavelet-based decompositions of ancestry variance along the genome and the correlation between ancestry and recombination reflect the joint effects of recombination, genetic drift, and genome-wide selection against introgressed alleles. Due to the clock-like effect of recombination in hybrids breaking up parental haplotypes, drift and selection produce predictable patterns of local ancestry variation at varying spatial genomic scales through time. Using wavelet approaches to identify the genomic scale of variance in ancestry and its correlates, we show that these methods can detect temporally localized effects of drift and selection. We apply these methods to previously published datasets from hybrid populations of swordtail fish (Xiphophorus) and baboons (Papio), and to inferred Neanderthal introgression in modern humans. Across systems, we find that upwards of 20% of the variation in local ancestry at the broadest genomic scales can be attributed to systematic selection against introgressed alleles, consistent with strong selection acting on early-generation hybrids. We also see signals of selection at fine genomic scales and much longer time scales. However, we show that our ability to confidently infer selection at fine scales is likely limited by inherent biases in current methods for estimating local ancestry from genomic similarity. Wavelet approaches will become widely applicable as genomic data from systems with introgression become increasingly available, and can help shed light on generalities of the genomic consequences of interspecific hybridization.

Age- and sex-related differences in baboon (Papio anubis) gray matter covariation

Age-related changes in cognition, brain morphology, and behavior are exhibited in several primate species. Baboons, like humans, naturally develop Alzheimer's disease-like pathology and cognitive declines with age and are an underutilized model for studies of aging. To determine age-related differences in gray matter covariation of 89 olive baboons (Papio anubis), we used source-based morphometry (SBM) to analyze data from magnetic resonance images. We hypothesized that we would find significant age effects in one or more SBM components, particularly those which include regions influenced by age in humans and other nonhuman primates (NHPs). A multivariate analysis of variance revealed that individual weighted gray matter covariation scores differed across the age classes. Elderly baboons contributed significantly less to gray matter covariation components including the brainstem, superior parietal cortex, thalamus, and pallidum compared to juveniles, and middle and superior frontal cortex compared to juveniles and young adults (p < 0.05). Future studies should examine the relationship between the changes in gray matter covariation reported here and age-related cognitive decline.

Blood pressure and the kidney cortex transcriptome response to high-sodium diet challenge in female nonhuman primates

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.

High-Resolution Estimates of Crossover and Noncrossover Recombination from a Captive Baboon Colony

Homologous recombination has been extensively studied in humans and a handful of model organisms. Much less is known about recombination in other species, including nonhuman primates. Here, we present a study of crossovers (COs) and noncrossover (NCO) recombination in olive baboons (Papio anubis) from two pedigrees containing a total of 20 paternal and 17 maternal meioses, and compare these results to linkage disequilibrium (LD) based recombination estimates from 36 unrelated olive baboons. We demonstrate how COs, combined with LD-based recombination estimates, can be used to identify genome assembly errors. We also quantify sex-specific differences in recombination rates, including elevated male CO and reduced female CO rates near telomeres. Finally, we add to the increasing body of evidence suggesting that while most NCO recombination tracts in mammals are short (e.g., <500 bp), there is a non-negligible fraction of longer (e.g., >1 kb) NCO tracts. For NCO tracts shorter than 10 kb, we fit a mixture of two (truncated) geometric distributions model to the NCO tract length distribution and estimate that >99% of all NCO tracts are very short (mean 24 bp), but the remaining tracts can be quite long (mean 4.3 kb). A single geometric distribution model for NCO tract lengths is incompatible with the data, suggesting that LD-based methods for estimating NCO recombination rates that make this assumption may need to be modified.

A Baboon Brain Atlas for Magnetic Resonance Imaging and Positron Emission Tomography Image Analysis

The olive baboon (Papio anubis) is phylogenetically proximal to humans. Investigation into the baboon brain has shed light on the function and organization of the human brain, as well as on the mechanistic insights of neurological disorders such as Alzheimer's and Parkinson's. Non-invasive brain imaging, including positron emission tomography (PET) and magnetic resonance imaging (MRI), are the primary outcome measures frequently used in baboon studies. PET functional imaging has long been used to study cerebral metabolic processes, though it lacks clear and reliable anatomical information. In contrast, MRI provides a clear definition of soft tissue with high resolution and contrast to distinguish brain pathology and anatomy, but lacks specific markers of neuroreceptors and/or neurometabolites. There is a need to create a brain atlas that combines the anatomical and functional/neurochemical data independently available from MRI and PET. For this purpose, a three-dimensional atlas of the olive baboon brain was developed to enable multimodal imaging analysis. The atlas was created on a population-representative template encompassing 89 baboon brains. The atlas defines 24 brain regions, including the thalamus, cerebral cortex, putamen, corpus callosum, and insula. The atlas was evaluated with four MRI images and 20 PET images employing the radiotracers for [11C]benzamide, [11C]metergoline, [18F]FAHA, and [11C]rolipram, with and without structural aids like [18F]flurodeoxyglycose images. The atlas-based analysis pipeline includes automated segmentation, registration, quantification of region volume, the volume of distribution, and standardized uptake value. Results showed that, in comparison to PET analysis utilizing the "gold standard" manual quantification by neuroscientists, the performance of the atlas-based analysis was at >80 and >70% agreement for MRI and PET, respectively. The atlas can serve as a foundation for further refinement, and incorporation into a high-throughput workflow of baboon PET and MRI data. The new atlas is freely available on the Figshare online repository (https://doi.org/10.6084/m9.figshare.16663339), and the template images are available from neuroImaging tools & resources collaboratory (NITRC) (https://www.nitrc.org/projects/haiko89/).

Research Relevant Conditions and Pathology in Nonhuman Primates

Biomedical research involving animal models continues to provide important insights into disease pathogenesis and treatment of diseases that impact human health. In particular, nonhuman primates (NHPs) have been used extensively in translational research due to their phylogenetic proximity to humans and similarities to disease pathogenesis and treatment responses as assessed in clinical trials. Microscopic changes in tissues remain a significant endpoint in studies involving these models. Spontaneous, expected (ie, incidental or background) histopathologic changes are commonly encountered and influenced by species, genetic variations, age, and geographical origin of animals, including exposure to infectious or parasitic agents. Often, the background findings confound study-related changes, because numbers of NHPs used in research are limited by animal welfare and other considerations. Moreover, background findings in NHPs can be exacerbated by experimental conditions such as treatment with xenobiotics (eg, infectious morphological changes related to immunosuppressive therapy). This review and summary of research-relevant conditions and pathology in rhesus and cynomolgus macaques, baboons, African green monkeys, common marmosets, tamarins, and squirrel and owl monkeys aims to improve the interpretation and validity of NHP studies.

A Fine-Scale Genetic Map for Vervet Monkeys

Despite its important biological role, the evolution of recombination rates remains relatively poorly characterized. This owes, in part, to the lack of high-quality genomic resources to address this question across diverse species. Humans and our closest evolutionary relatives, anthropoid apes, have remained a major focus of large-scale sequencing efforts, and thus recombination rate variation has been comparatively well studied in this group-with earlier work revealing a conservation at the broad- but not the fine-scale. However, in order to better understand the nature of this variation, and the time scales on which substantial modifications occur, it is necessary to take a broader phylogenetic perspective. I here present the first fine-scale genetic map for vervet monkeys based on whole-genome population genetic data from ten individuals and perform a series of comparative analyses with the great apes. The results reveal a number of striking features. First, owing to strong positive correlations with diversity and weak negative correlations with divergence, analyses suggest a dominant role for purifying and background selection in shaping patterns of variation in this species. Second, results support a generally reduced broad-scale recombination rate compared with the great apes, as well as a narrower fraction of the genome in which the majority of recombination events are observed to occur. Taken together, this data set highlights the great necessity of future research to identify genomic features and quantify evolutionary processes that are driving these rate changes across primates.

Dark Matter of Primate Genomes: Satellite DNA Repeats and Their Evolutionary Dynamics

A substantial portion of the primate genome is composed of non-coding regions, so-called "dark matter", which includes an abundance of tandemly repeated sequences called satellite DNA. Collectively known as the satellitome, this genomic component offers exciting evolutionary insights into aspects of primate genome biology that raise new questions and challenge existing paradigms. A complete human reference genome was recently reported with telomere-to-telomere human X chromosome assembly that resolved hundreds of dark regions, encompassing a 3.1 Mb centromeric satellite array that had not been identified previously. With the recent exponential increase in the availability of primate genomes, and the development of modern genomic and bioinformatics tools, extensive growth in our knowledge concerning the structure, function, and evolution of satellite elements is expected. The current state of knowledge on this topic is summarized, highlighting various types of primate-specific satellite repeats to compare their proportions across diverse lineages. Inter- and intraspecific variation of satellite repeats in the primate genome are reviewed. The functional significance of these sequences is discussed by describing how the transcriptional activity of satellite repeats can affect gene expression during different cellular processes. Sex-linked satellites are outlined, together with their respective genomic organization. Mechanisms are proposed whereby satellite repeats might have emerged as novel sequences during different evolutionary phases. Finally, the main challenges that hinder the detection of satellite DNA are outlined and an overview of the latest methodologies to address technological limitations is presented.

Hyperphosphorylation of fetal liver IGFBP-1 precedes slowing of fetal growth in nutrient-restricted baboons and may be a mechanism underlying IUGR

In cultured fetal liver cells, insulin-like growth factor (IGF) binding protein (IGFBP)-1 hyperphosphorylation in response to hypoxia and amino acid deprivation is mediated by inhibition of mechanistic target of rapamycin (mTOR) and activation of amino acid response (AAR) signaling and casein kinase (CK)2. We hypothesized that fetal liver mTOR inhibition, activation of AAR and CK2, and IGFBP-1 hyperphosphorylation occur before development of intrauterine growth restriction (IUGR). Pregnant baboons were fed a control (C) or a maternal nutrient restriction (MNR; 70% calories of control) diet starting at gestational day (GD) 30 (term GD 185). Umbilical blood and fetal liver tissue were obtained at GD 120 (C, n = 7; MNR, n = 10) and 165 (C, n = 7; MNR, n = 8). Fetal weights were unchanged at GD 120 but decreased at GD 165 in the MNR group (-13%, P = 0.03). IGFBP-1 phosphorylation, as determined by parallel reaction monitoring mass spectrometry (PRM-MS), immunohistochemistry, and/or Western blot, was enhanced in MNR fetal liver and umbilical plasma at GD 120 and 165. IGF-I receptor autophosphorylation^Tyr1135 (-64%, P = 0.05) was reduced in MNR fetal liver at GD 120. Furthermore, fetal liver CK2 (α/α'/β) expression, CK2β colocalization, proximity with IGFBP-1, and CK2 autophosphorylation^Tyr182 were greater at GD 120 and 165 in MNR vs. C. Additionally, mTOR complex (mTORC)1 (p-P70S6K^Thr389, -52%, P = 0.05) and mTORC2 (p-Akt^Ser473, -56%, P < 0.001) activity were decreased and AAR was activated (p-GCN2^Thr898, +117%, P = 0.02; p-eIF2α^Ser51, +294%, P = 0.002; p-ERK^Thr202, +111%, P = 0.03) in MNR liver at GD 120. Our data suggest that fetal liver IGFBP-1 hyperphosphorylation, mediated by mTOR inhibition and both AAR and CK2 activation, is a key link between restricted nutrient and oxygen availability and the development of IUGR.

A comparison of humans and baboons suggests germline mutation rates do not track cell divisions

In humans, most germline mutations are inherited from the father. This observation has been widely interpreted as reflecting the replication errors that accrue during spermatogenesis. If so, the male bias in mutation should be substantially lower in a closely related species with similar rates of spermatogonial stem cell divisions but a shorter mean age of reproduction. To test this hypothesis, we resequenced two 3-4 generation nuclear families (totaling 29 individuals) of olive baboons (Papio anubis), who reproduce at approximately 10 years of age on average, and analyzed the data in parallel with three 3-generation human pedigrees (26 individuals). We estimated a mutation rate per generation in baboons of 0.57×10-8 per base pair, approximately half that of humans. Strikingly, however, the degree of male bias in germline mutations is approximately 4:1, similar to that of humans-indeed, a similar male bias is seen across mammals that reproduce months, years, or decades after birth. These results mirror the finding in humans that the male mutation bias is stable with parental ages and cast further doubt on the assumption that germline mutations track cell divisions. Our mutation rate estimates for baboons raise a further puzzle, suggesting a divergence time between apes and Old World monkeys of 65 million years, too old to be consistent with the fossil record; reconciling them now requires not only a slowdown of the mutation rate per generation in humans but also in baboons.

Placental fatty acid transport across late gestation in a baboon model of intrauterine growth restriction

KEY POINTS

Intrauterine growth restriction (IUGR) is associated with perinatal morbidity and increased risk of lifelong disease, including neurodevelopmental impairment. Fatty acids (FA) are critical for normal brain development, although their transport across the placenta in IUGR pregnancies is poorly understood. The present study used a baboon model of IUGR (maternal nutrient restriction, MNR) to investigate placental expression of FA transport and binding proteins, and to determine gestational age-related changes in maternal and fetal plasma FA concentrations. We found MNR to be associated with increased placental expression of FA binding and transport proteins in late gestation, with fetal plasma FA concentrations that were similar to those of control animals. The present study is the first to report a profile of fetal and maternal plasma FA concentrations in a baboon model of growth restriction with data that suggest adaptation of placental transport to maintain delivery of critically needed FA.

ABSTRACT

Intrauterine growth restriction (IUGR) is associated with specific changes in placental transport of amino acids, folate and ions. However, little is known about placental fatty acid (FA) transport in IUGR. We hypothesized that placental FA transport proteins (FATP) and FA binding proteins (FABP) are up-regulated and fetal plasma FA concentrations are decreased at term in a baboon model of IUGR. Pregnant baboons were fed control or maternal nutrient restricted (MNR) diet (70% of control calories) from gestation day (GD) 30 (term 184 days). Plasma and placental samples were collected at GD120 (control n = 8, MNR n = 9), GD140 (control n = 6, MNR n = 7) and GD170 (control n = 6, MNR n = 6). Placentas were homogenized, and syncytiotrophoblast microvillous plasma membrane (MVM) and basal plasma membranes (BM) were isolated. Protein expression of FABP1, 3, 4 and 5 (homogenate) and FATP2, 4, and 6 (MVM, BM) was determined by Western blotting. FA content in maternal and umbilical vein plasma was measured by gas chromatography-mass spectrometry. Placental FABP1 and FABP5 expression was increased in MNR compared to controls at GD170, as was MVM FATP2 and FATP6 expression at GD140 and FATP2 expression at GD170. BM FATP4 and FATP6 expression was increased in MNR at GD140. Fetal plasma FA concentrations were similar in controls and MNR. These data suggest the adaptation of placental transport when aiming to maintain delivery of critically needed FAs for fetal growth and brain development.

Analysis of 100 high-coverage genomes from a pedigreed captive baboon colony

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.

Cardiac remodelling in a baboon model of intrauterine growth restriction mimics accelerated ageing

KEY POINTS

Rodent models of intrauterine growth restriction (IUGR) successfully identify mechanisms that can lead to short-term and long-term detrimental cardiomyopathies but differences between rodent and human cardiac physiology and placental-fetal development indicate a need for models in precocial species for translation to human development. We developed a baboon model for IUGR studies using a moderate 30% global calorie restriction of pregnant mothers and used cardiac magnetic resonance imaging to evaluate offspring heart function in early adulthood. Impaired diastolic and systolic cardiac function was observed in IUGR offspring with differences between male and female subjects, compared to their respective controls. Aspects of cardiac impairment found in the IUGR offspring were similar to those found in normal controls in a geriatric cohort. Understanding early cardiac biomarkers of IUGR using non-invasive imaging in this susceptible population, especially taking into account sexual dimorphisms, will aid recognition of the clinical presentation, development of biomarkers suitable for use in humans and management of treatment strategies.

ABSTRACT

Extensive rodent studies have shown that reduced perinatal nutrition programmes chronic cardiovascular disease. To enable translation to humans, we developed baboon offspring cohorts from mothers fed ad libitum (control) or 70% of the control ad libitum diet in pregnancy and lactation, which were growth restricted at birth. We hypothesized that intrauterine growth restriction (IUGR) offspring hearts would show impaired function and a premature ageing phenotype. We studied IUGR baboons (8 male, 8 female, 5.7 years), control offspring (8 male, 8 female, 5.6 years - human equivalent approximately 25 years), and normal elderly (OLD) baboons (6 male, 6 female, mean 15.9 years). Left ventricular (LV) morphology and systolic and diastolic function were evaluated with cardiac MRI and normalized to body surface area. Two-way ANOVA by group and sex (with P < 0.05) indicated ejection fraction, 3D sphericity indices, cardiac index, normalized systolic volume, normalized LV wall thickness, and average filling rate differed by group. Group and sex differences were found for normalized LV wall thickening and normalized myocardial mass, without interactions. Normalized peak LV filling rate and diastolic sphericity index were not correlated in control but strongly correlated in OLD and IUGR baboons. IUGR programming in baboons produces myocardial remodelling, reduces systolic and diastolic function, and results in the emergence of a premature ageing phenotype in the heart. To our knowledge, this is the first demonstration of the specific characteristics of cardiac programming and early life functional decline with ageing in an IUGR non-human primate model. Further studies across the life span will determine progression of cardiac dysfunction.

Population genetic insights into the social organization of Guinea baboons (Papio papio): Evidence for female-biased dispersal

Sex differences in philopatry and dispersal have important consequences on the genetic structure of populations, social groups, and social relationships within groups. Among mammals, male dispersal and female philopatry are most common and closely related taxa typically exhibit similar dispersal patterns. However, among four well-studied species of baboons, only hamadryas baboons exhibit female dispersal, thus differing from their congenerics, which show female philopatry and close-knit female social relationships. Until recently, knowledge of the Guinea baboon social system and dispersal pattern remained sparse. Previous observations suggested that the high degree of tolerance observed among male Guinea baboons could be due to kinship. This led us to hypothesize that this species exhibits male philopatry and female dispersal, conforming to the hamadryas pattern. We genotyped 165 individuals from five localities in the Niokolo-Koba National Park, Senegal, at 14 autosomal microsatellite loci and sequenced a fragment of the mitochondrial hypervariable region I (HVRI) of 55 individuals. We found evidence for higher population structuring in males than in females, as expected if males are the more philopatric sex. A comparison of relatedness between male-male and female-female dyads within and among communities did not yield conclusive results. HVRI diversity within communities was high and did not differ between the sexes, also suggesting female gene flow. Our study is the first comprehensive analysis of the genetic population structure in Guinea baboons and provides evidence for female-biased dispersal in this species. In conjunction with their multilevel social organization, this finding parallels the observations for human hunter-gatherers and strengthens baboons as an intriguing model to elucidate the processes that shaped the highly cooperative societies of Homo.

Cortical Folding of the Primate Brain: An Interdisciplinary Examination of the Genetic Architecture, Modularity, and Evolvability of a Significant Neurological Trait in Pedigreed Baboons (Genus Papio)

Folding of the primate brain cortex allows for improved neural processing power by increasing cortical surface area for the allocation of neurons. The arrangement of folds (sulci) and ridges (gyri) across the cerebral cortex is thought to reflect the underlying neural network. Gyrification, an adaptive trait with a unique evolutionary history, is affected by genetic factors different from those affecting brain volume. Using a large pedigreed population of ∼1000 Papio baboons, we address critical questions about the genetic architecture of primate brain folding, the interplay between genetics, brain anatomy, development, patterns of cortical-cortical connectivity, and gyrification's potential for future evolution. Through Mantel testing and cluster analyses, we find that the baboon cortex is quite evolvable, with high integration between the genotype and phenotype. We further find significantly similar partitioning of variation between cortical development, anatomy, and connectivity, supporting the predictions of tension-based models for sulcal development. We identify a significant, moderate degree of genetic control over variation in sulcal length, with gyrus-shape features being more susceptible to environmental effects. Finally, through QTL mapping, we identify novel chromosomal regions affecting variation in brain folding. The most significant QTL contain compelling candidate genes, including gene clusters associated with Williams and Down syndromes. The QTL distribution suggests a complex genetic architecture for gyrification with both polygeny and pleiotropy. Our results provide a solid preliminary characterization of the genetic basis of primate brain folding, a unique and biomedically relevant phenotype with significant implications in primate brain evolution.

New insights into samango monkey speciation in South Africa

The samango monkey is South Africa's only exclusively forest dwelling primate and represents the southernmost extent of the range of arboreal guenons in Africa. The main threats to South Africa's forests and thus to the samango are linked to increasing land-use pressure and increasing demands for forest resources, resulting in deforestation, degradation and further fragmentation of irreplaceable habitats. The species belongs to the highly polytypic Cercopithecus nictitans group which is sometimes divided into two species C. mitis and C. albogularis. The number of subspecies of C. albogularis is also under debate and is based only on differences in pelage colouration and thus far no genetic research has been undertaken on South African samango monkey populations. In this study we aim to further clarify the number of samango monkey subspecies, as well as their respective distributions in South Africa by combining molecular, morphometric and pelage data. Overall, our study provides the most comprehensive view to date into the taxonomic description of samango monkeys in South Africa. Our data supports the identification of three distinct genetic entities namely; C. a. labiatus, C. a. erythrarchus and C. a. schwarzi and argues for separate conservation management of the distinct genetic entities defined by this study.

Genetic influences on response to novel objects and dimensions of personality in Papio baboons

Behavioral variation within and between populations and species of the genus Papio has been studied extensively, but little is known about the genetic causes of individual- or population-level differences. This study investigates the influence of genetic variation on personality (sometimes referred to as temperament) in baboons and identifies a candidate gene partially responsible for the variation in that phenotype. To accomplish these goals, we examined individual variation in response to both novel objects and an apparent novel social partner (using a mirror test) among pedigreed baboons (n = 578) from the Southwest National Primate Research Center. We investigated the frequency and duration of individual behaviors in response to novel objects and used multivariate factor analysis to identify trait-like dimensions of personality. Exploratory factor analysis identified two distinct dimensions of personality within this population. Factor 1 accounts for 46.8 % of the variance within the behavioral matrix, and consists primarily of behaviors related to the "boldness" of the subject. Factor 2 accounts for 18.8 % of the variation, and contains several "anxiety" like behaviors. Several specific behaviors, and the two personality factors, were significantly heritable, with the factors showing higher heritability than most individual behaviors. Subsequent analyses show that the behavioral reactions observed in the test protocol are associated with animals' social behavior observed later in their home social groups. Finally we used linkage analysis to map quantitative trait loci for the measured phenotypes. Single nucleotide polymorphisms in a positional candidate gene (SNAP25) are associated with variation in one of the personality factors, and CSF levels of homovanillic acid and 3-methoxy-4-hydroxyphenylglycol. This study documents heritable variation in personality among baboons and suggests that sequence variation in SNAP25 may influence differences in behavior and neurochemistry in these nonhuman primates.

Fine-scale genetic assessment of sex-specific dispersal patterns in a multilevel primate society

Like humans, hamadryas baboons (Papio hamadryas) are unusual among primates in having a multilevel social system and stable pair bonds, and are thus a useful model for the evolution of human sociality. While the kinship structure and sex-biased dispersal patterns that underlie human social organization have been extensively elucidated, the impact of these factors on the social system of hamadryas baboons is currently unclear. Here we use genetic analysis of individuals to elucidate the patterns of male and female dispersal across multiple levels of society in a wild population of hamadryas baboons. We characterized 244 members of five hamadryas bands at Filoha, Ethiopia by genotyping one Y-linked and 23 autosomal microsatellite loci and sequencing part of the mitochondrial hypervariable control region I. We found both male and female dispersal to be limited at the level of the band, with more movement of females than males among bands. By integrating long-term behavioral data for Band 1, we also found evidence for male and female philopatry at the clan level. We speculate that male philopatry at the clan level and female dispersal across one-male units and clans may enable both kin-based cooperation among males and the maintenance of kin bonds among females after dispersal. This would mean that, as in humans, kin bonds within both sexes are a core feature of the hamadryas social system, thus contributing to our understanding of the evolution of social organization in humans.

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.

Genetic variations in the beta-tubulin gene and the internal transcribed spacer 2 region of Trichuris species from man and baboons

BACKGROUND

The whipworm Trichuris trichiura has been estimated to infect 604 - 795 million people worldwide. The current control strategy against trichuriasis using the benzimidazoles (BZs) albendazole (400 mg) or mebendazole (500 mg) as single-dose treatment is not satisfactory. The occurrence of single nucleotide polymorphisms (SNPs) in codons 167, 198 or 200 of the beta-tubulin gene has been reported to convey BZ-resistance in intestinal nematodes of veterinary importance. It was hypothesised that the low susceptibility of T. trichiura to BZ could be due to a natural occurrence of such SNPs. The aim of this study was to investigate whether these SNPs were present in the beta-tubulin gene of Trichuris spp. from humans and baboons. As a secondary objective, the degree of identity between T. trichiura from humans and Trichuris spp. from baboons was evaluated based on the beta-tubulin gene and the internal transcribed spacer 2 region (ITS2).

METHODS

Nucleotide sequences of the beta-tubulin gene were generated by PCR using degenerate primers, specific primers and DNA from worms and eggs of T. trichiura and worms of Trichuris spp. from baboons. The ITS2 region was amplified using adult Trichuris spp. from baboons. PCR products were sequenced and analysed. The beta-tubulin fragments were studied for SNPs in codons 167, 198 or 200 and the ITS2 amplicons were compared with GenBank records of T. trichiura.

RESULTS

No SNPs in codons 167, 198 or 200 were identified in any of the analysed Trichuris spp. from humans and baboons. Based on the ITS2 region, the similarity between Trichuris spp. from baboons and GenBank records of T. trichiura was found to be 98 - 99%.

CONCLUSIONS

Single nucleotide polymorphisms in codon 167, 198 and 200, known to confer BZ-resistance in other nematodes, were absent in the studied material. This study does not provide data that could explain previous reports of poor BZ treatment efficacy in terms of polymorphism in these codons of beta-tubulin. Based on a fragment of the beta-tubulin gene and the ITS2 region sequenced, it was found that T. trichiura from humans and Trichuris spp. isolated from baboons are closely related and may be the same species.

The baboon kidney transcriptome: analysis of transcript sequence, splice variants, and abundance

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.

Effects of maternal nutrient restriction, intrauterine growth restriction, and glucocorticoid exposure on phosphoenolpyruvate carboxykinase-1 expression in fetal baboon hepatocytes in vitro

BACKGROUND

The objective of this study was to develop a cell culture system for fetal baboon hepatocytes and to test the hypotheses that (i) expression of the gluconeogenic enzyme phosphoenolpyruvate carboxykinase-1 (PEPCK-1) is upregulated in hepatocytes isolated from fetuses of nutrient-restricted mothers (MNR) compared with ad libitum-fed controls (CTR), and (ii) glucocorticoids stimulate PEPCK-1 expression.

METHODS

Hepatocytes from 0.9G CTR and MNR fetuses were isolated and cultured. PEPCK-1 protein and mRNA levels in hepatocytes were determined by Western blot and quantitative PCR, respectively.

RESULTS

Fetuses of MNR mothers were intrauterine growth restricted (IUGR). Feasibility of culturing 0.9G fetal baboon hepatocytes was demonstrated. PEPCK-1 protein levels were increased in hepatocytes isolated from IUGR fetuses, and PEPCK-1 mRNA expression was stimulated by glucocorticoids in fetal hepatocytes.

CONCLUSIONS

Cultured fetal baboon hepatocytes that retain their in vivo phenotype provide powerful in vitro tools to investigate mechanisms that regulate normal and programmed hepatic function.

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.

Pleiotropic Effects of Genes for Insulin Resistance on Adiposity in Baboons

OBJECTIVE

Previous research has suggested a genetic contribution to the development of insulin resistance and obesity. We hypothesized that the same genes influencing insulin resistance might also contribute to the variation in adiposity.

RESEARCH METHODS AND PROCEDURES

A total of 601 (200 male, 401 female) adult baboons (Papio hamadryas) from nine families with pedigrees ranging in size from 43 to 121 were used in this study. Plasma insulin, glucose, C-peptide, and adiponectin were analyzed, and homeostasis model assessment of insulin resistance (HOMA IR) was calculated. Fat biopsies were collected from omental fat tissue, and triglyceride concentration per gram of fat tissue was determined. Body weight and length were measured, and BMI was derived. Univariate and bivariate quantitative genetic analyses were performed using SOLAR.

RESULTS

Insulin, glucose, C-peptide, and adiponectin levels, HOMA IR, triglyceride concentration of fat tissue, body weight, and BMI were all found to be significantly heritable, with heritabilities ranging from 0.15 to 0.80. Positive genetic correlations (r(G)s) were observed for HOMA IR with C-peptide (r(G) = 0.88 +/- 0.10, p = 0.01), triglyceride concentration in fat tissue (r(G) = 0.86 +/- 0.33, p = 0.02), weight (r(G) = 0.50 +/- 0.20, p = 0.03), and BMI (r(G) = 0.64 +/- 0.22, p = 0.02).

DISCUSSION

These results suggest that a set of genes contributing to insulin resistance also influence general and central adiposity phenotypes. Further genetic research in a larger sample size is needed to identify the common genes that constitute the genetic basis for the development of insulin resistance and obesity.

Adiponectin But Not Resistin Is Associated with Insulin Resistance-Related Phenotypes in Baboons

OBJECTIVE

The hormones adiponectin and resistin have been associated with insulin resistance. This paper analyzed the potential relationship between adiponectin and resistin and insulin resistance-related phenotypes in baboons.

RESEARCH METHODS AND PROCEDURES

One hundred eight adult baboons (84 female and 24 male) were studied. Weight was measured, and a blood sample was collected under fasting conditions for plasma and monocyte isolation. Fasting glucose, insulin, C-peptide, and adiponectin levels in plasma were measured by standard methods. Insulin resistance was calculated by the homeostasis model assessment index. Resistin mRNA abundance in monocytes was determined by real-time quantitative reverse transcription-polymerase chain reaction. Data were clustered by weight tertiles for statistical analysis.

RESULTS

As observed in humans, the insulin resistance-related phenotypes were related to weight, plasma levels of adiponectin, and C-peptide. No significant relationship between resistin circulating levels or expression in monocytes and insulin resistance-related phenotypes was found in baboons.

DISCUSSION

These findings suggest that resistin is not associated with insulin resistance. However, previous observations of relationships among weight, adiponectin, and insulin resistance are confirmed.

Trichuris sp. and Strongyloides sp. infections in a free-ranging baboon colony

We conducted cross-sectional surveys of parasites infecting a large free-living colony of baboons at the Southwest National Primate Research Center in San Antonio in October 2003 and April 2004, immediately before, and 6 mo after, treatment with ivermectin. Trichuris sp. was the predominant species present, infecting 79 and 69% of individual animals in the 2 surveys, with fecal egg counts (FEC) of up to 60,200 eggs per g (epg) (mean = 1,235 in October 2003 and 1,256 in April 2004). Prevalence remained fairly stable across age groups, and intensity was highest in animals <1 or >15 yr old, in contrast to patterns observed in humans, where school-age children show the heaviest infections. Strongyloides sp. was also identified, but the species identity remains uncertain. Small subunit ribosomal DNA sequences differed from published sequences of Strongyloides fuelleborni at multiple sites, but resided in a monophyletic group with other Strongyloides species with 92% bootstrap support. This may reflect a recent acquisition from a local host, or that the published sequence of S. fuelleborni is incorrect. Widespread infections with 2 nematode genera in a free-ranging baboon colony that are an important source of morbidity in human populations provide a useful model system for work on the epidemiology, control, pathology, and genetics of these parasites in a host species that is physiologically, immunologically, and genetically similar to humans.

A genetic linkage map and comparative mapping of the prairie vole (Microtus ochrogaster) genome

BACKGROUND

The prairie vole (Microtus ochrogaster) is an emerging rodent model for investigating the genetics, evolution and molecular mechanisms of social behavior. Though a karyotype for the prairie vole has been reported and low-resolution comparative cytogenetic analyses have been done in this species, other basic genetic resources for this species, such as a genetic linkage map, are lacking.

RESULTS

Here we report the construction of a genome-wide linkage map of the prairie vole. The linkage map consists of 406 markers that are spaced on average every 7 Mb and span an estimated ~90% of the genome. The sex average length of the linkage map is 1707 cM, which, like other Muroid rodent linkage maps, is on the lower end of the length distribution of linkage maps reported to date for placental mammals. Linkage groups were assigned to 19 out of the 26 prairie vole autosomes as well as the X chromosome. Comparative analyses of the prairie vole linkage map based on the location of 387 Type I markers identified 61 large blocks of synteny with the mouse genome. In addition, the results of the comparative analyses revealed a potential elevated rate of inversions in the prairie vole lineage compared to the laboratory mouse and rat.

CONCLUSIONS

A genetic linkage map of the prairie vole has been constructed and represents the fourth genome-wide high-resolution linkage map reported for Muroid rodents and the first for a member of the Arvicolinae sub-family. This resource will advance studies designed to dissect the genetic basis of a variety of social behaviors and other traits in the prairie vole as well as our understanding of genome evolution in the genus Microtus.

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.

Genetic and environmental contributions to variation in baboon cranial morphology

The development, function, and integration of morphological characteristics are all hypothesized to influence the utility of traits for phylogenetic reconstruction by affecting the way in which morphological characteristics evolve. We use a baboon model to test the hypotheses about phenotypic and quantitative genetic variation of traits in the cranium that bear on a phenotype's propensity to evolve. We test the hypotheses that: 1) individual traits in different functionally and developmentally defined regions of the cranium are differentially environmentally, genetically, and phenotypically variable; 2) genetic covariance with other traits constrains traits in one region of the cranium more than those in others; 3) and regions of the cranium subject to different levels of mechanical strain differ in the magnitude of variation in individual traits. We find that the levels of environmental and genetic variation in individual traits are randomly distributed across regions of the cranium rather than being structured by developmental origin or degree of exposure to strain. Individual traits in the cranial vault tend to be more constrained by covariance with other traits than those in other regions. Traits in regions subject to high degrees of strain during mastication are not any more variable at any level than other traits. If these results are generalizable to other populations, they indicate that there is no reason to suppose that individual traits from any one part of the cranium are intrinsically less useful for reconstructing patterns of evolution than those from any other part.

Evolutionary genetics in wild primates: combining genetic approaches with field studies of natural populations

Ecological and evolutionary studies of wild primates hold important keys to understanding both the shared characteristics of primate biology and the genetic and phenotypic differences that make specific lineages, including our own, unique. Although complementary genetic research on nonhuman primates has long been of interest, recent technological and methodological advances now enable functional and population genetic studies in an unprecedented manner. In the past several years, novel genetic data sets have revealed new information about the demographic history of primate populations and the genetics of adaptively important traits. In combination with the rich history of behavioral, ecological, and physiological work on natural primate populations, genetic approaches promise to provide a compelling picture of primate evolution in the past and in the present day.

Assisted Reproductive Technologies (ART) with baboons generate live offspring: a nonhuman primate model for ART and reproductive sciences

Human reproduction has benefited significantly by investigating nonhuman primate (NHP) models, especially rhesus macaques. To expand the Old World monkey species available for human reproductive studies, we present protocols in baboons, our closest Old World primate relatives, for assisted reproductive technologies (ART) leading to live born offspring. Baboons complement rhesus by confirming or modifying observations generated in humans often obtained by the study of clinically discarded specimens donated by anonymous infertility patient couples. Here, baboon ART protocols, including oocyte collection, in vitro fertilization (IVF), intracytoplasmic sperm injection (ICSI), preimplantation development to blastocyst stage, and embryo transfer techniques are described. With baboon ART methodologies in place, motility during baboon fertilization was investigated by time-lapse video microscopy (TLVM). The first ART baboons produced by ICSI, a pair of male twins, were delivered naturally at 165 days postgestation. Genetic testing of these twins confirmed their ART parental origins and demonstrated that they are unrelated fraternal twins not identicals. These results have implications for ART outcomes, embryonic stem cell (ESC) derivation, and reproductive sciences.

Identification of baboon microRNAs expressed in liver and lymphocytes

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.

Mapping primary gyrogenesis during fetal development in primate brains: high-resolution in utero structural MRI of fetal brain development in pregnant baboons

The global and regional changes in the fetal cerebral cortex in primates were mapped during primary gyrification (PG; weeks 17-25 of 26 weeks total gestation). Studying pregnant baboons using high-resolution MRI in utero, measurements included cerebral volume, cortical surface area, gyrification index and length and depth of 10 primary cortical sulci. Seven normally developing fetuses were imaged in two animals longitudinally and sequentially. We compared these results to those on PG that from the ferret studies and analyzed them in the context of our recent studies of phylogenetics of cerebral gyrification. We observed that in both primates and non-primates, the cerebrum undergoes a very rapid transformation into the gyrencephalic state, subsequently accompanied by an accelerated growth in brain volume and cortical surface area. However, PG trends in baboons exhibited some critical differences from those observed in ferrets. For example, in baboons, the growth along the long (length) axis of cortical sulci was unrelated to the growth along the short (depth) axis and far outpaced it. Additionally, the correlation between the rate of growth along the short sulcal axis and heritability of sulcal depth was negative and approached significance (r = -0.60; p < 0.10), while the same trend for long axis was positive and not significant (p = 0.3; p = 0.40). These findings, in an animal that shares a highly orchestrated pattern of PG with humans, suggest that ontogenic processes that influence changes in sulcal length and depth are diverse and possibly driven by different factors in primates than in non-primates.

The evolution of heterochiasmy: the role of sexual selection and sperm competition in determining sex-specific recombination rates in eutherian mammals

Early karyotypic work revealed that female and male recombination rates in many species show pronounced differences, and this pattern of heterochiasmy has also been observed in modern linkage mapping studies. Several hypotheses to explain this phenomenon have been offered, ranging from strictly biological mechanisms related to the gametic differences between the sexes, to more evolutionary models based on sexually antagonistic selection. However, despite the long history of interest in heterochiasmy, empirical data has failed to support any theory or pattern consistently. Here I test two alternative evolutionary hypotheses regarding heterochiasmy across the eutherian mammals, and show that sexual dimorphism, but not sperm competition, is strongly correlated with recombination rate, suggesting that sexual antagonism is an important influence. However, the observed relationship between heterochiasmy and sexual dimorphism runs counter to theoretical predictions, with male recombination higher in species with high levels of sexual dimorphism. This may be the response to male-biased dispersal, which, rather than the static male fitness landscape envisioned in the models tested here, could radically shift optimal male fitness parameters among generations.

Molecular cloning of the baboon UDP-glucuronosyltransferase 2B gene family and their activity in conjugating morphine

Glucuronidation by UDP-glucuronyltransferase 2B enzymes (UGT2Bs) is a major pathway for the elimination of endobiotics and xenobiotics, including therapeutic drugs. Morphine, a probe drug for UGT2B7, is metabolized to morphine-3-beta-glucuronide (M3G) and morphine-6-beta-glucuronide (M6G) in humans. Morphine has been used in a series of experiments in the baboon to characterize developmental changes in fetal glucuronidation. This study identifies the baboon UGT2B family of enzymes, compares them with that of the human and the monkey (Macaca fascicularis), and measures the activity of the individual baboon UGT2Bs toward morphine. UGT2B cDNAs were cloned from the liver of adult and newborn baboons and expressed in human embryonic kidney 293 cells. The UGT activity toward morphine was assessed by the rate of formation of M3G and M6G by high-performance liquid chromatography. Eight baboon UGT2Bs were cloned and identified: UGT2B41 and UGT2B42, which are 90% homologous to human UGT2B4; UGT2B43, which is 93% homologous to human UGT2B15; and UGT2B39, UGT2B40, UGT2B44, UGT2B45, and UGT2B46, which are 89 to 91% homologous to human UGT2B7. Homology between baboon and monkey UGT2B ranged from 92.6 to 99.1%, with the primary protein structure of UGT2B43 being 99.1% identical to monkey UGT2B20, including a unique R96I substitution. Gene conversion interfered with the phylogenetic signal in the baboon UGT2B7-like and the monkey UGT2B4-like groups and led to concerted evolution of these enzymes. All of the baboon UGT2Bs metabolized morphine to both M3G and M6G. This study lays the foundation for investigating the regulation of UGT2B enzymes during fetal and neonatal development in the baboon.

Quantitative loci regulating plasma levels of gamma glutamyl transferase and albumin and their genetic correlations with cardiovascular risk factors

gamma Glutamyl transferase (GGT) and albumin (ALB) are two markers of liver function. These two proteins have been associated with non-alcoholic fatty liver disease and cardiovascular disease. The objectives of this study were to explore the genetic factors that influence variation in the plasma levels of GGT and ALB and to evaluate their genetic correlations with cardiovascular risk factors. Baboons from the Southwest National Primate Research Center at the Southwest Foundation for Biomedical Research, San Antonio, TX, were used as an animal model. The baboons were fed a standard monkey chow diet ad libitum. Fasting plasma concentrations of GGT, ALB, triglycerides, total cholesterol and LDL cholesterol were measured in 350 pedigreed adult baboons by standard assay procedures. A maximum likelihood-based variance decomposition approach implemented in the computer program SOLAR was used to conduct genetic analyses. The heritabilities of GGT (h(2) = 0.55; P < 0.0001) and ALB (h(2) = 0.42; P < 0.01) were significant. No statistically significant associations were found between GGT and the cardiovascular-related phenotypes. Genetic correlations between ALB and total cholesterol, LDL cholesterol and triglycerides were significant. A QTL (LOD = 2.8) for GGT plasma levels was identified on the baboon homologue of human chromosome 22 between markers D22S304 and D22S280. A QTL (LOD = 2.3) near marker D10S1432 was detected on the baboon homologue of human chromosome 10 for ALB. These results imply that variations in the plasma levels of GGT and ALB are under significant genetic regulation and that a common genetic component influences ALB and cardiovascular risk factor phenotypes.

Integration of genetic and genomic methods for identification of genes and gene variants encoding QTLs in the nonhuman primate

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.

LRP5 sequence and polymorphisms in the baboon

BACKGROUND

LRP5 is known to have an important relationship with bone density and a variety of other biological processes. Mapping to human chromosome 11q13.2, LRP5 shows considerable evolutionary conservation. Orthologs of this gene exist in many species, although comparison of human LRP5 with other non-human primates has not been performed until now.

METHODS

We reported the complementary DNA (cDNA) sequence and deduced amino acid sequence for baboon LRP5, and compared the baboon and human sequences. cDNA sequences for 21 baboons were examined to identify single-nucleotide polymorphisms (SNPs).

RESULTS

Sequences of coding regions in human and baboon LRP5 showed 97- 99% homology. Twenty-five SNPs were identified in the coding region of baboon LRP5.

CONCLUSION

The observed degree of coding sequence homology in LRP5 led us to expect that the baboon may serve as a useful model for future research into the role(s) of this gene in primate metabolic diseases.

Chromosomal mapping of canine-derived BAC clones to the red fox and American mink genomes

High-quality sequencing of the dog (Canis lupus familiaris) genome has enabled enormous progress in genetic mapping of canine phenotypic variation. The red fox (Vulpes vulpes), another canid species, also exhibits a wide range of variation in coat color, morphology, and behavior. Although the fox genome has not yet been sequenced, canine genomic resources have been used to construct a meiotic linkage map of the red fox genome and begin genetic mapping in foxes. However, a more detailed gene-specific comparative map between the dog and fox genomes is required to establish gene order within homologous regions of dog and fox chromosomes and to refine breakpoints between homologous chromosomes of the 2 species. In the current study, we tested whether canine-derived gene-containing bacterial artificial chromosome (BAC) clones can be routinely used to build a gene-specific map of the red fox genome. Forty canine BAC clones were mapped to the red fox genome by fluorescence in situ hybridization (FISH). Each clone was uniquely assigned to a single fox chromosome, and the locations of 38 clones agreed with cytogenetic predictions. These results clearly demonstrate the utility of FISH mapping for construction of a whole-genome gene-specific map of the red fox. The further possibility of using canine BAC clones to map genes in the American mink (Mustela vison) genome was also explored. Much lower success was obtained for this more distantly related farm-bred species, although a few BAC clones were mapped to the predicted chromosomal locations.

Development of an integrative database with 499 novel microsatellite markers for Macaca fascicularis

BACKGROUND

Cynomolgus macaques (Macaca fascicularis) are a valuable resource for linkage studies of genetic disorders, but their microsatellite markers are not sufficient. In genetic studies, a prerequisite for mapping genes is development of a genome-wide set of microsatellite markers in target organisms. A whole genome sequence and its annotation also facilitate identification of markers for causative mutations. The aim of this study is to establish hundreds of microsatellite markers and to develop an integrative cynomolgus macaque genome database with a variety of datasets including marker and gene information that will be useful for further genetic analyses in this species.

RESULTS

We investigated the level of polymorphisms in cynomolgus monkeys for 671 microsatellite markers that are covered by our established Bacterial Artificial Chromosome (BAC) clones. Four hundred and ninety-nine (74.4%) of the markers were found to be polymorphic using standard PCR analysis. The average number of alleles and average expected heterozygosity at these polymorphic loci in ten cynomolgus macaques were 8.20 and 0.75, respectively.

CONCLUSION

BAC clones and novel microsatellite markers were assigned to the rhesus genome sequence and linked with our cynomolgus macaque cDNA database (QFbase). Our novel microsatellite marker set and genomic database will be valuable integrative resources in analyzing genetic disorders in cynomolgus macaques.

Localization of multiple pleiotropic genes for lipoprotein metabolism in baboons

We employed a novel approach to identify the key loci that harbor genes influencing lipoprotein metabolism in approximately 2,000 pedigreed baboons fed various diets differing in levels of fat and cholesterol. In this study, 126 overlapping traits related to both LDL and HDL metabolism were normalized and subjected to genome-wide linkage screening. As was expected, the traits were highly, but not completely, correlated. We exploited the information in these correlated traits by focusing on those genomic regions harboring quantitative trait loci (QTL) for multiple traits, reasoning that the more influential genes would impact a larger number of traits. This study identified five major QTL clusters (each with at least two significant logarithm of the odds scores >4.7), two of which had not been previously reported in baboons. One of these mapped to the baboon ortholog of human chromosome 1p32-p34 and influenced concentrations of LDL-cholesterol on Basal and high-fat, low-cholesterol diets. The other novel QTL cluster mapped to the baboon ortholog of human chromosome 12q13.13-q14.1 and influenced LDL size properties on high-fat, low-cholesterol and high-fat, high-cholesterol, but not Basal, diets. Confirming the value of this approach, three of the QTL clusters replicated published linkage findings for the same or similar traits.

Quantitative Genetics, Pleiotropy, and Morphological Integration in the Dentition of Papio hamadryas

Variation in the mammalian dentition is highly informative of adaptations and evolutionary relationships, and consequently has been the focus of considerable research. Much of the current research exploring the genetic underpinnings of dental variation can trace its roots to Olson and Miller's 1958 book Morphological Integration. These authors explored patterns of correlation in the post-canine dentitions of the owl monkey and Hyopsodus, an extinct condylarth from the Eocene. Their results were difficult to interpret, as was even noted by the authors, due to a lack of genetic information through which to view the patterns of correlation. Following in the spirit of Olson and Miller's research, we present a quantitative genetic analysis of dental variation in a pedigreed population of baboons. We identify patterns of genetic correlations that provide insight to the genetic architecture of the baboon dentition. This genetic architecture indicates the presence of at least three modules: an incisor module that is genetically independent of the post-canine dentition, and a premolar module that demonstrates incomplete pleiotropy with the molar module. We then compare this matrix of genetic correlations to matrices of phenotypic correlations between the same measurements made on museum specimens of another baboon subspecies and the Southeast Asian colobine Presbytis. We observe moderate significant correlations between the matrices from these three primate taxa. From these observations we infer similarity in modularity and hypothesize a common pattern of genetic integration across the dental arcade in the Cercopithecoidea.

Genetic variation in baboon craniofacial sexual dimorphism

Sexual dimorphism is a widespread phenomenon and contributes greatly to intraspecies variation. Despite a long history of active research, the genetic basis of dimorphism for complex traits remains unknown. Understanding the sex-specific differences in genetic architecture for cranial traits in a highly dimorphic species could identify possible mechanisms through which selection acts to produce dimorphism. Using distances calculated from three-dimensional landmark data from CT scans of 402 baboon skulls from a known genealogy, we estimated genetic variance parameters in both sexes to determine the presence of gene-by-sex (G x S) interactions and X-linked heritability. We hypothesize that traits exhibiting the greatest degree of sexual dimorphism (facial traits in baboons) will demonstrate either stronger G x S interactions or X-linked effects. We found G x S interactions and X-linked effects for a few measures that span the areas connecting the face to the neurocranium but for no traits restricted to the face. This finding suggests that facial traits will have a limited response to selection for further evolution of dimorphism in this population. We discuss the implications of our results with respect to the origins of cranial sexual dimorphism in this baboon sample, and how the genetic architecture of these traits affects their potential for future evolution.

Polymorphic microsatellite loci for the common marmoset (Callithrix jacchus) designed using a cost- and time-efficient method

We describe a cost- and time-efficient method for designing new microsatellite markers in any species with substantial genomic DNA sequence data available. Using this technique, we report 14 new polymorphic dinucleotide microsatellite loci isolated from the common marmoset. The relative yield of new polymorphisms was higher with less labor than described in previous marmoset studies. Of 20 loci initially evaluated, 14 were polymorphic and amplified reliably (70% success rate). The number of alleles ranged from 3 to 9 with heterozygosity varying from 0.48 to 0.83.

A genomewide linkage scan for quantitative trait loci influencing the craniofacial complex in baboons (Papio hamadryas spp.)

Numerous studies have detected significant contributions of genes to variation in development, size, and shape of craniofacial traits in a number of vertebrate taxa. This study examines 43 quantitative traits derived from lateral cephalographs of 830 baboons (Papio hamadryas) from the pedigreed population housed at the Southwest National Primate Research Center. Quantitative genetic analyses were conducted using the SOLAR 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, 14 significant quantitative trait loci (QTL) were identified for 12 traits from the three developmental components (basicranium, splanchnocranium, and neurocranium) of the craniofacial complex. These QTL were found on baboon chromosomes (and human orthologs) PHA1 (HSA1), PHA 2 (HSA3), PHA4 (HSA6), PHA11 (HSA12), PHA13 (HSA2), PHA16 (HSA17), and PHA17 (HSA13) (PHA, P. hamadryas; HSA, Homo sapiens). This study of the genetic architecture of the craniofacial complex in baboons provides the groundwork needed to establish the baboon as an animal model for the study of genetic and nongenetic influences on craniofacial variation.

Cross-species replication of a resistin mRNA QTL, but not QTLs for circulating levels of resistin, in human and baboon

Resistin has been associated with inflammation and risk for cardiovascular disease. We previously reported evidence of a QTL on chromosome 19p13 affecting the abundance of resistin (RETN) mRNA in the omental adipose tissue of baboons (L0D score 3.8). In this study, whole genome transcription levels were assessed in human lymphocyte samples from 1240 adults participating in the San Antonio Family Heart Study, using the Sentrix Human-6 Expression Beadchip. Lymphocytes were surveyed, as it has been proposed that their expression levels may reflect those in harder to ascertain tissues, such as adipose tissue, that are thought to be more directly relevant to disease procesn was conducted to detect loci affecting RETN mRNA levels. We obtained significant evidence for a QTL influencing the RETN expression (LOD score 10.7) on chromosome 19p. This region is orthologous/homologous to the region previously localized on baboon chromosome 19. The strongest positional candidate gene in this region is the structural gene for resistin, itself. We also found evidence for a QTL influencing resistin protein levels (LOD score 5.3) on chromosome 14q. This differs from our previously reported QTL on chromosome 18 in baboons. The different QTLs for circulating protein suggests that post-translational processing and turnover may be influenced by different or multiple genes in baboons and humans. The parallel findings of a cis-eQTL for RETN mRNA in baboon omental tissue and human lymphocytes lends support to the strategy of using lymphocyte gene expression levels as a surrogate for gene expression levels in other tissues.