Genetic characterization of a captive marmoset (Callithrix jacchus) colony using genotype-by-sequencing.

The marmoset is a fundamental nonhuman primate model for the study of aging, neurobiology, and many other topics. Genetic management of captive marmoset colonies is complicated by frequent chimerism in the blood and other tissues, a lack of tools to enable cost-effective, genome-wide interrogation of variation, and historic mergers and migrations of animals between colonies. We implemented genotype-by-sequencing (GBS) of hair follicle derived DNA (a minimally chimeric DNA source) of 82 marmosets housed at the Southwest National Primate Research Center (SNPRC). Our primary goals were the genetic characterization of our marmoset population for pedigree verification and colony management and to inform the scientific community of the functional genetic makeup of this valuable resource. We used the GBS data to reconstruct the genetic legacy of recent mergers between colonies, to identify genetically related animals whose relationships were previously unknown due to incomplete pedigree information, and to show that animals in the SNPRC colony appear to exhibit low levels of inbreeding. Of the >99,000 single-nucleotide variants (SNVs) that we characterized, >9800 are located within gene regions known to harbor pathogenic variants of clinical significance in humans. Overall, we show the combination of low-resolution (sparse) genotyping using hair follicle DNA is a powerful strategy for the genetic management of captive marmoset colonies and for identifying potential SNVs for the development of biomedical research models.

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.

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.

Genetic variant rs1205 is associated with COVID-19 outcomes: The Strong Heart Study and Strong Heart Family Study.

BACKGROUND: Although COVID-19 infection has been associated with a number of clinical and environmental risk factors, host genetic variation has also been associated with the incidence and morbidity of infection. The CRP gene codes for a critical component of the innate immune system and CRP variants have been reported associated with infectious disease and vaccination outcomes. We investigated possible associations between COVID-19 outcome and a limited number of candidate gene variants including rs1205. METHODOLOGY/PRINCIPAL FINDINGS: The Strong Heart and Strong Heart Family studies have accumulated detailed genetic, cardiovascular risk and event data in geographically dispersed American Indian communities since 1988. Genotypic data and 91 COVID-19 adjudicated deaths or hospitalizations from 2/1/20 through 3/1/23 were identified among 3,780 participants in two subsets. Among 21 candidate variants including genes in the interferon response pathway, APOE, TMPRSS2, TLR3, the HLA complex and the ABO blood group, only rs1205, a 3' untranslated region variant in the CRP gene, showed nominally significant association in T-dominant model analyses (odds ratio 1.859, 95%CI 1.001-3.453, p = 0.049) after adjustment for age, sex, center, body mass index, and a history of cardiovascular disease. Within the younger subset, association with the rs1205 T-Dom genotype was stronger, both in the same adjusted logistic model and in the SOLAR analysis also adjusting for other genetic relatedness. CONCLUSION: A T-dominant genotype of rs1205 in the CRP gene is associated with COVID-19 death or hospitalization, even after adjustment for relevant clinical factors and potential participant relatedness. Additional study of other populations and genetic variants of this gene are warranted.

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 (GWAS), 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 nonhuman primate models.

Hepatic transcript signatures predict atherosclerotic lesion burden prior to a 2-year high cholesterol, high fat diet challenge.

The purpose of this study was to identify molecular mechanisms by which the liver influences total lesion burden in a nonhuman primate model (NHP) of cardiovascular disease with acute and chronic feeding of a high cholesterol, high fat (HCHF) diet. Baboons (47 females, 64 males) were fed a HCHF diet for 2 years (y); liver biopsies were collected at baseline, 7 weeks (w) and 2y, and lesions were quantified in aortic arch, descending aorta, and common iliac at 2y. Unbiased weighted gene co-expression network analysis (WGCNA) revealed several modules of hepatic genes correlated with lesions at different time points of dietary challenge. Pathway and network analyses were performed to study the roles of hepatic module genes. More significant pathways were observed in males than females. In males, we found modules enriched for genes in oxidative phosphorylation at baseline, opioid signaling at 7w, and EIF2 signaling and HNF1A and HNF4A networks at baseline and 2y. One module enriched for fatty acid ß oxidation pathway genes was found in males and females at 2y. To our knowledge, this is the first study of a large NHP cohort to identify hepatic genes that correlate with lesion burden. Correlations of baseline and 7w module genes with lesions at 2y were observed in males but not in females. Pathway analyses of baseline and 7w module genes indicate EIF2 signaling, oxidative phosphorylation, and µ-opioid signaling are possible mechanisms that predict lesion formation induced by HCHF diet consumption in males. Our findings of coordinated hepatic transcriptional response in male baboons but not female baboons indicate underlying molecular mechanisms differ between female and male primate atherosclerosis.

Diet-induced leukocyte telomere shortening in a baboon model for early stage atherosclerosis.

Reported associations between leukocyte telomere length (LTL) attrition, diet and cardiovascular disease (CVD) are inconsistent. This study explores effects of prolonged exposure to a high cholesterol high fat (HCHF) diet on LTL in a baboon model of atherosclerosis. We measured LTL by qPCR in pedigreed baboons fed a chow (n = 105) or HCHF (n = 106) diet for 2 years, tested for effects of diet on LTL, and association between CVD risk factors and atherosclerotic lesions with LTL. Though not different at baseline, after 2 years median LTL is shorter in HCHF fed baboons (P < 0.0001). Diet predicts sex- and age-adjusted LTL and LTL attrition (P = 0.0009 and 0.0156, respectively). Serum concentrations of CVD biomarkers are associated with LTL at the 2-year endpoint and LTL accounts approximately 6% of the variance in aortic lesions (P = 0.04). Although heritable at baseline (h2 = 0.27, P = 0.027) and after 2 years (h2 = 0.46, P = 0.0038), baseline LTL does not predict lesion extent after 2 years. Atherogenic diet influences LTL, and LTL is a potential biomarker for early atherosclerosis. Prolonged exposure to an atherogenic diet decreases LTL and increases LTL attrition, and shortened LTL is associated with early-stage atherosclerosis in pedigreed baboons.

Toward Precision Medicine: TBC1D4 Disruption Is Common Among the Inuit and Leads to Underdiagnosis of Type 2 Diabetes.

OBJECTIVE: A common nonsense mutation in TBC1D4 was recently found to substantially increase the odds of type 2 diabetes in Greenlandic Inuit, leading to exclusively increased postprandial glucose. We investigated the frequency and effect of the TBC1D4 mutation on glucose metabolism and type 2 diabetes diagnosis among Canadian and Alaskan Inuit. RESEARCH DESIGN AND METHODS: Exome sequencing of the TBC1D4 variant was performed in 114 Inuit from Nunavik, Canada, and Sanger sequencing was undertaken in 1,027 Alaskan Inuit from the Genetics of Coronary Artery Disease in Alaskan Natives (GOCADAN) Study. Association testing evaluated the effect of the TBC1D4 variant on diabetes-related metabolic traits and diagnosis. RESULTS: The TBC1D4 mutation was present in 27% of Canadian and Alaskan Inuit. It was strongly associated with higher glucose (effect size +3.3 mmol/L; P = 2.5 x 10-6) and insulin (effect size +175 pmol/L; P = 0.04) 2 h after an oral glucose load in homozygote carriers. TBC1D4 carriers with prediabetes and type 2 diabetes had an increased risk of remaining undiagnosed unless postprandial glucose values were tested (odds ratio 5.4 [95% CI 2.5-12]) compared with noncarriers. Of carriers with prediabetes or type 2 diabetes, 32% would remain undiagnosed without an oral glucose tolerance test (OGTT). CONCLUSIONS: Disruption of TBC1D4 is common among North American Inuit, resulting in exclusively elevated postprandial glucose. This leads to underdiagnosis of type 2 diabetes, unless an OGTT is performed. Accounting for genetic factors in the care of Inuit with diabetes provides an opportunity to implement precision medicine in this population.

An initial genetic linkage map of the rhesus macaque (Macaca mulatta) genome using human microsatellite loci.

Rhesus macaques (Macaca mulatta) are the most widely used nonhuman primate species in biomedical research. To create new opportunities for genetic and genomic studies using rhesus monkeys, we constructed a genetic linkage map of the rhesus genome. This map consists of 241 microsatellite loci, all previously mapped in the human genome. These polymorphisms were genotyped in five pedigrees of rhesus monkeys totaling 865 animals. The resulting linkage map covers 2048 cM including all 20 rhesus autosomes, with average spacing between markers of 9.3 cM. Average heterozygosity among those markers is 0.73. This linkage map provides new comparative information concerning locus order and interlocus distances in humans and rhesus monkeys. The map will facilitate whole-genome linkage screens to locate quantitative trait loci (QTLs) that influence individual variation in phenotypic traits related to basic primate anatomy, physiology, and behavior, as well as QTLs relevant to risk factors for human disease.