Characterization, isolation, and culture of spermatogonial stem cells in Macaca fascicularis.

Spermatogonial stem cells (SSCs) have great applications in both reproductive and regenerative medicine. Primates including monkeys are very similar to humans with regard to physiology and pathology. Nevertheless, little is known about the isolation, the characteristics, and the culture of primate SSCs. This study was designed to identify, isolate, and culture monkey SSCs. Immunocytochemistry was used to identify markers for monkey SSCs. Glial cell line-derived neurotrophic factor family receptor alpha-1 (GFRA1)-enriched spermatogonia were isolated from monkeys, namely Macaca fascicularis (M. fascicularis), by two-step enzymatic digestion and magnetic-activated cell sorting, and they were cultured on precoated plates in the conditioned medium. Reverse transcription-polymerase chain reaction (RT-PCR), immunocytochemistry, and RNA sequencing were used to compare phenotype and transcriptomes in GFRA1-enriched spermatogonia between 0 day and 14 days of culture, and xenotransplantation was performed to evaluate the function of GFRA1-enriched spermatogonia. SSCs shared some phenotypes with rodent and human SSCs. GFRA1-enriched spermatogonia with high purity and viability were isolated from M. fascicularis testes. The freshly isolated cells expressed numerous markers for rodent SSCs, and they were cultured for 14 days. The expression of numerous SSC markers was maintained during the cultivation of GFRA1-enriched spermatogonia. RNA sequencing reflected a 97.3% similarity in global gene profiles between 0 day and 14 days of culture. The xenotransplantation assay indicated that the GFRA1-enriched spermatogonia formed colonies and proliferated in vivo in the recipient c-KitW/W (W) mutant mice. Collectively, GFRA1-enriched spermatogonia are monkey SSCs phenotypically both in vitro and in vivo. This study suggests that monkey might provide an alternative to human SSCs for basic research and application in human diseases.

Population genetics of wild Macaca fascicularis with low-coverage shotgun sequencing of museum specimens.

OBJECTIVES: Long-tailed macaques (Macaca fascicularis) are widely distributed throughout the mainland and islands of Southeast Asia, making them a useful model for understanding the complex biogeographical history resulting from drastic changes in sea levels throughout the Pleistocene. Past studies based on mitochondrial genomes (mitogenomes) of long-tailed macaque museum specimens have traced their colonization patterns throughout the archipelago, but mitogenomes trace only the maternal history. Here, our objectives were to trace phylogeographic patterns of long-tailed macaques using low-coverage nuclear DNA (nDNA) data from museum specimens. METHODS: We performed population genetic analyses and phylogenetic reconstruction on nuclear single nucleotide polymorphisms (SNPs) from shotgun sequencing of 75 long-tailed macaque museum specimens from localities throughout Southeast Asia. RESULTS: We show that shotgun sequencing of museum specimens yields sufficient genome coverage (average ~1.7%) for reconstructing population relationships using SNP data. Contrary to expectations of divergent results between nuclear and mitochondrial genomes for a female philopatric species, phylogeographical patterns based on nuclear SNPs proved to be closely similar to those found using mitogenomes. In particular, population genetic analyses and phylogenetic reconstruction from the nDNA identify two major clades within M. fascicularis: Clade A includes all individuals from the mainland along with individuals from northern Sumatra, while Clade B consists of the remaining island-living individuals, including those from southern Sumatra. CONCLUSIONS: Overall, we demonstrate that low-coverage sequencing of nDNA from museum specimens provides enough data for examining broad phylogeographic patterns, although greater genome coverage and sequencing depth would be needed to distinguish between very closely related populations, such as those throughout the Philippines.

Comparison of cellular localisation of the Ca2+ -binding proteins calbindin, calretinin and parvalbumin in the retina of four different Macaca species.

Ca2+ -binding proteins are differentially expressed in the nervous system; their functional role often remains unclear. This immunohistochemical study aimed at characterising and comparing the expression pattern of the Ca2+ -binding proteins calbindin (Calb), calretinin (Calr) and parvalbumin (Parv) in the retina of four species of macaque monkeys: Macaca fascicularis (cynomolgus macaque), M. mulatta (rhesus macaque), M. thibetana (Tibetan macaque) and M. fuscata (Japanese macaque). Calb was found in cone photoreceptors and in a subset of bipolar cells. Calr was expressed in a subpopulation of amacrine cells. Parv was present in horizontal and ganglion cells. In addition, Müller cells were stained using antibodies against the specific marker cellular retinaldehyde-binding protein (CRALBP). Immunostainings were used for calculation of the density of different cell populations. The expression pattern was similar between the examined species and between retinal regions.

Whole Mitochondrial Genomic and Y-Chromosomal Phylogenies of Burmese Long-Tailed Macaque (Macaca fascicularis aurea) Suggest Ancient Hybridization between fascicularis and sinica Species Groups.

Macaca fascicularis aurea (Burmese long-tailed macaque) is 1 of the 10 subspecies of Macaca fascicularis. Despite having few morphological differences from other subspecies, a recent phylogeographic study showed that M. f. aurea is clearly distinct genetically from Macaca fascicularis fascicularis (common long-tailed macaque) and suggests that M. f. aurea experienced a disparate evolutionary pathway versus other subspecies. To construct a detailed evolutionary history of M. f. aurea and its relationships with other macaque species, we performed phylogenetic analyses and divergence time estimation of whole mitochondrial genomes (2 M. f. aurea, 8 M. f. fascicularis, and 16 animals of 12 macaque species) and 2871 bp of the Y chromosome (1 M. f. aurea, 2 M. f. fascicularis, and 5 animals of 5 macaque species) and haplotype network analysis of 758 bp of the Y chromosome (1 M. f. aurea, 2 M. f. fascicularis, and 21 animals of 19 macaque species). Whereas the Y chromosome of M. f. aurea clustered with those of the fascicularis species group in the phylogenetic and haplotype network analyses, its mtDNA clustered within the clade of the sinica species group. Based on this phylogenetic incongruence and the estimated divergence times, we propose that proto-M. f. aurea underwent hybridization with a population of the sinica species group between 2.5 and 0.95 MYA after divergence from the common ancestor of M. fascicularis. Hybridization and introgression might have been central in the evolution of M. f. aurea, similar to what occurred in the evolution of other macaque species and subspecies.

Investigating biogeographic boundaries of the Sunda shelf: A phylogenetic analysis of two island populations of Macaca fascicularis.

OBJECTIVES: Cyclical submergence and re-emergence of the Sunda Shelf throughout the Pleistocene served as a dynamic biogeographic landscape, across which long-tailed macaques (Macaca fascicularis) have migrated and evolved. Here, we tested the integrity of the previously reported continental-insular haplotype divide reported among Y and mitochondrial DNA lineages across multiple studies. MATERIALS AND METHODS: The continental-insular haplotype divide was tested by heavily sampling wild macaques from two important biogeographic regions within Sundaland: (1) Singapore, the southernmost tip of continental Asia and (2) Bali, Indonesia, the southeastern edge of the Indonesian archipelago, immediately west of Wallace's line. Y DNA was haplotyped for samples from Bali, deep within the Indonesian archipelago. Mitochondrial D-loop from both islands was analyzed against existing data using Maximum Likelihood and Bayesian approaches. RESULTS: We uncovered both "continental" and "insular" Y DNA haplotypes in Bali. Between Singapore and Bali we found 52 unique mitochondrial haplotypes, none of which had been previously described. Phylogenetic analyses confirmed a major haplogroup division within Singapore and identified five new Singapore subclades and two primary subclades in Bali. DISCUSSION: While we confirmed the continental-insular divide among mtDNA haplotypes, maintenance of both Y DNA haplotypes on Bali, deep within the Indonesian archipelago calls into question the mechanism by which Y DNA diversity has been maintained. It also suggests the continental-insular designation is less appropriate for Y DNA, leading us to propose geographically neutral Y haplotype designations.

The whole mitochondrial genome of the Cynomolgus macaque (Macaca fascicularis).

Macaca fascicularis, known as the long-tailed macaque, is widely distributed in southern of East Asia and Southeast Asia. It was one of the most commonly used non-human primates in biomedical research. Thus, to illustrate the maternal phylogenetic status of M. fascicularis in primates based on the whole mitochondrial DNA (mtDNA) genome and determine a reference sequence for future population genetic studies by taking mtDNA as molecular marker, in this study, the high quality whole mtDNA genome of M. fascicularis was amplified and sequenced. Our data showed that the whole mtDNA genome of M. fascicularis includes 16,571 base pairs (bps). Further phylogenetic analyses of M. fascicularis were performed by incorporating the 83 available whole mtDNA genomes belonging to 77 primate species with Tupaia belangeri as out-group. Our result supported that M. fascicularis belongs to Macaca. Cercopithecinae. Cercopithecidae. Anthropoidea. Primates, which has the closest genetic affinity with Macaca mulatta. In addition, the ancestral divergence between the tarsier and other primate species was supported with evidence from the whole mtDNA genomes.

A genetic comparison of two alleged subspecies of Philippine cynomolgus macaques.

Two subspecies of cynomolgus macaques (Macaca fascicularis) are alleged to co-exist in the Philippines, M. f. philippensis in the north and M. f. fascicularis in the south. However, genetic differences between the cynomolgus macaques in the two regions have never been studied to document the propriety of their subspecies status. We genotyped samples of cynomolgus macaques from Batangas in southwestern Luzon and Zamboanga in southwestern Mindanao for 15 short tandem repeat (STR) loci and sequenced an 835 bp fragment of the mtDNA of these animals. The STR genotypes were compared with those of cynomolgus macaques from southern Sumatra, Singapore, Mauritius and Cambodia, and the mtDNA sequences of both Philippine populations were compared with those of cynomolgus macaques from southern Sumatra, Indonesia and Sarawak, Malaysia. We conducted STRUCTURE and PCA analyses based on the STRs and constructed a median joining network based on the mtDNA sequences. The Philippine population from Batangas exhibited much less genetic diversity and greater genetic divergence from all other populations, including the Philippine population from Zamboanga. Sequences from both Batangas and Zamboanga were most closely related to two different mtDNA haplotypes from Sarawak from which they are apparently derived. Those from Zamboanga were more recently derived than those from Batangas, consistent with their later arrival in the Philippines. However, clustering analyses do not support a sufficient genetic distinction of cynomolgus macaques from Batangas from other regional populations assigned to subspecies M. f. fascicularis to warrant the subspecies distinction M. f. philippensis.

Divergence and diversity of ULBP2 genes in rhesus and cynomolgus macaques.

Non-human primates such as rhesus macaque and cynomolgus macaque are important animals for medical research fields and they are classified as Old World monkey, in which genome structure is characterized by gene duplications. In the present study, we investigated polymorphisms in two genes for ULBP2 molecules that are ligands for NKG2D. A total of 15 and 11 ULBP2.1 alleles and 11 and 10 ULBP2.2 alleles were identified in rhesus macaques and cynomolgus macaques, respectively. Nucleotide sequences of exons for extra cellular domain were highly polymorphic and more than 70 % were non-synonymous variations in both ULBP2.1 and ULBP2.2. In addition, phylogenetic analyses revealed that the ULBP2.2 was diverged from a branch of ULBP2.1 along with ULBP2s of higher primates. Moreover, when 3D structural models were constructed for the rhesus ULBP2 molecules, residues at presumed contact sites with NKG2D were polymorphic in ULBP2.1 and ULBP2.2 in the rhesus macaque and cynomolgus macaque, respectively. These observations suggest that amino acid replacements at the interaction sites with NKG2D might shape a specific nature of ULBP2 molecules in the Old World monkeys.

Microsatellite variation in two subspecies of cynomolgus monkeys (Macaca fascicularis).

To estimate the genetic variability of two subspecies of cynomolgus monkeys (Macaca fascicularis fascicularis and M. f. aurea) using microsatellite markers, 26 microsatellite markers were selected from previous reports. Seventeen markers showed high polymorphism in a subset of monkeys and were used for the assessment of genetic diversity in the larger sample. The effective number of alleles, the polymorphism information content (PIC) and the expected heterozygosity of M. f. aurea monkeys were all statistically significantly higher than those of M. f. fascicularis monkeys (P < 0.05), suggesting the M. f. aurea monkeys had a higher degree of genetic variation than the M. f. fascicularis monkeys. Substantial differences in allele distribution were also detected between the two subspecies of cynomolgus monkeys. Private alleles restricted to the M. f. fascicularis or the M. f. aurea monkeys were found throughout the selected 17 loci. These private alleles may allow the discrimination of the two subspecies of cynomolgus monkeys. The selected markers could also be used to estimate the genetic variation for other subspecies of cynomolgus monkeys. Further work using additional animals obtained from native or independent sources will be important for a more complete understanding of the genetic differences between these two subgroups.

Comparative phylogenetics offer new insights into the biogeographic history of Macaca fascicularis and the origin of the Mauritian macaques.

We employ a comparative phylogenetic analysis to gain insight into the recent evolutionary history of Macaca fascicularis, the long-tailed macaque. Mitochondrial and Y-chromosomal topologies both show that, in general, the deepest intraspecific bifurcations separate Indochinese and Sundaic forms of this species. Sumatran populations, however, are an exception: they carry one Y-chromosomal lineage that clusters with continental populations, and another that clusters with insular stocks. This discovery provides insight into two events in the history of M. fascicularis. First, the presence of the 'continental' Y-lineage on Sumatra is one of the strongest lines of evidence to date for recent (Late Pleistocene) gene flow between Indochinese and Sundaic populations. Second, since Sumatra is the only region known to carry 'continental' YDNA and 'insular' mtDNA, it is considered the most likely source of the Mauritian macaques-an important biomedical research stock that appears to carry this mtDNA/YDNA combination exclusively.

Characterization of the ABO blood group genes in macaques: evidence for convergent evolution.

The ABO blood group system is known to act as a major transplantation barrier in primates. Different primate species share the presence of A and B antigens. The polymorphism of the macaque ABO blood group genes was analyzed by cloning and sequencing the exon 7 region. In the case of the rhesus macaque (Macaca mulatta) and cynomolgus monkey (Macaca fascicularis) we were able to identify ABO blood group gene segments which cluster into two lineages, namely: *A/*O1 and *B. In addition allelic variation was observed. The 2 amino acid replacements at positions 266 and 268, which are thought to be crucial for A or B transferase activity, could be confirmed for both macaque species. Comparison of primate sequences shows that A and B reactivity was generated independently from each other in the hominoids and Old World monkey lineages. Hence, the primate A and B blood group genes are subject to convergent evolution.