Population dynamics of GC-changing mutations in humans and great apes

The nucleotide composition of the genome is a balance between the origin and fixation rates of different mutations. For example, it is well-known that transitions occur more frequently than transversions, particularly at CpG sites. Differences in fixation rates of mutation types are less explored. Specifically, recombination-associated GC-biased gene conversion (gBGC) may differentially impact GC-changing mutations, due to differences in their genomic distributions and efficiency of mismatch repair mechanisms. Given that recombination evolves rapidly across species, we explore gBGC of different mutation types across human populations and great ape species. We report a stronger correlation between segregating GC frequency and recombination for transitions than for transversions. Notably, CpG transitions are most strongly affected by gBGC in humans and chimpanzees. We show that the overall strength of gBGC is generally correlated with effective population sizes in humans, with some notable exceptions, such as a stronger effect of gBGC on non-CpG transitions in populations of European descent. Furthermore, species of the Gorilla and Pongo genus have a greatly reduced gBGC effect on CpG sites. We also study the dependence of gBGC dynamics on flanking nucleotides and show that some mutation types evolve in opposition to the gBGC expectation, likely due to the hypermutability of specific nucleotide contexts. Our results highlight the importance of different gBGC dynamics experienced by GC-changing mutations and their impact on nucleotide composition evolution.     

Biased Gene Conversion Skews Allele Frequencies in Human Populations,Increasing the Disease Burden of Recessive Alleles

Gene conversion results in the nonreciprocal transfer of genetic information between two recombining sequences, and there is evidence that this process is biased toward G and C alleles. However, the strength of GC-biased gene conversion (gBGC) in human populations and its effects on hereditary disease have yet to be assessed on a genomic scale. Using high-coverage whole-genome sequences of African hunter-gatherers, agricultural populations, and primate outgroups, we quantified the effects of GC-biased gene conversion on population genomic data sets. We find that genetic distances (F_ST and population branch statistics) are modified by gBGC. In addition, the site frequency spectrum is left-shifted when ancestral alleles are favored by gBGC and right-shifted when derived alleles are favored by gBGC. Allele frequency shifts due to gBGC mimic the effects of natural selection. As expected, these effects are strongest in high-recombination regions of the human genome. By comparing the relative rates of fixation of unbiased and biased sites, the strength of gene conversion was estimated to be on the order of Nb ≈ 0.05 to 0.09. We also find that derived alleles favored by gBGC are much more likely to be homozygous than derived alleles at unbiased SNPs (+42.2% to 62.8%). This results in a curse of the converted, whereby gBGC causes substantial increases in hereditary disease risks. Taken together, our findings reveal that GC-biased gene conversion has important population genetic and public health implications.     

GC-biased gene conversion and selection affect GC content in the Oryza genus (rice)

Base composition varies among and within eukaryote genomes. Although mutational bias and selection have initially been invoked, more recently GC-biased gene conversion (gBGC) has been proposed to play a central role in shaping nucleotide landscapes, especially in yeast, mammals, and birds. gBGC is a kind of meiotic drive in favor of G and C alleles, associated with recombination. Previous studies have also suggested that gBGC could be at work in grass genomes. However, these studies were carried on third codon positions that can undergo selection on codon usage. As most preferred codons end in G or C in grasses, gBGC and selection can be confounded. Here we investigated further the forces that might drive GC content evolution in the rice genus using both coding and noncoding sequences. We found that recombination rates correlate positively with equilibrium GC content and that selfing species (Oryza sativa and O. glaberrima) have significantly lower equilibrium GC content compared with more outcrossing species. As recombination is less efficient in selfing species, these results suggest that recombination drives GC content. We also detected a positive relationship between expression levels and GC content in third codon positions, suggesting that selection favors codons ending with G or C bases. However, the correlation between GC content and recombination cannot be explained by selection on codon usage alone as it was also observed in noncoding positions. Finally, analyses of polymorphism data ruled out the hypothesis that genomic variation in GC content is due to mutational processes. Our results suggest that both gBGC and selection on codon usage affect GC content in the Oryza genus and likely in other grass species.     

A comparative analysis of regulatory regions of the transthyretin gene in the mouse, human, and chimpanzee genomes

A full genome analysis of differences between the gene expression in the human and chimpanzee brains revealed that the gene for transthyretin, the carrier of thyroid hormones, is differently transcribed in the cerebella of these species. A 7-kbp DNA fragment of chimpanzee was sequenced to identify possible regulatory sequences responsible for the differences in expression. One hundred and thirteen substitutions were found in the chimpanzee sequence in comparison with the human sequence. About 40% of the substitutions were revealed within the repeating elements of the genome; their location and sizes did not differ from those in the corresponding fragments of the human genome, and the nucleotide sequences had a high degree of identity. A comparison of nucleotide sequences of the transthyretin region of human, chimpanzee, and mouse genes revealed substantial differences in the distribution of G + C content along the examined fragment in the human (chimpanzee) and mouse genes and allowed us to localize three sequence tracts with a higher degree of identity in the three species. One of these tracts is located in the promoter region of the gene, and the other two probably determine the specificity of transthyretin gene expression in the liver and brain. One of the conserved tracts of the chimpanzee genome was found to have a single and a triple nucleotide substitution. The triple substitution distinguishes chimpanzees from humans and mice, which have identical sequences of this site. It is likely that these substitutions are responsible for the differences in the expression levels of the transthyretin gene in the human and chimpanzee brains.     

A Comparative Analysis of Regulatory Regions of the Transthyretin Gene in the Mouse, Human, and Chimpanzee Genomes

A full genome analysis of differences between the gene expression in the human and chimpanzee brains revealed that the gene for transthyretin, the carrier of thyroid hormones, is differently transcribed in the cerebella of these species. A 7-kbp DNA fragment of chimpanzee was sequenced to identify possible regulatory sequences responsible for the differences in expression. One hundred and thirteen substitutions were found in the chimpanzee sequence in comparison with the human sequence. About 40% of the substitutions were revealed within the repeating elements of the genome; their location and sizes did not differ from those in the corresponding fragments of the human genome, and the nucleotide sequences had a high degree of identity. A comparison of nucleotide sequences of the transthyretin region of human, chimpanzee, and mouse genes revealed substantial differences in the distribution of G + C content along the examined fragment in the human (chimpanzee) and mouse genes and allowed us to localize three sequence tracts with a higher degree of identity in the three species. One of these tracts was located in the promoter region of the gene, and the other two probably determine the specificity of transthyretin gene expression in the liver and brain. One of the conserved tracts of the chimpanzee genome was found to have a single and a triple nucleotide substitution. The triple substitution distinguishes chimpanzees from humans and mice, which have identical sequences of this site. It is likely that these substitutions are responsible for the differences in the expression levels of the transthyretin gene in the human and chimpanzee brains.     

Detecting positive selection within genomes: the problem of biased gene conversion

The identification of loci influenced by positive selection is a major goal of evolutionary genetics. A popular approach is to perform scans of alignments on a genome-wide scale in order to find regions evolving at accelerated rates on a particular branch of a phylogenetic tree. However, positive selection is not the only process that can lead to accelerated evolution. Notably, GC-biased gene conversion (gBGC) is a recombination-associated process that results in the biased fixation of G and C nucleotides. This process can potentially generate bursts of nucleotide substitutions within hotspots of meiotic recombination. Here, we analyse the results of a scan for positive selection on genes on branches across the primate phylogeny. We show that genes identified as targets of positive selection have a significant tendency to exhibit the genomic signature of gBGC. Using a maximum-likelihood framework, we estimate that more than 20 per cent of cases of significantly elevated non-synonymous to synonymous substitution rates ratio (d_N/d_S), particularly in shorter branches, could be due to gBGC. We demonstrate that in some cases, gBGC can lead to very high d_N/d_S (more than 2). Our results indicate that gBGC significantly affects the evolution of coding sequences in primates, often leading to patterns of evolution that can be mistaken for positive selection.     

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.     

Distribution of Recombination Hotspots in the Human Genome – A Comparison of Computer Simulations with Real Data

Recombination is the main cause of genetic diversity. Thus, errors in this process can lead to chromosomal abnormalities. Recombination events are confined to narrow chromosome regions called hotspots in which characteristic DNA motifs are found. Genomic analyses have shown that both recombination hotspots and DNA motifs are distributed unevenly along human chromosomes and are much more frequent in the subtelomeric regions of chromosomes than in their central parts. Clusters of motifs roughly follow the distribution of recombination hotspots whereas single motifs show a negative correlation with the hotspot distribution. To model the phenomena related to recombination, we carried out computer Monte Carlo simulations of genome evolution. Computer simulations generated uneven distribution of hotspots with their domination in the subtelomeric regions of chromosomes. They also revealed that purifying selection eliminating defective alleles is strong enough to cause such hotspot distribution. After sufficiently long time of simulations, the structure of chromosomes reached a dynamic equilibrium, in which number and global distribution of both hotspots and defective alleles remained statistically unchanged, while their precise positions were shifted. This resembles the dynamic structure of human and chimpanzee genomes, where hotspots change their exact locations but the global distributions of recombination events are very similar.     

The role of GC-biased gene conversion in shaping the fastest evolving regions of the human genome

GC-biased gene conversion (gBGC) is a recombination-associated evolutionary process that accelerates the fixation of guanine or cytosine alleles, regardless of their effects on fitness. gBGC can increase the overall rate of substitutions, a hallmark of positive selection. Many fast-evolving genes and noncoding sequences in the human genome have GC-biased substitution patterns, suggesting that gBGC-in contrast to adaptive processes-may have driven the human changes in these sequences. To investigate this hypothesis, we developed a substitution model for DNA sequence evolution that quantifies the nonlinear interacting effects of selection and gBGC on substitution rates and patterns. Based on this model, we used a series of lineage-specific likelihood ratio tests to evaluate sequence alignments for evidence of changes in mode of selection, action of gBGC, or both. With a false positive rate of less than 5% for individual tests, we found that the majority (76%) of previously identified human accelerated regions are best explained without gBGC, whereas a substantial minority (19%) are best explained by the action of gBGC alone. Further, more than half (55%) have substitution rates that significantly exceed local estimates of the neutral rate, suggesting that these regions may have been shaped by positive selection rather than by relaxation of constraint. By distinguishing the effects of gBGC, relaxation of constraint, and positive selection we provide an integrated analysis of the evolutionary forces that shaped the fastest evolving regions of the human genome, which facilitates the design of targeted functional studies of adaptation in humans.     

Surprising fitness consequences of GC-biased gene conversion. II. Heterosis

Heterosis is a widespread phenomenon corresponding to the increase in fitness following crosses between individuals from different populations or lines relative to their parents. Its genetic basis has been a topic of controversy since the early 20th century. The masking of recessive deleterious mutations in hybrids likely explains a substantial part of heterosis. The dynamics and consequences of these mutations have thus been studied in depth. Recently, it was suggested that GC-biased gene conversion (gBGC) might strongly affect the fate of deleterious mutations and may have significant fitness consequences. gBGC is a recombination-associated process mimicking selection in favor of G and C alleles, which can interfere with selection, for instance by increasing the frequency of GC deleterious mutations. I investigated how gBGC could affect the amount and genetic structure of heterosis through an analysis of the interaction between gBGC and selection in subdivided populations. To do so, I analyzed the infinite island model both by numerical computations and by analytical approximations. I showed that gBGC might have little impact on the total amount of heterosis but could greatly affect its genetic basis.     

High recombination rates and hotspots in a Plasmodium falciparum genetic cross

Background

The human malaria parasite Plasmodium falciparum survives pressures from the host immune system and antimalarial drugs by modifying its genome. Genetic recombination and nucleotide substitution are the two major mechanisms that the parasite employs to generate genome diversity. A better understanding of these mechanisms may provide important information for studying parasite evolution, immune evasion and drug resistance.

Results

Here, we used a high-density tiling array to estimate the genetic recombination rate among 32 progeny of a P. falciparum genetic cross (7G8 × GB4). We detected 638 recombination events and constructed a high-resolution genetic map. Comparing genetic and physical maps, we obtained an overall recombination rate of 9.6 kb per centimorgan and identified 54 candidate recombination hotspots. Similar to centromeres in other organisms, the sequences of P. falciparum centromeres are found in chromosome regions largely devoid of recombination activity. Motifs enriched in hotspots were also identified, including a 12-bp G/C-rich motif with 3-bp periodicity that may interact with a protein containing 11 predicted zinc finger arrays.

Conclusions

These results show that the P. falciparum genome has a high recombination rate, although it also follows the overall rule of meiosis in eukaryotes with an average of approximately one crossover per chromosome per meiosis. GC-rich repetitive motifs identified in the hotspot sequences may play a role in the high recombination rate observed. The lack of recombination activity in centromeric regions is consistent with the observations of reduced recombination near the centromeres of other organisms.     

Sex and recombination purge the genome of deleterious alleles: An Individual Based Modeling Approach

The prevalence of sexual reproduction in most animal species despite its considerable costs such as useless males, energy spent on mating, the cost of meiosis and genome dilution remains a puzzle in evolutionary theory. One prominent single factor attempt to solve this persistent puzzle is the claim that sexual reproduction is instrumental in eliminating deleterious alleles from the species genome by the mechanism of recombination. There are three major versions of the deleterious allele hypothesis: First, the mutational deterministic hypothesis (MDH), which rests on the assumption of negative epistasis, predicts that recombination will help to purge the species genome of deleterious alleles by breaking apart linkages between these alleles. The assumption is that the joint negative effects of linked deleterious alleles is sometimes greater than the effects of the alleles considered separately. Second, there is the hypothesis that sexual reproduction speeds up purifying (negative) selection, which purges the genome of deleterious alleles. Alleles that are less deleterious than the wild type are naturally selected. These alleles, attained via recombination, are sometimes ‘leaky’ mutations giving rise to reduced functionality of attendant proteins. This hypothesis does not necessarily rest on the assumption of negative epistasis, which some argue is relatively rare in nature (Kouyos, Silander and Bonhoeffer (2012)) and which arguably could be seen as a virtue of the purifying selection hypothesis vs. the MDH. Third, Muller's ratchet hypothesis predicts that recombination will help to prevent the buildup of deleterious mutations by the mechanism of recombination. In this study, we focus primarily on testing the purifying selection hypothesis. We performed an individual-based model computer simulation using the program EcoSim to test this hypothesis. The experimental runs for sexual reproduction, asexual reproduction and facultative reproduction involved introducing a deleterious allele into the genome, which exacts an intermediate-level energy penalty on individuals. It was found that whereas on average, deleteriousness consistently declined over 18,000 time-steps due to recombination in sexual reproduction, deleteriousness did not decline for asexual and facultative runs. These results corroborate the hypothesis that recombination due to sexual reproduction helps to eliminate deleterious alleles from the genome through the selection of reduced function mutations.     

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.     

Context-dependent mutation rates may cause spurious signatures of a fixation bias favoring higher GC-content in humans

Understanding the proximate and ultimate causes underlying the evolution of nucleotide composition in mammalian genomes is of fundamental interest to the study of molecular evolution. Comparative genomics studies have revealed that many more substitutions occur from G and C nucleotides to A and T nucleotides than the reverse, suggesting that mammalian genomes are not at equilibrium for base composition. Analysis of human polymorphism data suggests that mutations that increase GC-content tend to be at much higher frequencies than those that decrease or preserve GC-content when the ancestral allele is inferred via parsimony using the chimpanzee genome. These observations have been interpreted as evidence for a fixation bias in favor of G and C alleles due to either positive natural selection or biased gene conversion. Here, we test the robustness of this interpretation to violations of the parsimony assumption using a data set of 21,488 noncoding single nucleotide polymorphisms (SNPs) discovered by the National Institute of Environmental Health Sciences (NIEHS) SNPs project via direct resequencing of n = 95 individuals. Applying standard nonparametric and parametric population genetic approaches, we replicate the signatures of a fixation bias in favor of G and C alleles when the ancestral base is assumed to be the base found in the chimpanzee outgroup. However, upon taking into account the probability of misidentifying the ancestral state of each SNP using a context-dependent mutation model, the corrected distribution of SNP frequencies for GC-content increasing SNPs are nearly indistinguishable from the patterns observed for other types of mutations, suggesting that the signature of fixation bias is a spurious artifact of the parsimony assumption.     

Suppression of intrachromosomal gene conversion in mammalian cells by small degrees of sequence divergence

Pairs of closely linked defective herpes simplex virus (HSV) thymidine kinase (tk) gene sequences exhibiting various nucleotide heterologies were introduced into the genome of mouse Ltk- cells. Recombination events were recovered by selecting for the correction of a 16-bp insertion mutation in one of the tk sequences. We had previously shown that when two tk sequences shared a region of 232 bp of homology, interruption of the homology by two single nucleotide heterologies placed 19 bp apart reduced recombination nearly 20-fold. We now report that either one of the nucleotide heterologies alone reduces recombination only about 2.5-fold, indicating that the original pair of single nucleotide heterologies acted synergistically to inhibit recombination. We tested a variety of pairs of single nucleotide heterologies and determined that they reduced recombination from 7- to 175-fold. Substrates potentially leading to G-G or C-C mispairs in presumptive heteroduplex DNA (hDNA) intermediates displayed a particularly low rate of recombination. Additional experiments suggested that increased sequence divergence causes a shortening of gene conversion tracts. Collectively, our results suggest that suppression of recombination between diverged sequences is mediated via processing of a mispaired hDNA intermediate.     

The Red Queen Model of Recombination Hotspots Evolution in the Light of Archaic and Modern Human Genomes

Recombination is an essential process in eukaryotes, which increases diversity by disrupting genetic linkage between loci and ensures the proper segregation of chromosomes during meiosis. In the human genome, recombination events are clustered in hotspots, whose location is determined by the PRDM9 protein. There is evidence that the location of hotspots evolves rapidly, as a consequence of changes in PRDM9 DNA-binding domain. However, the reasons for these changes and the rate at which they occur are not known. In this study, we investigated the evolution of human hotspot loci and of PRDM9 target motifs, both in modern and archaic human lineages (Denisovan) to quantify the dynamic of hotspot turnover during the recent period of human evolution. We show that present-day human hotspots are young: they have been active only during the last 10% of the time since the divergence from chimpanzee, starting to be operating shortly before the split between Denisovans and modern humans. Surprisingly, however, our analyses indicate that Denisovan recombination hotspots did not overlap with modern human ones, despite sharing similar PRDM9 target motifs. We further show that high-affinity PRDM9 target motifs are subject to a strong self-destructive drive, known as biased gene conversion (BGC), which should lead to the loss of the majority of them in the next 3 MYR. This depletion of PRDM9 genomic targets is expected to decrease fitness, and thereby to favor new PRDM9 alleles binding different motifs. Our refined estimates of the age and life expectancy of human hotspots provide empirical evidence in support of the Red Queen hypothesis of recombination hotspots evolution.     

Extreme Recombination Frequencies Shape Genome Variation and Evolution in the Honeybee,Apis mellifera

Meiotic recombination is a fundamental cellular process, with important consequences for evolution and genome integrity. However, we know little about how recombination rates vary across the genomes of most species and the molecular and evolutionary determinants of this variation. The honeybee, Apis mellifera, has extremely high rates of meiotic recombination, although the evolutionary causes and consequences of this are unclear. Here we use patterns of linkage disequilibrium in whole genome resequencing data from 30 diploid honeybees to construct a fine-scale map of rates of crossing over in the genome. We find that, in contrast to vertebrate genomes, the recombination landscape is not strongly punctate. Crossover rates strongly correlate with levels of genetic variation, but not divergence, which indicates a pervasive impact of selection on the genome. Germ-line methylated genes have reduced crossover rate, which could indicate a role of methylation in suppressing recombination. Controlling for the effects of methylation, we do not infer a strong association between gene expression patterns and recombination. The site frequency spectrum is strongly skewed from neutral expectations in honeybees: rare variants are dominated by AT-biased mutations, whereas GC-biased mutations are found at higher frequencies, indicative of a major influence of GC-biased gene conversion (gBGC), which we infer to generate an allele fixation bias 5 – 50 times the genomic average estimated in humans. We uncover further evidence that this repair bias specifically affects transitions and favours fixation of CpG sites. Recombination, via gBGC, therefore appears to have profound consequences on genome evolution in honeybees and interferes with the process of natural selection. These findings have important implications for our understanding of the forces driving molecular evolution.     

Playing hide and seek with mammalian meiotic crossover hotspots

Crossovers (COs) are essential for meiosis and contribute to genome diversity by promoting the reassociation of alleles, and thus improve the efficiency of selection. COs are not randomly distributed but are found at specific regions, or CO hotspots. Recent results have revealed the historical recombination rates and the distribution of hotspots across the human genome. Surprisingly, CO hotspots are highly dynamic, as shown by differences in activity between individuals, populations and closely related species. We propose a role for DNA methylation in preventing the formation of COs, a regulation that might explain, in part, the correlation between recombination rates and GC content in mammals.