Genetic diversity and linkage disequilibrium in Chinese bread wheat (Triticum aestivum L.) revealed by SSR markers

Transcriptome Association Identifies Regulators of Wheat Spike Architecture

The architecture of wheat (Triticum aestivum) inflorescence and its complexity is among the most important agronomic traits that influence yield. For example, wheat spikes vary considerably in the number of spikelets, which are specialized reproductive branches, and the number of florets, which are spikelet branches that produce seeds. The large and repetitive nature of the three homologous and highly similar subgenomes of wheat has impeded attempts at using genetic approaches to uncover beneficial alleles that can be utilized for yield improvement. Using a population-associative transcriptomic approach, we analyzed the transcriptomes of developing spikes in 90 wheat lines comprising 74 landrace and 16 elite varieties and correlated expression with variations in spike complexity traits. In combination with coexpression network analysis, we inferred the identities of genes related to spike complexity. Importantly, further experimental studies identified regulatory genes whose expression is associated with and influences spike complexity. The associative transcriptomic approach utilized in this study allows rapid identification of the genetic basis of important agronomic traits in crops with complex genomes.

A genome-wide association study of 23 agronomic traits in Chinese wheat landraces

Uncovering the genetic basis of agronomic traits in wheat landraces is important for ensuring global food security via the development of improved varieties. Here, 723 wheat landraces from 10 Chinese agro-ecological zones were evaluated for 23 agronomic traits in six environments. All accessions could be clustered into five subgroups based on phenotypic data via discriminant function analysis, which was highly consistent with genotypic classification. A genome-wide association study was conducted for these traits using 52 303 DArT-seq markers to identify marker-trait associations and candidate genes. Using both the general linear model and the mixed linear model, 149 significant markers were identified for 21 agronomic traits based on best linear unbiased prediction values. Considering the linkage disequilibrium decay distance in this study, significant markers within 10 cM were combined as a quantitative trait locus (QTL), with a total of 29 QTL identified for 15 traits. Of these, five QTL for heading date, flag leaf width, peduncle length, and thousand kernel weight had been reported previously. Twenty-five candidate genes associated with significant markers were identified. These included the known vernalization genes VRN-B1 and vrn-B3 and the photoperiod response genes Ppd and PRR. Overall, this study should be helpful in elucidating the underlying genetic mechanisms of complex agronomic traits and performing marker-assisted selection in wheat.

Functional Conservation and Divergence among Homoeologs of TaSPL20 and TaSPL21, Two SBP-Box Genes Governing Yield-Related Traits in Hexaploid Wheat

Maintaining high and stable yields has become an increasing challenge in wheat breeding due to climate change. Although Squamosa-promoter binding protein (SBP)-box genes have important roles in plant development, very little is known about the actual biological functions of wheat SBP-box family members. Here, we dissect the functional conservation, divergence, and exploitation of homoeologs of two paralogous TaSPL wheat loci during domestication and breeding. TaSPL20 and TaSPL21 were highly expressed in the lemma and palea. Ectopic expressions of TaSPL20/21 in rice exhibited similar functions in terms of promoting panicle branching but had different functions during seed development. We characterized all six TaSPL20/21 genes located across the three homoeologous (A, B, and D) genomes. According to the functional analysis of naturally occurring variants in 20 environments, four favorable haplotypes were identified. Together, they reduced plant height by up to 27.5%, and TaSPL21-6D-HapII increased 1000-grain weight by 9.73%. Our study suggests that TaSPL20 and TaSPL21 homoeologs underwent diversification in function with each evolving its own distinctive characteristics. During domestication and breeding of wheat in China, favorable haplotypes of each set were selected and exploited to varying degrees due to their large effects on plant height and 1000-grain weight.

Elite Haplotypes of a Protein Kinase Gene TaSnRK2.3 Associated with Important Agronomic Traits in Common Wheat

Plant-specific protein kinase SnRK2s play crucial roles in response to various environmental stimuli. TaSnRK2.3, a SnRK2 member, was involved in the response to multiple abiotic stresses in wheat. To facilitate the use of TaSnRK2.3 in wheat breeding, the three genomic sequences of TaSnRK2.3, originating from the A, B, and D genomes of hexaploid wheat, were obtained. Sequence polymorphism assays showing 4 and 10 variations were detected at TaSnRK2.3-1A and at TaSnRK2.3-1B, respectively, yet no variation was identified at TaSnRK2.3-1D. Three haplotypes for A genome, and two main haplotypes for B genome of TaSnRK2.3 were identified in 32 genotypes. Functional markers (2.3AM1, 2.3AM2, 2.3BM1, 2.3BM2) were successfully developed to distinguish different haplotypes. Association analysis was performed with the general linear model in TASSEL 2.1. The results showed that both TaSnRK2.3-1A and TaSnRK2.3-1B were significantly associated with plant height (PH), length of peduncle and penultimate node, as well as 1,000-grain weight (TGW) under different environments. Additionally, TaSnRK2.3-1B was significantly associated with stem water-soluble carbohydrates at flowering and mid-grain filling stages. Hap-1A-1 had higher TGW and lower PH; Hap-1B-1 had higher TGW and stem water-soluble carbohydrates, as well as lower PH, thus the two haplotypes were considered as elite haplotypes. Geographic distribution and allelic frequencies indicated that the two preferred haplotypes Hap-1A-1 and Hap-1B-1 were positively selected in the process of Chinese wheat breeding. These results could be valuable for genetic improvement and germplasm enhancement using molecular marker assisted selection in wheat breeding.

The iSelect 9 K SNP analysis revealed polyploidization induced revolutionary changes and intense human selection causing strong haplotype blocks in wheat

A Chinese wheat mini core collection was genotyped using the wheat 9 K iSelect SNP array. Total 2420 and 2396 polymorphic SNPs were detected on the A and the B genome chromosomes, which formed 878 haplotype blocks. There were more blocks in the B genome, but the average block size was significantly (P < 0.05) smaller than those in the A genome. Intense selection (domestication and breeding) had a stronger effect on the A than on the B genome chromosomes. Based on the genetic pedigrees, many blocks can be traced back to a well-known Strampelli cross, which was made one century ago. Furthermore, polyploidization of wheat (both tetraploidization and hexaploidization) induced revolutionary changes in both the A and the B genomes, with a greater increase of gene diversity compared to their diploid ancestors. Modern breeding has dramatically increased diversity in the gene coding regions, though obvious blocks were formed on most of the chromosomes in both tetraploid and hexaploid wheats. Tag-SNP markers identified in this study can be used for marker assisted selection using haplotype blocks as a wheat breeding strategy. This strategy can also be employed to facilitate genome selection in other self-pollinating crop species.

Association analysis of resistance to cereal cyst nematodes (Heterodera avenae) and root lesion nematodes (Pratylenchus neglectus and P. thornei) in CIMMYT advanced spring wheat lines for semi-arid conditions

To identify loci linked to nematode resistance genes, a total of 126 of CIMMYT advanced spring wheat lines adapted to semi-arid conditions were screened for resistance to Heterodera avenae, Pratylenchus neglectus, and P. thornei, of which 107 lines were genotyped with 1,310 DArT. Association of DArT markers with nematode response was analyzed using the general linear model. Results showed that 11 markers were associated with resistance to H. avenae (pathotype Ha21), 25 markers with resistance to P. neglectus, and 9 significant markers were identified to be linked with resistance to P. thornei. In this work we confirmed that chromosome 4A (~90-105 cM) can be a source of resistance to P. thornei as has been recently reported. Other significant markers were also identified on chromosomal regions where no resistant genes have been reported for both nematodes species. These novel QTL were mapped to chromosomes 5A, 6A, and 7A for H. avenae; on chromosomes 1A, 1B, 3A, 3B, 6B, 7AS, and 7D for P. neglectus; and on chromosomes 1D, 2A, and 5B for P. thornei and represent potentially new loci linked to resistance that may be useful for selecting parents and deploying resistance into elite germplasm adapted to regions where nematodes are causing problem.

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.

Identification of SNP and SSR Markers in Finger Millet Using Next Generation Sequencing Technologies

Finger millet is an important cereal crop in eastern Africa and southern India with excellent grain storage quality and unique ability to thrive in extreme environmental conditions. Since negligible attention has been paid to improving this crop to date, the current study used Next Generation Sequencing (NGS) technologies to develop both Simple Sequence Repeat (SSR) and Single Nucleotide Polymorphism (SNP) markers. Genomic DNA from cultivated finger millet genotypes KNE755 and KNE796 was sequenced using both Roche 454 and Illumina technologies. Non-organelle sequencing reads were assembled into 207 Mbp representing approximately 13% of the finger millet genome. We identified 10,327 SSRs and 23,285 non-homeologous SNPs and tested 101 of each for polymorphism across a diverse set of wild and cultivated finger millet germplasm. For the 49 polymorphic SSRs, the mean polymorphism information content (PIC) was 0.42, ranging from 0.16 to 0.77. We also validated 92 SNP markers, 80 of which were polymorphic with a mean PIC of 0.29 across 30 wild and 59 cultivated accessions. Seventy-six of the 80 SNPs were polymorphic across 30 wild germplasm with a mean PIC of 0.30 while only 22 of the SNP markers showed polymorphism among the 59 cultivated accessions with an average PIC value of 0.15. Genetic diversity analysis using the polymorphic SNP markers revealed two major clusters; one of wild and another of cultivated accessions. Detailed STRUCTURE analysis confirmed this grouping pattern and further revealed 2 sub-populations within wild E. coracana subsp. africana. Both STRUCTURE and genetic diversity analysis assisted with the correct identification of the new germplasm collections. These polymorphic SSR and SNP markers are a significant addition to the existing 82 published SSRs, especially with regard to the previously reported low polymorphism levels in finger millet. Our results also reveal an unexploited finger millet genetic resource that can be included in the regional breeding programs in order to efficiently optimize productivity.

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.

The genetic diversity and structure of indica rice in China as detected by single nucleotide polymorphism analysis

Background

Rice (Oryza sativa L.) is the staple food of more than half of the world’s population. The identification of genetic diversity in local varieties of rice compared with that of improved or introduced varieties is important in breeding elite varieties for sustainable agriculture. Array-based single nucleotide polymorphism (SNP) detection is a useful technique for such studies and breeding applications.

Results

We developed a 5291-SNP genome-wide array and used it to genotype 471 indica rice accessions in China using Illumina’s Infinium technology. Local, introduced, and improved rice varieties were clustered into three sub-groups, with some overlapping shown in principal component analysis and neighbor-joining tree, also confirmed by model-based structure. A minor allele frequency ≥0.2 was observed in 72 % of polymorphic SNPs in local rice varieties, which was higher than that in other sub-groups. Local rice varieties also had the highest mean polymorphism information content (PIC) and genetic diversity. Analysis of molecular variance showed that 90.61 % of genetic variation was a result of differences within sub-groups.

Conclusions

Our results revealed that SNP analysis clustered local varieties, introduced varieties, and improved varieties into three clear sub-groups. The distribution of parameter PIC values on sub-group genomes revealed that genetic differentiation among them might not be on a genome-wide scale, but rather on selected loci or chromosomal intervals. The result of Gene Ontology enrichment analysis showed that genes nearby those selected SNPs associated different molecular functions or various traits among sub-groups.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-016-0361-x) contains supplementary material, which is available to authorized users.

Genome-Wide Association Studies of Anthracnose and Angular Leaf Spot Resistance in Common Bean (Phaseolus vulgaris L.)

The common bean (Phaseolus vulgaris L.) is the world’s most important legume for human consumption. Anthracnose (ANT; Colletotrichum lindemuthianum) and angular leaf spot (ALS; Pseudocercospora griseola) are complex diseases that cause major yield losses in common bean. Depending on the cultivar and environmental conditions, anthracnose and angular leaf spot infections can reduce crop yield drastically. This study aimed to estimate linkage disequilibrium levels and identify quantitative resistance loci (QRL) controlling resistance to both ANT and ALS diseases of 180 accessions of common bean using genome-wide association analysis. A randomized complete block design with four replicates was performed for the ANT and ALS experiments, with four plants per genotype in each replicate. Association mapping analyses were performed for ANT and ALS using a mixed linear model approach implemented in TASSEL. A total of 17 and 11 significant statistically associations involving SSRs were detected for ANT and ALS resistance loci, respectively. Using SNPs, 21 and 17 significant statistically associations were obtained for ANT and angular ALS, respectively, providing more associations with this marker. The SSR-IAC167 and PvM95 markers, both located on chromosome Pv03, and the SNP scaffold00021_89379, were associated with both diseases. The other markers were distributed across the entire common bean genome, with chromosomes Pv03 and Pv08 showing the greatest number of loci associated with ANT resistance. The chromosome Pv04 was the most saturated one, with six markers associated with ALS resistance. The telomeric region of this chromosome showed four markers located between approximately 2.5 Mb and 4.4 Mb. Our results demonstrate the great potential of genome-wide association studies to identify QRLs related to ANT and ALS in common bean. The results indicate a quantitative and complex inheritance pattern for both diseases in common bean. Our findings will contribute to more effective screening of elite germplasm to find resistance alleles for marker-assisted selection in breeding programs.

Two Novel AP2/EREBP Transcription Factor Genes TaPARG Have Pleiotropic Functions on Plant Architecture and Yield-Related Traits in Common Wheat

AP2/EREBPs play significant roles in plant growth and development. A novel, pleiotropic TaPARG (PLANT ARCHITECTURE-RELATED GENE), a member of the AP2/EREBP transcription factor gene family, and its flanking sequences were isolated in wheat (Triticum aestivum L.). Two TaPARG genes were identified and named as TaPARG-2A and TaPARG-2D. Their amino acid sequences were highly similar especially in the functional domains. TaPARG-2A on chromosome 2A was flanked by markers Xwmc63 and Xgwm372. TaPARG-2D was mapped to chromosome 2D. Subcellular localization revealed that TaPARG-2D was localized in the nucleus. The results of tissue expression pattern, overexpression in rice, association analysis and distinct population verification jointly revealed that TaPARG functions during the entire growth cycle of wheat. Its functions include regulation of plant architecture-related and yield-related traits. Association analysis, geographic distribution and allelic frequencies suggested that favored haplotypes Hap-2A-2 and Hap-2A-3 were selected in Chinese wheat breeding programs. Both favored haplotypes might be caused by a single amino acid substitution (His/Tyr). These results suggest that TaPARG is a regulatory factor in plant growth and development, and that the favored alleles might be useful for improving plant architecture and grain yield of wheat.

Molecular tagging of a new broad-spectrum powdery mildew resistance allele Pm2c in Chinese wheat landrace Niaomai

A new broad-spectrum powdery mildew resistance allele Pm2c was identified and mapped in Chinese wheat landrace Niaomai. Chinese wheat landrace Niaomai showed resistance to 27 of 28 Chinese Blumeria graminis f. sp tritici (Bgt) races. Genetic analysis of an F2 population and its derived F2:3 families from the cross Niaomai × Mingxian 169 and backcross population, Niaomai/2*Mingxian 169, indicated that the resistance of Niaomai to Bgt races was conferred by a single dominant resistance gene, temporarily designated PmNM. Molecular tagging showed that PmNM was located on chromosome 5DS and flanked by SSR markers Xcfd81 and Xcfd78 with the genetic distances of 0.1/0.4 cM and 4.9/7.5 cM, respectively. Niaomai showed a different array of responses compared to lines with Pm2a, Pm2b, PmD57-5D, PmLX66, PmX3986-2 and Pm48 genes, sharing the same Xcfd81 allele but differing from Xcfd78 allele for Pm2a and Pm2b lines. Allelism tests based on crosses of Niaomai with Ulka/8*Cc and KM2939 showed that PmNM is allelic to Pm2a and Pm2b. We concluded that PmNM is a new allele of Pm2, re-designated Pm2c. Pm2c could be transferred into wheat cultivars by marker-assisted selection to improve the powdery mildew resistance of breeding cultivars/lines.

A haplotype block associated with thousand-kernel weight on chromosome 5DS in common wheat (Triticum aestivum L.)

Spike number per unit area, number of grains per spike, and thousand-kernel weight (TKW) are important yield components for wheat (Triticum aestivum L.). TKW has the highest heritability among the three components. We validated 27 simple sequence repeat (SSR) loci associated with TKW in an F(2:5) breeding population grown in four environments. A cfd7(8265bp) marker on chromosome 5DS showed the strongest association with TKW and had a significantly positive effect on TKW compared to allele cfd7(8259bp), with mean increases of 5.17, 3.63, 4.11, and 5.16 g in the four environments. Markers cfd67 and cfd40 flanking cfd78 also showed significantly positive associations with TKW with increases of 5.11, 3.29, 4.31, and 4.50 g for cfd67(205), and 4.98, 3.49, 4.06, and 4.84 g for cfd40(187) compared with cfd67(203) and cfd40(190) in the four environments, respectively. A major quantitative trait locus for TKW spanning 2.94 cM on chromosome 5DS was detected by association mapping. Strong linkage disequilibrium (LD) (r(2)  > 0.2) was detected among the three linked markers, which formed three haplotype blocks in the F(2:5) breeding population. Mean TKW differences between HapB-I and HapB-II were 5.80, 4.41, 4.02, and 5.06 g in the four environments, respectively. Moreover, significant LD was detected only between cfd78 and cfd67 and between cfd67 and cfd40 in a germplasm collection. This study provides a base for cloning genes related to TKW on chromosome 5DS.

Development, identification and utilization of introgression lines using Chinese endemic and synthetic wheat as donors

Chromosome segmental introgression lines (ILs) are an effective way to utilize germplasm resources in crops. To improve agronomic traits of wheat cultivar (Triticum aestivum) Shi 4185, four sets of ILs were developed. The donors were Chinese endemic subspecies accessions Yunnan wheat (T. aestivum ssp. yunnanense) YN3, Tibetan semi-wild wheat (T. aestivum ssp. tibetanum) XZ-ZM19450, and Xinjiang wheat (T. aestivum ssp. petropavlovskyi) XJ5, and synthetic wheat HC-XM1620 derived from a cross between T. durum acc. D67.2/P66.270 with Aegilops tauschii acc. 218. Totals of 356, 366, 445 and 457 simple sequence repeat (SSR) markers were polymorphic between Shi 4185 and YN3, XZ-ZM19450, XJ5 and HC-XM1620, respectively. In total, 991 ILs were identified, including 300 derived from YN3, covering 95% of the genome of Shi 4185, 218 from XZ-ZM19450 (79%), 279 from XJ5 (97%), and 194 from HC-ZX1620 (84%). The sizes and locations of each introgression were determined from a consensus SSR linkage map. Using the ILs, 11 putative quantitative trait loci (QTLs) were identified for plant height (PH), spike length (SL) and grain number per spike (GNS). Comparative analyses of 24 elite ILs with the parents revealed that the four donor parents could be important resources to improve wheat SL and GNS. Our work offers a case for utilizing endemic landraces for QTL mapping and improvement of wheat cultivars using introgression lines.

Novel function of a putative MOC1 ortholog associated with spikelet number per spike in common wheat

Axillary meristems play an important role in determining final plant architecture and floral structures. Tomato Ls, Arabidopsis LAS and rice MOC1 are orthologous genes regulating axillary meristem initiation and outgrowth. Their functions are generally conserved but the functional specificities are divergent among species. Obvious differences between rice panicles and wheat spikes suggest the divergent functions of MOC1 and its wheat ortholog. We show that TaMOC1 might be involved in wheat spikelet development. TaMOC1 is a typical nucleus localized protein with transcriptional activation abilities. The variable N-termini of TaMOC1 protein is necessary for transcriptional activation. TaMOC1 is highly expressed in ears with length of 2, 3 and 6 cm. Significant associations between the TaMOC1-7A haplotype and spikelet number per spike were observed in ten environments over 3 years and 2 sites. TaMOC1-7A HapH, a favored haplotype acquired during wheat polyploidization, may make a positive contribution to spikelet number per spike. Based on evolutionary analysis, geographic distribution and frequency changes, TaMOC1-7A HapH might be associated with wheat domestication and Chinese wheat breeding history. The pyramiding favorable alleles of TaMOC1-7A HapH and TaSnRK2.10 (C, associated with higher TGW) can improve both spikelet number per spike and TGW simultaneously.

Subgenomic Diversity Patterns Caused by Directional Selection in Bread Wheat Gene Pools

Genetic diversity represents the fundamental key to breeding success, providing the basis for breeders to select varieties with constantly improving yield performance. On the other hand, strong selection during domestication and breeding have eliminated considerable genetic diversity in the breeding pools of major crops, causing erosion of genetic potential for adaptation to emerging challenges like climate change. High-throughput genomic technologies can address this dilemma by providing detailed knowledge to characterize and replenish genetic diversity in breeding programs. In hexaploid bread wheat (Triticum aestivum L.), the staple food for 35% of the world's population, bottlenecks during allopolyploidisation followed by strong artificial selection have considerably narrowed diversity to the extent that yields in many regions appear to be unexpectedly stagnating. In this study, we used a 90,000 single nucleotide polymorphism (SNP) wheat genotyping array to assay high-frequency, polymorphic SNP markers in 460 accessions representing different phenological diversity groups from Asian, Australian, European, and North American bread wheat breeding materials. Detailed analysis of subgroup diversity at the chromosome and subgenome scale revealed highly distinct patterns of conserved linkage disequilibrium between different gene pools. The data enable identification of genome regions in most need of rejuvenation with novel diversity and provide a high-resolution molecular basis for genomic-assisted introgression of new variation into chromosome segments surrounding directionally selected metaloci conferring important adaptation and quality traits.

Prospects and limits of marker imputation in quantitative genetic studies in European elite wheat (Triticum aestivum L.)

Background

The main goal of our study was to investigate the implementation, prospects, and limits of marker imputation for quantitative genetic studies contrasting map-independent and map-dependent algorithms. We used a diversity panel consisting of 372 European elite wheat (Triticum aestivum L.) varieties, which had been genotyped with SNP arrays, and performed intensive simulation studies.

Results

Our results clearly showed that imputation accuracy was substantially higher for map-dependent compared to map-independent methods. The accuracy of marker imputation depended strongly on the linkage disequilibrium between the markers in the reference panel and the markers to be imputed. For the decay of linkage disequilibrium present in European wheat, we concluded that around 45,000 markers are needed for low cost, low-density marker profiling. This will facilitate high imputation accuracy, also for rare alleles. Genomic selection and diversity studies profited only marginally from imputing missing values. In contrast, the power of association mapping increased substantially when missing values were imputed.

Conclusions

Imputing missing values is especially of interest for an economic implementation of association mapping in breeding populations.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1366-y) contains supplementary material, which is available to authorized users.

A yield-associated gene TaCWI, in wheat: its function, selection and evolution in global breeding revealed by haplotype analysis

Wheat anther-specific invertase genes were haplotyped in wheat. Strong allelic selection occurred during wheat polyploidization, domestication and breeding because of their association with yield traits. Plant invertase hydrolyzes sucrose into glucose and fructose. Cell wall invertase (CWI), one of the three types of invertase, is essential for plant development. Based on isolated TaCWI genes from chromosomes 4A, 5B and 5D, two SNPs were detected in the promoter region of TaCWI-4A, and four SNPs and two Indels were present in the TaCWI-5D gene. No polymorphism was detected in TaCWI-5B coding or promoter regions. CAPS markers caps4A and caps5D were developed to discriminate haplotypes of TaCWI-4A and TaCWI-5D. Marker/trait association analysis indicated that Hap-5D-C at TaCWI-5D was significantly associated with higher thousand kernel weight (TKW) in 348 Chinese modern cultivars grown in multiple environments. Geographic distributions and changes over time of favored haplotypes showed that Hap-5D-C was the most frequent haplotype in modern cultivars and was strongly positively selected in six major wheat production regions worldwide. However, selection for haplotypes at TaCWI-4A was not so evident, possibly due to balancing effects of the two haplotypes on TKW and grain number per spike (GN). In rainfed production regions, Hap-4A-C was favored because it brought more seeds, but in well irrigated conditions, Hap-4A-T was favored in modern breeding because of higher TKW. Evolutionary analysis among wheat and its relatives showed that genetic diversity of TaCWI genes on chromosomes 4A and 5D declined dramatically in progression from the diploid level to modern polyploid cultivars. There was strong allelic selection during polyploidization, domestication and breeding.

Molecular evaluation of orphan Afghan common wheat (Triticum aestivum L.) landraces collected by Dr. Kihara using single nucleotide polymorphic markers

BACKGROUND

Landraces are an important source of genetic diversity in common wheat, but archival collections of Afghan wheat landraces remain poorly characterised. The recent development of array based marker systems, particularly single nucleotide polymorphism (SNP) markers, provide an excellent tool for examining the genetic diversity of local populations. Here we used SNP analysis to demonstrate the importance of Afghan wheat landraces and found tremendous genetic diversity and province-specific characteristics unique to this geographic region.

RESULTS

A total of 446 Afghan wheat landraces were analysed using genotype by sequencing (GBS) arrays containing ~10 K unique markers. Pair-wise genetic distance analyses revealed significant genetic distances between landraces, particularly among those collected from distanced provinces. From these analyses, we were able to divide the landraces into 14 major classes, with the greatest degree of diversity evident among landraces isolated from Badakhshan province. Population-based analyses revealed an additional 15 sub-populations within our germplasm, and significant correlations were evident in both the provincial and botanical varieties. Genetic distance analysis was used to identify differences among provinces, with the strongest correlations seen between landraces from Herat and Ghor province, followed closely by those between Badakhshan and Takhar provinces. This result closely resembles existing agro-climatic zones within Afghanistan, as well as the wheat varieties commonly cultivated within these regions. Molecular variance analysis showed a higher proportion of intra-province variation among landraces compared with variation among all landraces as a whole.

CONCLUSION

The SNP analyses presented here highlight the importance and genetic diversity of Afghan wheat landraces. Furthermore, these data strongly refute a previous analysis that suggested low genetic diverse within this germplasm. Ongoing analyses include phenotypic characterisation of these landraces to identify functional traits associated with individual genotypes. Taken together, these analyses can be used to help improve wheat cultivation in Afghanistan, while providing insights into the evolution and selective pressures underlying these distinct landraces.

Association mapping and validation of QTLs for flour yield in the soft winter wheat variety Kitahonami

The winter wheat variety Kitahonami shows a superior flour yield in comparison to other Japanese soft wheat varieties. To map the quantitative trait loci (QTL) associated with this trait, association mapping was performed using a panel of lines from Kitahonami's pedigree, along with leading Japanese varieties and advanced breeding lines. Using a mixed linear model corrected for kernel types and familial relatedness, 62 marker-trait associations for flour yield were identified and classified into 21 QTLs. In eighteen of these, Kitahonami alleles showed positive effects. Pedigree analysis demonstrated that a continuous pyramiding of QTLs had occurred throughout the breeding history of Kitahonami. Linkage analyses using three sets of doubled haploid populations from crosses in which Kitahonami was used as a parent were performed, leading to the validation of five of the eight QTLs tested. Among these, QTLs on chromosomes 3B and 7A showed highly significant and consistent effects across the three populations. This study shows that pedigree-based association mapping using breeding materials can be a useful method for QTL identification at the early stages of breeding programs.

TEF-7A, a transcript elongation factor gene, influences yield-related traits in bread wheat (Triticum aestivum L.)

In this study, TaTEF-7A, a member of the transcript elongation factor gene family, and its flanking sequences were isolated. TaTEF-7A was located on chromosome 7A and was flanked by markers Xwmc83 and XP3156.3. Subcellular localization revealed that TaTEF-7A protein was localized in the nucleus. This gene was expressed in all organs, but the highest expression occurred in young spikes and developing seeds. Overexpression of TaTEF-7A in Arabidopsis thaliana produced pleiotropic effects on vegetative and reproductive development that enhanced grain length, silique number, and silique length. No diversity was found in the coding region of TaTEF-7A, but 16 single nucleotide polymorphisms and Indels were detected in the promoter regions of different cultivars. Markers based on sequence variations in the promoter regions (InDel-629 and InDel-604) were developed, and three haplotypes were identified based on those markers. Haplotype-trait association analysis of the Chinese wheat mini core collection revealed that TaTEF-7A was significantly associated with grain number per spike. Phenotyping of near-isogenic lines (NILs) confirmed that TaTEF-7A increases potential grain yield and yield-related traits. Frequency changes in favoured haplotypes gradually increased in cultivars released in China from the 1940s. Geographic distributions of favoured haplotypes were characterized in six major wheat production regions worldwide. The presence of Hap-7A-3, the favoured haplotype, showed a positive correlation with yield in a global set of breeding lines. These results suggest that TaTEF-7A is a functional regulatory factor for grain number per spike and provide a basis for marker-assisted selection.

TEF-7A, a transcript elongation factor gene, influences yield-related traits in bread wheat (Triticum aestivum L.)

In this study, TaTEF-7A, a member of the transcript elongation factor gene family, and its flanking sequences were isolated. TaTEF-7A was located on chromosome 7A and was flanked by markers Xwmc83 and XP3156.3. Subcellular localization revealed that TaTEF-7A protein was localized in the nucleus. This gene was expressed in all organs, but the highest expression occurred in young spikes and developing seeds. Overexpression of TaTEF-7A in Arabidopsis thaliana produced pleiotropic effects on vegetative and reproductive development that enhanced grain length, silique number, and silique length. No diversity was found in the coding region of TaTEF-7A, but 16 single nucleotide polymorphisms and Indels were detected in the promoter regions of different cultivars. Markers based on sequence variations in the promoter regions (InDel-629 and InDel-604) were developed, and three haplotypes were identified based on those markers. Haplotype-trait association analysis of the Chinese wheat mini core collection revealed that TaTEF-7A was significantly associated with grain number per spike. Phenotyping of near-isogenic lines (NILs) confirmed that TaTEF-7A increases potential grain yield and yield-related traits. Frequency changes in favoured haplotypes gradually increased in cultivars released in China from the 1940s. Geographic distributions of favoured haplotypes were characterized in six major wheat production regions worldwide. The presence of Hap-7A-3, the favoured haplotype, showed a positive correlation with yield in a global set of breeding lines. These results suggest that TaTEF-7A is a functional regulatory factor for grain number per spike and provide a basis for marker-assisted selection.

Application of next-generation sequencing technology to study genetic diversity and identify unique SNP markers in bread wheat from Kazakhstan

BACKGROUND

New SNP marker platforms offer the opportunity to investigate the relationships between wheat cultivars from different regions and assess the mechanism and processes that have led to adaptation to particular production environments. Wheat breeding has a long history in Kazakhstan and the aim of this study was to explore the relationship between key varieties from Kazakhstan and germplasm from breeding programs for other regions.

RESULTS

The study revealed 5,898 polymorphic markers amongst ten cultivars, of which 2,730 were mapped in the consensus genetic map. Mapped SNP markers were distributed almost equally across the A and B genomes, with between 279 and 484 markers assigned to each chromosome. Marker coverage was approximately 10-fold lower in the D genome. There were 863 SNP markers identified as unique to specific cultivars, and clusters of these markers (regions containing more than three closely mapped unique SNPs) showed specific patterns on the consensus genetic map for each cultivar. Significant intra-varietal genetic polymorphism was identified in three cultivars (Tzelinnaya 3C, Kazakhstanskaya rannespelaya and Kazakhstanskaya 15). Phylogenetic analysis based on inter-varietal polymorphism showed that the very old cultivar Erythrospermum 841 was the most genetically distinct from the other nine cultivars from Kazakhstan, falling in a clade together with the American cultivar Sonora and genotypes from Central and South Asia. The modern cultivar Kazakhstanskaya 19 also fell into a separate clade, together with the American cultivar Thatcher. The remaining eight cultivars shared a single sub-clade but were categorised into four clusters.

CONCLUSION

The accumulated data for SNP marker polymorphisms amongst bread wheat genotypes from Kazakhstan may be used for studying genetic diversity in bread wheat, with potential application for marker-assisted selection and the preparation of a set of genotype-specific markers.

Genome-wide linkage disequilibrium in nine-spined stickleback populations

Variation in the extent and magnitude of genome-wide linkage disequilibrium (LD) among populations residing in different habitats has seldom been studied in wild vertebrates. We used a total of 109 microsatellite markers to quantify the level and patterns of genome-wide LD in 13 Fennoscandian nine-spined stickleback (Pungitius pungitius) populations from four (viz. marine, lake, pond, and river) different habitat types. In general, high magnitude (D' > 0.5) of LD was found both in freshwater and marine populations, and the magnitude of LD was significantly greater in inland freshwater than in marine populations. Interestingly, three coastal freshwater populations located in close geographic proximity to the marine populations exhibited similar LD patterns and genetic diversity as their marine neighbors. The greater levels of LD in inland freshwater compared with marine and costal freshwater populations can be explained in terms of their contrasting demographic histories: founder events, long-term isolation, small effective sizes, and population bottlenecks are factors likely to have contributed to the high levels of LD in the inland freshwater populations. In general, these findings shed new light on the patterns and extent of variation in genome-wide LD, as well as the ecological and evolutionary factors driving them.

Genomic prediction for rust resistance in diverse wheat landraces

We have demonstrated that genomic selection in diverse wheat landraces for resistance to leaf, stem and strip rust is possible, as genomic breeding values were moderately accurate. Markers with large effects in the Bayesian analysis confirmed many known genes, while also discovering many previously uncharacterised genome regions associated with rust scores. Genomic selection, where selection decisions are based on genomic estimated breeding values (GEBVs) derived from genome-wide DNA markers, could accelerate genetic progress in plant breeding. In this study, we assessed the accuracy of GEBVs for rust resistance in 206 hexaploid wheat (Triticum aestivum) landraces from the Watkins collection of phenotypically diverse wheat genotypes from 32 countries. The landraces were genotyped for 5,568 SNPs using an Illumina iSelect 9 K bead chip assay and phenotyped for field-based leaf rust (Lr), stem rust (Sr) and stripe rust (Yr) responses across multiple years. Genomic Best Linear Unbiased Prediction (GBLUP) and a Bayesian Regression method (BayesR) were used to predict GEBVs. Based on fivefold cross-validation, the accuracy of genomic prediction averaged across years was 0.35, 0.27 and 0.44 for Lr, Sr and Yr using GBLUP and 0.33, 0.38 and 0.30 for Lr, Sr and Yr using BayesR, respectively. Inclusion of PCR-predicted genotypes for known rust resistance genes increased accuracy more substantially when the marker was diagnostic (Lr34/Sr57/Yr18) for the presence-absence of the gene rather than just linked (Sr2). Investigation of the impact of genetic relatedness between validation and reference lines on accuracy of genomic prediction showed that accuracy will be higher when each validation line had at least one close relationship to the reference lines. Overall, the prediction accuracies achieved in this study are encouraging, and confirm the feasibility of genomic selection in wheat. In several instances, estimated marker effects were confirmed by published literature and results of mapping experiments using Watkins accessions.

Homologous haplotypes, expression, genetic effects and geographic distribution of the wheat yield gene TaGW2

BACKGROUND

TaGW2-6A, cloned in earlier research, strongly influences wheat grain width and TKW. Here, we mainly analyzed haplotypes of TaGW2-6B and their effects on TKW and interaction with haplotypes at TaGW2-6A.

RESULTS

About 2.9 kb of the promoter sequences of TaGW2-6B and TaGW2-6D were cloned in 34 bread wheat cultivars. Eleven SNPs were detected in the promoter region of TaGW2-6B, forming 4 haplotypes, but no divergence was detected in the TaGW2-6D promoter or coding region. Three molecular markers including CAPS, dCAPS and ACAS, were developed to distinguish the TaGW2-6B haplotypes. Haplotype association analysis indicated that TaGW2-6B has a stronger influence than TaGW2-6A on TKW, and Hap-6B-1 was a favored haplotype increasing grain width and weight that had undergone strong positive selection in global wheat breeding. However, clear geographic distribution differences for TaGW2-6A haplotypes were found; Hap-6A-A was favored in Chinese, Australian and Russian cultivars, whereas Hap-6A-G was preferred in European, American and CIMMYT cultivars. This difference might be caused by a flowering and maturity time difference between the two haplotypes. Hap-6A-A is the earlier type. Haplotype interaction analysis between TaGW2-6A and TaGW2-6B showed additive effects between the favored haplotypes. Hap-6A-A/Hap-6B-1 was the best combination to increase TKW. Relative expression analysis of the three TaGW2 homoeologous genes in 22 cultivars revealed that TaGW2-6A underwent the highest expression. TaGW2-6D was the least expressed during grain development and TaGW2-6B was intermediate. Diversity of the three genes was negatively correlated with their effect on TKW.

CONCLUSIONS

Genetic effects, expression patterns and historic changes of haplotypes at three homoeologous genes of TaGW2 influencing yield were dissected in wheat cultivars. Strong and constant selection to favored haplotypes has been found in global wheat breeding during the past century. This research also provides a valuable case for understanding interaction of genes that control complex traits in polyploid species.

Genetic diversity and population structure analysis of European hexaploid bread wheat (Triticum aestivum L.) varieties

Progress in plant breeding is facilitated by accurate information about genetic structure and diversity. Here, Diversity Array Technology (DArT) was used to characterize a population of 94 bread wheat (Triticum aestivum L.) varieties of mainly European origin. In total, 1,849 of 7,000 tested markers were polymorphic and could be used for population structure analysis. Two major subgroups of wheat varieties, GrI and GrII, were identified using the program STRUCTURE, and confirmed by principal component analysis (PCA). These subgroups were largely separated according to origin; GrI comprised varieties from Southern and Eastern Europe, whereas GrII contained mostly modern varieties from Western and Northern Europe. A large proportion of the markers contributing most to the genetic separation of the subgroups were located on chromosome 2D near the Reduced height 8 (Rht8) locus, and PCR-based genotyping suggested that breeding for the Rht8 allele had a major impact on subgroup separation. Consistently, analysis of linkage disequilibrium (LD) suggested that different selective pressures had acted on chromosome 2D in the two subgroups. Our data provides an overview of the allele composition of bread wheat varieties anchored to DArT markers, which will facilitate targeted combination of alleles following DArT-based QTL studies. In addition, the genetic diversity and distance data combined with specific Rht8 genotypes can now be used by breeders to guide selection of crossing parents.

Global selection on sucrose synthase haplotypes during a century of wheat breeding

Spike number per unit area, number of grains per spike, and thousand kernel weight (TKW) are important yield components. In China, increases in wheat (Triticum aestivum) yields are mainly due to increases in grain number per spike and TKW. TKW mainly depends on starch content, as starch accounts for about 70% of the grain endosperm. Sucrose synthase catalysis is the first step in the conversion of sucrose to starch, that is, the conversion of sucrose to fructose and UDP-glucose by the wheat sucrose synthase genes (TaSus1 and TaSus2) that are located on chromosomes 7A/7B/7D and 2A/2B/2D, respectively. A total of 1,520 wheat accessions were genotyped at the six loci. Two, two, five, and two haplotypes were identified at the TaSus2-2A, TaSus2-2B, TaSus1-7A, and TaSus1-7B loci, respectively. Their main variations were detected within the introns. Significant differences between the haplotypes correlated with TKW differences among 348 modern Chinese cultivars from the core collection. Frequency changes for favored haplotypes showed gradual increases in cultivars released since beginning of the last century in China, Europe, and North America. Geographic distributions and time changes of favored haplotypes were characterized in six major wheat production regions worldwide. Strong selection bottlenecks to haplotype variations occurred at polyploidization and domestication and during breeding of wheat. Genetic-effect differences between haplotypes at the same locus influence the selection time and intensity. This work shows that the endosperm starch synthesis pathway is a major target of indirect selection in global wheat breeding for higher yield.

Linkage disequilibrium and population-structure analysis among Capsicum annuum L. cultivars for use in association mapping

Knowledge of population structure and linkage disequilibrium among the worldwide collections of peppers currently classified as hot, mild, sweet and ornamental types is indispensable for applying association mapping and genomic selection to improve pepper. The current study aimed to resolve the genetic diversity and relatedness of Capsicum annuum germplasm by use of simple sequence repeat (SSR) loci across all chromosomes in samples collected in 2011 and 2012. The physical distance covered by the entire set of SSRs used was 2,265.9 Mb from the 3.48-Gb hot-pepper genome size. The model-based program STRUCTURE was used to infer five clusters, which was further confirmed by classical molecular-genetic diversity analysis. Mean heterozygosity of various loci was estimated to be 0.15. Linkage disequilibrium (LD) was used to identify 17 LD blocks across various chromosomes with sizes from 0.154 Kb to 126.28 Mb. CAMS-142 of chromosome 1 was significantly associated with both capsaicin (CA) and dihydrocapsaicin (DCA) levels. Further, CAMS-142 was located in an LD block of 98.18 Mb. CAMS-142 amplified bands of 244, 268, 283 and 326 bp. Alleles 268 and 283 bp had positive effects on both CA and DCA levels, with an average R(2) of 12.15 % (CA) and 12.3 % (DCA). Eight markers from seven different chromosomes were significantly associated with fruit weight, contributing an average effect of 15 %. CAMS-199, HpmsE082 and CAMS-190 are the three major quantitative trait loci located on chromosomes 8, 9, and 10, respectively, and were associated with fruit weight in samples from both years of the study. This research demonstrates the effectiveness of using genome-wide SSR-based markers to assess features of LD and genetic diversity within C. annuum.

Genetic structure, linkage disequilibrium and association mapping of Verticillium wilt resistance in elite cotton (Gossypium hirsutum L.) germplasm population

Understanding the population structure and linkage disequilibrium in an association panel can effectively avoid spurious associations and improve the accuracy in association mapping. In this study, one hundred and fifty eight elite cotton (Gossypium hirsutum L.) germplasm from all over the world, which were genotyped with 212 whole genome-wide marker loci and phenotyped with an disease nursery and greenhouse screening method, were assayed for population structure, linkage disequilibrium, and association mapping of Verticillium wilt resistance. A total of 480 alleles ranging from 2 to 4 per locus were identified from all collections. Model-based analysis identified two groups (G1 and G2) and seven subgroups (G1a–c, G2a–d), and differentiation analysis showed that subgroup having a single origin or pedigree was apt to differentiate with those having a mixed origin. Only 8.12% linked marker pairs showed significant LD (P<0.001) in this association panel. The LD level for linked markers is significantly higher than that for unlinked markers, suggesting that physical linkage strongly influences LD in this panel, and LD level was elevated when the panel was classified into groups and subgroups. The LD decay analysis for several chromosomes showed that different chromosomes showed a notable change in LD decay distances for the same gene pool. Based on the disease nursery and greenhouse environment, 42 marker loci associated with Verticillium wilt resistance were identified through association mapping, which widely were distributed among 15 chromosomes. Among which 10 marker loci were found to be consistent with previously identified QTLs and 32 were new unreported marker loci, and QTL clusters for Verticillium wilt resistanc on Chr.16 were also proved in our study, which was consistent with the strong linkage in this chromosome. Our results would contribute to association mapping and supply the marker candidates for marker-assisted selection of Verticillium wilt resistance in cotton.

Next generation characterisation of cereal genomes for marker discovery

Cereal crops form the bulk of the world’s food sources, and thus their importance cannot be understated. Crop breeding programs increasingly rely on high-resolution molecular genetic markers to accelerate the breeding process. The development of these markers is hampered by the complexity of some of the major cereal crop genomes, as well as the time and cost required. In this review, we address current and future methods available for the characterisation of cereal genomes, with an emphasis on faster and more cost effective approaches for genome sequencing and the development of markers for trait association and marker assisted selection (MAS) in crop breeding programs.

A high density physical map of chromosome 1BL supports evolutionary studies, map-based cloning and sequencing in wheat

BACKGROUND

As for other major crops, achieving a complete wheat genome sequence is essential for the application of genomics to breeding new and improved varieties. To overcome the complexities of the large, highly repetitive and hexaploid wheat genome, the International Wheat Genome Sequencing Consortium established a chromosome-based strategy that was validated by the construction of the physical map of chromosome 3B. Here, we present improved strategies for the construction of highly integrated and ordered wheat physical maps, using chromosome 1BL as a template, and illustrate their potential for evolutionary studies and map-based cloning.

RESULTS

Using a combination of novel high throughput marker assays and an assembly program, we developed a high quality physical map representing 93% of wheat chromosome 1BL, anchored and ordered with 5,489 markers including 1,161 genes. Analysis of the gene space organization and evolution revealed that gene distribution and conservation along the chromosome results from the superimposition of the ancestral grass and recent wheat evolutionary patterns, leading to a peak of synteny in the central part of the chromosome arm and an increased density of non-collinear genes towards the telomere. With a density of about 11 markers per Mb, the 1BL physical map provides 916 markers, including 193 genes, for fine mapping the 40 QTLs mapped on this chromosome.

CONCLUSIONS

Here, we demonstrate that high marker density physical maps can be developed in complex genomes such as wheat to accelerate map-based cloning, gain new insights into genome evolution, and provide a foundation for reference sequencing.

Population structure, genetic diversity and linkage disequilibrium in elite winter wheat assessed with SNP and SSR markers

Modern genomics approaches rely on the availability of high-throughput and high-density genotyping platforms. A major breakthrough in wheat genotyping was the development of an SNP array. In this study, we used a diverse panel of 172 elite European winter wheat lines to evaluate the utility of the SNP array for genomic analyses in wheat germplasm derived from breeding programs. We investigated population structure and genetic relatedness and found that the results obtained with SNP and SSR markers differ. This suggests that additional research is required to determine the optimum approach for the investigation of population structure and kinship. Our analysis of linkage disequilibrium (LD) showed that LD decays within approximately 5-10 cM. Moreover, we found that LD is variable along chromosomes. Our results suggest that the number of SNPs needs to be increased further to obtain a higher coverage of the chromosomes. Taken together, SNPs can be a valuable tool for genomics approaches and for a knowledge-based improvement of wheat.

Polymorphism of TaSAP1-A1 and its association with agronomic traits in wheat

TaSAP1, a member of the stress association protein (SAP) gene family from wheat, is involved in response to several abiotic stresses, including drought, salt and cold. In this study, TaSAP1-A1, a TaSAP1 member on chromosome 7A, and its flanking sequences were isolated. Polymorphism analysis indicated that the average nucleotide diversity (π) of the whole region was 0.00296. The highest nucleotide diversity occurred in the promoter region (π = 0.00631) and no polymorphism was identified in the coding region. Three markers T7AM5, T7AM2606 and T7AM39 located in the promoter region, were developed from sequence variations (InDel5-1810, SNP-2606 and InDel39-1637). Six haplotypes were identified among 300 accessions based on the three markers. TaSAP1-A1 was located on chromosome 7A using marker T7AM39 and was flanked by markers Xwmc530 and Xbarc174. QTL for yield-related traits, including 1,000-grain weight, number of grains per spike and grain yield, were located in the same region. In marker- and haplotype-trait association analyses, TaSAP1-A1 was significantly associated with 1,000-grain weight, number of grains per spike, spike length, peduncle length and total number of spikelets per spike in multiple environments. These results provide useful information for marker-assisted selection for yield-related traits under well-watered and drought-stressed conditions.

Novel insights into the composition, variation, organization, and expression of the low-molecular-weight glutenin subunit gene family in common wheat

Low-molecular-weight glutenin subunits (LMW-GS), encoded by a complex multigene family, play an important role in the processing quality of wheat flour. Although members of this gene family have been identified in several wheat varieties, the allelic variation and composition of LMW-GS genes in common wheat are not well understood. In the present study, using the LMW-GS gene molecular marker system and the full-length gene cloning method, a comprehensive molecular analysis of LMW-GS genes was conducted in a representative population, the micro-core collections (MCC) of Chinese wheat germplasm. Generally, >15 LMW-GS genes were identified from individual MCC accessions, of which 4-6 were located at the Glu-A3 locus, 3-5 at the Glu-B3 locus, and eight at the Glu-D3 locus. LMW-GS genes at the Glu-A3 locus showed the highest allelic diversity, followed by the Glu-B3 genes, while the Glu-D3 genes were extremely conserved among MCC accessions. Expression and sequence analysis showed that 9-13 active LMW-GS genes were present in each accession. Sequence identity analysis showed that all i-type genes present at the Glu-A3 locus formed a single group, the s-type genes located at Glu-B3 and Glu-D3 loci comprised a unique group, while high-diversity m-type genes were classified into four groups and detected in all Glu-3 loci. These results contribute to the functional analysis of LMW-GS genes and facilitate improvement of bread-making quality by wheat molecular breeding programmes.

Association analysis of genomic loci important for grain weight control in elite common wheat varieties cultivated with variable water and fertiliser supply

Grain weight, an essential yield component, is under strong genetic control and markedly influenced by the environment. Here, by genome-wide association analysis with a panel of 94 elite common wheat varieties, 37 loci were found significantly associated with thousand-grain weight (TGW) in one or more environments differing in water and fertiliser levels. Five loci were stably associated with TGW under all 12 environments examined. Their elite alleles had positive effects on TGW. Four, two, three, and two loci were consistently associated with TGW in the irrigated and fertilised (IF), rainfed (RF), reduced nitrogen (RN), and reduced phosphorus (RP) environments. The elite alleles of the IF-specific loci enhanced TGW under well-resourced conditions, whereas those of the RF-, RN-, or RP-specific loci conferred tolerance to the TGW decrease when irrigation, nitrogen, or phosphorus were reduced. Moreover, the elite alleles of the environment-independent and -specific loci often acted additively to enhance TGW. Four additional loci were found associated with TGW in specific locations, one of which was shown to contribute to the TGW difference between two experimental sites. Further analysis of 14 associated loci revealed that nine affected both grain length and width, whereas the remaining loci influenced either grain length or width, indicating that these loci control grain weight by regulating kernel size. Finally, the elite allele of Xpsp3152 frequently co-segregated with the larger grain haplotype of TaGW2-6A, suggesting probable genetic and functional linkages between Xpsp3152 and GW2 that are important for grain weight control in cereal plants. Our study provides new knowledge on TGW control in elite common wheat lines, which may aid the improvement of wheat grain weight trait in further research.

Identifying loci influencing grain number by microsatellite screening in bread wheat (Triticum aestivum L.)

Grain number (GN) is one of three major yield-related components in wheat. We used the Chinese wheat mini core collection to undertake a genome-wide association analysis of grain number using 531 SSR markers randomly located on all 21 chromosomes. Grain numbers of all accessions were measured in four trials, i.e. two environments in four growing seasons. Association analysis based on a mixed linear model (MLM) revealed that 27 SSR loci were significantly associated with mean GN (MGN) estimated by the best linear unbiased predictor (BLUP) method. These included numerous breeder favorable alleles with strong positive effects at 23 loci. Significant or extremely significant differences were detected on MGN between varieties conveying favored allele and varieties with other alleles. Moreover, statistical simulation showed that the favored alleles have additive genetic effects. Although modern varieties combined larger numbers of favored alleles, the numbers of favored alleles were not significantly different from those in landraces, especially those alleles contributing mostly to the phenotypic variation. These results indicate that there is still considerable genetic potential for use of markers for genome selection of GN for high yield in wheat.

Association mapping and haplotype analysis of a 3.1-Mb genomic region involved in Fusarium head blight resistance on wheat chromosome 3BS

A previous study provided an in-depth understanding of molecular population genetics of European and Asian wheat gene pools using a sequenced 3.1-Mb contig (ctg954) on chromosome 3BS. This region is believed to carry the Fhb1 gene for response to Fusarium head blight. In this study, 266 wheat accessions were evaluated in three environments for Type II FHB response based on the single floret inoculation method. Hierarchical clustering (UPGMA) based on a Manhattan dissimilarity matrix divided the accessions into eight groups according to five FHB-related traits which have a high correlation between them; Group VIII comprised six accessions with FHB response levels similar to variety Sumai 3. Based on the compressed mixed linear model (MLM), association analysis between five FHB-related traits and 42 molecular markers along the 3.1-Mb region revealed 12 significant association signals at a threshold of P<0.05. The highest proportion of phenotypic variation (6.2%) in number of diseased spikelets (NDS) occurred at locus cfb6059, and the physical distance was about 2.9 Kb between umn10 and this marker. Haplotype block (HapB) analysis using a sliding window LD of 5 markers, detected six HapBs in the 3.1-Mb region at r(2)>0.1 and P<0.001 between random closely linked markers. F-tests among Haps with frequencies >0.05 within each HapB at r(2)>0.1 and P<0.001 showed significant differences between the Hap carried by FHB resistant resources, such as Sumai 3 and Wangshuibai, and susceptible genotypes in HapB3 and HapB6. These results suggest that Fhb1 is located within HapB6, with the possibility that another gene is located at or near HapB3. SSR markers and Haps detected in this study will be helpful in further understanding the genetic basis of FHB resistance, and provide useful information for marker-assisted selection of Fhb1 in wheat breeding.

Genome-wide evaluation of genetic diversity and linkage disequilibrium in winter and spring triticale (x Triticosecale Wittmack)

Background

Recent advances in genotyping with high-density markers nowadays enable genome-wide genomic analyses in crops. A detailed characterisation of the population structure and linkage disequilibrium (LD) is essential for the application of genomic approaches and consequently for knowledge-based breeding. In this study we used the triticale-specific DArT array to analyze population structure, genetic diversity, and LD in a worldwide set of 161 winter and spring triticale lines.

Results

The principal coordinate analysis revealed that the first principal coordinate divides the triticale population into two clusters according to their growth habit. The density distributions of the first ten principal coordinates revealed that several show a distribution indicative of population structure. In addition, we observed relatedness within growth habits which was higher among the spring types than among the winter types. The genome-wide analysis of polymorphic information content (PIC) showed that the PIC is variable among and along chromosomes and that especially the R genome of spring types possesses a reduced genetic diversity. We also found that several chromosomes showed regions of high genetic distance between the two growth habits, indicative of divergent selection. Regarding linkage disequilibrium, the A and B genomes showed a similar LD of 0.24 for closely linked markers and a decay within approximately 12 cM. LD in the R genome was lower with 0.19 and decayed within a shorter map distance of approximately 5 cM. The extent of LD was generally higher for the spring types compared to the winter types. In addition, we observed strong variability of LD along the chromosomes.

Conclusions

Our results confirm winter and spring growth habit are the major contributors to population structure in triticale, and a family structure exists in both growth types. The specific patterns of genetic diversity observed within these types, such as the low diversity on some rye chromosomes of spring habits, provide a basis for targeted broadening of the available breeding germplasm. In addition, the genome-wide analysis of the extent and the pattern of LD will assist scientists and breeders alike in the implementation and the interpretation of association mapping in triticale.

Genome wide linkage disequilibrium in Chinese asparagus bean (Vigna. unguiculata ssp. sesquipedialis) germplasm: implications for domestication history and genome wide association studies

Association mapping of important traits of crop plants relies on first understanding the extent and patterns of linkage disequilibrium (LD) in the particular germplasm being investigated. We characterize here the genetic diversity, population structure and genome wide LD patterns in a set of asparagus bean (Vigna. unguiculata ssp. sesquipedialis) germplasm from China. A diverse collection of 99 asparagus bean and normal cowpea accessions were genotyped with 1127 expressed sequence tag-derived single nucleotide polymorphism markers (SNPs). The proportion of polymorphic SNPs across the collection was relatively low (39%), with an average number of SNPs per locus of 1.33. Bayesian population structure analysis indicated two subdivisions within the collection sampled that generally represented the 'standard vegetable' type (subgroup SV) and the 'non-standard vegetable' type (subgroup NSV), respectively. Level of LD (r(2)) was higher and extent of LD persisted longer in subgroup SV than in subgroup NSV, whereas LD decayed rapidly (0-2 cM) in both subgroups. LD decay distance varied among chromosomes, with the longest (≈ 5 cM) five times longer than the shortest (≈ 1 cM). Partitioning of LD variance into within- and between-subgroup components coupled with comparative LD decay analysis suggested that linkage group 5, 7 and 10 may have undergone the most intensive epistatic selection toward traits favorable for vegetable use. This work provides a first population genetic insight into domestication history of asparagus bean and demonstrates the feasibility of mapping complex traits by genome wide association study in asparagus bean using a currently available cowpea SNPs marker platform.

Identifying loci influencing 1,000-kernel weight in wheat by microsatellite screening for evidence of selection during breeding

Chinese wheat mini core collection (262 accessions) was genotyped at 531 microsatellite loci representing a mean marker density of 5.1 cM. One-thousand-kernel weights (TKW) of lines were measured in five trials (three environments in four growing seasons). Structure analysis based on 42 unlinked SSR loci indicated that the materials formed two sub-populations, viz., landraces and modern varieties. A large difference in TKW (7.08 g, P<0.001) was found between the two sub-groups. Therefore, TKW is a major yield component that was improved in the past 6 decades; it increased from a mean 31.5 g in the 1940s to 44.64 g in the 2000s, representing a 2.19 g increase in each decade. Analyses based on a mixed linear model (MLM), population structure (Q) and relative kinship (K) revealed 22 SSR loci that were significantly associated with mean TKW (MTKW) of the five trials estimated by the best linear unbiased predictor (BLUP) method. They were mainly distributed on chromosomes of homoeologous groups 1, 2, 3, 5 and 7. Six loci, cfa2234-3A, gwm156-3B, barc56-5A, gwm234-5B, wmc17-7A and cfa2257-7A individually explained more than 11.84% of the total phenotypic variation. Favored alleles for breeding at the 22 loci were inferred according to their estimated effects on MTKW based on mean difference of varieties grouped by genotypes. Statistical simulation showed that these favored alleles have additive genetic effects. Frequency changes of alleles at loci associated with TKW are much more dramatic than those at neutral loci between the sub-groups. The numbers of favored alleles in modern varieties indicate there is still considerable genetic potential for their use as markers for genome selection of TKW in wheat breeding. Alleles that can be used globally to increase TKW were inferred according to their distribution by latitude and frequency of changes between landraces and the modern varieties.

Association mapping of dynamic developmental plant height in common wheat

Drought as a major abiotic stress often occurs from stem elongation to the grain filling stage of wheat in northern China. Plant height (PH) is a suitable trait to model the dissection of drought tolerance. The purposes of the present study were to validate molecular markers for PH developmental behavior and identify elite alleles of molecular markers. After the phenotyping of 154 accessions for PH dynamic development under well-watered (WW) and drought stressed (DS) conditions, and the genotyping of 60 SSR markers from six candidate chromosome regions related to PH found in our previous linkage mapping studies, both parameters PH and drought tolerance coefficient (DTC) calculated by the conditional analysis were used for association mapping. A total of 46 significant association signals (P < 0.01) were identified in 23 markers, and phenotypic variation ranged from 7 to 50%. Among them, four markers Xgwm261-2D, Xgwm495-4B, Xbarc109-4B and Xcfd23-4D were detected under both water regimes. Furthermore, 10 markers were associated with DTC, and four with both parameters PH and DTC at the same plant development stage. The results revealed different allelic effects of associated markers; for example, the 155 bp Xgwm495-4B allele was associated with a reduced height of -11.2 cm under DS and -15.3 cm under WW, whereas the 167 bp allele exhibited increased height effects of 3.9 and 8.1 cm, respectively. This study demonstrates a strong power of joint association analysis and linkage mapping for the identification of important genes in wheat.