Mitogenomics of macaques (Macaca) across Wallace's Line in the context of modern human dispersals

Phylogeography of the Sinica Group of Macaques in the Himalayas: Taxonomic and Evolutionary Implications

Owing to the taxonomic incongruence between the morphological features and genetic relationships of the sinica group of macaques (genus Macaca), the taxonomy of this macaque group has remained inconclusive. We aimed to resolve the taxonomic quandary and improve our understanding of the historical biogeography of the group by including macaque DNA samples from previously unsampled areas in the Himalayas. We sequenced and analyzed three mitochondrial DNA loci [cytochrome b (CYTB), cytochrome oxidase subunit 1 (COI) and D-loop; 2898 bp] for sequence polymorphism, phylogenetics, species delimitation, and ancestral area reconstruction. We confirmed the occurrence of Arunachal macaque (Macaca munzala) on the southern slopes of the Eastern Himalayas in the Xizang Zizhiqu (Tibet Autonomous Region) of China. The results revealed that the sinica group of macaques is a parapatric species group composed of seven distinct species. Phylogenetic and species delimitation analyses revealed that the two previously considered subspecies of Assamese macaques (the eastern subspecies M. assamensis assamensis and the western subspecies M. a. pelops) are two distinct species. The eastern Assamese macaque is a sister species to the Tibetan macaque, whereas the western Assamese macaque and Arunachal macaque are the closest genetic sister species. The sinica group of macaques underwent five vicariance and seven dispersal radiations in the past, which mainly coincided with the Quaternary climatic oscillations between the late Pliocene and the late Pleistocene. By integrating our phylogenetic and ancestral area reconstruction results with findings from previous paleontological and molecular studies, we propose a robust hypothesis about the phylogeography of the sinica group of macaques.

Phylogeography and demographic history of macaques, fascicularis species group, in East Asia: Inferred from multiple genomic markers

Climate changes at larger scales have influenced dispersal and range shifts of many taxa in East Asia. The fascicularis species group of macaques is composed of four species and is widely distributed in Southeast and East Asia. However, its phylogeography and demographic histories are currently poorly understood. Herein, we assembled autosomal, mitogenome, and Y-chromosome data for 106 individuals, and combined them with 174 mtDNA dloop haplotypes of this species group, with particular focus on the demographic histories and dispersal routes of Macaca fuscata, M. cyclopis, and M. mulatta. The results showed: (1) three monophyletic clades for M. fuscata, M. cyclopis, and M. mulatta based on the multiple genomics analyses; (2) the disparate demographic trajectories of the three species after their split ∼1.0 Ma revealed that M. cyclopis and M. fuscata were derived from an ancestral M. mulatta population; (3) the speciation time of M. cyclopis was later than that of M. fuscata, and their divergence time occurred at the beginning of "Ryukyu Coral Sea Stage" (1.0-0.2 Ma) when the East China Sea land bridge was completely submerged by the sea level rose; and (4) the three parallel rivers (Nujiang, Lancangjiang, and Jinshajiang) of Southwestern China divided M. mulatta into Indian and Chinese genetic populations ∼200 kya. These results shed light on understanding not only the evolutionary history of the fascicularis species group but also the formation mechanism of faunal diversity in East Asia during the Pleistocene.

Global phylogeographical distribution of Gloeoporus dichrous

Phylogeographic analyses are efficient in ecological and evolutionary studies to discover the origin of a lineage, its dispersal routes, and the divergence of ancestral traits. Studies on widespread wood-decay fungi have revealed the phylogenetic division of several polypores based on geographical distribution. In this study, specimens of Gloeoporus dichrous, a cosmopolitan polypore species, were collected globally and analyzed for their geographic distribution. Multi-marker Bayesian molecular clock and haplotype analyses revealed a clear division of G. dichrous populations by continent. The species diverged from its neighboring clades 10.3 (16.0-5.6) million years ago, with Asian and North American populations at the center of divergence. Possible dispersal mechanisms and pathways are predicted and discussed based on the evaluated transfer routes. The biogeography of G. dichrous analyzed in this study represents a fraction of the polypore evolution and may advance the understanding of the overall evolution of wood-decay fungi.

Pigs as Pets: Early Human Relations with the Sulawesi Warty Pig (Sus celebensis)

The Sulawesi warty pig (S. celebensis) is a wild and still-extant suid that is endemic to the Indonesian island of Sulawesi. It has long been theorised that S. celebensis was domesticated and/or deliberately introduced to other islands in Indonesia prior to the advent of the Neolithic farming transition in the region. Thus far, however, there has been no empirical support for this idea, nor have scientists critiqued the argument that S. celebensis was a pre-Neolithic domesticate in detail. Here, it is proposed that early foragers could have formed a relationship with S. celebensis that was similar in essence to the close association between Late Pleistocene foragers in Eurasia and the wild wolf ancestors of domestic dogs. That is, a longstanding practice of hunter-gatherers intensively socialising wild-caught S. celebensis piglets for adoption into human society as companion animals ('pets') may have altered the predator-prey dynamic, brought aspects of wild pig behaviour and reproduction under indirect human selection and control, and caused changes that differentiated human-associated pigs from their solely wild-living counterparts.

Wallacean and Melanesian Islands Promote Higher Rates of Diversification within the Global Passerine radiation Corvides

The complex island archipelagoes of Wallacea and Melanesia have provided empirical data behind integral theories in evolutionary biology, including allopatric speciation and island biogeography. Yet, questions regarding the relative impact of the layered biogeographic barriers, such as deep-water trenches and isolated island systems, on faunal diversification remain underexplored. One such barrier is Wallace's Line, a significant biogeographic boundary that largely separates Australian and Asian biodiversity. To assess the relative roles of biogeographic barriers-specifically isolated island systems and Wallace's Line-we investigated the tempo and mode of diversification in a diverse avian radiation, Corvides (Crows and Jays, Birds-of-paradise, Vangas, and allies). We combined a genus-level dataset of thousands of ultraconserved elements (UCEs) and a species-level, 12-gene Sanger sequence matrix to produce a well-resolved supermatrix tree that we leveraged to explore the group's historical biogeography and effects of biogeographic barriers on their macroevolutionary dynamics. The tree is well-resolved and differs substantially from what has been used extensively for past comparative analyses within this group. We confirmed that Corvides, and its major constituent clades, arose in Australia and that a burst of dispersals west across Wallace's Line occurred after the uplift of Wallacea during the mid-Miocene. We found that dispersal across this biogeographic barrier were generally rare, though westward dispersals were two times more frequent than eastward dispersals. Wallacea's central position between Sundaland and Sahul no doubt acted as a bridge for island-hopping dispersal out of Australia, across Wallace's Line, to colonize the rest of Earth. In addition, we found that the complex island archipelagoes east of Wallace's Line harbor the highest rates of net diversification and are a substantial source of colonists to continental systems on both sides of this biogeographic barrier. Our results support emerging evidence that island systems, particularly the geologically complex archipelagoes of the Indo-pacific, are drivers of species diversification.

Functional divergence of the pigmentation gene melanocortin-1 receptor (MC1R) in six endemic Macaca species on Sulawesi Island

Coat color is often highly variable within and between animal taxa. Among hundreds of pigmentation-related genes, melanocortin-1 receptor (MC1R) plays key roles in regulating the synthesis of the dark eumelanin and the red-yellow pheomelanin. The six species of macaques that inhabit Sulawesi Island diverged rapidly from their common ancestor, M. nemestrina. Unlike most macaques, Sulawesi macaques commonly have a dark coat color, with divergence in shade and color pattern. To clarify the genetic and evolutionary basis for coat color in Sulawesi macaques, we investigated the MC1R sequences and functional properties, including basal cAMP production and α-MSH-induced activity in vitro. We found fixed non-synonymous substitutions in MC1R in each species. Furthermore, we found that six species-specific variants corresponded with variation in agonist-induced and basal activity of MC1R. Inconsistent with the dark coat color, four substitutions independently caused decreases in the basal activity of MC1R in M. hecki, M. nigra, M. tonkeana, and M. ochreata. Selective analysis suggested MC1R of M. nigra and M. nigrescens underwent purifying selection. Overall, our results suggest that fixed differences in MC1R resulted in different functional characteristics and might contribute to divergence in color among the six Sulawesi macaque species.

Evolutionary history of orangutan plasmodia revealed by phylogenetic analysis of complete mtDNA genomes and new biogeographical divergence dating calibration models

During the past 15 years, researchers have shown a renewed interest in the study of the Plasmodium parasites that infect orangutans. Most recently, studies examined the phylogenetic relationships and divergence dates of these parasites in orangutans using complete mitochondrial DNA genomes. Questions regarding the dating of these parasites, however, remain. In the present study, we provide a new calibration model for dating the origins of Plasmodium parasites in orangutans using a modified date range for the origin of macaques in Asia. Our Bayesian phylogenetic analyses of complete Plasmodium sp. mitochondrial DNA genomes inferred two clades of plasmodia in orangutans (Pongo 1 and Pongo 2), and that these clades likely represent the previously identified species Plasmodium pitheci and Plasmodium silvaticum. However, we cannot identify which Pongo clade is representative of the morphologically described species. The most recent common ancestor of both Pongo sp. plasmodia, Plasmodium. hylobati, and Plasmodium. inui dates to 3-3.16 million years ago (mya) (95% highest posterior density [HPD]: 2.09-4.08 mya). The Pongo 1 parasite diversified 0.33-0.36 mya (95% HPD: 0.12-0.63), while the Pongo 2 parasite diversified 1.15-1.22 mya (95% HPD: 0.63-1.82 mya). It now seems likely that the monkey Plasmodium (P. inui) is the result of a host switch event from the Pongo 2 parasite to sympatric monkeys, or P. hylobati. Our new estimates for the divergence of orangutan malaria parasites, and subsequent diversification, are all several hundred thousand years later than previous Bayesian estimates.

Gene flow from rhesus (Macaca mulatta) to cynomolgus macaques (M. fascicularis) and effects of introgressive hybridization on reproduction in two biomedically relevant non-human primate species

BACKGROUND

We compared the reproductive patterns of wild Indochinese and Sundaic cynomolgus macaques (Mf) exhibiting different levels of genetic admixture with rhesus macaques (Mm).

METHODS

Ten adult females from each Indochinese (WHM) and Sundaic (KN/KTK) Mf populations, which exhibited 50% and 15% of Mm autosomal SNPs, were selected as focal animals. Animals were observed for 12 months, and the frequencies of sexual proceptivity, attractivity and receptivity, number of newborns, and changes in sex skin were recorded.

RESULTS

Both populations showed all three sexual behaviors throughout the year, but they were classified as moderately seasonal breeders because their 3-month birth counts were as high as ~50%. The fecundity of WHM was lower than the KN/KTK. Changes in sex skin of WHM were more prone to Mm's pattern than the KN/KTK.

CONCLUSION

The introgressive gene flow from Mm to Mf does not affect Mf's sexual behaviors; however, it can impact fecundity and physiological (sex skin) changes.

Genomic resources for rhesus macaques (Macaca mulatta)

Rhesus macaques (Macaca mulatta) are among the most extensively studied of nonhuman primates. This species has been the subject of many investigations concerning basic primate biology and behavior, including studies of social organization, developmental psychology, physiology, endocrinology, and neurodevelopment. Rhesus macaques are also critically important as a nonhuman primate model of human health and disease, including use in studies of infectious diseases, metabolic diseases, aging, and drug or alcohol abuse. Current research addressing fundamental biological and/or applied biomedical questions benefits from various genetic and genomic analyses. As a result, the genome of rhesus macaques has been the subject of more study than most nonhuman primates. This paper briefly discusses a number of information resources that can provide interested researchers with access to genetic and genomic data describing the content of the rhesus macaque genome, available information regarding genetic variation within the species, results from studies of gene expression, and other aspects of genomic analysis. Specific online databases are discussed, including the US National Center for Biotechnology Information, the University of California Santa Cruz genome browser, Ensembl genome browser, the Macaque Genotype and Phenotype database (mGAP), Rhesusbase, and others.

Mitonuclear interactions and introgression genomics of macaque monkeys (Macaca) highlight the influence of behaviour on genome evolution

In most macaques, females are philopatric and males migrate from their natal ranges, which results in pronounced divergence of mitochondrial genomes within and among species. We therefore predicted that some nuclear genes would have to acquire compensatory mutations to preserve compatibility with diverged interaction partners from the mitochondria. We additionally expected that these sex-differences would have distinctive effects on gene flow in the X and autosomes. Using new genomic data from 29 individuals from eight species of Southeast Asian macaque, we identified evidence of natural selection associated with mitonuclear interactions, including extreme outliers of interspecies differentiation and metrics of positive selection, low intraspecies polymorphism and atypically long runs of homozygosity associated with nuclear-encoded genes that interact with mitochondria-encoded genes. In one individual with introgressed mitochondria, we detected a small but significant enrichment of autosomal introgression blocks from the source species of her mitochondria that contained genes which interact with mitochondria-encoded loci. Our analyses also demonstrate that sex-specific demography sculpts genetic exchange across multiple species boundaries. These findings show that behaviour can have profound but indirect effects on genome evolution by influencing how interacting components of different genomic compartments (mitochondria, the autosomes and the sex chromosomes) move through time and space.

A new subspecies of Trypanosoma cyclops found in the Australian terrestrial leech Chtonobdella bilineata

Previously, it was suggested that haemadipsid leeches represent an important vector of trypanosomes amongst native animals in Australia. Consequently, Chtonobdella bilineata leeches were investigated for the presence of trypanosome species by polymerase chain reaction (PCR), DNA sequencing and in vitro isolation. Phylogenetic analysis ensued to further define the populations present. PCR targeting the 28S rDNA demonstrated that over 95% of C. bilineata contained trypanosomes; diversity profiling by deep amplicon sequencing of 18S rDNA indicated the presence of four different clusters related to the Trypanosoma (Megatrypanum) theileri. Novy–MacNeal–Nicolle slopes with liquid overlay were used to isolate trypanosomes into culture that proved similar in morphology to Trypanosoma cyclops in that they contained a large numbers of acidocalcisomes. Phylogeny of 18S rDNA/GAPDH/ND5 DNA sequences from primary cultures and subclones showed the trypanosomes were monophyletic, with T. cyclops as a sister group. Blood-meal analysis of leeches showed that leeches primarily contained blood from swamp wallaby (Wallabia bicolour), human (Homo sapiens) or horse (Equus sp.). The leech C. bilineata is a host for at least five lineages of Trypanosoma sp. and these are monophyletic with T. cyclops; we propose Trypanosoma cyclops australiensis as a subspecies of T. cyclops based on genetic similarity and biogeography considerations.

Reproductive seasonality in wild northern pig-tailed macaques (Macaca leonina)

Macaque reproductive patterns range from strictly seasonal breeding to non-seasonal breeding, but factors explaining this variation are not fully understood. Valid reproductive seasonality data are also still lacking for many wild macaque populations because the majority of birth data are from captive animals living outside of their geographic range. We evaluated whether the reproductive seasonality of wild northern pig-tailed macaques falls as expected by the ecological (latitude) or phylogenetic inertia hypotheses in comparison with other macaque species. We recorded monthly occurrences of births (N = 22), copulations (N = 563), and females exhibiting sex skin swellings (N = 18) in one group at Khao Yai National Park (KYNP), Thailand. Births, copulations, and females exhibiting sex skin swellings were significantly different from a random distribution. Using measures of circular statistics and the van Schaik and colleagues' (Schaik et al. Lee (ed), Comparative primate socioecology, Cambridge University Press, Cambridge, 1999) seasonality categorizations, the population at KYNP is best characterized as moderately seasonal breeding. Despite some inconsistency, macaque reproductive seasonality was significantly influenced by latitudinal location. We broadly found that: (1) non-seasonal breeding macaque populations (birth r-vector < 0.3) lived at latitudes close to the equator between 1°S and 3°N, (2) moderately seasonal breeding macaque populations (0.3 < birth r-vector < 0.7) were found between 3°N and 14°N and at 5°S, and (3) strictly seasonal breeding macaque populations (birth r-vector > 0.7) ranged ≥ 12°N. A strong phylogenetic signal in reproductive seasonality on the macaque phylogeny was also detected. However, further studies of wild macaque populations are still needed to better characterize reproductive seasonality in this taxon.

Multilocus phylogeny suggests a distinct species status for the Nepal population of Assam macaques ( Macaca assamensis): implications for evolution and conservation

Phylogenetic relationships within the sinica-group of macaques based on morphological, behavioral, and molecular characteristics have remained controversial. The Nepal population of Assam macaques ( Macaca assamensis) (NPAM), the westernmost population of the species, is morphologically distinct but has never been used in phylogenetic analyses. Here, the phylogenetic relationship of NPAM with other congeners was tested using multiple mitochondrial and Y-chromosomal loci. The divergence times and evolutionary genetic distances among macaques were also estimated. Results revealed two major mitochondrial DNA clades of macaques under the sinica-group: the first clade included M. thibetana, M. sinica, and eastern subspecies of Assam macaque ( M. assamensis assamensis); the second clade included M. radiata together with species from the eastern and central Himalaya, namely, M. leucogenys, M. munzala, and NPAM. Among the second-clade species, NPAM was the first to diverge from the other members of the clade around 1.9 million years ago. Our results revealed that NPAM is phylogenetically distinct from the eastern Assam macaques and closer to other species and hence may represent a separate species. Because of its phylogenetic distinctiveness, isolated distribution, and small population size, the Nepal population of sinica-group macaques warrants detailed taxonomic revision and high conservation priority.