Whole-genome sequencing and comparative genomics reveal candidate genes associated with quality traits in Dioscorea alata.

BACKGROUND: Quality traits are essential determinants of consumer preferences. Dioscorea alata (Greater Yam), is a starchy tuber crop in tropical regions. However, a comprehensive understanding of the genetic basis underlying yam tuber quality remains elusive. To address this knowledge gap, we employed population genomics and candidate gene association approaches to unravel the genetic factors influencing the quality attributes of boiled yam. METHODS AND RESULTS: Comparative genomics analysis of 45 plant species revealed numerous novel genes absent in the existing D. alata gene annotation. This approach, adding 48% more genes, significantly enhanced the functional annotation of three crucial metabolic pathways associated with boiled yam quality traits: pentose and glucuronate interconversions, starch and sucrose metabolism, and flavonoid biosynthesis. In addition, the whole-genome sequencing of 127 genotypes identified 27 genes under selection and 22 genes linked to texture, starch content, and color through a candidate gene association analysis. Notably, five genes involved in starch content and cell wall composition, including 1,3-beta Glucan synthase, ß-amylase, and Pectin methyl esterase, were common to both approaches and their expression levels were assessed by transcriptomic data. CONCLUSIONS: The analysis of the whole-genome of 127 genotypes of D. alata and the study of three specific pathways allowed the identification of important genes for tuber quality. Our findings provide insights into the genetic basis of yam quality traits and will help the enhancement of yam tuber quality through breeding programs.

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.

Transferability, development of simple sequence repeat (SSR) markers and application to the analysis of genetic diversity and population structure of the African fan palm (Borassus aethiopum Mart.) in Benin.

BACKGROUND: In Sub-Saharan Africa, Borassus aethiopum Mart. (African fan palm) is an important non-timber forest product-providing palm that faces multiple anthropogenic threats to its genetic diversity. However, this species is so far under-studied, which prevents its sustainable development as a resource. The present work is a first attempt at characterizing the genetic diversity and population structure of B. aethiopum across nine collection sites spanning the three climatic regions of Benin, West Africa, through the use of microsatellite markers. RESULTS: During a first phase we relied on the reported transferability of primers developed in other palm species. We find that, in disagreement with previously published results, only 22.5% of the markers tested enable amplification of B. aethiopum DNA and polymorphism detection is very low. In a second phase, we generated a B. aethiopum-specific genomic dataset through high-throughput sequencing and used it for the de novo detection of microsatellite loci. Among the primer pairs targeting these, 11 detected polymorphisms and were further used for analyzing genetic diversity. Across the nine sites, expected heterozygosity (He) ranges from 0.263 to 0.451 with an overall average of 0.354, showing a low genetic diversity. Analysis of molecular variance (AMOVA) shows that within-site variation accounts for 53% of the genetic variation. Accordingly, the low number of migrants and positive values of the fixation index (F) in sites from both the Central (Sudano-Guinean) and the Southern (Guinean) climatic regions suggest limited gene flow between sites. The global correlation between genetic and geographic distances is weak; however, our clustering analyses indicate that B. aethiopum palms from Savè (Center) are genetically more similar to those from the North than to samples from other Central sites. CONCLUSIONS: In the light of our results, we discuss the use of inter-species transfer vs. de novo development of microsatellite markers in genetic diversity analyses targeting under-studied species, and suggest future applications for our molecular resources. We propose that, while prominent short-range pollen and seed dispersal in Benin explain most of our results, gene flux between the Central and Northern regions, as a result of animal and/or human migrations, might underlie the Savè discrepancy.

Evolutionary forces affecting synonymous variations in plant genomes.

Base composition is highly variable among and within plant genomes, especially at third codon positions, ranging from GC-poor and homogeneous species to GC-rich and highly heterogeneous ones (particularly Monocots). Consequently, synonymous codon usage is biased in most species, even when base composition is relatively homogeneous. The causes of these variations are still under debate, with three main forces being possibly involved: mutational bias, selection and GC-biased gene conversion (gBGC). So far, both selection and gBGC have been detected in some species but how their relative strength varies among and within species remains unclear. Population genetics approaches allow to jointly estimating the intensity of selection, gBGC and mutational bias. We extended a recently developed method and applied it to a large population genomic dataset based on transcriptome sequencing of 11 angiosperm species spread across the phylogeny. We found that at synonymous positions, base composition is far from mutation-drift equilibrium in most genomes and that gBGC is a widespread and stronger process than selection. gBGC could strongly contribute to base composition variation among plant species, implying that it should be taken into account in plant genome analyses, especially for GC-rich ones.

Allele-specific expression and genetic determinants of transcriptomic variations in response to mild water deficit in tomato.

Characterizing the natural diversity of gene expression across environments is an important step in understanding how genotype-by-environment interactions shape phenotypes. Here, we analyzed the impact of water deficit onto gene expression levels in tomato at the genome-wide scale. We sequenced the transcriptome of growing leaves and fruit pericarps at cell expansion stage in a cherry and a large fruited accession and their F1 hybrid grown under two watering regimes. Gene expression levels were steadily affected by the genotype and the watering regime. Whereas phenotypes showed mostly additive inheritance, ~80% of the genes displayed non-additive inheritance. By comparing allele-specific expression (ASE) in the F1 hybrid to the allelic expression in both parental lines, respectively, 3005 genes in leaf and 2857 genes in fruit deviated from 1:1 ratio independently of the watering regime. Among these genes, ~55% were controlled by cis factors, ~25% by trans factors and ~20% by a combination of both types of factors. A total of 328 genes in leaf and 113 in fruit exhibited significant ASE-by-watering regime interaction, among which ~80% presented trans-by-watering regime interaction, suggesting a response to water deficit mediated through a majority of trans-acting loci in tomato. We cross-validated the expression levels of 274 transcripts in fruit and leaves of 124 recombinant inbred lines (RILs) and identified 163 expression quantitative trait loci (eQTLs) mostly confirming the divergences identified by ASE. Combining phenotypic and expression data, we observed a complex network of variation between genes encoding enzymes involved in the sugar metabolism.

A novel high-density grapevine (Vitis vinifera L.) integrated linkage map using GBS in a half-diallel population.

KEY MESSAGE: A half-diallel population involving five elite grapevine cultivars was generated and genotyped by GBS, and highly-informative segregation data was used to construct a high-density genetic map for Vitis vinifera L. Grapevine is one of the most relevant fruit crops in the world. Deeper genetic knowledge could assist modern grapevine breeding programs to develop new wine grape varieties able to face climate change effects. To assist in the rapid identification of markers for crop yield components, grape quality traits and adaptation potential, we generated a large Vitis vinifera L. population (N = 624) by crossing five red wine cultivars in a half-diallel scheme, which was subsequently sequenced by an efficient GBS procedure. A high number of fully informative genetic variants was detected using a novel mapping approach capable of reconstructing local haplotypes from adjacent biallelic SNPs, which were subsequently used to construct the densest consensus genetic map available for the cultivated grapevine to date. This 1378.3-cM map integrates 10 bi-parental consensus maps and orders 4437 markers in 3353 unique positions on 19 chromosomes. Markers are well distributed all along the grapevine reference genome, covering up to 98.8% of its genomic sequence. Additionally, a good agreement was observed between genetic and physical orders, adding confidence in the quality of this map. Collectively, our results pave the way for future genetic studies (such as fine QTL mapping) aimed to understand the complex relationship between genotypic and phenotypic variation in the cultivated grapevine. In addition, the method used (which efficiently delivers a high number of fully informative markers) could be of interest to other outbred organisms, notably perennial fruit crops.

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.

Epistatic determinism of durum wheat resistance to the wheat spindle streak mosaic virus.

KEY MESSAGE: The resistance of durum wheat to the Wheat spindle streak mosaic virus (WSSMV) is controlled by two main QTLs on chromosomes 7A and 7B, with a huge epistatic effect. Wheat spindle streak mosaic virus (WSSMV) is a major disease of durum wheat in Europe and North America. Breeding WSSMV-resistant cultivars is currently the only way to control the virus since no treatment is available. This paper reports studies of the inheritance of WSSMV resistance using two related durum wheat populations obtained by crossing two elite cultivars with a WSSMV-resistant emmer cultivar. In 2012 and 2015, 354 recombinant inbred lines (RIL) were phenotyped using visual notations, ELISA and qPCR and genotyped using locus targeted capture and sequencing. This allowed us to build a consensus genetic map of 8568 markers and identify three chromosomal regions involved in WSSMV resistance. Two major regions (located on chromosomes 7A and 7B) jointly explain, on the basis of epistatic interactions, up to 43% of the phenotypic variation. Flanking sequences of our genetic markers are provided to facilitate future marker-assisted selection of WSSMV-resistant cultivars.

A replicated climate change field experiment reveals rapid evolutionary response in an ecologically important soil invertebrate.

Whether species can respond evolutionarily to current climate change is crucial for the persistence of many species. Yet, very few studies have examined genetic responses to climate change in manipulated experiments carried out in natural field conditions. We examined the evolutionary response to climate change in a common annelid worm using a controlled replicated experiment where climatic conditions were manipulated in a natural setting. Analyzing the transcribed genome of 15 local populations, we found that about 12% of the genetic polymorphisms exhibit differences in allele frequencies associated to changes in soil temperature and soil moisture. This shows an evolutionary response to realistic climate change happening over short-time scale, and calls for incorporating evolution into models predicting future response of species to climate change. It also shows that designed climate change experiments coupled with genome sequencing offer great potential to test for the occurrence (or lack) of an evolutionary response.

Endemic insular and coastal Tunisian date palm genetic diversity.

The breeding of crop species relies on the valorisation of ancestral or wild varieties to enrich the cultivated germplasm. The Tunisian date palm genetic patrimony is being threatened by diversity loss and global climate change. We have conducted a genetic study to evaluate the potential of spontaneous coastal resources to improve the currently exploited Tunisian date palm genetic pool. Eighteen microsatellite loci of Phoenix dactylifera L. were used to compare the genetic diversity of coastal accessions from Kerkennah, Djerba, Gabès and continental date palm accessions from Tozeur. A collection of 105 date palms from the four regions was analysed. This study has provided us with an extensive understanding of the local genetic diversity and its distribution. The coastal date palm genotypes exhibit a high and specific genetic diversity. These genotypes are certainly an untapped reservoir of agronomically important genes to improve cultivated germplasm in continental date palm.

Genetic structure of the date palm (Phoenix dactylifera) in the Old World reveals a strong differentiation between eastern and western populations.

Since the publication of this paper, it has become apparent that an error was made in the legend to Fig. 3 and the colours referring to occidental and oriental are the wrong way round. The authors apologise for this error, and a correct version of the legend to Fig. 3 is given below.

Genetic structure of the date palm (Phoenix dactylifera) in the Old World reveals a strong differentiation between eastern and western populations.

BACKGROUND AND AIMS: Date palms (Phoenix dactylifera, Arecaceae) are of great economic and ecological value to the oasis agriculture of arid and semi-arid areas. However, despite the availability of a large date palm germplasm spreading from the Atlantic shores to Southern Asia, improvement of the species is being hampered by a lack of information on global genetic diversity and population structure. In order to contribute to the varietal improvement of date palms and to provide new insights on the influence of geographic origins and human activity on the genetic structure of the date palm, this study analysed the diversity of the species. METHODS: Genetic diversity levels and population genetic structure were investigated through the genotyping of a collection of 295 date palm accessions ranging from Mauritania to Pakistan using a set of 18 simple sequence repeat (SSR) markers and a plastid minisatellite. KEY RESULTS: Using a Bayesian clustering approach, the date palm genotypes can be structured into two different gene pools: the first, termed the Eastern pool, consists of accessions from Asia and Djibouti, whilst the second, termed the Western pool, consists of accessions from Africa. These results confirm the existence of two ancient gene pools that have contributed to the current date palm diversity. The presence of admixed genotypes is also noted, which points at gene flows between eastern and western origins, mostly from east to west, following a human-mediated diffusion of the species. CONCLUSIONS: This study assesses the distribution and level of genetic diversity of accessible date palm resources, provides new insights on the geographic origins and genetic history of the cultivated component of this species, and confirms the existence of at least two domestication origins. Furthermore, the strong genetic structure clearly established here is a prerequisite for any breeding programme exploiting the effective polymorphism related to each gene pool.

A small XY chromosomal region explains sex determination in wild dioecious V. vinifera and the reversal to hermaphroditism in domesticated grapevines.

BACKGROUND: In Vitis vinifera L., domestication induced a dramatic change in flower morphology: the wild sylvestris subspecies is dioecious while hermaphroditism is largely predominant in the domesticated subsp. V. v. vinifera. The characterisation of polymorphisms in genes underlying the sex-determining chromosomal region may help clarify the history of domestication in grapevine and the evolution of sex chromosomes in plants. In the genus Vitis, sex determination is putatively controlled by one major locus with three alleles, male M, hermaphrodite H and female F, with an allelic dominance M > H > F. Previous genetic studies located the sex locus on chromosome 2. We used DNA polymorphisms of geographically diverse V. vinifera genotypes to confirm the position of this locus, to characterise the genetic diversity and traces of selection in candidate genes, and to explore the origin of hermaphroditism. RESULTS: In V. v. sylvestris, a sex-determining region of 154.8 kb, also present in other Vitis species, spans less than 1% of chromosome 2. It displays haplotype diversity, linkage disequilibrium and differentiation that typically correspond to a small XY sex-determining region with XY males and XX females. In male alleles, traces of purifying selection were found for a trehalose phosphatase, an exostosin and a WRKY transcription factor, with strikingly low polymorphism levels between distant geographic regions. Both diversity and network analysis revealed that H alleles are more closely related to M than to F alleles. CONCLUSIONS: Hermaphrodite alleles appear to derive from male alleles of wild grapevines, with successive recombination events allowing import of diversity from the X into the Y chromosomal region and slowing down the expansion of the region into a full heteromorphic chromosome. Our data are consistent with multiple domestication events and show traces of introgression from other Asian Vitis species into the cultivated grapevine gene pool.

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

The African cultivated rice (Oryza glaberrima) was domesticated in West Africa 3000 years ago. Although less cultivated than the Asian rice (O. sativa), O. glaberrima landraces often display interesting adaptation to rustic environment (e.g. drought). Here, using RNA-seq technology, we were able to compare more than 12,000 transcripts between 9 O. glaberrima, 10 wild O. barthii and one O. meridionalis individuals. With a synonymous nucleotide diversity πs = 0.0006 per site, O. glaberrima appears as the least genetically diverse crop grass ever documented. Using approximate Bayesian computation, we estimated that O. glaberrima experienced a severe bottleneck during domestication. This demographic scenario almost fully accounts for the pattern of genetic diversity across O. glaberrima genome as we detected very few outliers regions where positive selection may have further impacted genetic diversity. Moreover, the large excess of derived nonsynonymous substitution that we detected suggests that the O. glaberrima population suffered from the 'cost of domestication'. In addition, we used this genome-scale data set to demonstrate that (i) O. barthii genetic diversity is positively correlated with recombination rate and negatively with gene density, (ii) expression level is negatively correlated with evolutionary constraint, and (iii) one region on chromosome 5 (position 4-6 Mb) exhibits a clear signature of introgression with a yet unidentified Oryza species. This work represents the first genome-wide survey of the African rice genetic diversity and paves the way for further comparison between the African and the Asian rice, notably regarding the genetics underlying domestication traits.

Developing nuclear DNA phylogenetic markers in the angiosperm genus Leucadendron (Proteaceae): a next-generation sequencing transcriptomic approach.

Despite the recent advances in generating molecular data, reconstructing species-level phylogenies for non-models groups remains a challenge. The use of a number of independent genes is required to resolve phylogenetic relationships, especially for groups displaying low polymorphism. In such cases, low-copy nuclear exons and non-coding regions, such as 3' untranslated regions (3'-UTRs) or introns, constitute a potentially interesting source of nuclear DNA variation. Here, we present a methodology meant to identify new nuclear orthologous markers using both public-nucleotide databases and transcriptomic data generated for the group of interest by using next generation sequencing technology. To identify PCR primers for a non-model group, the genus Leucadendron (Proteaceae), we adopted a framework aimed at minimizing the probability of paralogy and maximizing polymorphism. We anchored when possible the right-hand primer into the 3'-UTR and the left-hand primer into the coding region. Seven new nuclear markers emerged from this search strategy, three of those included 3'-UTRs. We further compared the phylogenetic potential between our new markers and the ribosomal internal transcribed spacer region (ITS). The sequenced 3'-UTRs yielded higher polymorphism rates than the ITS region did. We did not find strong incongruences with the phylogenetic signal contained in the ITS region and the seven new designed markers but they strongly improved the phylogeny of the genus Leucadendron. Overall, this methodology is efficient in isolating orthologous loci and is valid for any non-model group given the availability of transcriptomic data.

Primer pairs, PCR conditions, and peptide nucleic acid clamps affect fungal diversity assessment from plant root tissues.

High-throughput sequencing has become a prominent tool to assess plant-associated microbial diversity. Still, some technical challenges remain in characterising these communities, notably due to plant and fungal DNA co-amplification. Fungal-specific primers, Peptide Nucleic Acid (PNA) clamps, or adjusting PCR conditions are approaches to limit plant DNA contamination. However, a systematic comparison of these factors and their interactions, which could limit plant DNA contamination in the study of plant mycobiota, is still lacking. Here, three primers targeting the ITS2 region were evaluated alone or in combination with PNA clamps both on nettle (Urtica dioica) root DNA and a mock community. PNA clamps did not improve the richness or diversity of the fungal communities but increased the number of fungal reads. Among the tested factors, the most significant was the primer pair. Specifically, the 5.8S-Fun/ITS4-Fun pair exhibited a higher OTU richness but fewer fungal reads. Our study demonstrates that the choice of primers is critical for limiting plant and fungal DNA co-amplification. PNA clamps increase the number of fungal reads when ITS2 is targeted but do not result in higher fungal diversity recovery at high sequencing depth. At lower read depths, PNA clamps might enhance microbial diversity quantification for primer pairs lacking fungal specificity.

The homomorphic self-incompatibility system in Oleaceae is controlled by a hemizygous genomic region expressing a gibberellin pathway gene.

In flowering plants, outcrossing is commonly ensured by self-incompatibility (SI) systems. These can be homomorphic (typically with many different allelic specificities) or can accompany flower heteromorphism (mostly with just two specificities and corresponding floral types). The SI system of the Oleaceae family is unusual, with the long-term maintenance of only two specificities but often without flower morphology differences. To elucidate the genomic architecture and molecular basis of this SI system, we obtained chromosome-scale genome assemblies of Phillyrea angustifolia individuals and related them to a genetic map. The S-locus region proved to have a segregating 543-kb indel unique to one specificity, suggesting a hemizygous region, as observed in all distylous systems so far studied at the genomic level. Only one of the predicted genes in this indel region is found in the olive tree, Olea europaea, genome, also within a segregating indel. We describe complete association between the presence/absence of this gene and the SI types determined for individuals of seven distantly related Oleaceae species. This gene is predicted to be involved in catabolism of the gibberellic acid (GA) hormone, and experimental manipulation of GA levels in developing buds modified the male and female SI responses of the two specificities in different ways. Our results provide a unique example of a homomorphic SI system, where a single conserved gibberellin-related gene in a hemizygous indel underlies the long-term maintenance of two groups of reproductive compatibility.

Disentangling homeologous contigs in allo-tetraploid assembly: application to durum wheat.

BACKGROUND: Using Next Generation Sequencing, SNP discovery is relatively easy on diploid species and still hampered in polyploid species by the confusion due to homeology. We develop HomeoSplitter; a fast and effective solution to split original contigs obtained by RNAseq into two homeologous sequences. It uses the differential expression of the two homeologous genes in the RNA. We verify that the new sequences are closer to the diploid progenitors of the allopolyploid species than the original contig. By remapping original reads on these new sequences, we also verify that the number of valuable detected SNPs has significantly increased. RESULTS: HomeoSplitter is a fast and effective solution to disentangle homeologous sequences based on a maximum likelihood optimization. On a benchmark set of 2,505 clusters containing homologous sequences of urartu, speltoides and durum, HomeoSplitter was efficient to build sequences closer to the diploid references and increased the number of valuable SNPs from 188 out of 1,360 SNPs detected when mapping the reads on the de novo durum assembly to 762 out of 1,620 SNPs when mapping on HomeoSplitter contigs. CONCLUSIONS: The HomeoSplitter program is freely available at http://bioweb.supagro.inra.fr/homeoSplitter/. This work provides a practical solution to the complex problem of disentangling homeologous transcripts in allo-tetraploids, which further allows an improved SNP detection.

Male-specific DNA markers provide genetic evidence of an XY chromosome system, a recombination arrest and allow the tracing of paternal lineages in date palm.

Whether sex chromosomes are differentiated is an important aspect of our knowledge of dioecious plants, such as date palm (Phoenix dactylifera). In this crop plant, the female individuals produce dates, and are thus the more valuable sex. However, there is no way to identify the sex of date palm plants before reproductive age, and the sex-determining mechanism is still unclear. To identify sex-linked microsatellite markers, we surveyed a set of 52 male and 55 female genotypes representing the geographical diversity of the species. We found three genetically linked loci that are heterozygous only in males. Male-specific alleles allowed us to identify the gender in 100% of individuals. These results confirm the existence of an XY chromosomal system with a nonrecombining XY-like region in the date palm genome. The distribution of Y haplotypes in western and eastern haplogroups allowed us to trace two male ancestral paternal lineages that account for all known Y diversity in date palm. The very low diversity associated with Y haplotypes is consistent with clonal paternal transmission of a nonrecombining male-determining region. Our results establish the date palm as a biological model with one of the most ancient sex chromosomes in flowering plants.

Loss of genetic diversity as a signature of apricot domestication and diffusion into the Mediterranean Basin.

BACKGROUND: Domestication generally implies a loss of diversity in crop species relative to their wild ancestors because of genetic drift through bottleneck effects. Compared to native Mediterranean fruit species like olive and grape, the loss of genetic diversity is expected to be more substantial for fruit species introduced into Mediterranean areas such as apricot (Prunus armeniaca L.), which was probably primarily domesticated in China. By comparing genetic diversity among regional apricot gene pools in several Mediterranean areas, we investigated the loss of genetic diversity associated with apricot selection and diffusion into the Mediterranean Basin. RESULTS: According to the geographic origin of apricots and using Bayesian clustering of genotypes, Mediterranean apricot (207 genotypes) was structured into three main gene pools: 'Irano-Caucasian', 'North Mediterranean Basin' and 'South Mediterranean Basin'. Among the 25 microsatellite markers used, only one displayed deviations from the frequencies expected under neutrality. Similar genetic diversity parameters were obtained within each of the three main clusters using both all SSR loci and only 24 SSR loci based on the assumption of neutrality. A significant loss of genetic diversity, as assessed by the allelic richness and private allelic richness, was revealed from the 'Irano-Caucasian' gene pool, considered as a secondary centre of diversification, to the northern and southwestern Mediterranean Basin. A substantial proportion of shared alleles was specifically detected when comparing gene pools from the 'North Mediterranean Basin' and 'South Mediterranean Basin' to the secondary centre of diversification. CONCLUSIONS: A marked domestication bottleneck was detected with microsatellite markers in the Mediterranean apricot material, depicting a global image of two diffusion routes from the 'Irano-Caucasian' gene pool: North Mediterranean and Southwest Mediterranean. This study generated genetic insight that will be useful for management of Mediterranean apricot germplasm as well as genetic selection programs related to adaptive traits.