Deleterious Passengers in Adapting Populations

Most new mutations are deleterious and are eventually eliminated by natural selection. But in an adapting population, the rapid amplification of beneficial mutations can hinder the removal of deleterious variants in nearby regions of the genome, altering the patterns of sequence evolution. Here, we analyze the interactions between beneficial “driver” mutations and linked deleterious “passengers” during the course of adaptation. We derive analytical expressions for the substitution rate of a deleterious mutation as a function of its fitness cost, as well as the reduction in the beneficial substitution rate due to the genetic load of the passengers. We find that the fate of each deleterious mutation varies dramatically with the rate and spectrum of beneficial mutations and the deleterious substitution rate depends nonmonotonically on the population size and the rate of adaptation. By quantifying this dependence, our results allow us to estimate which deleterious mutations will be likely to fix and how many of these mutations must arise before the progress of adaptation is significantly reduced.     

Human DNA sequence variation in a 6.6-kb region containing the melanocortin 1 receptor promoter.

An approximately 6.6-kb region located upstream from the melanocortin 1 receptor (MC1R) gene and containing its promoter was sequenced in 54 humans (18 Africans, 18 Asians, and 18 Europeans) and in one chimpanzee, gorilla, and orangutan. Seventy-six polymorphic sites were found among the human sequences and the average nucleotide diversity (pi) was 0.141%, one of the highest among all studies of nuclear sequence variation in humans. Opposite to the pattern observed in the MC1R coding region, in the present region pi is highest in Africans (0.136%) compared to Asians (0.116%) and Europeans (0.122%). The distributions of pi, theta, and Fu and Li's F-statistic are nonuniform along the sequence and among continents. The pattern of genetic variation is consistent with a population expansion in Africans. We also suggest a possible phase of population size reduction in non-Africans and purifying selection acting in the middle subregion and parts of the 5' subregion in Africans. We hypothesize diversifying selection acting on some sites in the 5' and 3' subregions or in the MC1R coding region in Asians and Europeans, though we cannot reject the possibility of relaxation of functional constraints in the MC1R gene in Asians and Europeans. The mutation rate in the sequenced region is 1.65 x 10(-9) per site per year. The age of the most recent common ancestor for this region is similar to that for the other long noncoding regions studied to date, providing evidence for ancient gene genealogies. Our population screening and phylogenetic footprinting suggest potentially important sites for the MC1R promoter function.     

Inferring microevolutionary patterns from allele-size frequency distributions of minisatellite loci: a worldwide study of the APOB 3' hypervariable region polymorphism

The availability of numerous population and molecular data makes the apolipoprotein B 3' hypervariable region (APOB 3' HVR) polymorphism ideal for a pilot study of the relationships between the allele-size frequency distributions (referred to as allele-size distributions) of minisatellite loci and the microevolutionary processes underlying their present-day polymorphism in human populations. In this paper, we present a worldwide APOB 3' HVR study, based on published and unpublished data, which refers to 36 populations. We systematically compare APOB 3' HVR within-group diversity (in terms of heterozygosity, number of alleles, and allele-size variance) in numerous human populations, including African, European, Asian, Amerindian, Australomelanesian, and Polynesian groups. Overall, our analyses indicate a greater APOB 3' HVR diversity in Africans than non-Africans. Then, we compare APOB 3' HVR allele-size distributions. The APOB 3' HVR allele-size distribution is found to be quasi-unimodal in Africans and bimodal or nonunimodal in non-African populations. The analysis of the distribution of pairwise comparisons suggests that Africans expanded earlier and/or that their ancestral population was larger than other continental groups. As a final step, we examine APOB 3' HVR interpopulational relationships by using three genetic distances. The F(ST) genetic distance, which assumes genetic drift as being the agent that differentiates populations, provides results that are more congruent with established anthropological knowledge than mutation-based distances (D(SW) and R(ST)). We hypothesize that the ancestral population was characterized by a high heterozygosity, an extended range of allele size, and a quasi-unimodal allele-size distribution centered on allele *37, features persisting in examined African populations. Sampling processes during "out-of-Africa" migrations would be responsible for the decrease in APOB 3' HVR gene diversity and the nonunimodal allele-size distribution observed in non-Africans. Some possible confounding factors are discussed and a prospect of how the hypothesis could be refined and tested is given.     

Out of Africa? What do genes tell us?

Genetic diversity patterns in nuclear versus mitochondrial systems and in low versus high mutation rate systems do not support the hypothesis of a recent African origin for all of humanity following a split between Africans and non-Africans 100,000 years ago, nor do genetic distance data. Geographical analyses of nuclear and mitochondrial gene trees do not support the hypothesis of a recent global replacement of humans coming out of Africa, although a local replacement event in Europe is indicated by these analyses and recent studies on Neandertal DNA. The gene tree analyses instead indicate that genetic interchanges have ensured that all of humanity has evolved as a single evolutionary lineage with no major splits.     

Ancient structure in Africa unlikely to explain neanderthal and non-african genetic similarity

Neanderthals have been shown to share more genetic variants with present-day non-Africans than Africans. Recent admixture between Neanderthals and modern humans outside of Africa was proposed as the most parsimonious explanation for this observation. However, the hypothesis of ancient population structure within Africa could not be ruled out as an alternative explanation. We use simulations to test whether the site frequency spectrum, conditioned on a derived Neanderthal and an ancestral Yoruba (African) nucleotide (the doubly conditioned site frequency spectrum [dcfs]), can distinguish between models that assume recent admixture or ancient population structure. We compare the simulations to the dcfs calculated from data taken from populations of European, Chinese, and Japanese descent in the Complete Genomics Diversity Panel. Simulations under a variety of plausible demographic parameters were used to examine the shape of the dcfs for both models. The observed shape of the dcfs cannot be explained by any set of parameter values used in the simulations of the ancient structure model. The dcfs simulations for the recent admixture model provide a good fit to the observed dcfs for non-Africans, thereby supporting the hypothesis that recent admixture with Neanderthals accounts for the greater similarity of Neanderthals to non-Africans than Africans.     

Somatic mtDNA Mutation Spectra in the Aging Human Putamen

The accumulation of heteroplasmic mitochondrial DNA (mtDNA) deletions and single nucleotide variants (SNVs) is a well-accepted facet of the biology of aging, yet comprehensive mutation spectra have not been described. To address this, we have used next generation sequencing of mtDNA-enriched libraries (Mito-Seq) to investigate mtDNA mutation spectra of putamen from young and aged donors. Frequencies of the “common” deletion and other “major arc” deletions were significantly increased in the aged cohort with the fold increase in the frequency of the common deletion exceeding that of major arc deletions. SNVs also increased with age with the highest rate of accumulation in the non-coding control region which contains elements necessary for translation and replication. Examination of predicted amino acid changes revealed a skew towards pathogenic SNVs in the coding region driven by mutation bias. Levels of the pathogenic m.3243A>G tRNA mutation were also found to increase with age. Novel multimeric tandem duplications that resemble murine control region multimers and yeast ρ⁻ mtDNAs, were identified in both young and aged specimens. Clonal ∼50 bp deletions in the control region were found at high frequencies in aged specimens. Our results reveal the complex manner in which the mitochondrial genome alters with age and provides a foundation for studies of other tissues and disease states.     

X-Linked MTMR8 Diversity and Evolutionary History of Sub-Saharan Populations

The genetic diversity within an 11 kb segment of the MTMR8 gene in a sample of 111 sub-Saharan and 49 non-African X chromosomes was investigated to assess the early evolutionary history of sub-Saharan Africans and the out-of-Africa expansion. The analyses revealed a complex genetic structure of the Africans that contributed to the emergence of modern humans. We observed partitioning of two thirds of old lineages among southern, west/central and east African populations indicating ancient population stratification predating the out of Africa migration. Age estimates of these lineages, older than coalescence times of uniparentally inherited markers, raise the question whether contemporary humans originated from a single population or as an amalgamation of different populations separated by years of independent evolution, thus suggesting a greater antiquity of our species than generally assumed. While the oldest sub-Saharan lineages, ∼500 thousand years, are found among Khoe-San from southern-Africa, a distinct haplotype found among Biaka is likely due to admixture from an even older population. An East African population that gave rise to non-Africans underwent a selective sweep affecting the subcentromeric region where MTMR8 is located. This and similar sweeps in four other regions of the X chromosome, documented in the literature, effectively reduced genetic diversity of non-African chromosomes and therefore may have exacerbated the effect of the demographic bottleneck usually ascribed to the out of Africa migration. Our data is suggestive, however, that a bottleneck, occurred in Africa before range expansion.     

Searching for evidence of positive selection in the human genome using patterns of microsatellite variability

Both natural selection and nonequilibrium population-level processes can lead to a skew in the frequency distribution of polymorphisms. Population-level processes are expected to affect all loci in a roughly equal fashion, whereas selection will affect only some regions of the genome. We conducted a sliding-window analysis of the frequency distribution of microsatellite polymorphisms across the human genome to identify regions that may be under positive selection. The analysis was based on a published data set of 5,257 mapped microsatellites in individuals of European ancestry. Observed and expected numbers of alleles were compared under a stepwise mutation model (SMM) using analytical formulae. Observed and expected heterozygosities were compared under a SMM using coalescent simulations. The two sets of analyses gave similar results. Approximately one-fourth of all loci showed a significant deficit of heterozygosity, consistent with a recent population expansion. Forty-three windows were identified with extreme skews in the frequency distribution of polymorphisms (in the direction of a deficit of heterozygosity, given the number of alleles). If these extreme windows are tracking selection at linked sites, theory predicts that they should be more common in regions of the genome with less recombination. We tested this prediction by comparing recombination rates in these extreme windows and in other regions of the genome and found that extreme windows had a significantly lower recombination rate than the genomic average. The proportion of extreme windows was significantly higher on the X chromosome than on the autosomes. Moreover, all the windows with extreme skews on the X chromosome were found in two clusters near the centromere; both these clusters exhibit markedly reduced recombination rates. These analyses point to regions of the genome that may recently have been subject to positive selection. These results also suggest that the effects of positive selection may be more pronounced on the X chromosome than on the autosomes in humans.     

North African Populations Carry the Signature of Admixture with Neandertals

One of the main findings derived from the analysis of the Neandertal genome was the evidence for admixture between Neandertals and non-African modern humans. An alternative scenario is that the ancestral population of non-Africans was closer to Neandertals than to Africans because of ancient population substructure. Thus, the study of North African populations is crucial for testing both hypotheses. We analyzed a total of 780,000 SNPs in 125 individuals representing seven different North African locations and searched for their ancestral/derived state in comparison to different human populations and Neandertals. We found that North African populations have a significant excess of derived alleles shared with Neandertals, when compared to sub-Saharan Africans. This excess is similar to that found in non-African humans, a fact that can be interpreted as a sign of Neandertal admixture. Furthermore, the Neandertal''s genetic signal is higher in populations with a local, pre-Neolithic North African ancestry. Therefore, the detected ancient admixture is not due to recent Near Eastern or European migrations. Sub-Saharan populations are the only ones not affected by the admixture event with Neandertals.     

Extremely Rare Interbreeding Events Can Explain Neanderthal DNA in Living Humans

Considering the recent experimental discovery of Green et al that present-day non-Africans have 1 to of their nuclear DNA of Neanderthal origin, we propose here a model which is able to quantify the genetic interbreeding between two subpopulations with equal fitness, living in the same geographic region. The model consists of a solvable system of deterministic ordinary differential equations containing as a stochastic ingredient a realization of the neutral Wright-Fisher process. By simulating the stochastic part of the model we are able to apply it to the interbreeding ofthe African ancestors of Eurasians and Middle Eastern Neanderthal subpopulations and estimate the only parameter of the model, which is the number of individuals per generation exchanged between subpopulations. Our results indicate that the amount of Neanderthal DNA in living non-Africans can be explained with maximum probability by the exchange of a single pair of individuals between the subpopulations at each 77 generations, but larger exchange frequencies are also allowed with sizeable probability. The results are compatible with a long coexistence time of 130,000 years, a total interbreeding population of order individuals, and with all living humans being descendants of Africans both for mitochondrial DNA and Y chromosome.     

The Date of Interbreeding between Neandertals and Modern Humans

Comparisons of DNA sequences between Neandertals and present-day humans have shown that Neandertals share more genetic variants with non-Africans than with Africans. This could be due to interbreeding between Neandertals and modern humans when the two groups met subsequent to the emergence of modern humans outside Africa. However, it could also be due to population structure that antedates the origin of Neandertal ancestors in Africa. We measure the extent of linkage disequilibrium (LD) in the genomes of present-day Europeans and find that the last gene flow from Neandertals (or their relatives) into Europeans likely occurred 37,000–86,000 years before the present (BP), and most likely 47,000–65,000 years ago. This supports the recent interbreeding hypothesis and suggests that interbreeding may have occurred when modern humans carrying Upper Paleolithic technologies encountered Neandertals as they expanded out of Africa.     

Comparative nuclear and mitochondrial genome diversity in humans and chimpanzees

Restriction mapping and sequencing have shown that humans have substantially lower levels of mitochondrial genome diversity (d) than chimpanzees. In contrast, humans have substantially higher levels of heterozygosity (H) at protein-coding loci, suggesting a higher level of diversity in the nuclear genome. To investigate the discrepancy further, we sequenced a segment of the mitochondrial genome control region (CR) from 49 chimpanzees. The majority of these were from the Pan troglodytes versus subspecies, which was underrepresented in previous studies. We also estimated the average heterozygosity at 60 short tandem repeat (STR) loci in both species. For a total sample of 115 chimpanzees, d = 0.075 +/0 0.037, compared to 0.020 +/- 0.011 for a sample of 1,554 humans. The heterozygosity of human STR loci is significantly higher than that of chimpanzees. Thus, the higher level of nuclear genome diversity relative to mitochondrial genome diversity in humans is not restricted to protein-coding loci. It seems that humans, not chimpanzees, have an unusual d/H ratio, since the ratio in chimpanzees is similar to that in other catarrhines. This discrepancy in the relative levels of nuclear and mitochondrial genome diversity in the two species cannot be explained by differences in mutation rate. However, it may result from a combination of factors such as a difference in the extent of sex ratio disparity, the greater effect of population subdivision on mitochondrial than on nuclear genome diversity, a difference in the relative levels of male and female migration among subpopulations, diversifying selection acting to increase variation in the nuclear genome, and/or directional selection acting to reduce variation in the mitochondrial genome.      

The rate of adaptation in asexuals

I study the population genetics of adaptation in asexuals. I show that the rate of adaptive substitution in an asexual species or nonrecombining chromosome region is a bell-shaped function of the mutation rate: at some point, increasing the mutation rate decreases the rate of substitution. Curiously, the mutation rate that maximizes the rate of adaptation depends solely on the strength of selection against deleterious mutations. In particular, adaptation is fastest when the genomic rate of mutation, U, equals the harmonic mean of selection coefficients against deleterious mutations, where we assume that selection for favorable alleles is milder than that against deleterious ones. This simple result is independent of the shape of the distribution of effects among favorable and deleterious mutations, population size, and the action of clonal interference. In the course of this work, I derive an approximation to the probability of fixation of a favorable mutation in an asexual genome or nonrecombining chromosome region in which both favorable and deleterious mutations occur.     

Genetic Variation and Random Drift in Autotetraploid Populations

The rate of decay of genetic variation is determined for randomly mating autotetraploid populations of finite size, and the equilibrium homozygosity under mutation and random drift is calculated. It is shown that heterozygosity is lost at a slower rate than in diploid populations, and that the equilibrium heterozygosity with mutation and random drift is higher than for diploids. Outcrossing populations as well as populations that randomly self are analyzed. A method of comparing genetic variation between autotetraploid and diploid populations is proposed. Our treatment suggests that the ``gametic homozygosity' provides a unified approach for comparing genotypes within a population as well as comparing genetic variation between populations with different levels of ploidy.     

The Role of Viral Population Diversity in Adaptation of Bovine Coronavirus to New Host Environments

The high mutation rate of RNA viruses enables a diverse genetic population of viral genotypes to exist within a single infected host. In-host genetic diversity could better position the virus population to respond and adapt to a diverse array of selective pressures such as host-switching events. Multiple new coronaviruses, including SARS, have been identified in human samples just within the last ten years, demonstrating the potential of coronaviruses as emergent human pathogens. Deep sequencing was used to characterize genomic changes in coronavirus quasispecies during simulated host-switching. Three bovine nasal samples infected with bovine coronavirus were used to infect human and bovine macrophage and lung cell lines. The virus reproduced relatively well in macrophages, but the lung cell lines were not infected efficiently enough to allow passage of non lab-adapted samples. Approximately 12 kb of the genome was amplified before and after passage and sequenced at average coverages of nearly 950×(454 sequencing) and 38,000×(Illumina). The consensus sequence of many of the passaged samples had a 12 nucleotide insert in the consensus sequence of the spike gene, and multiple point mutations were associated with the presence of the insert. Deep sequencing revealed that the insert was present but very rare in the unpassaged samples and could quickly shift to dominate the population when placed in a different environment. The insert coded for three arginine residues, occurred in a region associated with fusion entry into host cells, and may allow infection of new cell types via heparin sulfate binding. Analysis of the deep sequencing data indicated that two distinct genotypes circulated at different frequency levels in each sample, and support the hypothesis that the mutations present in passaged strains were “selected” from a pre-existing pool rather than through de novo mutation and subsequent population fixation.     

The Oxytricha trifallax Macronuclear Genome: A Complex Eukaryotic Genome with 16,000 Tiny Chromosomes

The macronuclear genome of the ciliate Oxytricha trifallax displays an extreme and unique eukaryotic genome architecture with extensive genomic variation. During sexual genome development, the expressed, somatic macronuclear genome is whittled down to the genic portion of a small fraction (∼5%) of its precursor “silent” germline micronuclear genome by a process of “unscrambling” and fragmentation. The tiny macronuclear “nanochromosomes” typically encode single, protein-coding genes (a small portion, 10%, encode 2–8 genes), have minimal noncoding regions, and are differentially amplified to an average of ∼2,000 copies. We report the high-quality genome assembly of ∼16,000 complete nanochromosomes (∼50 Mb haploid genome size) that vary from 469 bp to 66 kb long (mean ∼3.2 kb) and encode ∼18,500 genes. Alternative DNA fragmentation processes ∼10% of the nanochromosomes into multiple isoforms that usually encode complete genes. Nucleotide diversity in the macronucleus is very high (SNP heterozygosity is ∼4.0%), suggesting that Oxytricha trifallax may have one of the largest known effective population sizes of eukaryotes. Comparison to other ciliates with nonscrambled genomes and long macronuclear chromosomes (on the order of 100 kb) suggests several candidate proteins that could be involved in genome rearrangement, including domesticated MULE and IS1595-like DDE transposases. The assembly of the highly fragmented Oxytricha macronuclear genome is the first completed genome with such an unusual architecture. This genome sequence provides tantalizing glimpses into novel molecular biology and evolution. For example, Oxytricha maintains tens of millions of telomeres per cell and has also evolved an intriguing expansion of telomere end-binding proteins. In conjunction with the micronuclear genome in progress, the O. trifallax macronuclear genome will provide an invaluable resource for investigating programmed genome rearrangements, complementing studies of rearrangements arising during evolution and disease.     

A High Resolution Genome-Wide Scan for Significant Selective Sweeps: An Application to Pooled Sequence Data in Laying Chickens

In most studies aimed at localizing footprints of past selection, outliers at tails of the empirical distribution of a given test statistic are assumed to reflect locus-specific selective forces. Significance cutoffs are subjectively determined, rather than being related to a clear set of hypotheses. Here, we define an empirical p-value for the summary statistic by means of a permutation method that uses the observed SNP structure in the real data. To illustrate the methodology, we applied our approach to a panel of 2.9 million autosomal SNPs identified from re-sequencing a pool of 15 individuals from a brown egg layer line. We scanned the genome for local reductions in heterozygosity, suggestive of selective sweeps. We also employed a modified sliding window approach that accounts for gaps in the sequence and increases scanning resolution by moving the overlapping windows by steps of one SNP only, and suggest to call this a “creeping window” strategy. The approach confirmed selective sweeps in the region of previously described candidate genes, i.e. TSHR, PRL, PRLHR, INSR, LEPR, IGF1, and NRAMP1 when used as positive controls. The genome scan revealed 82 distinct regions with strong evidence of selection (genome-wide p-value<0.001), including genes known to be associated with eggshell structure and immune system such as CALB1 and GAL cluster, respectively. A substantial proportion of signals was found in poor gene content regions including the most extreme signal on chromosome 1. The observation of multiple signals in a highly selected layer line of chicken is consistent with the hypothesis that egg production is a complex trait controlled by many genes.     

The early life social environment and DNA methylation: DNA methylation mediating the long-term impact of social environments early in life

Although epidemiological data provides evidence that there is an interaction between genetics (nature) and the social and physical environments (nurture) in human development; the main open question remains the mechanism. The pattern of distribution of methyl groups in DNA is different from cell-type to cell type and is conferring cell specific identity on DNA during cellular differentiation and organogenesis. This is an innate and highly programmed process. However, recent data suggests that DNA methylation is not only involved in cellular differentiation but that it is also involved in modulation of genome function in response to signals from the physical, biological and social environments. We propose that modulation of DNA methylation in response to environmental cues early in life serves as a mechanism of life-long genome “adaptation” that molecularly embeds the early experiences of a child (“nurture”) in the genome (“nature”). There is an emerging line of data supporting this hypothesis in rodents, non-human primates and humans that will be reviewed here. However, several critical questions remain including the identification of mechanisms that transmit the signals from the social environment to the DNA methylation/demethylation enzymes.Key words: DNA methylation, psychiatry, development, epidemiology, environment     

Mutations of different molecular origins exhibit contrasting patterns of regional substitution rate variation

Transitions at CpG dinucleotides, referred to as “CpG substitutions”, are a major mutational input into vertebrate genomes and a leading cause of human genetic disease. The prevalence of CpG substitutions is due to their mutational origin, which is dependent on DNA methylation. In comparison, other single nucleotide substitutions (for example those occurring at GpC dinucleotides) mainly arise from errors during DNA replication. Here we analyzed high quality BAC-based data from human, chimpanzee, and baboon to investigate regional variation of CpG substitution rates.

We show that CpG substitutions occur approximately 15 times more frequently than other single nucleotide substitutions in primate genomes, and that they exhibit substantial regional variation. Patterns of CpG rate variation are consistent with differences in methylation level and susceptibility to subsequent deamination. In particular, we propose a “distance-decaying” hypothesis, positing that due to the molecular mechanism of a CpG substitution, rates are correlated with the stability of double-stranded DNA surrounding each CpG dinucleotide, and the effect of local DNA stability may decrease with distance from the CpG dinucleotide.

Consistent with our “distance-decaying” hypothesis, rates of CpG substitution are strongly (negatively) correlated with regional G+C content. The influence of G+C content decays as the distance from the target CpG site increases. We estimate that the influence of local G+C content extends up to 1,500~2,000 bps centered on each CpG site. We also show that the distance-decaying relationship persisted when we controlled for the effect of long-range homogeneity of nucleotide composition. GpC sites, in contrast, do not exhibit such “distance-decaying” relationship. Our results highlight an example of the distinctive properties of methylation-dependent substitutions versus substitutions mostly arising from errors during DNA replication. Furthermore, the negative relationship between G+C content and CpG rates may provide an explanation for the observation that GC-rich SINEs show lower CpG rates than other repetitive elements.

     

Reverse Chemical Ecology Suggests Putative Primate Pheromones

Pheromonal communication is widespread among living organisms, but in apes and particularly in humans there is currently no strong evidence for such phenomenon. Among primates, lemurs use pheromones to communicate within members of the same species, whereas in some monkeys such capabilities seem to be lost. Chemical communication in humans appears to be impaired by the lack or malfunctioning of biochemical tools and anatomical structures mediating detection of pheromones. Here, we report on a pheromone-carrier protein (SAL) adopting a “reverse chemical ecology” approach to get insights on the structures of potential pheromones in a representative species of lemurs (Microcebus murinus) known to use pheromones, Old-World monkeys (Cercocebus atys) for which chemical communication has been observed, and humans (Homo sapiens), where pheromones and chemical communication are still questioned. We have expressed the SAL orthologous proteins of these primate species, after reconstructing the gene encoding the human SAL, which is disrupted due to a single base mutation preventing its translation into RNA. Ligand-binding experiments with the recombinant SALs revealed macrocyclic ketones and lactones as the best ligands for all three proteins, suggesting cyclopentadecanone, pentadecanolide, and closely related compounds as the best candidates for potential pheromones. Such hypothesis agrees with the presence of a chemical very similar to hexadecanolide in the gland secretions of Mandrillus sphinx, a species closely related to C. atys. Our results indicate that the function of this carrier protein has not changed much during evolution from lemurs to humans, although its physiological role has been certainly impaired in humans.