A second-generation genetic linkage map of the baboon (Papio hamadryas) genome

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.

Genetic regulatory effects in response to a high-cholesterol, high-fat diet in baboons

Steady-state expression quantitative trait loci (eQTLs) explain only a fraction of disease-associated loci identified through genome-wide association studies (GWASs), while eQTLs involved in gene-by-environment (GxE) interactions have rarely been characterized in humans due to experimental challenges. Using a baboon model, we found hundreds of eQTLs that emerge in adipose, liver, and muscle after prolonged exposure to high dietary fat and cholesterol. Diet-responsive eQTLs exhibit genomic localization and genic features that are distinct from steady-state eQTLs. Furthermore, the human orthologs associated with diet-responsive eQTLs are enriched for GWAS genes associated with human metabolic traits, suggesting that context-responsive eQTLs with more complex regulatory effects are likely to explain GWAS hits that do not seem to overlap with standard eQTLs. Our results highlight the complexity of genetic regulatory effects and the potential of eQTLs with disease-relevant GxE interactions in enhancing the understanding of GWAS signals for human complex disease using non-human primate models.

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.

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.

Selection against admixture and gene regulatory divergence in a long-term primate field study

Genetic admixture is central to primate evolution. We combined 50 years of field observations of immigration and group demography with genomic data from ~9 generations of hybrid baboons to investigate the consequences of admixture in the wild. Despite no obvious fitness costs to hybrids, we found signatures of selection against admixture similar to those described for archaic hominins. These patterns were concentrated near genes where ancestry is strongly associated with gene expression. Our analyses also show that introgression is partially predictable across the genome. This study demonstrates the value of integrating genomic and field data for revealing how "genomic signatures of selection" (e.g., reduced introgression in low-recombination regions) manifest in nature; moreover, it underscores the importance of other primates as living models for human evolution.

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.

Regulation of Uterine Spiral Artery Remodeling: a Review

Extravillous trophoblast remodeling of the uterine spiral arteries is essential for promoting blood flow to the placenta and fetal development, but little is known about the regulation of this process. A defect in spiral artery remodeling underpins adverse conditions of human pregnancy, notably early-onset preeclampsia and fetal growth restriction, which result in maternal and fetal morbidity and mortality. Many in vitro studies have been conducted to determine the ability of growth and other factors to stimulate trophoblast cells to migrate across a synthetic membrane. Clinical studies have investigated whether the maternal levels of various factors are altered during abnormal human pregnancy. Animal models have been established to assess the ability of various factors to recapitulate the pathophysiological symptoms of preeclampsia. This review analyzes the results of the in vitro, clinical, and animal studies and describes a nonhuman primate experimental paradigm of defective uterine artery remodeling to study the regulation of vessel remodeling.

Reduced meiotic recombination in rhesus macaques and the origin of the human recombination landscape

Characterizing meiotic recombination rates across the genomes of nonhuman primates is important for understanding the genetics of primate populations, performing genetic analyses of phenotypic variation and reconstructing the evolution of human recombination. Rhesus macaques (Macaca mulatta) are the most widely used nonhuman primates in biomedical research. We constructed a high-resolution genetic map of the rhesus genome based on whole genome sequence data from Indian-origin rhesus macaques. The genetic markers used were approximately 18 million SNPs, with marker density 6.93 per kb across the autosomes. We report that the genome-wide recombination rate in rhesus macaques is significantly lower than rates observed in apes or humans, while the distribution of recombination across the macaque genome is more uniform. These observations provide new comparative information regarding the evolution of recombination in primates.

Vascular Endothelial Growth Factor Delivery to Placental Basal Plate Promotes Uterine Artery Remodeling in the Primate

Extravillous trophoblast (EVT) uterine artery remodeling (UAR) promotes placental blood flow, but UAR regulation is unproven. Elevating estradiol (E2) in early baboon pregnancy suppressed UAR and EVT vascular endothelial growth factor (VEGF) expression, but this did not prove that VEGF mediated this process. Therefore, our primate model of prematurely elevating E2 and contrast-enhanced ultrasound cavitation of microbubble (MB) carriers was used to deliver VEGF DNA to the placental basal plate (PBP) to establish the role of VEGF in UAR. Baboons were treated on days 25 to 59 of gestation (term, 184 days) with E2 alone or with E2 plus VEGF DNA-conjugated MBs briefly infused via a maternal peripheral vein on days 25, 35, 45, and 55. At each of these times an ultrasound beam was directed to the PBP to collapse the MBs and release VEGF DNA. VEGF DNA-labeled MBs per contrast agent was localized in the PBP but not the fetus. Remodeling of uterine arteries >25 µm in diameter on day 60 was 75% lower (P < 0.001) in E2-treated (7% ± 2%) than in untreated baboons (30% ± 4%) and was restored to normal by E2/VEGF. VEGF protein levels (signals/nuclear area) within the PBP were twofold lower (P < 0.01) in E2-treated (4.2 ± 0.9) than in untreated (9.8 ± 2.8) baboons and restored to normal by E2/VEGF (11.9 ± 1.6), substantiating VEGF transfection. Thus, VEGF gene delivery selectively to the PBP prevented the decrease in UAR elicited by prematurely elevating E2 levels, establishing the role of VEGF in regulating UAR in vivo during primate pregnancy.

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.

Identification of coordinately regulated microRNA-gene networks that differ in baboons discordant for LDL-cholesterol

RATIONALE

Plasma 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 RESULTS

Using 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.

CONCLUSIONS

We 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.

Genomic Tools for the Use of Nonhuman Primates in Translational Research

Nonhuman primates (NHPs) are important preclinical models for understanding the etiology of human diseases and for developing therapies and vaccines to cure or eliminate disease. Most human diseases have genetic components. Therefore, to be of maximal utility, the NHP species used for translational science should be as well characterized in regard to their genome and transcriptome as possible. This article reviews the current status of genomic information for the five NHP species used most often in translational research: rhesus macaque, cynomolgus macaque, vervet (African green) monkey, baboon, and marmoset NHP. These species have published whole genome sequences (with the exception of the baboon) and relatively well-characterized transcriptomes. Some have also been characterized in regard to specific genetic loci that are particularly related to translational concerns, such as the major histocompatability complex and the cytochrome P40 genes. Genomic resources to aid in stratifying captive populations in regard to genetic and phenotypic characteristics have been developed as an aid to enhancing reproducibility and facilitating more efficient use of animals. Taken together, the current genomic resources and numerous studies currently underway to improve them should enhance the value of NHPs as preclinical models of human disease.

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.

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.

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.

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.

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.

Differential microRNA response to a high-cholesterol, high-fat diet in livers of low and high LDL-C baboons

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.

Genetical genomics of Th1 and Th2 immune response in a baboon model of atherosclerosis risk factors

OBJECTIVE

CD4(+) T-cells mediate inflammation in atherosclerosis, but additive genetic effects on associated pathways of Th1 and Th2 immune response have not been described. We sought to characterize heritability, pleiotropy, and QTL effects on the expression of genes implicated in Th1 and Th2 immune response in a baboon model of risk factors for atherosclerosis.

METHODS

We employed a maximum likelihood-based variance decomposition approach to estimate additive genetic effects on transcript levels generated from a gene expression profile of lymphocytes in 499 pedigreed baboons maintained on a basal diet. Transcript levels for 57 genes implicated in Th1 and Th2 immune response were selected for analysis based on significant heritability in this profile. Multipoint whole genome scans were conducted on heritable transcript levels to localize QTLs influencing these measures. To evaluate pleiotropic effects on transcript levels, we estimated genetic and phenotypic correlations among transcript measures, and assessed their correspondence using a Mantel test. Network analysis using GeneGo's MetaCore™ software was conducted to characterize known interaction among coded proteins.

RESULTS

Heritabilities for candidate gene transcript levels ranged from 0.092-0.786 (median h(2)=0.278, P=4.72×10(-4)). Linkage analyses yielded significant evidence (LOD≥2.73) for 14 eQTLs (LOD score range 2.76-14.87, genome-wide P=4.9×10(-2)-1.03×10(-14)). Estimates of genetic correlation supported shared additive genetic effects incorporating all 57 transcripts (null hypothesis of ρ(G)=0 rejected at FDR≤0.05 for 522 of 1596 estimates), and accounted for most of the observed phenotypic correlation among transcripts (Mantel test, r([ρP],)([ρG])=0.781, P<0.0001). Network analysis revealed direct interactions among 54 of the 57 coded proteins.

CONCLUSIONS

We conclude that major genetic effects influence expression levels of multiple genes implicated in Th1 and Th2 immune response. Additionally, we find that expression levels of these candidate genes are characterized by extensive pleiotropy, consistent with known interaction among their coded proteins, many of which are independently associated with atherosclerosis.

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.

Identification of a QTL for adipocyte volume and of shared genetic effects with aspartate aminotransferase

Plasma levels of aspartate aminotransferase (AST), a liver enzyme, are elevated in patients with visceral obesity. This study examined whether adipocyte volume is under the influence of genetic factors and evaluated its genetic correlations with AST. Fasting plasma levels of 344 pedigreed baboons from the Southwest National Primate Research Center in San Antonio, TX, USA, were assayed for AST. Adipocyte volume was measured using biopsies of omental adipose tissue. Adipocyte volume, body weight, and plasma AST were heritable. Genetic correlations between the measured adiposity-related phenotypes and AST were significant. A quantitative trait locus (LOD score 3.2) for adipocyte volume was identified on the baboon homolog of human chromosome 6 near marker D6S1028. These results suggest that omental adipocyte volume is under genetic regulation and that shared genetic factors influence adiposity-associated traits and AST.

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.

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.

Sequencing primate genomes: what have we learned?

We summarize the progress in whole-genome sequencing and analyses of primate genomes. These emerging genome datasets have broadened our understanding of primate genome evolution revealing unexpected and complex patterns of evolutionary change. This includes the characterization of genome structural variation, episodic changes in the repeat landscape, differences in gene expression, new models regarding speciation, and the ephemeral nature of the recombination landscape. The functional characterization of genomic differences important in primate speciation and adaptation remains a significant challenge. Limited access to biological materials, the lack of detailed phenotypic data and the endangered status of many critical primate species have significantly attenuated research into the genetic basis of primate evolution. Next-generation sequencing technologies promise to greatly expand the number of available primate genome sequences; however, such draft genome sequences will likely miss critical genetic differences within complex genomic regions unless dedicated efforts are put forward to understand the full spectrum of genetic variation.

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.

Cytological studies of recombination in rhesus males

An immunofluorescence approach was used to directly examine meiotic recombination events in 483 pachytene spermatocytes from 11 male rhesus monkeys. Specifically, we examined the nuclear localization patterns of the DNA mismatch repair protein MLH1, known from analyses of other mammalian species to be a useful marker of meiotic cross-overs. Our results indicated that rhesus pachytene spermatocytes contain approximately 40 cross-overs per cell, corresponding to about one cross-over per chromosome. The chromosomal distribution of these exchanges was consistent with data from human and mouse males but, surprisingly, the overall number of foci was lower, and the number of 'exchangeless' bivalents higher, than reported for either humans or mice.

Predictive models of insulin resistance derived from simple morphometric and biochemical indices related to obesity and the metabolic syndrome in baboons

Background

Non-human primates are valuable models for the study of insulin resistance and human obesity. In baboons, insulin sensitivity levels can be evaluated directly with the euglycemic clamp and is highly predicted by adiposity, metabolic markers of obesity and impaired glucose metabolism (i.e. percent body fat by DXA and HbA_1c). However, a simple method to screen and identify obese insulin resistant baboons for inclusion in interventional studies is not available.

Methods

We studied a population of twenty baboons with the euglycemic clamp technique to characterize a population of obese nondiabetic, insulin resistant baboons, and used a multivariate linear regression analysis (adjusted for gender) to test different predictive models of insulin sensitivity (insulin-stimulated glucose uptake = Rd) using abdominal circumference and fasting plasma insulin. Alternatively, we tested in a separate baboon population (n = 159), a simpler model based on body weight and fasting plasma glucose to predict the whole-body insulin sensitivity (Rd/SSPI) derived from the clamp.

Results

In the first model, abdominal circumference explained 59% of total insulin mediated glucose uptake (Rd). A second model, which included fasting plasma insulin (log transformed) and abdominal circumference, explained 64% of Rd. Finally, the model using body weight and fasting plasma glucose explained 51% of Rd/SSPI. Interestingly, we found that percent body fat was directly correlated with the adipocyte insulin resistance index (r = 0.755, p < 0.0001).

Conclusion

In baboons, simple morphometric measurements of adiposity/obesity, (i.e. abdominal circumference), plus baseline markers of glucose/lipid metabolism, (i.e. fasting plasma glucose and insulin) provide a feasible method to screen and identify overweight/obese insulin resistant baboons for inclusion in interventional studies aimed to study human obesity, insulin resistance and type 2 diabetes mellitus.

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.

A second generation genetic map for rainbow trout (Oncorhynchus mykiss)

Background

Genetic maps characterizing the inheritance patterns of traits and markers have been developed for a wide range of species and used to study questions in biomedicine, agriculture, ecology and evolutionary biology. The status of rainbow trout genetic maps has progressed significantly over the last decade due to interest in this species in aquaculture and sport fisheries, and as a model research organism for studies related to carcinogenesis, toxicology, comparative immunology, disease ecology, physiology and nutrition. We constructed a second generation genetic map for rainbow trout using microsatellite markers to facilitate the identification of quantitative trait loci for traits affecting aquaculture production efficiency and the extraction of comparative information from the genome sequences of model fish species.

Results

A genetic map ordering 1124 microsatellite loci spanning a sex-averaged distance of 2927.10 cM (Kosambi) and having 2.6 cM resolution was constructed by genotyping 10 parents and 150 offspring from the National Center for Cool and Cold Water Aquaculture (NCCCWA) reference family mapping panel. Microsatellite markers, representing pairs of loci resulting from an evolutionarily recent whole genome duplication event, identified 180 duplicated regions within the rainbow trout genome. Microsatellites associated with genes through expressed sequence tags or bacterial artificial chromosomes produced comparative assignments with tetraodon, zebrafish, fugu, and medaka resulting in assignments of homology for 199 loci.

Conclusion

The second generation NCCCWA genetic map provides an increased microsatellite marker density and quantifies differences in recombination rate between the sexes in outbred populations. It has the potential to integrate with cytogenetic and other physical maps, identifying paralogous regions of the rainbow trout genome arising from the evolutionarily recent genome duplication event, and anchoring a comparative map with the zebrafish, medaka, tetraodon, and fugu genomes. This resource will facilitate the identification of genes affecting traits of interest through fine mapping and positional cloning of candidate genes.

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.

Characterization of ghrelin in pedigreed baboons: evidence for heritability and pleiotropy

BACKGROUND

Ghrelin is an orexigenic hormone that is produced primarily in the stomach, and stimulates food intake via its receptors situated in the hypothalamus.

OBJECTIVE

The purpose of this study was to characterize baboon ghrelin cDNA and investigate the genetic influence on the variation in plasma ghrelin levels in baboons.

METHODS AND PROCEDURES

The sample consisted of 376 baboons (263 females, 113 males) from a pedigreed colony at the Southwest Foundation for Biomedical Research, San Antonio, Texas. Ghrelin cDNA was cloned by reverse-transcription polymerase chain reaction (RT-PCR) and sequenced. Real-time RT-PCR was performed to quantify mRNA from the collected tissues. Genetic contribution to plasma ghrelin was estimated using a variance components method implemented in SOLAR.

RESULTS

The baboon coding region and predicted amino acid sequence for ghrelin showed 97 and 96% sequence identity with humans, respectively. Maximum expression of ghrelin mRNA was detected in hypothalamus and stomach. Mean +/- s.e. plasma levels of ghrelin were 3,406 +/- 99 pg/ml. A significant heritability was observed for plasma ghrelin (h(2)= 0.25, P < 0.001). A genome-wide scan revealed the evidence of suggestive linkage for a locus affecting plasma ghrelin on chromosome 9q22 (between markers D9S910 and D9S261, logarithm of the odds (LOD) score = 2.3). Significant genetic correlations (P < 0.001) among ghrelin, body weight, and leptin were observed.

DISCUSSION

These results indicate a significant genetic component in the variation of plasma ghrelin in baboons and reveal a high degree of similarity between baboon and human ghrelin with respect to its cDNA and its correlation with other obesity traits.

Physiological and molecular determinants of insulin action in the baboon

OBJECTIVE

To quantitate insulin sensitivity in lean and obese nondiabetic baboons and examine the underlying cellular/molecular mechanisms responsible for impaired insulin action to characterize a baboon model of insulin resistance.

RESEARCH DESIGN AND METHODS

Twenty baboons received a hyperinsulinemic-euglycemic clamp with skeletal muscle and visceral adipose tissue biopsies at baseline and at 30 and 120 min after insulin. Genes and protein expression of key molecules involved in the insulin signaling cascade (insulin receptor, insulin receptor substrate-1, p85, phosphatidylinositol 3-kinase, Akt, and AS160) were sequenced, and insulin-mediated changes were analyzed.

RESULTS

Overall, baboons show a wide range of insulin sensitivity (6.2 +/- 4.8 mg x kg(-1) x min(-1)), and there is a strong inverse correlation between indexes of adiposity and insulin sensitivity (r = -0.946, P < 0.001 for % body fat; r = -0.72, P < 0.001 for waist circumference). The genes and protein sequences analyzed were found to have approximately 98% identity to those of man. Insulin-mediated changes in key signaling molecules were impaired both in muscle and adipose tissue in obese insulin-resistant compared with lean insulin-sensitive baboons.

CONCLUSIONS

The obese baboon is a pertinent nonhuman primate model to examine the underlying cellular/molecular mechanisms responsible for insulin resistance and eventual development of type 2 diabetes.

Genotype-by-diet effects on co-variation in Lp-PLA2 activity and LDL-cholesterol concentration in baboons fed an atherogenic diet

Both lipoprotein-associated phospholipase A(2) (Lp-PLA(2)) activity, a biomarker of inflammation, and concentration of its primary associated lipoprotein, LDL, are correlated with adverse coronary outcomes. We previously reported a quantitative trait locus (QTL) corresponding to HSA2p24.3-p23.2 with pleiotropic effects on Lp-PLA(2) activity and LDL-cholesterol (LDL-C) concentration in baboons fed a basal diet. Here, our goal was to locate pleiotropic QTLs influencing both traits in the same baboons fed a high-cholesterol, high-fat (HCHF) diet, and to assess whether shared genetic effects on these traits differ between diets. We assayed Lp-PLA(2) activity and LDL-C concentration in 683 baboons fed the HCHF diet. We used a bivariate maximum likelihood-based variance components approach in whole-genome linkage screens to locate a QTL [logarithm of odds (LOD) = 3.13, genome-wide P = 0.019] corresponding to HSA19q12-q13.2 with pleiotropic effects on Lp-PLA(2) activity and LDL-C levels in the HCHF diet. We additionally found significant evidence of genetic variance in response to diet for Lp-PLA(2) activity (P = 0.0017) and for LDL-C concentration (P = 0.00001), revealing a contribution of genotype-by-diet interaction to covariation in these two traits. We conclude that the pleiotropic QTLs detected at 2p24.3-p23.2 and 19q12-q13.2 on the basal and HCHF diets, respectively, exert diet-specific effects on covariation in Lp-PLA(2) activity and LDL-C concentration.

Adipokine expression is associated with adipocyte volume in baboons

Baboons show significant variation in body weight and composition, coupled with insulin resistance and phenotypes associated with the metabolic syndrome. An omental adipose tissue biopsy and a fasting blood sample were collected from 40 unrelated adult baboons from the colony at Southwest Foundation for Biomedical Research in San Antonio, TX. Serum was separated for analyses of circulating levels of glucose, insulin, adiponectin, resistin, interleukin 6 (IL-6) and monocyte chemoattractant protein-1 (MCP-1 or CCL-2). Adipose tissue biopsies were analyzed for cell volume and number. Total RNA was isolated from adipose tissue and adiponectin, resistin, delta-resistin, MCP-1 and IL-6 mRNA abundance were measured using real time, quantitative RT-PCR. Partial correlation coefficients were calculated among adipokine expression, fat tissue cell volume, and circulating levels of proteins. Cell volume was significantly correlated with expression of MCP-1 (r=0.44, p<0.05) and IL-6 mRNA (r=0.47, p<0.01). A step wise regression analysis was conducted with adipose tissue cell volume as dependent variable. The model identified IL-6 mRNA levels in adipose tissue as the only predictor. These observations support the role of IL-6 as a possible paracrine regulator in adipose tissue.

Genetics of variation in adiponectin in pedigreed baboons: evidence for pleiotropic effects on adipocyte volume and serum adiponectin

To detect and localize the effects of genes influencing variation in adiponectin mRNA and protein levels, we conducted statistical genetic analyses of circulating concentrations of adiponectin and adiponectin (ADIPOQ) mRNA expression in omental adipose tissue in adult, pedigreed baboons (Papio anubis). An omental adipose tissue biopsy and blood sample were collected from 427 baboons from the colony at the Southwest Foundation for Biomedical Research, San Antonio, TX. Total RNA was isolated from adipose tissue and adiponectin mRNA levels were assayed by real-time, quantitative reverse transcriptase-PCR. Adiponectin, insulin, glucose, cholesterol, high-density lipoproteins and triglycerides were measured in fasting serum. Quantitative genetic analyses were conducted for adiponectin mRNA and serum protein using a maximum likelihood-based variance decomposition approach. A genome-wide linkage analysis was conducted using adiponectin mRNA and protein levels as phenotypes. Significant heritability was estimated for ADIPOQ mRNA levels (h2=0.19+/-0.07, P=0.01) and protein levels (h2=0.28+/-0.14, P=0.003). Genetic correlations were found between adiponectin protein and body weight (rho(G)=-0.51, P=0.03), cell volume (rho(G)=-0.73, P=0.04), serum triglycerides (rho(G)=-0.67, P=0.03), and between adiponectin mRNA and glucose (rho(G)=0.93, P<0.01). A logarithm of odds score of 2.9 was found for ADIPOQ mRNA levels on baboon chromosome 4p, which is orthologous to human 6p21. There is a significant genetic component affecting variation in the analyzed traits, and common genes may be influencing adiponectin expression, adipocyte volume, body weight and circulating triglycerides. The region on 6p21 has been linked to diabetes-related phenotypes in human studies.

A pleiotropic QTL on 2p influences serum Lp-PLA2 activity and LDL cholesterol concentration in a baboon model for the genetics of atherosclerosis risk factors

Lipoprotein-associated phospholipase A(2) (Lp-PLA(2)), the major portion of which is bound to low-density lipoprotein, is an independent biomarker of cardiovascular disease risk. To search for common genetic determinants of variation in both Lp-PLA(2) activity and LDL cholesterol (LDL-C) concentration, we assayed these substances in serum from 679 pedigreed baboons. Using a maximum likelihood-based variance components approach, we detected significant evidence for a QTL affecting Lp-PLA(2) activity (LOD=2.79, genome-wide P=0.039) and suggestive evidence for a QTL affecting LDL-C levels (LOD=2.16) at the same location on the baboon ortholog of human chromosome 2p. Because we also found a significant genetic correlation between the two traits (rho(G)=0.50, P<0.00001), we conducted bivariate linkage analyses of Lp-PLA(2) activity and LDL-C concentration. These bivariate analyses improved the evidence (LOD=3.19, genome-wide P=0.015) for a QTL at the same location on 2p, corresponding to the human cytogenetic region 2p24.3-p23.2. The QTL-specific correlation between the traits (rho(Q)=0.62) was significantly different from both zero and 1 (P[rho(Q)=0]=0.047; P[rho(Q)=1]=0.022), rejecting the hypothesis of co-incident linkage and consistent with incomplete pleiotropy at this locus. We conclude that polymorphisms at the QTL described in this study exert some genetic effects that are shared between Lp-PLA(2) activity and LDL-C concentration.

A quantitative trait locus for variation in dopamine metabolism mapped in a primate model using reference sequences from related species

Non-human primates (NHP) provide crucial research models. Their strong similarities to humans make them particularly valuable for understanding complex behavioral traits and brain structure and function. We report here the genetic mapping of an NHP nervous system biologic trait, the cerebrospinal fluid (CSF) concentration of the dopamine metabolite homovanillic acid (HVA), in an extended inbred vervet monkey (Chlorocebus aethiops sabaeus) pedigree. CSF HVA is an index of CNS dopamine activity, which is hypothesized to contribute substantially to behavioral variations in NHP and humans. For quantitative trait locus (QTL) mapping, we carried out a two-stage procedure. We first scanned the genome using a first-generation genetic map of short tandem repeat markers. Subsequently, using >100 SNPs within the most promising region identified by the genome scan, we mapped a QTL for CSF HVA at a genome-wide level of significance (peak logarithm of odds score >4) to a narrow well delineated interval (<10 Mb). The SNP discovery exploited conserved segments between human and rhesus macaque reference genome sequences. Our findings demonstrate the potential of using existing primate reference genome sequences for designing high-resolution genetic analyses applicable across a wide range of NHP species, including the many for which full genome sequences are not yet available. Leveraging genomic information from sequenced to nonsequenced species should enable the utilization of the full range of NHP diversity in behavior and disease susceptibility to determine the genetic basis of specific biological and behavioral traits.

Identification of promoter variants in baboon endothelial lipase that regulate high-density lipoprotein cholesterol levels

BACKGROUND

High-density lipoprotein cholesterol (HDL) levels are a major risk factor for cardiovascular disease. Previously we identified a quantitative trait locus on baboon chromosome 18 that regulates HDL. From positional cloning studies and expression studies, we identified the endothelial lipase gene (LIPG) as the primary candidate gene for the quantitative trait locus. The mechanism by which LIPG variation influences HDL levels has not been determined.

METHODS AND RESULTS

We identified 164 LIPG polymorphisms in a panel of sibling baboons discordant for HDL1 and genotyped putative regulatory polymorphisms in a population of 951 pedigreed baboons. With the use of quantitative trait nucleotide analysis we identified 3 polymorphisms in the LIPG promoter associated with variation in serum HDL1 levels. In addition, we demonstrated that these 3 polymorphisms affect LIPG promoter activity in vitro. In silico analysis was used to identify putative transcription factors that differentially bind the functional promoter polymorphisms.

CONCLUSIONS

These results reveal LIPG variants that specifically contribute to HDL1 levels and demonstrate mechanisms by which these polymorphisms may regulate LIPG promoter activity. Results from the present study provide a mechanism, namely variation in LIPG promoter activity possibly caused by altered transcription factor binding, by which LIPG variation affects HDL levels.

A genetic linkage map of the vervet monkey (Chlorocebus aethiops sabaeus)

The spectacular progress in genomics increasingly highlights the importance of comparative biology in biomedical research. In particular, nonhuman primates, as model systems, provide a crucial intermediate between humans and mice. The close similarities between humans and other primates are stimulating primate studies in virtually every area of biomedical research, including development, anatomy, physiology, immunology, and behavior. The vervet monkey (Chlorocebus aethiops sabaeus) is an important model for studying human diseases and complex traits, especially behavior. We have developed a vervet genetic linkage map to enable mapping complex traits in this model organism and facilitate comparative genomic analysis between vervet and other primates. Here we report construction of an initial genetic map built with about 360 human orthologous short tandem repeats (STRs) that were genotyped in 434 members of an extended vervet pedigree. The map includes 226 markers mapped in a unique order with a resolution of 9.8 Kosambi centimorgans (cM) in the vervet monkey genome, and with a total length (including all 360 markers) of 2726 cM. At least one complex and 11 simple rearrangements in marker order distinguish vervet chromosomes from human homologs. While inversions and insertions can explain a similar number of changes in marker order between vervet and rhesus homologs, mostly inversions are observed when vervet chromosome organization is compared to that in human and chimpanzee. Our results support the notion that large inversions played a less prominent role in the evolution within the group of the Old World monkeys compared to the human and chimpanzee lineages.

Quantitative trait loci for peripheral blood cell counts: a study in baboons

Increasingly, baseline peripheral blood cell counts are implicated as risk factors for common complex diseases. While genetic influences on these hematologic parameters are firmly established, the genetic architecture of the blood counts is still poorly understood. In this article we used data from 582 healthy pedigreed baboons and variance components methods to localize quantitative trait loci (QTLs) influencing complete blood count variables. Besides performing genome-wide linkage scans for each trait individually, we conducted bivariate linkage analyses for all pairwise trait combinations to also identify pleiotropic QTLs influencing several blood counts. While significant and suggestive QTLs were localized throughout the genome (LOD range: 1.5-3.5), chromosomal regions associated with the expression of various hematologic parameters stand out. In particular, our results provide significant and consistent evidence for a QTL on the orthologous human chromosome 1p that is shared by several blood counts, mainly erythrocyte parameters. In addition, multiple suggestive evidence of linkage was detected on the orthologous human chromosomes 10 (near the q-terminus) and 19 (centromeric section). Future studies should help identify the genes responsible for these QTL and elucidate their role on baseline variation in hematologic indicators of health and disease.

Genetic correlations between mandibular molar cusp areas in baboons

Primate evolutionary studies rely significantly on dental variation given the large role that teeth play in how an organism interacts with its environment (animal and plant) and conspecifics. Variation in cusp size has been shown to vary among primate taxa, although most studies to date focused on extant and extinct hominoids. Here we test the assumed hypothesis that a significant proportion of this variation in baboons is due to the additive effects of genes. We perform quantitative genetic analyses on variation in two-dimensional (2-D) mandibular molar cusp size in a captive pedigreed breeding population of baboons (Papio hamadryas) from the Southwest National Primate Research Center. These analyses show that variation in cusp size is heritable and sexually dimorphic. Additionally, we tested for genetic correlations between cusps on the same crown, between morphological homologues along the tooth row, and between cusp area and crown buccolingual width. We find that four of the six cusp pairs on the first molar have a genetic correlation of one, save for the metaconid-hypoconid and entoconid-hypoconid, which are not statistically different from zero. The second and third molars have lower genetic correlations, although the metaconid-hypoconid correlation is similarly estimated at zero and the entoconid-protoconid correlation is estimated to be one. This cross pattern of genetic and no genetic correlation does not immediately accord with the known pattern of development and/or calcification. We propose two explanative hypotheses.

An evolutionary view of human recombination

Recombination has essential functions in mammalian meiosis, which impose several constraints on the recombination process. However, recent studies have shown that, in spite of these roles, recombination rates vary tremendously among humans, and show marked differences between humans and closely related species. These findings provide important insights into the determinants of recombination rates and raise new questions about the selective pressures that affect recombination over different genomic scales, with implications for human genetics and evolutionary biology.

QTL with pleiotropic effects on serum levels of bone-specific alkaline phosphatase and osteocalcin maps to the baboon ortholog of human chromosome 6p23-21.3

Bone ALP and OC are under partial genetic control. This study of 591 pedigreed baboons shows a QTL corresponding to human 6p23-21.3 that accounts for 25% (bone ALP) and 20% (OC) of the genetic variance. A gene affecting osteoblast activity, number, or recruitment likely resides in this area.

INTRODUCTION

Serum levels of bone alkaline phosphatase (ALP) and osteocalcin (OC) reflect osteoblast activity. Both of these measures are under partial genetic control. Genetic effects on bone ALP have not been previously localized to chromosomal regions in primates, nor has the degree to which genetic effects are shared (pleiotropic) between bone ALP and OC been studied.

MATERIALS AND METHODS

We applied variance components methods to a sample of 591 adult pedigreed baboons to detect and quantify effects of genes that influence bone ALP and that have pleiotropic effects on bone ALP and OC. A univariate linkage analysis was conducted for bone ALP. Bivariate linkage analyses were conducted in areas for which the bone ALP results presented here and a previous univariate OC linkage analysis showed evidence for linkage on the same chromosome for both bone ALP and OC.

RESULTS

A quantitative trait locus (QTL) for serum levels of bone ALP is evident on the baboon ortholog of human chromosomal region 6p (LOD 2.93). Thirty-seven percent (genetic correlation [rho(G)] = 0.61) of the genetic variance in bone ALP and OC is caused by pleiotropic effects of the same gene(s). Bivariate linkage analysis revealed a QTL in the region corresponding to human chromosome 6p23-21.3, with the strongest evidence for bivariate linkage near D6S422 (LOD = 2.97 at 22 cM from our pter-most marker). D6S422 maps to 20.4 Mb in the human genome. The QTL-specific heritability (h2) is 0.25 and 0.20 for bone ALP and OC, respectively.

CONCLUSIONS

This first formal test for shared genetic effects on two serum markers of osteoblast activity indicates that a significant pleiotropic effect on bone ALP and OC levels, and thus on bone formation, is detectible. The fact that this region corresponds to one on mouse chromosome 13 that has repeatedly yielded QTLs for BMD should encourage more intensive study of the effect of genes in this region on bone maintenance and turnover.