Genomewide comparison of DNA sequences between humans and chimpanzees

Oscillating evolution of a mammalian locus with overlapping reading frames: an XLalphas/ALEX relay

XLαs and ALEX are structurally unrelated mammalian proteins translated from alternative overlapping reading frames of a single transcript. Not only are they encoded by the same locus, but a specific XLαs/ALEX interaction is essential for G-protein signaling in neuroendocrine cells. A disruption of this interaction leads to abnormal human phenotypes, including mental retardation and growth deficiency. The region of overlap between the two reading frames evolves at a remarkable speed: the divergence between human and mouse ALEX polypeptides makes them virtually unalignable. To trace the evolution of this puzzling locus, we sequenced it in apes, Old World monkeys, and a New World monkey. We show that the overlap between the two reading frames and the physical interaction between the two proteins force the locus to evolve in an unprecedented way. Namely, to maintain two overlapping protein-coding regions the locus is forced to have high GC content, which significantly elevates its intrinsic evolutionary rate. However, the two encoded proteins cannot afford to change too quickly relative to each other as this may impair their interaction and lead to severe physiological consequences. As a result XLαs and ALEX evolve in an oscillating fashion constantly balancing the rates of amino acid replacements. This is the first example of a rapidly evolving locus encoding interacting proteins via overlapping reading frames, with a possible link to the origin of species-specific neurological differences.

One of the possible ways to achieve tight co-expression of two proteins is to encode them within a single mRNA. The GNAS1 gene in mammals does just that: it encodes two interacting signaling polypeptides within a single transcript using nested reading frames shifted one nucleotide relative to each other. The exceptionally high GC content of the region where the two reading frames overlap diminishes the probability of encountering stop codons but makes the locus highly mutable. To preserve their ability to interact functionally with each other despite the high mutation rate, the two polypeptides appear to evolve in an oscillating fashion, trying to maintain approximately equal rates of amino acid substitutions. This unexpected observation provides new insights into the evolution of mostly overlooked overlapping coding regions in eukaryotic genomes.

Trace amine-associated receptors form structurally and functionally distinct subfamilies of novel G protein-coupled receptors

Trace amines are endogenous compounds structurally related to classical biogenic amines that have been studied for decades, triggered by their link to psychiatric conditions of high epidemiological and economical relevance. The understanding of their pharmacology on the molecular level was hampered until the recent discovery of trace-amine-specific receptors. We completed the identification of all members of this novel GPCR family in human, chimpanzee, rat, and mouse and observed remarkable interspecies differences, even between human and chimpanzee. The analysis of the chromosomal localizations, phylogenetic relationships, and ligand pocket vectors reveals three distinct receptor subfamilies. As most of these receptors do not respond to trace amines, each subfamily will presumably have a distinct pharmacological profile, which remains to be identified. We propose a uniform nomenclature describing this novel GPCR family in all mammalian species as trace-amine-associated receptors (TAARs), which resolves the ambiguities and contradictions of the previous naming.

Why are young and old repetitive elements distributed differently in the human genome?

Alu elements are not distributed homogeneously throughout the human genome: old elements are preferentially found in the GC-rich parts of the genome, while young Alus are more often found in the GC-poor parts of the genome. The process giving rise to this differential distribution remains poorly understood. Here we investigate whether this pattern could be due to a preferential degradation of Alu elements integrated in GC-poor regions by small indel mutations. We aligned 5.1 Mb of human and chimpanzee sequences and examined whether the rate of insertion and deletion inside Alu elements differed according to the base composition surrounding them. We found that Alu elements are not preferentially degraded in GC-poor regions by indel events. We also looked at whether very young L1 elements show the same change in distribution compared to older ones. This analysis indicated that L1 elements also show a shift in their distribution, although we could not assess it as precisely as for Alu elements. We propose that the differential distribution of Alu elements is likely to be due to a change in their pattern of insertion or their probability of fixation through evolutionary time.

The scale of mutational variation in the murid genome

Mutation rates vary across mammalian genomes, but little is known about the scale over which this variation occurs. Knowledge of the magnitude and scale of mutational variation is required to understand the processes that drive mutation, and is essential in formulating a robust null hypothesis for comparative genomics studies. Here we estimate the scale of mutational variation in the murid genome by calculating the spatial autocorrelation of nucleotide substitution rates in ancestral repeats. Such transposable elements are good candidates for neutrally evolving sequence and therefore well suited for the study of mutation rate variation. We find that the autocorrelation coefficient decays to a value close to zero by approximately 15 Mb, with little apparent variation in mutation rate under 100 kb. We conclude that the primary scale over which mutation rates vary is subchromosomal. Furthermore, our analysis shows that within-chromosome mutational variability exceeds variation among chromosomes by approximately one order of magnitude. Thus, differences in mutation rate between different regions of the same chromosome frequently exceed differences both between whole autosomes and between autosomes and the X-chromosome. Our results indicate that factors other than the time spent in the male germ line are important in driving mutation rates. This raises questions about the biological mechanism(s) that produce new mutations and has implications for the study of male-driven evolution.

Male-driven evolution in closely related species of the mouse genus Mus

Recently, other researchers have found that closely related primate species had a lower male-to-female mutation rate ratio (alpha) than distantly related species. To determine if this is a general phenomenon affecting other mammalian orders, eleven species or subspecies of the rodent genus Mus and two outgroup species were compared. Intron sequences from a gene in the nonrecombining region of the Y chromosome Jarid1d (Smcy) and its X chromosomal gametolog, Jarid1c (Smcx), were analyzed in a phylogenetic context. The male-to-female mutation rate ratio for all thirteen taxa is approximately 2.5, which is similar to previous estimates in more distantly related rodents. However, when branches with lengths of more than 2.5% were removed from the analysis, the male-to-female mutation rate ratio dropped to 0.9. Thus, in closely related rodents, as in closely related primates, the male-to-female mutation rate ratio is lower than expected.

Multi-species microarrays reveal the effect of sequence divergence on gene expression profiles

Interspecies comparisons of gene expression levels will increase our understanding of the evolution of transcriptional mechanisms and help to identify targets of natural selection. This approach holds particular promise for apes, as many human-specific adaptations are thought to result from differences in gene expression rather than in coding sequence. To date, however, all studies directly comparing interspecies gene expression have been performed on single-species arrays, so that it has been impossible to distinguish differential hybridization due to sequence mismatches from underlying expression differences. To evaluate the severity of this potential problem, we constructed a new multiprimate cDNA array using probes from human, chimpanzee, orangutan, and rhesus. We find a large effect of sequence divergence on hybridization signal, even in the closest pair of species, human and chimpanzee. By comparing single-species array analyses with results from multispecies arrays, we examine how estimates of differential gene expression are affected by sequence divergence. Our results indicate that naive use of single-species arrays in direct interspecies comparisons can yield spurious results.

Nucleotide sequence comparison of a chromosome rearrangement on human chromosome 12 and the corresponding ape chromosomes

Chromosome rearrangement has been considered to be important in the evolutionary process. Here, we demonstrate the evolutionary relationship of the rearranged human chromosome 12 and the corresponding chromosome XII in apes (chimpanzee, bonobo, gorilla, orangutan, and gibbon) by examining PCR products derived from the breakpoints of inversions and by conducting shotgun sequencing of a gorilla fosmid clone containing the breakpoint and a "duplicated segment" (duplicon). We confirmed that a pair of 23-kb duplicons flank the breakpoints of inversions on the long and short arms of chimpanzee chromosome XII. Although only the 23-kb duplicon on the long arm of chimpanzee chromosome XII and its telomeric flanking sequence are found to be conserved among the hominoids (human, great apes, and gibbons), the duplicon on the short arm of chimpanzee chromosome XII is suggested to be the result of a duplication from that on the long arm. Furthermore, the shotgun sequencing of a gorilla fosmid indicated that the breakpoint on the long arm of the gorilla is located at a different position 1.9 kb from that of chimpanzee. The region is flanked by a sequence homologous to that of human chromosome 6q22. Our findings and sequence analysis suggest a close relationship between segmental duplication and chromosome rearrangement (or breakpoint of inversion) in Hominoidea. The role of the chromosome rearrangement in speciation is also discussed based on our new results.

Deep haplotype divergence and long-range linkage disequilibrium at xp21.1 provide evidence that humans descend from a structured ancestral population

Fossil evidence links human ancestry with populations that evolved from modern gracile morphology in Africa 130,000-160,000 years ago. Yet fossils alone do not provide clear answers to the question of whether the ancestors of all modern Homo sapiens comprised a single African population or an amalgamation of distinct archaic populations. DNA sequence data have consistently supported a single-origin model in which anatomically modern Africans expanded and completely replaced all other archaic hominin populations. Aided by a novel experimental design, we present the first genetic evidence that statistically rejects the null hypothesis that our species descends from a single, historically panmictic population. In a global sample of 42 X chromosomes, two African individuals carry a lineage of noncoding 17.5-kb sequence that has survived for >1 million years without any clear traces of ongoing recombination with other lineages at this locus. These patterns of deep haplotype divergence and long-range linkage disequilibrium are best explained by a prolonged period of ancestral population subdivision followed by relatively recent interbreeding. This inference supports human evolution models that incorporate admixture between divergent African branches of the genus Homo.

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.

Eighty percent of proteins are different between humans and chimpanzees

The chimpanzee is our closest living relative. The morphological differences between the two species are so large that there is no problem in distinguishing between them. However, the nucleotide difference between the two species is surprisingly small. The early genome comparison by DNA hybridization techniques suggested a nucleotide difference of 1-2%. Recently, direct nucleotide sequencing confirmed this estimate. These findings generated the common belief that the human is extremely close to the chimpanzee at the genetic level. However, if one looks at proteins, which are mainly responsible for phenotypic differences, the picture is quite different, and about 80% of proteins are different between the two species. Still, the number of proteins responsible for the phenotypic differences may be smaller since not all genes are directly responsible for phenotypic characters.

A comparison of the human and chimpanzee olfactory receptor gene repertoires

Olfactory receptor (OR) genes constitute the basis of the sense of smell and are encoded by the largest mammalian gene superfamily, with >1000 members. In humans, but not in mice or dogs, the majority of OR genes have become pseudogenes, suggesting that OR genes in humans evolve under different selection pressures than in other mammals. To explore this further, we compare the OR gene repertoire of human with its closest living evolutionary relative, by taking advantage of the recently sequenced genome of the chimpanzee. In agreement with previous reports based on a small number of ORs, we find that humans have a significantly higher proportion of OR pseudogenes than chimpanzees. Moreover, we can reject the possibility that humans have been accumulating OR pseudogenes at a constant neutral rate since the divergence of human and chimpanzee. The comparison of the two repertoires reveals two chimpanzee-specific OR subfamily expansions and three expansions specific to humans. It also suggests that a subset of OR genes are under positive selection in either the human or the chimpanzee lineage. Thus, although overall there is relaxed constraint on human olfaction relative to chimpanzee, species-specific sensory requirements appear to have shaped the evolution of the functional OR gene repertoires in both species.

Gene conversion and evolution of Xq28 duplicons involved in recurring inversions causing severe hemophilia A

Inversions breaking the 1041 bp int1h-1 or the 9.5-kb int22h-1 sequence of the F8 gene cause hemophilia A in 1/30,000 males. These inversions are due to homologous recombination between the above sequences and their inverted copies on the same DNA molecule, respectively, int1h-2 and int22h-2 or int22h-3. We find that (1) int1h and int22h duplicated more than 25 million years ago; (2) the identity of the copies (>99%) of these sequences in humans and other primates is due to gene conversion; (3) gene conversion is most frequent in the internal regions of int22h; (4) breakpoints of int22h-related inversions also tend to involve the internal regions of int22h; (5) sequence variations in a sample of human X chromosomes defined eight haplotypes of int22h-1 and 27 of int22h-2 plus int22h-3; (6) the latter two sequences, which lie, respectively, 500 and 600 kb telomeric to int22h-1 are five-fold more identical when in cis than when in trans, thus suggesting that gene conversion may be predominantly intrachromosomal; (7) int1h, int22h, and flanking sequences evolved at a rate of about 0.1% substitutions per million years during the divergence between humans and other primates, except for int1h during the human-chimpanzee divergence, when its rate of evolution was significantly lower. This is reminiscent of the slower evolution of palindrome arms in the male specific regions of the Y chromosome and we propose, as an explanation, that intrachromosomal gene conversion and cosegregation of the duplicated regions favors retention of the ancestral sequence and thus reduces the evolution rate.

Population genetics of CAPN10 and GPR35: implications for the evolution of type 2 diabetes variants

A positional cloning study of type 2 diabetes in Mexican Americans identified a region, termed "NIDDM1," on chromosome 2q37 with significant linkage evidence. Haplotype combinations at the calpain-10 gene (CAPN10) within this region were shown to increase diabetes risk in several populations. On the basis of the thrifty genotype hypothesis, variants that increase susceptibility to type 2 diabetes under modern lifestyle conditions provided a survival advantage in past environments by increasing the efficiency of energy use and storage. Here, our goal is to make inferences about the evolutionary forces shaping variation in genes in the NIDDM1 region and to investigate the population genetics models that may underlie the thrifty genotype hypothesis. To this end, we surveyed sequence variation in CAPN10 and in an adjacent gene, G-protein-coupled receptor 35 (GPR35), in four population samples from different ethnic groups. These data revealed two distinct deviations from the standard neutral model in CAPN10, whereas GPR35 variation was largely consistent with neutrality. CAPN10 showed a significant deficit of variation in the haplotype class defined by the derived allele at SNP44, a polymorphism that is significantly associated with diabetes in meta-analysis studies. This suggests that this haplotype class was quickly driven to high frequency by positive natural selection. Interestingly, the derived allele at SNP44 is protective against diabetes. CAPN10 also showed a local excess of polymorphism and linkage disequilibrium decay in intron 13. Simulations show that this pattern may be explained by long-standing balancing selection that maintains multiple selected alleles. Alternatively, it is possible that the local mutation and recombination rates changed since the divergence of human and chimpanzee; this scenario does not require the action of natural selection on intron 13 variation.

Further statistical analysis for genome-wide expression evolution in primate brain/liver/fibroblast tissues

In spite of only a 1-2 per cent genomic DNA sequence difference, humans and chimpanzees differ considerably in behaviour and cognition. Affymetrix microarray technology provides a novel approach to addressing a long-term debate on whether the difference between humans and chimpanzees results from the alteration of gene expressions. Here, we used several statistical methods (distance method, two-sample t-tests, regularised t-tests, ANOVA and bootstrapping) to detect the differential expression pattern between humans and great apes. Our analysis shows that the pattern we observed before is robust against various statistical methods; that is, the pronounced expression changes occurred on the human lineage after the split from chimpanzees, and that the dramatic brain expression alterations in humans may be mainly driven by a set of genes with increased expression (up-regulated) rather than decreased expression (down-regulated).

Clustering of genes coding for DNA binding proteins in a region of atypical evolution of the human genome

Comparison of the human and mouse genomes has revealed that significant variations in evolutionary rates exist among genomic regions and that a large part of this variation is interchromosomal. We confirm in this work, using a large collection of introns, that human chromosome 19 is the one that shows the highest divergence with respect to mouse. To search for other differences among chromosomes, we examine the distribution of gene functions in human and mouse chromosomes using the Gene Ontology definitions. We found by correspondence analysis that among the strongest clusterings of gene functions in human chromosomes is a group of genes coding for DNA binding proteins in chromosome 19. Interestingly, chromosome 19 also has a very high GC content, a feature that has been proposed to promote an opening of the chromatin, thereby facilitating binding of proteins to the DNA helix. In the mouse genome, however, a similar aggregation of genes coding for DNA binding proteins and high GC content cannot be found. This suggests that the distribution of genes coding for DNA binding proteins and the variations of the chromatin accessibility to these proteins are different in the human and mouse genomes. It is likely that the overall high synonymous and intron rates in chromosome 19 are a by-product of the high GC content of this chromosome.

Lineage-specific expansions of retroviral insertions within the genomes of African great apes but not humans and orangutans

Retroviral infections of the germline have the potential to episodically alter gene function and genome structure during the course of evolution. Horizontal transmissions between species have been proposed, but little evidence exists for such events in the human/great ape lineage of evolution. Based on analysis of finished BAC chimpanzee genome sequence, we characterize a retroviral element (Pan troglodytes endogenous retrovirus 1 [PTERV1]) that has become integrated in the germline of African great ape and Old World monkey species but is absent from humans and Asian ape genomes. We unambiguously map 287 retroviral integration sites and determine that approximately 95.8% of the insertions occur at non-orthologous regions between closely related species. Phylogenetic analysis of the endogenous retrovirus reveals that the gorilla and chimpanzee elements share a monophyletic origin with a subset of the Old World monkey retroviral elements, but that the average sequence divergence exceeds neutral expectation for a strictly nuclear inherited DNA molecule. Within the chimpanzee, there is a significant integration bias against genes, with only 14 of these insertions mapping within intronic regions. Six out of ten of these genes, for which there are expression data, show significant differences in transcript expression between human and chimpanzee. Our data are consistent with a retroviral infection that bombarded the genomes of chimpanzees and gorillas independently and concurrently, 3–4 million years ago. We speculate on the potential impact of such recent events on the evolution of humans and great apes.

Comparison of human and other primate genomes provides evidence for a retroviral infection that bombarded the genomes of chimpanzees and gorillas 3-4 million years ago

Fine-scale recombination patterns differ between chimpanzees and humans

Recombination rates seem to vary extensively along the human genome. Pedigree analysis suggests that rates vary by an order of magnitude when measured at the megabase scale, and at a finer scale, sperm typing studies point to the existence of recombination hotspots. These are short regions (1-2 kb) in which recombination rates are 10-1,000 times higher than the background rate. Less is known about how recombination rates change over time. Here we determined to what degree recombination rates are conserved among closely related species by estimating recombination rates from 14 Mb of linkage disequilibrium data in central chimpanzee and human populations. The results suggest that recombination hotspots are not conserved between the two species and that recombination rates in larger (50 kb) genomic regions are only weakly conserved. Therefore, the recombination landscape has changed markedly between the two species.

Comparison of fine-scale recombination rates in humans and chimpanzees

We compared fine-scale recombination rates at orthologous loci in humans and chimpanzees by analyzing polymorphism data in both species. Strong statistical evidence for hotspots of recombination was obtained in both species. Despite approximately 99% identity at the level of DNA sequence, however, recombination hotspots were found rarely (if at all) at the same positions in the two species, and no correlation was observed in estimates of fine-scale recombination rates. Thus, local patterns of recombination rate have evolved rapidly, in a manner disproportionate to the change in DNA sequence.

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.

A comparison of hominid and rodent lineages reveals that the gene control regions of hominids are not conserved and are accumulating mutations, suggesting widespread degradation of the hominid genome

Reduced variation on the chicken Z chromosome

Understanding the population genetic factors that shape genome variability is pivotal to the design and interpretation of studies using large-scale polymorphism data. We analyzed patterns of polymorphism and divergence at Z-linked and autosomal loci in the domestic chicken (Gallus gallus) to study the influence of mutation, effective population size, selection, and demography on levels of genetic diversity. A total of 14 autosomal introns (8316 bp) and 13 Z-linked introns (6856 bp) were sequenced in 50 chicken chromosomes from 10 highly divergent breeds. Genetic variation was significantly lower at Z-linked than at autosomal loci, with one segregating site every 39 bp at autosomal loci (theta(W) = 5.8 +/- 0.8 x 10(-3)) and one every 156 bp on the Z chromosome (theta(W) = 1.4 +/- 0.4 x 10(-3)). This difference may in part be due to a low male effective population size arising from skewed reproductive success among males, evident both in the wild ancestor-the red jungle fowl-and in poultry breeding. However, this effect cannot entirely explain the observed three- to fourfold reduction in Z chromosome diversity. Selection, in particular selective sweeps, may therefore have had an impact on reducing variation on the Z chromosome, a hypothesis supported by the observation of heterogeneity in diversity levels among loci on the Z chromosome and the lower recombination rate on Z than on autosomes. Selection on sex-linked genes may be particularly important in organisms with female heterogamety since the heritability of sex-linked sexually antagonistic alleles advantageous to males is improved when fathers pass a Z chromosome to their sons.

Are genes of human intelligence related to the metabolism of thyroid and steroids hormones? – Endocrine changes may explain human evolution and higher intelligence

We propose the hypothesis that genes of human intelligence are related with metabolism of thyroid and steroids hormones, which have a crucial role in brain development and function. First, there is evidence to support the idea that during hominid evolution small genetic differences were related with significant endocrine changes in thyroid and steroids hormones. Second, these neuroactive hormones are also related with unique features of human evolution such as body and brain size increase, penis and breast enlargement, pelvic sexual dimorphism, active sexuality, relative lack of hair and higher longevity. Besides underling many of the differences between humans and great apes, steroids hormones promote brain growth and development, are important in the myelination process, explain sexual dimorphisms in brain and intelligence and improve specific cognitive abilities in humans. Supporting our hypothesis, recent studies indicate differences in neuroactive hormones metabolism between humans and non-human primates. Furthermore, a link between X chromosome genes and sex steroids may explain why the frequency of genes affecting intelligence is so high on the X chromosome. This association suggests that, during hominid evolution, there was a positive feedback in both sexes on the same genes responsible for secondary sexual character development and intelligence. This interaction leads to acceleration of development of human brain and intelligence. Finally, we propose that neuroactive hormone therapy may provide significant improvement in some cognitive deficits in all stages of human life and in cases of neurodegenerative diseases. However, further investigation is needed, mainly in the enzymatic machinery, in order to understand the direct role of these hormones in intelligence.

CYP3A variation and the evolution of salt-sensitivity variants

Members of the cytochrome P450 3A subfamily catalyze the metabolism of endogenous substrates, environmental carcinogens, and clinically important exogenous compounds, such as prescription drugs and therapeutic agents. In particular, the CYP3A4 and CYP3A5 genes play an especially important role in pharmacogenetics, since they metabolize >50% of the drugs on the market. However, known genetic variants at these two loci are not sufficient to account for the observed phenotypic variability in drug response. We used a comparative genomics approach to identify conserved coding and noncoding regions at these genes and resequenced them in three ethnically diverse human populations. We show that remarkable interpopulation differences exist with regard to frequency spectrum and haplotype structure. The non-African samples are characterized by a marked excess of rare variants and the presence of a homogeneous group of long-range haplotypes at high frequency. The CYP3A5*1/*3 polymorphism, which is likely to influence salt and water retention and risk for salt-sensitive hypertension, was genotyped in >1,000 individuals from 52 worldwide population samples. The results reveal an unusual geographic pattern whereby the CYP3A5*3 frequency shows extreme variation across human populations and is significantly correlated with distance from the equator. Furthermore, we show that an unlinked variant, AGT M235T, previously implicated in hypertension and pre-eclampsia, exhibits a similar geographic distribution and is significantly correlated in frequency with CYP3A5*1/*3. Taken together, these results suggest that variants that influence salt homeostasis were the targets of a shared selective pressure that resulted from an environmental variable correlated with latitude.

Genomic data support the hominoid slowdown and an Early Oligocene estimate for the hominoid-cercopithecoid divergence

Several lines of indirect evidence suggest that hominoids (apes and humans) and cercopithecoids (Old World monkeys) diverged around 23-25 Mya. Importantly, although this range of dates has been used as both an initial assumption and as a confirmation of results in many molecular-clock analyses, it has not been critically assessed on its own merits. In this article we test the robusticity of the 23- to 25-Mya estimate with approximately 150,000 base pairs of orthologous DNA sequence data from two cercopithecoids and two hominoids by using quartet analysis. This method is an improvement over other estimates of the hominoid-cercopithecoid divergence because it incorporates two calibration points, one each within cercopithecoids and hominoids, and tests for a statistically appropriate model of molecular evolution. Most comparisons reject rate constancy in favor of a model incorporating two rates of evolution, supporting the "hominoid slowdown" hypothesis. By using this model of molecular evolution, the hominoid-cercopithecoid divergence is estimated to range from 29.2 to 34.5 Mya, significantly older than most previous analyses. Hominoid-cercopithecoid divergence dates of 23-25 Mya fall outside of the confidence intervals estimated, suggesting that as much as one-third of ape evolution has not been paleontologically sampled. Identifying stem cercopithecoids or hominoids from this period will be difficult because derived features that define crown catarrhines need not be present in early members of these lineages. More sites that sample primate habitats from the Oligocene of Africa are needed to better understand early ape and Old World monkey evolution.

Comparative primate genomics

With the completion of the human genome sequence and the advent of technologies to study functional aspects of genomes, molecular comparisons between humans and other primates have gained momentum. The comparison of the human genome to the genomes of species closely related to humans allows the identification of genomic features that set primates apart from other mammals and of features that set certain primates notably humans apart from other primates. In this article, we review recent progress in these areas with an emphasis on how comparative approaches may be used to identify functionally relevant features unique to the human genome.

Genetic and haplotype diversity among wild-derived mouse inbred strains

With the completion of the mouse genome sequence, it is possible to define the amount, type, and organization of the genetic variation in this species. Recent reports have provided an overview of the structure of genetic variation among classical laboratory mice. On the other hand, little is known about the structure of genetic variation among wild-derived strains with the exception of the presence of higher levels of diversity. We have estimated the sequence diversity due to substitutions and insertions/deletions among 20 inbred strains of Mus musculus, chosen to enable interpretation of the molecular variation within a clear evolutionary framework. Here, we show that the level of sequence diversity present among these strains is one to two orders of magnitude higher than the level of sequence diversity observed in the human population, and only a minor fraction of the sequence differences observed is found among classical laboratory strains. Our analyses also demonstrate that deletions are significantly more frequent than insertions. We estimate that 50% of the total variation identified in M. musculus may be recovered in intrasubspecific crosses. Alleles at variants positions can be classified into 164 strain distribution patterns, a number exceeding those reported and predicted in panels of classical inbred strains. The number of strains, the analysis of multiple loci scattered across the genome, and the mosaic nature of the genome in hybrid and classical strains contribute to the observed diversity of strain distribution patterns. However, phylogenetic analyses demonstrate that ancient polymorphisms that segregate across species and subspecies play a major role in the generation of strain distribution patterns.

Surveying phylogenetic footprints in large gene clusters: applications to Hox cluster duplications

Evolutionarily conserved non-coding genomic sequences represent a potentially rich source for the discovery of gene regulatory regions. Since these elements are subject to stabilizing selection they evolve much more slowly than adjacent non-functional DNA. These so-called phylogenetic footprints can be detected by comparison of the sequences surrounding orthologous genes in different species. Therefore the loss of phylogenetic footprints as well as the acquisition of conserved non-coding sequences in some lineages, but not in others, can provide evidence for the evolutionary modification of cis-regulatory elements. We introduce here a statistical model of footprint evolution that allows us to estimate the loss of sequence conservation that can be attributed to gene loss and other structural reasons. This approach to studying the pattern of cis-regulatory element evolution, however, requires the comparison of relatively long sequences from many species. We have therefore developed an efficient software tool for the identification of corresponding footprints in long sequences from multiple species. We apply this novel method to the published sequences of HoxA clusters of shark, human, and the duplicated zebrafish and Takifugu clusters as well as the published HoxB cluster sequences. We find that there is a massive loss of sequence conservation in the intergenic region of the HoxA clusters, consistent with the finding in [Chiu et al., PNAS 99 (2002) 5492]. The loss of conservation after cluster duplication is more extensive than expected from structural reasons. This suggests that binding site turnover and/or adaptive modification may also contribute to the loss of sequence conservation.

Comparative analyses reveal a complex history of molecular evolution for human MYH16

We describe the pattern of molecular evolution at a sarcomeric myosin gene, MYH16, using more than 30,000 bp of exon and intron sequence data from the chimpanzee and human genome sequencing projects to evaluate the timing and consequences of a human lineage-specific frameshift deletion. We estimate the age of the deletion at approximately 5.3 MYA. This estimate is consistent with the time of human and chimpanzee divergence and is significantly older than the first appearance of the genus Homo in the fossil record. We also find conflicting estimates of nonsynonymous fixation rates (d(N)) across different regions of this gene, revealing a complex pattern inconsistent with a simple model of pseudogene evolution for human MYH16.

Characterization of a common deletion polymorphism of the UGT2B17 gene linked to UGT2B15

Members of the human UDP-glucuronosyltransferase 2B family are located in a cluster on chromosome 4q13 and code for enzymes whose gene products are responsible for the normal catabolism of steroid hormones. Two members of this family, UGT2B15 and UGT2B17, share over 95% sequence identity. However, UGT2B17 exhibits broader substrate specificity due to a single amino acid difference. Using gene-specific primers to explore the genomic organization of these two genes, it was determined that UGT2B17 is absent in some human DNA samples. The gene-specific primers demonstrated the presence or absence of a 150 kb genomic interval spanning the entire UGT2B17 gene, revealing that UGT2B17 is present in the human genome as a deletion polymorphism linked to UGT2B15. Furthermore, it is shown that the UGT2B17 deletion polymorphism shows Mendelian segregation and allele frequencies that differ between African Americans and Caucasians.

Relaxation of selective constraint and loss of function in the evolution of human bitter taste receptor genes

Bitter taste perception prevents mammals from ingesting poisonous substances because many toxins taste bitter and cause aversion. We hypothesize that human bitter taste receptor (TAS2R) genes might be relaxed from selective constraints because of the change in diet, use of fire and reliance on other means of toxin avoidance that emerged in human evolution. Here, we examine the intra-specific variations of all 25 genes of the human TAS2R repertoire. Our data show hallmarks of neutral evolution, including similar rates of synonymous (d(S)) and non-synonymous (d(N)) nucleotide changes among rare polymorphisms, common polymorphisms and substitutions; no variation in d(N)/d(S) among functional domains; segregation of pseudogene alleles within species and fixation of loss-of-function mutations. These results, together with previous findings of large numbers of loss-of-function mutations in olfactory, pheromonal and visual sensory genes in humans, suggest surprisingly reduced sensory capabilities of humans in comparison with many other mammals.

Positive selection in MAOA gene is human exclusive: determination of the putative amino acid change selected in the human lineage

Monoamine oxidase A (MAOA) is the X-linked gene responsible for deamination and subsequent degradation of several neurotransmitters and other amines. Among other activities, the gene has been shown to play a role in locomotion, circadian rhythm, and pain sensitivity and to have a critical influence on behavior and cognition. Previous studies have reported a non-neutral evolution of the gene attributable to positive selection in the human lineage. To determine whether this selection was human-exclusive or shared with other species, we performed a population genetic analysis of the pattern of nucleotide variation in non-human species, including bonobo, chimpanzee, gorilla, and orangutan. Footprints of positive selection were absent in all analyzed species, suggesting that positive selection has been recent and unique to humans. To determine which human-unique genetic changes could have been responsible for this differential evolution, the coding region of the gene was compared between human, chimpanzee, and gorilla. Only one human exclusive non-conservative change is present in the gene: Glu151Lys. This human substitution affects protein dimerization according to a three-dimensional structural model that predicts a non-negligible functional shift. This is the only candidate position at present to have been selected to fixation in humans during an episode of positive selection. Divergence analysis among species has shown that, even under positive selection in the human lineage, the MAOA gene did not experience accelerated evolution in any of the analyzed lineages, and that tools such as K(a)/ K(s) would not have detected the selective history of the gene.

Testing the chromosomal speciation hypothesis for humans and chimpanzees

Fixed differences of chromosomal rearrangements between isolated populations may promote speciation by preventing between-population gene flow upon secondary contact, either because hybrids suffer from lowered fitness or, more likely, because recombination is reduced in rearranged chromosomal regions. This chromosomal speciation hypothesis thus predicts more rapid genetic divergence on rearranged than on colinear chromosomes because the former are less porous to gene flow. A number of studies of fungi, plants, and animals, including limited genetic data of humans and chimpanzees, support the hypothesis. Here we reexamine the hypothesis for humans and chimpanzees with substantially more genomic data than were used previously. No difference is observed between rearranged and colinear chromosomes in the level of genomic DNA sequence divergence between species. The same is also true for protein sequences. When the gorilla is used as an outgroup, no acceleration in protein sequence evolution associated with chromosomal rearrangements is found. Furthermore, divergence in expression pattern between orthologous genes is not significantly different for rearranged and colinear chromosomes. These results, showing that chromosomal rearrangements did not affect the rate of genetic divergence between humans and chimpanzees, are expected if incipient species on the evolutionary lineages separating humans and chimpanzees did not hybridize.

Evolutionary dynamics of duplicated microsatellites shared by sex chromosomes

Segmental duplications on sex chromosomes constitute an important proportion of recent duplications (approximately 30%). Among those, the evolution of duplicated noncoding DNA is still poorly investigated. We focus our work on repeated DNA sequences extensively used in population genetics and evolution: microsatellites. Six duplicated (CA), microsatellite loci, located on the homologous region of human sex chromosomes, were studied at the intraspecific level in Homo sapiens and by an orthologous comparison in eight primate species. At the intraspecific level, we evaluated the congruence in paralogous divergence between the flanking sequences of the six microsatellites and the approximately 2.2-kb surrounding sequences and observed that both phylogenies are congruent. At the interspecific level (8 species of primates: 54 individuals), we analyzed the sequence polymorphism and divergence of each orthologous locus for both the flanking sequence and the microsatellite. The results showed a lower divergence of flanking sequences than expected in noncoding DNA and a relative stability of the first nucleotides close to the microsatellite. The location of each CAIII locus in a Low Copy Repeated element containing duplicated VCX/Y genes (approximately 1 kb) suggested that direct or indirect selection could explain these results. Moreover, the substitution rates in the flanking sequences and in the microsatellites were correlated. Thus, the evolutionary dynamics of microsatellites seems closely linked to the variation of spontaneous mutations in the surrounding regions.

What makes us human?

The sequence of chimpanzee chromosome 22 is starting to help us to define the set of genetic attributes that are unique to humans, but interpreting the biological consequences of these remains a major challenge.

Genomic regionality in rates of evolution is not explained by clustering of genes of comparable expression profile

In mammalian genomes, linked genes show similar rates of evolution, both at fourfold degenerate synonymous sites (K4) and at nonsynonymous sites (KA). Although it has been suggested that the local similarity in the synonymous substitution rate is an artifact caused by the inclusion of disparately evolving gene pairs, we demonstrate here that this is not the case: after removal of disparately evolving genes, both (1) linked genes and (2) introns from the same gene have more similar silent substitution rates than expected by chance. What causes the local similarity in both synonymous and nonsynonymous substitution rates? One class of hypotheses argues that both may be related to the observed clustering of genes of comparable expression profile. We investigate these hypotheses using substitution rates from both human-mouse and mouse-rat comparisons, and employing three different methods to assay expression parameters. Although we confirm a negative correlation of expression breadth with both K4 and KA, we find no evidence that clustering of similarly expressed genes explains the clustering of genes of comparable substitution rates. If gene expression is not responsible, what about other causes? At least in the human-mouse comparison, the local similarity in KA can be explained by the covariation of KA and K4. As regards K4, our results appear consistent with the notion that local similarity is due to processes associated with meiotic recombination.

Origins of chromosomal rearrangement hotspots in the human genome: evidence from the AZFa deletion hotspots

BACKGROUND

The origins of the recombination hotspots that are a common feature of both allelic and non-allelic homologous recombination in the human genome are poorly understood. We have investigated, by comparative sequencing, the evolution of two hotspots of non-allelic homologous recombination on the Y chromosome that lie within paralogous sequences known to sponsor deletions resulting in male infertility.

RESULTS

These recombination hotspots are characterized by signatures of concerted evolution, which indicate that gene conversion between paralogs has been predominant in shaping their recent evolution. By contrast, the paralogous sequences that surround the hotspots exhibit little evidence of gene conversion. A second feature of these rearrangement hotspots is the extreme interspecific sequence divergence (around 2.5%) that places them among the most divergent orthologous sequences between humans and chimpanzees.

CONCLUSIONS

Several hominid-specific gene conversion events have rendered these hotspots better substrates for chromosomal rearrangements in humans than in chimpanzees or gorillas. Monte Carlo simulations of sequence evolution suggest that extreme sequence divergence is a direct consequence of gene conversion between paralogs. We propose that the coincidence of signatures of concerted evolution and recurrent breakpoints of chromosomal rearrangement (mapped at the sequence level) may enable the identification of putative rearrangement hotspots from analysis of comparative sequences from great apes.

Single nucleotide polymorphism (SNP) discovery in mammals: a targeted-gene approach

Single nucleotide polymorphisms (SNPs) have rarely been exploited in nonhuman and nonmodel organism genetic studies. This is due partly to difficulties in finding SNPs in species where little DNA sequence data exist, as well as to a lack of robust and inexpensive genotyping methods. We have explored one SNP discovery method for molecular ecology, evolution, and conservation studies to evaluate the method and its limitations for population genetics in mammals. We made use of 'CATS' (or 'EPIC') primers to screen for novel SNPs in mammals. Most of these primer sets were designed from primates and/or rodents, for amplifying intron regions from conserved genes. We have screened 202 loci in 16 representatives of the major mammalian clades. Polymerase chain reaction (PCR) success correlated with phylogenetic distance from the human and mouse sequences used to design most primers; for example, specific PCR products from primates and the mouse amplified the most consistently and the marsupial and armadillo amplifications were least successful. Approximately 24% (opossum) to 65% (chimpanzee) of primers produced usable PCR product(s) in the mammals tested. Products produced generally high but variable levels of readable sequence and similarity to the expected genes. In a preliminary screen of chimpanzee DNA, 12 SNPs were identified from six (of 11) sequenced regions, yielding a SNP on average every 400 base pairs (bp). Given the progress in genome sequencing, and the large numbers of CATS-like primers published to date, this approach may yield sufficient SNPs per species for population and conservation genetic studies in nonmodel mammals and other organisms.

Mutation rate variation in the mammalian genome

Recent advances in the large-scale sequencing of mammalian genomes have provided a means to study divergence in not only genic sequences but also in the non-coding bulk of DNA. There is evidence of significant variation in the levels of divergence between presumably neutral regions, pointing at an underlying variation in the rate of mutation across the genome. Apparently, such variation occurs on different scales, including sequence context effects (the influence of neighboring nucleotides on the rate of mutation at individual sites), variation within chromosomes (on the scales of kilobases as well as megabases), and between chromosomes (among autosomes as well as between autosomes and sex chromosomes). An important aspect for further research in this area is to study whether there is an ultimate evolutionary explanation for mutation rate variation within mammalian genomes.

DNA sequence and comparative analysis of chimpanzee chromosome 22

Human-chimpanzee comparative genome research is essential for narrowing down genetic changes involved in the acquisition of unique human features, such as highly developed cognitive functions, bipedalism or the use of complex language. Here, we report the high-quality DNA sequence of 33.3 megabases of chimpanzee chromosome 22. By comparing the whole sequence with the human counterpart, chromosome 21, we found that 1.44% of the chromosome consists of single-base substitutions in addition to nearly 68,000 insertions or deletions. These differences are sufficient to generate changes in most of the proteins. Indeed, 83% of the 231 coding sequences, including functionally important genes, show differences at the amino acid sequence level. Furthermore, we demonstrate different expansion of particular subfamilies of retrotransposons between the lineages, suggesting different impacts of retrotranspositions on human and chimpanzee evolution. The genomic changes after speciation and their biological consequences seem more complex than originally hypothesized.

Absence of the TAP2 human recombination hotspot in chimpanzees

Recent experiments using sperm typing have demonstrated that, in several regions of the human genome, recombination does not occur uniformly but instead is concentrated in “hotspots” of 1–2 kb. Moreover, the crossover asymmetry observed in a subset of these has led to the suggestion that hotspots may be short-lived on an evolutionary time scale. To test this possibility, we focused on a region known to contain a recombination hotspot in humans, TAP2, and asked whether chimpanzees, the closest living evolutionary relatives of humans, harbor a hotspot in a similar location. Specifically, we used a new statistical approach to estimate recombination rate variation from patterns of linkage disequilibrium in a sample of 24 western chimpanzees (Pan troglodytes verus). This method has been shown to produce reliable results on simulated data and on human data from the TAP2 region. Strikingly, however, it finds very little support for recombination rate variation at TAP2 in the western chimpanzee data. Moreover, simulations suggest that there should be stronger support if there were a hotspot similar to the one characterized in humans. Thus, it appears that the human TAP2 recombination hotspot is not shared by western chimpanzees. These findings demonstrate that fine-scale recombination rates can change between very closely related species and raise the possibility that rates differ among human populations, with important implications for linkage-disequilibrium based association studies.

The human TAP2 recombination hotspot is absent from the homologous region in western chimpanzees, with important implications for association studies, the HapMap project and understanding fine-scale variation in recombination rates

Selective pressures on the olfactory receptor repertoire since the human-chimpanzee divergence

The availability of the sequence of the chimpanzee genome provides an opportunity to examine human genes and their chimpanzee orthologs and to analyze selective pressures that have been shaping the olfactory receptor repertoire since the human-chimpanzee divergence. We determined the ratio of nonsynonymous to synonymous changes for each of 186 orthologous pairs and then examined how the distribution of these ratios compares with the distribution expected under neutral drift. Consistent with the diminishing importance of olfaction for these species, we find no evidence for positive selection and we find evidence of weak purifying selection affecting over half of the repertoire.

Male-biased mutation rate and divergence in autosomal, z-linked and w-linked introns of chicken and Turkey

To investigate mutation-rate variation between autosomes and sex chromosomes in the avian genome, we have analyzed divergence between chicken (Gallus gallus) and turkey (Meleagris galopavo) sequences from 33 autosomal, 28 Z-linked, and 14 W-linked introns with a total ungapped alignment length of approximately 43,000 bp. There are pronounced differences in the mean divergence among autosomes and sex chromosomes (autosomes [A] = 10.08%, Z chromosome = 10.99%, and W chromosome = 5.74%), and we use these data to estimate the male-to-female mutation-rate ratio (alpha(m)) from Z/A, Z/W, and A/W comparisons at 1.71, 2.37, and 2.52, respectively. Because the alpha(m) estimates of the three comparisons do not differ significantly, we find no statistical support for a specific reduction in the Z chromosome mutation rate (Z reduction estimated at 4.89%, P = 0.286). The idea of mutation-rate reduction in the sex chromosome hemizygous in one sex (i.e., X in mammals, Z in birds) has been suggested on the basis of theory on adaptive mutation-rate evolution. If it exists in birds, the effect would, thus, seem to be weak; a preliminary power analysis suggests that it is significantly less than 18%. Because divergence may vary within chromosomal classes as a result of variation in mutation and/or selection, we developed a novel double-bootstrapping method, bootstrapping both by introns and sites from concatenated alignments, to estimate confidence intervals for chromosomal class rates and for alpha(m). The narrowest interval for the alpha(m) estimate is 1.88 to 2.97 from the Z/W comparison. We also estimated alpha(m) using maximum likelihood on data from all three chromosome classes; this method yielded alpha(m) = 2.47 and approximate 95% confidence intervals of 2.27 to 2.68. Our data are broadly consistent with the idea that mutation-rate differences between chromosomal classes can be explained by the male mutation bias alone.

Insertions and deletions are male biased too: a whole-genome analysis in rodents

It is presently accepted that, in mammals, due to the greater number of cell divisions in the male germline than in the female germline, nucleotide substitutions occur more frequently in males. The data on mutation bias in insertions and deletions (indels) are contradictory, with some studies indicating no sex bias and others indicating either female or male bias. The sequenced rat and mouse genomes provide a unique opportunity to investigate a potential sex bias for different types of mutations. Indeed, mutation rates can be accurately estimated from a large number of orthologous loci in organisms similar in generation time and in the number of germline cell divisions. Here we compare the mutation rates between chromosome X and autosomes for likely neutral sites in eutherian ancestral interspersed repetitive elements present at orthologous locations in the rat and mouse genomes. We find that small indels are male biased: The male-to-female mutation rate ratio (alpha) for indels in rodents is approximately 2. Similarly, our whole-genome analysis in rodents indicates an approximately twofold excess of nucleotide substitutions originating in males over that in females. This is the same as the male-to-female ratio of the number of germline cell divisions in rat and mouse. Thus, this is consistent with nucleotide substitutions and small indels occurring primarily during DNA replication.

Pash: efficient genome-scale sequence anchoring by Positional Hashing

Pash is a computer program for efficient, parallel, all-against-all comparison of very long DNA sequences. Pash implements Positional Hashing, a novel parallelizable method for sequence comparison based on k-mer representation of sequences. The Positional Hashing method breaks the comparison problem in a unique way that avoids the quadratic penalty encountered with other sensitive methods and confers inherent low-level parallelism. Furthermore, Positional Hashing allows one to readily and predictably trade between sensitivity and speed. In a simulated comparison task, anchoring computationally mutated reads onto a genome, the sensitivity of Pash was equal to or greater than that of BLAST and BLAT, with Pash outperforming these programs as the reads became shorter and less similar to the genome. Using modest computing resources, we employed Pash for two large-scale sequence comparison tasks: comparison of three mammalian genomes, and anchoring millions of chimpanzee whole-genome shotgun sequencing reads onto the human genome. The results of these comparisons by Pash agree with those computed by other methods that use more than an order of magnitude more computing resources. These results confirm the sensitivity of Positional Hashing.

Inversion, duplication, and changes in gene context are associated with human chromosome 18 evolution

Human chromosome 18 differs from its homologues in the great apes by a pericentric inversion. We have identified a chimpanzee bacterial artificial chromosome that spans a region where a break is likely to have occurred in a human progenitor and have characterized the corresponding regions in both chimpanzees and humans. Interspecies sequence comparisons indicate that the ancestral break occurred between the genes ROCK1 and USP14. In humans, the inversion places ROCK1 near centromeric heterochromatin and USP14 adjacent to highly repetitive subtelomeric repeats. In addition, we provide evidence for a human segmental duplication that may have provided a mechanism for the inversion.

Inferring nonneutral evolution from human-chimp-mouse orthologous gene trios

Even though human and chimpanzee gene sequences are nearly 99% identical, sequence comparisons can nevertheless be highly informative in identifying biologically important changes that have occurred since our ancestral lineages diverged. We analyzed alignments of 7645 chimpanzee gene sequences to their human and mouse orthologs. These three-species sequence alignments allowed us to identify genes undergoing natural selection along the human and chimp lineage by fitting models that include parameters specifying rates of synonymous and nonsynonymous nucleotide substitution. This evolutionary approach revealed an informative set of genes with significantly different patterns of substitution on the human lineage compared with the chimpanzee and mouse lineages. Partitions of genes into inferred biological classes identified accelerated evolution in several functional classes, including olfaction and nuclear transport. In addition to suggesting adaptive physiological differences between chimps and humans, human-accelerated genes are significantly more likely to underlie major known Mendelian disorders.