Combining sperm typing and linkage disequilibrium analyses reveals differences in selective pressures or recombination rates across human populations

A previous polymorphism survey of the type 2 diabetes gene CAPN10 identified a segment showing an excess of polymorphism levels in all population samples, coinciding with localized breakdown of linkage disequilibrium (LD) in a sample of Hausa from Cameroon, but not in non-African samples. This raised the possibility that a recombination hotspot is present in all populations and we had insufficient power to detect it in the non-African data. To test this possibility, we estimated the crossover rate by sperm typing in five non-African men; these estimates were consistent with the LD decay in the non-African, but not in the Hausa data. Moreover, resequencing the orthologous region in a sample of Western chimpanzees did not show either an excess of polymorphism level or rapid LD decay, suggesting that the processes underlying the patterns observed in humans operated only on the human lineage. These results suggest that a hotspot of recombination has recently arisen in humans and has reached higher frequency in the Hausa than in non-Africans, or that there is no elevation in crossover rate in any human population, and the observed variation results from long-standing balancing selection.     

Bayesian Inference of Shared Recombination Hotspots Between Humans and Chimpanzees

Recombination generates variation and facilitates evolution. Recombination (or lack thereof) also contributes to human genetic disease. Methods for mapping genes influencing complex genetic diseases via association rely on linkage disequilibrium (LD) in human populations, which is influenced by rates of recombination across the genome. Comparative population genomic analyses of recombination using related primate species can identify factors influencing rates of recombination in humans. Such studies can indicate how variable hotspots for recombination may be both among individuals (or populations) and over evolutionary timescales. Previous studies have suggested that locations of recombination hotspots are not conserved between humans and chimpanzees. We made use of the data sets from recent resequencing projects and applied a Bayesian method for identifying hotspots and estimating recombination rates. We also reanalyzed SNP data sets for regions with known hotspots in humans using samples from the human and chimpanzee. The Bayes factors (BF) of shared recombination hotspots between human and chimpanzee across regions were obtained. Based on the analysis of the aligned regions of human chromosome 21, locations where the two species show evidence of shared recombination hotspots (with high BFs) were identified. Interestingly, previous comparative studies of human and chimpanzee that focused on the known human recombination hotspots within the β-globin and HLA regions did not find overlapping of hotspots. Our results show high BFs of shared hotspots at locations within both regions, and the estimated locations of shared hotspots overlap with the locations of human recombination hotspots obtained from sperm-typing studies.     

Accelerated Evolution of the ASPM Gene Controlling Brain Size Begins Prior to Human Brain Expansion

Primary microcephaly (MCPH) is a neurodevelopmental disorder characterized by global reduction in cerebral cortical volume. The microcephalic brain has a volume comparable to that of early hominids, raising the possibility that some MCPH genes may have been evolutionary targets in the expansion of the cerebral cortex in mammals and especially primates. Mutations in ASPM, which encodes the human homologue of a fly protein essential for spindle function, are the most common known cause of MCPH. Here we have isolated large genomic clones containing the complete ASPM gene, including promoter regions and introns, from chimpanzee, gorilla, orangutan, and rhesus macaque by transformation-associated recombination cloning in yeast. We have sequenced these clones and show that whereas much of the sequence of ASPM is substantially conserved among primates, specific segments are subject to high Ka/Ks ratios (nonsynonymous/synonymous DNA changes) consistent with strong positive selection for evolutionary change. The ASPM gene sequence shows accelerated evolution in the African hominoid clade, and this precedes hominid brain expansion by several million years. Gorilla and human lineages show particularly accelerated evolution in the IQ domain of ASPM. Moreover, ASPM regions under positive selection in primates are also the most highly diverged regions between primates and nonprimate mammals. We report the first direct application of TAR cloning technology to the study of human evolution. Our data suggest that evolutionary selection of specific segments of the ASPM sequence strongly relates to differences in cerebral cortical size.     

Recombination hotspots: Models and tools for detection

Recombination hotspots are the regions within the genome where the rate, and the frequency of recombination are optimum with a size varying from 1 to 2 kb. The recombination event is mediated by the double-stranded break formation, guided by the combined enzymatic action of DNA topoisomerase and Spo 11 endonuclease. These regions are distributed non-uniformly throughout the human genome and cause distortions in the genetic map. Numerous lines of evidence suggest that the number of hotspots known in humans has increased manifold in recent years. A few facts about the hotspot evolutions were also put forward, indicating the differences in the hotspot position between chimpanzees and humans. In mice, recombination hot spots were found to be clustered within the major histocompatibility complex (MHC) region. Several models, that help explain meiotic recombination has been proposed. Moreover, scientists also developed some computational tools to locate the hotspot position and estimate their recombination rate in humans is of great interest to population and medical geneticists. Here we reviewed the molecular mechanisms, models and in silico prediction techniques of hot spot residues.     

Recombination rate variation in closely related species

Despite their importance to successful meiosis and various evolutionary processes, meiotic recombination rates sometimes vary within species or between closely related species. For example, humans and chimpanzees share virtually no recombination hotspot locations in the surveyed portion of the genomes. However, conservation of recombination rates between closely related species has also been documented, raising an apparent contradiction. Here, we evaluate how and why conflicting patterns of recombination rate conservation and divergence may be observed, with particular emphasis on features that affect recombination, and the scale and method with which recombination is surveyed. Additionally, we review recent studies identifying features influencing fine-scale and broad-scale recombination patterns and informing how quickly recombination rates evolve, how changes in recombination impact selection and evolution in natural populations, and more broadly, which forces influence genome evolution.     

Similarity in recombination rate estimates highly correlates with genetic differentiation in humans

Recombination varies greatly among species, as illustrated by the poor conservation of the recombination landscape between humans and chimpanzees. Thus, shorter evolutionary time frames are needed to understand the evolution of recombination. Here, we analyze its recent evolution in humans. We calculated the recombination rates between adjacent pairs of 636,933 common single-nucleotide polymorphism loci in 28 worldwide human populations and analyzed them in relation to genetic distances between populations. We found a strong and highly significant correlation between similarity in the recombination rates corrected for effective population size and genetic differentiation between populations. This correlation is observed at the genome-wide level, but also for each chromosome and when genetic distances and recombination similarities are calculated independently from different parts of the genome. Moreover, and more relevant, this relationship is robustly maintained when considering presence/absence of recombination hotspots. Simulations show that this correlation cannot be explained by biases in the inference of recombination rates caused by haplotype sharing among similar populations. This result indicates a rapid pace of evolution of recombination, within the time span of differentiation of modern humans.     

Variation in estimated recombination rates across human populations

Recently it has been reported that recombination hotspots appear to be highly variable between humans and chimpanzees, and there is evidence for between-person variability in hotspots, and evolutionary transience. To understand the nature of variation in human recombination rates, it is important to describe patterns of variability across populations. Direct measurement of recombination rates remains infeasible on a large scale, and population-genetic approaches can be imprecise, and are affected by demographic history. Reports to date have suggested broad similarity in recombination rates at large genomic scales and across human populations. Here, we examine recombination rate estimates at a finer population and genomic scale: 28 worldwide populations and 107 SNPs in a 1 Mb stretch of chromosome 22q. We employ analysis of variance of recombination rate estimates, corrected for differences in effective population size using genome-wide microsatellite mutation rate estimates. We find substantial variation in fine-scale rates between populations, but reduced variation within continental groups. All effects examined (SNP-pair, region, population and interactions) were highly significant. Adjustment for effective population size made little difference to the conclusions. Observed hotspots tended to be conserved across populations, albeit at varying intensities. This holds particularly for populations from the same region, and also to a considerable degree across geographical regions. However, some hotspots appear to be population-specific. Several results from studies on the population history of humans are in accordance with our analysis. Our results suggest that between-population variation in DNA sequences may underly recombination rate variation.     

Evidence for a complex demographic history of chimpanzees

To characterize patterns of genomic variation in central chimpanzees(Pan troglodytes troglodytes) and gain insight into their evolution,we sequenced nine unlinked, intergenic regions, representinga total of 19,000 base pairs, in 14 individuals. When theseDNA sequences are compared with homologous sequences previouslycollected in humans and in western chimpanzees (Pan troglodytesverus), nucleotide diversity is higher in central chimpanzeesthan in western chimpanzees or in humans. Consistent with alarger effective population size of central chimpanzees, levelsof linkage disequilibrium are lower than in humans. Patternsof linkage disequilibrium further suggest that homologous geneconversion may be an important contributor to genetic exchangeat short distances, in agreement with a previous study of thesame DNA sequences in humans. In central chimpanzees, but notin western chimpanzees, the allele frequency spectrum is significantlyskewed towards rare alleles, pointing to population size changesor fine-scale population structure. Strikingly, the extent ofgenetic differentiation between western and central chimpanzeesis much stronger than what is seen between human populations.This suggests that careful attention should be paid to geographicsampling in studies of chimpanzee genetic variation.     

Molecular mapping of a recombination hotspot located in the second intron of the human TAP2 locus.

Recombination across the HLA class II region is not randomly distributed, as indicated by both strong linkage disequilibrium within the 100 kb encompassing the DRB1-DQA1-DQB1 loci and complete equilibrium between TAP1 and TAP2, the closest variant sites of which are < 15 kb. In an attempt to explain these observations, 39 novel polymorphic markers in a region encompassing the TAP, LMP, and DOB genes were used to delineate the site of crossover in 11 class II recombinant chromosomes. SSCP demonstrated that two recombination events occurred within an 850-bp interval in the second intron of TAP2, which separates the variant sites of TAP1 and TAP2. These data indicate the presence of a recombination hotspot, the first to be identified from the analysis of familial transmission in the human major histocompatibility complex. The region of crossover was cloned and sequenced from one of the recombinants, further defining the crossover site to a 138-bp segment nested within the 850-bp region. This represents the most precisely defined region of recombination in the human genome.     

Comparative linkage-disequilibrium analysis of the beta-globin hotspot in primates

Recombination rates vary both across the genome and between different species, but little information is available about the temporal and physical scales over which such rates change. To shed light on these questions, we performed a high-resolution analysis of a genomic region within the beta-globin gene cluster that is known to experience elevated recombination rates in humans. For this purpose, we developed new linkage disequilibrium-based methods that thoroughly search for subsets of the data with unusually high or unusually low estimated values of the population-recombination parameter (4Nr, where N is the effective population size and r is the crossover rate between adjacent base pairs). By resequencing a 15-kb segment in a human population sample, we were able to narrow the recombinational hotspot to a segment <2 kb in length that coincides with the beta-globin replication origin. In addition, we analyzed the orthologous region in samples of rhesus macaques and common chimpanzees. Whereas the analysis of the chimpanzee data is complicated by the sample structure, the macaque data imply that this region may not be a hotspot in that species. These results suggest a time scale for the evolution of hotspots in primates. Furthermore, they allow us to propose diverged sequence elements that may contribute to the differences in the recombinational landscape in the two species.     

A coalescent model of recombination hotspots

Recent experimental findings suggest that the assumption of a homogeneous recombination rate along the human genome is too naive. These findings point to block-structured recombination rates; certain regions (called hotspots) are more prone than other regions to recombination. In this report a coalescent model incorporating hotspot or block-structured recombination is developed and investigated analytically as well as by simulation. Our main results can be summarized as follows: (1) The expected number of recombination events is much lower in a model with pure hotspot recombination than in a model with pure homogeneous recombination, (2) hotspots give rise to large variation in recombination rates along the genome as well as in the number of historical recombination events, and (3) the size of a (nonrecombining) block in the hotspot model is likely to be overestimated grossly when estimated from SNP data. The results are discussed with reference to the current debate about block-structured recombination and, in addition, the results are compared to genome-wide variation in recombination rates. A number of new analytical results about the model are derived.     

Application of coalescent methods to reveal fine-scale rate variation and recombination hotspots

There has been considerable recent interest in understanding the way in which recombination rates vary over small physical distances, and the extent of recombination hotspots, in various genomes. Here we adapt, apply, and assess the power of recently developed coalescent-based approaches to estimating recombination rates from sequence polymorphism data. We apply full-likelihood estimation to study rate variation in and around a well-characterized recombination hotspot in humans, in the beta-globin gene cluster, and show that it provides similar estimates, consistent with those from sperm studies, from two populations deliberately chosen to have different demographic and selectional histories. We also demonstrate how approximate-likelihood methods can be used to detect local recombination hotspots from genomic-scale SNP data. In a simulation study based on 80 100-kb regions, these methods detect 43 out of 60 hotspots (ranging from 1 to 2 kb in size), with only two false positives out of 2000 subregions that were tested for the presence of a hotspot. Our study suggests that new computational tools for sophisticated analysis of population diversity data are valuable for hotspot detection and fine-scale mapping of local recombination rates.     

Molecular evolution of recombination hotspots and highly recombining pseudoautosomal regions in hominoids

We examined the effects of recombination on the molecular evolution of noncoding regions in pseudoautosomal regions (PARs) and recombination hotspots in hominoids. The PAR-linked regions analyzed had on average longer branch lengths than those of the recombination hotspots. Moreover, contrary to previous observations, we found no correlation between recombination rate and silent site divergence in our data set and little change in the GC content during recent hominoid evolution. This suggests that the current rate of recombination is not a good indicator of the past rates of recombination for these highly recombining regions. Furthermore, human recombination hotspots show increased AT to GC substitutions in the human lineage, while no such pattern is detected for PAR-linked regions. Together, these observations suggest that recombination hotspots in hominoids are transient in the evolutionary time-scale. Interestingly, the 16p13.3 recombination hotspot locus violates a local molecular clock, though the locus appears to be noncoding and should evolve neutrally. We hypothesize that sudden changes in recombination rate have caused the changes in substitution rate at this locus.     

Number of ancestral human species: a molecular perspective

Despite the remarkable developments in molecular biology over the past three decades, anthropological genetics has had only limited impact on systematics in human evolution. Genetics offers the opportunity to objectively test taxonomies based on morphology and may be used to supplement conventional approaches to hominid systematics. Our analyses, examining chromosomes and 46 estimates of genetic distance, indicate there may have been only around 4 species on the direct line to modern humans and 5 species in total. This contrasts with current taxonomies recognising up to 23 species.

The genetic proximity of humans and chimpanzees has been used to suggest these species are congeneric. Our analysis of genetic distances between them is consistent with this proposal. It is time that chimpanzees, living humans and all fossil humans be classified in Homo. The creation of new genera can no longer be a solution to the complexities of fossil morphologies. Published genetic distances between common chimpanzees and bonobos, along with evidence for interbreeding, suggest they should be assigned to a single species.

The short distance between humans and chimpanzees also places a strict limit on the number of possible evolutionary side branches that might be recognised on the human lineage. All fossil taxa were genetically very close to each other and likely to have been below congeneric genetic distances seen for many mammals.

Our estimates of genetic divergence suggest that periods of around 2 million years are required to produce sufficient genetic distance to represent speciation. Therefore, Neanderthals and so-called H. erectus were genetically so close to contemporary H. sapiens they were unlikely to have been separate species. Thus, it is likely there was only one species of human (H. sapiens) for most of the last 2 million years. We estimate the divergence time of H. sapiensfrom 16 genetic distances to be around 1.7 Ma which is consistent with evidence for the earliest migration out of Africa. These findings call into question the mitochondrial «African Eve» hypothesis based on a far more recent origin for H. sapiens and show that humans did not go through a bottleneck in their recent evolutionary history.

Given the large offset in evolutionary rates of molecules and morphology seen in human evolution, Homo species are likely to be characterised by high levels of morphological variation and low levels of genetic variability. Thus, molecular data suggest the limits for intraspecific morphological variation used by many palaeoanthropologists have been set too low. The role of phenotypic plasticity has been greatly underestimated in human evolution.

We call into question the use of mtDNA for studies of human evolution. This DNA is under strong selection, which violates the assumption of selective neutrality. This issue should be addressed by geneticists, including a reassessment of its use for molecular clocks. There is a need for greater cooperation between palaeoanthropologists and anthropological geneticists to better understand human evolution and to bring palaeoanthropology into the mainstream of evolutionary biology.     

High diversity at PRDM9 in chimpanzees and bonobos

Background

The PRDM9 locus in mammals has increasingly attracted research attention due to its role in mediating chromosomal recombination and possible involvement in hybrid sterility and hence speciation processes. The aim of this study was to characterize sequence variation at the PRDM9 locus in a sample of our closest living relatives, the chimpanzees and bonobos.

Methodology/Principal Findings

PRDM9 contains a highly variable and repetitive zinc finger array. We amplified this domain using long-range PCR and determined the DNA sequences using conventional Sanger sequencing. From 17 chimpanzees representing three subspecies and five bonobos we obtained a total of 12 alleles differing at the nucleotide level. Based on a data set consisting of our data and recently published Pan PRDM9 sequences, we found that at the subspecies level, diversity levels did not differ among chimpanzee subspecies or between chimpanzee subspecies and bonobos. In contrast, the sample of chimpanzees harbors significantly more diversity at PRDM9 than samples of humans. Pan PRDM9 shows signs of rapid evolution including no alleles or ZnFs in common with humans as well as signals of positive selection in the residues responsible for DNA binding.

Conclusions and Significance

The high number of alleles specific to the genus Pan, signs of positive selection in the DNA binding residues, and reported lack of conservation of recombination hotspots between chimpanzees and humans suggest that PRDM9 could be active in hotspot recruitment in the genus Pan. Chimpanzees and bonobos are considered separate species and do not have overlapping ranges in the wild, making the presence of shared alleles at the amino acid level between the chimpanzee and bonobo species interesting in view of the hypothesis that PRDM9 plays a universal role in interspecific hybrid sterility.     

The Case for Selection at CCR5-Δ32

The C-C chemokine receptor 5, 32 base-pair deletion (CCR5-Δ32) allele confers strong resistance to infection by the AIDS virus HIV. Previous studies have suggested that CCR5-Δ32 arose within the past 1,000 y and rose to its present high frequency (5%–14%) in Europe as a result of strong positive selection, perhaps by such selective agents as the bubonic plague or smallpox during the Middle Ages. This hypothesis was based on several lines of evidence, including the absence of the allele outside of Europe and long-range linkage disequilibrium at the locus. We reevaluated this evidence with the benefit of much denser genetic maps and extensive control data. We find that the pattern of genetic variation at CCR5-Δ32 does not stand out as exceptional relative to other loci across the genome. Moreover using newer genetic maps, we estimated that the CCR5-Δ32 allele is likely to have arisen more than 5,000 y ago. While such results can not rule out the possibility that some selection may have occurred at C-C chemokine receptor 5 (CCR5), they imply that the pattern of genetic variation seen atCCR5-Δ32 is consistent with neutral evolution. More broadly, the results have general implications for the design of future studies to detect the signs of positive selection in the human genome.     

Demographic history and genetic differentiation in apes

Comparisons of genetic variation between humans and great apes are hampered by the fact that we still know little about the demographics and evolutionary history of the latter species. In addition, characterizing ape genetic variation is important because they are threatened with extinction, and knowledge about genetic differentiation among groups may guide conservation efforts. We sequenced multiple intergenic autosomal regions totaling 22,400 base pairs (bp) in ten individuals each from western, central, and eastern chimpanzee groups and in nine bonobos, and 16,000 bp in ten Bornean and six Sumatran orangutans. These regions are analyzed together with homologous information from three human populations and gorillas. We find that whereas orangutans have the highest diversity, western chimpanzees have the lowest, and that the demographic histories of most groups differ drastically. Special attention should therefore be paid to sampling strategies and the statistics chosen when comparing levels of variation within and among groups. Finally, we find that the extent of genetic differentiation among "subspecies" of chimpanzees and orangutans is comparable to that seen among human populations, calling the validity of the "subspecies" concept in apes into question.     

Live hot, die young: transmission distortion in recombination hotspots

There is strong evidence that hotspots of meiotic recombination in humans are transient features of the genome. For example, hotspot locations are not shared between human and chimpanzee. Biased gene conversion in favor of alleles that locally disrupt hotspots is a possible explanation of the short lifespan of hotspots. We investigate the implications of such a bias on human hotspots and their evolution. Our results demonstrate that gene conversion bias is a sufficiently strong force to produce the observed lack of sharing of intense hotspots between species, although sharing may be much more common for weaker hotspots. We investigate models of how hotspots arise, and find that only models in which hotspot alleles do not initially experience drive are consistent with observations of rather hot hotspots in the human genome. Mutations acting against drive cannot successfully introduce such hotspots into the population, even if there is direct selection for higher recombination rates, such as to ensure correct segregation during meiosis. We explore the impact of hotspot alleles on patterns of haplotype variation, and show that such alleles mask their presence in population genetic data, making them difficult to detect.     

Similar numbers but different repertoires of olfactory receptor genes in humans and chimpanzees

Animals recognize their external world through the detection of tens of thousands of chemical odorants. Olfactory receptor (OR) genes encode proteins for detecting odorant molecules and form the largest multigene family in mammals. It is known that humans have fewer OR genes and a higher fraction of OR pseudogenes than mice or dogs. To investigate whether these features are human specific or common to all higher primates, we identified nearly complete sets of OR genes from the chimpanzee and macaque genomes and compared them with the human OR genes. In contrast to previous studies, here we show that the number of OR genes ( approximately 810) and the fraction of pseudogenes (51%) in chimpanzees are very similar to those in humans, though macaques have considerably fewer OR genes. The pseudogenization rates and the numbers of genes affected by positive selection are also similar between humans and chimpanzees. Moreover, the most recent common ancestor between humans and chimpanzees had a larger number of functional OR genes (>500) and a lower fraction of pseudogenes (41%) than its descendents, suggesting that the OR gene repertoires are in a phase of deterioration in both lineages. Interestingly, despite the close evolutionary relationship between the 2 species, approximately 25% of their functional gene repertoires are species specific due to massive gene losses. These findings suggest that the tempo of evolution of OR genes is similar between humans and chimpanzees, but the OR gene repertoires are quite different between them. This difference might be responsible for the species-specific ability of odor perception.     

Evidence that a RecQ helicase slows senescence by resolving recombining telomeres

RecQ helicases, including Saccharomyces cerevisiae Sgs1p and the human Werner syndrome protein, are important for telomere maintenance in cells lacking telomerase activity. How maintenance is accomplished is only partly understood, although there is evidence that RecQ helicases function in telomere replication and recombination. Here we use two-dimensional gel electrophoresis (2DGE) and telomere sequence analysis to explore why cells lacking telomerase and Sgs1p (tlc1 sgs1 mutants) senesce more rapidly than tlc1 mutants with functional Sgs1p. We find that apparent X-shaped structures accumulate at telomeres in senescing tlc1 sgs1 mutants in a RAD52- and RAD53-dependent fashion. The X-structures are neither Holliday junctions nor convergent replication forks, but instead may be recombination intermediates related to hemicatenanes. Direct sequencing of examples of telomere I-L in senescing cells reveals a reduced recombination frequency in tlc1 sgs1 compared with tlc1 mutants, indicating that Sgs1p is needed for tlc1 mutants to complete telomere recombination. The reduction in recombinants is most prominent at longer telomeres, consistent with a requirement for Sgs1p to generate viable progeny following telomere recombination. We therefore suggest that Sgs1p may be required for efficient resolution of telomere recombination intermediates, and that resolution failure contributes to the premature senescence of tlc1 sgs1 mutants.