Search

The genus Macaca is widely distributed, occupies a variety of habitats, shows diverse phenotypic characteristics, and is one of the best-studied genera of nonhuman primates. Here, we reported five re-sequencing Macaca genomes, including one M. cyclopis, one M. fuscata, one M. thibetana, one M. silenus, and one M. sylvanus. Together with published genomes of other macaque species, we combined 20 genome sequences of 10 macaque species to investigate the gene introgression and genetic differences among the species. The network analysis of the SNV-fragment trees indicates a reticular phylogeny of macaque species. Combining the results from various analytical methods, we identified extensive ancient introgression events among macaque species. The multiple introgression signals between different species groups were also observed, such as between fascicularis group species and silenus group species. However, gene flow signals between fascicularis and sinica group were not as strong as those between fascicularis group and silenus group. On the other hand, the unidirect gene flow in M. arctoides probably occurred between the progenitor of M. arctoides and the common ancestor of fascicularis group. Our study also shows that the genetic backgrounds and genetic diversity of different macaques vary dramatically among species, even among populations of the same species. In conclusion, using whole genome sequences and multiple methods, we have studied the evolutionary history of the genus Macaca and provided evidence for extensive introgression among the species., (© 2024. The Author(s), under exclusive licence to Springer Science+Business Media, LLC, part of Springer Nature.)

Cis - and trans -regulatory mutations are important contributors to transcriptome evolution. Quantifying their relative contributions to intraspecific variation in gene expression is essential for understanding the population genetic processes that underlie evolutionary changes in gene expression. Here, we have examined this issue by quantifying genome-wide, allele-specific expression (ASE) variation using a crossing scheme that produces F 1 hybrids between 18 different Drosophila melanogaster strains sampled from the Drosophila Genetic Reference Panel and a reference strain from another population. Head and body samples from F 1 adult females were subjected to RNA sequencing and the subsequent ASE quantification. Cis - and trans -regulatory effects on expression variation were estimated from these data. A higher proportion of genes showed significant cis -regulatory variation (∼28%) than those that showed significant trans -regulatory variation (∼9%). The sizes of cis -regulatory effects on expression variation were 1.98 and 1.88 times larger than trans -regulatory effects in heads and bodies, respectively. A generalized linear model analysis revealed that both cis - and trans -regulated expression variation was strongly associated with nonsynonymous nucleotide diversity and tissue specificity. Interestingly, trans -regulated variation showed a negative correlation with local recombination rate. Also, our analysis on proximal transposable element (TE) insertions suggested that they affect transcription levels of ovary-expressed genes more pronouncedly than genes not expressed in the ovary, possibly due to defense mechanisms against TE mobility in the germline. Collectively, our detailed quantification of ASE variations from a natural population has revealed a number of new relationships between genomic factors and the effects of cis - and trans -regulatory factors on expression variation., (Copyright © 2017 by the Genetics Society of America.)

The chicken domestication process represents a typical model of artificial selection, and gives significant insight into the general understanding of the influence of artificial selection on recognizable phenotypes. Two Japanese domesticated chicken varieties, the fighting cock (Shamo) and the long-crowing chicken (Naganakidori), have been selectively bred for dramatically different phenotypes. The former has been selected exclusively for aggressiveness and the latter for long crowing with an obedient sitting posture. To understand the particular mechanism behind these genetic changes during domestication, we investigated the degree of genetic differentiation in the aforementioned chickens, focusing on dopamine receptor D2, D3, and D4 genes. We studied other ornamental chickens such as Chabo chickens as a reference for comparison. When genetic differentiation was measured by an index of nucleotide differentiation (NST) newly devised in this study, we found that the NST value of DRD4 for Shamo (0.072) was distinctively larger than those of the other genes among the three populations, suggesting that aggressiveness has been selected for in Shamo by collecting a variety of single nucleotide polymorphisms. In addition, we found that in DRD4 in Naganakidori, there is a deletion variant of one proline at the 24th residue in the repeat of nine prolines of exon 1. We thus conclude that artificial selection has operated on these different kinds of genetic variation in the DRD4 genes of Shamo and Naganakidori so strongly that the two domesticated varieties have differentiated to obtain their present opposite features in a relatively short period of time.

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., (© The Author 2014. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

PR/SET domain containing 9 (Prdm9) mediates histone modifications such as H3K4me3 and marks hotspots of meiotic recombination. In many mammalian species, the Prdm9 gene is highly polymorphic. Prdm9 polymorphism is assumed to play two critical roles in evolution: to diversify the spectrum of meiotic recombination hotspots and to cause male hybrid sterility, leading to reproductive isolation and speciation. Nevertheless, information about Prdm9 sequences in natural populations is very limited. In this study, we conducted a comprehensive population survey on Prdm9 polymorphism in the house mouse, Mus musculus. Overall M. musculus Prdm9 displays an extraordinarily high level of polymorphism, particularly in regions encoding zinc finger repeats, which recognize recombination hotspots. Prdm9 alleles specific to various M. musculus subspecies dominate in subspecies territories. Moreover, introgression into other subspecies territories was found for highly divergent Prdm9 alleles associated with t-haplotype. The results of our phylogeographical analysis suggest that the requirement for hotspot diversity depends on geographical range and time span in mouse evolution, and that Prdm9 polymorphism has not been maintained by a simple balanced selection in the population of each subspecies., (© The Author 2014. Published by Oxford University Press on behalf of Kazusa DNA Research Institute.)

The ratio of genetic diversity on X chromosomes relative to autosomes in organisms with XX/XY sex chromosomes could provide fundamental insight into the process of genome evolution. Here we report this ratio for 24 cynomolgus monkeys (Macaca fascicularis) originating in Indonesia, Malaysia, and the Philippines. The average X/A diversity ratios in these samples was 0.34 and 0.20 in the Indonesian-Malaysian and Philippine populations, respectively, considerably lower than the null expectation of 0.75. A Philippine population supposed to derive from an ancestral population by founding events showed a significantly lower ratio than the parental population, suggesting a demographic effect for the reduction. Taking sex-specific mutation rate bias and demographic effect into account, expected X/A diversity ratios generated by computer simulations roughly agreed with the observed data in the intergenic regions. In contrast, silent sites in genic regions on X chromosomes showed strong reduction in genetic diversity and the observed X/A diversity ratio in the genic regions cannot be explained by mutation rate bias and demography, indicating that natural selection also reduces the level of polymorphism near genes. Whole-genome analysis of a female cynomolgus monkey also supported the notion of stronger reduction of genetic diversity near genes on the X chromosome.

The genes involved in host defences are known to undergo rapid evolution. Therefore, it is often difficult to assign orthologs in multigene families among various vertebrate species. Chemokines are a large family of small cytokines that orchestrate cell migration in health and disease. Herein, we have surveyed the genomes of 18 representative vertebrate species for chemokine genes and identified a total of 553 genes. We have determined their orthologous relationships and classified them in accordance with the current systematic chemokine nomenclature system. Our study reveals an interesting evolutionary history that gave origin and diversification to the vertebrate chemokine superfamily., (© 2012 The Authors Genes to Cells © 2012 by the Molecular Biology Society of Japan and Wiley Publishing Asia Pty Ltd.)

The "nearly neutral" theory of molecular evolution proposes that many features of genomes arise from the interaction of three weak evolutionary forces: mutation, genetic drift, and natural selection acting at its limit of efficacy. Such forces generally have little impact on allele frequencies within populations from generation to generation but can have substantial effects on long-term evolution. The evolutionary dynamics of weakly selected mutations are highly sensitive to population size, and near neutrality was initially proposed as an adjustment to the neutral theory to account for general patterns in available protein and DNA variation data. Here, we review the motivation for the nearly neutral theory, discuss the structure of the model and its predictions, and evaluate current empirical support for interactions among weak evolutionary forces in protein evolution. Near neutrality may be a prevalent mode of evolution across a range of functional categories of mutations and taxa. However, multiple evolutionary mechanisms (including adaptive evolution, linked selection, changes in fitness-effect distributions, and weak selection) can often explain the same patterns of genome variation. Strong parameter sensitivity remains a limitation of the nearly neutral model, and we discuss concave fitness functions as a plausible underlying basis for weak selection.

Recovering the genetic divergence between species is one of the major interests in the evolutionary biology. It requires accurate estimation of the neutral substitution rates. Arabidopsis thaliana, the first whole-genome sequenced plant, and its out-crossing relatives provide an ideal model for examining the split between sister species. In the study, rates of molecular evolution at markers frequently used for systematics and population genetics, including 14 nuclear genes spanning most chromosomes, three noncoding regions of chloroplast genome, and one intron of mitochondrial genome, between A. thaliana and four relatives were estimated. No deviation from neutrality was detected in the genes examined. Based on the known divergence between A. thaliana and its sisters about 8.0-17.6 MYA, evolutionary rates of the eighteen genes were estimated. Accordingly, the ratio of rates of synonymous substitutions among mitochondrial, chloroplast and nuclear genes was calculated with an average and 95% confidence interval of 1 (0.25-1.75): 15.77 (7.48-114.09): 74.79 (36.27-534.61). Molecular evolutionary rates of nuclear genes varied, with a range of 0.383-0.856×10(-8) for synonymous substitutions per site per year and 0.036-0.081×10(-9) for nonsynonymous substitutions per site per year. Compared with orthologs in Populus, a long life-span tree, genes in Arabidopsis evolved faster in an order of magnitude at the gene level, agreeing with a generation time hypothesis. The estimated substitution rates of these genes can be used as a reference for molecular dating., (Copyright © 2012 Elsevier B.V. All rights reserved.)

Accelerated rates of mitochondrial protein evolution have been proposed to reflect Darwinian coadaptation for efficient energy production for mammalian flight and brain activity. However, several features of mammalian mtDNA (absence of recombination, small effective population size, and high mutation rate) promote genome degradation through the accumulation of weakly deleterious mutations. Here, we present evidence for "compensatory" adaptive substitutions in nuclear DNA- (nDNA) encoded mitochondrial proteins to prevent fitness decline in primate mitochondrial protein complexes. We show that high mutation rate and small effective population size, key features of primate mitochondrial genomes, can accelerate compensatory adaptive evolution in nDNA-encoded genes. We combine phylogenetic information and the 3D structure of the cytochrome c oxidase (COX) complex to test for accelerated compensatory changes among interacting sites. Physical interactions among mtDNA- and nDNA-encoded components are critical in COX evolution; amino acids in close physical proximity in the 3D structure show a strong tendency for correlated evolution among lineages. Only nuclear-encoded components of COX show evidence for positive selection and adaptive nDNA-encoded changes tend to follow mtDNA-encoded amino acid changes at nearby sites in the 3D structure. This bias in the temporal order of substitutions supports compensatory weak selection as a major factor in accelerated primate COX evolution.

Background: Elucidating the pattern of evolutionary changes in drug-metabolizing genes is an important subject not only for evolutionary but for biomedical research. We investigated the pattern of divergence and polymorphisms of macaque CYP1A1 and CYP1A2 genes, which are major drug-metabolizing genes in humans. In humans, CYP1A2 is specifically expressed in livers while CYP1A1 has a wider gene expression pattern in extrahepatic tissues. In contrast, macaque CYP1A2 is expressed at a much lower level than CYP1A1 in livers. Interestingly, a previous study has shown that Macaca fascicularis CYP1A2 harbored unusually high genetic diversity within species. Genomic regions showing high genetic diversity within species is occasionally interpreted as a result of balancing selection, where natural selection maintains highly diverged alleles with different functions. Nevertheless many other forces could create such signatures., Results: We found that the CYP1A1/2 gene copy number and orientation has been highly conserved among mammalian genomes. The signature of gene conversion between CYP1A1 and CYP1A2 was detected, but the last gene conversion event in the simian primate lineage occurred before the Catarrhini-Platyrrhini divergence. The high genetic diversity of macaque CYP1A2 therefore cannot be explained by gene conversion between CYP1A1 and CYP1A2. By surveying CYP1A2 polymorphisms in total 91 M. fascicularis and M. mulatta, we found several null alleles segregating in these species, indicating functional constraint on CYP1A2 in macaques may have weakened after the divergence between humans and macaques. We propose that the high genetic diversity in macaque CYP1A2 is partly due to the degeneration of CpG sites, which had been maintained at a high level by purifying selection, and the rapid degeneration process was initiated by the loss of functional constraint on macaque CYP1A2., Conclusions: Our findings show that the highly polymorphic CYP1A2 gene in macaques has not been created by balancing selection but by the burst of CpG site degeneration after loss of functional constraint. Because the functional importance of CYP1A1/2 genes is different between humans and macaques, we have to be cautious in extrapolating a drug-testing data using substrates metabolized by CYP1A genes from macaques to humans, despite of their somewhat overlapping substrate specificity.

Population genetics of species living in marginal habitats could be particularly informative about the genetics of adaptation, but such analyses have not been readily feasible until recently. Sonneratia alba, a mangrove species widely distributed in the Indo-West Pacific, provides a very suitable system for the study of local adaptation. In this study, we analyzed DNA variation by pooling 71 genes from 85-100 individuals for DNA sequencing. For each of the two nearby S. alba populations, we obtained ~2,500 × coverage on the Illumina GA platform and for the Sanya population, an additional 5,400 × coverage on the AB SOLiD platform. For the Sanya sample, although each sequencing method called many putative single nucleotide polymorphisms, the two sets of calls did not overlap, suggesting platform-dependent errors. Conventional sequencing corroborated that each population is monomorphic. The two populations differ by 54 bp of 79,000 sites, but 90% of the variants are found in 10% of the genes. Strong local adaptation and high migration may help to explain the extensive monomorphism shared by the two populations in the presence of a small number of highly differentiated loci.

Chemokines play an important role in orchestrating cell recruitment and localization in both physiological and pathological conditions. More than 44 ligands have been identified in the human genome. A significantly different set of chemokines, however, is found in the mouse genome, suggesting a rapid evolution of the chemokine system in mammalian genomes. Thus, there are lineage and even individual-specific differences in chemokine genes in mammals. Differences in the expression and function between even recently duplicated genes are also evident. In this review, we discuss how evolutionary events such as gene duplication and gene conversion have shaped the diverse arrays of chemokines in mammalian genomes., (Copyright © 2010 Elsevier Ltd. All rights reserved.)

Background: Cynomolgus macaques (Macaca fascicularis) are widely used as experimental animals in biomedical research and are closely related to other laboratory macaques, such as rhesus macaques (M. mulatta). We isolated 85,721 clones and determined 9407 full-insert sequences from cynomolgus monkey brain, testis, and liver. These sequences were annotated based on homology to human genes and stored in a database, QFbase http://genebank.nibio.go.jp/qfbase/., Results: We found that 1024 transcripts did not represent any public human cDNA sequence and examined their expression using M. fascicularis oligonucleotide microarrays. Significant expression was detected for 544 (51%) of the unidentified transcripts. Moreover, we identified 226 genes containing exon alterations in the untranslated regions of the macaque transcripts, despite the highly conserved structure of the coding regions. Considering the polymorphism in the common ancestor of cynomolgus and rhesus macaques and the rate of PCR errors, the divergence time between the two species was estimated to be around 0.9 million years ago., Conclusion: Transcript data from Old World monkeys provide a means not only to determine the evolutionary difference between human and non-human primates but also to unveil hidden transcripts in the human genome. Increasing the genomic resources and information of macaque monkeys will greatly contribute to the development of evolutionary biology and biomedical sciences.

It is possible to assess positive selection by using the ratio of K(a) (nonsynonymous substitutions per plausible nonsynonymous sites) to K(s) (synonymous substitutions per plausible synonymous sites). We have searched candidate genes positively selected during primate evolution by using 5'-end sequences of 21,302 clones derived from cynomolgus monkey (Macaca fascicularis) brain cDNA libraries. Among these candidates, 10 genes that had not been shown by previous studies to undergo positive selection exhibited a K(a)/K(s) ratio > 1. Of the 10 candidate genes we found, 5 were included in the mitochondrial respiratory enzyme complexes, suggesting that these nuclear-encoded genes coevolved with mitochondrial-encoded genes, which have high mutation rates. The products of other candidate genes consisted of a cell-surface protein, a member of the lipocalin family, a nuclear transcription factor, and hypothetical proteins.