A Comparison of Brain Gene Expression Levels in Domesticated and Wild Animals

Domestication has led to similar changes in morphology and behavior in several animal species, raising the question whether similarities between different domestication events also exist at the molecular level. We used mRNA sequencing to analyze genome-wide gene expression patterns in brain frontal cortex in three pairs of domesticated and wild species (dogs and wolves, pigs and wild boars, and domesticated and wild rabbits). We compared the expression differences with those between domesticated guinea pigs and a distant wild relative (Cavia aperea) as well as between two lines of rats selected for tameness or aggression towards humans. There were few gene expression differences between domesticated and wild dogs, pigs, and rabbits (30–75 genes (less than 1%) of expressed genes were differentially expressed), while guinea pigs and C. aperea differed more strongly. Almost no overlap was found between the genes with differential expression in the different domestication events. In addition, joint analyses of all domesticated and wild samples provided only suggestive evidence for the existence of a small group of genes that changed their expression in a similar fashion in different domesticated species. The most extreme of these shared expression changes include up-regulation in domesticates of SOX6 and PROM1, two modulators of brain development. There was almost no overlap between gene expression in domesticated animals and the tame and aggressive rats. However, two of the genes with the strongest expression differences between the rats (DLL3 and DHDH) were located in a genomic region associated with tameness and aggression, suggesting a role in influencing tameness. In summary, the majority of brain gene expression changes in domesticated animals are specific to the given domestication event, suggesting that the causative variants of behavioral domestication traits may likewise be different.     

Regions of Homozygosity in the Porcine Genome: Consequence of Demography and the Recombination Landscape

Inbreeding has long been recognized as a primary cause of fitness reduction in both wild and domesticated populations. Consanguineous matings cause inheritance of haplotypes that are identical by descent (IBD) and result in homozygous stretches along the genome of the offspring. Size and position of regions of homozygosity (ROHs) are expected to correlate with genomic features such as GC content and recombination rate, but also direction of selection. Thus, ROHs should be non-randomly distributed across the genome. Therefore, demographic history may not fully predict the effects of inbreeding. The porcine genome has a relatively heterogeneous distribution of recombination rate, making Sus scrofa an excellent model to study the influence of both recombination landscape and demography on genomic variation. This study utilizes next-generation sequencing data for the analysis of genomic ROH patterns, using a comparative sliding window approach. We present an in-depth study of genomic variation based on three different parameters: nucleotide diversity outside ROHs, the number of ROHs in the genome, and the average ROH size. We identified an abundance of ROHs in all genomes of multiple pigs from commercial breeds and wild populations from Eurasia. Size and number of ROHs are in agreement with known demography of the populations, with population bottlenecks highly increasing ROH occurrence. Nucleotide diversity outside ROHs is high in populations derived from a large ancient population, regardless of current population size. In addition, we show an unequal genomic ROH distribution, with strong correlations of ROH size and abundance with recombination rate and GC content. Global gene content does not correlate with ROH frequency, but some ROH hotspots do contain positive selected genes in commercial lines and wild populations. This study highlights the importance of the influence of demography and recombination on homozygosity in the genome to understand the effects of inbreeding.     

Male and female recombination landscapes of diploid Arabidopsis arenosa

The number and placement of meiotic crossover events during meiosis have important implications for the fidelity of chromosome segregation as well as patterns of inheritance. Despite the functional importance of recombination, recombination landscapes vary widely among and within species, and this can have a strong impact on evolutionary processes. A good knowledge of recombination landscapes is important for model systems in evolutionary and ecological genetics, since it can improve interpretation of genomic patterns of differentiation and genome evolution, and provides an important starting point for understanding the causes and consequences of recombination rate variation. Arabidopsis arenosa is a powerful evolutionary genetic model for studying the molecular basis of adaptation and recombination rate evolution. Here, we generate genetic maps for 2 diploid A. arenosa individuals from distinct genetic lineages where we have prior knowledge that meiotic genes show evidence of selection. We complement the genetic maps with cytological approaches to map and quantify recombination rates, and test the idea that these populations might have distinct patterns of recombination. We explore how recombination differs at the level of populations, individuals, sexes and genomic regions. We show that the positioning of crossovers along a chromosome correlates with their number, presumably a consequence of crossover interference, and discuss how this effect can cause differences in recombination landscape among sexes or species. We identify several instances of female segregation distortion. We found that averaged genome-wide recombination rate is lower and sex differences subtler in A. arenosa than in Arabidopsis thaliana.     

Study of genetic relationships between wild and domesticated grapevine distributed from Middle East Regions to European countries

Archaeobotanical-archaeological, cultural and historical data indicate that grapevine domestication can be dated back from 6000 to 7000 years ago and that it took place in the Caucasian and Middle East Regions. However, events leading to the domestication of this crop species are still an open issue. In this paper, 6 chloroplast microsatellites have been used to assess genetic similarities among, and within, domesticated and wild grapevine accessions representative of 7 distinct geographical regions from the Middle-East to Western Europe. Results show that 2 out of the 6 analyzed chloroplast loci are polymorphic within the 193 domesticated individuals and the 387 samples of 69 wild populations. Allele variants of the Cp-SSR loci combine in a total of 6 different haplotypes. The data show that the haplotype distribution is not homogeneous: the 6 haplotypes are present in the domesticated varieties, but only 5 (haplotype VI is absent) are observed in wild populations. The analysis of haplotype distribution allows discussion of the relationships between the two grape subspecies. The contribution of the wild grape germplasm to the domesticated gene pool still growing in different geographical regions can be, in cases, made evident, suggesting that beside domestication, gene introgression has also played a role in shaping the current varietal landscape of the European viticulture.

---

     

Domestication Shapes Recombination Patterns in Tomato

Meiotic recombination is a biological process of key importance in breeding, to generate genetic diversity and develop novel or agronomically relevant haplotypes. In crop tomato, recombination is curtailed as manifested by linkage disequilibrium decay over a longer distance and reduced diversity compared with wild relatives. Here, we compared domesticated and wild populations of tomato and found an overall conserved recombination landscape, with local changes in effective recombination rate in specific genomic regions. We also studied the dynamics of recombination hotspots resulting from domestication and found that loss of such hotspots is associated with selective sweeps, most notably in the pericentromeric heterochromatin. We detected footprints of genetic changes and structural variants, among them associated with transposable elements, linked with hotspot divergence during domestication, likely causing fine-scale alterations to recombination patterns and resulting in linkage drag.     

The complex history of the domestication of rice

BACKGROUND: Rice has been found in archaeological sites dating to 8000 bc, although the date of rice domestication is a matter of continuing debate. Two species of domesticated rice, Oryza sativa (Asian) and Oryza glaberrima (African) are grown globally. Numerous traits separate wild and domesticated rices including changes in: pericarp colour, dormancy, shattering, panicle architecture, tiller number, mating type and number and size of seeds. SCOPE: Genetic studies using diverse methodologies have uncovered a deep population structure within domesticated rice. Two main groups, the indica and japonica subspecies, have been identified with several subpopulations existing within each group. The antiquity of the divide has been estimated at more than 100 000 years ago. This date far precedes domestication, supporting independent domestications of indica and japonica from pre-differentiated pools of the wild ancestor. Crosses between subspecies display sterility and segregate for domestication traits, indicating that different populations are fixed for different networks of alleles conditioning these traits. Numerous domestication QTLs have been identified in crosses between the subspecies and in crosses between wild and domesticated accessions of rice. Many of the QTLs cluster in the same genomic regions, suggesting that a single gene with pleiotropic effects or that closely linked clusters of genes underlie these QTL. Recently, several domestication loci have been cloned from rice, including the gene controlling pericarp colour and two loci for shattering. The distribution and evolutionary history of these genes gives insight into the domestication process and the relationship between the subspecies. CONCLUSIONS: The evolutionary history of rice is complex, but recent work has shed light on the genetics of the transition from wild (O. rufipogon and O. nivara) to domesticated (O. sativa) rice. The types of genes involved and the geographic and genetic distribution of alleles will allow scientists to better understand our ancestors and breed better rice for our descendents.     

Geography and end use drive the diversification of worldwide winter rye populations

To meet the current challenges in human food production, improved understanding of the genetic diversity of crop species that maximizes the selection efficacy in breeding programs is needed. The present study offers new insights into the diversity, genetic structure and demographic history of cultivated rye (Secale cereale L.). We genotyped 620 individuals from 14 global rye populations with a different end use (grain or forage) at 32 genome‐wide simple sequence repeat markers. We reveal the relationships among these populations, their sizes and the timing of domestication events using population genetics and model‐based inference with approximate Bayesian computation. Our main results demonstrate (i) a high within‐population variation and genetic diversity, (ii) an unexpected absence of reduction in diversity with an increasing improvement level and (iii) patterns suggestive of multiple domestication events. We suggest that the main drivers of diversification of winter rye are the end use of rye in two early regions of cultivation: rye forage in the Mediterranean area and grain in northeast Europe. The lower diversity and stronger differentiation of eastern European populations were most likely due to more intensive cultivation and breeding of rye in this region, in contrast to the Mediterranean region where it was considered a secondary crop or even a weed. We discuss the relevance of our results for the management of gene bank resources and the pitfalls of inference methods applied to crop domestication due to violation of model assumptions and model complexity.     

Genetic mapping of a new race specific resistance allele effective to Puccinia hordei at the Rph9/Rph12locus on chromosome 5HL in barley

Background

Tea is one of the most popular beverages in the world. Many species in the Thea section of the Camellia genus can be processed for drinking and have been domesticated. However, few investigations have focused on the genetic consequence of domestication and geographic origin of landraces on tea plants using credible wild and planted populations of a single species. Here, C. taliensis provides us with a unique opportunity to explore these issues.

Results

Fourteen nuclear microsatellite loci were employed to determine the genetic diversity and domestication origin of C. taliensis, which were represented by 587 individuals from 25 wild, planted and recently domesticated populations. C. taliensis showed a moderate high level of overall genetic diversity. The greater reduction of genetic diversity and stronger genetic drift were detected in the wild group than in the recently domesticated group, indicating the loss of genetic diversity of wild populations due to overexploitation and habitat fragmentation. Instead of the endangered wild trees, recently domesticated individuals were used to compare with the planted trees for detecting the genetic consequence of domestication. A little and non-significant reduction in genetic diversity was found during domestication. The long life cycle, selection for leaf traits and gene flow between populations will delay the emergence of bottleneck in planted trees. Both phylogenetic and assignment analyses suggested that planted trees may have been domesticated from the adjacent central forest of western Yunnan and dispersed artificially to distant places.

Conclusions

This study contributes to the knowledge about levels and distribution of genetic diversity of C. taliensis and provides new insights into genetic consequence of domestication and geographic origin of planted trees of this species. As an endemic tea source plant, wild, planted and recently domesticated C. taliensis trees should all be protected for their unique genetic characteristics, which are valuable for tea breeding.     

Population genomics reveals a fine‐scale recombination landscape for genetic improvement of cotton

Recombination breaks up ancestral linkage disequilibrium, creates combinations of alleles, affects the efficiency of natural selection, and plays a major role in crop domestication and improvement. However, there is little knowledge regarding the variation in the population‐scaled recombination rate in cotton. We constructed recombination maps and characterized the difference in the genomic landscape of the population‐scaled recombination rate between Gossypium hirsutum and G. arboreum and sub‐genomes based on the 381 sequenced G. hirsutum and 215 G. arboreum accessions. Comparative genomics identified large structural variations and syntenic genes in the recombination regions, suggesting that recombination was related to structural variation and occurred preferentially in the distal chromosomal regions. Correlation analysis indicated that recombination was only slightly affected by geographical distribution and breeding period. A genome‐wide association study (GWAS) was performed with 15 agronomic traits using 267 cotton accessions and identified 163 quantitative trait loci (QTL) and an important candidate gene (Ghir_COL2) for early maturity traits. Comparative analysis of recombination and a GWAS revealed that the QTL of fibre quality traits tended to be more common in high‐recombination regions than were those of yield and early maturity traits. These results provide insights into the population‐scaled recombination landscape, suggesting that recombination contributed to the domestication and improvement of cotton, which provides a useful reference for studying recombination in other species.     

Evaluating sample allocation and effort in detecting population differentiation for discrete and continuously distributed individuals

One of the most pressing issues in spatial genetics concerns sampling. Traditionally, substructure and gene flow are estimated for individuals sampled within discrete populations. Because many species may be continuously distributed across a landscape without discrete boundaries, understanding sampling issues becomes paramount. Given large-scale, geographically broad conservation efforts, researchers are looking for guidance as to the trade-offs between sampling more individuals within a population versus few individuals scattered across more populations. Here, we conducted simulations that address these issues. We first established two archetypical patterns of dispersion: (1) individuals within discrete populations, and (2) continuously distributed individuals with limited dispersal. We used genotypes generated from a spatially-explicit, individual-based program and simulated genetic structure in individuals from nine different population sizes across a landscape that either had barriers to movement (defining discrete populations) or isolation-by-distance patterns (defining continuously distributed individuals). Then, given each pattern of dispersion, we allocated samples across four different sampling strategies for each of the nine population sizes in various configurations for sampling more individuals within a population versus fewer individuals scattered across more populations. We assessed the population genetic substructure with both the population-based metric, F _ST, and an individual-based metric, D _PS regardless of the true pattern of dispersion to allow us to better understand the effect of incorrectly matching the metric and the distribution (e.g., F _ST with continuously distributed individuals, and vice versa). We show that sampling many subpopulations (or sampling areas), thus sampling fewer individuals per subpopulation, overestimates measures of population subdivision with the population-based metric for both patterns of dispersion. In contrast, using the individual-based metric gives the opposite results: sampling too few subpopulations, and many individuals per subpopulation, produces an underestimate of the strength of isolation-by-distance. By comparing all results, we were able to suggest a strong predictive model of a chosen genetic structure metric for elucidating the sampling design trade-offs given each pattern of dispersion and configuration on the landscape.     

Natural selection in gene-dense regions shapes the genomic pattern of polymorphism in wild and domesticated rice

Levels of nucleotide variability are frequently positively correlated with recombination rate and negatively associated with gene density due to the effects of selection on linked variation. These relationships are determined by properties that frequently differ among species, including the mating system, and aspects of genome organization such as how genes are distributed along chromosomes. In rice, genes are found at highest density in regions with frequent crossing-over. This association between gene density and recombination rate provides an opportunity to evaluate the effects of selection in a genomic context that differs from other model organisms. Using single-nucleotide polymorphism data from Asian domesticated rice Oryza sativa ssp. japonica and ssp. indica and their progenitor species O. rufipogon, we observe a significant negative association between levels of polymorphism and both gene and coding site density, but either no association, or a negative correlation, between nucleotide variability and recombination rate. We establish that these patterns are unlikely to be explained by neutral mutation rate biases and demonstrate that a model of background selection with variable rates of deleterious mutation is sufficient to account for the gene density effect in O. rufipogon. In O. sativa ssp. japonica, we report a strong negative correlation between polymorphism and recombination rate and greater losses of variation during domestication in the euchromatic chromosome arms than heterochromatin. This is consistent with Hill-Robertson interference in low-recombination regions, which may limit the efficacy of selection for domestication traits. Our results suggest that the physical distribution of selected mutations is a primary factor that determines the genomic pattern of polymorphism in wild and domesticated rice species.     

A Genome-Wide mQTL-seq Scan Identifies Potential Molecular Signatures Regulating Plant Height in Chickpea

The current study employed a high-throughput genome-wide next-generation sequencing-led multiple QTL-seq (mQTL-seq) strategy in two inter- and intra-specific recombinant inbred line (RIL) mapping populations to identify the major genomic regions underlying robust quantitative trait loci (QTLs) regulating plant height in chickpea. The whole genome resequencing discovered 446,475 and 150,434 high-quality homozygous single nucleotide polymorphisms (SNPs) exhibiting polymorphism between tall and dwarf/semi-dwarf mapping parents and bulk/homozygous individuals selected from each of two chickpea RIL populations. These SNP-led mQTL-seq assays in RIL mapping populations scaled-down two longer major genomic regions (1.26–1.34 Mb) underlying robust plant height QTLs into the shorter high-resolution QTL intervals (653.2–756.3 kb) on chickpea chromosomes 3 and 8. This essentially delineated regulatory novel natural SNP allelic variants from brassinosteroid insensitive 1-receptor kinase 1 (BAK1) and gibberellin (GA) 20-oxidase genes governing plant height in chickpea. A strong impact of evolutionary bottlenecks including strong artificial/natural selection on two plant height gene loci during chickpea domestication was observed. The shoot apical meristem-specific expression aside from down-regulation of two plant height genes especially in dwarf/semi-dwarf as compared to tall parents and homozygous mapping individuals of two aforementioned RIL populations was apparent. The integrated genomics-assisted breeding strategy combining mQTL-seq with differential gene expression profiling and functional allelic diversity-based trait domestication study collectively identified potential natural allelic variants of candidate genes underlying major plant height QTLs in chickpea. These functionally relevant molecular signatures can be of immense use for marker-aided genetic enhancement to develop high seed- and pod-yielding non-lodging cultivars restructured with desirable plant height in chickpea.     

Dopamine Receptor Genes and Evolutionary Differentiation in the Domestication of Fighting Cocks and Long-Crowing Chickens

The chicken domestication process represents a typical model of artificial selection, and gives significant insight into the general understanding of the influence of artificial selection on recognizable phenotypes. Two Japanese domesticated chicken varieties, the fighting cock (Shamo) and the long-crowing chicken (Naganakidori), have been selectively bred for dramatically different phenotypes. The former has been selected exclusively for aggressiveness and the latter for long crowing with an obedient sitting posture. To understand the particular mechanism behind these genetic changes during domestication, we investigated the degree of genetic differentiation in the aforementioned chickens, focusing on dopamine receptor D2, D3, and D4 genes. We studied other ornamental chickens such as Chabo chickens as a reference for comparison. When genetic differentiation was measured by an index of nucleotide differentiation (N _ST) newly devised in this study, we found that the N_ST value of DRD4 for Shamo (0.072) was distinctively larger than those of the other genes among the three populations, suggesting that aggressiveness has been selected for in Shamo by collecting a variety of single nucleotide polymorphisms. In addition, we found that in DRD4 in Naganakidori, there is a deletion variant of one proline at the 24^th residue in the repeat of nine prolines of exon 1. We thus conclude that artificial selection has operated on these different kinds of genetic variation in the DRD4 genes of Shamo and Naganakidori so strongly that the two domesticated varieties have differentiated to obtain their present opposite features in a relatively short period of time.     

Largely unlinked gene sets targeted by selection for domestication syndrome phenotypes in maize and sorghum

The domestication of diverse grain crops from wild grasses was a result of artificial selection for a suite of overlapping traits producing changes referred to in aggregate as ‘domestication syndrome’. Parallel phenotypic change can be accomplished by either selection on orthologous genes or selection on non‐orthologous genes with parallel phenotypic effects. To determine how often artificial selection for domestication traits in the grasses targeted orthologous genes, we employed resequencing data from wild and domesticated accessions of Zea (maize) and Sorghum (sorghum). Many ‘classic’ domestication genes identified through quantitative trait locus mapping in populations resulting from wild/domesticated crosses indeed show signatures of parallel selection in both maize and sorghum. However, the overall number of genes showing signatures of parallel selection in both species is not significantly different from that expected by chance. This suggests that while a small number of genes will extremely large phenotypic effects have been targeted repeatedly by artificial selection during domestication, the optimization part of domestication targeted small and largely non‐overlapping subsets of all possible genes which could produce equivalent phenotypic alterations.     

Divergent and linked selection shape patterns of genomic differentiation between European and North American Atlantic salmon (Salmo salar)

As populations diverge many processes can shape genomic patterns of differentiation. Regions of high differentiation can arise due to divergent selection acting on selected loci, genetic hitchhiking of nearby loci, or through repeated selection against deleterious alleles (linked background selection); this divergence may then be further elevated in regions of reduced recombination. Atlantic salmon (Salmo salar) from Europe and North America diverged >600,000 years ago and despite some evidence of secondary contact, the majority of genetic data indicate substantial divergence between lineages. This deep divergence with potential gene flow provides an opportunity to investigate the role of different mechanisms that shape the genomic landscape during early speciation. Here, using 184,295 single nucleotide polymorphisms (SNPs) and 80 populations, we investigate the genomic landscape of differentiation across the Atlantic Ocean with a focus on highly differentiated regions and the processes shaping them. We found evidence of high (mean F_ST = 0.26) and heterogeneous genomic differentiation between continents. Genomic regions associated with high trans‐Atlantic differentiation ranged in size from single loci (SNPs) within important genes to large regions (1–3 Mbp ) on four chromosomes (Ssa06, Ssa13, Ssa16 and Ssa19). These regions showed signatures consistent with selection, including high linkage disequilibrium, despite no significant reduction in recombination. Genes and functional enrichment of processes associated with differentiated regions may highlight continental differences in ocean navigation and parasite resistance. Our results provide insight into potential mechanisms underlying differences between continents, and evidence of near‐fixed and potentially adaptive trans‐Atlantic differences concurrent with a background of high genome‐wide differentiation supports subspecies designation in Atlantic salmon.     

Analysis of a slow-growing line reveals wide genetic variability of carcass and meat quality-related traits

Background

Extensive genome-wide analyses of many human populations, using microarrays containing hundreds of thousands of single-nucleotide polymorphisms, have provided us with abundant information about global genomic diversity. However, these data can also be used to analyze local variability in individual genomic regions. In this study, we analyzed the variability in two genomic regions carrying the genes of the GSTA and GSTM subfamilies, located on different chromosomes.

Results

Analysis of the polymorphisms in GSTA and GSTM gene clusters showed similarities in their allelic and haplotype diversities. These patterns were similar in three Russian populations and the CEU population of European origin. There were statistically significant differences in all the haploblocks of both the GSTM and GSTA regions when the Russian populations were compared with populations from China and Japan. Most haploblocks also differed between the Russians and Nigerians from Yoruba, but, some of them had similar allelic frequencies. Special attention was paid to SNP rs4986947 from the intron of the GSTA4 gene, which is represented in apes by an A nucleotide. In the Asian and African samples, it was represented only by a G allele, and both allelic variants (G/A) occurred in the Russian and European populations.

Conclusions

The results obtained suggest the presence of common features in the evolutionary histories of the GSTA and GSTM gene regions, and that African subpopulations were involved differently in the formation of the European and Asian human lineages.     

Multiple domestication events explain the origin of Gossypium hirsutum landraces in Mexico

Several Mesoamerican crops constitute wild‐to‐domesticated complexes generated by multiple initial domestication events, and continuous gene flow among crop populations and between these populations and their wild relatives. It has been suggested that the domestication of cotton (Gossypium hirsutum) started in the northwest of the Yucatán Peninsula, from where it spread to other regions inside and outside of Mexico. We tested this hypothesis by assembling chloroplast genomes of 23 wild, landraces, and breeding lines (transgene‐introgressed and conventional). The phylogenetic analysis showed that the evolutionary history of cotton in Mexico involves multiple events of introgression and genetic divergence. From this, we conclude that Mexican landraces arose from multiple wild populations. Our results also revealed that their structural and functional chloroplast organizations had been preserved. However, genetic diversity decreases as a consequence of domestication, mainly in transgene‐introgressed (TI) individuals (π = 0.00020, 0.00001, 0.00016, 0, and 0, of wild, TI‐wild, landraces, TI‐landraces, and breeding lines, respectively). We identified homologous regions that differentiate wild from domesticated plants and indicate a relationship among the samples. A decrease in genetic diversity associated with transgene introgression in cotton was identified for the first time, and our outcomes are therefore relevant to both biosecurity and agrobiodiversity conservation.     

A 500-year tale of co-evolution,adaptation, and virulence: Helicobacter pylori in the Americas

Helicobacter pylori is a common component of the human stomach microbiota, possibly dating back to the speciation of Homo sapiens. A history of pathogen evolution in allopatry has led to the development of genetically distinct H. pylori subpopulations, associated with different human populations, and more recent admixture among H. pylori subpopulations can provide information about human migrations. However, little is known about the degree to which some H. pylori genes are conserved in the face of admixture, potentially indicating host adaptation, or how virulence genes spread among different populations. We analyzed H. pylori genomes from 14 countries in the Americas, strains from the Iberian Peninsula, and public genomes from Europe, Africa, and Asia, to investigate how admixture varies across different regions and gene families. Whole-genome analyses of 723 H. pylori strains from around the world showed evidence of frequent admixture in the American strains with a complex mosaic of contributions from H. pylori populations originating in the Americas as well as other continents. Despite the complex admixture, distinctive genomic fingerprints were identified for each region, revealing novel American H. pylori subpopulations. A pan-genome Fst analysis showed that variation in virulence genes had the strongest fixation in America, compared with non-American populations, and that much of the variation constituted non-synonymous substitutions in functional domains. Network analyses suggest that these virulence genes have followed unique evolutionary paths in the American populations, spreading into different genetic backgrounds, potentially contributing to the high risk of gastric cancer in the region.Subject terms: Population genetics, Microbial genetics     

Unraveling the Complex Hybrid Ancestry and Domestication History of Cultivated Strawberry

Cultivated strawberry (Fragaria × ananassa) is one of our youngest domesticates, originating in early eighteenth-century Europe from spontaneous hybrids between wild allo-octoploid species (Fragaria chiloensis and Fragaria virginiana). The improvement of horticultural traits by 300 years of breeding has enabled the global expansion of strawberry production. Here, we describe the genomic history of strawberry domestication from the earliest hybrids to modern cultivars. We observed a significant increase in heterozygosity among interspecific hybrids and a decrease in heterozygosity among domesticated descendants of those hybrids. Selective sweeps were found across the genome in early and modern phases of domestication—59–76% of the selectively swept genes originated in the three less dominant ancestral subgenomes. Contrary to the tenet that genetic diversity is limited in cultivated strawberry, we found that the octoploid species harbor massive allelic diversity and that F. × ananassa harbors as much allelic diversity as either wild founder. We identified 41.8 M subgenome-specific DNA variants among resequenced wild and domesticated individuals. Strikingly, 98% of common alleles and 73% of total alleles were shared between wild and domesticated populations. Moreover, genome-wide estimates of nucleotide diversity were virtually identical in F. chiloensis,F. virginiana, and F. × ananassa (π = 0.0059–0.0060). We found, however, that nucleotide diversity and heterozygosity were significantly lower in modern F. × ananassa populations that have experienced significant genetic gains and have produced numerous agriculturally important cultivars.