Genome-wide marker development for the wheat D genome based on single nucleotide polymorphisms identified from transcripts in the wild wheat progenitor Aegilops tauschii

RNA-seq analysis revealed considerable genetic diversity and enabled the development of specific KASP markers for Psathyrostachys huashanica

Psathyrostachys huashanica, which grows exclusively in Huashan, China, is an important wild relative of common wheat that has many desirable traits relevant for wheat breeding. However, the poorly characterized interspecific phylogeny and genomic variations and the relative lack of species-specific molecular markers have limited the utility of P. huashanica as a genetic resource for enhancing wheat germplasm. In this study, we sequenced the P. huashanica transcriptome, resulting in 50,337,570 clean reads that were assembled into 65,617 unigenes, of which 38,428 (58.56%) matched at least one sequence in public databases. The phylogenetic analysis of P. huashanica, Triticeae species, and Poaceae species was conducted using 68 putative orthologous gene clusters. The data revealed the distant evolutionary relationship between P. huashanica and common wheat as well as the substantial diversity between the P. huashanica genome and the wheat D genome. By comparing the transcriptomes of P. huashanica and Chinese Spring, 750,759 candidate SNPs between P. huashanica Ns genes and their common wheat orthologs were identified. Among the 90 SNPs in the exon regions with different functional annotations, 58 (64.4%) were validated as Ns genome-specific SNPs in the common wheat background by KASP genotyping assays. Marker validation analyses indicated that six specific markers can discriminate between P. huashanica and the other wheat-related species. In addition, five markers are unique to P. huashanica, P. juncea, and Leymus species, which carry the Ns genome. The Ns genome-specific markers in a wheat background were also validated regarding their specificity and stability for detecting P. huashanica chromosomes in four wheat-P. huashanica addition lines. Four and eight SNP markers were detected in wheat-P. huashanica 2Ns and 7Ns addition lines, respectively, and one marker was specific to both wheat-P. huashanica 3Ns, 4Ns, and 7Ns addition lines. These markers developed using transcriptome data may be used to elucidate the genetic relationships among Psathyrostachys, Leymus, and other closely-related species. They may also facilitate precise introgressions and the high-throughput monitoring of P. huashanica exogenous chromosomes or segments in future crop breeding programs.

Exploring the Genetic Diversity and Population Structure of Wheat Landrace Population Conserved at ICARDA Genebank

Landraces are considered a valuable source of potential genetic diversity that could be used in the selection process in any plant breeding program. Here, we assembled a population of 600 bread wheat landraces collected from eight different countries, conserved at the ICARDA's genebank, and evaluated the genetic diversity and the population structure of the landraces using single nucleotide polymorphism (SNP) markers. A total of 11,830 high-quality SNPs distributed across the genomes A (40.5%), B (45.9%), and D (13.6%) were used for the final analysis. The population structure analysis was evaluated using the model-based method (STRUCTURE) and distance-based methods [discriminant analysis of principal components (DAPC) and principal component analysis (PCA)]. The STRUCTURE method grouped the landraces into two major clusters, with the landraces from Syria and Turkey forming two clusters with high proportions of admixture, whereas the DAPC and PCA analysis grouped the population into three subpopulations mostly according to the geographical information of the landraces, i.e., Syria, Iran, and Turkey with admixture. The analysis of molecular variance revealed that the majority of the variation was due to genetic differences within the populations as compared with between subpopulations, and it was the same for both the cluster-based and distance-based methods. Genetic distance analysis was also studied to estimate the differences between the landraces from different countries, and it was observed that the maximum genetic distance (0.389) was between the landraces from Spain and Palestine, whereas the minimum genetic distance (0.013) was observed between the landraces from Syria and Turkey. It was concluded from the study that the model-based methods (DAPC and PCA) could dissect the population structure more precisely when compared with the STRUCTURE method. The population structure and genetic diversity analysis of the bread wheat landraces presented here highlight the complex genetic architecture of the landraces native to the Fertile Crescent region. The results of this study provide useful information for the genetic improvement of hexaploid wheat and facilitate the use of landraces in wheat breeding programs.

High genetic diversity in Aegilops tauschii Coss. accessions from North Iran as revealed by IRAP and REMAP markers

BACKGROUND

Aegilops tauschii Coss. as a donor of wheat D genome has an important role in wheat breeding programs. Genetic and phylogeographic diversity of 79 Ae. tauschii accessions collected from north and northwest of Iran were analyzed based on retroelement insertional polymorphisms using inter-retrotransposon amplified polymorphism (IRAP) and retrotransposon-microsatellite amplified polymorphism (REMAP) markers.

RESULTS

In total, 306 and 151 polymorphic bands were amplified in IRAP and REMAP analyses, respectively. As a result, a high level of polymorphism was observed among the studied accessions as revealed by an average of 25.5 bands per primer/primer combination and mean PIC value of 0.47 in IRAP and an average of 25.16 bands per primer combination and mean PIC value of 0.47 in REMAP. Genetic relationships of the accessions were analyzed using distance- and model-based cluster analyses.

CONCLUSION

The result showed that genetic distance did not seem to be related to geographic distribution, and the accessions could be divided into three groups, which was further supported by principal coordinate analysis. These results on genetic diversity and population structure of Ae. tauschii in Iran should provide important knowledge on genetic resources and their applications in wheat breeding programs.

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 sequencing-based coverage analyses facilitate high-resolution detection of deletions linked to phenotypes of gamma-irradiated wheat mutants

BACKGROUND

Gamma-irradiated mutants of Triticum aestivum L., hexaploid wheat, provide novel and agriculturally important traits and are used as breeding materials. However, the identification of causative genomic regions of mutant phenotypes is challenging because of the large and complicated genome of hexaploid wheat. Recently, the combined use of high-quality reference genome sequences of common wheat and cost-effective resequencing technologies has made it possible to evaluate genome-wide polymorphisms, even in complex genomes.

RESULTS

To investigate whether the genome sequencing approach can effectively detect structural variations, such as deletions, frequently caused by gamma irradiation, we selected a grain-hardness mutant from the gamma-irradiated population of Japanese elite wheat cultivar "Kitahonami." The Hardness (Ha) locus, including the puroindoline protein-encoding genes Pina-D1 and Pinb-D1 on the short arm of chromosome 5D, primarily regulates the grain hardness variation in common wheat. We performed short-read genome sequencing of wild-type and grain-hardness mutant plants, and subsequently aligned their short reads to the reference genome of the wheat cultivar "Chinese Spring." Genome-wide comparisons of depth-of-coverage between wild-type and mutant strains detected ~ 130 Mbp deletion on the short arm of chromosome 5D in the mutant genome. Molecular markers for this deletion were applied to the progeny populations generated by a cross between the wild-type and the mutant. A large deletion in the region including the Ha locus was associated with the mutant phenotype, indicating that the genome sequencing is a powerful and efficient approach for detecting a deletion marker of a gamma-irradiated mutant phenotype. In addition, we investigated a pre-harvest sprouting tolerance mutant and identified a 67.8 Mbp deletion on chromosome 3B where Viviparous-B1 and GRAS family transcription factors are located. Co-dominant markers designed to detect the deletion-polymorphism confirmed the association with low germination rate, leading to pre-harvest sprouting tolerance.

CONCLUSIONS

Short read-based genome sequencing of gamma-irradiated mutants facilitates the identification of large deletions linked to mutant phenotypes when combined with segregation analyses in progeny populations. This method allows effective application of mutants with agriculturally important traits in breeding using marker-assisted selection.

RNA-Seq-based DNA marker analysis of the genetics and molecular evolution of Triticeae species

The release of high-quality chromosome-level genome sequences of members of the Triticeae tribe has greatly facilitated genetic and genomic analyses of important crops such as wheat (Triticum aestivum) and barley (Hordeum vulgare). Due to the large diploid genome size of Triticeae plants (ca. 5 Gbp), transcript analysis is an important method for identifying genetic and genomic differences among Triticeae species. In this review, we summarize our results of RNA-Seq analyses of diploid wheat accessions belonging to the genera Aegilops and Triticum. We also describe studies of the molecular relationships among these accessions and provide insight into the evolution of common hexaploid wheat. DNA markers based on polymorphisms within species can be used to map loci of interest. Even though the genome sequence of diploid Aegilops tauschii, the D-genome donor of common wheat, has been released, the diploid barley genome continues to provide key information about the physical structures of diploid wheat genomes. We describe how a series of RNA-Seq analyses of wheat relatives has helped uncover the structural and evolutionary features of genomic and genetic systems in wild and cultivated Triticeae species.

Genetic Diversity and Population Structure Analysis of Triticum aestivum L. Landrace Panel from Afghanistan

Landraces are a potential source of genetic diversity and provide useful genetic resources to cope with the current and future challenges in crop breeding. Afghanistan is located close to the centre of origin of hexaploid wheat. Therefore, understanding the population structure and genetic diversity of Afghan wheat landraces is of enormous importance in breeding programmes for the development of high-yielding cultivars as well as broadening the genetic base of bread wheat. Here, a panel of 363 bread wheat landraces collected from seven north and north-eastern provinces of Afghanistan were evaluated for population structure and genetic diversity using single nucleotide polymorphic markers (SNPs). The genotyping-by-sequencing of studied landraces after quality control provided 4897 high-quality SNPs distributed across the genomes A (33.75%), B (38.73%), and D (27.50%). The population structure analysis was carried out by two methods using model-based STRUCTURE analysis and cluster-based discriminant analysis of principal components (DAPC). The analysis of molecular variance showed a higher proportion of variation within the sub-populations compared with the variation observed as a whole between sub-populations. STRUCTURE and DAPC analysis grouped the majority of the landraces from Badakhshan and Takhar together in one cluster and the landraces from Baghlan and Kunduz in a second cluster, which is in accordance with the micro-climatic conditions prevalent within the north-eastern agro-ecological zone. Genetic distance analysis was also studied to identify differences among the Afghan regions; the strongest correlation was observed for the Badakhshan and Takhar (0.003), whereas Samangan and Konarha (0.399) showed the highest genetic distance. The population structure and genetic diversity analysis highlighted the complex genetic variation present in the landraces which were highly correlated to the geographic origin and micro-climatic conditions within the agro-climatic zones of the landraces. The higher proportions of admixture could be attributed to historical unsupervised exchanges of seeds between the farmers of the central and north-eastern provinces of Afghanistan. The results of this study will provide useful information for genetic improvement in wheat and is essential for association mapping and genomic prediction studies to identify novel sources for resistance to abiotic and biotic stresses.

Development of RNA-seq-based molecular markers for characterizing Thinopyrum bessarabicum and Secale introgressions in wheat

Sequence-based markers have added a new dimension in the efficiency of identifying alien introgressions in wheat. Expressed sequence tag-sequence tagged sites (EST-STS) markers have proved useful in tracing alien chromatin. In this study, we report the development of Thinopyrum bessarabicum- and Secale anatolicum-specific EST-STS markers and their application in tracing respective alien chromatin introgressions in wheat. The parental lines, Chinese Spring (CS), ISR991.1 (CS/Th. bessarabicum amphidiploid), and ISR1049.2 (CS/Secale anatolicum amphidiploid), were used as core experimental materials. Using comparative analysis of RNA-Seq data, 10 903 and 10 660 candidate sequences specific to Th. bessarabicum and S. anatolicum, respectively, were assembled and identified. To validate the genome specificity of these candidate sequences, 68 and 64 EST-STS markers were developed from randomly selected candidate sequences of Th. bessarabicum and S. anatolicum, respectively, and tested on sets of alien addition lines. Fifty-five and 53 markers for Th. bessarabicum and S. anatolicum chromatin, respectively, were assigned to chromosomal location(s), covering all seven chromosomes. Approximately 83% of S. anatolicum-specific markers were transferable to S. cereale. The genome-specific candidate sequences identified and the EST-STS markers developed will be valuable resources for exploitation of Th. bessarabicum and Secale species diversity in wheat and triticale breeding.

Appraisal of wheat genomics for gene discovery and breeding applications: a special emphasis on advances in Asia

We discussed the most recent efforts in wheat functional genomics to discover new genes and their deployment in breeding with special emphasis on advances in Asian countries. Wheat research community is making significant progress to bridge genotype-to-phenotype gap and then applying this knowledge in genetic improvement. The advances in genomics and phenomics have intrigued wheat researchers in Asia to make best use of this knowledge in gene and trait discovery. These advancements include, but not limited to, map-based gene cloning, translational genomics, gene mapping, association genetics, gene editing and genomic selection. We reviewed more than 57 homeologous genes discovered underpinning important traits and multiple strategies used for their discovery. Further, the complementary advancements in wheat phenomics and analytical approaches to understand the genetics of wheat adaptability, resilience to climate extremes and resistance to pest and diseases were discussed. The challenge to build a gold standard reference genome sequence of bread wheat is now achieved and several de novo reference sequences from the cultivars representing different gene pools will be available soon. New pan-genome sequencing resources of wheat will strengthen the foundation required for accelerated gene discovery and provide more opportunities to practice the knowledge-based breeding.

Phenotypic effects of the U-genome variation in nascent synthetic hexaploids derived from interspecific crosses between durum wheat and its diploid relative Aegilops umbellulata

Aegilops umbellulata is a wild diploid wheat species with the UU genome that is an important genetic resource for wheat breeding. To exploit new synthetic allohexaploid lines available as bridges for wheat breeding, a total of 26 synthetic hexaploid lines were generated through crossing between the durum wheat cultivar Langdon and 26 accessions of Ae. umbellulata. In nascent synthetic hexaploids with the AABBUU genome, the presence of the set of seven U-genome chromosomes was confirmed with U-genome chromosome-specific markers developed based on RNA-seq-derived data from Ae. umbellulata. The AABBUU synthetic hexaploids showed large variations in flowering- and morphology-related traits, and these large variations transmitted well from the parental Ae. umbellulata accessions. However, the variation ranges in most traits examined were reduced under the AABBUU hexaploid background compared with under the diploid parents. The AABBUU and AABBDD synthetic hexaploids were clearly discriminated by several morphological traits, and an increase of plant height and in the number of spikes and a decrease of spike length were commonly observed in the AABBUU synthetics. Thus, interspecific differences in several morphological traits between Ae. umbellulata and A. tauschii largely affected the basic plant architecture of the synthetic hexaploids. In conclusion, the AABBUU synthetic hexaploid lines produced in the present study are useful resources for the introgression of desirable genes from Ae. umbellulata to common wheat.

Diploid genome differentiation conferred by RNA sequencing-based survey of genome-wide polymorphisms throughout homoeologous loci in Triticum and Aegilops

Background

Triticum and Aegilops diploid species have morphological and genetic diversity and are crucial genetic resources for wheat breeding. According to the chromosomal pairing-affinity of these species, their genome nomenclatures have been defined. However, evaluations of genome differentiation based on genome-wide nucleotide variations are still limited, especially in the three genomes of the genus Aegilops: Ae. caudata L. (CC genome), Ae. comosa Sibth. et Sm. (MM genome), and Ae. uniaristata Vis. (NN genome). To reveal the genome differentiation of these diploid species, we first performed RNA-seq-based polymorphic analyses for C, M, and N genomes, and then expanded the analysis to include the 12 diploid species of Triticum and Aegilops.

Results

Genetic divergence of the exon regions throughout the entire chromosomes in the M and N genomes was larger than that between A- and A^m-genomes. Ae. caudata had the second highest genetic diversity following Ae. speltoides, the putative B genome donor of common wheat. In the phylogenetic trees derived from the nuclear and chloroplast genome-wide polymorphism data, the C, D, M, N, U, and S genome species were connected with short internal branches, suggesting that these diploid species emerged during a relatively short period in the evolutionary process. The highly consistent nuclear and chloroplast phylogenetic topologies indicated that nuclear and chloroplast genomes of the diploid Triticum and Aegilops species coevolved after their diversification into each genome, accounting for most of the genome differentiation among the diploid species.

Conclusions

RNA-sequencing-based analyses successfully evaluated genome differentiation among the diploid Triticum and Aegilops species and supported the chromosome-pairing-based genome nomenclature system, except for the position of Ae. speltoides. Phylogenomic and epigenetic analyses of intergenic and centromeric regions could be essential for clarifying the mechanisms behind this inconsistency.

Genotyping-by-Sequencing Derived Single Nucleotide Polymorphisms Provide the First Well-Resolved Phylogeny for the Genus Triticum (Poaceae)

Wheat (Triticum spp.) has been an important staple food crop for mankind since the beginning of agriculture. The genus Triticum L. is composed of diploid, tetraploid, and hexaploid species, majority of which have not yet been discriminated clearly, and hence their phylogeny and classification remain unresolved. Genotyping-by-sequencing (GBS) is an easy and affordable method that allows us to generate genome-wide single nucleotide polymorphism (SNP) markers. In this study, we used GBS to obtain SNPs covering all seven chromosomes from 283 accessions of Triticum-related genera. After filtering low-quality and redundant SNPs based on haplotype information, the GBS assay provided 14,188 high-quality SNPs that were distributed across the A (71%), B (26%), and D (2.4%) genomes. Cluster analysis and discriminant analysis of principal components (DAPC) allowed us to distinguish six distinct groups that matched well with Triticum species complexity. We constructed a Bayesian phylogenetic tree using 14,188 SNPs, in which 17 Triticum species and subspecies were discriminated. Dendrogram analysis revealed that the polyploid wheat species could be divided into groups according to the presence of A, B, D, and G genomes with strong nodal support and provided new insight into the evolution of spelt wheat. A total of 2,692 species-specific SNPs were identified to discriminate the common (T. aestivum) and durum (T. turgidum) wheat cultivar and landraces. In principal component analysis grouping, the two wheat species formed individual clusters and the SNPs were able to distinguish up to nine groups of 10 subspecies. This study demonstrated that GBS-derived SNPs could be used efficiently in genebank management to classify Triticum species and subspecies that are very difficult to distinguish by their morphological characters.

Characterization of genetic diversity and population structure in wheat using array based SNP markers

Genetic diversity is crucial for successful adaptation and sustained improvement in crops. India is bestowed with diverse agro-climatic conditions which makes it rich in wheat germplasm adapted to various niches. Germplasm repository consists of local landraces, trait specific genetic stocks including introgressions from wild relatives, exotic collections, released varieties, and improved germplasm. Characterization of genetic diversity is done using morpho-physiological characters as well as by analyzing variations at DNA level. However, there are not many reports on array based high throughput SNP markers having characteristics of genome wide coverage employed in Indian spring wheat germplasm. Amongst wheat SNP arrays, 35K Axiom Wheat Breeder's Array has the highest SNP polymorphism efficiency suitable for genetic mapping and genetic diversity characterization. Therefore, genotyping was done using 35K in 483 wheat genotypes resulting in 14,650 quality filtered SNPs, that were distributed across the B (~ 50%), A (~ 39%), and D (~ 10%) genomes. The total genetic distance coverage was 4477.85 cM with 3.27 SNP/cM and 0.49 cM/SNP as average marker density and average inter-marker distance, respectively. The PIC ranged from 0.09 to 0.38 with an average of 0.29 across genomes. Population structure and Principal Coordinate Analysis resulted in two subpopulations (SP1 and SP2). The analysis of molecular variance revealed the genetic variation of 2% among and 98% within subpopulations indicating high gene flow between SP1 and SP2. The subpopulation SP2 showed high level of genetic diversity based on genetic diversity indices viz. Shannon's information index (I) = 0.648, expected heterozygosity (He) = 0.456 and unbiased expected heterozygosity (uHe) = 0.456. To the best of our knowledge, this study is the first to include the largest set of Indian wheat genotypes studied exclusively for genetic diversity. These findings may serve as a potential source for the identification of uncharacterized QTL/gene using genome wide association studies and marker assisted selection in wheat breeding programs.

Genetic Contribution of Synthetic Hexaploid Wheat to CIMMYT's Spring Bread Wheat Breeding Germplasm

Synthetic hexaploid (SH) wheat (AABBD'D') is developed by artificially generating a fertile hybrid between tetraploid durum wheat (Triticum turgidum, AABB) and diploid wild goat grass (Aegilops tauschii, D'D'). Over three decades, the International Maize and Wheat Improvement Center (CIMMYT) has developed and utilized SH wheat to bridge gene transfer from Ae. tauschii and durum wheat to hexaploid bread wheat. This is a unique example of success utilizing wild relatives in mainstream breeding at large scale worldwide. Our study aimed to determine the genetic contribution of SH wheat to CIMMYT's global spring bread wheat breeding program. We estimated the theoretical and empirical contribution of D' to synthetic derivative lines using the ancestral pedigree and marker information using over 1,600 advanced lines and their parents. The average marker-estimated D' contribution was 17.5% with difference in genome segments suggesting application of differential selection pressure. The pedigree-based contribution was correlated with marker-based estimates without providing chromosome segment specific variation. Results from international yield trials showed that 20% of the lines were synthetic derived with an average D' contribution of 15.6%. Our results underline the importance of SH wheat in maintaining and enhancing genetic diversity and genetic gain over years and is important for development of a more targeted introgression strategy. The study provides retrospective view into development and utilization of SH in the CIMMYT Global Wheat Program.

Genetic diversity and population structure analysis of synthetic and bread wheat accessions in Western Siberia

Recurrent selection and intercrossing between best of the best parents in each generation of breeding cycle resulted in a narrower genetic diversity in elite wheat (Triticum aestivum L.) germplasm. Therefore, we investigated diverse source of 143 synthetic and bread wheat accessions for identifying potentially rich genetic resources for improving the genetic diversity in wheat. This study identified 47,526 genotyping-by-sequencing-derived SNP markers that were nearly evenly distributed across three genomes of wheat. The population structure analysis identified three distinct clusters (Japan synthetics, CIMMYT synthetics, and bread wheat) of wheat genotypes on the basis of type and geographical origin of wheat accessions. Population differentiation using analysis of molecular variance indicated 21% of the total genetic variance among subgroups and the remainder within subgroups. This study also identified that the Japan synthetic group was the most divergent group compared with other subgroups. The genetic diversity comparisons between synthetic and bread wheat lines showed that the gene diversity of synthetic wheat was 33% higher than bread wheat accessions, indicating the potential use of these lines for broadening the genetic diversity of modern wheat cultivars. The results from this study will be helpful in further understanding genomic features of wheat and facilitate their use in wheat breeding programs.

Genome-wide polymorphisms from RNA sequencing assembly of leaf transcripts facilitate phylogenetic analysis and molecular marker development in wild einkorn wheat

A survey of genome-wide polymorphisms between closely related species is required to understand the molecular basis of the evolutionary differentiation of their genomes. Two wild diploid wheat species, namely Triticum monococcum ssp. aegilopoides and T. urartu, are closely related and harbour the A^m and A genomes, respectively. The A-genome donor of tetraploid and common wheat is T. urartu, and T. monococcum ssp. monococcum is the cultivated form derived from the wild einkorn wheat subspecies aegilopoides. Although subspecies aegilopoides has been a useful genetic resource in wheat breeding, genome-wide molecular markers for this subspecies have not been sufficiently developed. Here, we describe the detection of genome-wide polymorphisms such as single-nucleotide polymorphisms (SNPs) and insertions/deletions (indels) from RNA sequencing (RNA-seq) data of leaf transcripts in 15 accessions of the two diploid wheat species. The SNPs and indels, detected using the A genome of common wheat as the reference genome, covered the entire chromosomes of these species. The polymorphism information facilitated a comparison of the genetic diversity of einkorn wheat with that of two related diploid Aegilops species, namely, Ae. tauschii and Ae. umbellulata. Cleaved amplified polymorphic sequence (CAPS) markers converted from the SNP data were efficiently developed to confirm the addition of aegilopoides subspecies chromosomes to tetraploid wheat in nascent allohexaploid lines with AABBA^mA^m genomes. In addition, the CAPS markers permitted linkage map construction in mapping populations of aegilopoides subspecies accessions. Therefore, these RNA-seq data provide information for further breeding of closely related species with no reference genome sequence data.

Origin of wheat B-genome chromosomes inferred from RNA sequencing analysis of leaf transcripts from section Sitopsis species of Aegilops

Dramatic changes occasionally occur in intergenic regions leading to genomic alterations during speciation and will consequently obscure the ancestral species that have contributed to the formation of allopolyploid organisms. The S genome of five species of section Sitopsis of genus Aegilops is considered to be an origin of B-genome in cultivated tetraploid and hexaploid wheat species, although its actual donor is still unclear. Here, we attempted to elucidate phylogenetic relationship among Sitopsis species by performing RNA sequencing of the coding regions of each chromosome. Thus, genome-wide polymorphisms were extensively analyzed in 19 accessions of the Sitopsis species in reference to the tetraploid and hexaploid wheat B genome sequences and consequently were efficiently anchored to the B-genome chromosomes. The results of our genome-wide exon sequencing and resultant phylogenetic analysis indicate that Ae. speltoides is likely to be the direct donor of all chromosomes of the wheat B genome. Our results also indicate that the genome differentiation during wheat allopolyploidization from S to B proceeds at different speeds over the chromosomes rather than at constant rate and recombination could be a factor determining the speed. This observation is potentially generalized to genome differentiation during plant allopolyploid evolution.

RNA Sequencing-Based Bulked Segregant Analysis Facilitates Efficient D-genome Marker Development for a Specific Chromosomal Region of Synthetic Hexaploid Wheat

Common wheat originated from interspecific hybridization between cultivated tetraploid wheat and its wild diploid relative Aegilops tauschii followed by amphidiploidization. This evolutionary process can be reproduced artificially, resulting in synthetic hexaploid wheat lines. Here we performed RNA sequencing (RNA-seq)-based bulked segregant analysis (BSA) using a bi-parental mapping population of two synthetic hexaploid wheat lines that shared identical A and B genomes but included with D-genomes of distinct origins. This analysis permitted identification of D-genome-specific polymorphisms around the Net2 gene, a causative locus to hybrid necrosis. The resulting single nucleotide polymorphisms (SNPs) were classified into homoeologous polymorphisms and D-genome allelic variations, based on the RNA-seq results of a parental tetraploid and two Ae. tauschii accessions. The difference in allele frequency at the D-genome-specific SNP sites between the contrasting bulks (ΔSNP-index) was higher on the target chromosome than on the other chromosomes. Several SNPs with the highest ΔSNP-indices were converted into molecular markers and assigned to the Net2 chromosomal region. These results indicated that RNA-seq-based BSA can be applied efficiently to a synthetic hexaploid wheat population to permit molecular marker development in a specific chromosomal region of the D genome.

RNA-seq analysis reveals considerable genetic diversity and provides genetic markers saturating all chromosomes in the diploid wild wheat relative Aegilops umbellulata

BACKGROUND

Aegilops umbellulata Zhuk. (2n = 14), a wild diploid wheat relative, has been the source of trait improvement in wheat breeding. Intraspecific genetic variation of Ae. umbellulata, however, has not been well studied and the genomic information in this species is limited.

RESULTS

To develop novel genetic markers distributed over all chromosomes of Ae. umbellulata and to evaluate its genetic diversity, we performed RNA sequencing of 12 representative accessions and reconstructed transcripts by de novo assembly of reads for each accession. A large number of single nucleotide polymorphisms (SNPs) and insertions/deletions (indels) were obtained and anchored to the pseudomolecules of Ae. tauschii and barley (Hordeum vulgare L.), which were regarded as virtual chromosomes of Ae. umbellulata. Interestingly, genetic diversity in Ae. umbellulata was higher than in Ae. tauschii, despite the narrow habitat of Ae. umbellulata. Comparative analyses of nucleotide polymorphisms between Ae. umbellulata and Ae. tauschii revealed no clear lineage differentiation and existence of alleles with rarer frequencies predominantly in Ae. umbellulata, with patterns clearly distinct from those in Ae. tauschii.

CONCLUSIONS

The anchored SNPs, covering all chromosomes, provide sufficient genetic markers between Ae. umbellulata accessions. The alleles with rarer frequencies might be the main source of the high genetic diversity in Ae. umbellulata.

Unlocking the novel genetic diversity and population structure of synthetic Hexaploid wheat

Background

Synthetic hexaploid wheat (SHW) is a reconstitution of hexaploid wheat from its progenitors (Triticum turgidum ssp. durum L.; AABB x Aegilops tauschii Coss.; DD) and has novel sources of genetic diversity for broadening the genetic base of elite bread wheat (BW) germplasm (T. aestivum L). Understanding the diversity and population structure of SHWs will facilitate their use in wheat breeding programs. Our objectives were to understand the genetic diversity and population structure of SHWs and compare the genetic diversity of SHWs with elite BW cultivars and demonstrate the potential of SHWs to broaden the genetic base of modern wheat germplasm.

Results

The genotyping-by-sequencing of SHW provided 35,939 high-quality single nucleotide polymorphisms (SNPs) that were distributed across the A (33%), B (36%), and D (31%) genomes. The percentage of SNPs on the D genome was nearly same as the other two genomes, unlike in BW cultivars where the D genome polymorphism is generally much lower than the A and B genomes. This indicates the presence of high variation in the D genome in the SHWs. The D genome gene diversity of SHWs was 88.2% higher than that found in a sample of elite BW cultivars. Population structure analysis revealed that SHWs could be separated into two subgroups, mainly differentiated by geographical location of durum parents and growth habit of the crop (spring and winter type). Further population structure analysis of durum and Ae. parents separately identified two subgroups, mainly based on type of parents used. Although Ae. tauschii parents were divided into two sub-species: Ae. tauschii ssp. tauschii and ssp. strangulate, they were not clearly distinguished in the diversity analysis outcome. Population differentiation between SHWs (Spring_SHW and Winter_SHW) samples using analysis of molecular variance indicated 17.43% of genetic variance between populations and the remainder within populations.

Conclusions

SHWs were diverse and had a clearly distinguished population structure identified through GBS-derived SNPs. The results of this study will provide valuable information for wheat genetic improvement through inclusion of novel genetic variation and is a prerequisite for association mapping and genomic selection to unravel economically important marker-trait associations and for cultivar development.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4969-2) contains supplementary material, which is available to authorized users.

Adaptive selection of founder segments and epistatic control of plant height in the MAGIC winter wheat population WM-800

BACKGROUND

Multi-parent advanced generation intercross (MAGIC) populations are a newly established tool to dissect quantitative traits. We developed the high resolution MAGIC wheat population WM-800, consisting of 910 F4:6 lines derived from intercrossing eight recently released European winter wheat cultivars.

RESULTS

Genotyping WM-800 with 7849 SNPs revealed a low mean genetic similarity of 59.7% between MAGIC lines. WM-800 harbours distinct genomic regions exposed to segregation distortion. These are mainly located on chromosomes 2 to 6 of the wheat B genome where founder specific DNA segments were positively or negatively selected. This suggests adaptive selection of individual founder alleles during population development. The application of a genome-wide association study identified 14 quantitative trait loci (QTL) controlling plant height in WM-800, including the known semi-dwarf genes Rht-B1 and Rht-D1 and a potentially novel QTL on chromosome 5A. Additionally, epistatic effects controlled plant height. For example, two loci on chromosomes 2B and 7B gave rise to an additive epistatic effect of 13.7 cm.

CONCLUSION

The present study demonstrates that plant height in the MAGIC-WHEAT population WM-800 is mainly determined by large-effect QTL and di-genic epistatic interactions. As a proof of concept, our study confirms that WM-800 is a valuable tool to dissect the genetic architecture of important agronomic traits.

Wheat genetic resources in the post-genomics era: promise and challenges

Background

Wheat genetic resources have been used for genetic improvement since 1876, when Stephen Wilson (Transactions and Proceedings of the Botanical Society of Edinburgh 12: 286) consciously made the first wide hybrid involving wheat and rye in Scotland. Wide crossing continued with sporadic attempts in the first half of 19th century and became a sophisticated scientific discipline during the last few decades with considerable impact in farmers' fields. However, a large diversity of untapped genetic resources could contribute in meeting future wheat production challenges.

Perspectives and Conclusion

Recently the complete reference genome of hexaploid (Chinese Spring) and tetraploid (Triticum turgidum ssp. dicoccoides) wheat became publicly available coupled with on-going international efforts on wheat pan-genome sequencing. We anticipate that an objective appraisal is required in the post-genomics era to prioritize genetic resources for use in the improvement of wheat production if the goal of doubling yield by 2050 is to be met. Advances in genomics have resulted in the development of high-throughput genotyping arrays, improved and efficient methods of gene discovery, genomics-assisted selection and gene editing using endonucleases. Likewise, ongoing advances in rapid generation turnover, improved phenotyping, envirotyping and analytical methods will significantly accelerate exploitation of exotic genes and increase the rate of genetic gain in breeding. We argue that the integration of these advances will significantly improve the precision and targeted identification of potentially useful variation in the wild relatives of wheat, providing new opportunities to contribute to yield and quality improvement, tolerance to abiotic stresses, resistance to emerging biotic stresses and resilience to weather extremes.

Genetic mapping reveals a dominant awn-inhibiting gene related to differentiation of the variety anathera in the wild diploid wheat Aegilops tauschii

Aegilops tauschii, a wild wheat relative, is the D-genome donor of common wheat. Subspecies and varieties of Ae. tauschii are traditionally classified based on differences in their inflorescence architecture. However, the genetic information for their diversification has been quite limited in the wild wheat relatives. The variety anathera has no awn on the lemma, but the genetic basis for this diagnostic character is unknown. Wide variations in awn length traits at the top and middle spikes were found in the Ae. tauschii core collection, and the awn length at the middle spike was significantly smaller in the eastward-dispersed sublineage than in those in other sublineages. To clarify loci controlling the awnless phenotype of var. anathera, we measured awn length of an intervariety F₂ mapping population, and found that the F₂ individuals could be divided into two groups mainly based on the awn length at the middle of spike, namely short and long awn groups, significantly fitting a 3:1 segregation ratio, which indicated that a single locus controls the awnless phenotype. The awnless locus, Anathera (Antr), was assigned to the distal region of the short arm of chromosome 5D. Quantitative trait locus analysis using the awn length data of each F₂ individual showed that only one major locus was at the same chromosomal position as Antr. These results suggest that a single dominant allele determines the awnless diagnostic character in the variety anathera. The Antr dominant allele is a novel gene inhibiting awn elongation in wheat and its relatives.

Genotyping-by-Sequencing (GBS) Revealed Molecular Genetic Diversity of Iranian Wheat Landraces and Cultivars

Background: Genetic diversity is an essential resource for breeders to improve new cultivars with desirable characteristics. Recently, genotyping-by-sequencing (GBS), a next-generation sequencing (NGS) technology that can simplify complex genomes, has now be used as a high-throughput and cost-effective molecular tool for routine breeding and screening in many crop species, including the species with a large genome. Results: We genotyped a diversity panel of 369 Iranian hexaploid wheat accessions including 270 landraces collected between 1931 and 1968 in different climate zones and 99 cultivars released between 1942 to 2014 using 16,506 GBS-based single nucleotide polymorphism (GBS-SNP) markers. The B genome had the highest number of mapped SNPs while the D genome had the lowest on both the Chinese Spring and W7984 references. Structure and cluster analyses divided the panel into three groups with two landrace groups and one cultivar group, suggesting a high differentiation between landraces and cultivars and between landraces. The cultivar group can be further divided into four subgroups with one subgroup was mostly derived from Iranian ancestor(s). Similarly, landrace groups can be further divided based on years of collection and climate zones where the accessions were collected. Molecular analysis of variance indicated that the genetic variation was larger between groups than within group. Conclusion: Obvious genetic diversity in Iranian wheat was revealed by analysis of GBS-SNPs and thus breeders can select genetically distant parents for crossing in breeding. The diverse Iranian landraces provide rich genetic sources of tolerance to biotic and abiotic stresses, and they can be useful resources for the improvement of wheat production in Iran and other countries.

Candidate loci involved in domestication and improvement detected by a published 90K wheat SNP array

Selection is one of the most important forces in crop evolution. Common wheat is a major world food crop and a typical allopolyploid with a huge and complex genome. We applied four approaches to detect loci selected in wheat during domestication and improvement. A total of 7,984 candidate loci were detected, accounting for 23.3% of all 34,317 SNPs analysed, a much higher proportion than estimated in previous reports. We constructed a first generation wheat selection map which revealed the following new insights on genome-wide selection: (1) diversifying selection acted by increasing, decreasing or not affecting gene frequencies; (2) the number of loci under selection during domestication was much higher than that during improvement; (3) the contribution to wheat improvement by the D sub-genome was relatively small due to the bottleneck of hexaploidisation and diversity can be expanded by using synthetic wheat and introgression lines; and (4) clustered selection regions occur throughout the wheat genome, including the centromere regions. This study will not only help future wheat breeding and evolutionary studies, but will also accelerate study of other crops, especially polyploids.

Fine mapping and genetic association analysis of Net2, the causative D-genome locus of low temperature-induced hybrid necrosis in interspecific crosses between tetraploid wheat and Aegilops tauschii

Hybrid necrosis has been observed in many interspecific hybrids from crosses between tetraploid wheat and the wheat D-genome donor Aegilops tauschii. Type II necrosis is a kind of hybrid incompatibility that is specifically characterized by low-temperature induction and growth suppression. Two complementary genes, Net1 on the AB genome and Net2 on the D genome, putatively control type II necrosis in ABD triploids and synthetic hexaploid wheat. Toward map-based cloning of Net2, a fine map around the Net2 region on 2DS was constructed in this study. Using the draft genome sequence of Ae. tauschii and the physical map of the barley genome, the Net2 locus was mapped within a 0.6 cM interval between two closely linked markers. Although local chromosomal rearrangements were observed in the Net2-corresponding region between the barley/Brachypodium and Ae. tauschii genomes, the two closely linked markers were significantly associated with type II necrosis in Ae. tauschii. These results suggest that these markers will aid efficient selection of Net2 non-carrier individuals from the Ae. tauschii population and intraspecific progeny, and could help with introgression of agriculturally important genes from Ae. tauschii to common wheat.

Quantitative trait locus analysis for flowering-related traits using two F2 populations derived from crosses between Japanese common wheat cultivars and synthetic hexaploids

Flowering time is an important trait for Japanese wheat breeding. Aegilops tauschii, the D-genome donor of hexaploid wheat, is a useful resource to enlarge the D-genome diversity of common wheat. Previously, we identified flowering-related QTLs in F2 populations of synthetic hexaploid wheat lines between the tetraploid wheat cultivar Langdon and Ae. tauschii accessions. Here, to evaluate the usefulness of the early-flowering alleles from Ae. tauschii for Japanese wheat breeding, QTL analyses were conducted in two F2 populations derived from crosses between Japanese wheat cultivars and early-flowering lines of synthetic hexaploid wheat. Only two chromosomal regions controlling flowering-related traits were identified, on chromosomes 2DS and 5AL in the mapping populations, and no previously identified QTLs were found in the synthetic hexaploid lines. The strong effect of the 2DS QTL, putatively corresponding to Ppd-D1, was considered to hide any significant expression of other QTLs with small effects on flowering-related traits. When F2 individuals carrying Ae. tauschii-homozygous alleles around the 2DS QTL region were selected, the Ae. tauschii-derived alleles of the previously identified flowering QTLs partly showed an early-flowering phenotype compared with the Japanese wheat-derived alleles. Thus, some early-flowering alleles from Ae. tauschii may be useful for production of early-flowering Japanese wheat cultivars.

Genomics of crop wild relatives: expanding the gene pool for crop improvement

Plant breeders require access to new genetic diversity to satisfy the demands of a growing human population for more food that can be produced in a variable or changing climate and to deliver the high-quality food with nutritional and health benefits demanded by consumers. The close relatives of domesticated plants, crop wild relatives (CWRs), represent a practical gene pool for use by plant breeders. Genomics of CWR generates data that support the use of CWR to expand the genetic diversity of crop plants. Advances in DNA sequencing technology are enabling the efficient sequencing of CWR and their increased use in crop improvement. As the sequencing of genomes of major crop species is completed, attention has shifted to analysis of the wider gene pool of major crops including CWR. A combination of de novo sequencing and resequencing is required to efficiently explore useful genetic variation in CWR. Analysis of the nuclear genome, transcriptome and maternal (chloroplast and mitochondrial) genome of CWR is facilitating their use in crop improvement. Genome analysis results in discovery of useful alleles in CWR and identification of regions of the genome in which diversity has been lost in domestication bottlenecks. Targeting of high priority CWR for sequencing will maximize the contribution of genome sequencing of CWR. Coordination of global efforts to apply genomics has the potential to accelerate access to and conservation of the biodiversity essential to the sustainability of agriculture and food production.

Application of Population Sequencing (POPSEQ) for Ordering and Imputing Genotyping-by-Sequencing Markers in Hexaploid Wheat

The advancement of next-generation sequencing technologies in conjunction with new bioinformatics tools enabled fine-tuning of sequence-based, high-resolution mapping strategies for complex genomes. Although genotyping-by-sequencing (GBS) provides a large number of markers, its application for association mapping and genomics-assisted breeding is limited by a large proportion of missing data per marker. For species with a reference genomic sequence, markers can be ordered on the physical map. However, in the absence of reference marker order, the use and imputation of GBS markers is challenging. Here, we demonstrate how the population sequencing (POPSEQ) approach can be used to provide marker context for GBS in wheat. The utility of a POPSEQ-based genetic map as a reference map to create genetically ordered markers on a chromosome for hexaploid wheat was validated by constructing an independent de novo linkage map of GBS markers from a Synthetic W7984 × Opata M85 recombinant inbred line (SynOpRIL) population. The results indicated that there is strong agreement between the independent de novo linkage map and the POPSEQ mapping approach in mapping and ordering GBS markers for hexaploid wheat. After ordering, a large number of GBS markers were imputed, thus providing a high-quality reference map that can be used for QTL mapping for different traits. The POPSEQ-based reference map and whole-genome sequence assemblies are valuable resources that can be used to order GBS markers and enable the application of highly accurate imputation methods to leverage the application GBS markers in wheat.

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.

The cuticular wax inhibitor locus Iw2 in wild diploid wheat Aegilops tauschii: phenotypic survey, genetic analysis, and implications for the evolution of common wheat

BACKGROUND

Cuticular wax production on plant surfaces confers a glaucous appearance and plays important roles in plant stress tolerance. Most common wheat cultivars, which are hexaploid, and most tetraploid wheat cultivars are glaucous; in contrast, a wild wheat progenitor, Aegilops tauschii, can be glaucous or non-glaucous. A dominant non-glaucous allele, Iw2, resides on the short arm of chromosome 2D, which was inherited from Ae. tauschii through polyploidization. Iw2 is one of the major causal genes related to variation in glaucousness among hexaploid wheat. Detailed genetic and phylogeographic knowledge of the Iw2 locus in Ae. tauschii may provide important information and lead to a better understanding of the evolution of common wheat.

RESULTS

Glaucous Ae. tauschii accessions were collected from a broad area ranging from Armenia to the southwestern coastal part of the Caspian Sea. Linkage analyses with five mapping populations showed that the glaucous versus non-glaucous difference was mainly controlled by the Iw2 locus in Ae. tauschii. Comparative genomic analysis of barley and Ae. tauschii was then used to develop molecular markers tightly linked with Ae. tauschii Iw2. Chromosomal synteny around the orthologous Iw2 regions indicated that some chromosomal rearrangement had occurred during the genetic divergence leading to Ae. tauschii, barley, and Brachypodium. Genetic associations between specific Iw2-linked markers and respective glaucous phenotypes in Ae. tauschii indicated that at least two non-glaucous accessions might carry other glaucousness-determining loci outside of the Iw2 locus.

CONCLUSION

Allelic differences at the Iw2 locus were the main contributors to the phenotypic difference between the glaucous and non-glaucous accessions of Ae. tauschii. Our results supported the previous assumption that the D-genome donor of common wheat could have been any Ae. tauschii variant that carried the recessive iw2 allele.

Identification of quantitative trait loci for abscisic acid responsiveness in the D-genome of hexaploid wheat

In crop species such as wheat, abiotic stresses and preharvest sprouting reduce grain yield and quality. The plant hormone abscisic acid (ABA) plays important roles in abiotic stress tolerance and seed dormancy. In previous studies, we evaluated ABA responsiveness of 67 Aegilops tauschii accessions and their synthetic hexaploid wheat lines, finding wide variation that was due to the D-genome. In this study, quantitative trait locus (QTL) analysis was performed using an F2 population derived from crosses of highly ABA-responsive and less-responsive synthetic wheat lines. A significant QTL was detected on chromosome 6D, in a similar location to that reported for ABA responsiveness using recombinant inbred lines derived from common wheat cultivars Mironovskaya 808 and Chinese Spring. A comparative map and physiological and expression analyses of the 6D QTL suggested that this locus involved in line differences among wheat synthetics is different from that involved in cultivar differences in common wheat. The common wheat 6D QTL was found to affect seed dormancy and the regulation of cold-responsive/late embryogenesis abundant genes during dehydration. However, in synthetic wheat, we failed to detect any association of ABA responsiveness with abiotic stress tolerance or seed dormancy, at least under our experimental conditions. Development of near-isogenic lines will be important for functional analyses of the synthetic wheat 6D QTL.