Genetic diversity, population structure and linkage disequilibrium in elite Chinese winter wheat investigated with SSR markers

Genetic diversity, population structure analysis and codon substitutions of Indicine Badri cattle using ddRAD sequencing

The present study was carried out on 96 animals representing three distinct colour variants of Badri cattle to investigate the genetic diversity, population structure and substitution mutations in the genetic codons due to single nucleotide variations. The DNA samples of 96 Badri cows were genotyped using a double digestion restriction associated DNA (ddRAD) sequencing approach. A standardized bioinformatics pipeline was employed to identify single nucleotide polymorphisms (SNPs), initially detecting 7,168,552 SNPs through alignment with the Bos indicus reference genome assembly. Subsequent stringent quality filtration yielded 65,483 high-confidence SNPs for downstream analysis. Genetic diversity analysis of the Badri cattle population resulted in average values of 0.145, 0.088, and 0.091 for Shannon's diversity Index (I), Simpson's Diversity (h), and Simpson's Unbiased Diversity (uh), respectively. Genetic similarities between the black and brown, black and grey, and brown and grey Badri variants were found to be 0.9972, 0.9980 and 0.9970, respectively. Tajima's D diversity value was observed to be significant and positive for 99.29% of high-confidence SNPs (65,483). STRUCTURE analysis showed admixture among the three Badri colour variants, suggesting a lack of genetic differentiation. Annotation of high-confidence SNPs regarding genetic codon changes indicated maximum substitutions in the GGC with GGT (22 occurrences), followed by AAC to AGC (20 occurrences), GAA to TAA (19 occurrences) and CAA to CAG (19 occurrences). The study concludes there are genetic similarities among colour variants, lack of rare alleles, balancing selection, sudden population contraction and genetic codon substitutions within the Badri cattle population. Insights derived from SNP data analysis hold potential significance for conservation initiatives and breed improvement programs for indicine cattle.

Genetic Diversity and Population Structure in Türkiye Bread Wheat Genotypes Revealed by Simple Sequence Repeats (SSR) Markers

Wheat genotypes should be improved through available germplasm genetic diversity to ensure food security. This study investigated the molecular diversity and population structure of a set of Türkiye bread wheat genotypes using 120 microsatellite markers. Based on the results, 651 polymorphic alleles were evaluated to determine genetic diversity and population structure. The number of alleles ranged from 2 to 19, with an average of 5.44 alleles per locus. Polymorphic information content (PIC) ranged from 0.031 to 0.915 with a mean of 0.43. In addition, the gene diversity index ranged from 0.03 to 0.92 with an average of 0.46. The expected heterozygosity ranged from 0.00 to 0.359 with a mean of 0.124. The unbiased expected heterozygosity ranged from 0.00 to 0.319 with an average of 0.112. The mean values of the number of effective alleles (Ne), genetic diversity of Nei (H) and Shannon's information index (I) were estimated at 1.190, 1.049 and 0.168, respectively. The highest genetic diversity (GD) was estimated between genotypes G1 and G27. In the UPGMA dendrogram, the 63 genotypes were grouped into three clusters. The three main coordinates were able to explain 12.64, 6.38 and 4.90% of genetic diversity, respectively. AMOVA revealed diversity within populations at 78% and between populations at 22%. The current populations were found to be highly structured. Model-based cluster analyses classified the 63 genotypes studied into three subpopulations. The values of F-statistic (Fst) for the identified subpopulations were 0.253, 0.330 and 0.244, respectively. In addition, the expected values of heterozygosity (He) for these sub-populations were recorded as 0.45, 0.46 and 0.44, respectively. Therefore, SSR markers can be useful not only in genetic diversity and association analysis of wheat but also in its germplasm for various agronomic traits or mechanisms of tolerance to environmental stresses.

Widely untargeted metabolomic profiling unearths metabolites and pathways involved in leaf senescence and N remobilization in spring-cultivated wheat under different N regimes

Progression of leaf senescence consists of both degenerative and nutrient recycling processes in crops including wheat. However, the levels of metabolites in flag leaves in spring-cultivated wheat, as well as biosynthetic pathways involved under different nitrogen fertilization regimes, are largely unknown. Therefore, the present study employed a widely untargeted metabolomic profiling strategy to identify metabolites and biosynthetic pathways that could be used in a wheat improvement program aimed at manipulating the rate and onset of senescence by handling spring wheat (Dingxi 38) flag leaves sampled from no-, low-, and high-nitrogen (N) conditions (designated Groups 1, 2, and 3, respectively) across three sampling times: anthesis, grain filling, and end grain filling stages. Through ultrahigh-performance liquid chromatography-tandem mass spectrometry, a total of 826 metabolites comprising 107 flavonoids, 51 phenol lipids, 37 fatty acyls, 37 organooxygen compounds, 31 steroids and steroid derivatives, 18 phenols, and several unknown compounds were detected. Upon the application of the stringent screening criteria for differentially accumulated metabolites (DAMs), 28 and 23 metabolites were differentially accumulated in Group 1_vs_Group 2 and Group 1_vs_Group 3, respectively. From these, 1-O-Caffeoylglucose, Rhoifolin, Eurycomalactone;Ingenol, 4-Methoxyphenyl beta-D-glucopyranoside, and Baldrinal were detected as core conserved DAMs among the three groups with all accumulated higher in Group 1 than in the other two groups. Kyoto Encyclopedia of Genes and Genomes pathway analysis revealed that tropane, piperidine, and pyridine alkaloid biosynthesis; acarbose and validamycin biosynthesis; lysine degradation; and biosynthesis of alkaloids derived from ornithine, lysine, and nicotinic acid pathways were the most significantly (p < 0.05) enriched in Group 1_vs_Group 2, while flavone and flavonol as well as anthocyanins biosynthetic pathways were the most significantly (p < 0.05) enriched in Group 1_vs_Group 3. The results from this study provide a foundation for the manipulation of the onset and rate of leaf senescence and N remobilization in wheat.

Genetic Relationships of 118 Castanea Specific Germplasms and Construction of Their Molecular ID Based on Morphological Characteristics and SSR Markers

To understand the genetic relationships of Castanea species, 16 phenotypic traits were measured, simple sequence repeat (SSR) markers were analyzed, and molecular identity cards (IDs) were constructed for 118 Castanea materials using fluorescent capillary electrophoresis. The coefficient of variation values of the 16 morphological traits of the test materials ranged from 11.11% to 60.38%. A total of 58 alleles were detected using six pairs of SSR core primers, with an average number of 9.7 alleles per locus. The average number of valid alleles per locus was 3.9419 and the proportion of valid alleles was 40.78%. A total of 105 genotypes were detected, and the number of genotypic species that could be amplified per primer pair ranged from 8 to 26. The mean value of the observed heterozygosity was 0.4986. The variation in the He, H, and PIC values was similar; the size of I value was approximately 2.21 times larger, and its mean number of variations was 0.7390, 0.7359, 0.6985, and 1.6015, respectively. The classification of 118 Castanea species was performed using three analytical methods: structure analysis, neighbor-joining (NJ) cluster analysis, and principal coordinate analysis (PCoA), and the results of the three methods were in high agreement. Six pairs of SSR core primers with high polymorphism and strong discriminatory properties were used to identify 118 Castanea plants, and a unique molecular ID card was constructed for each material. These results provide insight into the genetic diversity and population structure of Castanea plants and a theoretical basis for improving the phenomenon of mixed varieties and substandard plants in the Castanea plant market.

Genotyping-by-Sequencing Based Molecular Genetic Diversity of Pakistani Bread Wheat (Triticum aestivum L.) Accessions

Spring wheat (Triticum aestivum L.) is one of the most imperative staple food crops, with an annual production of 765 million tons globally to feed ∼40% world population. Genetic diversity in available germplasm is crucial for sustainable wheat improvement to ensure global food security. A diversity panel of 184 Pakistani wheat accessions was genotyped using 123,596 high-quality single nucleotide polymorphism (SNP) markers generated by genotyping-by-sequencing with 42% of the SNPs mapped on B, 36% on A, and 22% on D sub-genomes of wheat. Chromosome 2B contains the most SNPs (9,126), whereas 4D has the least (2,660) markers. The mean polymorphic information content, genetic diversity, and major allele frequency of the population were 0.157, 0.1844, and 0.87, respectively. Analysis of molecular variance revealed a higher genetic diversity (80%) within the sub-population than among the sub-populations (20%). The genome-wide linkage disequilibrium was 0.34 Mbp for the whole wheat genome. Among the three subgenomes, A has the highest LD decay value (0.29 Mbp), followed by B (0.2 Mbp) and D (0.07 Mbp) genomes, respectively. The results of population structure, principal coordinate analysis, phylogenetic tree, and kinship analysis also divided the whole population into three clusters comprising 31, 33, and 120 accessions in group 1, group 2, and group 3, respectively. All groups were dominated by the local wheat accessions. Estimation of genetic diversity will be a baseline for the selection of breeding parents for mutations and the genome-wide association and marker-assisted selection studies.

Genome-Wide Association Study of Root-Lesion Nematodes Pratylenchus Species and Crown Rot Fusarium culmorum in Bread Wheat

Triticum aestivum L., also known as common wheat, is affected by many biotic stresses. Root diseases are the most difficult to tackle due to the complexity of phenotypic evaluation and the lack of resistant sources compared to other biotic stress factors. Soil-borne pathogens such as the root-lesion nematodes caused by the Pratylenchus species and crown rot caused by various Fusarium species are major wheat root diseases, causing substantial yield losses globally. A set of 189 advanced spring bread wheat lines obtained from the International Maize and Wheat Improvement Center (CIMMYT) were genotyped with 4056 single nucleotide polymorphisms (SNP) markers and screened for root-lesion nematodes and crown rot resistance. Population structure revealed that the genotypes could be divided into five subpopulations. Genome-Wide Association Studies were carried out for both resistances to Pratylenchus and Fusarium species. Based on our results, 11 different SNPs on chromosomes 1A, 1B, 2A, 3A, 4A, 5B, and 5D were significantly associated with root-lesion nematode resistance. Seven markers demonstrated association with P. neglectus, while the remaining four were linked to P. thornei resistance. In the case of crown rot, eight different markers on chromosomes 1A, 2B, 3A, 4B, 5B, and 7D were associated with Fusarium crown rot resistance. Identification and screening of root diseases is a challenging task; therefore, the newly identified resistant sources/genotypes could be exploited by breeders to be incorporated in breeding programs. The use of the identified markers in marker-assisted selection could enhance the selection process and cultivar development with root-lesion nematode and crown rot resistance.

Discovery of miRNAs and Development of Heat-Responsive miRNA-SSR Markers for Characterization of Wheat Germplasm for Terminal Heat Tolerance Breeding

A large proportion of the Asian population fulfills their energy requirements from wheat (Triticum aestivum L.). Wheat quality and yield are critically affected by the terminal heat stress across the globe. It affects approximately 40% of the wheat-cultivating regions of the world. Therefore, there is a critical need to develop improved terminal heat-tolerant wheat varieties. Marker-assisted breeding with genic simple sequence repeats (SSR) markers have been used for developing terminal heat-tolerant wheat varieties; however, only few studies involved the use of microRNA (miRNA)-based SSR markers (miRNA-SSRs) in wheat, which were found as key players in various abiotic stresses. In the present study, we identified 104 heat-stress-responsive miRNAs reported in various crops. Out of these, 70 miRNA-SSR markers have been validated on a set of 20 terminal heat-tolerant and heat-susceptible wheat genotypes. Among these, only 19 miRNA-SSR markers were found to be polymorphic, which were further used to study the genetic diversity and population structure. The polymorphic miRNA-SSRs amplified 61 SSR loci with an average of 2.9 alleles per locus. The polymorphic information content (PIC) value of polymorphic miRNA-SSRs ranged from 0.10 to 0.87 with a mean value of 0.48. The dendrogram constructed using unweighted neighbor-joining method and population structure analysis clustered these 20 wheat genotypes into 3 clusters. The target genes of these miRNAs are involved either directly or indirectly in providing tolerance to heat stress. Furthermore, two polymorphic markers miR159c and miR165b were declared as very promising diagnostic markers, since these markers showed specific alleles and discriminated terminal heat-tolerant genotypes from the susceptible genotypes. Thus, these identified miRNA-SSR markers will prove useful in the characterization of wheat germplasm through the study of genetic diversity and population structural analysis and in wheat molecular breeding programs aimed at terminal heat tolerance of wheat varieties.

Genetic Diversity, Linkage Disequilibrium and Population Structure of Bulgarian Bread Wheat Assessed by Genome-Wide Distributed SNP Markers: From Old Germplasm to Semi-Dwarf Cultivars

Genetic diversity and population structure are key resources for breeding purposes and genetic studies of important agronomic traits in crops. In this study, we described SNP-based genetic diversity, linkage disequilibrium and population structure in a panel of 179 bread wheat advanced cultivars and old accessions from Bulgaria, using an optimized wheat 25K Infinium iSelect array. Out of 19,019 polymorphic SNPs, 17,968 had a known chromosome position on the A (41%), B (42%) and D (11%) genome, and 6% were not assigned to any chromosome. Homoeologous group 4, in particular chromosome 4D, was the least polymorphic. In the total population, the Nei’s gene diversity was within the range 0.1–0.5, and the polymorphism information content ranged from 0.1 to 0.4. Significant differences between the old and modern collections were revealed with respect to the linkage disequilibrium (LD): the average values for LD (r²), the percentage of the locus pairs in LD and the LD decay were 0.64, 16% and 3.3 for the old germplasm, and 0.43, 30% and 4.1 for the modern releases, respectively. Structure and k-means clustering algorithm divided the panel into three groups. The old accessions formed a distinct subpopulation. The cluster analysis further distinguished the modern releases according to the geographic region and genealogy. Gene exchange was evidenced mainly between the subpopulations of contemporary cultivars. The achieved understanding of the genetic diversity and structure of the Bulgarian wheat population and distinctiveness of the old germplasm could be of interest for breeders developing cultivars with improved characteristics. The obtained knowledge about SNP informativeness and the LD estimation are worthwhile for selecting markers and for considering the composition of a population in association mapping studies of traits of interest.

Population structure and genome-wide association studies in bread wheat for phosphorus efficiency traits using 35 K Wheat Breeder's Affymetrix array

Soil bioavailability of phosphorus (P) is a major concern for crop productivity worldwide. As phosphatic fertilizers are a non-renewable resource associated with economic and environmental issues so, the sustainable option is to develop P use efficient crop varieties. We phenotyped 82 diverse wheat (Triticum aestivum L.) accessions in soil and hydroponics at low and sufficient P. To identify the genic regions for P efficiency traits, the accessions were genotyped using the 35 K-SNP array and genome-wide association study (GWAS) was performed. The high-quality SNPs across the genomes were evenly distributed with polymorphic information content values varying between 0.090 and 0.375. Structure analysis revealed three subpopulations (C1, C2, C3) and the phenotypic responses of these subpopulations were assessed for P efficiency traits. The C2 subpopulation showed the highest genetic variance and heritability values for numerous agronomically important traits as well as strong correlation under both P levels in soil and hydroponics. GWAS revealed 78 marker-trait associations (MTAs) but only 35 MTAs passed Bonferroni Correction. A total of 297 candidate genes were identified for these MTAs and their annotation suggested their involvement in several biological process. Out of 35, nine (9) MTAs were controlling polygenic trait (two controlling four traits, one controlling three traits and six controlling two traits). These multi-trait MTAs (each controlling two or more than two correlated traits) could be utilized for improving bread wheat to tolerate low P stress through marker-assisted selection (MAS).

Genome-wide association analysis of Mexican bread wheat landraces for resistance to yellow and stem rust

Deploying under-utilized landraces in wheat breeding has been advocated to accelerate genetic gains in current era of genomics assisted breeding. Mexican bread wheat landraces (Creole wheats) represent an important resource for the discovery of novel alleles including disease resistance. A core set of 1,098 Mexican landraces was subjected to multi-location testing for rust diseases in India, Mexico and Kenya. The landrace core set showed a continuous variation for yellow (YR) and stem rust (SR) disease severity. Principal component analysis differentiated Mexican landraces into three groups based on their respective collection sites. Linkage disequilibrium (LD) decay varied from 10 to 32 Mb across chromosomes with an averge of 23Mb across whole genome. Genome-wide association analysis revealed marker-trait associations for YR resistance in India and Mexico as well as for SR resistance in Kenya. In addition, significant additive-additive interaction effects were observed for both YR and SR resistance including genomic regions on chromosomes 1BL and 3BS, which co-locate with pleiotropic genes Yr29/Lr46/Sr58/Pm39/Ltn2 and Sr2/Yr30/Lr27, respectively. Study reports novel genomic associations for YR (chromosomes 1AL, 2BS, and 3BL) and SR (chromosomes 2AL, 4DS, and 5DS). The novel findings in Creole wheat landraces can be efficiently utilized for the wheat genetic improvement.

Genome-Wide Association Study of Tan Spot Resistance in a Hexaploid Wheat Collection From Kazakhstan

Tan spot, caused by Pyrenophora tritici-repentis, is a serious foliar disease of wheat in Kazakhstan with reported yield losses as high as 50% during epidemic years. Here, we report the evaluation of a collection of 191 hexaploid spring and winter wheat lines for tan spot resistance and its underlying genetic architecture using genome-wide association study (GWAS). Our wheat collection comprised candidate varieties from Kazakhstan, Russia, and CIMMYT. It was genotyped using the DArTseq technology and phenotyped for resistance to tan spot at seedling and adult plant stages in Kazakhstan. DArTseq SNPs revealed high genetic diversity (average polymorphic information content = 0.33) in the panel and genome-wide linkage disequilibrium decay at 22 Mb (threshold r² = 0.1). Principal component analysis revealed a clear separation of Eurasian germplasm from CIMMYT and IWWIP lines. GWAS identified 34 marker-trait associations (MTA) for resistance to tan spot and the amount of phenotypic variation explained by these MTA ranged from 4% to 13.7%. Our results suggest the existence of novel valuable resistant alleles on chromosomes 3BS, and 5DL and 6AL for resistance to Race 1 and Race 5, respectively, in addition to known genes tsn1 and tsc2. On chromosome 6AL, a genomic region spanning 3 Mb was identified conferring resistance to both Race 1 and Race 5. Epistatic interaction of associated loci was revealed on chromosomes 1B, 5B, 7B, 5A, and 6A contributing to additional variation of 3.2–11.7%. Twenty-five lines with the best allele combinations of SNPs associated with resistance to both races have been identified as candidates for future variety release and breeding. The results of the present study will be further validated in other independent genetic backgrounds to be able to use markers in breeding.

Population structure of Nepali spring wheat (Triticum aestivum L.) germplasm

BACKGROUND

Appropriate information about genetic diversity and population structure of germplasm improves the efficiency of plant breeding. The low productivity of Nepali bread wheat (Triticum aestivum L.) is a major concern particularly since Nepal is ranked the 4th most vulnerable nation globally to climate change. The genetic diversity and population structure of Nepali spring wheat have not been reported. This study aims to improve the exploitation of more diverse and under-utilized genetic resources to contribute to current and future breeding efforts for global food security.

RESULTS

We used genotyping-by-sequencing (GBS) to characterize a panel of 318 spring wheat accessions from Nepal including 166 landraces, 115 CIMMYT advanced lines, and 34 Nepali released varieties. We identified 95 K high-quality SNPs. The greatest genetic diversity was observed among the landraces, followed by CIMMYT lines, and released varieties. Though we expected only 3 groupings corresponding to these 3 seed origins, the population structure revealed two large, distinct subpopulations along with two smaller and scattered subpopulations in between, with significant admixture. This result was confirmed by principal component analysis (PCA) and UPGMA distance-based clustering. The pattern of LD decay differed between subpopulations, ranging from 60 to 150 Kb. We discuss the possibility that germplasm explorations during the 1970s-1990s may have mistakenly collected exotic germplasm instead of local landraces and/or collected materials that had already cross-hybridized since exotic germplasm was introduced starting in the 1950s.

CONCLUSION

We suggest that only a subset of wheat "landraces" in Nepal are authentic which this study has identified. Targeting these authentic landraces may accelerate local breeding programs to improve the food security of this climate-vulnerable nation. Overall, this study provides a novel understanding of the genetic diversity of wheat in Nepal and this may contribute to global wheat breeding initiatives.

Phenotypic and genetic variability of a tetraploid wheat collection grown in Kazakhstan

New cultivars adapted to major durum wheat growing environments are essential for the cultivation of this crop. The development of new cultivars has required the availability of diverse genetic material and their extensive field trials. In this work, a collection of tetraploid wheat consisting of 85 accessions was tested in the field conditions of Almaty region during 2018 and 2019. The accessions were ranged according to nine agronomic traits studied, and accessions with the highest yield performance for Almaty region of Kazakhstan were revealed. The ANOVA suggested that the performance of agronomic traits were influenced both by Environment and Genotype. Also, the collection was analyzed using seven SSR (simple sequence repeats) markers. From 3 to 6 alleles per locus were revealed, with an average of 4.6, while the effective number of alleles was 2.8. Nei's genetic diversity was in the range of 0.45-0.69. The results showed high values of polymorphism index content (PIC) in the range of 0.46-0.70, with an average of 0.62, suggesting that 6 out of 7 SSRs were highly informative (PIC > 0.5). Phylogenetic analysis of the collection has allowed the separation of accessions into six clusters. The local accessions were presented in all six clusters with the majority of them grouped in the first three clusters designated as A, B, and C, respectively. The relations between SSR markers and agronomic traits in the collection were studied. The results can be efficiently used for the enhancement of local breeding projects for the improvement of yield productivity in durum wheat.

Genome-wide association analysis of stem water-soluble carbohydrate content in bread wheat

GWAS identified 36 potentially new loci for wheat stem water-soluble carbohydrate (WSC) contents and 13 pleiotropic loci affecting WSC and thousand-kernel weight. Five KASP markers were developed and validated. Water-soluble carbohydrates (WSC) reserved in stems contribute significantly to grain yield (GY) in wheat. However, knowledge of the genetic architecture underlying stem WSC content (SWSCC) is limited. In the present study, 166 diverse wheat accessions from the Yellow and Huai Valleys Winter Wheat Zone of China and five other countries were grown in four well-watered environments. SWSCC at 10 days post-anthesis (10DPA), 20DPA and 30DPA, referred as WSC10, WSC20 and WSC30, respectively, and thousand-kernel weight (TKW) were assessed. Correlation analysis showed that TKW was significantly and positively correlated with WSC10 and WSC20. Genome-wide association study was performed on SWSCC and TKW with 373,106 markers from the wheat 660 K and 90 K SNP arrays. Totally, 62 stable loci were detected for SWSCC, with 36, 24 and 19 loci for WSC10, WSC20 and WSC30, respectively; among these, 36 are potentially new, 16 affected SWSCC at two or three time-points, and 13 showed pleiotropic effects on both SWSCC and TKW. Linear regression showed clear cumulative effects of favorable alleles for increasing SWSCC and TKW. Genetic gain analyses indicated that pyramiding favorable alleles of SWSCC had simultaneously improved TKW. Kompetitive allele-specific PCR markers for five pleiotropic loci associated with both SWSCC and TKW were developed and validated. This study provided a genome-wide landscape of the genetic architecture of SWSCC, gave a perspective for understanding the relationship between WSC and GY and explored the theoretical basis for co-improvement of WSC and GY. It also provided valuable loci and markers for future breeding.

Genetic Diversity and Population Structure of Asian and European Common Wheat Accessions Based on Genotyping-By-Sequencing

Obtaining information on the genetic diversity and population structure of germplasm facilitates its use in wheat breeding programs. Recently, with the development of next-generation sequencing technology, genotyping-by-sequencing (GBS) has been used as a high-throughput and cost-effective molecular tool for examination of the genetic diversity of wheat breeding lines. In this study, GBS was used to characterize a population of 180 accessions of common wheat originating from Asia and Europe between the latitudes 30° and 45°N. In total, 24,767 high-quality single-nucleotide polymorphism (SNP) markers were used for analysis of genetic diversity and population structure. The B genome contained the highest number of SNPs, followed by the A and D genomes. The polymorphism information content was in the range of 0.1 to 0.4, with a mean of 0.26. The distribution of SNPs markers on the 21 chromosomes ranged from 243 on chromosome 4D to 2,337 on chromosome 3B. Structure and cluster analyses divided the panel of accessions into two subgroups (G1 and G2). G1 principally consisted of European and partial Asian accessions, and G2 comprised mainly accessions from the Middle East and partial Asia. Molecular analysis of variance showed that the genetic variation was greater within groups (99%) than between groups (1%). Comparison of the two subgroups indicated that G1 and G2 contained a high level of genetic diversity. The genetic diversity of G2 was slightly higher as indicated by the observed heterozygosity (H o) = 0.23, and unbiased diversity index (uh) = 0.34. The present results will not only help breeders to understand the genetic diversity of wheat germplasm on the Eurasian continent between the latitudes of 30° and 45°N, but also provide valuable information for wheat genetic improvement through introgression of novel genetic variation in this region.

Selection criteria for high-yielding and early-flowering bread wheat hybrids under heat stress

Hybrid performance during wheat breeding can be improved by analyzing genetic distance (GD) among wheat genotypes and determining its correlation with heterosis. This study evaluated the GD between 16 wheat genotypes by using 60 simple sequence repeat (SSR) markers to classify them according to their relationships and select those with greater genetic diversity, evaluate the correlation of the SSR marker distance with heterotic performance and specific combining ability (SCA) for heat stress tolerance, and identify traits that most influence grain yield (GY). Eight parental genotypes with greater genetic diversity and their 28 F1 hybrids generated using diallel crossing were evaluated for 12 measured traits in two seasons. The GD varied from 0.235 to 0.911 across the 16 genotypes. Cluster analysis based on the GD estimated using SSRs classified the genotypes into three major groups and six sub-groups, almost consistent with the results of principal coordinate analysis. The combined data indicated that five hybrids showed 20% greater yield than mid-parent or better-parent. Two hybrids (P2 × P4) and (P2 × P5), which showed the highest performance of days to heading (DH), grain filling duration (GFD), and GY, and had large genetic diversity among themselves (0.883 and 0.911, respectively), were deemed as promising heat-tolerant hybrids. They showed the best mid-parent heterosis and better-parent heterosis (BPH) for DH (-11.57 and -7.65%; -13.39 and -8.36%, respectively), GFD (12.74 and 12.17%; 12.09 and 10.59%, respectively), and GY (36.04 and 20.04%; 44.06 and 37.73%, respectively). Correlation between GD and each of BPH and SCA effects based on SSR markers was significantly positive for GFD, hundred kernel weight, number of kernels per spike, harvest index, GY, and grain filling rate and was significantly negative for DH. These correlations indicate that the performance of wheat hybrids with high GY and earliness could be predicted by determining the GD of the parents by using SSR markers. Multivariate analysis (stepwise regression and path coefficient) suggested that GFD, hundred kernel weight, days to maturity, and number of kernels per spike had the highest influence on GY.

Longin CF, Ganal M, Himmelbach A, Reif JC, Schulthess AW. Exome association analysis sheds light onto leaf rust (Puccinia triticina) resistance genes currently used in wheat breeding (Triticum aestivum L.)

Resistance breeding is crucial for a sustainable control of leaf rust (Puccinia triticina) in wheat (Triticum aestivum L.) while directly targeting functional variants is the Holy Grail for efficient marker-assisted selection and map-based cloning. We assessed the limits and prospects of exome association analysis for severity of leaf rust in a large hybrid wheat population of 1574 single-crosses plus their 133 parents. After imputation and quality control, exome sequencing revealed 202 875 single-nucleotide polymorphisms (SNPs) covering 19.7% of the high-confidence annotated gene space. We performed intensive data mining and found significant associations for 2171 SNPs corresponding to 50 different loci. Some of these associations mapped in the proximity of the already known resistance genes Lr21, Lr34-B, Lr1 and Lr10, while other associated genomic regions, such as those on chromosomes 1A and 3D, harboured several annotated genes putatively involved in resistance. Validation with an independent population helped to narrow down the list of putative resistance genes that should be targeted by fine-mapping. We expect that the proposed strategy of intensive data mining coupled with validation will significantly influence research in plant genetics and breeding.

Exploring the genetic base of the soybean germplasm from Africa, America and Asia as well as mining of beneficial allele for flowering and seed weight

Genetic diversity is the foundation for any breeding program. The present study analyzed the genetic base of 163 soybean genotypes from three continents viz. Africa, America and Asia using 68 trait-linked simple sequence repeats (SSR) markers. The average number of alleles among the germplasm from the three continents followed the trend as Asia (9) > America (8) > Africa (7). Similar trends were observed for gene diversity (0.76 > 0.74 > 0.71) and polymorphism information content (PIC) (0.73 > 0.71 > 0.68). These findings revealed that soybean germplasm from Asia has wider genetic base followed by America, and least in Africa. The 163 genotypes were grouped into 4 clusters by phylogenetic analysis, whereas model-based population structure analysis also divided them into 4 subpopulations comprising 80.61% pure lines and 19.39% admixtures. The genotypes from Africa were easily distinguished from those of other two continents using phylogenetic analysis, indicating important role of geographyical differentiation for this genetic variability. Our results indicated that soybean germplasm has moved from Asia to America, and from America to Africa. Analysis of molecular variance (AMOVA) showed 8.41% variation among the four subpopulations, whereas 63.12% and 28.47% variation existed among and within individuals in the four subpopulations, respectively. Based on the association mapping, a total of 21 SSR markers showed significant association with days to flowering (DoF) and 100-seed weight (HSW). Two markers Satt365 and Satt581 on chromosome 6 and 10, respectively, showed pleiotropic effect or linkage on both traits. Genotype A50 (Gakuran Daizu/PI 506679) from Japan has 8 out of the 13 beneficial alleles for increased HSW. The diverse genotypes, polymorphic SSR markers and desirable alleles identified for DoF and HSW will be used in future breeding programs to improve reproductive, yield and quality traits.

Molecular Markers Associated with Agro-Physiological Traits under Terminal Drought Conditions in Bread Wheat

Terminal drought stress poses a big challenge to sustain wheat grain production in rain-fed environments. This study aimed to utilize the genetically diverse pre-breeding lines for identification of genomic regions associated with agro-physiological traits at terminal stage drought stress in wheat. A total of 339 pre-breeding lines panel derived from three-way crosses of 'exotics × elite × elite' lines were evaluated in field conditions at Obregon, Mexico for two years under well irrigated as well as drought stress environments. Drought stress was imposed at flowering by skipping the irrigations at pre and post anthesis stage. Results revealed that drought significantly reduced grain yield (Y), spike length (SL), number of grains spikes-1 (NGS) and thousand kernel weight (TKW), while kernel abortion (KA) was increased. Population structure analysis in this panel uncovered three sub-populations. Genome wide linkage disequilibrium (LD) decay was observed at 2.5 centimorgan (cM). The haplotypes-based genome wide association study (GWAS) identified significant associations of Y, SL, and TKW on three chromosomes; 4A (HB10.7), 2D (HB6.10) and 3B (HB8.12), respectively. Likewise, associations on chromosomes 6B (HB17.1) and 3A (HB7.11) were found for NGS while on chromosome 3A (HB7.12) for KA. The genomic analysis information generated in the study can be efficiently utilized to improve Y and/or related parameters under terminal stage drought stress through marker-assisted breeding.

Morphological and Genetic Diversity within Salt Tolerance Detection in Eighteen Wheat Genotypes

Salinity is a major obstacle to wheat production worldwide. Salt-affected soils could be used by improving salt-tolerant genotypes depending upon the genetic variation and salt stress response of adapted and donor wheat germplasm. We used a comprehensive set of morpho-physiological and biochemical parameters and simple sequence repeat (SSR) marker technique with multivariate analysis to accurately demonstrate the phenotypic and genetic variation of 18 wheat genotypes under salinity stress. All genotypes were evaluated without NaCl as a control and with 150 mM NaCl, until the onset of symptoms of death in the sensitive plant (after 43 days of salinity treatment). The results showed that the relative change of the genetic variation was high for all parameters, heritability (> 60%), and genetic gain (> 20%). Stepwise regression analysis, noting the importance of the root dry matter, relative turgidity, and their respective contributions to the shoot dry matter, indicated their relevance in improving and evaluating the salt-tolerant genotypes of breeding programs. The relative change of the genotypes in terms of the relative turgidity and shoot dry matter during salt stress was verified using clustering methods. For cluster analysis, the genotypes were classified into three groups: tolerant, intermediate, and sensitive, representing five, six, and seven genotypes, respectively. The morphological and genetic distances were significantly correlated based on the Mantel test. Of the 23 SSR markers that showed polymorphism, 17 were associated with almost all examined parameters. Therefore, based on the observed molecular marker-phenotypic trait association, the markers were highly useful in detecting tolerant and sensitive genotypes. Thus, it considers a helpful tool for salt tolerance through marker-assisted selection.

Whole-genome diversity, population structure and linkage disequilibrium analysis of globally diverse wheat genotypes using genotyping-by-sequencing DArTseq platform

In this study, 129 wheat genotypes from globally diverse origins were genotyped using DArTseq (SilicoDArT and SNP) markers. After filtering markers for quality-filtering, 14,270 SilicoDArTs and 6484 SNPs were retained and used for genetic diversity, population structure and linkage disequilibrium analyses. The highest number of SilicoDArT and SNP markers mapped on genome A and B compared to genome D. In both marker types, polymorphism information content (PIC) values ranged from 0.1 to 0.5, while > 0.80% of SilicoDArTs and > 0.44% SNPs showed PIC value more than median (0.25%). Un-weighted Neighbor Joining cluster analysis and Bayesian-based model population structure grouped wheat genotypes into three and four clusters, respectively. Principal component analysis and discriminant analysis of principal component results showed highly match with cluster and population structure analysis. Linkage disequilibrium (LD) was more extensive in both marker types, while graphical display of LD decay for both marker types showed that LD declined in the region close to 15 kbp, where r ²-values corresponded to r ² = 0.16. Overall, our genetic diversity analysis showed high level of variation in studied wheat genotypes, even though there was no relationship between wheat grouping and origins. This might be attributed to admixture level that occurred during long-term natural selection of wheat genotypes in different parts of the world. Highly diverse wheat genotypes used in this study may possess unique genes and are useful sources in breeding programs to improve grain yield and quality.

GWAS to Identify Genetic Loci for Resistance to Yellow Rust in Wheat Pre-Breeding Lines Derived From Diverse Exotic Crosses

Yellow rust (YR) or stripe rust, caused by Puccinia striformis f. sp tritici Eriks (Pst), is a major challenge to resistance breeding in wheat. A genome wide association study (GWAS) was performed using 22,415 single nucleotide polymorphism (SNP) markers and 591 haplotypes to identify genomic regions associated with resistance to YR in a subset panel of 419 pre-breeding lines (PBLs) developed at International Center for Maize and Wheat Improvement (CIMMYT). The 419 PBLs were derived from an initial set of 984 PBLs generated by a three-way crossing scheme (exotic/elite1//elite2) among 25 best elites and 244 exotics (synthetics, landraces) from CIMMYT's germplasm bank. For the study, 419 PBLs were characterized with 22,415 high-quality DArTseq-SNPs and phenotyped for severity of YR disease at five locations in Mexico. A population structure was evident in the panel with three distinct subpopulations, and a genome-wide linkage disequilibrium (LD) decay of 2.5 cM was obtained. Across all five locations, 14 SNPs and 7 haplotype blocks were significantly (P < 0.001) associated with the disease severity explaining 6.0 to 14.1% and 7.9 to 19.9% of variation, respectively. Based on average LD decay of 2.5 cM, identified 14 SNP-trait associations were delimited to seven quantitative trait loci in total. Seven SNPs were part of the two haplotype blocks on chromosome 2A identified in haplotypes-based GWAS. In silico analysis of the identified SNPs showed hits with interesting candidate genes, which are related to pathogenic process or known to regulate induction of genes related to pathogenesis such as those coding for glunolactone oxidase, quinate O-hydroxycinnamoyl transferase, or two-component histidine kinase. The two-component histidine kinase, for example, acts as a sensor in the perception of phytohormones ethylene and cytokinin. Ethylene plays a very important role in regulation of multiple metabolic processes of plants, including induction of defense mechanisms mediated by jasmonate. The SNPs linked to the promising genes identified in the study can be used for marker-assisted selection.

Genome-wide variation patterns between landraces and cultivars uncover divergent selection during modern wheat breeding

Genetic diversity, population structure, LD decay, and selective sweeps in 687 wheat accessions were analyzed, providing relevant guidelines to facilitate the use of the germplasm in wheat breeding. Common wheat (Triticum aestivum L.) is one of the most widely grown crops in the world. Landraces were subjected to strong human-mediated selection in developing high-yielding, good quality, and widely adapted cultivars. To investigate the genome-wide patterns of allelic variation, population structure and patterns of selective sweeps during modern wheat breeding, we tested 687 wheat accessions, including landraces (148) and cultivars (539) mainly from China and Pakistan in a wheat 90 K single nucleotide polymorphism array. Population structure analysis revealed that cultivars and landraces from China and Pakistan comprised three relatively independent genetic clusters. Cultivars displayed lower nucleotide diversity and a wider average LD decay across whole genome, indicating allelic erosion and a diversity bottleneck due to the modern breeding. Analysis of genetic differentiation between landraces and cultivars from China and Pakistan identified allelic variants subjected to selection during modern breeding. In total, 477 unique genome regions showed signatures of selection, where 109 were identified in both China and Pakistan germplasm. The majority of genomic regions were located in the B genome (225), followed by the A genome (175), and only 77 regions were located in the D genome. EigenGWAS was further used to identify key selection loci in modern wheat cultivars from China and Pakistan by comparing with global winter wheat and spring wheat diversity panels, respectively. A few known functional genes or loci found within these genome regions corresponded to known phenotypes for disease resistance, vernalization, quality, adaptability and yield-related traits. This study uncovered molecular footprints of modern wheat breeding and explained the genetic basis of polygenic adaptation in wheat. The results will be useful for understanding targets of modern wheat breeding, and in devising future breeding strategies to target beneficial alleles currently not pursued.

From landraces to improved cultivars: Assessment of genetic diversity and population structure of Mediterranean wheat using SNP markers

Assessment of genetic diversity and population structure in crops is essential for breeding and germplasm conservation. A collection of 354 bread wheat genotypes, including Mediterranean landraces and modern cultivars representative of the ones most widely grown in the Mediterranean Basin, were characterized with 11196 single nucleotide polymorphism (SNP) markers. Total genetic diversity (H_T) and polymorphic information content (PIC) were 0.36 and 0.30 respectively for both landraces and modern cultivars. Linkage disequilibrium for the modern cultivars was higher than for the landraces (0.18 and 0.12, respectively). Analysis of the genetic structure showed a clear geographical pattern for the landraces, which were clustered into three subpopulations (SPs) representing the western, northern and eastern Mediterranean, whereas the modern cultivars were structured according to the breeding programmes that developed them: CIMMYT/ICARDA, France/Italy, and Balkan/eastern European countries. The modern cultivars showed higher genetic differentiation (G_ST) and lower gene flow (0.1673 and 2.49, respectively) than the landraces (0.1198 and 3.67, respectively), indicating a better distinction between subpopulations. The maximum gene flow was observed between landraces from the northern Mediterranean SPs and the modern cultivars released mainly by French and Italian breeding programmes.

Genetic Structure Analysis of a Collection of Tunisian Durum Wheat Germplasm

The Tunisian durum wheat germplasm includes modern cultivars and traditional varieties that are still cultivated in areas where elite cultivars or intensive cultivation systems are not suitable. Within the frame of a collection program of the National Gene Bank of Tunisia (NGBT), durum wheat germplasm was collected from different Tunisian agro-ecological zones. The collected samples were studied using simple sequence repeats (SSRs) markers to explore the genetic diversity and evaluate the genetic structure in Tunisian germplasm. The results demonstrated significant diversity in the Tunisian durum wheat germplasm, with clear differentiation between traditional varieties and modern cultivars. The population structure analysis allowed the identification of five subpopulations, two of which appear to be more strongly represented in germplasm collected in central and southern Tunisia, where environmental conditions at critical development phases of the plant are harsher. Moreover these subpopulations are underrepresented in modern varieties, suggesting that traits of adaptation useful for breeding more resilient varieties might be present in central and southern germplasm. Moreover, our results will support, the activity of in situ on farm conservation of Tunisian durum wheat germplasm started by the National Gene Bank of Tunisia along with the ex situ approach.

Broadening the bread wheat D genome

Although Ae. tauschii has been extensively utilised for wheat breeding, the D-genome-containing allopolyploids have largely remained unexploited. In this review, we discuss approaches that can be used to exploit the D genomes of the different Aegilops species for the improvement of bread wheat. The D genome of allohexaploid bread wheat (Triticum aestivum, 2n = AABBDD) is the least diverse of the three wheat genomes and is unarguably less diverse than that of diploid progenitor Aegilops tauschii (2n = DD). Useful genetic variation and phenotypic traits also exist within each of the wheat group species containing a copy of the D genome: allopolyploid Aegilops species Ae. cylindrica (2n = D^cD^cC^cC^c), Ae. crassa 4x (2n = D¹D¹X^crX^cr), Ae. crassa 6x (2n = D¹D¹X^crX^crD^crD^cr), Ae. ventricosa (2n = D^vD^vN^vN^v), Ae. vavilovii (2n = D¹D¹X^crX^crS^vS^v) and Ae. juvenalis (2n = D¹D¹X^crX^crU^jU^j). Although Ae. tauschii has been extensively utilised for wheat breeding, the D-genome-containing allopolyploids have largely remained unexploited. Some of these D genomes appear to be modified relative to the bread wheat and Ae. tauschii D genomes, and others present in the allopolyploids may also contain useful variation as a result of adaptation to an allopolyploid, multi-genome environment. We summarise the genetic relationships, karyotypic variation and phenotypic traits known to be present in each of the D genome species that could be of relevance for bread wheat improvement and discuss approaches that can be used to exploit the D genomes of the different Aegilops species for the improvement of bread wheat. Better understanding of factors controlling chromosome inheritance and recombination in wheat group interspecific hybrids, as well as effective utilisation of new and developing genetics and genomics technologies, have great potential to improve the agronomic potential of the bread wheat D genome.

Imputation accuracy of wheat genotyping-by-sequencing (GBS) data using barley and wheat genome references

Genotyping-by-sequencing (GBS) provides high SNP coverage and has recently emerged as a popular technology for genetic and breeding applications in bread wheat (Triticum aestivum L.) and many other plant species. Although GBS can discover millions of SNPs, a high rate of missing data is a major concern for many applications. Accurate imputation of those missing data can significantly improve the utility of GBS data. This study compared imputation accuracies among four genome references including three wheat references (Chinese Spring survey sequence, W7984, and IWGSC RefSeq v1.0) and one barley reference genome by comparing imputed data derived from low-depth sequencing to actual data from high-depth sequencing. After imputation, the average number of imputed data points was the highest in the B genome (~48.99%). The D genome had the lowest imputed data points (~15.02%) but the highest imputation accuracy. Among the four reference genomes, IWGSC RefSeq v1.0 reference provided the most imputed data points, but the lowest imputation accuracy for the SNPs with < 10% minor allele frequency (MAF). The W7984 reference, however, provided the highest imputation accuracy for the SNPs with < 10% MAF.

Resequencing 200 Flax Cultivated Accessions Identifies Candidate Genes Related to Seed Size and Weight and Reveals Signatures of Artificial Selection

Seed size and weight are key traits determining crop yield, which often undergo strongly artificial selection during crop domestication. Although seed sizes differ significantly between oil flax and fiber flax, the genetic basis of morphological differences and artificial selection characteristics in seed size remains largely unclear. Here we re-sequenced 200 flax cultivated accessions to generate a genome variation map based on chromosome assembly reference genomes. We provide evidence that oil flax group is the ancestor of cultivated flax, and the oil-fiber dual purpose group (OF) is the evolutionary intermediate transition state between oil and fiber flax. Genome-wide association studies (GWAS) were combined with LD Heatmap to identify candidate regions related to seed size and weight, then candidate genes were screened based on detailed functional annotations and estimation of nucleotide polymorphism effects. Using this strategy, we obtained 13 candidate genes related to seed size and weight. Selective sweeps analysis indicates human-involved selection of small seeds during the oil to fiber flax transition. Our study shows the existence of elite alleles for seed size and weight in flax germplasm and provides molecular insights into approaches for further improvement.

Identification of Novel Quantitative Trait Loci Linked to Crown Rot Resistance in Spring Wheat

Crown rot (CR), caused by various Fusarium species, is a major disease in many cereal-growing regions worldwide. Fusarium culmorum is one of the most important species, which can cause significant yield losses in wheat. A set of 126 advanced International Maize and Wheat Improvement Center (CIMMYT) spring bread wheat lines were phenotyped against CR for field crown, greenhouse crown and stem, and growth room crown resistance scores. Of these, 107 lines were genotyped using Diversity Array Technology (DArT) markers to identify quantitative trait loci linked to CR resistance by genome-wide association study. Results of the population structure analysis grouped the accessions into three sub-groups. Genome wide linkage disequilibrium was large and declined on average within 20 cM (centi-Morgan) in the panel. General linear model (GLM), mixed linear model (MLM), and naïve models were tested for each CR score and the best model was selected based on quarantine-quarantine plots. Three marker-trait associations (MTAs) were identified linked to CR resistance; two of these on chromosome 3B were associated with field crown scores, each explaining 11.4% of the phenotypic variation and the third MTA on chromosome 2D was associated with greenhouse stem score and explained 11.6% of the phenotypic variation. Together, these newly identified loci provide opportunity for wheat breeders to exploit in enhancing CR resistance via marker-assisted selection or deployment in genomic selection in wheat breeding programs.

Genetic Diversity and Molecular Characterization of Iranian Durum Wheat Landraces (Triticum turgidum durum (Desf.) Husn.) Using DArT Markers

Durum wheat is grown mainly in rain-fed regions of Iran and the Mediterranean district under stressful conditions. Different environmental conditions and agricultural practices among ancient communities have led to the development of locally adapted genotypes known as landraces. Landraces are a valued source of genetic variety and show definite adaptation to local environmental conditions according to their home of origin. This study aimed to explore linkage disequilibrium (LD) analysis and the population structure and genetic diversity of Iranian durum wheat landraces. In this study, population structure and genome-wide LD were investigated in 129 durum landrace accessions using 1500 DArT markers. Both structure and discriminant analysis of principal components obviously subdivided the sample collection into seven distinct groups centered on key ancestors and regions of origin of the germplasm. Genetic diversity among the populations was primarily within population (68 vs. 32%). Mean LD values across the entire population sample decayed below r² of 0.11 after 1 cM. LD decay of genomes A and B of Iranian durum wheat landrace is approximately 2-3 cM (r² = 0.11) and approximately 0.5 cM (r² = 0.12), respectively. Altogether, low LD decay, a high number of subpopulations, and the high existence of genetic diversity among and within populations were characteristics of the Iranian durum landrace collection. Hence, the existing genetic diversity within the population can be associated with the very long evolutionary history of plants in Iran. The populations we studied are hence presented as a valuable resource that can be used in basic and applied research in durum wheat breeding.

Molecular analysis of buckwheat using gene specific markers

Buckwheat (Fagopyrium esculentum) is a pseudo-cereal which has spread troughout the world and nowadays it represents cultural, economic and nutritionally important pseudocereal. It´s enviromentally friendly, characterized by high fiber, routine, protein and B vitamins, and is general-purpose. The goal of the present study was to analyze 17 genotypes of buckwheat by using 7 SCoT markers. In total, 52 fragments were detected, of which 38 were polymorphic. The average number of polymorphic fragments was 5.43. The most polymorphic fragments were detected in SCoT 26 and SCoT 29 markers, and the average percentage of polymorphism was 73.36 %. SCoT 29 reached the highest percentage of polymorphism (87.5 %) and SCoT 36 was lowest (60 %). The DI values ”‹”‹ranged from 0.625 (SCoT 36) to 0.887 (SCoT 26) and the average DI value was 0.749. The average PIC value was 0.729 with PIC values ranging from 0.386 (SCoT 36) to 0.831 (SCoT 26). To determine the genetic diversity of 17 genotypes of the buckwheat, a dendrogram was created using the hierarchical cluster analysis. The genotypes were divided into two major clusters (I and II). Cluster I was divided into three other subgroups. Sixteen genotypes were included in cluster I and the genotype of Madawaska (USA) was genetically the farthest in cluster II. Genetically the closest were the varieties of Ballada (Russia) and Bamby (Austria). Used SCoT markers were sufficiently polymorphic, were able identify and differentiate chosen set of buckwheat genotypes.

Genome-Wide Association Study for Spot Blotch Resistance in Hard Winter Wheat

Spot blotch (SB) caused by Cochliobolus sativus (anamorph: Bipolaris sorokiniana) is an economically important disease of wheat worldwide. Under a severe epidemic condition, the disease can cause yield losses up to 70%. Previous approaches like bi-parental mapping for identifying SB resistant genes/QTLs exploited only a limited portion of the available genetic diversity with a lower capacity to detect polygenic traits, and had a lower marker density. In this study, we performed genome-wide association study (GWAS) for SB resistance in hard winter wheat association mapping panel (HWWAMP) of 294 genotypes. The HWWAMP was evaluated for response to B. sorokiniana (isolate SD40), and a range of reactions was observed with 10 resistant, 38 moderately resistant, 120 moderately resistant- moderately susceptible, 111 moderately susceptible, and 15 susceptible genotypes. GWAS using 15,590 high-quality SNPs and 294 genotypes we identified six QTLs (p = <0.001) on chromosomes 2D, 3A, 4A, 4B, 5A, and 7B that collectively explained 30% of the total variation for SB resistance. Highly associated SNPs were identified for all six QTLs, QSb.sdsu-2D.1 (SNP: Kukri_c31121_1460, R2 = 4%), QSb.sdsu-3A.1 (SNP: Excalibur_c46082_440, R2 = 4%), QSb.sdsu-4A.1 (SNP: IWA8475, R2 = 5.5%), QSb.sdsu-4B.1 (SNP: Excalibur_rep_c79414_306, R2 = 4%), QSb.sdsu-5A.1 (SNP: Kukri_rep_c104877_2166, R2 = 6%), and QSb.sdsu-7B.1 (SNP: TA005844-0160, R2 = 6%). Our study not only validates three (2D, 5A, and 7B) genomic regions identified in previous studies but also provides highly associated SNP markers for marker assisted selection. In addition, we identified three novel QTLs (QSb.sdsu-3A.1, QSb.sdsu-4A.1, and QSb.sdsu-4B.1) for SB resistance in wheat. Gene annotation analysis of the candidate regions identified nine NBS-LRR and 38 other plant defense-related protein families across multiple QTLs, and these could be used for fine mapping and further characterization of SB resistance in wheat. Comparative analysis with barley indicated the SB resistance locus on wheat chromosomes 2D, 3A, 5A, and 7B identified in our study are syntenic to the previously identified SB resistance locus on chromosomes 2H, 3H, 5H, and 7H in barley. The 10 highly resistant genotypes and SNP markers identified in our study could be very useful resources for breeding of SB resistance in wheat.

Variability Assessment of Aromatic Rice Germplasm by Pheno-Genomic traits and Population Structure Analysis

While the pleasant scent of aromatic rice is making it more popular, with demand for aromatic rice expected to rise in future, varieties of this have low yield potential. Genetic diversity and population structure of aromatic germplasm provide valuable information for yield improvement which has potential market value and farm profit. Here, we show diversity and population structure of 113 rice germplasm based on phenotypic and genotypic traits. Phenotypic traits showed that considerable variation existed across the germplasm. Based on Shannon-Weaver index, the most variable phenotypic trait was lemma-palea color. Detecting 140 alleles, 11 were unique and suitable as a germplasm diagnostic tool. Phylogenetic cluster analysis using genotypic traits classified germplasm into three major groups. Moreover, model-based population structure analysis divided all germplasm into three groups, confirmed by principal component and neighbors joining tree analyses. An analysis of molecular variance (AMOVA) and pairwise FST test showed significant differentiation among all population pairs, ranging from 0.023 to 0.068, suggesting that all three groups differed. Significant correlation coefficient was detected between phenotypic and genotypic traits which could be valuable to select further improvement of germplasm. Findings from this study have the potential for future use in aromatic rice molecular breeding programs.

Genetic Diversity and Population Structure of F3:6 Nebraska Winter Wheat Genotypes Using Genotyping-By-Sequencing

The availability of information on the genetic diversity and population structure in wheat (Triticum aestivum L.) breeding lines will help wheat breeders to better use their genetic resources and manage genetic variation in their breeding program. The recent advances in sequencing technology provide the opportunity to identify tens or hundreds of thousands of single nucleotide polymorphism (SNPs) in large genome species (e.g., wheat). These SNPs can be utilized for understanding genetic diversity and performing genome wide association studies (GWAS) for complex traits. In this study, the genetic diversity and population structure were investigated in a set of 230 genotypes (F_3:6) derived from various crosses as a prerequisite for GWAS and genomic selection. Genotyping-by-sequencing provided 25,566 high-quality SNPs. The polymorphism information content (PIC) across chromosomes ranged from 0.09 to 0.37 with an average of 0.23. The distribution of SNPs markers on the 21 chromosomes ranged from 319 on chromosome 3D to 2,370 on chromosome 3B. The analysis of population structure revealed three subpopulations (G1, G2, and G3). Analysis of molecular variance identified 8% variance among and 92% within subpopulations. Of the three subpopulations, G2 had the highest level of genetic diversity based on three genetic diversity indices: Shannon’s information index (I) = 0.494, diversity index (h) = 0.328 and unbiased diversity index (uh) = 0.331, while G3 had lowest level of genetic diversity (I = 0.348, h = 0.226 and uh = 0.236). This high genetic diversity identified among the subpopulations can be used to develop new wheat cultivars.

Three endoplasmic reticulum-associated fatty acyl-coenzyme a reductases were involved in the production of primary alcohols in hexaploid wheat (Triticum aestivum L.)

BACKGROUND

The cuticle covers the surface of the polysaccharide cell wall of leaf epidermal cells and forms an essential diffusion barrier between the plant and the environment. The cuticle is composed of cutin and wax. Cuticular wax plays an important role in the survival of plants by serving as the interface between plants and their biotic and abiotic environments, especially restricting nonstomatal water loss. Leaf cuticular waxes of hexaploid wheat at the seedling stage mainly consist of primary alcohols, aldehydes, fatty acids, alkane and esters. Primary alcohols account for more than 80% of the total wax load. Therefore, we cloned several genes encoding fatty acyl-coenzyme A reductases from wheat and analyzed their function in yeast and plants. We propose the potential use of these genes in wheat genetic breeding.

RESULTS

We reported the cloning and characterization of three TaFARs, namely TaFAR6, TaFAR7 and TaFAR8, encoding fatty acyl-coenzyme A reductases (FAR) in wheat leaf cuticle. Expression analysis revealed that TaFAR6, TaFAR7 and TaFAR8 were expressed at the higher levels in the seedling leaf blades, and were expressed moderately or weakly in stamen, glumes, peduncle, flag leaf blade, sheath, spike, and pistil. The heterologous expression of three TaFARs in yeast (Saccharomyces cerevisiae) led to the production of C24:0 and C26:0 primary alcohols. Transgenic expression of the three TaFARs in tomato (Solanum lycopersicum) and rice (Oryza sativa) led to increased accumulation of C24:0-C30:0 primary alcohols. Transient expression of GFP protein-tagged TaFARs revealed that the three TaFAR proteins were localized to the endoplasmic reticulum (ER), the site of wax biosynthesis. The three TaFAR genes were transcriptionally induced by drought, cold, heat, powdery mildew (Blumeria graminis) infection, abscisic acid (ABA) and methyl jasmonate (MeJa) treatments.

CONCLUSIONS

These results indicated that wheat TaFAR6, TaFAR7 and TaFAR8 are involved in biosynthesis of very-long-chain primary alcohols in hexaploid wheat and in response to multiple environmental stresses.

Analysis of the Genetic Diversity and Population Structure of Austrian and Belgian Wheat Germplasm within a Regional Context Based on DArT Markers

Analysis of crop genetic diversity and structure provides valuable information needed to broaden the narrow genetic base as well as to enhance the breeding and conservation strategies of crops. In this study, 95 Austrian and Belgian wheat cultivars maintained at the Centre for Genetic Resources (CGN) in the Netherlands were characterised using 1052 diversity array technology (DArT) markers to evaluate their genetic diversity, relationships and population structure. The rarefacted allelic richness recorded in the Austrian and Belgian breeding pools (A₂₅ = 1.396 and 1.341, respectively) indicated that the Austrian germplasm contained a higher genetic diversity than the Belgian pool. The expected heterozygosity (H_E) values of the Austrian and Belgian pools were 0.411 and 0.375, respectively. Moreover, the values of the polymorphic information content (PIC) of the Austrian and Belgian pools were 0.337 and 0.298, respectively. Neighbour-joining tree divided each of the Austrian and Belgian germplasm pools into two genetically distinct groups. The structure analyses of the Austrian and Belgian pools were in a complete concordance with their neighbour-joining trees. Furthermore, the 95 cultivars were compared to 618 wheat genotypes from nine European countries based on a total of 141 common DArT markers in order to place the Austrian and Belgian wheat germplasm in a wider European context. The rarefacted allelic richness (A₁₀) varied from 1.224 (Denmark) to 1.397 (Austria). Cluster and principal coordinates (PCoA) analyses divided the wheat genotypes of the nine European countries into two main clusters. The first cluster comprised the Northern and Western European wheat genotypes, whereas the second included the Central European cultivars. The structure analysis of the 618 European wheat genotypes was in a complete concordance with the results of cluster and PCoA analyses. Interestingly, a highly significant difference was recorded between regions (26.53%). In conclusion, this is the first study to reveal the high diversity levels and structure of the uncharacterised Austrian and Belgian wheat germplasm maintained at the CGN as well as place them in a wider European context. The results should help plant breeders to utilise the most promising wheat genotypes of this study in future breeding programmes for enhancing wheat cultivars.

Genome-Wide Association Study for Spot Blotch Resistance in Hard Winter Wheat

Spot blotch (SB) caused by Cochliobolus sativus (anamorph: Bipolaris sorokiniana) is an economically important disease of wheat worldwide. Under a severe epidemic condition, the disease can cause yield losses up to 70%. Previous approaches like bi-parental mapping for identifying SB resistant genes/QTLs exploited only a limited portion of the available genetic diversity with a lower capacity to detect polygenic traits, and had a lower marker density. In this study, we performed genome-wide association study (GWAS) for SB resistance in hard winter wheat association mapping panel (HWWAMP) of 294 genotypes. The HWWAMP was evaluated for response to B. sorokiniana (isolate SD40), and a range of reactions was observed with 10 resistant, 38 moderately resistant, 120 moderately resistant- moderately susceptible, 111 moderately susceptible, and 15 susceptible genotypes. GWAS using 15,590 high-quality SNPs and 294 genotypes we identified six QTLs (p = <0.001) on chromosomes 2D, 3A, 4A, 4B, 5A, and 7B that collectively explained 30% of the total variation for SB resistance. Highly associated SNPs were identified for all six QTLs, QSb.sdsu-2D.1 (SNP: Kukri_c31121_1460, R² = 4%), QSb.sdsu-3A.1 (SNP: Excalibur_c46082_440, R² = 4%), QSb.sdsu-4A.1 (SNP: IWA8475, R² = 5.5%), QSb.sdsu-4B.1 (SNP: Excalibur_rep_c79414_306, R² = 4%), QSb.sdsu-5A.1 (SNP: Kukri_rep_c104877_2166, R² = 6%), and QSb.sdsu-7B.1 (SNP: TA005844-0160, R² = 6%). Our study not only validates three (2D, 5A, and 7B) genomic regions identified in previous studies but also provides highly associated SNP markers for marker assisted selection. In addition, we identified three novel QTLs (QSb.sdsu-3A.1, QSb.sdsu-4A.1, and QSb.sdsu-4B.1) for SB resistance in wheat. Gene annotation analysis of the candidate regions identified nine NBS-LRR and 38 other plant defense-related protein families across multiple QTLs, and these could be used for fine mapping and further characterization of SB resistance in wheat. Comparative analysis with barley indicated the SB resistance locus on wheat chromosomes 2D, 3A, 5A, and 7B identified in our study are syntenic to the previously identified SB resistance locus on chromosomes 2H, 3H, 5H, and 7H in barley. The 10 highly resistant genotypes and SNP markers identified in our study could be very useful resources for breeding of SB resistance in wheat.

A novel linkage-disequilibrium corrected genomic relationship matrix for SNP-heritability estimation and genomic prediction

Single nucleotide polymorphism (SNP)-heritability estimation is an important topic in several research fields, including animal, plant and human genetics, as well as in ecology. Linear mixed model estimation of SNP-heritability uses the structures of genomic relationships between individuals, which is constructed from genome-wide sets of SNP-markers that are generally weighted equally in their contributions. Proposed methods to handle dependence between SNPs include, “thinning” the marker set by linkage disequilibrium (LD)-pruning, the use of haplotype-tagging of SNPs, and LD-weighting of the SNP-contributions. For improved estimation, we propose a new conceptual framework for genomic relationship matrix, in which Mahalanobis distance-based LD-correction is used in a linear mixed model estimation of SNP-heritability. The superiority of the presented method is illustrated and compared to mixed-model analyses using a VanRaden genomic relationship matrix, a matrix used by GCTA and a matrix employing LD-weighting (as implemented in the LDAK software) in simulated (using real human, rice and cattle genotypes) and real (maize, rice and mice) datasets. Despite of the computational difficulties, our results suggest that by using the proposed method one can improve the accuracy of SNP-heritability estimates in datasets with high LD.

Genome-wide association mapping of black point reaction in common wheat (Triticum aestivum L.)

BACKGROUND

Black point is a serious threat to wheat production and can be managed by host resistance. Marker-assisted selection (MAS) has the potential to accelerate genetic improvement of black point resistance in wheat breeding. We performed a genome-wide association study (GWAS) using the high-density wheat 90 K and 660 K single nucleotide polymorphism (SNP) assays to better understand the genetic basis of black point resistance and identify associated molecular markers.

RESULTS

Black point reactions were evaluated in 166 elite wheat cultivars in five environments. Twenty-five unique loci were identified on chromosomes 2A, 2B, 3A, 3B (2), 3D, 4B (2), 5A (3), 5B (3), 6A, 6B, 6D, 7A (5), 7B and 7D (2), respectively, explaining phenotypic variation ranging from 7.9 to 18.0%. The highest number of loci was detected in the A genome (11), followed by the B (10) and D (4) genomes. Among these, 13 were identified in two or more environments. Seven loci coincided with known genes or quantitative trait locus (QTL), whereas the other 18 were potentially novel loci. Linear regression showed a clear dependence of black point scores on the number of favorable alleles, suggesting that QTL pyramiding will be an effective approach to increase resistance. In silico analysis of sequences of resistance-associated SNPs identified 6 genes possibly involved in oxidase, signal transduction and stress resistance as candidate genes involved in black point reaction.

CONCLUSION

SNP markers significantly associated with black point resistance and accessions with a larger number of resistance alleles can be used to further enhance black point resistance in breeding. This study provides new insights into the genetic architecture of black point reaction.

Genetic diversity and structure of elite cotton germplasm (Gossypium hirsutum L.) using genome-wide SNP data

Cotton (Gossypium spp.) is the most important natural textile fiber crop, and Gossypium hirsutum L. is responsible for 90% of the annual cotton crop in the world. Information on cotton genetic diversity and population structure is essential for new breeding lines. In this study, we analyzed population structure and genetic diversity of 288 elite Gossypium hirsutum cultivar accessions collected from around the world, and especially from China, using genome-wide single nucleotide polymorphisms (SNP) markers. The average polymorphsim information content (PIC) was 0.25, indicating a relatively low degree of genetic diversity. Population structure analysis revealed extensive admixture and identified three subgroups. Phylogenetic analysis supported the subgroups identified by STRUCTURE. The results from both population structure and phylogenetic analysis were, for the most part, in agreement with pedigree information. Analysis of molecular variance revealed a larger amount of variation was due to diversity within the groups. Establishment of genetic diversity and population structure from this study could be useful for genetic and genomic analysis and systematic utilization of the standing genetic variation in upland cotton.

Genomewide association study for seeding emergence and tiller number using SNP markers in an elite winter wheat population

Seeding emergence and tiller number are the most important traits for wheat (Triticum aestivum L.) yield, but the inheritance of seeding emergence and tillering is poorly understood. We conducted a genomewide association study focussing on seeding emergence and tiller number at different growth stages with a panel of 205 elite winter wheat accessions. The population was genotyped with a high-density Illumina iSelect 90K SNPs assay. A total of 31 loci were found to be associated with seeding emergence rate (SER) and tiller number in different growth stages. Loci distributed among 12 chromosomes accounted for 5.35 to 11.33% of the observed phenotypic variation. With this information, 10 stable SNPs were identified for eventual development of cleaved amplified polymorphic sequence markers for SER and tiller number in different growth stages. Additionally, a set of elite alleles were identified, such as Ra_c14761_1348-T, which may increase SER by 13.35%, and Excalibur_c11045_236-A and BobWhite_c8436_391-T, which may increase the rate of available tillering by 14.78 and 8.47%, respectively. These results should provide valuable information for marker-assisted selection and parental selection in wheat breeding programmes.

SNP-based association analysis for seedling traits in durum wheat (Triticum turgidum L. durum (Desf.))

In the present study, 150 accessions of worldwide originated durum wheat germplasm (Triticum turgidum spp. durum) were observed for major seedling traits and their growth. The accessions were evaluated for major seedling traits under controlled conditions of hydroponics at the 13th, 20th, 27th and 34th day-after germination. Biomass traits were measured at the 34th day-after germination. Correlation analysis was conducted among the seedling traits and three field traits at maturity, plant height, grain weight and 1000-grain weight observed in four consecutive years. Associations of the measured seedling traits and SNP markers were analyzed based on the mixed linear model (MLM). The results indicated that highly significant genetic variation and robust heritability were found for the seedling and field mature traits. In total, 259 significant associations were detected for all the traits and four growth stages. The phenotypic variation explained (R2) by a single SNP marker is higher than 10% for most (84%) of the significant SNP markers. Forty-six SNP markers associated with multiple traits, indicating non-neglectable pleiotropy in seedling stage. The associated SNP markers could be helpful for genetic analysis of seedling traits, and marker-assisted breeding of new wheat varieties with strong seedling vigor.

Genetic Structure of Modern Durum Wheat Cultivars and Mediterranean Landraces Matches with Their Agronomic Performance

A collection of 172 durum wheat landraces from 21 Mediterranean countries and 20 modern cultivars were phenotyped in 6 environments for 14 traits including phenology, biomass, yield and yield components. The genetic structure of the collection was ascertained with 44 simple sequence repeat markers that identified 448 alleles, 226 of them with a frequency lower than 5%, and 10 alleles per locus on average. In the modern cultivars all the alleles were fixed in 59% of the markers. Total genetic diversity was H_T = 0.7080 and the genetic differentiation value was G_ST = 0.1730. STRUCTURE software allocated 90.1% of the accessions in five subpopulations, one including all modern cultivars, and the four containing landrace related to their geographic origin: eastern Mediterranean, eastern Balkans and Turkey, western Balkans and Egypt, and western Mediterranean. Mean yield of subpopulations ranged from 2.6 t ha^-1 for the western Balkan and Egyptian landraces to 4.0 t ha^-1 for modern cultivars, with the remaining three subpopulations showing similar values of 3.1 t ha^-1. Modern cultivars had the highest number of grains m^-2 and harvest index, and the shortest cycle length. The diversity was lowest in modern cultivars (H_T = 0.4835) and highest in landraces from the western Balkans and Egypt (H_T = 0.6979). Genetic diversity and AMOVA indicated that variability between subpopulations was much lower (17%) than variability within them (83%), though all subpopulations had similar biomass values in all growth stages. A dendrogram based on simple sequence repeat data matched with the clusters obtained by STRUCTURE, improving this classification for some accessions that have a large admixture. landraces included in the subpopulation from the eastern Balkans and Turkey were separated into two branches in the dendrogram drawn with phenotypic data, suggesting a different origin for the landraces collected in Serbia and Macedonia. The current study shows a reliable relationship between genetic and phenotypic population structures, and the connection of both with the geographic origin of the landraces.

Genetic Architecture of Resistance to Stripe Rust in a Global Winter Wheat Germplasm Collection

Virulence shifts in populations of Puccinia striiformis f. sp. tritici (Pst), the causal pathogen of wheat stripe rust, are a major challenge to resistance breeding. The majority of known resistance genes are already ineffective against current races of Pst, necessitating the identification and introgression of new sources of resistance. Germplasm core collections that reflect the range of genetic and phenotypic diversity of crop species are ideal platforms for examining the genetic architecture of complex traits such as resistance to stripe rust. We report the results of genetic characterization and genome-wide association analysis (GWAS) for resistance to stripe rust in a core subset of 1175 accessions in the National Small Grains Collection (NSGC) winter wheat germplasm collection, based on genotyping with the wheat 9K single nucleotide polymorphism (SNP) iSelect assay and phenotyping of seedling and adult plants under natural disease epidemics in four environments. High correlations among the field data translated into high heritability values within and across locations. Population structure was evident when accessions were grouped by stripe rust reaction. GWAS identified 127 resistance loci that were effective across at least two environments, including 20 with significant genome-wide adjusted P-values. Based on relative map positions of previously reported genes and QTL, five of the QTL with significant genome-wide adjusted P-values in this study represent potentially new loci. This study provides an overview of the diversity of Pst resistance in the NSGC winter wheat germplasm core collection, which can be exploited for diversification of stripe rust resistance in breeding programs.

Molecular and Morpho-Agronomical Characterization of Root Architecture at Seedling and Reproductive Stages for Drought Tolerance in Wheat

Water availability is a major limiting factor for wheat (Triticum aestivum L.) production in rain-fed agricultural systems worldwide. Root architecture is important for water and nutrition acquisition for all crops, including wheat. A set of 158 diverse wheat genotypes of Australian (72) and Indian (86) origin were studied for morpho-agronomical traits in field under irrigated and drought stress conditions during 2010-11 and 2011-12.Out of these 31 Indian wheat genotypes comprising 28 hexaploid (Triticum aestivum L.) and 3 tetraploid (T. durum) were characterized for root traits at reproductive stage in polyvinyl chloride (PVC) pipes. Roots of drought tolerant genotypes grew upto137cm (C306) as compared to sensitive one of 63cm with a mean value of 94.8cm. Root architecture traits of four drought tolerant (C306, HW2004, HD2888 and NI5439) and drought sensitive (HD2877, HD2012, HD2851 and MACS2496) genotypes were also observed at 6 and 9 days old seedling stage. The genotypes did not show any significant variation for root traits except for longer coleoptiles and shoot and higher absorptive surface area in drought tolerant genotypes. The visible evaluation of root images using WinRhizo Tron root scanner of drought tolerant genotype HW2004 indicated compact root system with longer depth while drought sensitive genotype HD2877 exhibited higher horizontal root spread and less depth at reproductive stage. Thirty SSR markers were used to study genetic variation which ranged from 0.12 to 0.77 with an average value of 0.57. The genotypes were categorized into three subgroups as highly tolerant, sensitive, moderately sensitive and tolerant as intermediate group based on UPGMA cluster, STRUCTURE and principal coordinate analyses. The genotypic clustering was positively correlated to grouping based on root and morpho-agronomical traits. The genetic variability identified in current study demonstrated these traits can be used to improve drought tolerance and association mapping.

A unified framework for hybrid breeding and the establishment of heterotic groups in wheat

Global wheat genetic diversity can be used in a unified framework to support and accelerate hybrid breeding and the development of heterotic groups in wheat. Hybrid wheat breeding has great potential to increase the global wheat grain yield level particularly in view of the increasing abiotic and biotic stress challenges as well as variable climatic conditions. For the long-term success of hybrid wheat breeding and the maximum exploitation of heterosis, high-yielding heterotic patterns must be established. Here, we propose a unified framework for hybrid breeding and the establishment of heterotic groups in autogamous crops and exemplify it for hybrid wheat breeding in Germany. A key component is the establishment of genetic distance between heterotic groups and in this context, we assessed genetic diversity in a global collection of 1110 winter wheat varieties released during the past decades in 35 countries but with a focus on European origin. Our analyses revealed the absence of major population structure but nevertheless suggest genetically distinct subgroups with potential for hybrid wheat breeding. Taking our molecular results and additional phenotypic data together, we propose how global genetic diversity can be used to accelerate and support reciprocal recurrent selection for the development of genetically distinct heterotic groups in hybrid wheat breeding.

Genetic Diversity, Population Structure, and Linkage Disequilibrium in Bread Wheat (Triticum aestivum L.)

Bread wheat (Triticum aestivum L.) gene pool was analyzed with 117 microsatellite markers scattered throughout A, B, and D genomes. Ninety microsatellite markers were giving 1620 polymorphic alleles in 55 different bread wheat genotypes. These genotypes were found to be divided into three subgroups based on Bayesian model and Principal component analysis. The highest polymorphism information content value for the markers resides on A genome was estimated for wmc262 marker located on 4A chromosome with the polymorphism information content value of 0.960. The highest polymorphism information content value (0.954) among the markers known to be located on B genome was realized for wmc44 marker located on 1B chromosome. The highest polymorphism information content value for the markers specific to D genome was found in gwm174 marker located on 5D chromosome with the polymorphism information content value of 0.948. The presence of linkage disequilibrium between 81 pairwise SSR markers reside on the same chromosome was tested and very limited linkage disequilibrium was observed. The results confirmed that the most distant genotype pairs were as follows Ceyhan-99-Behoth 6, Gerek 79-Douma 40989, and Karahan-99-Douma 48114.

Development and Identification of SSR Markers Associated with Starch Properties and β-Carotene Content in the Storage Root of Sweet Potato (Ipomoea batatas L.)

Sweet potato (Ipomoea batatas L.) is a nutritious food crop and, based on the high starch content of its storage root, a potential bioethanol feedstock. Enhancing the nutritional value and starch quantity of storage roots are important goals of sweet potato breeding programs aimed at developing improved varieties for direct consumption, processing, and industrial uses. However, developing improved lines of sweet potato is challenging due to the genetic complexity of this plant and the lack of genome information. Short sequence repeat (SSR) markers are powerful molecular tools for tracking important loci in crops and for molecular-based breeding strategies; however, few SSR markers and marker-trait associations have hitherto been identified in sweet potato. In this study, we identified 1824 SSRs by using a de novo assembly of publicly available ESTs and mRNAs in sweet potato, and designed 1476 primer pairs based on SSR-containing sequences. We mapped 214 pairs of primers in a natural population comprised of 239 germplasms, and identified 1278 alleles with an average of 5.972 alleles per locus and a major allele frequency of 0.7702. Population structure analysis revealed two subpopulations in this panel of germplasms, and phenotypic characterization demonstrated that this panel is suitable for association mapping of starch-related traits. We identified 32, 16, and 17 SSR markers associated with starch content, β-carotene content, and starch composition in the storage root, respectively, using association analysis and further evaluation of a subset of sweet potato genotypes with various characteristics. The SSR markers identified here can be used to select varieties with desired traits and to investigate the genetic mechanism underlying starch and carotenoid formation in the starchy roots of sweet potato.