Chloroplast and nuclear DNA variation in common wheat: insight into the origin and evolution of common wheat

Paternally imprinted LATE-FLOWERING2 transcription factor contributes to paternal-excess interploidy hybridization barriers in wheat

Interploidy hybridization between hexaploid and tetraploid genotypes occurred repeatedly during genomic introgression events throughout wheat evolution, and is commonly employed in wheat breeding programs. Hexaploid wheat usually serves as maternal parent because the reciprocal cross generates progeny with severe defects and poor seed germination, but the underlying mechanism is poorly understood. Here, we performed detailed analysis of phenotypic variation in endosperm between two interploidy reciprocal crosses arising from tetraploid (Triticum durum, AABB) and hexaploid wheat (Triticum aestivum, AABBDD). In the paternal- versus the maternal-excess cross, the timing of endosperm cellularization was delayed and starch granule accumulation in the endosperm was repressed, causing reduced germination percentage. The expression profiles of genes involved in nutrient metabolism differed strongly between these endosperm types. Furthermore, expression patterns of parental alleles were dramatically disturbed in interploidy versus intraploidy crosses, leading to increased number of imprinted genes. The endosperm-specific TaLFL2 showed a paternally imprinted expression pattern in interploidy crosses partially due to allele-specific DNA methylation. Paternal TaLFL2 binds to and represses a nutrient accumulation regulator TaNAC019, leading to reduced storage protein and starch accumulation during endosperm development in paternal-excess cross, as confirmed by interploidy crosses between tetraploid wild-type and clustered regularly interspaced palindromic repeats (CRISPR) - CRISPR-associated protein 9 generated hexaploid mutants. These findings reveal a contribution of genomic imprinting to paternal-excess interploidy hybridization barriers during wheat evolution history and explains why experienced breeders preferentially exploit maternal-excess interploidy crosses in wheat breeding programs.

Genetic Analysis of Hexaploid Wheat (Triticum aestivum L.) Using the Complete Sequencing of Chloroplast DNA and Haplotype Analysis of the Wknox1 Gene

The aim of the presented study is a genetic characterization of the hexaploid wheat Triticum aestivum L. Two approaches were used for the genealogical study of hexaploid wheats—the complete sequencing of chloroplast DNA and PCR-based haplotype analysis of the fourth intron of Wknox1d and of the fifth-to-sixth-exon region of Wknox1b. The complete chloroplast DNA sequences of 13 hexaploid wheat samples were determined: Free-threshing—T. aestivum subsp. aestivum, one sample; T. aestivum subsp. compactum, two samples; T. aestivum subsp. sphaerococcum, one sample; T. aestivum subsp. carthlicoides, four samples. Hulled—T. aestivum subsp. spelta, three samples; T. aestivum subsp. vavilovii jakubz., two samples. The comparative analysis of complete cpDNA sequences of 20 hexaploid wheat samples (13 samples in this article plus 7 samples sequenced in this laboratory in 2018) was carried out. PCR-based haplotype analysis of the fourth intron of Wknox1d and of the fifth-to-sixth exon region of Wknox1b of all 20 hexaploid wheat samples was carried out. The 20 hexaploid wheat samples (13 samples in this article plus 7 samples in 2018) can be divided into two groups—T. aestivum subsp. spelta, three samples and T. aestivum subsp. vavilovii collected in Armenia, and the remaining 16 samples, including T. aestivum subsp. vavilovii collected in Europe (Sweden). If we take the cpDNA of Chinese Spring as a reference, 25 SNPs can be identified. Furthermore, 13–14 SNPs can be identified in T. aestivum subsp. spelta and subsp. vavilovii (Vav1). In the other samples up to 11 SNPs were detected. 22 SNPs are found in the intergenic regions, 2 found in introns, and 10 SNPs were found in the genes, of which seven are synonymous. PCR-based haplotype analysis of the fourth intron of Wknox1d and the fifth-to-sixth-exon region of Wknox1b provides an opportunity to make an assumption that hexaploid wheats T. aestivum subsp. macha var. palaeocolchicum and var. letshckumicum differ from other macha samples by the absence of a 42 bp insertion in the fourth intron of Wknox1d. One possible explanation for this observation would be that two Aegilops tauschii Coss. (A) and (B) participated in the formation of hexaploids through the D genome: Ae. tauschii (A)—macha (1–5, 7, 8, 10–12), and Ae. tauschii (B)—macha M6, M9, T. aestivum subsp. aestivum cv. ‘Chinese Spring’ and cv. ‘Red Doly’.

Geographical distribution and adaptive variation of VRN-A3 alleles in worldwide polyploid wheat (Triticum spp.) species collection

The distribution of early flowering alleles of VRN-A3 was found to be biased to low latitudes, and these alleles may contribute to environmental adaptability to low latitudes in cultivated emmer wheat. In wheat (Triticum spp.), the flowering time is an important trait for successful seed production and yield by adapting to the regional environment. An early flowering allele of VRN-A3 with 7- and 25-bp insertions in the promoter region (Vrn-A3a-h1) has recently been reported from the analysis of an emmer wheat (Triticum turgidum L. ssp. dicoccum) accession, TN26. This early flowering allele of VRN-A3 might be associated with the regional adaptation of wheat. In this study, we elucidated its geographic distribution to assess the importance of the early flowering allele of VRN-A3 in worldwide wheat collection. From sequence analysis, we identified six VRN-A3 alleles with the 7- and 25-bp insertions, namely, Vrn-A3a-h2, Vrn-A3a-h3, Vrn-A3a-h4, Vrn-A3a-h5, Vrn-A3a-h6, and Vrn-A3c-h2 from wild emmer wheat, while we identified two VRN-A3 alleles with these insertions, Vrn-A3a-h2 and Vrn-A3c-h1 from cultivated tetraploid and hexaploid wheat species in addition to Vrn-A3a-h1. Among VRN-A3 alleles distributed in cultivated wheat, we found that Vrn-A3a-h2 promoted early heading, whereas Vrn-A3c-h1 did not affect heading time. Our analysis showed that the distribution of early flowering alleles of VRN-A3 dominated in cultivated emmer wheat in Ethiopia and India, which actually showed an early flowering phenotype. This implied that the early flowering alleles of VRN-A3 contribute to adaptability to a low-latitude environment in cultivated emmer wheat. We could not find durum (T. turgidum L. ssp. durum) and bread wheat (T. aestivum L. ssp. aestivum) accessions with these early flowering alleles. Our findings indicated that Vrn-A3a-h1 and Vrn-A3a-h2 were useful for breeding of early flowering cultivars in durum and bread wheat varieties.

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.

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.

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.

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

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

Analysis of ATP6 sequence diversity in the Triticum-Aegilops species group reveals the crucial role of rearrangement in mitochondrial genome evolution

Mutation and chromosomal rearrangements are the two main forces of increasing genetic diversity for natural selection to act upon, and ultimately drive the evolutionary process. Although genome evolution is a function of both forces, simultaneously, the ratio of each can be varied among different genomes and genomic regions. It is believed that in plant mitochondrial genome, rearrangements play a more important role than point mutations, but relatively few studies have directly addressed this phenomenon. To address this issue, we isolated and sequenced the ATP6-1 and ATP6-2 genes from 46 different euplasmic and alloplasmic wheat lines. Four different ATP6-1 orthologs were detected, two of them reported for the first time. Expression analysis revealed that all four orthologs are transcriptionally active. Results also indicated that both point mutation and genomic rearrangement are involved in the evolution of ATP6. However, rearrangement is the predominant force that triggers drastic variation. Data also indicated that speciation of domesticated wheat cultivars were simultaneous with the duplication of this gene. These results directly support the notion that rearrangement plays a significant role in driving plant mitochondrial genome evolution.

Sequencing of chloroplast genomes from wheat, barley, rye and their relatives provides a detailed insight into the evolution of the Triticeae tribe

Using Roche/454 technology, we sequenced the chloroplast genomes of 12 Triticeae species, including bread wheat, barley and rye, as well as the diploid progenitors and relatives of bread wheat Triticum urartu, Aegilops speltoides and Ae. tauschii. Two wild tetraploid taxa, Ae. cylindrica and Ae. geniculata, were also included. Additionally, we incorporated wild Einkorn wheat Triticum boeoticum and its domesticated form T. monococcum and two Hordeum spontaneum (wild barley) genotypes. Chloroplast genomes were used for overall sequence comparison, phylogenetic analysis and dating of divergence times. We estimate that barley diverged from rye and wheat approximately 8–9 million years ago (MYA). The genome donors of hexaploid wheat diverged between 2.1–2.9 MYA, while rye diverged from Triticum aestivum approximately 3–4 MYA, more recently than previously estimated. Interestingly, the A genome taxa T. boeoticum and T. urartu were estimated to have diverged approximately 570,000 years ago. As these two have a reproductive barrier, the divergence time estimate also provides an upper limit for the time required for the formation of a species boundary between the two. Furthermore, we conclusively show that the chloroplast genome of hexaploid wheat was contributed by the B genome donor and that this unknown species diverged from Ae. speltoides about 980,000 years ago. Additionally, sequence alignments identified a translocation of a chloroplast segment to the nuclear genome which is specific to the rye/wheat lineage. We propose the presented phylogeny and divergence time estimates as a reference framework for future studies on Triticeae.

A possible breakage of linkage disequilibrium between mitochondrial and chloroplast genomes during Emmer and Dinkel wheat evolution

In wheat (Triticum) and Aegilops, chloroplast and mitochondrial genomes have been studied for over three decades to clarify the phylogenetic relationships among species, and most of the maternal lineages of polyploid species have been clarified. Mitochondrial genomes of Emmer (tetraploid with nuclear genome AABB) and Dinkel (hexaploid with AABBDD) wheat are classified into two different types, VIIa and VIIb, by the presence–absence of the third largest HindIII fragment (named H3) in the mitochondrial DNA. Although the mitochondrial genome in the genera often provides useful information to clarify the phylogenetic relationship among closely related species, the phylogenetic significance of this dimorphism has yet not been clarified. In this study, to facilitate analysis using a large number of accessions, a sequence characterized amplified region (SCAR) marker that distinguishes the type VIIb mitochondrial genome from type VIIa was first developed. Mitochondrial genome type was determined for each of 30 accessions of wild and cultivated Emmer wheat and 25 accessions of Dinkel wheat. The mitochondrial genome type for each accession was compared with the plastogroup that had been determined using chloroplast microsatellite markers. Unexpectedly, the distribution of mitochondrial genome type was not in accordance with that of the plastogroups, suggesting occasional paternal leakage of either the mitochondrial or chloroplast genome during speciation and differentiation of Emmer and Dinkel wheat. An alternative possibility that substoichiometric shifting is involved in the observed dimorphism of the mitochondrial genome is also discussed.

Evolution of polyploid triticum wheats under cultivation: the role of domestication, natural hybridization and allopolyploid speciation in their diversification

The evolution of the polyploid Triticum wheats is distinctive in that domestication, natural hybridization and allopolyploid speciation have all had significant impacts on their diversification. In this review, I outline the phylogenetic relationships of cultivated wheats and their wild relatives and provide an overview of the recent progress and remaining issues in understanding the genetic and ecological factors that favored their evolution. An attempt is made to view the evolution of the polyploid Triticum wheats as a continuous process of diversification that was initiated by domestication of tetraploid emmer wheat and driven by various natural events ranging from interploidy introgression via hybridization to allopolyploid speciation of hexaploid common wheat, instead of viewing it as a group of discrete evolutionary processes that separately proceeded at the tetraploid and hexaploid levels. This standpoint underscores the important role of natural hybridization in the reticulate diversification of the tetraploid-hexaploid Triticum wheat complex and highlights critical, but underappreciated, issues that concern the allopolyploid speciation of common wheat.

Repeated polyploidy drove different levels of crossover suppression between homoeologous chromosomes in Brassica napus allohaploids

Allopolyploid species contain more than two sets of related chromosomes (homoeologs) that must be sorted during meiosis to ensure fertility. As polyploid species usually have multiple origins, one intriguing, yet largely underexplored, question is whether different mechanisms suppressing crossovers between homoeologs may coexist within the same polyphyletic species. We addressed this question using Brassica napus, a young polyphyletic allopolyploid species. We first analyzed the meiotic behavior of 363 allohaploids produced from 29 accessions, which represent a large part of B. napus genetic diversity. Two main clear-cut meiotic phenotypes were observed, encompassing a twofold difference in the number of univalents at metaphase I. We then sequenced two chloroplast intergenic regions to gain insight into the maternal origins of the same 29 accessions; only two plastid haplotypes were found, and these correlated with the dichotomy of meiotic phenotypes. Finally, we analyzed genetic diversity at the PrBn locus, which was shown to determine meiotic behavior in a segregating population of B. napus allohaploids. We observed that segregation of two alleles at PrBn could adequately explain a large part of the variation in meiotic behavior found among B. napus allohaploids. Overall, our results suggest that repeated polyploidy resulted in different levels of crossover suppression between homoeologs in B. napus allohaploids.

Natural variation for fertile triploid F1 hybrid formation in allohexaploid wheat speciation

The tempo, mode, and geography of allopolyploid speciation are influenced by natural variation in the ability of parental species to express postzygotic reproductive phenotypes that affect hybrid fertility. To shed light on the impact of such natural variations, we used allohexaploid Triticum aestivum wheats' evolution as a model and analyzed the geographic and phylogenetic distributions of Aegilops tauschii (diploid progenitor) accessions involved in the expression of abnormality and fertility in triploid F(1) hybrids with Triticum turgidum (tetraploid progenitor). Artificial-cross experiments and chloroplast-DNA-based evolutionary analyses showed that hybrid-abnormality-causing accessions had limited geographic and phylogenetic distributions, indicative that postzygotic hybridization barriers are underdeveloped between these species. In contrast, accessions that are involved with fertile triploid F(1) hybrid formation have wide geographic and phylogenetic distributions, indicative of a deep evolutionary origin. Wide-spread hybrid-fertilizing accessions support the theory that T. aestivum speciation occurred at multiple sites within the species range of Ae. tauschii, in which existing conditions enabled natural hybridization with T. turgidum. Implications of our findings on how natural variation in the ability of Ae. tauschii to express those postzygotic reproductive phenotypes diversified and contributed to the speciation of T. aestivum are discussed.