Sex-linked inheritance in macaque monkeys: implications for effective population size and dispersal to Sulawesi

Taxonomy, Evolutionary and Dispersal Events of Pig-Tailed Macaque, Macaca nemestrina (Linnaeus, 1766) in Southeast Asia with Description of a New Subspecies, Macaca nemestrina perakensis in Malaysia

The pig-tailed macaque, Macaca nemestrina, which is distributed in Peninsular Malaysia, Thailand, Borneo, and Sumatra, has been the subject of unstable and changing taxonomic entity in the M. nemestrina group. This species is involved with a human-macaque conflict in Malaysia and at the same time played an important role in the ethnozoological culture of Malaysian. Even so, comprehensive phylogenetic, population genetics and biogeographical analysis of M. nemestrina in Malaysia are non-existent after decades of intensive research on the genus itself. Thus, we conducted the first comprehensive genetic study of M. nemestrina in Malaysia, based on three mitochondrial loci-Cytochrome b (567 bp), D-loop (398 bp), and COI (577 bp)-from 27 individuals representing Malaysia, plus an additional 26 sequences of Southeast Asian macaques from Genbank. Comparative biogeographical analysis in this study supports the positions of M. nemestrina in M. nemestrina groups as opposed to the silenus or Sulawesi groups. Results from this study also indicate that Bornean populations are the first extant lineages to separate from the other examined lineages of M. nemestrina, M. leonina, M. pagensis, and M. siberu in Southeast Asia. Molecular clock analysis suggested that M. nemestrina arrived in the Malay Peninsula about 0.32 million years ago (MYA). Our results indicate that the population of pig-tailed macaque from Perak (west Peninsular Malaysia) differs genetically based on all phylogenetic and population genetic analyses. Morphologically, Perak's pig-tailed macaque shows brighter coloration than M. n. nemestrina. Thus, we proposed a new subspecies for Perak's pig-tailed macaque as Macaca nemestrina perakensis distributed in the state of Perak, Peninsular Malaysia. This research helps resolve the taxonomic position and population genetics of pig-tailed macaque in Malaysia, which contribute directly to conservation and management of the species in Malaysia.

Testing Sex-Biased Admixture Origin of Macaque Species Using Autosomal and X-Chromosomal Genomic Sequences

The role of sex-specific demography in hybridization and admixture of genetically diverged species and populations is essential to understand the origins of the genomic diversity of sexually reproducing organisms. In order to infer how sex-linked loci have been differentiated undergoing frequent hybridization and admixture, we examined 17 whole-genome sequences of seven species representing the genus Macaca, which shows frequent inter-specific hybridization and predominantly female philopatry. We found that hybridization and admixture were prevalent within these species. For three cases of suggested hybrid origin of species/subspecies, Macaca arctoides, Macaca fascicularis ssp. aurea, and Chinese Macaca mulatta, we examined the level of admixture of X chromosomes, which is less affected by male-biased migration than that of autosomes. In one case, we found that Macaca cyclopis and Macaca fuscata was genetically closer to Chinese M. mulatta than to the Indian M. mulatta, and the admixture level of Chinese M. mulatta and M. fuscata/cyclopis was more pronounced on the X chromosome than on autosomes. Since the mitochondrial genomes of Chinese M. mulatta, M. cyclopis, and M. fuscata were found to cluster together, and the mitochondrial genome of Indian M. mulatta is more distantly related, the observed pattern of genetic differentiation on X-chromosomal loci is consistent with the nuclear swamping hypothesis, in which strong, continuous male-biased introgression from the ancestral Chinese M. mulatta population to a population related to M. fuscata and M. cyclopis generated incongruencies between the genealogies of the mitochondrial and nuclear genomes.

Evolutionary history of mountain voles of the subgenus Aschizomys (Cricetidae, Rodentia), inferred from mitochondrial and nuclear markers

In this study, we present an assessment of the evolutionary history and phylogenetic relationships of Asian mountain voles of the subgenus Aschizomys, genus Alticola, based on extensive sampling and phylogenetic analyses of data from mitochondrial and nuclear markers. Two species of this subgenus are widespread in the mountain areas of north-eastern Asia. However, both their distribution and taxonomic borders remained questionable for more than 100 years. Our study showed discordance in the phylogenetic patterns between nuclear and mtDNA markers. We found that mtDNA in A. lemminus is paraphyletic relative to A. macrotis, but nuclear markers demonstrated reciprocal monophyly. According to species distribution modeling, ranges of A. macrotis and A. lemminus experienced a secondary contact during the Last Glacial Maximum (approximately 22 kyr BP), and thus a hybridization event seems plausible during that period. Species tree analyses recovered a sister group relationship between the two species of the Aschizomys subgenus, with an estimated divergence date of around 0.8 Ma. Our results provided good support for currently recognized subspecies within both A. macrotis and A. lemminus based on mitochondrial and nuclear datasets. A new, yet undescribed form, supposedly of a subspecific status within A. lemminus, was found in the Bureinskiy Range in the Khabarovsk area. This finding expands the current species distribution range further to the southeast.

Speciation over the edge: gene flow among non-human primate species across a formidable biogeographic barrier

Many genera of terrestrial vertebrates diversified exclusively on one or the other side of Wallace's Line, which lies between Borneo and Sulawesi islands in Southeast Asia, and demarcates one of the sharpest biogeographic transition zones in the world. Macaque monkeys are unusual among vertebrate genera in that they are distributed on both sides of Wallace's Line, raising the question of whether dispersal across this barrier was an evolutionary one-off or a more protracted exchange-and if the latter, what were the genomic consequences. To explore the nature of speciation over the edge of this biogeographic divide, we used genomic data to test for evidence of gene flow between macaque species across Wallace's Line after macaques colonized Sulawesi. We recovered evidence of post-colonization gene flow, most prominently on the X chromosome. These results are consistent with the proposal that gene flow is a pervasive component of speciation-even when barriers to gene flow seem almost insurmountable.

Multicopy gene family evolution on primate Y chromosomes

Background

The primate Y chromosome is distinguished by a lack of inter-chromosomal recombination along most of its length, extensive gene loss, and a prevalence of repetitive elements. A group of genes on the male-specific portion of the Y chromosome known as the “ampliconic genes” are present in multiple copies that are sometimes part of palindromes, and that undergo a form of intra-chromosomal recombination called gene conversion, wherein the nucleotides of one copy are homogenized by those of another. With the aim of further understanding gene family evolution of these genes, we collected nucleotide sequence and gene copy number information for several species of papionin monkey. We then tested for evidence of gene conversion, and developed a novel statistical framework to evaluate alternative models of gene family evolution using our data combined with other information from a human, a chimpanzee, and a rhesus macaque.

Results

Our results (i) recovered evidence for several novel examples of gene conversion in papionin monkeys and indicate that (ii) ampliconic gene families evolve faster than autosomal gene families and than single-copy genes on the Y chromosome and that (iii) Y-linked singleton and autosomal gene families evolved faster in humans and chimps than they do in the other Old World Monkey lineages we studied.

Conclusions

Rapid evolution of ampliconic genes cannot be attributed solely to residence on the Y chromosome, nor to variation between primate lineages in the rate of gene family evolution. Instead other factors, such as natural selection and gene conversion, appear to play a role in driving temporal and genomic evolutionary heterogeneity in primate gene families.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2187-8) contains supplementary material, which is available to authorized users.

Complete mitochondrial genome sequence of the Eastern gorilla (Gorilla beringei) and implications for african ape biogeography

The Western and Eastern species of gorillas (Gorilla gorilla and Gorilla beringei) began diverging in the mid-Pleistocene, but in a complex pattern with ongoing gene flow following their initial split. We sequenced the complete mitochondrial genomes of 1 Eastern and 1 Western gorilla to provide the most accurate date for their mitochondrial divergence, and to analyze patterns of nucleotide substitutions. The most recent common ancestor of these genomes existed about 1.9 million years ago, slightly more recent than that of chimpanzee and bonobo. We in turn use this date as a calibration to reanalyze sequences from the Eastern lowland and mountain gorilla subspecies to estimate their mitochondrial divergence at approximately 380000 years ago. These dates help frame a hypothesis whereby populations became isolated nearly 2 million years ago with restricted maternal gene flow, followed by ongoing male migration until the recent past. This process of divergence with prolonged hybridization occurred against the backdrop of the African Pleistocene, characterized by intense fluctuations in temperature and aridity, while at the same time experiencing tectonic uplifting and consequent shifts in the drainage of major river systems. Interestingly, this same pattern of introgression following divergence and discrepancies between mitochondrial and nuclear loci is seen in fossil hominins from Eurasia, suggesting that such processes may be common in hominids and that living gorillas may provide a useful model for understanding isolation and migration in our extinct relatives.

Reduced representation genome sequencing suggests low diversity on the sex chromosomes of tonkean macaque monkeys

In species with separate sexes, social systems can differ in the relative variances of male versus female reproductive success. Papionin monkeys (macaques, mangabeys, mandrills, drills, baboons, and geladas) exhibit hallmarks of a high variance in male reproductive success, including a female-biased adult sex ratio and prominent sexual dimorphism. To explore the potential genomic consequences of such sex differences, we used a reduced representation genome sequencing approach to quantifying polymorphism at sites on autosomes and sex chromosomes of the tonkean macaque (Macaca tonkeana), a species endemic to the Indonesian island of Sulawesi. The ratio of nucleotide diversity of the X chromosome to that of the autosomes was less than the value (0.75) expected with a 1:1 sex ratio and no sex differences in the variance in reproductive success. However, the significance of this difference was dependent on which outgroup was used to standardize diversity levels. Using a new model that includes the effects of varying population size, sex differences in mutation rate between the autosomes and X chromosome, and GC-biased gene conversion (gBGC) or selection on GC content, we found that the maximum-likelihood estimate of the ratio of effective population size of the X chromosome to that of the autosomes was 0.68, which did not differ significantly from 0.75. We also found evidence for 1) a higher level of purifying selection on genic than nongenic regions, 2) gBGC or natural selection favoring increased GC content, 3) a dynamic demography characterized by population growth and contraction, 4) a higher mutation rate in males than females, and 5) a very low polymorphism level on the Y chromosome. These findings shed light on the population genomic consequences of sex differences in the variance in reproductive success, which appear to be modest in the tonkean macaque; they also suggest the occurrence of hitchhiking on the Y chromosome.

The condition dependency of fitness in males and females: the fitness consequences of juvenile diet assessed in environments differing in key adult resources

Variation in environmental or genetic quality leads to phenotypic variation in condition, but how much variation in fitness is created by this variation in condition? Using Drosophila melanogaster, we manipulated condition via alternative larval diets and then tested several key factors predicted to influence how much variation in fitness results from differences in condition. Specifically, we were interested in whether male and female fitness are affected equally by condition and whether the strength of selection on condition depends on the abundance of key resources limiting the reproductive output of each sex. We measured selection on condition in alternative assay contexts that varied in the abundance of adult food (a key resource for females) or in the abundance of females (a key resource for males). Overall, selection tended to be stronger on males than females. However, selection on males was weakened when the abundance of their key resource (females) was elevated. Increasing the abundance of the key resource for females (live yeast) elevated their reproductive output as expected but did not change the strength of selection in this sex. Instead, this manipulation increased selection on males, suggesting that this environmental factor indirectly affects selection on males via their interaction with females.

The effect of nonindependent mate pairing on the effective population size

The effective population size (N(e)) quantifies the effectiveness of genetic drift in finite populations. When generations overlap, theoretical expectations for N(e) typically assume that the sampling of offspring genotypes from a given individual is independent among successive breeding events, even though this is not true in many species, including humans. To explore the effects on N(e) of nonindependent mate pairing across breeding events, we simulated the genetic drift of mitochondrial DNA, autosomal DNA, and sex chromosome DNA under three mating systems. Nonindependent mate pairing across breeding seasons has no effect when all adults mate pair for life, a small or moderate effect when females reuse stored sperm, and a large effect when there is intense male-male competition for reproduction in a harem social system. If adult females reproduce at a constant rate irrespective of the type of mate pairing, the general effect of nonindependent mate pairing is to decrease N(e) for paternally inherited components of the genome. These findings have significant implications for the relative N(e) values of different genomic regions, and hence for the expected levels of DNA sequence diversity in these regions.

Phylogenetic inconsistency caused by ancient sex-biased gene migration

Inferences of ancient sex-biased migration patterns using sex-linked genetic markers are usually difficult because of a stochastic process of allele fixation. Nevertheless, incongruent phylogenetic trees between different sex-linked markers and between sex-linked and autosomal markers are frequently interpreted as a signature of sex-biased migration without further statistical evaluation. I investigated the types of incongruent phylogenetic trees from which past sex-biased migration events can be statistically supported under the coalescent model. In the case of mammals, detecting a sex-biased migration pattern is not guaranteed by comparing the phylogenetic pattern of mitochondrial and Y-chromosomal loci. Likewise, evidence of introgression at a mitochondrial locus, but not at autosomal loci, does not support the hypothesis of an ancient female-biased migration pattern with statistical significance. In contrast, evidence of introgression at ≥ 5 unlinked autosomal loci, but not at a Y-chromosomal locus, would reject the null hypothesis of a sexually equal migration rate with statistical significance. A similar argument can be made to infer a male-biased migration pattern. Furthermore, the investigation of many recombining sex-biased markers such as X-chromosomal loci in mammals has the potential to efficiently detect ancient sex-biased demographic patterns.

The role of the effective population size in compensatory evolution

The impact of the effective population size (N_e) on the efficacy of selection has been the focus of many theoretical and empirical studies over the recent years. Yet, the effect of N_e on evolution under epistatic fitness interactions is not well understood. In this study, we compare selective constraints at independently evolving (unpaired) and coevolving (paired) sites in orthologous transfer RNAs (tRNA molecules for vertebrate and drosophilid species pairs of different N_e. We show that patterns of nucleotide variation for the two classes of sites are explained well by Kimura's one- and two-locus models of sequence evolution under mutational pressure. We find that constraints in orthologous tRNAs increase with increasing N_e of the investigated species pair. Thereby, the effect of N_e on the efficacy of selection is stronger at unpaired sites than at paired sites. Furthermore, we identify a “core” set of tRNAs with high structural similarity to tRNAs from all major kingdoms of life and a “peripheral” set with lower similarity. We observe that tRNAs in the former set are subject to higher constraints and less prone to the effect of N_e, whereas constraints in tRNAs of the latter set show a large influence of N_e. Finally, we are able to demonstrate that constraints are relaxed in X-linked drosophilid tRNAs compared with autosomal tRNAs and suggest that N_e is responsible for this difference. The observed effects of N_e are consistent with the hypothesis that evolution of most tRNAs is governed by slightly to moderately deleterious mutations (i.e., |N_es| ≤ 5).