Microsatellites uncover extraordinary diversity in native American land races and wild populations of cultivated sunflower

The contemporary oilseed sunflower (Helianthus annuus L.) gene pool is a product of multiple breeding and domestication bottlenecks. Despite substantial phenotypic diversity, modest differences in molecular genetic diversity have been uncovered in anciently and recently domesticated sunflowers. The paucity of molecular marker polymorphisms in early analyses led to the hypothesis of a single domestication origin. Phylogenetic analyses were performed on 47 domesticated and wild germplasm accessions using 122 microsatellite loci distributed throughout the sunflower genome. Extraordinary allelic diversity was found in the Native American land races and wild populations, and progressively less allelic diversity was found in germplasm produced by successive cycles of domestication and breeding. Of 1,341 microsatellite alleles, 489 were unique to land races, exotic domesticates and wild populations, whereas only 15 were unique to elite inbred lines. The number of taxon-specific alleles was 35-fold greater among wild populations (26.27) than elite inbred lines (0.75). Microsatellite genotyping uncovered the possibility of multiple domestication origins. Land races domesticated by Native Americans of the southwestern US (Hopi and Havasupai) formed a clade independent of land races domesticated by Native Americans of the Great Plains and eastern US (Arikara and Seneca). Predictably, domestication and breeding have ratcheted genetic diversity down in sunflower. The contemporary oilseed sunflower gene pool, while not imperiled, could profit from an infusion of novel alleles from the reservoir of latent genetic diversity present in wild populations and Native American land races.     

The USDA Barley Core Collection: Genetic Diversity,Population Structure,and Potential for Genome-Wide Association Studies

New sources of genetic diversity must be incorporated into plant breeding programs if they are to continue increasing grain yield and quality, and tolerance to abiotic and biotic stresses. Germplasm collections provide a source of genetic and phenotypic diversity, but characterization of these resources is required to increase their utility for breeding programs. We used a barley SNP iSelect platform with 7,842 SNPs to genotype 2,417 barley accessions sampled from the USDA National Small Grains Collection of 33,176 accessions. Most of the accessions in this core collection are categorized as landraces or cultivars/breeding lines and were obtained from more than 100 countries. Both STRUCTURE and principal component analysis identified five major subpopulations within the core collection, mainly differentiated by geographical origin and spike row number (an inflorescence architecture trait). Different patterns of linkage disequilibrium (LD) were found across the barley genome and many regions of high LD contained traits involved in domestication and breeding selection. The genotype data were used to define ‘mini-core’ sets of accessions capturing the majority of the allelic diversity present in the core collection. These ‘mini-core’ sets can be used for evaluating traits that are difficult or expensive to score. Genome-wide association studies (GWAS) of ‘hull cover’, ‘spike row number’, and ‘heading date’ demonstrate the utility of the core collection for locating genetic factors determining important phenotypes. The GWAS results were referenced to a new barley consensus map containing 5,665 SNPs. Our results demonstrate that GWAS and high-density SNP genotyping are effective tools for plant breeders interested in accessing genetic diversity in large germplasm collections.     

Genetic Variation and Breeding Signature in Mass Selection Lines of the Pacific Oyster (Crassostrea gigas) Assessed by SNP Markers

In breeding industries, a challenging problem is how to keep genetic diversity over generations. To investigate genetic variation and identify breeding signatures in mass selected lines of Pacific oyster (Crassostrea gigas), three sixth-generation selected lines and four wild populations were assessed using 103 single nucleotide polymorphism (SNP) markers. The genetic diversity data indicated that the selected lines exhibited a significant reduction in the observed heterozygosity and observed number of alleles per locus compared with the wild populations (P≤0.05), indicating the selected lines tended to lose genetic diversity contrasted with the wild populations. The unweighted pair-group method with arithmetic mean (UPGMA) analysis showed that the wild populations and selected lines were not separated into two groups. Using four outlier tests, a total of 17 loci were found under selection at two levels. The global outlier detection suggested that 4 common outlier loci were subject to selection using both the hierarchical island model and Bayesian likelihood approaches. At regional level, 3 SNPs were detected as outlier using at least two outlier tests and one outlier SNP (CgSNP309) was overlapped in the two wild-selected population comparisons. The candidate outlier SNPs provide valuable resources for future association studies in C. gigas.     

Resequencing Reveals Different Domestication Rate for BADH1 and BADH2 in Rice (Oryza sativa)

BADH1 and BADH2 are two homologous genes, encoding betaine aldehyde dehydrogenase in rice. In the present study, we scanned BADHs sequences of 295 rice cultivars, and 10 wild rice accessions to determine the polymorphisms, gene functions and domestication of these two genes. A total of 16 alleles for BADH1 and 10 alleles for BADH2 were detected in transcribed region of cultivars and wild species. Association study showed that BADH1 has significant correlation with salt tolerance in rice during germination stage, the SNP (T/A) in exon 4 is highly correlated with salt tolerance index (STI) (P<10⁻⁴). While, BADH2 was only responsible for rice fragrance, of which two BADH2 alleles (8 bp deletion in exon 7 and C/T SNP in exon 13) explain 97% of aroma variation in our germplasm. Theses indicate that there are no overlapping functions between the two homologous genes. In addition, a large LD block was detected in BADH2 region, however, there was no large LD blocks in a 4-Mb region of BADH1. We found that BADH2 region only showed significant bias in Tajima’s D value from the balance. Extended haplotype homozygosity study revealed fragrant accessions had a large LD block that extended around the mutation site (8 bp deletion in exon 7) of BADH2, while both of the BADH1 alleles (T/A in exon 4) did not show large extended LD block. All these results suggested that BADH2 was domesticated during rice evolution, while BADH1 was not selected by human beings.     

Wheat domestication: lessons for the future

Wheat was one of the first crops to be domesticated more than 10,000 years ago in the Middle East. Molecular genetics and archaeological data have allowed the reconstruction of plausible domestication scenarios leading to modern cultivars. For diploid einkorn and tetraploid durum wheat, a single domestication event has likely occurred in the Karacadag Mountains, Turkey. Following a cross between tetraploid durum and diploid T.?tauschii, the resultant hexaploid bread wheat was domesticated and disseminated around the Caucasian region. These polyploidisation events facilitated wheat domestication and created genetic bottlenecks, which excluded potentially adaptive alleles. With the urgent need to accelerate genetic progress to confront the challenges of climate change and sustainable agriculture, wild ancestors and old landraces represent a reservoir of underexploited genetic diversity that may be utilized through modern breeding methods. Understanding domestication processes may thus help identifying new strategies.     

SNP frequency and allelic haplotype structure of Beta vulgaris expressed genes

Based on EST sequences, fragments of 37 genes have been amplified and sequenced in two inbred lines of sugar beet. The rate of single nucleotide polymorphisms (SNP) corresponded to 1 every 130 bp, with an average (nucleotide diversity) value of 7.6×10^–3. When extrapolated to the whole sugar beet genome, randomly compared lines differ at 5.4×10⁶ SNPs in the genetic pool considered. In a wider search for SNP-related polymorphisms, 96 fragments of expressed genes were scanned with SSCP (single-strand conformation polymorphism) and heteroduplex (HA) analyses in 8 inbred lines. One SSCP or HA polymorphism was found every 1,470 bp of amplified DNA, corresponding to 5×10⁵ SSCP or HA loci in the whole genome. This frequency, 11 times lower than the SNP rate, was attributed to the high frequency of base pair substitution along the amplified fragment analysed electrophoretically. Therefore nucleotide variability was further studied by sequencing fragments of 10 genes in the same 8 lines. The results indicate that sugar beet alleles of expressed genes are very frequently organized as robust intragene haplotypes. In the 8 lines analysed, two haplotypes were identified for each of three gene fragments, three haplotypes for six gene fragments and four haplotypes for one gene fragment which is in good correspondence with the number of alleles detected by SSCP and HA analysis. In a cross between two lines, SSCP or HA alleles of expressed genes have 54% probability to be different.     

Development and Evaluation of SoySNP50K,a High-Density Genotyping Array for Soybean

The objective of this research was to identify single nucleotide polymorphisms (SNPs) and to develop an Illumina Infinium BeadChip that contained over 50,000 SNPs from soybean (Glycine max L. Merr.). A total of 498,921,777 reads 35–45bp in length were obtained from DNA sequence analysis of reduced representation libraries from several soybean accessions which included six cultivated and two wild soybean (G. soja Sieb. et Zucc.) genotypes. These reads were mapped to the soybean whole genome sequence and 209,903 SNPs were identified. After applying several filters, a total of 146,161 of the 209,903 SNPs were determined to be ideal candidates for Illumina Infinium II BeadChip design. To equalize the distance between selected SNPs, increase assay success rate, and minimize the number of SNPs with low minor allele frequency, an iteration algorithm based on a selection index was developed and used to select 60,800 SNPs for Infinium BeadChip design. Of the 60,800 SNPs, 50,701 were targeted to euchromatic regions and 10,000 to heterochromatic regions of the 20 soybean chromosomes. In addition, 99 SNPs were targeted to unanchored sequence scaffolds. Of the 60,800 SNPs, a total of 52,041 passed Illumina’s manufacturing phase to produce the SoySNP50K iSelect BeadChip. Validation of the SoySNP50K chip with 96 landrace genotypes, 96 elite cultivars and 96 wild soybean accessions showed that 47,337 SNPs were polymorphic and generated successful SNP allele calls. In addition, 40,841 of the 47,337 SNPs (86%) had minor allele frequencies ≥10% among the landraces, elite cultivars and the wild soybean accessions. A total of 620 and 42 candidate regions which may be associated with domestication and recent selection were identified, respectively. The SoySNP50K iSelect SNP beadchip will be a powerful tool for characterizing soybean genetic diversity and linkage disequilibrium, and for constructing high resolution linkage maps to improve the soybean whole genome sequence assembly.     

Genetic Characterization and Linkage Disequilibrium Estimation of a Global Maize Collection Using SNP Markers

A newly developed maize Illumina GoldenGate Assay with 1536 SNPs from 582 loci was used to genotype a highly diverse global maize collection of 632 inbred lines from temperate, tropical, and subtropical public breeding programs. A total of 1229 informative SNPs and 1749 haplotypes within 327 loci was used to estimate the genetic diversity, population structure, and familial relatedness. Population structure identified tropical and temperate subgroups, and complex familial relationships were identified within the global collection. Linkage disequilibrium (LD) was measured overall and within chromosomes, allelic frequency groups, subgroups related by geographic origin, and subgroups of different sample sizes. The LD decay distance differed among chromosomes and ranged between 1 to 10 kb. The LD distance increased with the increase of minor allelic frequency (MAF), and with smaller sample sizes, encouraging caution when using too few lines in a study. The LD decay distance was much higher in temperate than in tropical and subtropical lines, because tropical and subtropical lines are more diverse and contain more rare alleles than temperate lines. A core set of inbreds was defined based on haplotypes, and 60 lines capture 90% of the haplotype diversity of the entire panel. The defined core sets and the entire collection can be used widely for different research targets.     

Applications of single nucleotide polymorphisms in crop genetics

The discovery of single nucleotide polymorphisms (SNPs) and insertions/deletions, which are the basis of most differences between alleles, has been simplified by recent developments in sequencing technology. SNP discovery in many crop species, such as corn and soybean, is relatively straightforward because of their high level of intraspecific nucleotide diversity, and the availability of many gene and expressed sequence tag (EST) sequences. For these species, direct readout of SNP haplotypes is possible. Haplotype-based analysis is more informative than analysis based on individual SNPs, and has more power in analyzing association with phenotypes. The elite germplasm of some crops may have been subjected to bottlenecks relatively recently, increasing the amount of linkage disequilibrium (LD) present and facilitating the association of SNP haplotypes at candidate gene loci with phenotypes. Whole-genome scans may help identify genome regions that are associated with interesting phenotypes if sufficient LD is present. Technological improvements make the use of SNP and indel markers attractive for high-throughput use in marker-assisted breeding, EST mapping and the integration of genetic and physical maps.     

Linkage Disequilibrium and Demographic History of Wild and Domestic Canids

Assessing the extent of linkage disequilibrium (LD) in natural populations of a nonmodel species has been difficult due to the lack of available genomic markers. However, with advances in genotyping and genome sequencing, genomic characterization of natural populations has become feasible. Using sequence data and SNP genotypes, we measured LD and modeled the demographic history of wild canid populations and domestic dog breeds. In 11 gray wolf populations and one coyote population, we find that the extent of LD as measured by the distance at which r² = 0.2 extends <10 kb in outbred populations to >1.7 Mb in populations that have experienced significant founder events and bottlenecks. This large range in the extent of LD parallels that observed in 18 dog breeds where the r² value varies from ~20 kb to >5 Mb. Furthermore, in modeling demographic history under a composite-likelihood framework, we find that two of five wild canid populations exhibit evidence of a historical population contraction. Five domestic dog breeds display evidence for a minor population contraction during domestication and a more severe contraction during breed formation. Only a 5% reduction in nucleotide diversity was observed as a result of domestication, whereas the loss of nucleotide diversity with breed formation averaged 35%.     

Analysis of gene-derived SNP marker polymorphism in US wheat (Triticum aestivum L.) cultivars

In this study, we developed 359 detection primers for single nucleotide polymorphisms (SNPs) previously discovered within intron sequences of wheat genes and used them to evaluate SNP polymorphism in common wheat (Triticum aestivum L.). These SNPs showed an average polymorphism information content (PIC) of 0.18 among 20 US elite wheat cultivars, representing seven market classes. This value increased to 0.23 when SNPs were pre-selected for polymorphisms among a diverse set of 13 hexaploid wheat accessions (excluding synthetic wheats) used in the wheat SNP discovery project (). PIC values for SNP markers in the D genome were approximately half of those for the A and B genomes. D genome SNPs also showed a larger PIC reduction relative to the other genomes (P < 0.05) when US cultivars were compared with the more diverse set of 13 wheat accessions. Within those accessions, D genome SNPs show a higher proportion of alleles with low minor allele frequencies (<0.125) than found in the other two genomes. These data suggest that the reduction of PIC values in the D genome was caused by differential loss of low frequency alleles during the population size bottleneck that accompanied the development of modern commercial cultivars. Additional SNP discovery efforts targeted to the D genome in elite wheat germplasm will likely be required to offset the lower diversity of this genome. With increasing SNP discovery projects and the development of high-throughput SNP assay technologies, it is anticipated that SNP markers will play an increasingly important role in wheat genetics and breeding applications. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Comprehensive genotyping of the USA national maize inbred seed bank

Background

Genotyping by sequencing, a new low-cost, high-throughput sequencing technology was used to genotype 2,815 maize inbred accessions, preserved mostly at the National Plant Germplasm System in the USA. The collection includes inbred lines from breeding programs all over the world.

Results

The method produced 681,257 single-nucleotide polymorphism (SNP) markers distributed across the entire genome, with the ability to detect rare alleles at high confidence levels. More than half of the SNPs in the collection are rare. Although most rare alleles have been incorporated into public temperate breeding programs, only a modest amount of the available diversity is present in the commercial germplasm. Analysis of genetic distances shows population stratification, including a small number of large clusters centered on key lines. Nevertheless, an average fixation index of 0.06 indicates moderate differentiation between the three major maize subpopulations. Linkage disequilibrium (LD) decays very rapidly, but the extent of LD is highly dependent on the particular group of germplasm and region of the genome. The utility of these data for performing genome-wide association studies was tested with two simply inherited traits and one complex trait. We identified trait associations at SNPs very close to known candidate genes for kernel color, sweet corn, and flowering time; however, results suggest that more SNPs are needed to better explore the genetic architecture of complex traits.

Conclusions

The genotypic information described here allows this publicly available panel to be exploited by researchers facing the challenges of sustainable agriculture through better knowledge of the nature of genetic diversity.     

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.     

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.     

Genome-wide haplotype changes produced by artificial selection during modern rice breeding in Japan

During the last 90 years, the breeding of rice has delivered cultivars with improved agronomic and economic characteristics. Crossing of different lines and successive artificial selection of progeny based on their phenotypes have changed the chromosomal constitution of the ancestors of modern rice; however, the nature of these changes is unclear. The recent accumulation of data for genome-wide single-nucleotide polymorphisms (SNPs) in rice has allowed us to investigate the change in haplotype structure and composition. To assess the impact of these changes during modern breeding, we studied 177 Japanese rice accessions, which were categorized into three groups: landraces, improved cultivars developed from 1931 to 1974 (the early breeding phase), and improved cultivars developed from 1975 to 2005 (the late breeding phase). Phylogenetic tree and structure analysis indicated genetic differentiation between non-irrigated (upland) and irrigated (lowland) rice groups as well as genetic structuring within the irrigated rice group that corresponded to the existence of three subgroups. Pedigree analysis revealed that a limited number of landraces and cultivars was used for breeding at the beginning of the period of systematic breeding and that 11 landraces accounted for 70% of the ancestors of the modern improved cultivars. The values for linkage disequilibrium estimated from SNP alleles and the haplotype diversity determined from consecutive alleles in five-SNP windows indicated that haplotype blocks became less diverse over time as a result of the breeding process. A decrease in haplotype diversity, caused by a reduced number of polymorphisms in the haplotype blocks, was observed in several chromosomal regions. However, our results also indicate that new haplotype polymorphisms have been generated across the genome during the breeding process. These findings will facilitate our understanding of the association between particular haplotypes and desirable phenotypes in modern Japanese rice cultivars.     

Molecular characterization of global maize breeding germplasm based on genome-wide single nucleotide polymorphisms

Characterization of genetic diversity is of great value to assist breeders in parental line selection and breeding system design. We screened 770 maize inbred lines with 1,034 single nucleotide polymorphism (SNP) markers and identified 449 high-quality markers with no germplasm-specific biasing effects. Pairwise comparisons across three distinct sets of germplasm, CIMMYT (394), China (282), and Brazil (94), showed that the elite lines from these diverse breeding pools have been developed with only limited utilization of genetic diversity existing in the center of origin. Temperate and tropical/subtropical germplasm clearly clustered into two separate groups. The temperate germplasm could be further divided into six groups consistent with known heterotic patterns. The greatest genetic divergence was observed between temperate and tropical/subtropical lines, followed by the divergence between yellow and white kernel lines, whereas the least divergence was observed between dent and flint lines. Long-term selection for hybrid performance has contributed to significant allele differentiation between heterotic groups at 20% of the SNP loci. There appeared to be substantial levels of genetic variation between different breeding pools as revealed by missing and unique alleles. Two SNPs developed from the same candidate gene were associated with the divergence between two opposite Chinese heterotic groups. Associated allele frequency change at two SNPs and their allele missing in Brazilian germplasm indicated a linkage disequilibrium block of 142 kb. These results confirm the power of SNP markers for diversity analysis and provide a feasible approach to unique allele discovery and use in maize breeding programs.     

Population structure and genetic diversity in a commercial maize breeding program assessed with SSR and SNP markers

Information about the genetic diversity and population structure in elite breeding material is of fundamental importance for the improvement of crops. The objectives of our study were to (a) examine the population structure and the genetic diversity in elite maize germplasm based on simple sequence repeat (SSR) markers, (b) compare these results with those obtained from single nucleotide polymorphism (SNP) markers, and (c) compare the coancestry coefficient calculated from pedigree records with genetic distance estimates calculated from SSR and SNP markers. Our study was based on 1,537 elite maize inbred lines genotyped with 359 SSR and 8,244 SNP markers. The average number of alleles per locus, of group specific alleles, and the gene diversity (D) were higher for SSRs than for SNPs. Modified Roger’s distance (MRD) estimates and membership probabilities of the STRUCTURE matrices were higher for SSR than for SNP markers but the germplasm organization in four heterotic pools was consistent with STRUCTURE results based on SSRs and SNPs. MRD estimates calculated for the two marker systems were highly correlated (0.87). Our results suggested that the same conclusions regarding the structure and the diversity of heterotic pools could be drawn from both markers types. Furthermore, although our results suggested that the ratio of the number of SSRs and SNPs required to obtain MRD or D estimates with similar precision is not constant across the various precision levels, we propose that between 7 and 11 times more SNPs than SSRs should be used for analyzing population structure and genetic diversity.