Variation in repeated nucleotide sequences sheds light on the phylogeny of the wheat B and G genomes

Chromosome-scale genome assembly of bread wheat's wild relative Triticum timopheevii

Wheat (Triticum aestivum) is one of the most important food crops with an urgent need for increase in its production to feed the growing world. Triticum timopheevii (2n = 4x = 28) is an allotetraploid wheat wild relative species containing the At and G genomes that has been exploited in many pre-breeding programmes for wheat improvement. In this study, we report the generation of a chromosome-scale reference genome assembly of T. timopheevii accession PI 94760 based on PacBio HiFi reads and chromosome conformation capture (Hi-C). The assembly comprised a total size of 9.35 Gb, featuring a contig N50 of 42.4 Mb and included the mitochondrial and plastid genome sequences. Genome annotation predicted 166,325 gene models including 70,365 genes with high confidence. DNA methylation analysis showed that the G genome had on average more methylated bases than the At genome. In summary, the T. timopheevii genome assembly provides a valuable resource for genome-informed discovery of agronomically important genes for food security.

Genetic dissection of the roles of β-hydroxylases in carotenoid metabolism, photosynthesis, and plant growth in tetraploid wheat (Triticum turgidum L.)

Functional redundancy and subfunctionalization of β-hydroxylases in tetraploid wheat tissues open up opportunities for manipulation of carotenoid metabolism for trait improvement. The genetic diversity provided by subgenome homoeologs in allopolyploid wheat can be leveraged for developing improved wheat varieties with modified chemical traits, including profiles of carotenoids, which play critical roles in photosynthesis, photoprotection, and growth regulation. Carotenoid profiles are greatly influenced by hydroxylation catalyzed by β-hydroxylases (HYDs). To genetically dissect the contribution of HYDs to carotenoid metabolism and wheat growth and yield, we isolated loss-of-function mutants of the two homoeologs of HYD1 (HYD-A1 and HYD-B1) and HYD2 (HYD-A2 and HYD-B2) from the sequenced ethyl methanesulfonate mutant population of the tetraploid wheat cultivar Kronos, and generated various mutant combinations. Although functional redundancy between HYD1 and HYD2 paralogs was observed in leaves, HYD1 homoeologs contributed more than HYD2 homoeologs to carotenoid β-ring hydroxylation in this tissue. By contrast, the HYD2 homoeologs functioned toward production of lutein, the major carotenoid in mature grains, whereas HYD1 homoeologs had a limited role. These results collectively suggested subfunctionalization of HYD genes and homoeologs in different tissues of tetraploid wheat. Despite reduced photoprotective responses observed in the triple hyd-A1 hyd-B1 hyd-A2 and the quadruple hyd-A1 hyd-B1 hyd-A2 hyd-B2 combinatorial mutants, comprehensive plant phenotyping analysis revealed that all mutants analyzed were comparable to the control for growth, yield, and fertility, except for a slight delay in anthesis and senescence as well as accelerated germination in the quadruple mutant. Overall, this research takes steps toward untangling the functions of HYDs in wheat and has implications for improving performance and consumer traits of this economically important global crop.

Identification of Fusarium Head Blight Resistance in Triticum timopheevii Accessions and Characterization of Wheat-T. timopheevii Introgression Lines for Enhanced Resistance

A diverse panel of wheat wild relative species was screened for resistance to Fusarium head blight (FHB) by spray inoculation. The great majority of species and accessions were susceptible or highly susceptible to FHB. Accessions of Triticum timopheevii (P95-99.1-1), Agropyron desertorum (9439957), and Elymus vaillantianus (531552) were highly resistant to FHB while additional accessions of T. timopheevii were found to be susceptible to FHB. A combination of spray and point inoculation assessments over two consecutive seasons indicated that the resistance in accession P95-99.1-1 was due to enhanced resistance to initial infection of the fungus (type 1 resistance), and not to reduction in spread (type 2 resistance). A panel of wheat-T. timopheevii (accession P95-99.1-1) introgression lines was screened for FHB resistance over two consecutive seasons using spray inoculation. Most introgression lines were similar in susceptibility to FHB as the wheat recipient (Paragon) but substitution of the terminal portion of chromosome 3BS of wheat with a similar-sized portion of 3G of T. timopheevii significantly enhanced FHB resistance in the wheat background.

Assessing the Heat Tolerance of Meiosis in Spanish Landraces of Tetraploid Wheat Triticum turgidum

Heat stress alters the number and distribution of meiotic crossovers in wild and cultivated plant species. Hence, global warming may have a negative impact on meiosis, fertility, and crop productions. Assessment of germplasm collections to identify heat-tolerant genotypes is a priority for future crop improvement. Durum wheat, Triticum turgidum, is an important cultivated cereal worldwide and given the genetic diversity of the durum wheat Spanish landraces core collection, we decided to analyse the heat stress effect on chiasma formation in a sample of 16 landraces of T. turgidum ssp. turgidum and T. turgidum ssp. durum, from localities with variable climate conditions. Plants of each landrace were grown at 18–22 °C and at 30 °C during the premeiotic temperature-sensitive stage. The number of chiasmata was not affected by heat stress in three genotypes, but decreased by 0.3–2 chiasmata in ten genotypes and more than two chiasmata in the remaining three ones. Both thermotolerant and temperature-sensitive genotypes were found in the two subspecies, and in some of the agroecological zones studied, which supports that genotypes conferring a heat tolerant meiotic phenotype are not dependent on subspecies or geographical origin. Implications of heat adaptive genotypes in future research and breeding are discussed.

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.

Introgression of the Triticum timopheevii Genome Into Wheat Detected by Chromosome-Specific Kompetitive Allele Specific PCR Markers

Triticum timopheevii (2n = 28, A ^t A ^t GG) is a tetraploid wild relative species with great potential to increase the genetic diversity of hexaploid wheat Triticum aestivum (2n = 42, AABBDD) for various important agronomic traits. A breeding scheme that propagated advanced backcrossed populations of wheat-T. timopheevii introgression lines through further backcrossing and self-fertilisation resulted in the generation of 99 introgression lines (ILs) that carried 309 homozygous segments from the A ^t and G subgenomes of T. timopheevii. These introgressions contained 89 and 74 unique segments from the A ^t and G subgenomes, respectively. These overlapping segments covered 98.9% of the T. timopheevii genome that has now been introgressed into bread wheat cv. Paragon including the entirety of all T. timopheevii chromosomes via varying sized segments except for chromosomes 3A ^t , 4G, and 6G. Homozygous ILs contained between one and eight of these introgressions with an average of three per introgression line. These homozygous introgressions were detected through the development of a set of 480 chromosome-specific Kompetitive allele specific PCR (KASP) markers that are well-distributed across the wheat genome. Of these, 149 were developed in this study based on single nucleotide polymorphisms (SNPs) discovered through whole genome sequencing of T. timopheevii. A majority of these KASP markers were also found to be T. timopheevii subgenome specific with 182 detecting A ^t subgenome and 275 detecting G subgenome segments. These markers showed that 98% of the A ^t segments had recombined with the A genome of wheat and 74% of the G genome segments had recombined with the B genome of wheat with the rest recombining with the D genome of wheat. These results were validated through multi-colour in situ hybridisation analysis. Together these homozygous wheat-T. timopheevii ILs and chromosome-specific KASP markers provide an invaluable resource to wheat breeders for trait discovery to combat biotic and abiotic stress factors affecting wheat production due to climate change.

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.

The transcription factor TaLAX1 interacts with Q to antagonistically regulate grain threshability and spike morphogenesis in bread wheat

The domestication gene Q is largely responsible for the widespread cultivation of wheat because it confers multiple domestication traits. However, the underlying molecular mechanisms of how Q regulates these domestication traits remain unclear. In this study, we identify a Q-interacting protein TaLAX1, a basic helix-loop-helix transcription factor, through yeast two-hybrid assays. Using biochemical and genetic approaches, we explore the roles of TaLAX1 in regulating wheat domestication traits. Overexpression of TaLAX1 produces phenotypes, reminiscent of the q allele; loss-of-function Talax1 mutations confer compact spikes, largely similar to the Q-overexpression wheat lines. The two transcription factors TaLAX1 and Q disturb each other's activity to antagonistically regulate the expression of the lignin biosynthesis-related gene TaKNAT7-4D. More interestingly, a natural variation (InDel, +/- TATA), which occurs in the promoter of TaLAX1, is associated with the promoter activity difference between the D subgenome of bread wheat and its ancestor Aegilops tauschii accession T093. This study reveals that the transcription factor TaLAX1 physically interacts with Q to antagonistically regulate wheat domestication traits and a natural variation (InDel, +/- TATA) is associated with the diversification of TaLAX1 promoter activity.

Preferential Subgenome Elimination and Chromosomal Structural Changes Occurring in Newly Formed Tetraploid Wheat-Aegilops ventricosa Amphiploid (AABBDvDvNvNv)

Artificial allopolyploids derived from the genera Triticum and Aegilops have been used as genetic resources for wheat improvement and are a classic example of evolution via allopolyploidization. In this study, we investigated chromosomes and subgenome transmission behavior in the newly formed allopolyploid of wheat group via multicolor Fluorescence in situ hybridization (mc-FISH), using pSc119.2, pTa535, and (GAA)₇ as probe combinations, to enabled us to precisely identify individual chromosomes in 381 S₃ and S₄ generations plants derived from reciprocal crosses between Ae. ventricosa (D^vD^vN^vN^v) and T. turgidum (AABB). A higher rate of aneuploidy, constituting 66.04–86.41% individuals, was observed in these two early generations. Of the four constituent subgenomes, D^v showed the highest frequency of elimination, followed by N^v and B, while A was the most stable. In addition, structural chromosomal changes occurred ubiquitously in the selfed progenies of allopolyploids. Among the constituent subgenomes, B showed the highest number of aberrations. In terms of chromosomal dynamics, there was no significant association between the chromosomal behavior model and the cytoplasm, with the exception of chromosomal loss in the D^v subgenome. The chromosome loss frequency in the D^v subgenome was significantly higher in the T. turgidum × Ae. ventricosa cross than in the Ae. ventricosa × T. turgidum cross. This result indicates that, although the D subgenome showed great instability, allopolyploids containing D subgenome could probably be maintained after a certain hybridization in which the D subgenome donor was used as the maternal parent at its onset stage. Our findings provide valuable information pertaining to the behavior patterns of subgenomes during allopolyploidization. Moreover, the allopolyploids developed here could be used as potential resources for the genetic improvement of wheat.

Genome-wide introgression from a bread wheat × Lophopyrum elongatum amphiploid into wheat

We introgressed wheatgrass germplasm from the octoploid amphiploid Triticum aestivum× Lophopyrum elongatum into wheat by manipulating the wheat Ph1 gene and discovered and characterized 130 introgression lines harboring single or, in various combinations, complete and recombined L. elongatum chromosomes. Diploid wheatgrass Lophopyrum elongatum (genomes EE) possesses valuable traits for wheat genetics and breeding. We evaluated several strategies for introgression of this germplasm into wheat. To detect it, we developed and validated multiplexed sets of Sequenom MassARRAY single nucleotide polymorphism (SNP) markers, which differentiated disomic and monosomic L. elongatum chromosomes from wheat chromosomes. We identified 130 introgression lines (ILs), which harbored 108 complete and 89 recombined L. elongatum chromosomes. Of the latter, 59 chromosomes were recombined by one or more crossovers and 30 were involved in centromeric (Robertsonian) translocations or were telocentric. To identify wheat chromosomes substituted for or recombined with L. elongatum chromosomes, we genotyped the ILs with the wheat 90-K Infinium SNP array. We found that most of the wheat 90-K probes correctly detected their targets in the L. elongatum genome and showed that some wheat SNPs are ancient and had originated prior to the divergence of the wheat and L. elongatum lineages. Of the 130 ILs, 52% were homozygous for Ph1 deletion and thus are staged to be recombined further. We failed to detect in the L. elongatum genome the 4/5 reciprocal translocation that has been reported in Thinopyrum bessarabicum and several other Triticeae genomes.

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.

Recombination between homoeologous chromosomes induced in durum wheat by the Aegilops speltoides Su1-Ph1 suppressor

Su1-Ph1, which we previously introgressed into wheat from Aegilops speltoides, is a potent suppressor of Ph1 and a valuable tool for gene introgression in tetraploid wheat. We previously introgressed Su1-Ph1, a suppressor of the wheat Ph1 gene, from Aegilops speltoides into durum wheat cv Langdon (LDN). Here, we evaluated the utility of the introgressed suppressor for inducing introgression of alien germplasm into durum wheat. We built LDN plants heterozygous for Su1-Ph1 that simultaneously contained a single LDN chromosome 5B and a single Ae. searsii chromosome 5S^se, which targeted them for recombination. We genotyped 28 BC₁F₁ and 84 F₂ progeny with the wheat 90-K Illumina single-nucleotide polymorphism assay and detected extensive recombination between the two chromosomes, which we confirmed by non-denaturing fluorescence in situ hybridization (ND-FISH). We constructed BC₁F₁ and F₂ genetic maps that were 65.31 and 63.71 cM long, respectively. Recombination rates between the 5B and 5S^se chromosomes were double the expected rate computed from their meiotic pairing, which we attributed to selection against aneuploid gametes. Recombination rate between 5B and 5S^se was depressed compared to that between 5B chromosomes in the proximal region of the long arm. We integrated ND-FISH signals into the genetic map and constructed a physical map, which we used to map a 172,188,453-bp Ph1 region. Despite the location of the region in a low-recombination region of the 5B chromosome, we detected three crossovers in it. Our data show that Su1-Ph1 is a valuable tool for gene introgression and gene mapping based on recombination between homoeologous chromosomes in wheat.

The novel function of the Ph1 gene to differentiate homologs from homoeologs evolved in Triticum turgidum ssp. dicoccoides via a dramatic meiosis-specific increase in the expression of the 5B copy of the C-Ph1 gene

The Ph1 gene is the principal regulator of homoeologous chromosome pairing control (HECP) that ensures the diploid-like meiotic chromosome pairing behavior of polyploid wheat. The HECP control was speculated to have evolved after the first event of polyploidization. With the objective to accurately understand the evolution of the HECP control, wild emmer wheat accessions previously known to differ for HECP control were characterized for the structure and expression of the candidate Ph1 gene, C-Ph1. The C-TdPh1-5A and 5B gene copies of emmer wheat showed 98 and 99% DNA sequence similarity respectively with the corresponding hexaploid wheat copies. Further, the C-TdPh1-5B carried the C-Ph1-5B specific structural changes and transcribed three splice variants as observed in the hexaploid wheat. Further, single nucleotide changes differentiating accessions varying for HECP control were identified. Analyzed by quantitative expression analysis, the wild emmer accessions with HECP control showed ~ 10,000-fold higher transcript abundance of the C-TdPh1-5B copy during prophase-I compared to accessions lacking the control. Differential transcriptional regulation of C-TdPh1-5B splice variants further revealed that C-Ph1-5B^alt1 variant is mainly responsible for differential accumulation of C-Ph1-5B copy in accessions with HECP control. Taken together, these results showed that the HECP control evolved via transcriptional regulation of splice variants during meiosis.

Alien chromosome segment from Aegilops speltoides and Dasypyrum villosum increases drought tolerance in wheat via profuse and deep root system

BACKGROUND

Recurrent drought associated with climate change is a major constraint to wheat (Triticum aestivum L.) productivity. This study aimed to (i) quantify the effects of addition/substitution/translocation of chromosome segments from wild relatives of wheat on the root, physiological and yield traits of hexaploid wheat under drought, and (ii) understand the mechanism(s) associated with drought tolerance or susceptibility in wheat-alien chromosome lines.

METHODS

A set of 48 wheat-alien chromosome lines (addition/substitution/translocation lines) with Chinese Spring background were used. Seedling root traits were studied on solid agar medium. To understand the influence of drought on the root system of adult plants, these 48 lines were grown in 150-cm columns for 65 d under full irrigation or withholding water for 58 d. To quantify the effect of drought on physiological and yield traits, the 48 lines were grown in pots under full irrigation until anthesis; after that, half of the plants were drought stressed by withholding water for 16 d before recording physiological and yield-associated traits.

RESULTS

The alien chromosome lines exhibited altered root architecture and decreased photochemical efficiency and seed yield and its components under drought. The wheat-alien chromosome lines T5DS·5S#3L (TA5088) with a chromosome segment from Aegilops speltoides (5S) and T5DL.5 V#3S (TA5638) with a chromosome segment from Dasypyrum villosum (5 V) were identified as drought tolerant, and the drought tolerance mechanism was associated with a deep, thin and profuse root system.

CONCLUSIONS

The two germplasm lines (TA5088 and TA5638) could be used in wheat breeding programs to improve drought tolerance in wheat and understand the underlying molecular genetic mechanisms of root architecture and drought tolerance.

Btr1-A Induces Grain Shattering and Affects Spike Morphology and Yield-Related Traits in Wheat

Spike brittleness represents an important domestication trait in crops. Although the brittle rachis of wild wheat was cloned, however, the molecular mechanism underlying spike brittleness is yet to be elucidated. Here, we identified a single dominant brittle rachis gene Br-Ab on chromosome arm 3AbS using an F2 population of diploid wheat and designated Btr1-Ab. Sequence analysis of the Btr1-A gene in 40 diploid wheat accessions, 80 tetraploid wheat accessions and 38 hexaploid wheat accessions showed that two independent mutations (Ala119Thr for diploid and Gly97* for polyploids) in the Btr1-A coding region resulting in the nonbrittle rachis allele. Overexpression of Btr1-Ab in nonbrittle hexaploid wheat led to brittle rachis in transgenic plants. RNA-Seq analysis revealed that Btr1-A represses the expression of cell wall biosynthesis genes during wheat rachis development. In addition, we found that Btr1-A can modify spike morphology and reduce threshability, grain size and thousand grain weight in transgenic wheat. These results demonstrated that Btr1-A reduces cell wall synthesis in rachis nodes, resulting in natural spikelet shattering, and that the transition from Btr1-A to btr1-A during wheat domestication had profound effects on evolution of spike morphology and yield-related traits.

Nucleotide diversity patterns at the DREB1 transcriptional factor gene in the genome donor species of wheat (Triticum aestivum L)

Bread wheat (AABBDD) originated from the diploid progenitor Triticum urartu (AA), a relative of Aegilops speltoides (BB), and Ae. tauschii (DD). The DREB1 transcriptional factor plays key regulatory role in low-temperature tolerance. The modern breeding strategies resulted in serious decrease of the agricultural biodiversity, which led to a loss of elite genes underlying abiotic stress tolerance in crops. However, knowledge of this gene's natural diversity is largely unknown in the genome donor species of wheat. We characterized the dehydration response element binding protein 1 (DREB1) gene-diversity pattern in Ae. speltoides, Ae. tauschii, T. monococcum and T. urartu. The highest nucleotide diversity value was detected in Ae. speltoides, followed by Ae. tauschii and T. monococcum. The lowest nucleotide diversity value was observed in T. urartu. Nucleotide diversity and haplotype data might suggest no reduction of nucleotide diversity during T. monococcum domestication. Alignment of the 68 DREB1 sequences found a large-size (70 bp) insertion/deletion in the accession PI486264 of Ae. speltoides, which was different from the copy of sequences from other accessions of Ae. speltoides, suggesting a likely existence of two different ancestral Ae. speltoides forms. Implication of sequences variation of Ae. speltoides on origination of B genome in wheat was discussed.

Broadening the bread wheat D genome

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

Aegilops tauschii Genome Sequence: A Framework for Meta-analysis of Wheat QTLs

Numerous quantitative trait loci (QTL) have been mapped in tetraploid and hexaploid wheat and wheat relatives, mostly with simple sequence repeat (SSR) or single nucleotide polymorphism (SNP) markers. To conduct meta-analysis of QTL requires projecting them onto a common genomic framework, either a consensus genetic map or genomic sequence. The latter strategy is pursued here. Of 774 QTL mapped in wheat and wheat relatives found in the literature, 585 (75.6%) were successfully projected onto the Aegilops tauschii pseudomolecules. QTL mapped with SNP markers were more successfully projected (92.2%) than those mapped with SSR markers (66.2%). The QTL were not distributed homogeneously along chromosome arms. Their frequencies increased in the proximal-to-distal direction but declined in the most distal regions and were weakly correlated with recombination rates along the chromosome arms. Databases for projected SSR markers and QTL were constructed and incorporated into the Ae. tauschii JBrowse. To facilitate meta-QTL analysis, eight clusters of QTL were used to estimate standard deviations ([Formula: see text]) of independently mapped QTL projected onto the Ae. tauschii genome sequence. The standard deviations [Formula: see text] were modeled as an exponential decay function of recombination rates along the Ae. tauschii chromosomes. We implemented four hypothesis tests for determining the membership of query QTL. The hypothesis tests and estimation procedure for [Formula: see text] were implemented in a web portal for meta-analysis of projected QTL. Twenty-one QTL for Fusarium head blight resistance mapped on wheat chromosomes 3A, 3B, and 3D were analyzed to illustrate the use of the portal for meta-QTL analyses.

Reassessment of the evolution of wheat chromosomes 4A, 5A, and 7B

Comparison of genome sequences of wild emmer wheat and Aegilops tauschii suggests a novel scenario of the evolution of rearranged wheat chromosomes 4A, 5A, and 7B. Past research suggested that wheat chromosome 4A was subjected to a reciprocal translocation T(4AL;5AL)1 that occurred in the diploid progenitor of the wheat A subgenome and to three major rearrangements that occurred in polyploid wheat: pericentric inversion Inv(4AS;4AL)1, paracentric inversion Inv(4AL;4AL)1, and reciprocal translocation T(4AL;7BS)1. Gene collinearity along the pseudomolecules of tetraploid wild emmer wheat (Triticum turgidum ssp. dicoccoides, subgenomes AABB) and diploid Aegilops tauschii (genomes DD) was employed to confirm these rearrangements and to analyze the breakpoints. The exchange of distal regions of chromosome arms 4AS and 4AL due to pericentric inversion Inv(4AS;4AL)1 was detected, and breakpoints were validated with an optical Bionano genome map. Both breakpoints contained satellite DNA. The breakpoints of reciprocal translocation T(4AL;7BS)1 were also found. However, the breakpoints that generated paracentric inversion Inv(4AL;4AL)1 appeared to be collocated with the 4AL breakpoints that had produced Inv(4AS;4AL)1 and T(4AL;7BS)1. Inv(4AS;4AL)1, Inv(4AL;4AL)1, and T(4AL;7BS)1 either originated sequentially, and Inv(4AL;4AL)1 was produced by recurrent chromosome breaks at the same breakpoints that generated Inv(4AS;4AL)1 and T(4AL;7BS)1, or Inv(4AS;4AL)1, Inv(4AL;4AL)1, and T(4AL;7BS)1 originated simultaneously. We prefer the latter hypothesis since it makes fewer assumptions about the sequence of events that produced these chromosome rearrangements.

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.

Introgression of Aegilops speltoides segments in Triticum aestivum and the effect of the gametocidal genes

BACKGROUND AND AIMS

Bread wheat (Triticum aestivum) has been through a severe genetic bottleneck as a result of its evolution and domestication. It is therefore essential that new sources of genetic variation are generated and utilized. This study aimed to generate genome-wide introgressed segments from Aegilops speltoides. Introgressions generated from this research will be made available for phenotypic analysis.

METHODS

Aegilops speltoides was crossed as the male parent to T. aestivum 'Paragon'. The interspecific hybrids were then backcrossed to Paragon. Introgressions were detected and characterized using the Affymetrix Axiom Array and genomic in situ hybridization (GISH).

KEY RESULTS

Recombination in the gametes of the F1 hybrids was at a level where it was possible to generate a genetic linkage map of Ae. speltoides. This was used to identify 294 wheat/Ae. speltoides introgressions. Introgressions from all seven linkage groups of Ae. speltoides were found, including both large and small segments. Comparative analysis showed that overall macro-synteny is conserved between Ae. speltoides and T. aestivum, but that Ae. speltoides does not contain the 4A/5A/7B translocations present in wheat. Aegilops speltoides has been reported to carry gametocidal genes, i.e. genes that ensure their transmission through the gametes to the next generation. Transmission rates of the seven Ae. speltoides linkage groups introgressed into wheat varied. A 100 % transmission rate of linkage group 2 demonstrates the presence of the gametocidal genes on this chromosome.

CONCLUSIONS

A high level of recombination occurs between the chromosomes of wheat and Ae. speltoides, leading to the generation of large numbers of introgressions with the potential for exploitation in breeding programmes. Due to the gametocidal genes, all germplasm developed will always contain a segment from Ae. speltoides linkage group 2S, in addition to an introgression from any other linkage group.

Molecular cytogenetic and genomic analyses reveal new insights into the origin of the wheat B genome

This work pinpointed the goatgrass chromosomal segment in the wheat B genome using modern cytogenetic and genomic technologies, and provided novel insights into the origin of the wheat B genome. Wheat is a typical allopolyploid with three homoeologous subgenomes (A, B, and D). The donors of the subgenomes A and D had been identified, but not for the subgenome B. The goatgrass Aegilops speltoides (genome SS) has been controversially considered a possible candidate for the donor of the wheat B genome. However, the relationship of the Ae. speltoides S genome with the wheat B genome remains largely obscure. The present study assessed the homology of the B and S genomes using an integrative cytogenetic and genomic approach, and revealed the contribution of Ae. speltoides to the origin of the wheat B genome. We discovered noticeable homology between wheat chromosome 1B and Ae. speltoides chromosome 1S, but not between other chromosomes in the B and S genomes. An Ae. speltoides-originated segment spanning a genomic region of approximately 10.46 Mb was detected on the long arm of wheat chromosome 1B (1BL). The Ae. speltoides-originated segment on 1BL was found to co-evolve with the rest of the B genome. Evidently, Ae. speltoides had been involved in the origin of the wheat B genome, but should not be considered an exclusive donor of this genome. The wheat B genome might have a polyphyletic origin with multiple ancestors involved, including Ae. speltoides. These novel findings will facilitate genome studies in wheat and other polyploids.

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.

Development of a set of PCR markers specific to Aegilops longissima chromosome arms and application in breeding a translocation line

Transcriptome data were used to develop 134 Aegilops longissima specific PCR markers and their comparative maps were constructed by contrasting with the homologous genes in the wheat B genome. Three wheat- Ae. longissima 1BL·1S ^l S translocation lines were identified using the correspondence markers. Aegilops longissima is an important wild species of common wheat that harbors many genes that can be used to improve various traits of common wheat (Triticum aestivum L.). To efficiently transfer the traits conferred by these Ae. longissima genes into wheat, we sequenced the whole expression transcript of Ae. longissima. Using the transcriptome data, we developed 134 specific polymerase chain reaction markers located on the 14 chromosome arms of Ae. longissima. These novel molecular markers were assigned to specific chromosome locations based on a comparison with the homologous genes in the B genome of wheat. Annotation of these genes showed that most had functions related to metabolic processes, hydrolase activity, or catalytic activity. Additionally, we used these markers to identify three wheat-Ae. longissima 1BL·1S^lS translocation lines in somatic variation populations resulting from a cross between wheat cultivar Westonia and a wheat-Ae. longissima substitution line 1S^l(1B). The translocation lines had several low molecular weight glutenin subunits encoding genes beneficial to flour processing quality that came from Ae. longissima 1S^lS. The three translocation lines were also confirmed by genomic in situ hybridization. These translocation lines will be further evaluated for potential quality improvement of bread-making properties of wheat.

Introgression of the Aegilops speltoides Su1-Ph1 Suppressor into Wheat

Meiotic pairing between homoeologous chromosomes in polyploid wheat is inhibited by the Ph1 locus on the long arm of chromosome 5 in the B genome. Aegilops speltoides (genomes SS), the closest relative of the progenitor of the wheat B genome, is polymorphic for genetic suppression of Ph1. Using this polymorphism, two major suppressor loci, Su1-Ph1 and Su2-Ph1, have been mapped in Ae. speltoides. Su1-Ph1 is located in the distal, high-recombination region of the long arm of the Ae. speltoides chromosome 3S. Its location and tight linkage to marker Xpsr1205-3S makes Su1-Ph1 a suitable target for introgression into wheat. Here, Xpsr1205-3S was introgressed into hexaploid bread wheat cv. Chinese Spring (CS) and from there into tetraploid durum wheat cv. Langdon (LDN). Sequential fluorescence in situ hybridization and genomic in situ hybridization showed that an Ae. speltoides segment with Xpsr1205-3S replaced the distal end of the long arm of chromosome 3A. In the CS genetic background, the chromosome induced homoeologous chromosome pairing in interspecific hybrids with Ae. peregrina but not in progenies from crosses involving alien disomic substitution lines. In the LDN genetic background, the chromosome induced homoeologous chromosome pairing in both interspecific hybrids and progenies from crosses involving alien disomic substitution lines. We conclude that the recombined chromosome harbors Su1-Ph1 but its expression requires expression of complementary gene that is present in LDN but absent in CS. We suggest that it is unlikely that Su1-Ph1 and ZIP4-1, a paralog of Ph1 located on wheat chromosomes 3A and 3B and Ae. tauschii chromosome 3D, are equivalent. The utility of Su1-Ph1 for induction of recombination between homoeologous chromosomes in wheat is illustrated.

Fine organization of genomic regions tagged to the 5S rDNA locus of the bread wheat 5B chromosome

BACKGROUND

The multigene family encoding the 5S rRNA, one of the most important structurally-functional part of the large ribosomal subunit, is an obligate component of all eukaryotic genomes. 5S rDNA has long been a favored target for cytological and phylogenetic studies due to the inherent peculiarities of its structural organization, such as the tandem arrays of repetitive units and their high interspecific divergence. The complex polyploid nature of the genome of bread wheat, Triticum aestivum, and the technically difficult task of sequencing clusters of tandem repeats mean that the detailed organization of extended genomic regions containing 5S rRNA genes remains unclear. This is despite the recent progress made in wheat genomic sequencing. Using pyrosequencing of BAC clones, in this work we studied the organization of two distinct 5S rDNA-tagged regions of the 5BS chromosome of bread wheat.

RESULTS

Three BAC-clones containing 5S rDNA were identified in the 5BS chromosome-specific BAC-library of Triticum aestivum. Using the results of pyrosequencing and assembling, we obtained six 5S rDNA- containing contigs with a total length of 140,417 bp, and two sets (pools) of individual 5S rDNA sequences belonging to separate, but closely located genomic regions on the 5BS chromosome. Both regions are characterized by the presence of approximately 70-80 copies of 5S rDNA, however, they are completely different in their structural organization. The first region contained highly diverged short-type 5S rDNA units that were disrupted by multiple insertions of transposable elements. The second region contained the more conserved long-type 5S rDNA, organized as a single tandem array. FISH using probes specific to both 5S rDNA unit types showed differences in the distribution and intensity of signals on the chromosomes of polyploid wheat species and their diploid progenitors.

CONCLUSION

A detailed structural organization of two closely located 5S rDNA-tagged genomic regions on the 5BS chromosome of bread wheat has been established. These two regions differ in the organization of both 5S rDNA and the neighboring sequences comprised of transposable elements, implying different modes of evolution for these regions.

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.

Major structural genomic alterations can be associated with hybrid speciation in Aegilops markgrafii (Triticeae)

During evolutionary history many grasses from the tribe Triticeae have undergone interspecific hybridization, resulting in allopolyploidy; whereas homoploid hybrid speciation was found only in rye. Homoeologous chromosomes within the Triticeae preserved cross-species macrocolinearity, except for a few species with rearranged genomes. Aegilops markgrafii, a diploid wild relative of wheat (2n = 2x = 14), has a highly asymmetrical karyotype that is indicative of chromosome rearrangements. Molecular cytogenetics and next-generation sequencing were used to explore the genome organization. Fluorescence in situ hybridization with a set of wheat cDNAs allowed the macrostructure and cross-genome homoeology of the Ae. markgrafii chromosomes to be established. Two chromosomes maintained colinearity, whereas the remaining were highly rearranged as a result of inversions and inter- and intrachromosomal translocations. We used sets of barley and wheat orthologous gene sequences to compare discrete parts of the Ae. markgrafii genome involved in the rearrangements. Analysis of sequence identity profiles and phylogenic relationships grouped chromosome blocks into two distinct clusters. Chromosome painting revealed the distribution of transposable elements and differentiated chromosome blocks into two groups consistent with the sequence analyses. These data suggest that introgressive hybridization accompanied by gross chromosome rearrangements might have had an impact on karyotype evolution and homoploid speciation in Ae. markgrafii.

Physical Map of the Short Arm of Bread Wheat Chromosome 3D

Bread wheat ( L.) is one of the most important crops worldwide. Although a reference genome sequence would represent a valuable resource for wheat improvement through genomics-assisted breeding and gene cloning, its generation has long been hampered by its allohexaploidy, high repeat content, and large size. As a part of a project coordinated by the International Wheat Genome Sequencing Consortium (IWGSC), a physical map of the short arm of wheat chromosome 3D (3DS) was prepared to facilitate reference genome assembly and positional gene cloning. It comprises 869 contigs with a cumulative length of 274.5 Mbp and represents 85.5% of the estimated chromosome arm size. Eighty-six Mbp of survey sequences from chromosome arm 3DS were assigned in silico to physical map contigs via next-generation sequencing of bacterial artificial chromosome pools, thus providing a high-density framework for physical map ordering along the chromosome arm. About 60% of the physical map was anchored in this single experiment. Finally, 1393 high-confidence genes were anchored to the physical map. Comparisons of gene space of the chromosome arm 3DS with genomes of closely related species [ (L.) P.Beauv., rice ( L.), and sorghum [ (L.) Moench] and homeologous wheat chromosomes provided information about gene movement on the chromosome arm.

Solar Radiation-Associated Adaptive SNP Genetic Differentiation in Wild Emmer Wheat, Triticum dicoccoides

Whole-genome scans with large number of genetic markers provide the opportunity to investigate local adaptation in natural populations and identify candidate genes under positive selection. In the present study, adaptation genetic differentiation associated with solar radiation was investigated using 695 polymorphic SNP markers in wild emmer wheat originated in a micro-site at Yehudiyya, Israel. The test involved two solar radiation niches: (1) sun, in-between trees; and (2) shade, under tree canopy, separated apart by a distance of 2-4 m. Analysis of molecular variance showed a small (0.53%) but significant portion of overall variation between the sun and shade micro-niches, indicating a non-ignorable genetic differentiation between sun and shade habitats. Fifty SNP markers showed a medium (0.05 ≤ FST ≤ 0.15) or high genetic differentiation (FST > 0.15). A total of 21 outlier loci under positive selection were identified by using four different FST -outlier testing algorithms. The markers and genome locations under positive selection are consistent with the known patterns of selection. These results suggested that genetic differentiation between sun and shade habitats is substantial, radiation-associated, and therefore ecologically determined. Hence, the results of this study reflected effects of natural selection through solar radiation on EST-related SNP genetic diversity, resulting presumably in different adaptive complexes at a micro-scale divergence. The present work highlights the evolutionary theory and application significance of solar radiation-driven natural selection in wheat improvement.

A step change in the transfer of interspecific variation into wheat from Amblyopyrum muticum

Despite some notable successes, only a fraction of the genetic variation available in wild relatives has been utilized to produce superior wheat varieties. This is as a direct result of the lack of availability of suitable high-throughput technologies to detect wheat/wild relative introgressions when they occur. Here, we report on the use of a new SNP array to detect wheat/wild relative introgressions in backcross progenies derived from interspecific hexaploid wheat/Ambylopyrum muticum F₁ hybrids. The array enabled the detection and characterization of 218 genomewide wheat/Am. muticum introgressions, that is a significant step change in the generation and detection of introgressions compared to previous work in the field. Furthermore, the frequency of introgressions detected was sufficiently high to enable the construction of seven linkage groups of the Am. muticum genome, thus enabling the syntenic relationship between the wild relative and hexaploid wheat to be determined. The importance of the genetic variation from Am. muticum introduced into wheat for the development of superior varieties is discussed.

Reconciling the evolutionary origin of bread wheat (Triticum aestivum)

The origin of bread wheat (Triticum aestivum; AABBDD) has been a subject of controversy and of intense debate in the scientific community over the last few decades. In 2015, three articles published in New Phytologist discussed the origin of hexaploid bread wheat (AABBDD) from the diploid progenitors Triticum urartu (AA), a relative of Aegilops speltoides (BB) and Triticum tauschii (DD). Access to new genomic resources since 2013 has offered the opportunity to gain novel insights into the paleohistory of modern bread wheat, allowing characterization of its origin from its diploid progenitors at unprecedented resolution. We propose a reconciled evolutionary scenario for the modern bread wheat genome based on the complementary investigation of transposable element and mutation dynamics between diploid, tetraploid and hexaploid wheat. In this scenario, the structural asymmetry observed between the A, B and D subgenomes in hexaploid bread wheat derives from the cumulative effect of diploid progenitor divergence, the hybrid origin of the D subgenome, and subgenome partitioning following the polyploidization events.

TaGW2, a Good Reflection of Wheat Polyploidization and Evolution

Hexaploid wheat consists of three subgenomes, namely, A, B, and D. These well-characterized ancestral genomes also exist at the diploid and tetraploid levels, thereby rendering wheat as a good model species for studying polyploidization. Here, we performed intra- and inter-species comparative analyses of wheat and its relatives to dissect polymorphism and differentiation of the TaGW2 genes. Our results showed that genetic diversity of TaGW2 decreased with progression from the diploids to tetraploids and hexaploids. The strongest selection occurred in the promoter regions of TaGW2-6A and TaGW2-6B. Phylogenetic trees clearly indicated that Triticum urartu and Ae. speltoides were the donors of the A and B genomes in tetraploid and hexaploid wheats. Haplotypes detected among hexaploid genotypes traced back to the tetraploid level. Fst and π values revealed that the strongest selection on TaGW2 occurred at the tetraploid level rather than in hexaploid wheat. This infers that grain size enlargement, especially increased kernel width, mainly occurred in tetraploid genotypes. In addition, relative expression levels of TaGW2s significantly declined from the diploid level to tetraploids and hexaploids, further indicating that these genes negatively regulate kernel size. Our results also revealed that the polyploidization events possibly caused much stronger differentiation than domestication and breeding.

A set of Triticum aestivum-Aegilops speltoides Robertsonian translocation lines

Here we report the production of a set of wheat - Aegilops speltoides Robertsonian translocations covering all Ae. speltoides chromosome arms except the long arm of the homoeologous group 4 chromosome. Aegilops speltoides of the Poaceae family is the most probable donor of the B and G genomes of polyploid Triticum species and also an important source of resistance to diseases and pests of wheat. Previously, we reported the production of a complete set of T aestivum-Ae. speltoides chromosome addition lines and a set of disomic S(B/A)-genome chromosome substitution lines. The isolation of compensating Robertsonian translocations (RobTs) composed of alien chromosome arms translocated to homoeologous wheat chromosome arms is the important next step to exploit the genetic variation of a wild relative of wheat. Here, we report the development of molecular markers specific for the S-genome chromosomes and their use in the isolation of a set of 13 compensating wheat-Ae. speltoides RobTs covering the S genome of Ae. speltoides except for the long arm of chromosome 4S. Most of the RobTs were fully fertile and will facilitate mapping of genes to specific chromosome arms and also will accelerate the introgression of agronomically useful traits from Ae. speltoides into wheat by homologous recombination.

The occurrence of spring forms in tetraploid Timopheevi wheat is associated with variation in the first intron of the VRN-A1 gene

BACKGROUND

Triticum araraticum and Triticum timopheevii are tetraploid species of the Timopheevi group. The former includes both winter and spring forms with a predominance of winter forms, whereas T. timopheevii is considered a spring species. In order to clarify the origin of the spring growth habit in T. timopheevii, allelic variability of the VRN-1 gene was investigated in a set of accessions of both tetraploid species, together with the diploid species Ae. speltoides, presumed donor of the G genome to these tetraploids.

RESULTS

The promoter region of the VRN-A1 locus in all studied tetraploid accessions of both T. araraticum and T. timopheevii represents the previously described allele VRN-A1f with a 50 bp deletion near the start codon. Three additional alleles were identified namely, VRN-A1f-del, VRN-A1f-ins and VRN-A1f-del/ins, which contained large mutations in the first (1^st) intron of VRN-A1. The first allele, carrying a deletion of 2.7 kb in a central part of intron 1, occurred in a few accessions of T. araraticum and no accessions of T. timopheevii. The VRN-A1f-ins allele, containing the insertion of a 0.4 kb MITE element about 0.4 kb upstream from the start of intron 1, and allele VRN-A1f-del/ins having this insertion coupled with a deletion of 2.7 kb are characteristic only for T. timopheevii. Allelic variation at the VRN-G1 locus includes the previously described allele VRN-G1a (with the insertion of a 0.2 kb MITE in the promoter) found in a few accessions of both tetraploid species. We showed that alleles VRN-A1f-del and VRN-G1a have no association with the spring growth habit, while in all accessions of T. timopheevii this habit was associated with the dominant VRN-A1f-ins and VRN-A1f-del/ins alleles. None of the Ae. speltoides accessions included in this study had changes in the promoter or 1^st intron regions of VRN-1 which might confer a spring growth habit. The VRN-1 promoter sequences analyzed herein and downloaded from databases have been used to construct a phylogram to assess the time of divergence of Ae. speltoides in relation to other wheat species.

CONCLUSIONS

Among accessions of T. araraticum, the preferentially winter predecessor of T. timopheevii, two large mutations were found in both VRN-A1 and VRN-G1 loci (VRN-A1f-del and VRN-G1a) that were found to have no effect on vernalization requirements. Spring tetraploid T. timopheevii had one VRN-1 allele in common for two species (VRN-G1a), and two that were specific (VRN-A1f-ins, VRN-A1f-del/ins). The latter alleles include mutations in the 1^st intron of VRN-A1 and also share a 0.4 kb MITE insertion near the start of intron 1. We suggested that this insertion resulted in a spring growth habit in a progenitor of T. timopheevii which has probably been selected during subsequent domestication. The phylogram constructed on the basis of the VRN-1 promoter sequences confirmed the early divergence (~3.5 MYA) of the ancestor(s) of the B/G genomes from Ae. speltoides.

Sequencing and comparative analyses of Aegilops tauschii chromosome arm 3DS reveal rapid evolution of Triticeae genomes

Bread wheat (Triticum aestivum, AABBDD) is an allohexaploid species derived from two rounds of interspecific hybridizations. A high-quality genome sequence assembly of diploid Aegilops tauschii, the donor of the wheat D genome, will provide a useful platform to study polyploid wheat evolution. A combined approach of BAC pooling and next-generation sequencing technology was employed to sequence the minimum tiling path (MTP) of 3176 BAC clones from the short arm of Ae. tauschii chromosome 3 (At3DS). The final assembly of 135 super-scaffolds with an N50 of 4.2 Mb was used to build a 247-Mb pseudomolecule with a total of 2222 predicted protein-coding genes. Compared with the orthologous regions of rice, Brachypodium, and sorghum, At3DS contains 38.67% more genes. In comparison to At3DS, the short arm sequence of wheat chromosome 3B (Ta3BS) is 95-Mb large in size, which is primarily due to the expansion of the non-centromeric region, suggesting that transposable element (TE) bursts in Ta3B likely occurred there. Also, the size increase is accompanied by a proportional increase in gene number in Ta3BS. We found that in the sequence of short arm of wheat chromosome 3D (Ta3DS), there was only less than 0.27% gene loss compared to At3DS. Our study reveals divergent evolution of grass genomes and provides new insights into sequence changes in the polyploid wheat genome.

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.

Diversification of the Homoeologous Lr34 Sequences in Polyploid Wheat Species and Their Diploid Progenitors

Allopolyploidization induces a multiple processes of genomic reorganization, including the structurally functional diversification of the homoeologous genes. An example of such diversification is the appearance of the Lr34 gene on chromosome 7D of bread wheat T. aestivum (BAD), the gene conferring durable, race non-specific protection against three fungal pathogens. In this study, we focused on the variability of a functionally critical region between exons 10-12 of Lr34 among diploid progenitors of wheat genomes and their respective polyploids. In the diploid A-genome species, two basic forms of the studied region have been revealed: (1) non-functional forms containing stop codons, or/and frameshifts (T. monococcum/T. urartu) and (2) forms with no such a mutations (T. boeoticum). The Lr34 sequence of T. urartu containing a TGA stop codon was inherited by the first tetraploid T. dicoccoides (BA), and then reorganized in some accessions of this species due to the insertion of an LTR retroelement in exon 10. Besides T. boeoticum, the second form of the Lr34 sequence is also characteristic of A. speltoides, which presumably donated this form to all polyploid descendants bearing B-genome. No differences were found between the D-genome-specific Lr34 sequences studied here and downloaded from databases, implying the highest level of conservation of the Lr34 predecessor throughout evolution. The sequence data were later used to construct phylograms, and apparent peculiarities in the evolution of the studied region of Lr34 genes discussed.

To Have and to Hold: Selection for Seed and Fruit Retention During Crop Domestication

Crop domestication provides a useful model system to characterize the molecular and developmental bases of morphological variation in plants. Among the most universal changes resulting from selection during crop domestication is the loss of seed and fruit dispersal mechanisms, which greatly facilitates harvesting efficiency. In this review, we consider the molecular genetic and developmental bases of the loss of seed shattering and fruit dispersal in six major crop plant families, three of which are primarily associated with seed crops (Poaceae, Brassicaceae, Fabaceae) and three of which are associated with fleshy-fruited crops (Solanaceae, Rosaceae, Rutaceae). We find that the developmental basis of the loss of seed/fruit dispersal is conserved in a number of independently domesticated crops, indicating the widespread occurrence of developmentally convergent evolution in response to human selection. With regard to the molecular genetic approaches used to characterize the basis of this trait, traditional biparental quantitative trait loci mapping remains the most commonly used strategy; however, recent advances in next-generation sequencing technologies are now providing new avenues to map and characterize loss of shattering/dispersal alleles. We anticipate that continued application of these approaches, together with candidate gene analyses informed by known shattering candidate genes from other crops, will lead to a rapid expansion of our understanding of this critical domestication trait.

Genome-Wide Comparative Analysis of Flowering-Related Genes in Arabidopsis, Wheat, and Barley

Early flowering is an important trait influencing grain yield and quality in wheat (Triticum aestivum L.) and barley (Hordeum vulgare L.) in short-season cropping regions. However, due to large and complex genomes of these species, direct identification of flowering genes and their molecular characterization remain challenging. Here, we used a bioinformatic approach to predict flowering-related genes in wheat and barley from 190 known Arabidopsis (Arabidopsis thaliana (L.) Heynh.) flowering genes. We identified 900 and 275 putative orthologs in wheat and barley, respectively. The annotated flowering-related genes were clustered into 144 orthologous groups with one-to-one, one-to-many, many-to-one, and many-to-many orthology relationships. Our approach was further validated by domain and phylogenetic analyses of flowering-related proteins and comparative analysis of publicly available microarray data sets for in silico expression profiling of flowering-related genes in 13 different developmental stages of wheat and barley. These further analyses showed that orthologous gene pairs in three critical flowering gene families (PEBP, MADS, and BBX) exhibited similar expression patterns among 13 developmental stages in wheat and barley, suggesting similar functions among the orthologous genes with sequence and expression similarities. The predicted candidate flowering genes can be confirmed and incorporated into molecular breeding for early flowering wheat and barley in short-season cropping regions.

Molecular characterization of two isoforms of a farnesyl pyrophosphate synthase gene in wheat and their roles in sesquiterpene synthesis and inducible defence against aphid infestation

Aphids are important pests of wheat (Triticum aestivum) that affect crop production globally. Herbivore-induced emission of sesquiterpenes can repel pests, and farnesyl pyrophosphate synthase (FPS) is a key enzyme involved in sesquiterpene biosynthesis. However, fps orthologues in wheat and their functional roles in sesquiterpene synthesis and defence against aphid infestation are unknown. Here, two fps isoforms, Tafps1 and Tafps2, were identified in wheat. Quantitative real-time polymerase chain reaction (qRT-PCR) and in vitro catalytic activity analyses were conducted to investigate expression patterns and activity. Heterologous expression of these isoforms in Arabidopsis thaliana, virus-induced gene silencing (VIGS) in wheat and aphid behavioural assays were performed to understand the functional roles of these two isoforms. We demonstrated that Tafps1 and Tafps2 played different roles in induced responses to aphid infestation and in sesquiterpene synthesis. Heterologous expression in A. thaliana resulted in repulsion of the peach aphid (Myzus persicae). Wheat plants with these two isoforms transiently silenced were significantly attractive to grain aphid (Sitobion avenae). Our results provide new insights into induced defence against aphid herbivory in wheat, in particular, the different roles of the two Tafps isoforms in both sesquiterpene biosynthesis and defence against aphid infestation.

Phylogenetic analysis of Gossypium L. using restriction fragment length polymorphism of repeated sequences

Cotton is the world's leading textile fiber crop and is also grown as a bioenergy and food crop. Knowledge of the phylogeny of closely related species and the genome origin and evolution of polyploid species is significant for advanced genomics research and breeding. We have reconstructed the phylogeny of the cotton genus, Gossypium L., and deciphered the genome origin and evolution of its five polyploid species by restriction fragment analysis of repeated sequences. Nuclear DNA of 84 accessions representing 35 species and all eight genomes of the genus were analyzed. The phylogenetic tree of the genus was reconstructed using the parsimony method on 1033 polymorphic repeated sequence restriction fragments. The genome origin of its polyploids was determined by calculating the diploid-polyploid restriction fragment correspondence (RFC). The tree is consistent with the morphological classification, genome designation and geographic distribution of the species at subgenus, section and subsection levels. Gossypium lobatum (D7) was unambiguously shown to have the highest RFC with the D-subgenomes of all five polyploids of the genus, while the common ancestor of Gossypium herbaceum (A1) and Gossypium arboreum (A2) likely contributed to the A-subgenomes of the polyploids. These results provide a comprehensive phylogenetic tree of the cotton genus and new insights into the genome origin and evolution of its polyploid species. The results also further demonstrate a simple, rapid and inexpensive method suitable for phylogenetic analysis of closely related species, especially congeneric species, and the inference of genome origin of polyploids that constitute over 70 % of flowering plants.

VRN-1 gene- associated prerequisites of spring growth habit in wild tetraploid wheat T. dicoccoides and the diploid A genome species

BACKGROUND

In order to clarify the origin of spring growth habit in modern domesticated wheat, allelic variability of the VRN-1 gene was investigated in a wide set of accessions of the wild tetraploid species Triticum dicoccoides (BBAA), together with diploid species T. monococcum, T. boeoticum and T. urartu, presumable donors of the A genome to polyploid wheats.

RESULTS

No significant variation was found at the VRN-B1 locus of T. dicoccoides, whereas at VRN-A1 a number of previously described alleles were found with small deletions in the promoter (VRN-A1b, VRN-A1d) or a large deletion in the first (1st) intron (VRN-A1L). The diploid A genome species were characterized by their own set of VRN-1 alleles including previously described VRN-A1f and VRN-A1h alleles with deletions in the promoter region and the VRN-A1ins allele containing a 0.5 kb insertion in the 1st intron. Based on the CAPS screening data, alleles VRN-A1f and VRN-A1ins were species-specific for T. monococcum, while allele VRN-A1h was specific for T. boeoticum. Different indels were revealed in both the promoter and 1(st) intron of the recessive VRN-A1u allele providing specific identification of T. urartu, the proposed donor of the A genome to modern wheat. We found that alleles VRN-A1b and VRN-A1h, previously described as dominant, have either no or weak association with spring growth habit, while in some diploid accessions this habit was associated with the recessive VRN-A1 allele.

CONCLUSIONS

Spring growth habit in diploid wheats was only partially associated with indels in regulatory regions of the VRN-1 gene. An exception is T. monococcum where dominant mutations in both the promoter region and, especially, the 1st intron were selected during domestication resulting in a greater variety of spring forms. The wild tetraploid T. dicoccoides had a distinct set of VRN-A1 alleles compared to the diploids in this study, indicating an independent origin of spring tetraploid forms that likely occurred after combining of diploid genomes. These alleles were subsequently inherited by cultivated polyploid (tetraploid and hexaploid) descendants.

Analysis of agronomic and domestication traits in a durum × cultivated emmer wheat population using a high-density single nucleotide polymorphism-based linkage map

Development of a high-density SNP map and evaluation of QTL shed light on domestication events in tetraploid wheat and the potential utility of cultivated emmer wheat for durum wheat improvement. Cultivated emmer wheat (Triticum turgidum ssp. dicoccum) is tetraploid and considered as one of the eight founder crops that spawned the Agricultural Revolution about 10,000 years ago. Cultivated emmer has non-free-threshing seed and a somewhat fragile rachis, but mutations in genes governing these and other agronomic traits occurred that led to the formation of today's fully domesticated durum wheat (T. turgidum ssp. durum). Here, we evaluated a population of recombinant inbred lines (RILs) derived from a cross between a cultivated emmer accession and a durum wheat variety. A high-density single nucleotide polymorphism (SNP)-based genetic linkage map consisting of 2,593 markers was developed for the identification of quantitative trait loci. The major domestication gene Q had profound effects on spike length and compactness, rachis fragility, and threshability as expected. The cultivated emmer parent contributed increased spikelets per spike, and the durum parent contributed higher kernel weight, which led to the identification of some RILs that had significantly higher grain weight per spike than either parent. Threshability was governed not only by the Q locus, but other loci as well including Tg-B1 on chromosome 2B and a putative Tg-A1 locus on chromosome 2A indicating that mutations in the Tg loci occurred during the transition of cultivated emmer to the fully domesticated tetraploid. These results not only shed light on the events that shaped wheat domestication, but also demonstrate that cultivated emmer is a useful source of genetic variation for the enhancement of durum varieties.

The chloroplast view of the evolution of polyploid wheat

Polyploid wheats comprise four species: Triticum turgidum (AABB genomes) and T. aestivum (AABBDD) in the Emmer lineage, and T. timopheevii (AAGG) and T. zhukovskyi (AAGGA(m) A(m) ) in the Timopheevi lineage. Genetic relationships between chloroplast genomes were studied to trace the evolutionary history of the species. Twenty-five chloroplast genomes were sequenced, and 1127 plant accessions were genotyped, representing 13 Triticum and Aegilops species. The A. speltoides (SS genome) diverged before the divergence of T. urartu (AA), A. tauschii (DD) and the Aegilops species of the Sitopsis section. Aegilops speltoides forms a monophyletic clade with the polyploid Emmer and Timopheevi wheats, which originated within the last 0.7 and 0.4 Myr, respectively. The geographic distribution of chloroplast haplotypes of the wild tetraploid wheats and A. speltoides illustrates the possible geographic origin of the Emmer lineage in the southern Levant and the Timopheevi lineage in northern Iraq. Aegilops speltoides is the closest relative of the diploid donor of the chloroplast (cytoplasm), as well as the B and G genomes to Timopheevi and Emmer lineages. Chloroplast haplotypes were often shared by species or subspecies within major lineages and between the lineages, indicating the contribution of introgression to the evolution and domestication of polyploid wheats.

Mapping resistance to the bird cherry-oat aphid and the greenbug in wheat using sequence-based genotyping

Identification of novel resistance QTL against wheat aphids. First QTL-resistance report for R. padi in wheat and chromosome 2DL for S. graminum . These sources have potential use in wheat breeding. The aphids Rhopalosiphum padi and Schizaphis graminum are important pests of common wheat (Triticum aestivum L.). Characterization of the genetic bases of resistance sources is crucial to facilitate the development of resistant wheat cultivars to these insects. We examined 140 recombinant inbred lines (RILs) from the cross of Seri M82 wheat (susceptible) with the synthetic hexaploid wheat CWI76364 (resistant). RILs were phenotyped for R. padi antibiosis and tolerance traits. Phenotyping of S. graminum resistance was based on leaf chlorosis in a greenhouse screening and the number of S. graminum/tiller in the field. RILs were also scored for pubescence. Using a sequence-based genotyping method, we located genomic regions associated with these resistance traits. A quantitative trait locus (QTL) for R. padi antibiosis (QRp.slu.4BL) that explained 10.2 % of phenotypic variation was found in chromosome 4BL and located 14.6 cM apart from the pubescence locus. We found no association between plant pubescence and the resistance traits. We found two QTLs for R. padi tolerance (QRp.slu.5AL and QRp.slu.5BL) in chromosomes 5AL and 5BL, with an epistatic interaction between a locus in chromosome 3AL (EnQRp.slu.5AL) and QRp.slu.5AL. These genomic regions explained about 35 % of the phenotypic variation. We re-mapped a previously reported gene for S. graminum resistance (putatively Gba) in 7DL and found a novel QTL associated with the number of aphids/tiller (QGb.slu-2DL) in chromosome 2DL. This is the first report on the genetic mapping of R. padi resistance in wheat and the first report where chromosome 2DL is shown to be associated with S. graminum resistance.

SNP Discovery for mapping alien introgressions in wheat

Background

Monitoring alien introgressions in crop plants is difficult due to the lack of genetic and molecular mapping information on the wild crop relatives. The tertiary gene pool of wheat is a very important source of genetic variability for wheat improvement against biotic and abiotic stresses. By exploring the 5M^g short arm (5M^gS) of Aegilops geniculata, we can apply chromosome genomics for the discovery of SNP markers and their use for monitoring alien introgressions in wheat (Triticum aestivum L).

Results

The short arm of chromosome 5M^g of Ae. geniculata Roth (syn. Ae. ovata L.; 2n = 4x = 28, U^gU^gM^gM^g) was flow-sorted from a wheat line in which it is maintained as a telocentric chromosome. DNA of the sorted arm was amplified and sequenced using an Illumina Hiseq 2000 with ~45x coverage. The sequence data was used for SNP discovery against wheat homoeologous group-5 assemblies. A total of 2,178 unique, 5M^gS-specific SNPs were discovered. Randomly selected samples of 59 5M^gS-specific SNPs were tested (44 by KASPar assay and 15 by Sanger sequencing) and 84% were validated. Of the selected SNPs, 97% mapped to a chromosome 5M^g addition to wheat (the source of t5M^gS), and 94% to 5M^g introgressed from a different accession of Ae. geniculata substituting for chromosome 5D of wheat. The validated SNPs also identified chromosome segments of 5M^gS origin in a set of T5D-5M^g translocation lines; eight SNPs (25%) mapped to TA5601 [T5DL · 5DS-5M^gS(0.75)] and three (8%) to TA5602 [T5DL · 5DS-5M^gS (0.95)]. SNPs (gsnp_5ms83 and gsnp_5ms94), tagging chromosome T5DL · 5DS-5M^gS(0.95) with the smallest introgression carrying resistance to leaf rust (Lr57) and stripe rust (Yr40), were validated in two released germplasm lines with Lr57 and Yr40 genes.

Conclusion

This approach should be widely applicable for the identification of species/genome-specific SNPs. The development of a large number of SNP markers will facilitate the precise introgression and monitoring of alien segments in crop breeding programs and further enable mapping and cloning novel genes from the wild relatives of crop plants.

Increased genomic prediction accuracy in wheat breeding through spatial adjustment of field trial data

In crop breeding, the interest of predicting the performance of candidate cultivars in the field has increased due to recent advances in molecular breeding technologies. However, the complexity of the wheat genome presents some challenges for applying new technologies in molecular marker identification with next-generation sequencing. We applied genotyping-by-sequencing, a recently developed method to identify single-nucleotide polymorphisms, in the genomes of 384 wheat (Triticum aestivum) genotypes that were field tested under three different water regimes in Mediterranean climatic conditions: rain-fed only, mild water stress, and fully irrigated. We identified 102,324 single-nucleotide polymorphisms in these genotypes, and the phenotypic data were used to train and test genomic selection models intended to predict yield, thousand-kernel weight, number of kernels per spike, and heading date. Phenotypic data showed marked spatial variation. Therefore, different models were tested to correct the trends observed in the field. A mixed-model using moving-means as a covariate was found to best fit the data. When we applied the genomic selection models, the accuracy of predicted traits increased with spatial adjustment. Multiple genomic selection models were tested, and a Gaussian kernel model was determined to give the highest accuracy. The best predictions between environments were obtained when data from different years were used to train the model. Our results confirm that genotyping-by-sequencing is an effective tool to obtain genome-wide information for crops with complex genomes, that these data are efficient for predicting traits, and that correction of spatial variation is a crucial ingredient to increase prediction accuracy in genomic selection models.