Search

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.)

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.

There is a mounting evidence for the correlation between the gene expression pattern and sequence divergence. However, little is known about the relationship between the gene expression pattern and polymorphism. We compiled the gene expression, polymorphism, and divergence data from the public databases of the human genome. The ratios of nonsynonymous (A) to synonymous (S) substitutions in polymorphism and divergence in the human genome were strongly influenced by the expression pattern and breadth of genes and showed strong correlations. Among the tissues we analyzed, the brain-expressed genes have the smallest and the liver-expressed genes have the largest proportion of amino acid changes both in polymorphism and divergence. The analysis implies that negative selection is the primary factor affecting expression-dependent gene evolution and the prevalent but nonuniform distribution of slightly deleterious mutations in the genome. Although the genes under relaxed negative selection evolved faster than the other genes, these genes are even more liable to slightly deleterious mutations in the population. On the other hand, nonneutral mutations in the highly conservative genes, such as brain-expressed and housekeeping genes, are largely deleterious and eliminated before they enter the population.

Nucleotide sequence polymorphisms affecting gene expression occur in the regulatory region of genes (in cis) and elsewhere in the genome (in trans). Further study is required to weigh the relative importance of cis- and trans-acting mutations in mediating gene expression differences within and between species. Here, microarray hybridization experiments were used to isolate 363 gene expression differences between the female fly head transcriptomes of 2 Drosophila melanogaster strains. One strain (French) represented the cosmopolitan M mating race and the other strain (ZS30) represented the Z mating race derived from Zimbabwe, Africa. From chromosomal substitution strains engineered from the 2 strains, we inferred that the expression differences between M and Z alleles largely could be attributed to the genotype of the chromosomes where the differentially expressed genes were located, that is, cis-regulatory polymorphisms prominently influence gene expression differences between M and Z. The effects of trans-regulatory polymorphisms were apparent yet difficult to quantify. Results have implications for models of gene regulatory evolution as well as experimental studies trying to identify the nucleotide sequence polymorphisms underlying gene expression differences between Drosophila strains.

Retrotransposable element-1 (RTE-1) is a class of long interspersed nucleotide elements that contain in its open reading frame an apurinic/apyrimidinic endonuclease domain (AP-END) and a reverse transcriptase domain. Ruminants have a clade-specific RTE-1 (BovB/RTE). The bovine bcnt gene (bucentaur or craniofacial developmental protein 1) has a duplicated paralog (bcntp97) in tandem that recruited an AP-END of BovB/RTE as a coding exon (RTE exon). We obtained sequence of the bcnt region from several animals and showed that other ruminants also have the bcntp97 with a conserved RTE exon while camels and pigs do not. Genomic Southern analysis showed that camels and pigs have multiple bcnt-related sequences but not BovB/RTE which bovines and lesser mouse deer have abundantly. These results indicate that the bcnt gene duplication followed by the creation of bcntp97 including recruitment of the RTE exon occurred in the ancestral ruminant about 55 MYA. The indication of time frame is supported by a phylogenetic analysis. Taken together with a result of differential tissue expression of the two bcnt paralogs, we conclude that bcntp97 was created concurrently with the early radiation of BovB/RTE in an ancestral ruminant and then acquired a novel function.

The substitution rate and structural divergence in the 5'-untranslated region (UTR) were investigated by using human and cynomolgus monkey cDNA sequences. Due to the weaker functional constraint in the UTR than in the coding sequence, the divergence between humans and macaques would provide a good estimate of the nucleotide substitution rate and structural divergence in the 5'UTR. We found that the substitution rate in the 5'UTR (K5UTR) averaged approximately 10%-20% lower than the synonymous substitution rate (Ks). However, both the K5UTR and nonsynonymous substitution rate (Ka) were significantly higher in the testicular cDNAs than in the brain cDNAs, whereas the Ks did not differ. Further, an in silico analysis revealed that 27% (169/622) of macaque testicular cDNAs had an altered exon-intron structure in the 5'UTR compared with the human cDNAs. The fraction of cDNAs with an exon alteration was significantly higher in the testicular cDNAs than in the brain cDNAs. We confirmed by using reverse transcriptase-polymerase chain reaction that about one-third (6/16) of in silico "macaque-specific" exons in the 5'UTR were actually macaque specific in the testis. The results imply that positive selection increased K5UTR and structural alteration rate of a certain fraction of genes as well as Ka. We found that both positive and negative selection can act on the 5'UTR sequences.

Ruminant Bcnt protein with a molecular mass of 97 kDa (designated p97Bcnt) includes a region derived from the endonuclease domain of a retrotransposable element RTE-1. Human and mouse Bcnt proteins lack the corresponding region but have a highly conserved 82-amino acid region at the C-terminus that is not present in p97Bcnt. By screening a bovine BAC library, we found two more bcnt-related genes: human-type bcnt (h-type bcnt) and its processed pseudogene. Whereas the pseudogene is localized on chromosome 26, both bcntp97 and the h-type bcnt genes are found on bovine chromosome 18, a synteny region of human chromosome 16 on which human BCNT is localized. Complete nucleotide sequencing of the BAC clone reveals that the bcntp97 and h-type bcnt genes are located just 6 kb apart in a tandem manner. The two h-type bcnt and bcntp97genes are active at both the transcriptional level and the protein level. H-type bovine Bcnt is more like human BCNT than p97Bcnt, when compared at their N-terminal regions. However, phylogenetic analysis using the N-terminal region of the bcnt gene family revealed that the duplication of bovine genes occurred within the bovine lineage with significantly accelerated substitution in bcntp97. This acceleration was not ascribed definitely to positive selection. After duplication, one of the bovine bcnt genes recruited the endonuclease domain of an intronic RTE-1 repeat accompanied by the accelerated substitution at the 5'-ORF, resulting in creation of a novel type of Bcnt protein in bovine.