The chloroplast view of the evolution of polyploid wheat

Phylogeny and taxonomy of the polyploid species that contain St genome (Triticeae; Poaceae) based on four nuclear DNA and three chloroplast genes

BACKGROUND

The genus Pseudoroegneria (Nevski) Á.Löve contributes the St genome for more than 60% of perennial Triticeae species. However, the strong dominant character of the St genome makes it challenging to distinguish each species and/or even genus based on single or combined morphological traits. Moreover, the phylogeny and taxonomy of the St-genome containing polyploid genera remain controversial.

RESULTS

In this study, we used nuclear and chloroplast DNA-based phylogenetic analyses to reveal the systematic relationships between the St-genome containing polyploid species. The maximum likelihood (ML) tree based on nuclear ribosomal internal transcribed spacer region (nrITS) and three single-copy nuclear genes data (Acc1 + Pgk1 + DMC1) showed that polyploid species with the St genome were separated into seven genera with StStHH, StStYY, StStYYHH, StStYYPP, StStYYWW, StStPP, and StStEE genome constitutions, moreover, the polyploid species in Caucasus, America, and Australia have independent polyploidization events. The ML tree for the chloroplast DNA fragments (matK + rbcL + trnL-trnF) displayed that the P genome served as a maternal donor of Kengyilia melanthera and K. dingqinensis from the Hengduan Mountains region, while the St or StY genome served as the maternal donor of other St-genome containing species. Herein, we reported the genomic constitution of Kengyilia tibetica, K. changduensis, and K. dingqinensis with the StStYYPP genome for the first time.

CONCLUSIONS

The St-genome-containing polyploid species should be treated as distinct genera according to different genome constitutions, and those species experienced independent allo-polyploidization events in different distribution regions and had two relatively independent maternal origins from the P or St/StY genomes. Besides, the Xp genome might have contributed to the unknown Y genome formation.

Characterization of the complete chloroplast genome of Chasmanthium latifolium (Michx.) H.O.Yates, 1966 (Poaceae)

Chasmanthium latifolium (Michx.) H.O.Yates is a popular ornamental plant native to southeastern North America. Genomic data and genetic studies related to Chasmanthium latifolium are limited. Therefore, the complete chloroplast genome of Chasmanthium latifolium was sequenced, assembled, and characterized in this study. The complete chloroplast genome was 138,934 bp in length and contained 105 unique genes (77 protein-coding genes, 24 tRNA genes, and 4 rRNA genes). Phylogenetic analyses showed that Chasmanthium latifolium and Chasmanthium laxum clustered into a separate clade with the closest affinity to the clade comprising Zeugites pittieri Hack and Lophatherum gracile Brongn. In conclusion, our study describes the complete chloroplast genome of Chasmanthium latifolium for the first time, contributing to a better understanding of its taxonomy and evolution.

Identifying novel sources of resistance to wheat stem sawfly in five wild wheat species

BACKGROUND

The wheat stem sawfly (WSS, Cephus cinctus) is a major pest of wheat (Triticum aestivum) and can cause significant yield losses. WSS damage results from stem boring and/or cutting, leading to the lodging of wheat plants. Although solid-stem wheat genotypes can effectively reduce larval survival, they may have lower yields than hollow-stem genotypes and show inconsistent solidness expression. Because of limited resistance sources to WSS, evaluating diverse wheat germplasm for novel resistance genes is crucial. We evaluated 91 accessions across five wild wheat species (Triticum monococcum, T. urartu, T. turgidum, T. timopheevii, and Aegilops tauschii) and common wheat cultivars (T. aestivum) for antixenosis (host selection) and antibiosis (host suitability) to WSS. Host selection was measured as the number of eggs after adult oviposition, and host suitability was determined by examining the presence or absence of larval infestation within the stem. The plants were grown in the greenhouse and brought to the field for WSS infestation. In addition, a phylogenetic analysis was performed to determine the relationship between the WSS traits and phylogenetic clustering.

RESULTS

Overall, Ae. tauschii, T. turgidum and T. urartu had lower egg counts and larval infestation than T. monococcum, and T. timopheevii. T. monococcum, T. timopheevii, T. turgidum, and T. urartu had lower larval weights compared with T. aestivum.

CONCLUSION

This study shows that wild relatives of wheat could be a valuable source of alleles for enhancing resistance to WSS and identifies specific germplasm resources that may be useful for breeding. © 2024 The Authors. Pest Management Science published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

Genomic insights into the origin and evolution of spelt (Triticum spelta L.) as a valuable gene pool for modern wheat breeding

Spelt (Triticum aestivum ssp. spelta) is an important wheat subspecies mainly cultivated in Europe before the 20th century that has contributed to modern wheat breeding as a valuable genetic resource. However, relatively little is known about the origins and maintenance of spelt populations. Here, using resequencing data from 416 worldwide wheat accessions, including representative spelt wheat, we demonstrate that European spelt emerged when primitive hexaploid wheat spread to the west and hybridized with pre-settled domesticated emmer, the putative maternal donor. Genomic introgression regions from domesticated emmer confer spelt's primitive morphological characters used for species taxonomy, such as tenacious glumes and later flowering. We propose a haplotype-based "spelt index" to identify spelt-type wheat varieties and to quantify utilization of the spelt gene pool in modern wheat cultivars. This study reveals the genetic basis for the establishment of the spelt wheat subspecies in a specific ecological niche and the vital role of the spelt gene pool as a unique germplasm resource in modern wheat breeding.

Exploring genetic diversity of wild and related tetraploid wheat species Triticum turgidum and Triticum timopheevii

INTRODUCTION

The domestication bottleneck has reduced genetic diversity inwheat, necessitating the use of wild relatives in breeding programs. Wild tetraploid wheat are widely used in the breeding programs but with morphological characters, it is difficult to distinguish these, resulting in misclassification/mislabeling or duplication of accessions in the Gene bank.

OBJECTIVES

The study aims to exploreGenotyping by sequencing (GBS) to characterize wild and domesticated tetraploid wheat accessions to generate a core set of accessions to be used in the breeding program.

METHODS

TASSEL-GBS pipeline was used for SNP discovery, fastStructure was used to determine the population structure and PowerCore was used to generate a core sets. Nucleotide diversity matrices of Nie's and F-statistics (FST) index were used to determine the center of genetic diversity.

RESULTS

We found 65 % and 47 % duplicated accessions in Triticum timopheevii and T. turgidum respectively. Genome-wide nucleotide diversity and FST scan uncovered a lower intra and higher inter-species differentiation. Distinct FST regions were identified in genomic regions belonging to domestication genes: non-brittle rachis (Btr1) and vernalization (VRN-1).Our results suggest that Israel, Jordan, Syria, and Lebanonas the hub of genetic diversity of wild emmer;Turkey, and Georgia for T. durum; and Iraq, Azerbaijan, and Armenia for theT. timopheevii. Identified core set accessions preserved more than 93 % of the available genetic diversity. Genome wide association study (GWAS) indicated the potential chromosomal segment for resistance to leaf rust in T. timopheevii.

CONCLUSION

The present study explored the potential of GBS technology in data reduction while maintaining the significant genetic diversity of the species. Wild germplasm showed more differentiation than domesticated accessions, indicating the availability of sufficient diversity for crop improvement. With reduced complexity, the core set preserves the genetic diversity of the gene bank collections and will aid in a more robust characterization of wild germplasm.

G-SAIP: Graphical Sequence Alignment Through Parallel Programming in the Post-Genomic Era

A common task in bioinformatics is to compare DNA sequences to identify similarities between organisms at the sequence level. An approach to such comparison is the dot-plots, a 2-dimensional graphical representation to analyze DNA or protein alignments. Dot-plots alignment software existed before the sequencing revolution, and now there is an ongoing limitation when dealing with large-size sequences, resulting in very long execution times. High-Performance Computing (HPC) techniques have been successfully used in many applications to reduce computing times, but so far, very few applications for graphical sequence alignment using HPC have been reported. Here, we present G-SAIP (Graphical Sequence Alignment in Parallel), a software capable of spawning multiple distributed processes on CPUs, over a supercomputing infrastructure to speed up the execution time for dot-plot generation up to 1.68× compared with other current fastest tools, improve the efficiency for comparative structural genomic analysis, phylogenetics because the benefits of pairwise alignments for comparison between genomes, repetitive structure identification, and assembly quality checking.

The grain quality of wheat wild relatives in the evolutionary context

KEY MESSAGE

We evaluated the potential of wheat wild relatives for the improvement in grain quality characteristics including micronutrients (Fe, Zn) and gluten and identified diploid wheats and the timopheevii lineage as the most promising resources. Domestication enabled the advancement of civilization through modification of plants according to human requirements. Continuous selection and cultivation of domesticated plants induced genetic bottlenecks. However, ancient diversity has been conserved in crop wild relatives. Wheat (Triticum aestivum L.; Triticum durum Desf.) is one of the most important staple foods and was among the first domesticated crop species. Its evolutionary diversity includes diploid, tetraploid and hexaploid species from the Triticum and Aegilops taxa and different genomes, generating an AA, BBAA/GGAA and BBAADD/GGAAAmAm genepool, respectively. Breeding and improvement in wheat altered its grain quality. In this review, we identified evolutionary patterns and the potential of wheat wild relatives for quality improvement regarding the micronutrients Iron (Fe) and Zinc (Zn), the gluten storage proteins α-gliadins and high molecular weight glutenin subunits (HMW-GS), and the secondary metabolite phenolics. Generally, the timopheevii lineage has been neglected to date regarding grain quality studies. Thus, the timopheevii lineage should be subject to grain quality research to explore the full diversity of the wheat gene pool.

A phylogenetic analysis of Bromus (Poaceae: Pooideae: Bromeae) based on nuclear ribosomal and plastid data, with a focus on Bromus sect. Bromus

To investigate phylogenetic relationships among and within major lineages of Bromus, with focus on Bromus sect. Bromus, we analyzed DNA sequences from two nuclear ribosomal (ITS, ETS) and two plastid (rpl32-trnL_UAG , matK) regions. We sampled 103 ingroup accessions representing 26 taxa of B. section Bromus and 15 species of other Bromus sections. Our analyses confirm the monophyly of Bromus s.l. and identify incongruence between nuclear ribosomal and plastid data partitions for relationships within and among major Bromus lineages. Results support classification of B. pumilio and B. gracillimus within B. sect. Boissiera and B. sect. Nevskiella, respectively. These species are sister groups and are closely related to B. densus (B. sect. Mexibromus) in nrDNA trees and Bromus sect. Ceratochloa in plastid trees. Bromus sect. Bromopsis is paraphyletic. In nrDNA trees, species of Bromus sects. Bromopsis, Ceratochloa, Neobromus, and Genea plus B. rechingeri of B. sect. Bromus form a clade, in which B. tomentellus is sister to a B. sect. Genea-B. rechingeri clade. In the plastid trees, by contrast, B. sect. Bromopsis species except B. tomentosus form a clade, and B. tomentosus is sister to a clade comprising B. sect. Bromus and B. sect. Genea species. Affinities of B. gedrosianus, B. pulchellus, and B. rechingeri (members of the B. pectinatus complex), as well as B. oxyodon and B. sewerzowii, are discordant between nrDNA and plastid trees. We infer these species may have obtained their plastomes via chloroplast capture from species of B. sect. Bromus and B. sect. Genea. Within B. sect. Bromus, B. alopecuros subsp. caroli-henrici, a clade comprising B. hordeaceus and B. interruptus, and B. scoparius are successive sister groups to the rest of the section in the nrDNA phylogeny. Most relationships among the remaining species of B. sect. Bromus are unresolved in the nrDNA and plastid trees. Given these results, we infer that most B. sect. Bromus species likely diversified relatively recently. None of the subdivisional taxa proposed for Bromus sect. Bromus over the last century correspond to natural groups identified in our phylogenetic analyses except for a group including B. hordeaceus and B. interruptus.

Comparative and Phylogenetic Analysis of Complete Chloroplast Genomes in Leymus (Triticodae, Poaceae)

Leymus is a perennial genus that belongs to the tribe Triticeae (Poaceae) which has an adaptive capacity to ecological conditions and strong resistance to cold, drought, and salinity. Most Leymus species are fine herbs that can be used for agriculture, conservation, and landscaping. Due to confusion taxonomy within genera, the complete chloroplast (cp) genome of 13 Leymus species was sequenced, assembled, and compared with those of three other previously published Leymus species (Leymus condensatus, Leymus angustus, and Leymus mollis) to clarify the issue. Overall, the whole cp genome size ranged between 135,057 (L. condensatus) and 136,906 bp (Leymus coreanus) and showed a typical quadripartite structure. All studied species had 129 genes, including 83 protein-coding genes, 38 transfer RNAs, and 8 ribosomal RNAs. In total, 800 tandem repeats and 707 SSR loci were detected, most of which were distributed in the large single-copy region, followed by the inverted repeat (IR) and small single-copy regions. The sequence identity of all sequences was highly similar, especially concerning the protein-coding and IR regions; in particular, the protein-coding regions were significantly similar to those in the IR regions, regardless of small sequence differences in the whole cp genome. Moreover, the coding regions were more conserved than the non-coding regions. Comparisons of the IR boundaries showed that IR contraction and expansion events were reflected in different locations of rpl22, rps19, ndhH, and psbA genes. The close phylogenetic relationship of Leymus and Psathyrostachys indicated that Psathyrostachys possibly is the donor of the Ns genome sequence identified in Leymus. Altogether, the complete cp genome sequence of Leymus will lay a solid foundation for future population genetics and phylogeography studies, as well as for the analysis of the evolution of economically valuable plants.

Dissection of Structural Reorganization of Wheat 5B Chromosome Associated With Interspecies Recombination Suppression

Chromosomal rearrangements that lead to recombination suppression can have a significant impact on speciation, and they are also important for breeding. The regions of recombination suppression in wheat chromosome 5B were identified based on comparisons of the 5B map of a cross between the Chinese Spring (CS) variety of hexaploid wheat and CS-5Bdic (genotype CS with 5B substituted with its homologue from tetraploid Triticum dicoccoides) with several 5B maps of tetraploid and hexaploid wheat. In total, two regions were selected in which recombination suppression occurred in cross CS × CS-5Bdic when compared with other maps: one on the short arm, 5BS_RS, limited by markers BS00009810/BS00022336, and the second on the long arm, 5BL_RS, between markers Ra_c10633_2155 and BS00087043. The regions marked as 5BS_RS and 5BL_RS, with lengths of 5 Mb and 3.6 Mb, respectively, were mined from the 5B pseudomolecule of CS and compared to the homoeologous regions (7.6 and 3.8 Mb, respectively) of the 5B pseudomolecule of Zavitan (T. dicoccoides). It was shown that, in the case of 5BS_RS, the local heterochromatin islands determined by the satellite DNA (119.2) and transposable element arrays, as well as the dissimilarity caused by large insertions/deletions (chromosome rearrangements) between 5BSs aestivum/dicoccoides, are likely the key determinants of recombination suppression in the region. Two major and two minor segments with significant loss of similarity were recognized within the 5BL_RS region. It was shown that the loss of similarity, which can lead to suppression of recombination in the 5BL_RS region, is caused by chromosomal rearrangements, driven by the activity of mobile genetic elements (both DNA transposons and long terminal repeat retrotransposons) and their divergence during evolution. It was noted that the regions marked as 5BS_RS and 5BL_RS are associated with chromosomal rearrangements identified earlier by С-banding analysis of intraspecific polymorphism of tetraploid emmer wheat. The revealed divergence in 5BS_RS and 5BL_RS may be a consequence of interspecific hybridization, plant genetic adaptation, or both.

Global Patterns of Subgenome Evolution in Organelle-Targeted Genes of Six Allotetraploid Angiosperms

Whole-genome duplications (WGDs) are a prominent process of diversification in eukaryotes. The genetic and evolutionary forces that WGD imposes on cytoplasmic genomes are not well understood, despite the central role that cytonuclear interactions play in eukaryotic function and fitness. Cellular respiration and photosynthesis depend on successful interaction between the 3,000+ nuclear-encoded proteins destined for the mitochondria or plastids and the gene products of cytoplasmic genomes in multi-subunit complexes such as OXPHOS, organellar ribosomes, Photosystems I and II, and Rubisco. Allopolyploids are thus faced with the critical task of coordinating interactions between the nuclear and cytoplasmic genes that were inherited from different species. Because the cytoplasmic genomes share a more recent history of common descent with the maternal nuclear subgenome than the paternal subgenome, evolutionary "mismatches" between the paternal subgenome and the cytoplasmic genomes in allopolyploids might lead to the accelerated rates of evolution in the paternal homoeologs of allopolyploids, either through relaxed purifying selection or strong directional selection to rectify these mismatches. We report evidence from six independently formed allotetraploids that the subgenomes exhibit unequal rates of protein-sequence evolution, but we found no evidence that cytonuclear incompatibilities result in altered evolutionary trajectories of the paternal homoeologs of organelle-targeted genes. The analyses of gene content revealed mixed evidence for whether the organelle-targeted genes are lost more rapidly than the non-organelle-targeted genes. Together, these global analyses provide insights into the complex evolutionary dynamics of allopolyploids, showing that the allopolyploid subgenomes have separate evolutionary trajectories despite sharing the same nucleus, generation time, and ecological context.

Salinity Tolerance in a Synthetic Allotetraploid Wheat (SlSlAA) Is Similar to Its Higher Tolerant Parent Aegilops longissima (SlSl) and Linked to Flavonoids Metabolism

Allotetraploidization between A and S (closely related to B) genome species led to the speciation of allotetraploid wheat (genome BBAA). However, the immediate metabolic outcomes and adaptive changes caused by the allotetraploidization event are poorly understood. Here, we investigated how allotetraploidization affected salinity tolerance using a synthetic allotetraploid wheat line (genome S^lS^lAA, labeled as 4x), its Aegilops longissima (genome S^lS^l, labeled as S^lS^l) and Triticum urartu (AA genome, labeled as AA) parents. We found that the degree of salinity tolerance of 4x was similar to its S^lS^l parent, and both were substantially more tolerant to salinity stress than AA. This suggests that the S^lS^l subgenome exerts a dominant effect for this trait in 4x. Compared with S^lS^l and 4x, the salinity-stressed AA plants did not accumulate a higher concentration of Na⁺ in leaves, but showed severe membrane peroxidation and accumulated a higher concentration of ROS (H₂O₂ and O₂ ^⋅⁣-) and a lesser concentration of flavonoids, indicating that ROS metabolism plays a key role in saline sensitivity. Exogenous flavonoid application to roots of AA plants significantly relieved salinity-caused injury. Our results suggest that the higher accumulation of flavonoids in S^lS^l may contribute to ROS scavenging and salinity tolerance, and these physiological properties were stably inherited by the nascent allotetraploid S^lS^lAA.

Genome sequences of five Sitopsis species of Aegilops and the origin of polyploid wheat B subgenome

Common wheat (Triticum aestivum, BBAADD) is a major staple food crop worldwide. The diploid progenitors of the A and D subgenomes have been unequivocally identified; that of B, however, remains ambiguous and controversial but is suspected to be related to species of Aegilops, section Sitopsis. Here, we report the assembly of chromosome-level genome sequences of all five Sitopsis species, namely Aegilops bicornis, Ae. longissima, Ae. searsii, Ae. sharonensis, and Ae. speltoides, as well as the partial assembly of the Amblyopyrum muticum (synonym Aegilops mutica) genome for phylogenetic analysis. Our results reveal that the donor of the common wheat B subgenome is a distinct, and most probably extinct, diploid species that diverged from an ancestral progenitor of the B lineage to which the still extant Ae. speltoides and Am. muticum belong. In addition, we identified interspecific genetic introgressions throughout the evolution of the Triticum/Aegilops species complex. The five Sitopsis species have various assembled genome sizes (4.11-5.89 Gb) with high proportions of repetitive sequences (85.99%-89.81%); nonetheless, they retain high collinearity with other genomes or subgenomes of species in the Triticum/Aegilops complex. Differences in genome size were primarily due to independent post-speciation amplification of transposons. We also identified a set of Sitopsis genes pertinent to important agronomic traits that can be harnessed for wheat breeding. These newly assembled genome resources provide a new roadmap for evolutionary and genetic studies of the Triticum/Aegilops complex, as well as for wheat improvement.

The Plastome Sequences of Triticum sphaerococcum (ABD) and Triticum turgidum subsp. durum (AB) Exhibit Evolutionary Changes, Structural Characterization, Comparative Analysis, Phylogenomics and Time Divergence

The mechanism and course of Triticum plastome evolution is currently unknown; thus, it remains unclear how Triticum plastomes evolved during recent polyploidization. Here, we report the complete plastomes of two polyploid wheat species, Triticum sphaerococcum (AABBDD) and Triticum turgidum subsp. durum (AABB), and compare them with 19 available and complete Triticum plastomes to create the first map of genomic structural variation. Both T. sphaerococcum and T. turgidum subsp. durum plastomes were found to have a quadripartite structure, with plastome lengths of 134,531 bp and 134,015 bp, respectively. Furthermore, diploid (AA), tetraploid (AB, AG) and hexaploid (ABD, AGA^m) Triticum species plastomes displayed a conserved gene content and commonly harbored an identical set of annotated unique genes. Overall, there was a positive correlation between the number of repeats and plastome size. In all plastomes, the number of tandem repeats was higher than the number of palindromic and forward repeats. We constructed a Triticum phylogeny based on the complete plastomes and 42 shared genes from 71 plastomes. We estimated the divergence of Hordeum vulgare from wheat around 11.04-11.9 million years ago (mya) using a well-resolved plastome tree. Similarly, Sitopsis species diverged 2.8-2.9 mya before Triticum urartu (AA) and Triticum monococcum (AA). Aegilops speltoides was shown to be the maternal donor of polyploid wheat genomes and diverged ~0.2-0.9 mya. The phylogeny and divergence time estimates presented here can act as a reference framework for future studies of Triticum evolution.

Transposable Element Populations Shed Light on the Evolutionary History of Wheat and the Complex Co-Evolution of Autonomous and Non-Autonomous Retrotransposons

Wheat has one of the largest and most repetitive genomes among major crop plants, containing over 85% transposable elements (TEs). TEs populate genomes much in the way that individuals populate ecosystems, diversifying into different lineages, sub-families and sub-populations. The recent availability of high-quality, chromosome-scale genome sequences from ten wheat lines enables a detailed analysis how TEs evolved in allohexaploid wheat, its diploids progenitors, and in various chromosomal haplotype segments. LTR retrotransposon families evolved into distinct sub-populations and sub-families that were active in waves lasting several hundred thousand years. Furthermore, It is shown that different retrotransposon sub-families were active in the three wheat sub-genomes, making them useful markers to study and date polyploidization events and chromosomal rearrangements. Additionally, haplotype-specific TE sub-families are used to characterize chromosomal introgressions in different wheat lines. Additionally, populations of non-autonomous TEs co-evolved over millions of years with their autonomous partners, leading to complex systems with multiple types of autonomous, semi-autonomous and non-autonomous elements. Phylogenetic and TE population analyses revealed the relationships between non-autonomous elements and their mobilizing autonomous partners. TE population analysis provided insights into genome evolution of allohexaploid wheat and genetic diversity of species, and may have implication for future crop breeding.

Genetic diversity, distribution and domestication history of the neglected GGAtAt genepool of wheat

We present a comprehensive survey of cytogenetic and genomic diversity of the GGA^tA^t genepool of wheat, thereby unlocking these plant genetic resources for wheat improvement. Wheat yields are stagnating around the world and new sources of genes for resistance or tolerances to abiotic traits are required. In this context, the tetraploid wheat wild relatives are among the key candidates for wheat improvement. Despite its potential huge value for wheat breeding, the tetraploid GGA^tA^t genepool is largely neglected. Understanding the population structure, native distribution range, intraspecific variation of the entire tetraploid GGA^tA^t genepool and its domestication history would further its use for wheat improvement. The paper provides the first comprehensive survey of genomic and cytogenetic diversity sampling the full breadth and depth of the tetraploid GGA^tA^t genepool. According to the results obtained, the extant GGA^tA^t genepool consists of three distinct lineages. We provide detailed insights into the cytogenetic composition of GGA^tA^t wheats, revealed group- and population-specific markers and show that chromosomal rearrangements play an important role in intraspecific diversity of T. araraticum. The origin and domestication history of the GGA^tA^t lineages is discussed in the context of state-of-the-art archaeobotanical finds. We shed new light on the complex evolutionary history of the GGA^tA^t wheat genepool and provide the basis for an increased use of the GGA^tA^t wheat genepool for wheat improvement. The findings have implications for our understanding of the origins of agriculture in southwest Asia.

Introducing Beneficial Alleles from Plant Genetic Resources into the Wheat Germplasm

Wheat (Triticum sp.) is one of the world's most important crops, and constantly increasing its productivity is crucial to the livelihoods of millions of people. However, more than a century of intensive breeding and selection processes have eroded genetic diversity in the elite genepool, making new genetic gains difficult. Therefore, the need to introduce novel genetic diversity into modern wheat has become increasingly important. This review provides an overview of the plant genetic resources (PGR) available for wheat. We describe the most important taxonomic and phylogenetic relationships of these PGR to guide their use in wheat breeding. In addition, we present the status of the use of some of these resources in wheat breeding programs. We propose several introgression schemes that allow the transfer of qualitative and quantitative alleles from PGR into elite germplasm. With this in mind, we propose the use of a stage-gate approach to align the pre-breeding with main breeding programs to meet the needs of breeders, farmers, and end-users. Overall, this review provides a clear starting point to guide the introgression of useful alleles over the next decade.

Chloroplast Phylogenomic Analyses Resolve Multiple Origins of the Kengyilia Species (Poaceae: Triticeae) via Independent Polyploidization Events

Kengyilia is a group of allohexaploid species that arose from two hybridization events followed by genome doubling of three ancestral diploid species with different genomes St, Y, and P in the Triticeae. Estimating the phylogenetic relationship in resolution of the maternal lineages has been difficult, owing to the extremely low rate of sequence divergence. Here, phylogenetic reconstructions based on the plastome sequences were used to explore the role of maternal progenitors in the establishment of Kengyilia polyploid species. The plastome sequences of 11 Kengyilia species were analyzed together with 12 tetraploid species (PP, StP, and StY) and 33 diploid taxa representing 20 basic genomes in the Triticeae. Phylogenomic analysis and genetic divergence patterns suggested that (1) Kengyilia is closely related to Roegneria, Pseudoroegneria, Agropyron, Lophopyrum, Thinopyrum, and Dasypyrum; (2) both the StY genome Roegneria tetraploids and the PP genome Agropyron tetraploids served as the maternal donors during the speciation of Kengyilia species; (3) the different Kengyilia species derived their StY genome from different Roegneria species. Multiple origins of species via independent polyploidization events have occurred in the genus Kengyilia, resulting in a maternal haplotype polymorphism. This helps explain the rich diversity and wide adaptation of polyploid species in the genus Kengyilia.

Characterizing chloroplast genomes and inferring maternal divergence of the Triticum-Aegilops complex

The Triticum (wheat)-Aegilops (goatgrass) complex has been extensively studied, but the evolutionary history of polyploid wheats has not been fully elucidated. The chloroplast (cp) with maternal inheritance and homoplasy can simplify the sequence-based evolutionary inferences, but informative inferences would require a complete and accurate cp genome sequence. In this study, 16 cp genomes representing five Aegilops and 11 Triticum species and subspecies were sequenced, assembled and annotated, yielding five novel circular cp genome sequences. Analyzing the assembled cp genomes revealed no marked differences in genome structure and gene arrangement across the assayed species. A polymorphism analysis of 72 published cp genome sequences representing 10 Aegilops and 15 Triticum species and subspecies detected 1183 SNPs and 1881 SSRs. More than 80% SNPs detected resided on the downstream and upstream gene regions and only 2.78% or less SNPs were predicted to be deleterious. The largest nucleotide diversity was observed in the short single-copy genomic region. Relatively weak selection pressure on cp coding genes was detected. Different phylogenetic analyses confirmed that the maternal divergence of the Triticum-Aegilops complex had three deep lineages each representing a diploid species with nuclear A, B, or D genome. Dating the maternal divergence yielded age estimates of divergence that matched well with those reported previously. The divergence between emmer and bread wheats occurred at 8200-11,200 years ago. These findings are useful for further genomic studies, provide insight into cp genome evolvability and allow for better understanding of the maternal divergence of the Triticum-Aegilops complex.

Chromosome and Molecular Analyses Reveal Significant Karyotype Diversity and Provide New Evidence on the Origin of Aegilops columnaris

Aegilops columnaris Zhuk. is tetraploid grass species (2n = 4x = 28, UcUcXcXc) closely related to Ae. neglecta and growing in Western Asia and a western part of the Fertile Crescent. Genetic diversity of Ae. columnaris was assessed using C-banding, FISH, nuclear and chloroplast (cp) DNA analyses, and gliadin electrophoresis. Cytogenetically Ae. columnaris was subdivided into two groups, C-I and C-II, showing different karyotype structure, C-banding, and FISH patterns. C-I group was more similar to Ae. neglecta. All types of markers revealed significant heterogeneity in C-II group, although group C-I was also polymorphic. Two chromosomal groups were consistent with plastogroups identified in a current study based on sequencing of three chloroplast intergenic spacer regions. The similarity of group C-I of Ae. columnaris with Ae. neglecta and their distinctness from C-II indicate that divergence of the C-I group was associated with minor genome modifications. Group C-II could emerge from C-I relatively recently, probably due to introgression from another Aegilops species followed by a reorganization of the parental genomes. Most C-II accessions were collected from a very narrow geographic region, and they might originate from a common ancestor. We suggest that the C-II group is at the initial stage of species divergence and undergoing an extensive speciation process.

Variation in Plastome Sizes Accompanied by Evolutionary History in Monogenomic Triticeae (Poaceae: Triticeae)

To investigate the pattern of chloroplast genome variation in Triticeae, we comprehensively analyzed the indels in protein-coding genes and intergenic sequence, gene loss/pseudonization, intron variation, expansion/contraction in inverted repeat regions, and the relationship between sequence characteristics and chloroplast genome size in 34 monogenomic Triticeae plants. Ancestral genome reconstruction suggests that major length variations occurred in four-stem branches of monogenomic Triticeae followed by independent changes in each genus. It was shown that the chloroplast genome sizes of monogenomic Triticeae were highly variable. The chloroplast genome of Pseudoroegneria, Dasypyrum, Lophopyrum, Thinopyrum, Eremopyrum, Agropyron, Australopyrum, and Henradia in Triticeae had evolved toward size reduction largely because of pseudogenes elimination events and length deletion fragments in intergenic. The Aegilops/Triticum complex, Taeniatherum, Secale, Crithopsis, Herteranthelium, and Hordeum in Triticeae had a larger chloroplast genome size. The large size variation in major lineages and their subclades are most likely consequences of adaptive processes since these variations were significantly correlated with divergence time and historical climatic changes. We also found that several intergenic regions, such as petN-trnC and psbE-petL containing unique genetic information, which can be used as important tools to identify the maternal relationship among Triticeae species. Our results contribute to the novel knowledge of plastid genome evolution in Triticeae.

Comparison of the cytoplastic genomes by resequencing: insights into the genetic diversity and the phylogeny of the agriculturally important genus Brassica

Background

The genus Brassica mainly comprises three diploid and three recently derived allotetraploid species, most of which are highly important vegetable, oil or ornamental crops cultivated worldwide. Despite being extensively studied, the origination of B. napus and certain detailed interspecific relationships within Brassica genus remains undetermined and somewhere confused. In the current high-throughput sequencing era, a systemic comparative genomic study based on a large population is necessary and would be crucial to resolve these questions.

Results

The chloroplast DNA and mitochondrial DNA were synchronously resequenced in a selected set of Brassica materials, which contain 72 accessions and maximally integrated the known Brassica species. The Brassica genomewide cpDNA and mtDNA variations have been identified. Detailed phylogenetic relationships inside and around Brassica genus have been delineated by the cpDNA- and mtDNA- variation derived phylogenies. Different from B. juncea and B. carinata, the natural B. napus contains three major cytoplasmic haplotypes: the cam-type which directly inherited from B. rapa, polima-type which is close to cam-type as a sister, and the mysterious but predominant nap-type. Certain sparse C-genome wild species might have primarily contributed the nap-type cytoplasm and the corresponding C subgenome to B. napus, implied by their con-clustering in both phylogenies. The strictly concurrent inheritance of mtDNA and cpDNA were dramatically disturbed in the B. napus cytoplasmic male sterile lines (e.g., mori and nsa). The genera Raphanus, Sinapis, Eruca, Moricandia show a strong parallel evolutional relationships with Brassica.

Conclusions

The overall variation data and elaborated phylogenetic relationships provide further insights into genetic understanding of Brassica, which can substantially facilitate the development of novel Brassica germplasms.

Homoeologous exchanges occur through intragenic recombination generating novel transcripts and proteins in wheat and other polyploids

Recombination between homeologous chromosomes, also known as homeologous exchange (HE), plays a significant role in shaping genome structure and gene expression in interspecific hybrids and allopolyploids of several plant species. However, the molecular mechanisms that govern HEs are not well understood. Here, we studied HE events in the progeny of a nascent allotetraploid (genome AADD) derived from two diploid progenitors of hexaploid bread wheat using cytological and whole-genome sequence analyses. In total, 37 HEs were identified and HE junctions were mapped precisely. HEs exhibit typical patterns of homologous recombination hotspots, being biased toward low-copy, subtelomeric regions of chromosome arms and showing association with known recombination hotspot motifs. But, strikingly, while homologous recombination preferentially takes place upstream and downstream of coding regions, HEs are highly enriched within gene bodies, giving rise to novel recombinant transcripts, which in turn are predicted to generate new protein fusion variants. To test whether this is a widespread phenomenon, a dataset of high-resolution HE junctions was analyzed for allopolyploid Brassica, rice, Arabidopsis suecica, banana, and peanut. Intragenic recombination and formation of chimeric genes was detected in HEs of all species and was prominent in most of them. HE thus provides a mechanism for evolutionary novelty in transcript and protein sequences in nascent allopolyploids.

Evolutionary patterns of plastome uncover diploid-polyploid maternal relationships in Triticeae

To investigate the diploid-polyploid relationships and the role of maternal progenitors in establishment of polyploid richness in Triticeae, 35 polyploids representing almost all genomic constitutions together with 48 diploid taxa representing 20 basic genomes in the tribe were analyzed. Phylogenomic reconstruction, genetic distance matrix, and nucleotide diversity patterns of plastome sequences indicated that (1) The maternal donor of the annual polyploid species with the U- and D-genome are related to extant Ae. umbellulata and Ae. tauschii, respectively. The maternal donor to the annual polyploid species with the S-, G-, and B-genome originated from the species of Sitopsis section of the genus Aegilops. The annual species with the Xe-containing polyploids were donated by Eremopyrum as the female parent; (2) Pseudoroegneria and Psathyrostachys were the maternal donor of perennial species with the St- and Ns-containing polyploids, respectively; (3) The Lophopyrum, Thinopyrum and Dasypyrum genomes contributed cytoplasm genome to Pseudoroegneria species as a result of incomplete lineage sorting and/or chloroplast captures, and these lineages were genetically transmitted to the St-containing polyploid species via polyploidization; (4) There is a reticulate relationship among the St-containing polyploid species. It can be suggested that genetic heterogeneity might associate with the richness of the polyploids in Triticeae.

The study of organelle DNA variability in alloplasmic barley lines in the NGS era

Alloplasmic lines are a suitable model for studying molecular coevolution and interrelations between genetic systems of plant cells. Whole chloroplast (cp) and mitochondrial (mt) genome sequences were obtained by the MiSeq System (Illumina). Organelle DNA samples were prepared from a set of 12 alloplasmic barley lines with different cytoplasms of Hordeum vulgare ssp. spontaneum and H. vulgare ssp. vulgare, as well as from their paternal varieties. A bioinformatic approach for analysis of NGS data obtained on an organellar DNA mix has been developed and verified. A comparative study of Hordeum organelle genomes’ variability and disposition of polymorphic loci was conducted. Eight types of chloroplast DNA and 5 types of mitochondrial DNA were distinguished for the barley sample set examined. These results were compared with the previous data of a restriction fragment length polymorphism (RFLP) study of organelle DNAs for the same material. Formerly established data about a field evaluation of alloplasmic barley lines were revised in the light of information about organelle genomes gained after NGS. Totally 17 polymorphic loci were found at exons of chloroplast genomes. Seven of the SNPs were located in the genes of the Ndh complex. The nonsynonymous changes of nucleotides were detected in the matK, rpoC1, ndhK, ndhG and infA genes. Some of the SNPs detected are very similar in codon position and in the type of amino acid substitution to the places where RNA editing can occur. Thus, these results outline new perspectives for the future study of nuclear-cytoplasmic interactions in alloplasmic lines.

The complete chloroplast genomes of seventeen Aegilops tauschii: genome comparative analysis and phylogenetic inference

The D genome progenitor of bread wheat, Aegilops tauschii Cosson (DD, 2n = 2x = 14), which is naturally distributed in Central Eurasia, ranging from northern Syria and Turkey to western China, is considered a potential genetic resource for improving bread wheat. In this study, the chloroplast (cp) genomes of 17 Ae. tauschii accessions were reconstructed. The cp genome sizes ranged from 135,551 bp to 136,009 bp and contained a typical quadripartite structure of angiosperms. Within these genomes, we identified a total of 124 functional genes, including 82 protein-coding genes, 34 transfer RNA genes and eight ribosomal RNA genes, with 17 duplicated genes in the IRs. Although the comparative analysis revealed that the genomic structure (gene order, gene number and IR/SC boundary regions) is conserved, a few variant loci were detected, predominantly in the non-coding regions (intergenic spacer regions). The phylogenetic relationships determined based on the complete genome sequences were consistent with the hypothesis that Ae. tauschii populations in the Yellow River region of China originated in South Asia not Xinjiang province or Iran, which could contribute to more effective utilization of wild germplasm resources. Furthermore, we confirmed that Ae. tauschii was derived from monophyletic speciation rather than hybrid speciation at the cp genome level. We also identified four variable genomic regions, rpl32-trnL-UAG, ccsA-ndhD, rbcL-psaI and rps18-rpl20, showing high levels of nucleotide polymorphisms, which may accordingly prove useful as cpDNA markers in studying the intraspecific genetic structure and diversity of Ae. tauschii.

Multiregional origins of the domesticated tetraploid wheats

We used genotyping-by-sequencing (GBS) to investigate the evolutionary history of domesticated tetraploid wheats. With a panel of 189 wild and domesticated wheats, we identified 1,172,469 single nucleotide polymorphisms (SNPs) with a read depth ≥3. Principal component analyses (PCAs) separated the Triticum turgidum and Triticum timopheevii accessions, as well as wild T. turgidum from the domesticated emmers and the naked wheats, showing that SNP typing by GBS is capable of providing robust information on the genetic relationships between wheat species and subspecies. The PCAs and a neighbour-joining analysis suggested that domesticated tetraploid wheats have closest affinity with wild emmers from the northern Fertile Crescent, consistent with the results of previous genetic studies on the origins of domesticated wheat. However, a more detailed examination of admixture and allele sharing between domesticates and different wild populations, along with genome-wide association studies (GWAS), showed that the domesticated tetraploid wheats have also received a substantial genetic input from wild emmers from the southern Levant. Taking account of archaeological evidence that tetraploid wheats were first cultivated in the southern Levant, we suggest that a pre-domesticated crop spread from this region to southeast Turkey and became mixed with a wild emmer population from the northern Fertile Crescent. Fixation of the domestication traits in this mixed population would account for the allele sharing and GWAS results that we report. We also propose that feralization of the component of the pre-domesticated population that did not acquire domestication traits has resulted in the modern wild population from southeast Turkey displaying features of both the domesticates and wild emmer from the southern Levant, and hence appearing to be the sole progenitor of domesticated tetraploids when the phylogenetic relationships are studied by methods that assume a treelike pattern of evolution.

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.

Cytonuclear Coevolution following Homoploid Hybrid Speciation in Aegilops tauschii

The diploid D-genome lineage of the Triticum/Aegilops complex has an evolutionary history involving genomic contributions from ancient A- and B/S-genome species. We explored here the possible cytonuclear evolutionary responses to this history of hybridization. Phylogenetic analysis of chloroplast DNAs indicates that the D-genome lineage has a maternal origin of the A-genome or some other closely allied lineage. Analyses of the nuclear genome in the D-genome species Aegilops tauschii indicate that accompanying and/or following this ancient hybridization, there has been biased maintenance of maternal A-genome ancestry in nuclear genes encoding cytonuclear enzyme complexes (CECs). Our study provides insights into mechanisms of cytonuclear coevolution accompanying the evolution and eventual stabilization of homoploid hybrid species. We suggest that this coevolutionary process includes likely rapid fixation of A-genome CEC orthologs as well as biased retention of A-genome nucleotides in CEC homologs following population level recombination during the initial generations.

Pervasive hybridizations in the history of wheat relatives

Cultivated wheats are derived from an intricate history of three genomes, A, B, and D, present in both diploid and polyploid species. It was recently proposed that the D genome originated from an ancient hybridization between the A and B lineages. However, this result has been questioned, and a robust phylogeny of wheat relatives is still lacking. Using transcriptome data from all diploid species and a new methodological approach, our comprehensive phylogenomic analysis revealed that more than half of the species descend from an ancient hybridization event but with a more complex scenario involving a different parent than previously thought-Aegilops mutica, an overlooked wild species-instead of the B genome. We also detected other extensive gene flow events that could explain long-standing controversies in the classification of wheat relatives.

A discriminatory test for the wheat B and G genomes reveals misclassified accessions of Triticum timopheevii and Triticum turgidum

The tetraploid wheat species Triticum turgidum and Triticum timopheevii are morphologically similar, and misidentification of material collected from the wild is possible. We compared published sequences for the Ppd-A1, Ppd-B1 and Ppd-G1 genes from multiple accessions of T. turgidum and T. timopheevii and devised a set of four polymerase chain reactions (PCRs), two specific for Ppd-B1 and two for Ppd-G1. We used these PCRs with 51 accessions of T. timopheevii and 20 of T. turgidum. Sixty of these accessions gave PCR products consistent with their taxon identifications, but the other eleven accessions gave anomalous results: ten accessions that were classified as T. turgidum were identified as T. timopheevii by the PCRs, and one T. timopheevii accession was typed as T. turgidum. We believe that these anomalies are not due to errors in the PCR tests because the results agree with a more comprehensive analysis of genome-wide single nucleotide polymorphisms, which similarly suggest that these eleven accessions have been misclassified. Our results therefore show that the accepted morphological tests for discrimination between T. turgidum and T. timopheevii might not be entirely robust, but that species identification can be made cheaply and quickly by PCRs directed at the Ppd-1 gene.

Divergence between bread wheat and Triticum militinae in the powdery mildew resistance QPm.tut-4A locus and its implications for cloning of the resistance gene

A segment of Triticum militinae chromosome 7G harbors a gene(s) conferring powdery mildew resistance which is effective at both the seedling and the adult plant stages when transferred into bread wheat (T. aestivum). The introgressed segment replaces a piece of wheat chromosome arm 4AL. An analysis of segregating materials generated to positionally clone the gene highlighted that in a plant heterozygous for the introgression segment, only limited recombination occurs between the introgressed region and bread wheat 4A. Nevertheless, 75 genetic markers were successfully placed within the region, thereby confining the gene to a 0.012 cM window along the 4AL arm. In a background lacking the Ph1 locus, the localized rate of recombination was raised 33-fold, enabling the reduction in the length of the region containing the resistance gene to a 480 kbp stretch harboring 12 predicted genes. The substituted segment in the reference sequence of bread wheat cv. Chinese Spring is longer (640 kbp) and harbors 16 genes. A comparison of the segments' sequences revealed a high degree of divergence with respect to both their gene content and nucleotide sequence. Of the 12 T. militinae genes, only four have a homolog in cv. Chinese Spring. Possible candidate genes for the resistance have been identified based on function predicted from their sequence.

Experimental evolutionary studies on the genetic autonomy of the cytoplasmic genome "plasmon" in the Triticum (wheat)-Aegilops complex

The term "plasmon" is used to indicate the whole cytoplasmic genetic system, whereas "genome" refers to the whole nuclear genetic system. Although maternal inheritance of the plasmon is well documented in angiosperms, its genetic autonomy from the coexisting nuclear genome still awaits critical examination. We tested this autonomy in two related studies: One was to determine the persistence of the genetic effect of the plasmon of Aegilops caudata (genome CC) on the phenotype of common wheat, Triticum aestivum strain "Tve" (genome AABBDD), during 63 y (one generation per year) of repeated backcrosses of Ae. caudata and its offspring with pollen of the same Tve wheat, and the second was to reconstruct an Ae. caudata strain from the genome of this strain and its plasmon that had been resident in Tve wheat for 50 generations, and to compare the phenotypic and organellar DNA characteristics between the native and reconstructed strains. Results indicated no change in the effect of Ae. caudata plasmon on Tve wheat during its stay in wheat for more than half a century, and no difference between the native and reconstructed caudata strains in their phenotype and simple sequence repeats in their organellar DNAs, thus demonstrating the prolonged genetic autonomy of the plasmon from the coexisting genomes of wheat and several other species that were used in the reconstruction of Ae. caudata The relationship between the proven genetic autonomy of the plasmon under changing nuclear conditions and its diversification during evolution of the Triticum-Aegilops complex is discussed.

DArTseq-based analysis of genomic relationships among species of tribe Triticeae

Precise utilization of wild genetic resources to improve the resistance of their cultivated relatives to environmental growth limiting factors, such as salinity stress and diseases, requires a clear understanding of their genomic relationships. Although seriously criticized, analyzing these relationships in tribe Triticeae has largely been based on meiotic chromosome pairing in hybrids of wide crosses, a specialized and labourious strategy. In this study, DArTseq, an efficient genotyping-by-sequencing platform, was applied to analyze the genomes of 34 Triticeae species. We reconstructed the phylogenetic relationships among diploid and polyploid Aegilops and Triticum species, including hexaploid wheat. Tentatively, we have identified the diploid genomes that are likely to have been involved in the evolution of five polyploid species of Aegilops, which have remained unresolved for decades. Explanations which cast light on the progenitor of the A genomes and the complex genomic status of the B/G genomes of polyploid Triticum species in the Emmer and Timopheevi lineages of wheat have also been provided. This study has, therefore, demonstrated that DArTseq genotyping can be effectively applied to analyze the genomes of plants, especially where their genome sequence information are not available.

Structural variation and rates of genome evolution in the grass family seen through comparison of sequences of genomes greatly differing in size

Homology was searched with genes annotated in the Aegilops tauschii pseudomolecules against genes annotated in the pseudomolecules of tetraploid wild emmer wheat, Brachypodium distachyon, sorghum and rice. Similar searches were performed with genes annotated in the rice pseudomolecules. Matrices of collinear genes and rearrangements in their order were constructed. Optical BioNano genome maps were constructed and used to validate rearrangements unique to the wild emmer and Ae. tauschii genomes. Most common rearrangements were short paracentric inversions and short intrachromosomal translocations. Intrachromosomal translocations outnumbered segmental intrachromosomal duplications. The densities of paracentric inversion lengths were approximated by exponential distributions in all six genomes. Densities of collinear genes along the Ae. tauschii chromosomes were highly correlated with meiotic recombination rates but those of rearrangements were not, suggesting different causes of the erosion of gene collinearity and evolution of major chromosome rearrangements. Frequent rearrangements sharing breakpoints suggested that chromosomes have been rearranged recurrently at some sites. The distal 4 Mb of the short arms of rice chromosomes Os11 and Os12 and corresponding regions in the sorghum, B. distachyon and Triticeae genomes contain clusters of interstitial translocations including from 1 to 7 collinear genes. The rates of acquisition of major rearrangements were greater in the large wild emmer wheat and Ae. tauschii genomes than in the lineage preceding their divergence or in the B. distachyon, rice and sorghum lineages. It is suggested that synergy between large quantities of dynamic transposable elements and annual growth habit have been the primary causes of the fast evolution of the Triticeae genomes.

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.

Asymmetrical changes of gene expression, small RNAs and chromatin in two resynthesized wheat allotetraploids

Polyploidy occurs in some animals and all flowering plants, including important crops such as wheat. The consequences of polyploidy in crops remain elusive, partly because their progenitors are unknown. Using two resynthesized wheat allotetraploids S^l S^l AA and AADD with known diploid progenitors, we analyzed mRNA and small RNA transcriptomes in the endosperm, compared transcriptomes between endosperm and root in AADD, and examined chromatin changes in the allotetraploids. In the endosperm, there were more non-additively expressed genes in S^l S^l AA than in AADD. In AADD, non-additively expressed genes were developmentally regulated, and the majority (62-70%) were repressed. The repressed genes in AADD included a group of histone methyltransferase gene homologs, which correlated with reduced histone H3K9me2 levels and activation of various transposable elements in AADD. In S^l S^l AA, there was a tendency for expression dominance of S^l over A homoeologs, but the histone methyltransferase gene homologs were additively expressed, correlating with insignificant changes in histone H3K9me2 levels. Moreover, more 24-nucleotide small inferring RNAs (siRNAs) in the A subgenome were disrupted in AADD than in S^l S^l AA, which were associated with expression changes of siRNA-associated genes. Our results indicate that asymmetrical changes in siRNAs, chromatin modifications, transposons and gene expression coincide with unstable AADD genomes and stable S^l S^l AA genomes, which could help explain the evolutionary trajectories of wheat allotetraploids formed by different progenitors.

A 250 plastome phylogeny of the grass family (Poaceae): topological support under different data partitions

The systematics of grasses has advanced through applications of plastome phylogenomics, although studies have been largely limited to subfamilies or other subgroups of Poaceae. Here we present a plastome phylogenomic analysis of 250 complete plastomes (179 genera) sampled from 44 of the 52 tribes of Poaceae. Plastome sequences were determined from high throughput sequencing libraries and the assemblies represent over 28.7 Mbases of sequence data. Phylogenetic signal was characterized in 14 partitions, including (1) complete plastomes; (2) protein coding regions; (3) noncoding regions; and (4) three loci commonly used in single and multi-gene studies of grasses. Each of the four main partitions was further refined, alternatively including or excluding positively selected codons and also the gaps introduced by the alignment. All 76 protein coding plastome loci were found to be predominantly under purifying selection, but specific codons were found to be under positive selection in 65 loci. The loci that have been widely used in multi-gene phylogenetic studies had among the highest proportions of positively selected codons, suggesting caution in the interpretation of these earlier results. Plastome phylogenomic analyses confirmed the backbone topology for Poaceae with maximum bootstrap support (BP). Among the 14 analyses, 82 clades out of 309 resolved were maximally supported in all trees. Analyses of newly sequenced plastomes were in agreement with current classifications. Five of seven partitions in which alignment gaps were removed retrieved Panicoideae as sister to the remaining PACMAD subfamilies. Alternative topologies were recovered in trees from partitions that included alignment gaps. This suggests that ambiguities in aligning these uncertain regions might introduce a false signal. Resolution of these and other critical branch points in the phylogeny of Poaceae will help to better understand the selective forces that drove the radiation of the BOP and PACMAD clades comprising more than 99.9% of grass diversity.

Complete structure and variation of the chloroplast genome of Agropyron cristatum (L.) Gaertn

Agropyron cristatum (L.) Gaertner, a perennial grass in the tribe Triticeae (Poaceae), is a wild relative of cereal crops that is suitable for genetic improvement. In this study, we first sequenced the complete chloroplast (cp) genome of Ag. cristatum using Hiseq4000 PE150. The Ag. cristatum chloroplast genome is 135,554bp in length, has a typical quadripartite structure and contains 76 protein-coding genes, 29 tRNA genes and four rRNA genes. The cp genome of Ag. cristatum was used for comparison with other seven Triticeae species. One large variable region (800bp), which primarily contained the rpl23 (non-reciprocally translocated from IRs) and accD genes, was detected between rbcL gene and psaI gene within LSC region. The deletion of the accD and translocated rpl23 genes in Ag. cristatum indicated an independent gene-loss events or an additional divergence in Triticeae. Analyses of the dn/ds ratio and K2-P's genetic distance for 76 protein-coding genes showed that genes with evolutionary divergence might suffer from the effect of sequence regional constraints or gene functional constraints in Triticeae species. Our research will generally contribute to the knowledge of plastid genome evolution in Triticeae.

Evolution of the S-Genomes in Triticum-Aegilops Alliance: Evidences From Chromosome Analysis

Five diploid Aegilops species of the Sitopsis section: Ae. speltoides, Ae. longissima, Ae. sharonensis, Ae. searsii, and Ae. bicornis, two tetraploid species Ae. peregrina (= Ae. variabilis) and Ae. kotschyi (Aegilops section) and hexaploid Ae. vavilovii (Vertebrata section) carry the S-genomes. The B- and G-genomes of polyploid wheat are also the derivatives of the S-genome. Evolution of the S-genome species was studied using Giemsa C-banding and fluorescence in situ hybridization (FISH) with DNA probes representing 5S (pTa794) and 18S-5.8S-26S (pTa71) rDNAs as well as nine tandem repeats: pSc119.2, pAesp_SAT86, Spelt-1, Spelt-52, pAs1, pTa-535, and pTa-s53. To correlate the C-banding and FISH patterns we used the microsatellites (CTT)₁₀ and (GTT)₉, which are major components of the C-banding positive heterochromatin in wheat. According to the results obtained, diploid species split into two groups corresponding to Emarginata and Truncata sub-sections, which differ in the C-banding patterns, distribution of rDNA and other repeats. The B- and G-genomes of polyploid wheat are most closely related to the S-genome of Ae. speltoides. The genomes of allopolyploid wheat have been evolved as a result of different species-specific chromosome translocations, sequence amplification, elimination and re-patterning of repetitive DNA sequences. These events occurred independently in different wheat species and in Ae. speltoides _. The 5S rDNA locus of chromosome 1S was probably lost in ancient Ae. speltoides prior to formation of Timopheevii wheat, but after the emergence of ancient emmer. Evolution of Emarginata species was associated with an increase of C-banding and (CTT)₁₀-positive heterochromatin, amplification of Spelt-52, re-pattering of the pAesp_SAT86, and a gradual decrease in the amount of the D-genome-specific repeats pAs1, pTa-535, and pTa-s53. The emergence of Ae. peregrina and Ae. kotschyi did not lead to significant changes of the S^*-genomes. However, partial elimination of 45S rDNA repeats from 5S^* and 6S^* chromosomes and alterations of C-banding and FISH-patterns have been detected. Similarity of the S^v-genome of Ae. vavilovii with the S^s genome of diploid Ae. searsii confirmed the origin of this hexaploid. A model of the S-genome evolution is suggested.

A High-Density Genetic Map of Wild Emmer Wheat from the Karaca Dağ Region Provides New Evidence on the Structure and Evolution of Wheat Chromosomes

Wild emmer (Triticum turgidum ssp. dicoccoides) is a progenitor of all cultivated wheat grown today. It has been hypothesized that emmer was domesticated in the Karaca Dağ region in southeastern Turkey. A total of 445 recombinant inbred lines of T. turgidum ssp. durum cv. 'Langdon' x wild emmer accession PI 428082 from this region was developed and genotyped with the Illumina 90K single nucleotide polymorphism Infinium assay. A genetic map comprising 2,650 segregating markers was constructed. The order of the segregating markers and an additional 8,264 co-segregating markers in the Aegilops tauschii reference genome sequence was used to compare synteny of the tetraploid wheat with the Brachypodium distachyon, rice, and sorghum. These comparisons revealed the presence of 15 structural chromosome rearrangements, in addition to the already known 4A-5A-7B rearrangements. The most common type was an intra-chromosomal translocation in which the translocated segment was short and was translocated only a short distance along the chromosome. A large reciprocal translocation, one small non-reciprocal translocation, and three large and one small paracentric inversions were also discovered. The use of inversions for a phylogeny reconstruction in the Triticum-Aegilops alliance was illustrated. The genetic map was inconsistent with the current model of evolution of the rearranged chromosomes 4A-5A-7B. Genetic diversity in the rearranged chromosome 4A showed that the rearrangements might have been contemporary with wild emmer speciation. A selective sweep was found in the centromeric region of chromosome 4A in Karaca Dağ wild emmer but not in 4A of T. aestivum. The absence of diversity from a large portion of chromosome 4A of wild emmer, believed to be ancestral to all domesticated wheat, is puzzling.

Revisiting Pivotal-Differential Genome Evolution in Wheat

An interesting and possibly unique pattern of genome evolution following polyploidy can be observed among allopolyploids of the Triticum and Aegilops genera (wheat group). Most polyploids in this group are presumed to share a common unaltered (pivotal) subgenome (U, D, or A) together with one or two modified (differential) subgenomes, a status that has been referred to as 'pivotal-differential' genome evolution. In this review we discuss various mechanisms that could be responsible for this evolutionary pattern, as well as evidence for and against the putative evolutionary mechanisms involved. We suggest that, in light of recent advances in genome sequencing and related technologies in the wheat group, the time has come to reopen the investigation into pivotal-differential genome evolution.

Dated tribe-wide whole chloroplast genome phylogeny indicates recurrent hybridizations within Triticeae

BACKGROUND

Triticeae, the tribe of wheat grasses, harbours the cereals barley, rye and wheat and their wild relatives. Although economically important, relationships within the tribe are still not understood. We analysed the phylogeny of chloroplast lineages among nearly all monogenomic Triticeae taxa and polyploid wheat species aiming at a deeper understanding of the tribe's evolution. We used on- and off-target reads of a target-enrichment experiment followed by Illumina sequencing.

RESULTS

The read data was used to assemble the plastid locus ndhF for 194 individuals and the whole chloroplast genome for 183 individuals, representing 53 Triticeae species and 15 genera. We conducted Bayesian and multispecies coalescent analyses to infer relationships and estimate divergence times of the taxa. We present the most comprehensive dated Triticeae chloroplast phylogeny and review previous hypotheses in the framework of our results. Monophyly of Triticeae chloroplasts could not be confirmed, as either Bromus or Psathyrostachys captured a chloroplast from a lineage closely related to a Bromus-Triticeae ancestor. The most recent common ancestor of Triticeae occurred approximately between ten and 19 million years ago.

CONCLUSIONS

The comparison of the chloroplast phylogeny with available nuclear data in several cases revealed incongruences indicating past hybridizations. Recent events of chloroplast capture were detected as individuals grouped apart from con-specific accessions in otherwise monopyhletic groups.

Genomic history of the origin and domestication of common bean unveils its closest sister species

Background

Modern civilization depends on only a few plant species for its nourishment. These crops were derived via several thousands of years of human selection that transformed wild ancestors into high-yielding domesticated descendants. Among cultivated plants, common bean (Phaseolus vulgaris L.) is the most important grain legume. Yet, our understanding of the origins and concurrent shaping of the genome of this crop plant is limited.

Results

We sequenced the genomes of 29 accessions representing 12 Phaseolus species. Single nucleotide polymorphism-based phylogenomic analyses, using both the nuclear and chloroplast genomes, allowed us to detect a speciation event, a finding further supported by metabolite profiling. In addition, we identified ~1200 protein coding genes (PCGs) and ~100 long non-coding RNAs with domestication-associated haplotypes. Finally, we describe asymmetric introgression events occurring among common bean subpopulations in Mesoamerica and across hemispheres.

Conclusions

We uncover an unpredicted speciation event in the tropical Andes that gave rise to a sibling species, formerly considered the “wild ancestor” of P. vulgaris, which diverged before the split of the Mesoamerican and Andean P. vulgaris gene pools. Further, we identify haplotypes strongly associated with genes underlying the emergence of domestication traits. Our findings also reveal the capacity of a predominantly autogamous plant to outcross and fix loci from different populations, even from distant species, which led to the acquisition by domesticated beans of adaptive traits from wild relatives. The occurrence of such adaptive introgressions should be exploited to accelerate breeding programs in the near future.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-017-1190-6) contains supplementary material, which is available to authorized users.

Identification and Analysis of RNA Editing Sites in the Chloroplast Transcripts of Aegilops tauschii L

RNA editing is an important way to convert cytidine (C) to uridine (U) at specific sites within RNA molecules at a post-transcriptional level in the chloroplasts of higher plants. Although it has been systematically studied in many plants, little is known about RNA editing in the wheat D genome donor Aegilops tauschii L. Here, we investigated the chloroplast RNA editing of Ae. tauschii and compared it with other wheat relatives to trace the evolution of wheat. Through bioinformatics prediction, a total of 34 C-to-U editing sites were identified, 17 of which were validated using RT-PCR product sequencing. Furthermore, 60 sites were found by the RNA-Seq read mapping approach, 24 of which agreed with the prediction and six were validated experimentally. The editing sites were biased toward tCn or nCa trinucleotides and 5′-pyrimidines, which were consistent with the flanking bases of editing sites of other seed plants. Furthermore, the editing events could result in the alteration of the secondary structures and topologies of the corresponding proteins, suggesting that RNA editing might impact the function of target genes. Finally, comparative analysis found some evolutionarily conserved editing sites in wheat and two species-specific sites were also obtained. This study is the first to report on RNA editing in Aegilops tauschii L, which not only sheds light on the evolution of wheat from the point of view of RNA editing, but also lays a foundation for further studies to identify the mechanisms of C-to-U alterations.

Recurrence of Chromosome Rearrangements and Reuse of DNA Breakpoints in the Evolution of the Triticeae Genomes

Chromosomal rearrangements (CRs) play important roles in karyotype diversity and speciation. While many CR breakpoints have been characterized at the sequence level in yeast, insects, and primates, little is known about the structure of evolutionary CR breakpoints in plant genomes, which are much more dynamic in genome size and sequence organization. Here, we report identification of breakpoints of a translocation between chromosome arms 4L and 5L of Triticeae, which is fixed in several species, including diploid wheat and rye, by comparative mapping and analysis of the draft genome and chromosome survey sequences of the Triticeae species. The wheat translocation joined the ends of breakpoints downstream of a WD40 gene on 4AL and a gene of the PMEI family on 5AL. A basic helix-loop-helix transcription factor gene in 5AL junction was significantly restructured. Rye and wheat share the same position for the 4L breakpoint, but the 5L breakpoint positions are not identical, although very close in these two species, indicating the recurrence of 4L/5L translocations in the Triticeae. Although barley does not carry the translocation, collinearity across the breakpoints was violated by putative inversions and/or transpositions. Alignment with model grass genomes indicated that the translocation breakpoints coincided with ancient inversion junctions in the Triticeae ancestor. Our results show that the 4L/5L translocation breakpoints represent two CR hotspots reused during Triticeae evolution, and support breakpoint reuse as a widespread mechanism in all eukaryotes. The mechanisms of the recurrent translocation and its role in Triticeae evolution are also discussed.

Genome-wide characterization of microsatellites in Triticeae species: abundance, distribution and evolution

Microsatellites are an important constituent of plant genome and distributed across entire genome. In this study, genome-wide analysis of microsatellites in 8 Triticeae species and 9 model plants revealed that microsatellite characteristics were similar among the Triticeae species. Furthermore, genome-wide microsatellite markers were designed in wheat and then used to analyze the evolutionary relationship of wheat and other Triticeae species. Results displayed that Aegilops tauschii was found to be the closest species to Triticum aestivum, followed by Triticum urartu, Triticum turgidum and Aegilops speltoides, while Triticum monococcum, Aegilops sharonensis and Hordeum vulgare showed a relatively lower PCR amplification effectivity. Additionally, a significantly higher PCR amplification effectivity was found in chromosomes at the same subgenome than its homoeologous when these markers were subjected to search against different chromosomes in wheat. After a rigorous screening process, a total of 20,666 markers showed high amplification and polymorphic potential in wheat and its relatives, which were integrated with the public available wheat markers and then anchored to the genome of wheat (CS). This study not only provided the useful resource for SSR markers development in Triticeae species, but also shed light on the evolution of polyploid wheat from the perspective of microsatellites.

Molecular evolution of Wcor15 gene enhanced our understanding of the origin of A, B and D genomes in Triticum aestivum

The allohexaploid bread wheat originally derived from three closely related species with A, B and D genome. Although numerous studies were performed to elucidate its origin and phylogeny, no consensus conclusion has reached. In this study, we cloned and sequenced the genes Wcor15-2A, Wcor15-2B and Wcor15-2D in 23 diploid, 10 tetraploid and 106 hexaploid wheat varieties and analyzed their molecular evolution to reveal the origin of the A, B and D genome in Triticum aestivum. Comparative analyses of sequences in diploid, tetraploid and hexaploid wheats suggest that T. urartu, Ae. speltoides and Ae. tauschii subsp. strangulata are most likely the donors of the Wcor15-2A, Wcor15-2B and Wcor15-2D locus in common wheat, respectively. The Wcor15 genes from subgenomes A and D were very conservative without insertion and deletion of bases during evolution of diploid, tetraploid and hexaploid. Non-coding region of Wcor15-2B gene from B genome might mutate during the first polyploidization from Ae. speltoides to tetraploid wheat, however, no change has occurred for this gene during the second allopolyploidization from tetraploid to hexaploid. Comparison of the Wcor15 gene shed light on understanding of the origin of the A, B and D genome of common wheat.