Transcriptome population genomics reveals severe bottleneck and domestication cost in the African rice (Oryza glaberrima)

The genomic history and global migration of a windborne pest

Many insect pests, including the brown planthopper (BPH), undergo windborne migration that is challenging to observe and track. It remains controversial about their migration patterns and largely unknown regarding the underlying genetic basis. By analyzing 360 whole genomes from around the globe, we clarify the genetic sources of worldwide BPHs and illuminate a landscape of BPH migration showing that East Asian populations perform closed-circuit journeys between Indochina and the Far East, while populations of Malay Archipelago and South Asia undergo one-way migration to Indochina. We further find round-trip migration accelerates population differentiation, with highly diverged regions enriching in a gene desert chromosome that is simultaneously the speciation hotspot between BPH and related species. This study not only shows the power of applying genomic approaches to demystify the migration in windborne migrants but also enhances our understanding of how seasonal movements affect speciation and evolution in insects.

Porous borders at the wild-crop interface promote weed adaptation in Southeast Asia

High reproductive compatibility between crops and their wild relatives can provide benefits for crop breeding but also poses risks for agricultural weed evolution. Weedy rice is a feral relative of rice that infests paddies and causes severe crop losses worldwide. In regions of tropical Asia where the wild progenitor of rice occurs, weedy rice could be influenced by hybridization with the wild species. Genomic analysis of this phenomenon has been very limited. Here we use whole genome sequence analyses of 217 wild, weedy and cultivated rice samples to show that wild rice hybridization has contributed substantially to the evolution of Southeast Asian weedy rice, with some strains acquiring weed-adaptive traits through introgression from the wild progenitor. Our study highlights how adaptive introgression from wild species can contribute to agricultural weed evolution, and it provides a case study of parallel evolution of weediness in independently-evolved strains of a weedy crop relative.

A comprehensive framework for detecting copy number variants from single nucleotide polymorphism data: 'rCNV', a versatile r package for paralogue and CNV detection

Recent studies have highlighted the significant role of copy number variants (CNVs) in phenotypic diversity, environmental adaptation and species divergence across eukaryotes. The presence of CNVs also has the potential to introduce genotyping biases, which can pose challenges to accurate population and quantitative genetic analyses. However, detecting CNVs in genomes, particularly in non-model organisms, presents a formidable challenge. To address this issue, we have developed a statistical framework and an accompanying r software package that leverage allelic-read depth from single nucleotide polymorphism (SNP) data for accurate CNV detection. Our framework capitalises on two key principles. First, it exploits the distribution of allelic-read depth ratios in heterozygotes for individual SNPs by comparing it against an expected distribution based on binomial sampling. Second, it identifies SNPs exhibiting an apparent excess of heterozygotes under Hardy-Weinberg equilibrium. By employing multiple statistical tests, our method not only enhances sensitivity to sampling effects but also effectively addresses reference biases, resulting in optimised SNP classification. Our framework is compatible with various NGS technologies (e.g. RADseq, Exome-capture). This versatility enables CNV calling from genomes of diverse complexities. To streamline the analysis process, we have implemented our framework in the user-friendly r package 'rCNV', which automates the entire workflow seamlessly. We trained our models using simulated data and validated their performance on four datasets derived from different sequencing technologies, including RADseq (Chinook salmon-Oncorhynchus tshawytscha), Rapture (American lobster-Homarus americanus), Exome-capture (Norway spruce-Picea abies) and WGS (Malaria mosquito-Anopheles gambiae).

Genetic Incompatibilities and Evolutionary Rescue by Wild Relatives Shaped Grain Amaranth Domestication

Crop domestication and the subsequent expansion of crops have long been thought of as a linear process from a wild ancestor to a domesticate. However, evidence of gene flow from locally adapted wild relatives that provided adaptive alleles into crops has been identified in multiple species. Yet, little is known about the evolutionary consequences of gene flow during domestication and the interaction of gene flow and genetic load in crop populations. We study the pseudo-cereal grain amaranth that has been domesticated three times in different geographic regions of the Americas. We quantify the amount and distribution of gene flow and genetic load along the genome of the three grain amaranth species and their two wild relatives. Our results show ample gene flow between crop species and between crops and their wild relatives. Gene flow from wild relatives decreased genetic load in the three crop species. This suggests that wild relatives could provide evolutionary rescue by replacing deleterious alleles in crops. We assess experimental hybrids between the three crop species and found genetic incompatibilities between one Central American grain amaranth and the other two crop species. These incompatibilities might have created recent reproductive barriers and maintained species integrity today. Together, our results show that gene flow played an important role in the domestication and expansion of grain amaranth, despite genetic species barriers. The domestication of plants was likely not linear and created a genomic mosaic by multiple contributors with varying fitness effects for today's crops.

Deleterious Variation in Natural Populations and Implications for Conservation Genetics

Deleterious mutations decrease reproductive fitness and are ubiquitous in genomes. Given that many organisms face ongoing threats of extinction, there is interest in elucidating the impact of deleterious variation on extinction risk and optimizing management strategies accounting for such mutations. Quantifying deleterious variation and understanding the effects of population history on deleterious variation are complex endeavors because we do not know the strength of selection acting on each mutation. Further, the effect of demographic history on deleterious mutations depends on the strength of selection against the mutation and the degree of dominance. Here we clarify how deleterious variation can be quantified and studied in natural populations. We then discuss how different demographic factors, such as small population size, nonequilibrium population size changes, inbreeding, and gene flow, affect deleterious variation. Lastly, we provide guidance on studying deleterious variation in nonmodel populations of conservation concern.

Interspecific Hybridization Is an Important Driving Force for Origin and Diversification of Asian Cultivated Rice Oryza sativa L

As one of the most important crops, Asian cultivated rice has evolved into a complex group including several subgroups adapting various eco-climate-systems around the globe. Here, we pictured a comprehensive view of its original domestication, divergences, and the origin of different subgroups by integrating agriculture, archeology, genetics, nuclear, and cytoplasm genome results. Then, it was highlighted that interspecific hybridization-introgression has played important role in improving the genetic diversity and adaptation of Oryza sativa during its evolution process. Natural hybridization-introgression led to the origin of indica, aus, and basmatic subgroups, which adapted to changing cultivated environments, and produced feral weedy rice coexisting and competing with cultivars under production management. Artificial interspecific hybridization-introgression gained several breakthroughs in rice breeding, such as developing three-line hybrid rice, new rice for Africa (NERICA), and some important pest and disease resistance genes in rice genetic improvement, contributing to the stable increase of rice production to meet the expanding human population. We proposed a series to exploit the virtues of hybridization-introgression in the genetic improvement of Asian cultivated rice. But some key issues such as reproductive barriers especially hybrid sterility should be investigated further, which are conducive to gene exchange between cultivated rice and its relatives, and even is beneficial to exploiting interspecific hybrid vigor. New technologies help introduce favorable genes from distant wild species to Asian cultivated rice, such as transgenic and genome editing systems. Rising introgression lines in a wider range with multi-donor benefits allele mining, understanding genetic network of rice growth and development, yield formation, and environmental adaptation. Then, integration of new tools and interspecific hybridization can be a future direction to develop more usable breeding populations which can make Asian cultivated rice more resilient to the changing climate and world.

Gene fusion as an important mechanism to generate new genes in the genus Oryza

Background

Events of gene fusion have been reported in several organisms. However, the general role of gene fusion as part of new gene origination remains unknown.

Results

We conduct genome-wide interrogations of four Oryza genomes by designing and implementing novel pipelines to detect fusion genes. Based on the phylogeny of ten plant species, we detect 310 fusion genes across four Oryza species. The estimated rate of origination of fusion genes in the Oryza genus is as high as 63 fusion genes per species per million years, which is fixed at 16 fusion genes per species per million years and much higher than that in flies. By RNA sequencing analysis, we find more than 44% of the fusion genes are expressed and 90% of gene pairs show strong signals of purifying selection. Further analysis of CRISPR/Cas9 knockout lines indicates that newly formed fusion genes regulate phenotype traits including seed germination, shoot length and root length, suggesting the functional significance of these genes.

Conclusions

We detect new fusion genes that may drive phenotype evolution in Oryza. This study provides novel insights into the genome evolution of Oryza.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02696-w.

Gene network simulations provide testable predictions for the molecular domestication syndrome

The domestication of plant species lead to repeatable morphological evolution, often referred to as the phenotypic domestication syndrome. Domestication is also associated with important genomic changes, such as the loss of genetic diversity compared to adequately large wild populations, and modifications of gene expression patterns. Here, we explored theoretically the effect of a domestication-like scenario on the evolution of gene regulatory networks. We ran population genetics simulations in which individuals were featured by their genotype (an interaction matrix encoding a gene regulatory network) and their gene expressions, representing the phenotypic level. Our domestication scenario included a population bottleneck and a selection switch mimicking human-mediated directional and canalizing selection, i.e., change in the optimal gene expression level and selection towards more stable expression across environments. We showed that domestication profoundly alters genetic architectures. Based on four examples of plant domestication scenarios, our simulations predict (i) a drop in neutral allelic diversity, (ii) a change in gene expression variance that depends upon the domestication scenario, (iii) transient maladaptive plasticity, (iv) a deep rewiring of the gene regulatory networks, with a trend towards gain of regulatory interactions, and (v) a global increase in the genetic correlations among gene expressions, with a loss of modularity in the resulting coexpression patterns and in the underlying networks. We provide empirically testable predictions on the differences of genetic architectures between wild and domesticated forms. The characterization of such systematic evolutionary changes in the genetic architecture of traits contributes to define a molecular domestication syndrome.

Comparative transcriptome analysis and identification of candidate adaptive evolution genes of Miscanthus lutarioriparius and Miscanthus sacchariflorus

Miscanthus species are perennial C4 grasses that are considered promising energy crops because of their high biomass yields, excellent adaptability and low management costs. Miscanthus lutarioriparius and Miscanthus sacchariflorus are closely related subspecies that are distributed in different habitats. However, there are only a few reports on the mechanisms by which Miscanthus adapts to different environments. Here, comparative transcriptomic and morphological analyses were used to study the evolutionary adaptation of M. lutarioriparius and M. sacchariflorus to different habitats. In total, among 7586 identified orthologs, 2060 orthologs involved in phenylpropanoid biosynthesis and plant hormones were differentially expressed between the two species. Through an analysis of the Ka/Ks ratios of the orthologs, we estimated that the divergence time between the two species was approximately 4.37 Mya. In addition, 37 candidate positively selected orthologs (PSGs) that played important roles in the adaptation of these species to different habitats were identified. Then, the expression levels of 20 PSGs in response to flooding and drought stress were analyzed, and the analysis revealed significant changes in their expression levels. These results facilitate our understanding of the evolutionary adaptation to habitats and the speciation of M. lutarioriparius and M. sacchariflorus. We hypothesise that lignin synthesis genes are the main cause of the morphological differences between the two species. In summary, the plant nonspecific phospholipase C gene family and the receptor-like protein kinase gene family played important roles in the evolution of these two species.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01030-1.

Large-scale genomic analysis reveals the genetic cost of chicken domestication

BACKGROUND

Species domestication is generally characterized by the exploitation of high-impact mutations through processes that involve complex shifting demographics of domesticated species. These include not only inbreeding and artificial selection that may lead to the emergence of evolutionary bottlenecks, but also post-divergence gene flow and introgression. Although domestication potentially affects the occurrence of both desired and undesired mutations, the way wild relatives of domesticated species evolve and how expensive the genetic cost underlying domestication is remain poorly understood. Here, we investigated the demographic history and genetic load of chicken domestication.

RESULTS

We analyzed a dataset comprising over 800 whole genomes from both indigenous chickens and wild jungle fowls. We show that despite having a higher genetic diversity than their wild counterparts (average π, 0.00326 vs. 0.00316), the red jungle fowls, the present-day domestic chickens experienced a dramatic population size decline during their early domestication. Our analyses suggest that the concomitant bottleneck induced 2.95% more deleterious mutations across chicken genomes compared with red jungle fowls, supporting the "cost of domestication" hypothesis. Particularly, we find that 62.4% of deleterious SNPs in domestic chickens are maintained in heterozygous states and masked as recessive alleles, challenging the power of modern breeding programs to effectively eliminate these genetic loads. Finally, we suggest that positive selection decreases the incidence but increases the frequency of deleterious SNPs in domestic chicken genomes.

CONCLUSION

This study reveals a new landscape of demographic history and genomic changes associated with chicken domestication and provides insight into the evolutionary genomic profiles of domesticated animals managed under modern human selection.

Dioecy Is Associated with High Genetic Diversity and Adaptation Rates in the Plant Genus Silene

About 15,000 angiosperm species (∼6%) have separate sexes, a phenomenon known as dioecy. Why dioecious taxa are so rare is still an open question. Early work reported lower species richness in dioecious compared with nondioecious sister clades, raising the hypothesis that dioecy may be an evolutionary dead-end. This hypothesis has been recently challenged by macroevolutionary analyses that detected no or even positive effect of dioecy on diversification. However, the possible genetic consequences of dioecy at the population level, which could drive the long-term fate of dioecious lineages, have not been tested so far. Here, we used a population genomics approach in the Silene genus to look for possible effects of dioecy, especially for potential evidence of evolutionary handicaps of dioecy underlying the dead-end hypothesis. We collected individual-based RNA-seq data from several populations in 13 closely related species with different sexual systems: seven dioecious, three hermaphroditic, and three gynodioecious species. We show that dioecy is associated with increased genetic diversity, as well as higher selection efficacy both against deleterious mutations and for beneficial mutations. The results hold after controlling for phylogenetic inertia, differences in species census population sizes and geographic ranges. We conclude that dioecious Silene species neither show signs of increased mutational load nor genetic evidence for extinction risk. We discuss these observations in the light of the possible demographic differences between dioecious and self-compatible hermaphroditic species and how this could be related to alternatives to the dead-end hypothesis to explain the rarity of dioecy.

Genetics and Genomics of African Rice (Oryza glaberrima Steud) Domestication

African rice (Oryza glaberrima Steud) is one of the two independently domesticated rice species, the other one being Asian rice (Oryza sativa L.). Despite major progress being made in understanding the evolutionary and domestication history of African rice, key outstanding issues remain controversial. There appears to be an underlying difficulty in identifying the domestication centre and number of times the crop has been domesticated. Advances in genomics have provided unprecedented opportunities for understanding the genetic architecture of domestication related traits. For most of the domestication traits, the underlying genes and mutations have been identified. Comparative analysis of domestication genes between Asian and African rice has revealed that the two species went through an independent but convergent evolution process. The genetic and developmental basis of some of the domestic traits are conserved not only between Asian and African rice but also with other domesticated crop species. Analysis of genome data and its interpretation is emerging as a major challenge facing studies of domestication in African rice as key studies continue giving contradictory findings and conclusions. Insights obtained on the domestication of this species are vital for guiding crop improvement efforts.

Harnessing Knowledge from Maize and Rice Domestication for New Crop Breeding

Crop domestication has fundamentally altered the course of human history, causing a shift from hunter-gatherer to agricultural societies and stimulating the rise of modern civilization. A greater understanding of crop domestication would provide a theoretical basis for how we could improve current crops and develop new crops to deal with environmental challenges in a sustainable manner. Here, we provide a comprehensive summary of the similarities and differences in the domestication processes of maize and rice, two major staple food crops that feed the world. We propose that maize and rice might have evolved distinct genetic solutions toward domestication. Maize and rice domestication appears to be associated with distinct regulatory and evolutionary mechanisms. Rice domestication tended to select de novo, loss-of-function, coding variation, while maize domestication more frequently favored standing, gain-of-function, regulatory variation. At the gene network level, distinct genetic paths were used to acquire convergent phenotypes in maize and rice domestication, during which different central genes were utilized, orthologous genes played different evolutionary roles, and unique genes or regulatory modules were acquired for establishing new traits. Finally, we discuss how the knowledge gained from past domestication processes, together with emerging technologies, could be exploited to improve modern crop breeding and domesticate new crops to meet increasing human demands.

The Road to Sorghum Domestication: Evidence From Nucleotide Diversity and Gene Expression Patterns

Native African cereals (sorghum, millets) ensure food security to millions of low-income people from low fertility and drought-prone regions of Africa and Asia. In spite of their agronomic importance, the genetic bases of their phenotype and adaptations are still not well-understood. Here we focus on Sorghum bicolor, which is the fifth cereal worldwide for grain production and constitutes the staple food for around 500 million people. We leverage transcriptomic resources to address the adaptive consequences of the domestication process. Gene expression and nucleotide variability were analyzed in 11 domesticated and nine wild accessions. We documented a downregulation of expression and a reduction of diversity both in nucleotide polymorphism (30%) and gene expression levels (18%) in domesticated sorghum. These findings at the genome-wide level support the occurrence of a global reduction of diversity during the domestication process, although several genes also showed patterns consistent with the action of selection. Nine hundred and forty-nine genes were significantly differentially expressed between wild and domesticated gene pools. Their functional annotation points to metabolic pathways most likely contributing to the sorghum domestication syndrome, such as photosynthesis and auxin metabolism. Coexpression network analyzes revealed 21 clusters of genes sharing similar expression patterns. Four clusters (totaling 2,449 genes) were significantly enriched in differentially expressed genes between the wild and domesticated pools and two were also enriched in domestication and improvement genes previously identified in sorghum. These findings reinforce the evidence that the combined and intricated effects of the domestication and improvement processes do not only affect the behaviors of a few genes but led to a large rewiring of the transcriptome. Overall, these analyzes pave the way toward the identification of key domestication genes valuable for genetic resources characterization and breeding purposes.

Genomic inference of complex domestication histories in three Solanaceae species

Domestication is a human-induced selection process that imprints the genomes of domesticated populations over a short evolutionary time scale and that occurs in a given demographic context. Reconstructing historical gene flow, effective population size changes and their timing is therefore of fundamental interest to understand how plant demography and human selection jointly shape genomic divergence during domestication. Yet, the comparison under a single statistical framework of independent domestication histories across different crop species has been little evaluated so far. Thus, it is unclear whether domestication leads to convergent demographic changes that similarly affect crop genomes. To address this question, we used existing and new transcriptome data on three crop species of Solanaceae (eggplant, pepper and tomato), together with their close wild relatives. We fitted twelve demographic models of increasing complexity on the unfolded joint allele frequency spectrum for each wild/crop pair, and we found evidence for both shared and species-specific demographic processes between species. A convergent history of domestication with gene flow was inferred for all three species, along with evidence of strong reduction in the effective population size during the cultivation stage of tomato and pepper. The absence of any reduction in size of the crop in eggplant stands out from the classical view of the domestication process; as does the existence of a "protracted period" of management before cultivation. Our results also suggest divergent management strategies of modern cultivars among species as their current demography substantially differs. Finally, the timing of domestication is species-specific and supported by the few historical records available.

Evolutionary transcriptomics reveals the origins of olives and the genomic changes associated with their domestication

The olive (Olea europaea L. subsp. europaea) is one of the oldest and most socio-economically important cultivated perennial crop in the Mediterranean region. Yet, its origins are still under debate and the genetic bases of the phenotypic changes associated with its domestication are unknown. We generated RNA-sequencing data for 68 wild and cultivated olive trees to study the genetic diversity and structure both at the transcription and sequence levels. To localize putative genes or expression pathways targeted by artificial selection during domestication, we employed a two-step approach in which we identified differentially expressed genes and screened the transcriptome for signatures of selection. Our analyses support a major domestication event in the eastern part of the Mediterranean basin followed by dispersion towards the West and subsequent admixture with western wild olives. While we found large changes in gene expression when comparing cultivated and wild olives, we found no major signature of selection on coding variants and weak signals primarily affected transcription factors. Our results indicated that the domestication of olives resulted in only moderate genomic consequences and that the domestication syndrome is mainly related to changes in gene expression, consistent with its evolutionary history and life history traits.

Advances in Molecular Genetics and Genomics of African Rice (Oryza glaberrima Steud)

African rice (Oryza glaberrima) has a pool of genes for resistance to diverse biotic and abiotic stresses, making it an important genetic resource for rice improvement. African rice has potential for breeding for climate resilience and adapting rice cultivation to climate change. Over the last decade, there have been tremendous technological and analytical advances in genomics that have dramatically altered the landscape of rice research. Here we review the remarkable advances in knowledge that have been witnessed in the last few years in the area of genetics and genomics of African rice. Advances in cheap DNA sequencing technologies have fuelled development of numerous genomic and transcriptomic resources. Genomics has been pivotal in elucidating the genetic architecture of important traits thereby providing a basis for unlocking important trait variation. Whole genome re-sequencing studies have provided great insights on the domestication process, though key studies continue giving conflicting conclusions and theories. However, the genomic resources of African rice appear to be under-utilized as there seems to be little evidence that these vast resources are being productively exploited for example in practical rice improvement programmes. Challenges in deploying African rice genetic resources in rice improvement and the genomics efforts made in addressing them are highlighted.

Comparisons of molecular diversity indices, selective sweeps and population structure of African rice with its wild progenitor and Asian rice

The extent of molecular diversity parameters across three rice species was compared using large germplasm collection genotyped with genomewide SNPs and SNPs that fell within selective sweep regions. Previous studies conducted on limited number of accessions have reported very low genetic variation in African rice (Oryza glaberrima Steud.) as compared to its wild progenitor (O. barthii A. Chev.) and to Asian rice (O. sativa L.). Here, we characterized a large collection of African rice and compared its molecular diversity indices and population structure with the two other species using genomewide single nucleotide polymorphisms (SNPs) and SNPs that mapped within selective sweeps. A total of 3245 samples representing African rice (2358), Asian rice (772) and O. barthii (115) were genotyped with 26,073 physically mapped SNPs. Using all SNPs, the level of marker polymorphism, average genetic distance and nucleotide diversity in African rice accounted for 59.1%, 63.2% and 37.1% of that of O. barthii, respectively. SNP polymorphism and overall nucleotide diversity of the African rice accounted for 20.1-32.1 and 16.3-37.3% of that of the Asian rice, respectively. We identified 780 SNPs that fell within 37 candidate selective sweeps in African rice, which were distributed across all 12 rice chromosomes. Nucleotide diversity of the African rice estimated from the 780 SNPs was 8.3 × 10^-4, which is not only 20-fold smaller than the value estimated from all genomewide SNPs (π = 1.6 × 10^-2), but also accounted for just 4.1%, 0.9% and 2.1% of that of O. barthii, lowland Asian rice and upland Asian rice, respectively. The genotype data generated for a large collection of rice accessions conserved at the AfricaRice genebank will be highly useful for the global rice community and promote germplasm use.

Origins and geographic diversification of African rice (Oryza glaberrima)

Rice is a staple food for the majority of the world’s population. Whereas Asian rice (Oryza sativa) has been extensively studied, the exact origins of African rice (Oryza glaberrima) are still contested. Previous studies have supported either a centric or a non-centric geographic origin of African rice domestication. Here we review the evidence for both scenarios through a critical reassessment of 206 whole genome sequences of domesticated and wild African rice. While genetic diversity analyses support a severe bottleneck caused by domestication, signatures of recent and strong positive selection do not unequivocally point to candidate domestication genes, suggesting that domestication proceeded differently than in Asian rice–either by selection on different alleles, or different modes of selection. Population structure analysis revealed five genetic clusters localising to different geographic regions. Isolation by distance was identified in the coastal populations, which could account for parallel adaptation in geographically separated demes. Although genome-wide phylogenetic relationships support an origin in the eastern cultivation range followed by diversification along the Atlantic coast, further analysis of domestication genes shows distinct haplotypes in the southwest—suggesting that at least one of several key domestication traits might have originated there. These findings shed new light on an old controversy concerning plant domestication in Africa by highlighting the divergent roots of African rice cultivation, including a separate centre of domestication activity in the Guinea Highlands. We thus suggest that the commonly accepted centric origin of African rice must be reconsidered in favour of a non-centric or polycentric view.

On the Road to Breeding 4.0: Unraveling the Good, the Bad, and the Boring of Crop Quantitative Genomics

Understanding the quantitative genetics of crops has been and will continue to be central to maintaining and improving global food security. We outline four stages that plant breeding either has already achieved or will probably soon achieve. Top-of-the-line breeding programs are currently in Breeding 3.0, where inexpensive, genome-wide data coupled with powerful algorithms allow us to start breeding on predicted instead of measured phenotypes. We focus on three major questions that must be answered to move from current Breeding 3.0 practices to Breeding 4.0: ( a) How do we adapt crops to better fit agricultural environments? ( b) What is the nature of the diversity upon which breeding can act? ( c) How do we deal with deleterious variants? Answering these questions and then translating them to actual gains for farmers will be a significant part of achieving global food security in the twenty-first century.

Using Genomic Data to Infer Historic Population Dynamics of Nonmodel Organisms

Genome sequence data are now being routinely obtained from many nonmodel organisms. These data contain a wealth of information about the demographic history of the populations from which they originate. Many sophisticated statistical inference procedures have been developed to infer the demographic history of populations from this type of genomic data. In this review, we discuss the different statistical methods available for inference of demography, providing an overview of the underlying theory and logic behind each approach. We also discuss the types of data required and the pros and cons of each method. We then discuss how these methods have been applied to a variety of nonmodel organisms. We conclude by presenting some recommendations for researchers looking to use genomic data to infer demographic history.

Genomic approaches for studying crop evolution

Understanding how crop plants evolved from their wild relatives and spread around the world can inform about the origins of agriculture. Here, we review how the rapid development of genomic resources and tools has made it possible to conduct genetic mapping and population genetic studies to unravel the molecular underpinnings of domestication and crop evolution in diverse crop species. We propose three future avenues for the study of crop evolution: establishment of high-quality reference genomes for crops and their wild relatives; genomic characterization of germplasm collections; and the adoption of novel methodologies such as archaeogenetics, epigenomics, and genome editing.

Next-generation metrics for monitoring genetic erosion within populations of conservation concern

Genetic erosion is a major threat to biodiversity because it can reduce fitness and ultimately contribute to the extinction of populations. Here, we explore the use of quantitative metrics to detect and monitor genetic erosion. Monitoring systems should not only characterize the mechanisms and drivers of genetic erosion (inbreeding, genetic drift, demographic instability, population fragmentation, introgressive hybridization, selection) but also its consequences (inbreeding and outbreeding depression, emergence of large-effect detrimental alleles, maladaptation and loss of adaptability). Technological advances in genomics now allow the production of data the can be measured by new metrics with improved precision, increased efficiency and the potential to discriminate between neutral diversity (shaped mainly by population size and gene flow) and functional/adaptive diversity (shaped mainly by selection), allowing the assessment of management-relevant genetic markers. The requirements of such studies in terms of sample size and marker density largely depend on the kind of population monitored, the questions to be answered and the metrics employed. We discuss prospects for the integration of this new information and metrics into conservation monitoring programmes.

De Novo Assemblies of Three Oryza glaberrima Accessions Provide First Insights about Pan-Genome of African Rices

Oryza glaberrima is one of the two cultivated species of rice, and harbors various interesting agronomic traits, especially in biotic and abiotic resistance, compared with its Asian cousin O. sativa. A previous reference genome was published but newer studies highlighted some missing parts. Moreover, global species diversity is known nowadays to be represented by more than one single individual. For that purpose, we sequenced, assembled and annotated de novo three different cultivars from O. glaberrima. After validating our assemblies, we were able to better solve complex regions than the previous assembly and to provide a first insight in pan-genomic divergence between individuals. The three assemblies shown large common regions, but almost 25% of the genome present collinearity breakpoints or are even individual specific.

Elevated Proportions of Deleterious Genetic Variation in Domestic Animals and Plants

A fraction of genetic variants segregating in any population are deleterious, which negatively impacts individual fitness. The domestication of animals and plants is associated with population bottlenecks and artificial selection, which are predicted to increase the proportion of deleterious variants. However, the extent to which this is a general feature of domestic species is unclear. Here, we examine the effects of domestication on the prevalence of deleterious variation using pooled whole-genome resequencing data from five domestic animal species (dog, pig, rabbit, chicken, and silkworm) and two domestic plant species (rice and soybean) compared with their wild ancestors. We find significantly reduced genetic variation and increased proportion of nonsynonymous amino acid changes in all but one of the domestic species. These differences are observable across a range of allele frequencies, both common and rare. We find proportionally more single nucleotide polymorphisms in highly conserved elements in domestic species and a tendency for domestic species to harbor a higher proportion of changes classified as damaging. Our findings most likely reflect an increased incidence of deleterious variants in domestic species, which is most likely attributable to population bottlenecks that lead to a reduction in the efficacy of selection. An exception to this pattern is displayed by European domestic pigs, which do not show traces of a strong population bottleneck and probably continued to exchange genes with wild boar populations after domestication. The results presented here indicate that an elevated proportion of deleterious variants is a common, but not ubiquitous, feature of domestic species.

Sequencing of bulks of segregants allows dissection of genetic control of amylose content in rice

Amylose content (AC) is a key quality trait in rice. A cross between Oryza glaberrima (African rice) and Oryza sativa (Asian rice) segregating for AC was analysed by sequencing bulks of individuals with high and low AC. SNP associated with the granule bound starch synthase (GBSS1) locus on chromosome 6 were polymorphic between the bulks. In particular, a G/A SNP that would result in an Asp to Asn mutation was identified. This amino acid substitution may be responsible for differences in GBSS activity as it is adjacent to a disulphide linkage conserved in all grass GBSS proteins. Other polymorphisms in genomic regions closely surrounding this variation may be the result of linkage drag. In addition to the variant in the starch biosynthesis gene, SNP on chromosomes 1 and 11 linked to AC was also identified. SNP was found in the genes encoding the NAC and CCAAT-HAP5 transcription factors that have previously been linked to starch biosynthesis. This study has demonstrated that the approach of sequencing bulks was able to identify genes on different chromosomes associated with this complex trait.

Population transcriptomic characterization of the genetic and expression variation of a candidate progenitor of Miscanthus energy crops

The use of transcriptome data in the study of the population genetics of a species can capture faint signals of both genetic variation and expression variation and can provide a broad picture of a species' genomic response to environmental conditions. In this study, we characterized the genetic and expression diversity of Miscanthus lutarioriparius by comparing more than 16,225 transcripts obtained from 78 individuals, belonging to 10 populations distributed across the species' entire geographic range. We only observed a low level of nucleotide diversity (π = 0.000434) among the transcriptome data of these populations, which is consistent with highly conserved sequences of functional elements and protein-coding genes captured with this method. Tests of population divergence using the transcriptome data were consistent with previous microsatellite data but proved to be more sensitive, particularly if gene expression variation was considered as well. For example, the analysis of expression data showed that genes involved in photosynthetic processes and responses to temperature or reactive oxygen species stimuli were significantly enriched in certain populations. This differential gene expression was primarily observed among populations and not within populations. Interestingly, nucleotide diversity was significantly negatively correlated with expression diversity within populations, while this correlation was positive among populations. This suggests that genetic and expression variation play separate roles in adaptation and population persistence. Combining analyses of genetic and gene expression variation represents a promising approach for studying the population genetics of wild species and may uncover both adaptive and nonadaptive processes.

Domestication rewired gene expression and nucleotide diversity patterns in tomato

Plant domestication has led to considerable phenotypic modifications from wild species to modern varieties. However, although changes in key traits have been well documented, less is known about the underlying molecular mechanisms, such as the reduction of molecular diversity or global gene co-expression patterns. In this study, we used a combination of gene expression and population genetics in wild and crop tomato to decipher the footprints of domestication. We found a set of 1729 differentially expressed genes (DEG) between the two genetic groups, belonging to 17 clusters of co-expressed DEG, suggesting that domestication affected not only individual genes but also regulatory networks. Five co-expression clusters were enriched in functional terms involving carbohydrate metabolism or epigenetic regulation of gene expression. We detected differences in nucleotide diversity between the crop and wild groups specific to DEG. Our study provides an extensive profiling of the rewiring of gene co-expression induced by the domestication syndrome in one of the main crop species.

Plastic potential: how the phenotypes and adaptations of pathogens are influenced by microbial interactions within plants

Predicting the effects of plant-associated microbes on emergence, spread, and evolution of plant pathogens demands an understanding of how pathogens respond to these microbes at two levels of biological organization: that of an individual pathogen and that of a pathogen population across multiple individual plants. We first examine the plastic responses of individual plant pathogens to microbes within a shared host, as seen through changes in pathogen growth and multiplication. We then explore the limited understanding of how within-plant microbial interactions affect pathogen populations and discuss the need to incorporate population-level observations with population genomic techniques. Finally, we suggest that integrating across levels will further our understanding of the ecological and evolutionary impacts of within-plant microbial interactions on pathogens.

Differences in meristem size and expression of branching genes are associated with variation in panicle phenotype in wild and domesticated African rice

Background

The African rice Oryza glaberrima was domesticated from its wild relative Oryza barthii about 3000 years ago. During the domestication process, panicle complexity changed from a panicle with low complexity in O. barthii, to a highly branched panicle carrying more seeds in O. glaberrima. To understand the basis of this differential panicle development between the two species, we conducted morphological and molecular analyses of early panicle development.

Results

Using X-ray tomography, we analyzed the morphological basis of early developmental stages of panicle development. We uncovered evidence for a wider rachis meristem in O. glaberrima than in O. barthii. At the molecular level, spatial and temporal expression profiles of orthologs of O. sativa genes related to meristem activity and meristem fate control were obtained using in situ hybridization and qRT-PCR. Despite highly conserved spatial expression patterns between O. glaberrima and O. barthii, differences in the expression levels of these early acting genes were detected.

Conclusion

The higher complexity of the O. glaberrima panicle compared to that of its wild relative O. barthii is associated with a wider rachis meristem and a modification of expression of branching-related genes. Our study indicates that the expression of genes in the miR156/miR529/SPL and TAW1 pathways, along with that of their target genes, is altered from the unbranched stage of development. This suggests that differences in panicle complexity between the two African rice species result from early alterations to gene expression during reproductive development.

Electronic supplementary material

The online version of this article (doi:10.1186/s13227-017-0065-y) contains supplementary material, which is available to authorized users.

High Level of Nonsynonymous Changes in Common Bean Suggests That Selection under Domestication Increased Functional Diversity at Target Traits

Crop species have been deeply affected by the domestication process, and there have been many efforts to identify selection signatures at the genome level. This knowledge will help geneticists to better understand the evolution of organisms, and at the same time, help breeders to implement successful breeding strategies. Here, we focused on domestication in the Mesoamerican gene pool of Phaseolus vulgaris by sequencing 49 gene fragments from a sample of 45 P. vulgaris wild and domesticated accessions, and as controls, two accessions each of the closely related species Phaseolus coccineus and Phaseolus dumosus. An excess of nonsynonymous mutations within the domesticated germplasm was found. Our data suggest that the cost of domestication alone cannot explain fully this finding. Indeed, the significantly higher frequency of polymorphisms in the coding regions observed only in the domesticated plants (compared to noncoding regions), the fact that these mutations were mostly nonsynonymous and appear to be recently derived mutations, and the investigations into the functions of their relative genes (responses to biotic and abiotic stresses), support a scenario that involves new functional mutations selected for adaptation during domestication. Moreover, consistent with this hypothesis, selection analysis and the possibility to compare data obtained for the same genes in different studies of varying sizes, data types, and methodologies allowed us to identify four genes that were strongly selected during domestication. Each selection candidate is involved in plant resistance/tolerance to abiotic stresses, such as heat, drought, and salinity. Overall, our study suggests that domestication acted to increase functional diversity at target loci, which probably controlled traits related to expansion and adaptation to new agro-ecological growing conditions.

miR2118-triggered phased siRNAs are differentially expressed during the panicle development of wild and domesticated African rice species

BACKGROUND

Rice exhibits a wide range of panicle structures. To explain these variations, much emphasis has been placed on changes in transcriptional regulation, but no large-scale study has yet reported on changes in small RNA regulation in the various rice species. To evaluate this aspect, we performed deep sequencing and expression profiling of small RNAs from two closely related species with contrasting panicle development: the cultivated African rice Oryza glaberrima and its wild relative Oryza barthii.

RESULTS

Our RNA-seq analysis revealed a dramatic difference between the two species in the 21 nucleotide small RNA population, corresponding mainly to miR2118-triggered phased siRNAs. A detailed expression profiling during the panicle development of O. glaberrima and O. barthii using qRT-PCRs and in situ hybridization, confirmed a delayed expression of the phased siRNAs as well as their lncRNA precursors and regulators (miR2118 and MEL1 gene) in O. glaberrima compared to O. barthii. We provide evidence that the 21-nt phasiRNA pathway in rice is associated with male-gametogenesis but is initiated in spikelet meristems.

CONCLUSION

Differential expression of the miR2118-triggered 21-nt phasiRNA pathway between the two African rice species reflects differential rates of determinate fate acquisition of panicle meristems between the two species.

Taming the Past: Ancient DNA and the Study of Animal Domestication

During the last decade, ancient DNA research has been revolutionized by the availability of increasingly powerful DNA sequencing and ancillary genomics technologies, giving rise to the new field of paleogenomics. In this review, we show how our understanding of the genetic basis of animal domestication and the origins and dispersal of livestock and companion animals during the Upper Paleolithic and Neolithic periods is being rapidly transformed through new scientific knowledge generated with paleogenomic methods. These techniques have been particularly informative in revealing high-resolution patterns of artificial and natural selection and evidence for significant admixture between early domestic animal populations and their wild congeners.

A new eIF4E1 allele characterized by RNAseq data mining is associated with resistance to potato virus Y in tomato albeit with a low durability

Allele mining on susceptibility factors offers opportunities to find new sources of resistance among crop wild relatives for breeding purposes. As a proof of concept, we used available RNAseq data to investigate polymorphisms among the four tomato genes encoding translation initiation factors [eIF4E1 and eIF4E2, eIFiso4E and the related gene new cap-binding protein(nCBP)] to look for new potential resistance alleles to potyviruses. By analysing polymorphism among RNAseq data obtained for 20 tomato accessions, 10 belonging to the cultivated type Solanum lycopersicum and 10 belonging to the closest related wild species Solanum pimpinellifolium, we isolated one new eIF4E1 allele, in the S. pimpinellifolium LA0411 accession, which encodes a potential new resistance allele, mainly due to a polymorphism associated with an amino acid change within eIF4E1 region II. We confirmed that this new allele, pot12, is indeed associated with resistance to potato virus Y, although with a restricted resistance spectrum and a very low durability potential. This suggests that mutations occurring in eIF4E region II only may not be sufficient to provide efficient and durable resistance in plants. However, our study emphasizes the opportunity brought by RNAseq data to mine for new resistance alleles. Moreover, this approach could be extended to seek for putative new resistance alleles by screening for variant forms of susceptibility genes encoding plant host proteins known to interact with viral proteins.

A reverse genetic approach identifies an ancestral frameshift mutation in RP1 causing recessive progressive retinal degeneration in European cattle breeds

Background

Domestication and artificial selection have resulted in strong genetic drift, relaxation of purifying selection and accumulation of deleterious mutations. As a consequence, bovine breeds experience regular outbreaks of recessive genetic defects which might represent only the tip of the iceberg since their detection depends on the observation of affected animals with distinctive symptoms. Thus, recessive mutations resulting in embryonic mortality or in non-specific symptoms are likely to be missed. The increasing availability of whole-genome sequences has opened new research avenues such as reverse genetics for their investigation. Our aim was to characterize the genetic load of 15 European breeds using data from the 1000 bull genomes consortium and prove that widespread harmful mutations remain to be detected.

Results

We listed 2489 putative deleterious variants (in 1923 genes) segregating at a minimal frequency of 5 % in at least one of the breeds studied. Gene enrichment analysis showed major enrichment for genes related to nervous, visual and auditory systems, and moderate enrichment for genes related to cardiovascular and musculoskeletal systems. For verification purposes, we investigated the phenotypic consequences of a frameshift variant in the retinitis pigmentosa-1 gene segregating in several breeds and at a high frequency (27 %) in Normande cattle. As described in certain human patients, clinical and histological examination revealed that this mutation causes progressive degeneration of photoreceptors leading to complete blindness in homozygotes. We established that the deleterious allele was even more frequent in the Normande breed before 1975 (>40 %) and has been progressively counter-selected likely because of its associated negative effect on udder morphology. Finally, using identity-by-descent analysis we demonstrated that this mutation resulted from a unique ancestral event that dates back to ~2800 to 4000 years.

Conclusions

We provide a list of mutations that likely represent a substantial part of the genetic load of domestication in European cattle. We demonstrate that they accumulated non-randomly and that genes related to cognition and sensory functions are particularly affected. Finally, we describe an ancestral deleterious variant segregating in different breeds causing progressive retinal degeneration and irreversible blindness in adult animals.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-016-0232-y) contains supplementary material, which is available to authorized users.

Domestication history and geographical adaptation inferred from a SNP map of African rice

African rice (Oryza glaberrima Steud.) is a cereal crop species closely related to Asian rice (Oryza sativa L.) but was independently domesticated in West Africa ∼3,000 years ago. African rice is rarely grown outside sub-Saharan Africa but is of global interest because of its tolerance to abiotic stresses. Here we describe a map of 2.32 million SNPs of African rice from whole-genome resequencing of 93 landraces. Population genomic analysis shows a population bottleneck in this species that began ∼13,000-15,000 years ago with effective population size reaching its minimum value ∼3,500 years ago, suggesting a protracted period of population size reduction likely commencing with predomestication management and/or cultivation. Genome-wide association studies (GWAS) for six salt tolerance traits identify 11 significant loci, 4 of which are within ∼300 kb of genomic regions that possess signatures of positive selection, suggesting adaptive geographical divergence for salt tolerance in this species.

Superheroes and masterminds of plant domestication

Domestication is one of the most fundamental changes in the evolution of human societies. The geographical origins of domesticated plants are inferred from archaeology, ecology and genetic data. Scenarios vary among species and include single, diffuse or multiple independent domestications. Cultivated plants present a panel of traits, the "domestication syndrome" that distinguish them from their wild relatives. It encompasses yield-, food usage-, and cultivation-related traits. Most genes underlying those traits are "masterminds" affecting the regulation of gene networks. Phenotypic convergence of domestication traits across species or within species between independently domesticated forms rarely coincides with convergence at the gene level. We review here current data/models that propose a protracted transition model for domestication and investigate the impact of mating system, life cycle and gene flow on the pace of domestication. Finally, we discuss the cost of domestication, pointing to the importance of characterizing adaptive functional variation in wild resources.

Determinants of genetic diversity

Genetic polymorphism varies among species and within genomes, and has important implications for the evolution and conservation of species. The determinants of this variation have been poorly understood, but population genomic data from a wide range of organisms now make it possible to delineate the underlying evolutionary processes, notably how variation in the effective population size (Ne) governs genetic diversity. Comparative population genomics is on its way to providing a solution to 'Lewontin's paradox' - the discrepancy between the many orders of magnitude of variation in population size and the much narrower distribution of diversity levels. It seems that linked selection plays an important part both in the overall genetic diversity of a species and in the variation in diversity within the genome. Genetic diversity also seems to be predictable from the life history of a species.

The Limits of Natural Selection in a Nonequilibrium World

Evolutionary theory predicts that factors such as a small population size or low recombination rate can limit the action of natural selection. The emerging field of comparative population genomics offers an opportunity to evaluate these hypotheses. However, classical theoretical predictions assume that populations are at demographic equilibrium. This assumption is likely to be violated in the very populations researchers use to evaluate selection's limits: populations that have experienced a recent shift in population size and/or effective recombination rates. Here we highlight theory and data analyses concerning limitations on the action of natural selection in nonequilibrial populations and argue that substantial care is needed to appropriately test whether species and populations show meaningful differences in selection efficacy. A move toward model-based inferences that explicitly incorporate nonequilibrium dynamics provides a promising approach to more accurately contrast selection efficacy across populations and interpret its significance.

Comparative transcriptome analysis highlights the crucial roles of photosynthetic system in drought stress adaptation in upland rice

Drought stress is one of the major adverse environmental factors reducing plant growth. With the aim to elucidate the underlying molecular basis of rice response to drought stress, comparative transcriptome analysis was conducted between drought susceptible rice cultivar Zhenshan97 and tolerant cultivar IRAT109 at the seedling stage. 436 genes showed differential expression and mainly enriched in the Gene Ontology (GO) terms of stress defence. A large number of variations exist between these two genotypes including 2564 high-quality insertion and deletions (INDELs) and 70,264 single nucleotide polymorphism (SNPs). 1041 orthologous gene pairs show the ratio of nonsynonymous nucleotide substitution rate to synonymous nucleotide substitutions rate (Ka/Ks) larger than 1.5, indicating the rapid adaptation to different environments during domestication. GO and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of positive selection genes suggested that photosynthesis represents the most significant category. The collocation of positively selected genes with the QTLs of photosynthesis and the different photosynthesis performance of these two cultivars further illuminate the crucial function of photosynthesis in rice adaptation to drought stress. Our results also provide fruitful functional markers and candidate genes for future genetic research and improvement of drought tolerance in rice.

Widespread and Adaptive Alterations in Genome-Wide Gene Expression Associated with Ecological Divergence of Two Oryza Species

Ecological speciation is a common mechanism by which new species arise. Despite great efforts, the role of gene expression in ecological divergence and speciation is poorly understood. Here, we conducted a genome-wide gene expression investigation of two Oryza species that are evolutionarily young and distinct in ecology and morphology. Using digital gene expression technology and the paired-end RNA sequencing method, we obtained 21,415 expressed genes across three reproduction-related tissues. Of them, approximately 8% (1,717) differed significantly in expression levels between the two species and these differentially expressed genes are randomly distributed across the genome. Moreover, 62% (1,064) of the differentially expressed genes exhibited a signature of directional selection in at least one species. Importantly, the genes with differential expression between species evolved more rapidly at the 5' flanking sequences than the genes without differential expression relative to coding sequences, suggesting that cis-regulatory changes are likely adaptive and play an important role in the ecological divergence of the two species. Finally, we showed evidence of significant differentiation between species in phenotype traits and observed that genes with differential expression were overrepresented with functional terms involving phenotypic and ecological differentiation between the two species, including reproduction- and stress-related characteristics. Our findings demonstrate that ecological speciation is associated with widespread and adaptive alterations in genome-wide gene expression and provide new insights into the importance of regulatory evolution in ecological speciation in plants.

Chromosomal passports provide new insights into diffusion of emmer wheat

Emmer wheat, Triticum dicoccon schrank (syn. T. dicoccum (schrank) schÜbl.), is one of the earliest domesticated crops, harboring a wide range of genetic diversity and agronomically valuable traits. The crop, however, is currently largely neglected. We provide a wealth of karyotypic information from a comprehensive collection of emmer wheat and related taxa. In addition to C-banding polymorphisms, we identified 43 variants of chromosomal rearrangements in T. dicoccon; among them 26 (60.4%) were novel. The T7A:5B translocation was most abundant in Western Europe and the Mediterranean. The plant genetic resources investigated here might become important in the future for wheat improvement. Based on cluster analysis four major karyotypic groups were discriminated within the T. dicoccon genepool, each harboring characteristic C-banding patterns and translocation spectra: the balkan, asian, european and ethiopian groups. We postulate four major diffusion routes of the crop and discuss their migration out of the Fertile Crescent considering latest archaeobotanical findings.

Molecular Evolution of Freshwater Snails with Contrasting Mating Systems

Because mating systems affect population genetics and ecology, they are expected to impact the molecular evolution of species. Self-fertilizing species experience reduced effective population size, recombination rates, and heterozygosity, which in turn should decrease the efficacy of natural selection, both adaptive and purifying, and the strength of meiotic drive processes such as GC-biased gene conversion. The empirical evidence is only partly congruent with these predictions, depending on the analyzed species, some, but not all, of the expected effects have been observed. One possible reason is that self-fertilization is an evolutionary dead-end, so that most current selfers recently evolved self-fertilization, and their genome has not yet been strongly impacted by selfing. Here, we investigate the molecular evolution of two groups of freshwater snails in which mating systems have likely been stable for several millions of years. Analyzing coding sequence polymorphism, divergence, and expression levels, we report a strongly reduced genetic diversity, decreased efficacy of purifying selection, slower rate of adaptive evolution, and weakened codon usage bias/GC-biased gene conversion in the selfer Galba compared with the outcrosser Physa, in full agreement with theoretical expectations. Our results demonstrate that self-fertilization, when effective in the long run, is a major driver of population genomic and molecular evolutionary processes. Despite the genomic effects of selfing, Galba truncatula seems to escape the demographic consequences of the genetic load. We suggest that the particular ecology of the species may buffer the negative consequences of selfing, shedding new light on the dead-end hypothesis.

Prehistoric genomes reveal the genetic foundation and cost of horse domestication

The domestication of the horse ∼ 5.5 kya and the emergence of mounted riding, chariotry, and cavalry dramatically transformed human civilization. However, the genetics underlying horse domestication are difficult to reconstruct, given the near extinction of wild horses. We therefore sequenced two ancient horse genomes from Taymyr, Russia (at 7.4- and 24.3-fold coverage), both predating the earliest archeological evidence of domestication. We compared these genomes with genomes of domesticated horses and the wild Przewalski's horse and found genetic structure within Eurasia in the Late Pleistocene, with the ancient population contributing significantly to the genetic variation of domesticated breeds. We furthermore identified a conservative set of 125 potential domestication targets using four complementary scans for genes that have undergone positive selection. One group of genes is involved in muscular and limb development, articular junctions, and the cardiac system, and may represent physiological adaptations to human utilization. A second group consists of genes with cognitive functions, including social behavior, learning capabilities, fear response, and agreeableness, which may have been key for taming horses. We also found that domestication is associated with inbreeding and an excess of deleterious mutations. This genetic load is in line with the "cost of domestication" hypothesis also reported for rice, tomatoes, and dogs, and it is generally attributed to the relaxation of purifying selection resulting from the strong demographic bottlenecks accompanying domestication. Our work demonstrates the power of ancient genomes to reconstruct the complex genetic changes that transformed wild animals into their domesticated forms, and the population context in which this process took place.