Adaptive Evolution in Linked Genomes

Adaptive evolution is governed by various forces: Mutations occur randomly in the genome and generate variability in the individuals’reproductive success; natural selection shifts this variability in the population towards individuals with high fitness; genetic drift introduces random fluctuations in the number of offspring of an individual and affects weakly selected or neutral mutations. On top of these, genetic linkage can be an important evolutionary force. Linkage generates interference interactions, by which simultaneously occurring mutations affect each other’s chance of fixation. Here, we develop a comprehensive model of adaptive evolution in linked genomes, which integrates interference interactions between multiple beneficial and deleterious mutations into a unified framework. By an approximate analytical solution, we predict the fixation rates of these mutations, as well as the probabilities of beneficial and deleterious alleles at fixed genomic sites. We find that interference interactions generate a regime of emergent neutrality: all genomic sites with selection coefficients smaller in magnitude than a characteristic threshold have nearly random fixed alleles, and both beneficial and deleterious mutations at these sites have nearly neutral fixation rates. We show that this dynamics limits not only the speed of adaptation, but also a population’s degree of adaptation in its current environment. We apply the model to different scenarios: stationary adaptation in a time-dependent environment, and approach to equilibrium in a fixed environment. In both cases, the analytical predictions are in good agreement with numerical simulations. Our results suggest that interference can severely compromise biological functions in an adapting population, which sets viability limits on adaptive evolution under linkage. We furthermore develop a likelihood-based inference method for genomic data, which explicitly takes into account genetic linkage. Tests with simulated datasets show that our method correctly predicts the amount of positive selection in linked sequence. In contrast, many existing tests falsely interpret traces from linkage as spurious positive selection. We apply our method to fruit fly genome data (Drosophila melanogaster), and find that a substantial fraction of sequence differences between two related fly species is in fact caused by linkage instead of natural selection

Rate and cost of adaptation in the Drosophila Genome

Recent studies have consistently inferred high rates of adaptive molecular evolution between Drosophila species. At the same time, the Drosophila genome evolves under different rates of recombination, which results in partial genetic linkage between alleles at neighboring genomic loci. Here we analyze how linkage correlations affect adaptive evolution. We develop a new inference method for adaptation that takes into account the effect on an allele at a focal site caused by neighboring deleterious alleles (background selection) and by neighboring adaptive substitutions (hitchhiking). Using complete genome sequence data and fine-scale recombination maps, we infer a highly heterogeneous scenario of adaptation in Drosophila. In high-recombining regions, about 50% of all amino acid substitutions are adaptive, together with about 20% of all substitutions in proximal intergenic regions. In low-recombining regions, only a small fraction of the amino acid substitutions are adaptive, while hitchhiking accounts for the majority of these changes. Hitchhiking of deleterious alleles generates a substantial collateral cost of adaptation, leading to a fitness decline of about 30/2N per gene and per million years in the lowest-recombining regions. Our results show how recombination shapes rate and efficacy of the adaptive dynamics in eukaryotic genomes

The asexual genome of Drosophila

The rate of recombination affects the mode of molecular evolution. In high-recombining sequence, the targets of selection are individual genetic loci; under low recombination, selection collectively acts on large, genetically linked genomic segments. Selection under linkage can induce clonal interference, a specific mode of evolution by competition of genetic clades within a population. This mode is well known in asexually evolving microbes, but has not been traced systematically in an obligate sexual organism. Here we show that the Drosophila genome is partitioned into two modes of evolution: a local interference regime with limited effects of genetic linkage, and an interference condensate with clonal competition. We map these modes by differences in mutation frequency spectra, and we show that the transition between them occurs at a threshold recombination rate that is predictable from genomic summary statistics. We find the interference condensate in segments of low-recombining sequence that are located primarily in chromosomal regions flanking the centromeres and cover about 20% of the Drosophila genome. Condensate regions have characteristics of asexual evolution that impact gene function: the efficacy of selection and the speed of evolution are lower and the genetic load is higher than in regions of local interference. Our results suggest that multicellular eukaryotes can harbor heterogeneous modes and tempi of evolution within one genome. We argue that this variation generates selection on genome architecture

High-coverage genome of the Tyrolean Iceman reveals unusually high Anatolian farmer ancestry

The Tyrolean Iceman is known as one of the oldest human glacier mummies, directly dated to 3350-3120 calibrated BCE. A previously published low-coverage genome provided novel insights into European prehistory, despite high present-day DNA contamination. Here, we generate a high-coverage genome with low contamination (15.3×) to gain further insights into the genetic history and phenotype of this individual. Contrary to previous studies, we found no detectable Steppe-related ancestry in the Iceman. Instead, he retained the highest Anatolian-farmer-related ancestry among contemporaneous European populations, indicating a rather isolated Alpine population with limited gene flow from hunter-gatherer-ancestry-related populations. Phenotypic analysis revealed that the Iceman likely had darker skin than present-day Europeans and carried risk alleles associated with male-pattern baldness, type 2 diabetes, and obesity-related metabolic syndrome. These results corroborate phenotypic observations of the preserved mummified body, such as high pigmentation of his skin and the absence of hair on his head

Genomic islands of speciation separate cichlid ecomorphs in an East African crater lake

The genomic causes and effects of divergent ecological selection during speciation are still poorly understood. Here, we report the discovery and detailed characterization of early-stage adaptive divergence of two cichlid fish ecomorphs in a small (700m diameter) isolated crater lake in Tanzania. The ecomorphs differ in depth preference, male breeding color, body shape, diet and trophic morphology. With whole genome sequences of 146 fish, we identify 98 clearly demarcated genomic ‘islands’ of high differentiation and demonstrate association of genotypes across these islands to divergent mate preferences. The islands contain candidate adaptive genes enriched for functions in sensory perception (including rhodopsin and other twilight vision associated genes), hormone signaling and morphogenesis. Our study suggests mechanisms and genomic regions that may play a role in the closely related mega-radiation of Lake Malawi.The work was funded by Royal Society-Leverhulme Trust Africa Awards AA100023 and AA130107 (M.J.G., B.P.N. and G.F.T.), a Wellcome Trust PhD studentship grant 097677/Z/11/Z (M.M.), Wellcome Trust grant WT098051 (S.S. and R.D.), Wellcome Trust and Cancer Research UK core support and a Wellcome Trust Senior Investigator Award (E.A.M.), a Leverhulme Trust Research Fellowship RF-2014-686 (M.J.G.), a University of Bristol Research Committee award (M.G.), a Bangor University Anniversary PhD studentship (to A.M.T.) and a Fisheries Society of the British Isles award (G.F.T.). Raw sequencing reads are in the SRA nucleotide archive: RAD sequencing (BioProject: PRJNA286304; accessions SAMN03768857 to SAMN03768912) and whole genome sequencing (BioProject PRJEB1254: sample accessions listed in Table S16). The RAD based phylogeny and alignments have been deposited in TreeBase (TB2:S18241). Whole genome variant calls in the VCF format, phylogenetic trees, and primer sequences for Sequenom genotyping are available from the Dryad Digital Repository (http://dx.doi.org/10.5061/dryad.770mc). RD declares his interests as a founder and non-executive director of Congenica Ltd., that he owns stock in Illumina from previous consulting, and is a scientific advisory board member of Dovetail Inc. We thank R. Schley for generating pharyngeal jaw data; S. Mzighani, J. Kihedu and staff of the Tanzanian Fisheries Research Institute for logistical support; A. Smith, H. Sungani, A. Shechonge, P. Parsons, J. Swanstrom, G. Cooke and J. Bridle for contributions to sampling and aquarium maintenance, the Sanger Institute sequencing core for DNA sequencing and Dr. H. Imai (Kyoto University) for the use of spectrometer in his laboratory.This is the author accepted manuscript. The final version is available from AAAS via http://dx.doi.org/10.1126/science.aac992

Intramural child burials in Iron Age Navarra: How ancient DNA can contribute to household archaeology

The transition from the Late Bronze to the Iron Age on the Iberian Peninsula saw a shift in mortuary customs from mainly inhumation to cremation of the deceased. The poor preservation characteristic of cremated skeletal remains has hindered molecular analyses (isotope analyses, ancient DNA) of the Iberian Final Bronze and Iron Age communities of Iberia. Incidentally, a limited number of young children, often newborns, were exempt from the predominant cremation ritual, in favour of intramural inhumations inside buildings at certain settlements. The discourse surrounding the mean- ing and interpretation of this particular burial rite has developed over a long time in Iberian archaeology but has always been hampered by the limited anthropological, archaeological, and molecular data from these intramural inhumations. Here, we study the genomes of 37 intramurally buried children found in three Early Iron Age settlements, dated between c. 800–450 BC. Population genetic analyses on the newly reported individuals extend our understanding of ancient Iberia by revealing previously unsampled genetic diversity as well as showing a lesser influence of Mediterranean ancestry than on previously published Iron Age individuals from northern Spain. We also provide insights into the sex and biological relatedness of the children, and in so doing, elucidate differ- ent aspects of the intramural burial ritual and building use in settlements. More broadly, the genetic data from these individuals fill an important gap in the archaeogenetic record of northern Spain and offer a unique opportunity to study the genetic makeup and population changes from the Bronze Age to Antiquity.This project has received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (grant agreement number 851511). It has also been supported by the research project »Convergence and interaction between complex Bronze Age societies« from the Academia program of the Institució Catalana de Recerca i Estudis Avançats (ICREA) of the Catalan Government and the Spanish Ministry for Science and Innovation (PID2020-112909GB-100)

Genomic evidence for the Pleistocene and recent population history of Native Americans

This is the author’s version of the work. It is posted here by permission of the AAAS for personal use, not for redistribution. The definitive version was published in Science on 2015 August 21; 349(6250), DOI: 10.1126/science.aab3884.How and when the Americas were populated remains contentious. Using ancient and modern genome-wide data, we find that the ancestors of all present-day Native Americans, including Athabascans and Amerindians, entered the Americas as a single migration wave from Siberia no earlier than 23 thousand years ago (KYA), and after no more than 8,000-year isolation period in Beringia. Following their arrival to the Americas, ancestral Native Americans diversified into two basal genetic branches around 13 KYA, one that is now dispersed across North and South America and the other is restricted to North America. Subsequent gene flow resulted in some Native Americans sharing ancestry with present-day East Asians (including Siberians) and, more distantly, Australo-Melanesians. Putative ‘Paleoamerican’ relict populations, including the historical Mexican Pericúes and South American Fuego-Patagonians, are not directly related to modern Australo-Melanesians as suggested by the Paleoamerican Model

Genomic islands of speciation separate cichlid ecomorphs in an East African crater lake.

The genomic causes and effects of divergent ecological selection during speciation are still poorly understood. Here we report the discovery and detailed characterization of early-stage adaptive divergence of two cichlid fish ecomorphs in a small (700 meters in diameter) isolated crater lake in Tanzania. The ecomorphs differ in depth preference, male breeding color, body shape, diet, and trophic morphology. With whole-genome sequences of 146 fish, we identified 98 clearly demarcated genomic "islands" of high differentiation and demonstrated the association of genotypes across these islands with divergent mate preferences. The islands contain candidate adaptive genes enriched for functions in sensory perception (including rhodopsin and other twilight-vision-associated genes), hormone signaling, and morphogenesis. Our study suggests mechanisms and genomic regions that may play a role in the closely related mega-radiation of Lake Malawi.The work was funded by Royal Society-Leverhulme Trust Africa Awards AA100023 and AA130107 (M.J.G., B.P.N. and G.F.T.), a Wellcome Trust PhD studentship grant 097677/Z/11/Z (M.M.), Wellcome Trust grant WT098051 (S.S. and R.D.), Wellcome Trust and Cancer Research UK core support and a Wellcome Trust Senior Investigator Award (E.A.M.), a Leverhulme Trust Research Fellowship RF-2014-686 (M.J.G.), a University of Bristol Research Committee award (M.G.), a Bangor University Anniversary PhD studentship (to A.M.T.) and a Fisheries Society of the British Isles award (G.F.T.). Raw sequencing reads are in the SRA nucleotide archive: RAD sequencing (BioProject: PRJNA286304; accessions SAMN03768857 to SAMN03768912) and whole genome sequencing (BioProject PRJEB1254: sample accessions listed in Table S16). The RAD based phylogeny and alignments have been deposited in TreeBase (TB2:S18241). Whole genome variant calls in the VCF format, phylogenetic trees, and primer sequences for Sequenom genotyping are available from the Dryad Digital Repository (http://dx.doi.org/10.5061/dryad.770mc). RD declares his interests as a founder and non-executive director of Congenica Ltd., that he owns stock in Illumina from previous consulting, and is a scientific advisory board member of Dovetail Inc. We thank R. Schley for generating pharyngeal jaw data; S. Mzighani, J. Kihedu and staff of the Tanzanian Fisheries Research Institute for logistical support; A. Smith, H. Sungani, A. Shechonge, P. Parsons, J. Swanstrom, G. Cooke and J. Bridle for contributions to sampling and aquarium maintenance, the Sanger Institute sequencing core for DNA sequencing and Dr. H. Imai (Kyoto University) for the use of spectrometer in his laboratory.This is the author accepted manuscript. The final version is available from AAAS via http://dx.doi.org/10.1126/science.aac992

Tracing the Route of Modern Humans out of Africa by Using 225 Human Genome Sequences from Ethiopians and Egyptians

The predominantly African origin of all modern human populations is well established, but the route taken out of Africa is still unclear. Two alternative routes, via Egypt and Sinai or across the Bab el Mandeb strait into Arabia, have traditionally been proposed as feasible gateways in light of geographic, paleoclimatic, archaeological, and genetic evidence. Distinguishing among these alternatives has been difficult. We generated 225 whole-genome sequences (225 at 8× depth, of which 8 were increased to 30×; Illumina HiSeq 2000) from six modern Northeast African populations (100 Egyptians and five Ethiopian populations each represented by 25 individuals). West Eurasian components were masked out, and the remaining African haplotypes were compared with a panel of sub-Saharan African and non-African genomes. We showed that masked Northeast African haplotypes overall were more similar to non-African haplotypes and more frequently present outside Africa than were any sets of haplotypes derived from a West African population. Furthermore, the masked Egyptian haplotypes showed these properties more markedly than the masked Ethiopian haplotypes, pointing to Egypt as the more likely gateway in the exodus to the rest of the world. Using five Ethiopian and three Egyptian high-coverage masked genomes and the multiple sequentially Markovian coalescent (MSMC) approach, we estimated the genetic split times of Egyptians and Ethiopians from non-African populations at 55,000 and 65,000 years ago, respectively, whereas that of West Africans was estimated to be 75,000 years ago. Both the haplotype and MSMC analyses thus suggest a predominant northern route out of Africa via Egypt

Ancient genomes reveal social and genetic structure of Late Neolithic Switzerland

Genetic studies of Neolithic and Bronze Age skeletons from Europe have provided evidence for strong population genetic changes at the beginning and the end of the Neolithic period. To further understand the implications of these in Southern Central Europe, we analyze 96 ancient genomes from Switzerland, Southern Germany, and the Alsace region in France, covering the Middle/Late Neolithic to Early Bronze Age. Similar to previously described genetic changes in other parts of Europe from the early 3rd millennium BCE, we detect an arrival of ancestry related to Late Neolithic pastoralists from the Pontic-Caspian steppe in Switzerland as early as 2860-2460 calBCE. Our analyses suggest that this genetic turnover was a complex process lasting almost 1000 years and involved highly genetically structured populations in this region