Single nucleotide polymorphisms and linkage disequilibrium in sunflower

Genetic diversity in modern sunflower (Helianthus annuus L.) cultivars (elite oilseed inbred lines) has been shaped by domestication and breeding bottlenecks and wild and exotic allele introgression⁻the former narrowing and the latter broadening genetic diversity. To assess single nucleotide polymorphism (SNP) frequencies, nucleotide diversity, and linkage disequilibrium (LD) in modern cultivars, alleles were resequenced from 81 genic loci distributed throughout the sunflower genome. DNA polymorphisms were abundant; 1078 SNPs (1/45.7 bp) and 178 insertions-deletions (INDELs) (1/277.0 bp) were identified in 49.4 kbp of DNA/genotype. SNPs were twofold more frequent in noncoding (1/32.1 bp) than coding (1/62.8 bp) sequences. Nucleotide diversity was only slightly lower in inbred lines (θ = 0.0094) than wild populations (θ = 0.0128). Mean haplotype diversity was 0.74. When extraploted across the genome (~3500 Mbp), sunflower was predicted to harbor at least 76.4 million common SNPs among modern cultivar alleles. LD decayed more slowly in inbred lines than wild populations (mean LD declined to 0.32 by 5.5 kbp in the former, the maximum physical distance surveyed), a difference attributed to domestication and breeding bottlenecks. SNP frequencies and LD decay are sufficient in modern sunflower cultivars for very high-density genetic mapping and high-resolution association mapping.     

Carotenoid biosynthesis genes provide evidence of geographical subdivision and extensive linkage disequilibrium in the carrot

According to the history of the cultivated carrot, root colour can be considered as a structural factor of carrot germplasm. Therefore, molecular variations of carotenoid biosynthesis genes, these being involved in colour traits, represent a good putative source of polymorphism related to diversity structure. Seven candidate genes involved in the carotenoid biosynthesis pathway have been analysed from a sample of 48 individual plants, each one from a different cultivar of carrot (Daucus carota L. ssp. sativus). The cultivars were chosen to represent a large diversity and a wide range of root colour. A high single nucleotide polymorphism (SNP) frequency of 1 SNP per 22 bp (mean π _sil = 0.020) was found on average within these genes. The analysis of genetic structure from carotenoid biosynthesis gene sequences and 17 putatively neutral microsatellites showed moderate genetic differentiation between cultivars originating from the West and the East (F _ST = 0.072), this being consistent with breeding history, but not previously evidenced by molecular tools. Surprisingly, carotenoid biosynthesis genes did not exhibit decay of LD (mean r ²  = 0.635) within the 700–1,000 bp analysed, even though a fast decay level of LD is expected in outcrossing species. The high level of intralocus LD found for carotenoid biosynthesis genes implies that candidate-gene association mapping for carrot root colour should be useful to validate gene function, but may be unable to identify precisely the causative variations involved in trait determinism. Finally this study affords the first molecular evidence of a genetic structure in cultivated carrot germplasm related to phylogeography.     

Grinding up wheat: a massive loss of nucleotide diversity since domestication

Several demographic and selective events occurred during the domestication of wheat from the allotetraploid wild emmer (Triticum turgidum ssp. dicoccoides). Cultivated wheat has since been affected by other historical events. We analyzed nucleotide diversity at 21 loci in a sample of 101 individuals representing 4 taxa corresponding to representative steps in the recent evolution of wheat (wild, domesticated, cultivated durum, and bread wheats) to unravel the evolutionary history of cultivated wheats and to quantify its impact on genetic diversity. Sequence relationships are consistent with a single domestication event and identify 2 genetically different groups of bread wheat. The wild group is not highly polymorphic, with only 212 polymorphic sites among the 21,720 bp sequenced, and, during domestication, diversity was further reduced in cultivated forms--by 69% in bread wheat and 84% in durum wheat--with considerable differences between loci, some retaining no polymorphism at all. Coalescent simulations were performed and compared with our data to estimate the intensity of the bottlenecks associated with domestication and subsequent selection. Based on our 21-locus analysis, the average intensity of domestication bottleneck was estimated at about 3--giving a population size for the domesticated form about one third that of wild dicoccoides. The most severe bottleneck, with an intensity of about 6, occurred in the evolution of durum wheat. We investigated whether some of the genes departed from the empirical distribution of most loci, suggesting that they might have been selected during domestication or breeding. We detected a departure from the null model of demographic bottleneck for the hypothetical gene HgA. However, the atypical pattern of polymorphism at this locus might reveal selection on the linked locus Gsp1A, which may affect grain softness--an important trait for end-use quality in wheat.     

Nucleotide diversity and linkage disequilibrium in wild avocado (Persea americana Mill.)

Resequencing studies provide the ultimate resolution of genetic diversity because they identify all mutations in a gene that are present within the sampled individuals. We report a resequencing study of Persea americana, a subtropical tree species native to Meso- and Central America and the progenitor of cultivated avocado. The sample includes 21 wild accessions from Mexico, Costa Rica, Ecuador, and the Dominican Republic. Estimated levels of nucleotide polymorphism and linkage disequilibrium (LD) are obtained from fully resolved haplotype data from 4 nuclear loci that span 5960 nucleotide sites. Results show that, although avocado is a subtropical tree crop and a predominantly outcrossing plant, the overall level of genetic variation is not exceptionally high (nucleotide diversity at silent sites, pi(sil) = 0.0102) compared with available estimates from temperate plant species. Intralocus LD decays rapidly to half the initial value within about 1 kb. Estimates of recombination rate (based on the sequence data) show that the rate is not exceptionally high when compared with annual plants such as wild barley or maize. Interlocus LD is significant owing to substantial population structure induced by mixing of the 3 botanical races of avocado.     

Association Mapping and the Genomic Consequences of Selection in Sunflower

The combination of large-scale population genomic analyses and trait-based mapping approaches has the potential to provide novel insights into the evolutionary history and genome organization of crop plants. Here, we describe the detailed genotypic and phenotypic analysis of a sunflower (Helianthus annuus L.) association mapping population that captures nearly 90% of the allelic diversity present within the cultivated sunflower germplasm collection. We used these data to characterize overall patterns of genomic diversity and to perform association analyses on plant architecture (i.e., branching) and flowering time, successfully identifying numerous associations underlying these agronomically and evolutionarily important traits. Overall, we found variable levels of linkage disequilibrium (LD) across the genome. In general, islands of elevated LD correspond to genomic regions underlying traits that are known to have been targeted by selection during the evolution of cultivated sunflower. In many cases, these regions also showed significantly elevated levels of differentiation between the two major sunflower breeding groups, consistent with the occurrence of divergence due to strong selection. One of these regions, which harbors a major branching locus, spans a surprisingly long genetic interval (ca. 25 cM), indicating the occurrence of an extended selective sweep in an otherwise recombinogenic interval.     

Patterns of nucleotide diversity in wild and cultivated rice

There are few reports of the patterns of polymorphism in the non-coding regions of plant genomes. In this study, we explored nucleotide diversity and linkage disequilibrium (LD) in 47 non-coding regions on chromosome 4 of wild and cultivated rice. The cultivated rice retained about 70% of the diversity of wild rice, which was verified by coalescent simulations with one population bottleneck for 198 combinations of duration and population sizes. Multi-locus likelihood analysis showed that the severity of the bottleneck ranged from 2.25 to 3.33, with an average value of 2.70; i.e., the diversity found in the cultivated rice could be explained by a founding population of 2,700 individuals if the initial domestication event occurred over a period of 1,000 years. LD decreased more rapidly in wild rice than in cultivated rice within 10 kb, and the LD observed in cultivated rice was increased at 100–140 kb by comparison with wild rice. The patterns of LD indicated the possibility of a haplotype block in cultivated rice but not in wild rice.     

A comparison of sequence-based polymorphism and haplotype content in transcribed and anonymous regions of the barley genome.

Direct estimates of sequence diversity provides an abundant source of DNA polymorphisms based on single nucleotide polymorphisms (SNPs). The frequency and distribution of nucleotide diversity within 23 genes associated with grain germination in barley were determined in a sample of accessions representing European cultivars, landraces, and wild barley accessions from throughout the fertile crescent. The overall nucleotide diversity ranged from 0.0021 to 0.0189 with a single nucleotide change being detected every 78 bp and insertion-deletion events being observed every 680 bp. Within the cultivated (H. vulgare) genepool, a small number of haplotypes were detected, the total number of haplotypes observed in H. spontaneum was almost double that detected in H. vulgare (46 and 26, respectively). Distinct haplotypes were observed in the H. spontaneum and landrace genepools, which are highly divergent from H. vulgare. A comparison of SNP-based haplotype data with EST-derived SSRs and genomic SSRs revealed a similar trend of decreasing variability in the cultivated genepool. However, the number of unique alleles identified in the cultivated sample was much greater with genomic SSRs (18%) compared with only 2.1% for SNPs and 3.8% for EST-derived SSRs. The potential utility of SNPs and EST-derived SSRs for association mapping in barley is discussed.     

Multilocus analysis of nucleotide variation of Oryza sativa and its wild relatives: severe bottleneck during domestication of rice

Varying degrees of reduction of genetic diversity in crops relative to their wild progenitors occurred during the process of domestication. Such information, however, has not been available for the Asian cultivated rice (Oryza sativa) despite its importance as a staple food and a model organism. To reveal levels and patterns of nucleotide diversity and to elucidate the genetic relationship and demographic history of O. sativa and its close relatives (Oryza rufipogon and Oryza nivara), we investigated nucleotide diversity data from 10 unlinked nuclear loci in species-wide samples of these species. The results indicated that O. rufipogon and O. nivara possessed comparable levels of nucleotide variation ((sil) = 0.0077 approximately 0.0095) compared with the relatives of other crops. In contrast, nucleotide diversity of O. sativa was as low as (sil) = 0.0024 and even lower ((sil) = 0.0021 for indica and 0.0011 for japonica), if we consider the 2 subspecies separately. Overall, only 20-10% of the diversity in the wild species was retained in 2 subspecies of the cultivated rice (indica and japonica), respectively. Because statistic tests did not reject the assumption of neutrality for all 10 loci, we further used coalescent to simulate bottlenecks under various lengths and population sizes to better understand the domestication process. Consistent with the dramatic reduction in nucleotide diversity, we detected a severe domestication bottleneck and demonstrated that the sequence diversity currently found in the rice genome could be explained by a founding population of 1,500 individuals if the initial domestication event occurred over a 3,000-year period. Phylogenetic analyses revealed close genetic relationships and ambiguous species boundary of O. rufipogon and O. nivara, providing additional evidence to treat them as 2 ecotypes of a single species. Lowest linkage disequilibrium (LD) was found in the perennial O. rufipogon where the r(2) value dropped to a negligible level within 400 bp, and the highest in the japonica rice where LD extended to the entirely sequenced region ( approximately 900 bp), implying that LD mapping by genome scans may not be feasible in wild rice due to the high density of markers needed.     

A first insight into peach [Prunus persica (L.) Batsch] SNP variability

Three factors may have reduced the diversity at both individual gene and whole genome levels in cultivated peach: its self-compatible mating system, the narrow genetic basis of most commercial cultivars, and the recent strong selection towards agronomically interesting traits. Previous diversity analyses with markers such as simple sequence repeats (SSRs) have revealed low levels of genetic variability. Here, we sequenced 23 genome-wide distributed DNA fragments in 47 occidental peach varieties, also observing reduced variability levels. On average, there was one single nucleotide polymorphism (SNP) every 598 bp and one indel every 4,189 bp. As expected, variability was higher in non-coding than in coding regions (one SNP every 390 non-coding bp versus one in 1,850 bp in coding DNA). In general, SNPs were observed at relatively high frequency, mean minor allele frequency?=?0.225, meaning that a large proportion of the SNPs discovered by sequencing similar germplasm will be useful for other purposes, such as association mapping. The average heterozygosity of the varieties was 0.28, with a low correlation between SSR and SNP heterozygosity. The whole sequence of two candidate genes, a pectate lyase 1 candidate for fruit firmness (CGPAA2668) and a sucrose synthase 1 candidate for sugar content (CGPPB6189), in the 47 varieties revealed that they both may have suffered a process of balancing selection.     

Molecular characterization and evolutionary pattern of the 9- cis -epoxycarotenoid dioxygenase NCED1 gene in grapevine

This study provides the first analysis of the level and patterns of nucleotide polymorphism of the NCED1 gene in grapevine (Vitis vinifera L.). A total of 123 sequences of the gene were analyzed to give a sample of 50 wild accessions and 73 cultivars. A high single nucleotide polymorphism and haplotype diversity was revealed in the cultivars studied, especially Tunisian germplasms which present an important and diverse reservoir of genetic diversity for grape breeding and conservation. The haplotype distribution highlights two origins of the cultivars studied: one may be related to primary grapevine gene pool domestication while the second seems to be more recent. Thus, besides domestication, gene introgression has also played a role in shaping the current varietal landscape of grape cultivars. Higher nucleotide and haplotype polymorphism was recorded for cultivars. This was accompanied by a higher recombination rate in cultivated grapevines for this gene, a recent selective sweep in wild samples and a balancing selection in cultivars. The conservation of genetic diversity of the endangered wild germplasm is important to ensure that the wild population can be used in future breeding programs of the domesticated cultivars. The high number of alleles discovered can be used as a valuable source for association studies between allele frequencies and phenotypic variations in this gene. In addition to natural selection, molecular evidence shows that genetic variation in this locus appears to be shaped by a combination of mutation and recombination events.     

Molecular genetic diversity and association mapping of nut and kernel traits in Slovenian hazelnut (<Emphasis Type="Italic">Corylus avellana</Emphasis>) germplasm

European hazelnut (Corylus avellana L.), cultivated in several areas of the world including Europe, Anatolia, and the USA, is an economically important nut crop due to its high mineral, oleic acid, amino acid, and phenolic compound content and pleasant flavor. This study examined molecular genetic diversity and population structure of 54 wild accessions and 48 cultivars from the Slovenian national hazelnut collection using amplified fragment length polymorphism (AFLP) and simple sequence repeat (SSR) markers. Eleven AFLP primer combinations and 49 SSR markers yielded 532 and 504 polymorphic fragments, respectively. As expected for a wind-pollinated, self-incompatible species, levels of genetic diversity were high with cultivars and wild accessions having mean dissimilarity values of 0.50 and 0.60, respectively. In general, cultivars and wild accessions clustered separately in dendrogram, principal coordinate, and population structure analyses with regional clustering of the wild material. The accessions were also characterized for ten nut and seven kernel traits and some wild accessions were shown to have breeding potential. Morphological principal component analysis showed distinct clustering of cultivars and wild accessions. An association mapping panel composed of 64 hazelnut cultivars and wild accessions had considerable variation for the nut and kernel quality traits. Morphological and molecular data were associated to identify markers controlling the traits. In all, 49 SSR markers were significantly associated with nut and kernel traits [P < 0.0001 and LD value (r ²) = 0.15–0.50]. This work is the first use of association mapping in hazelnut and has identified molecular markers associated with important quality parameters in this important nut crop.     

Reduced representation genome sequencing reveals patterns of genetic diversity and selection in apple

Identifying DNA sequence variations is a fundamental step towards deciphering the genetic basis of traits of interest.Here,a total of 20 cultivated and 10 wild apples were genotyped using specific-locus amplified fragment sequencing,and 39,635 single nucleotide polymorphisms with no missing genotypes and evenly distributed along the genome were selected to investigate patterns of genome-wide genetic variations between cultivated and wild apples.Overall,wild apples displayed higher levels of genetic diversity than cultivated apples.Linkage disequilibrium(LD) decays were observed quite rapidly in cultivated and wild apples,with an r~2-value below 0.2 at 440 and 280 bp,respectively.Moreover,bidirectional gene flow and different distribution patterns of LD blocks were detected between domesticated and wild apples.Most LD blocks unique to cultivated apples were located within QTL regions controlling fruit quality,thus suggesting that fruit quality had probably undergone selection during apple domestication.The genome of the earliest cultivated apple in China,Nai,was highly similar to that of Malus sieversii,and contained a small portion of genetic material from other wild apple species.This suggested that introgression could have been an important driving force during initial domestication of apple.These findings will facilitate future breeding and genetic dissection of complex traits in apple.     

Assessment of linkage disequilibrium in potato genome with single nucleotide polymorphism markers

The extent of linkage disequilibrium (LD) is an important factor in designing association mapping experiments. Unlike other plant species that have been analyzed so far for the extent of LD, cultivated potato (Solanum tuberosum L.), an outcrossing species, is a highly heterozygous autotetraploid. The favored genotypes of modern cultivars are maintained by vegetative propagation through tubers. As a first step in the LD analysis, we surveyed both coding and noncoding regions of 66 DNA fragments from 47 accessions for single nucleotide polymorphism (SNP). In the process, we combined information from the potato SNP database with experimental SNP detection. The total length of all analyzed fragments was >25 kb, and the number of screened sequence bases reached almost 1.4 million. Average nucleotide polymorphism (=11.5x10(-3)) and diversity (pi=14.6x10(-3)) was high compared to the other plant species. The overall Tajima's D value (0.5) was not significant, but indicates a deficit of low-frequency alleles relative to expectation. To eliminate the possibility that an elevated D value occurs due to population subdivision, we assessed the population structure with probabilistic statistics. The analysis did not reveal any significant subdivision, indicating a relatively homogenous population structure. However, the analysis of individual fragments revealed the presence of subgroups in the fragment closely linked to the R1 resistance gene. Data pooled from all fragments show relatively fast decay of LD in the short range (r2=0.208 at 1 kb) but slow decay afterward (r2=0.137 at approximately 70 kb). The estimate from our data indicates that LD in potato declines below 0.10 at a distance of approximately 10 cM. We speculate that two conflicting factors play a vital role in shaping LD in potato: the outcrossing mating type and the very limited number of meiotic generations.     

Genetic diversity and population structure in cultivated sunflower and a comparison to its wild progenitor, <Emphasis Type="Italic">Helianthus annuus</Emphasis> L

Crop germplasm collections are valuable resources for ongoing plant breeding efforts. To fully utilize such collections, however, researchers need detailed information about the amount and distribution of genetic diversity present within collections. Here, we report the results of a population genetic analysis of the primary gene pool of sunflower (Helianthus annuus L.) based on a broad sampling of 433 cultivated accessions from North America and Europe, as well as a range-wide collection of 24 wild sunflower populations. Gene diversity across the cultivars was 0.47, as compared with 0.70 in the wilds, indicating that cultivated sunflower harbors roughly two-thirds of the total genetic diversity present in wild sunflower. Population structure analyses revealed that wild sunflower can be subdivided into four genetically distinct population clusters throughout its North American range, whereas the cultivated sunflower gene pool could be split into two main clusters separating restorer lines from the balance of the gene pool. Use of a maximum likelihood method to estimate the contribution of the wild gene pool to the cultivated sunflower germplasm revealed that the bulk of the cultivar diversity is derived from two wild sunflower population genetic clusters that are primarily composed of individuals from the east-central United States, the same general region in which sunflower domestication is believed to have occurred. We also identified a nested subset of accessions that capture as much of the allelic diversity present within the sampled cultivated sunflower germplasm collection as possible. At the high end, a core set of 288 captured nearly 90% of the alleles present in the full set of 433, whereas a core set of just 12 accessions was sufficient to capture nearly 50% of the total allelic diversity present within this sample of cultivated sunflower.     

Diversity,genetic mapping,and signatures of domestication in the carrot (Daucus carota L.) genome,as revealed by Diversity Arrays Technology (DArT) markers

Carrot is one of the most economically important vegetables worldwide, but genetic and genomic resources supporting carrot breeding remain limited. We developed a Diversity Arrays Technology (DArT) platform for wild and cultivated carrot and used it to investigate genetic diversity and to develop a saturated genetic linkage map of carrot. We analyzed a set of 900 DArT markers in a collection of plant materials comprising 94 cultivated and 65 wild carrot accessions. The accessions were attributed to three separate groups: wild, Eastern cultivated and Western cultivated. Twenty-seven markers showing signatures for selection were identified. They showed a directional shift in frequency from the wild to the cultivated, likely reflecting diversifying selection imposed in the course of domestication. A genetic linkage map constructed using 188 F2 plants comprised 431 markers with an average distance of 1.1 cM, divided into nine linkage groups. Using previously anchored single nucleotide polymorphisms, the linkage groups were physically attributed to the nine carrot chromosomes. A cluster of markers mapping to chromosome 8 showed significant segregation distortion. Two of the 27 DArT markers with signatures for selection were segregating in the mapping population and were localized on chromosomes 2 and 6. Chromosome 2 was previously shown to carry the Vrn1 gene governing the biennial growth habit essential for cultivated carrot. The results reported here provide background for further research on the history of carrot domestication and identify genomic regions potentially important for modern carrot breeding.     

Using multilocus sequence data to assess population structure, natural selection, and linkage disequilibrium in wild tomatoes

We employed a multilocus approach to examine the effects of population subdivision and natural selection on DNA polymorphism in 2 closely related wild tomato species (Solanum peruvianum and Solanum chilense), using sequence data for 8 nuclear loci from populations across much of the species' range. Both species exhibit substantial levels of nucleotide variation. The species-wide level of silent nucleotide diversity is 18% higher in S. peruvianum (pi(sil) approximately 2.50%) than in S. chilense (pi(sil) approximately 2.12%). One of the loci deviates from neutral expectations, showing a clinal pattern of nucleotide diversity and haplotype structure in S. chilense. This geographic pattern of variation is suggestive of an incomplete (ongoing) selective sweep, but neutral explanations cannot be entirely dismissed. Both wild tomato species exhibit moderate levels of population differentiation (average F(ST) approximately 0.20). Interestingly, the pooled samples (across different demes) exhibit more negative Tajima's D and Fu and Li's D values; this marked excess of low-frequency polymorphism can only be explained by population (or range) expansion and is unlikely to be due to population structure per se. We thus propose that population structure and population/range expansion are among the most important evolutionary forces shaping patterns of nucleotide diversity within and among demes in these wild tomatoes. Patterns of population differentiation may also be impacted by soil seed banks and historical associations mediated by climatic cycles. Intragenic linkage disequilibrium (LD) decays very rapidly with physical distance, suggesting high recombination rates and effective population sizes in both species. The rapid decline of LD seems very promising for future association studies with the purpose of mapping functional variation in wild tomatoes.     

Historical Divergence and Gene Flow in the Genus Zea

Gene flow plays a fundamental role in plant evolutionary history, yet its role in population divergence—and ultimately speciation—remains poorly understood. We investigated gene flow and the modalities of divergence in the domesticate Zea mays ssp. mays and three wild Zea taxa using sequence polymorphism data from 26 nuclear loci. We described diversity across loci and assessed evidence for adaptive and purifying selection at nonsynonymous sites. For each of three divergence events in the history of these taxa, we used approximate Bayesian simulation to estimate population sizes and divergence times and explicitly compare among alternative models of divergence. Our estimates of divergence times are surprisingly consistent with previous data from other markers and suggest rapid diversification of lineages within Zea in the last ~150,000 years. We found widespread evidence of historical gene flow, including evidence for divergence in the face of gene flow. We speculate that cultivated maize may serve as a bridge for gene flow among otherwise allopatric wild taxa.     

Evolution of genetic diversity during the domestication of common-bean (Phaseolus vulgaris L.)

M13 DNA fingerprinting was used to determine evolutionary changes that occurred in Latin American germ plasm and USA cultivars of commonbean (Phaseolus vulgaris L.) during domestication. Linkage mapping experiments showed that M13-related sequences in the common-bean genome were either located at the distal ends of linkage groups or that they were unlinked to each other or to any previously mapped markers. Levels of polymorphism observed by hybridization with M13 (1 probe-enzyme combination) were comparable to those observed by hybridization with single-copy random PstI genomic probes (36 enzyme-probe combinations) but were higher than those observed for isozymes (10 loci). Results indicated that the wild ancestor had diverged into two taxa, one distributed in Middle America (Mexico, Central America, and Colombia) and the other in the Andes (Peru and Argentina); they also suggested separate domestications in the two areas leading to two cultivated gene pools. Domestication in both areas led to pronounced reductions in diversity in cultivated descendants in Middle America and the Andes. The marked lack of polymorphism within commercial classes of USA cultivars suggests that the dispersal of cultivars from the centers of origin and subsequent breeding of improved cultivars led to high levels of genetic uniformity. To our knowledge, this is the first crop for which this reduction in diversity has been documented with a single type of marker in lineages that span the evolution between wild ancestor and advanced cultivars.