Abstract:Macaques are the most widely distributed nonhuman primates and used as animal models in biomedical research. The availability of full-genome sequences from them would be essential to both biomedical and primate evolutionary studies. Previous studies have reported whole-genome sequences from rhesus macaque (Macaca mulatta) and cynomolgus macaque (M. fascicularis, CE), both of which belong to the fascicularis group. Here, we present a 37-fold coverage genome sequence of the Tibetan macaque (M. thibetana; TM). TM is an endemic species to China belonging to the sinica group. On the basis of mapping to the rhesus macaque genome, we identified approximately 11.9 million single-nucleotide variants), of which 3.9 million were TM specific, as assessed by comparison two Chinese rhesus macaques (CR) and two CE genomes. Some genes carried TM-specific homozygous nonsynonymous variants (TSHNVs), which were scored as deleterious in human by both PolyPhen-2 and SIFT (Sorting Tolerant From Intolerant) and were enriched in the eye disease genes. In total, 273 immune response and disease-related genes carried at least one TSHNV. The heterozygosity rates of two CRs (0.002617 and 0.002612) and two CEs (0.003004 and 0.003179) were approximately three times higher than that of TM (0.000898). Polymerase chain reaction resequencing of 18 TM individuals showed that 29 TSHNVs exhibited high allele frequencies, thus confirming their low heterozygosity. Genome-wide genetic divergence analysis demonstrated that TM was more closely related to CR than to CE. We further detected unusual low divergence regions between TM and CR. In addition, after applying statistical criteria to detect putative introgression regions (PIRs) in the TM genome, up to 239,620 kb PIRs (8.84% of the genome) were identified. Given that TM and CR have overlapping geographical distributions, had the same refuge during the Middle Pleistocene, and show similar mating behaviors, it is highly likely that there was an ancient introgression event between them. Moreover, demographic inferences revealed that TM exhibited a similar demographic history as other macaques until 0.5 Ma, but then it maintained a lower effective population size until present time. Our study has provided new insight into the macaque evolutionary history, confirming hybridization events between macaque species groups based on genome-wide data.

Identification and Characterization of Polymorphic Alu Insertions in the Tibetan Macaque ( Macaca Thibetana )

Development of eighteen tetranucleotide microsatellite markers in Tibetan Macaque ( Macaca thibetana ) and genetic diversity analysis of captive population

Identification and Characterization of Short Tandem Repeats in the Tibetan Macaque Genome Based on Resequencing Data.

Isolation and characterization of twenty-one polymorphic microsatellite loci in the Tibetan macaque ( Macaca thibetana )

An Effective Microsatellite Marker Panel for Noninvasive Samples in Tibetan Macaques (macaca Thibetana)

Development and Characterization of 29 SNP Markers for the Tibetan Macaque (macaca Thibetana)

Phylogeny of the Macaques (cercopithecidae: Macaca) Based on Alu Elements

Identification of MHC-DMA and -DMB alleles in Tibetan macaques (Macaca thibetana).

Identification of 17 Novel Major Histocompatibility Complex-a Alleles in a Population of Chinese-origin Rhesus Macaques.

Identification of Major Histocompatibility Complex Class I Alleles in Chinese Rhesus Macaques

Identification of MHC‐DMA and ‐DMB Alleles in Tibetan Macaques (macaca Thibetana)

Sixteen Novel Microsatellite Loci Developed for the Giant Panda (ailuropoda Melanoleuca)

Specific Detection Method to Distinguish Tibetan Macaque from Contaminants

Simple Sequence Repeat (ssr) Polymorphisms and Population Genetics in Sichuan Wild Rhesus Macaques

Low Major Histocompatibility Complex Class II DQA Diversity in the Giant Panda (ailuropoda Melanoleuca)

Mitochondrial Genome and Nuclear Markers Provide New Insight into the Evolutionary History of Macaques.

Whole-Genome Sequencing of Tibetan Macaque (macaca Thibetana) Provides New Insight into the Macaque Evolutionary History

Major histocompatibility complex class II variation in the giant panda (Ailuropoda melanoleuca): MHC VARIATION IN THE GIANT PANDA

Characterization of the Major Histocompatibility Complex Class II DOB, DPB1, and DQB1 Alleles in Cynomolgus Macaques of Vietnamese Origin

Comprehensive identification of high-frequency and co-occurring Mafa-B, Mafa-DQB1, and Mafa-DRB alleles in cynomolgus macaques of Vietnamese origin.

Molecular Characterization of Major Histocompatibility Complex Class I Genes from the Giant Panda (ailuropoda Melanoleuca).