A genome-wide study of DNA methylation patterns and gene expression levels in multiple human and chimpanzee tissues

The modification of DNA by methylation is an important epigenetic mechanism that affects the spatial and temporal regulation of gene expression. Methylation patterns have been described in many contexts within and across a range of species. However, the extent to which changes in methylation might underlie inter-species differences in gene regulation, in particular between humans and other primates, has not yet been studied. To this end, we studied DNA methylation patterns in livers, hearts, and kidneys from multiple humans and chimpanzees, using tissue samples for which genome-wide gene expression data were also available. Using the multi-species gene expression and methylation data for 7,723 genes, we were able to study the role of promoter DNA methylation in the evolution of gene regulation across tissues and species. We found that inter-tissue methylation patterns are often conserved between humans and chimpanzees. However, we also found a large number of gene expression differences between species that might be explained, at least in part, by corresponding differences in methylation levels. In particular, we estimate that, in the tissues we studied, inter-species differences in promoter methylation might underlie as much as 12%-18% of differences in gene expression levels between humans and chimpanzees.     

The evolutionary history of human and chimpanzee Y-chromosome gene loss

Recent studies have suggested that gene gain and loss may contribute significantly to the divergence between humans and chimpanzees. Initial comparisons of the human and chimpanzee Y-chromosomes indicate that chimpanzees have a disproportionate loss of Y-chromosome genes, which may have implications for the adaptive evolution of sex-specific as well as reproductive traits, especially because one of the genes lost in chimpanzees is critically involved in spermatogenesis in humans. Here we have characterized Y-chromosome sequences in gorilla, bonobo, and several chimpanzee subspecies for 7 chimpanzee gene-disruptive mutations. Our analyses show that 6 of these gene-disruptive mutations predate chimpanzee-bonobo divergence at approximately 1.8 MYA, which indicates significant Y-chromosome change in the chimpanzee lineage relatively early in the evolutionary divergence of humans and chimpanzees.     

MicroRNA expression and regulation in human, chimpanzee, and macaque brains

Among other factors, changes in gene expression on the human evolutionary lineage have been suggested to play an important role in the establishment of human-specific phenotypes. However, the molecular mechanisms underlying these expression changes are largely unknown. Here, we have explored the role of microRNA (miRNA) in the regulation of gene expression divergence among adult humans, chimpanzees, and rhesus macaques, in two brain regions: prefrontal cortex and cerebellum. Using a combination of high-throughput sequencing, miRNA microarrays, and Q-PCR, we have shown that up to 11% of the 325 expressed miRNA diverged significantly between humans and chimpanzees and up to 31% between humans and macaques. Measuring mRNA and protein expression in human and chimpanzee brains, we found a significant inverse relationship between the miRNA and the target genes expression divergence, explaining 2%-4% of mRNA and 4%-6% of protein expression differences. Notably, miRNA showing human-specific expression localize in neurons and target genes that are involved in neural functions. Enrichment in neural functions, as well as miRNA-driven regulation on the human evolutionary lineage, was further confirmed by experimental validation of predicted miRNA targets in two neuroblastoma cell lines. Finally, we identified a signature of positive selection in the upstream region of one of the five miRNA with human-specific expression, miR-34c-5p. This suggests that miR-34c-5p expression change took place after the split of the human and the Neanderthal lineages and had adaptive significance. Taken together these results indicate that changes in miRNA expression might have contributed to evolution of human cognitive functions.     

Divergence of Gene Body DNA Methylation and Evolution of Plant Duplicate Genes

It has been shown that gene body DNA methylation is associated with gene expression. However, whether and how deviation of gene body DNA methylation between duplicate genes can influence their divergence remains largely unexplored. Here, we aim to elucidate the potential role of gene body DNA methylation in the fate of duplicate genes. We identified paralogous gene pairs from Arabidopsis and rice (Oryza sativa ssp. japonica) genomes and reprocessed their single-base resolution methylome data. We show that methylation in paralogous genes nonlinearly correlates with several gene properties including exon number/gene length, expression level and mutation rate. Further, we demonstrated that divergence of methylation level and pattern in paralogs indeed positively correlate with their sequence and expression divergences. This result held even after controlling for other confounding factors known to influence the divergence of paralogs. We observed that methylation level divergence might be more relevant to the expression divergence of paralogs than methylation pattern divergence. Finally, we explored the mechanisms that might give rise to the divergence of gene body methylation in paralogs. We found that exonic methylation divergence more closely correlates with expression divergence than intronic methylation divergence. We show that genomic environments (e.g., flanked by transposable elements and repetitive sequences) of paralogs generated by various duplication mechanisms are associated with the methylation divergence of paralogs. Overall, our results suggest that the changes in gene body DNA methylation could provide another avenue for duplicate genes to develop differential expression patterns and undergo different evolutionary fates in plant genomes.     

Complex events in the evolution of the haptoglobin gene cluster in primates

Southern blot analyses of genomic DNA show that new world monkeys have only one haptoglobin gene but that chimpanzees, gorillas, orangutans, and old world monkeys have three. Humans have two: haptoglobin (Hp) and haptoglobin-related (Hpr). These observations suggest that a triplication of the haptoglobin locus occurred after the divergence of the new world monkeys, followed by a deletion of one locus in humans. To investigate these events, we have cloned the haptoglobin gene cluster in chimpanzee. The organization of the Hp and Hpr genes in chimpanzees is the same as in humans, including a retrovirus-like sequence in the first intron of Hpr. The third gene, which we name Hpp for haptoglobin primate, is 16 kilobases downstream of Hpr. A second copy of the retrovirus-like sequence occurs between Hpr and Hpp. The nucleotide sequence of the chimpanzee Hpp gene suggests that it may code for a functional protein, but the chimpanzee Hpr gene has a single base deletion in exon 5 that causes a frameshift. Comparison of the human and chimpanzee sequences suggests that the human Hpr gene was generated by a homologous unequal crossover between ancestral Hpr and Hpp genes. The crossover point lies within a 1.3-kilobase region containing exon 5 and 500 nucleotides 3' to the genes, but the exact point is obscured by a subsequent gene conversion event.     

Human and chimpanzee gene expression differences replicated in mice fed different diets

Although the human diet is markedly different from the diets of closely related primate species, the influence of diet on phenotypic and genetic differences between humans and other primates is unknown. In this study, we analyzed gene expression in laboratory mice fed diets typical of humans and of chimpanzees. The effects of human diets were found to be significantly different from that of a chimpanzee diet in the mouse liver, but not in the brain. Importantly, 10% of the genes that differ in their expression between humans and chimpanzee livers differed also between the livers of mice fed the human and chimpanzee diets. Furthermore, both the promoter sequences and the amino acid sequences of these diet-related genes carry more differences between humans and chimpanzees than random genes. Our results suggest that the mouse can be used to study at least some aspects of human-specific traits.     

Chromosomal rearrangements and the genomic distribution of gene-expression divergence in humans and chimpanzees

The genomic DNA sequences of humans and chimpanzees differ by only 1.24%. Recently, however, substantial differences in gene-expression patterns between the two species have been revealed. In this article, we investigate the genomic distribution of such differences. Besides confirming previous findings about the evolution of sex chromosomes and duplications, we show that chromosomal rearrangements are associated with increased gene-expression differences in the brain and that rearrangements can have both direct and indirect effects on the expression of linked genes. In addition, our results are consistent with a role for some rearrangements in the original speciation events that separated the human and chimpanzee lineages.     

Reexamination of the African hominoid trichotomy with additional sequences from the primate beta-globin gene cluster.

Additional DNA sequence information from a range of primates, including 13.7 kb from pygmy chimpanzee (Pan paniscus), was added to data sets of beta-globin gene cluster sequence alignments that span the gamma 1, gamma 2, and psi eta loci and their flanking and intergenic regions. This enlarged body of data was used to address the issue of whether the ancestral separations of gorilla, chimpanzee, and human lineages resulted from only one trichotomous branching or from two dichotomous branching events. The degree of divergence, corrected for superimposed substitutions, seen in the beta-globin gene cluster between human alleles is about a third to a half that observed between two species of chimpanzee and about a fourth that between human and chimpanzee. The divergence either between chimpanzee and gorilla or between human and gorilla is slightly greater than that between human and chimpanzee, suggesting that the ancestral separations resulted from two closely spaced dichotomous branchings. Maximum parsimony analysis further strengthened the evidence that humans and chimpanzees share the longest common ancestry. Support for this human-chimpanzee clade is statistically significant at P = 0.002 over a human-gorilla clade or a chimpanzee-gorilla clade. An analysis of expected and observed homoplasy revealed that the number of sequence changes uniquely shared by human and chimpanzee lineages is too large to be attributed to homoplasy. Molecular clock calculations that accommodated lineage variations in rates of molecular evolution yielded hominoid branching times that ranged from 17-19 million years ago (MYA) for the separation of gibbon from the other hominoids to 5-7 MYA for the separation of chimpanzees from humans. Based on the relatively late dates and mounting corroborative evidence from unlinked nuclear genes and mitochondrial DNA for the close sister grouping of humans and chimpanzees, a cladistic classification would place all apes and humans in the same family. Within this family, gibbons would be placed in one subfamily and all other extant hominoids in another subfamily. The later subfamily would be divided into a tribe for orangutans and another tribe for gorillas, chimpanzees, and humans. Finally, gorillas would be placed in one subtribe with chimpanzees and humans in another, although this last division is not as strongly supported as the other divisions.     

Evidence for widespread degradation of gene control regions in hominid genomes

Although sequences containing regulatory elements located close to protein-coding genes are often only weakly conserved during evolution, comparisons of rodent genomes have implied that these sequences are subject to some selective constraints. Evolutionary conservation is particularly apparent upstream of coding sequences and in first introns, regions that are enriched for regulatory elements. By comparing the human and chimpanzee genomes, we show here that there is almost no evidence for conservation in these regions in hominids. Furthermore, we show that gene expression is diverging more rapidly in hominids than in murids per unit of neutral sequence divergence. By combining data on polymorphism levels in human noncoding DNA and the corresponding human–chimpanzee divergence, we show that the proportion of adaptive substitutions in these regions in hominids is very low. It therefore seems likely that the lack of conservation and increased rate of gene expression divergence are caused by a reduction in the effectiveness of natural selection against deleterious mutations because of the low effective population sizes of hominids. This has resulted in the accumulation of a large number of deleterious mutations in sequences containing gene control elements and hence a widespread degradation of the genome during the evolution of humans and chimpanzees.     

Contrasting effects of natural selection on human and chimpanzee CC chemokine receptor 5

Human immunodeficiency virus type 1 (HIV-1) evolved via cross-species transmission of simian immunodeficiency virus (SIVcpz) from chimpanzees (Pan troglodytes). Chimpanzees, like humans, are susceptible to infection by HIV-1. However, unlike humans, infected chimpanzees seldom develop immunodeficiency when infected with SIVcpz or HIV-1. SIVcpz and most strains of HIV-1 require the cell-surface receptor CC chemokine receptor 5 (CCR5) to infect specific leukocyte subsets, and, subsequent to infection, the level of CCR5 expression influences the amount of HIV-1 entry and the rate of HIV-1 replication. Evidence that variants in the 5' cis-regulatory region of CCR5 (5'CCR5) affect disease progression in humans suggests that variation in CCR5 might also influence the response of chimpanzees to HIV-1/SIVcpz. To determine whether patterns of genetic variation at 5'CCR5 in chimpanzees are similar to those in humans, we analyzed patterns of DNA sequence variation in 37 wild-born chimpanzees (26 P. t. verus, 9 P. t. troglodytes, and 2 P. t. schweinfurthii), along with previously published 5'CCR5 data from 112 humans and 50 noncoding regions in the human and chimpanzee genomes. These analyses revealed that patterns of variation in 5'CCR5 differ dramatically between chimpanzees and humans. In chimpanzees, 5'CCR5 was less diverse than 80% of noncoding regions and was characterized by an excess of rare variants. In humans, 5'CCR5 was more diverse than 90% of noncoding regions and had an excess of common variants. Under a wide range of demographic histories, these patterns suggest that, whereas human 5'CCR5 has been subject to balancing selection, chimpanzee 5'CCR5 has been influenced by a selective sweep. This result suggests that chimpanzee 5'CCR5 might harbor or be linked to functional variants that influence chimpanzee resistance to disease caused by SIVcpz/HIV-1.     

Genome-wide DNA methylation analyses in the brain reveal four differentially methylated regions between humans and non-human primates

ABSTRACT: BACKGROUND: The highly improved cognitive function is the most significant change in human evolutionary history. Recently, several large-scale studies reported the evolutionary roles of DNA methylation; however, the role of DNA methylation on brain evolution is largely unknown. RESULTS: To test if DNA methylation has contributed to the evolution of human brain, with the use of MeDIP-Chip and SEQUENOM MassARRAY, we conducted a genome-wide analysis to identify differentially methylated regions (DMRs) in the brain between humans and rhesus macaques. We first identified a total of 150 candidate DMRs by the MeDIP-Chip method, among which 4DMRs were confirmed by the MassARRAY analysis. All 4 DMRs are within or close to the CpG islands, and a MIR3 repeat element was identified in one DMR, but no repeat sequence was observed in the other 3 DMRs. For the 4 DMR genes, their proteins tend to be conserved and two genes have neural related functions. Bisulfite sequencing and phylogenetic comparison among human, chimpanzee, rhesus macaque and rat suggested several regions of lineage specific DNA methylation, including a human specific hypomethylated region in the promoter of K6IRS2 gene. CONCLUSIONS: Our study provides a new angle of studying human brain evolution and understanding the evolutionary role of DNA methylation in the central nervous system. The results suggest that the patterns of DNA methylation in the brain are in general similar between humans and nonhuman primates, and only a few DMRs were identified.     

A scan for positively selected genes in the genomes of humans and chimpanzees

Since the divergence of humans and chimpanzees about 5 million years ago, these species have undergone a remarkable evolution with drastic divergence in anatomy and cognitive abilities. At the molecular level, despite the small overall magnitude of DNA sequence divergence, we might expect such evolutionary changes to leave a noticeable signature throughout the genome. We here compare 13,731 annotated genes from humans to their chimpanzee orthologs to identify genes that show evidence of positive selection. Many of the genes that present a signature of positive selection tend to be involved in sensory perception or immune defenses. However, the group of genes that show the strongest evidence for positive selection also includes a surprising number of genes involved in tumor suppression and apoptosis, and of genes involved in spermatogenesis. We hypothesize that positive selection in some of these genes may be driven by genomic conflict due to apoptosis during spermatogenesis. Genes with maximal expression in the brain show little or no evidence for positive selection, while genes with maximal expression in the testis tend to be enriched with positively selected genes. Genes on the X chromosome also tend to show an elevated tendency for positive selection. We also present polymorphism data from 20 Caucasian Americans and 19 African Americans for the 50 annotated genes showing the strongest evidence for positive selection. The polymorphism analysis further supports the presence of positive selection in these genes by showing an excess of high-frequency derived nonsynonymous mutations.     

Intragene higher order repeats in neuroblastoma breakpoint family genes distinguish humans from chimpanzees

Much attention has been devoted to identifying genomic patterns underlying the evolution of the human brain and its emergent advanced cognitive capabilities, which lie at the heart of differences distinguishing humans from chimpanzees, our closest living relatives. Here, we identify two particular intragene repeat structures of noncoding human DNA, spanning as much as a hundred kilobases, that are present in human genome but are absent from the chimpanzee genome and other nonhuman primates. Using our novel computational method Global Repeat Map, we examine tandem repeat structure in human and chimpanzee chromosome 1. In human chromosome 1, we find three higher order repeats (HORs), two of them novel, not reported previously, whereas in chimpanzee chromosome 1, we find only one HOR, a 2mer alphoid HOR instead of human alphoid 11mer HOR. In human chromosome 1, we identify an HOR based on 39-bp primary repeat unit, with secondary, tertiary, and quartic repeat units, fully embedded in human hornerin gene, related to regenerating and psoriatric skin. Such an HOR is not found in chimpanzee chromosome 1. We find a remarkable human 3mer HOR organization based on the ~1.6-kb primary repeat unit, fully embedded within the neuroblastoma breakpoint family genes, which is related to the function of the human brain. Such HORs are not present in chimpanzees. In general, we find that human-chimpanzee differences are much larger for tandem repeats, in particularly for HORs, than for gene sequences. This may be of great significance in light of recent studies that are beginning to reveal the large-scale regulatory architecture of the human genome, in particular the role of noncoding sequences. We hypothesize about the possible importance of human accelerated HOR patterns as components in the gene expression multilayered regulatory network.