A re-evaluation of the homoploid hybrid origin of Aegilops tauschii, the donor of the wheat D-subgenome

The phylogeny of the Triticeae: Resolution and phylogenetic conflict based on genomewide nuclear loci

PREMISE

The wheat tribe, Triticeae, has been the subject of molecular phylogenetic analyses for nearly three decades, and extensive phylogenetic conflict has been apparent from the earliest comparisons among DNA-based data sets. While most previous analyses focused primarily on nuclear vs. chloroplast DNA conflict, the present analysis provides a broader picture of conflict among nuclear loci throughout the tribe.

METHODS

Exon data were generated from over 1000 nuclear loci using targeted sequence capture with custom baits, and nearly complete chloroplast genome sequences were recovered. Phylogenetic conflict was assessed among the trees from the chloroplast genomes, the concatenated nuclear loci, and a series of nuclear-locus subsets guided by Hordeum chromosome gene maps.

RESULTS

At the intergeneric level, the analyses collectively revealed a few broadly consistent relationships. However, the prevailing pattern was one of extensive phylogenetic conflict throughout the tribe, among both deep and shallow branches, and with the extent of the conflict varying among data subsets.

CONCLUSIONS

The results suggest continual introgression or lineage sorting within and among the named lineages of the Triticeae, shaping both deep and shallow relationships in the tribe.

Characterization of the complete chloroplast genome of Elymus alashanicus (Keng) S. L. Chen, and its phylogenetic analysis

Elymus alashanicus (Keng) S. L. Chen, a herbaceous plant endemic to China, plays a crucial role in the local ecosystems. In this study, we sequenced and characterized the complete chloroplast (cp) genome of E. alashanicus, which is 135,072 bp in length and arranged in a circular form. The cp genome includes a pair of inverted repeats (IRa and IRb) of 20,813 bp each, separated by a large single-copy (LSC) region of 80,678 bp and a small single-copy (SSC) region of 12,768 bp. The cp genome contains 130 genes, including 83 protein-coding genes, 39 tRNA genes, and eight rRNA genes. Phylogenetic analysis revealed that E. alashanicus is closely related to Elymus breviaristatus and Campeiostachys dahurica var. tangutorum in current sampling. Our findings provide valuable insights into the cp genome of E. alashanicus, which could contribute to further studies on the evolution and conservation of this species.

Genome shock in a synthetic allotetraploid wheat invokes subgenome-partitioned gene regulation, meiotic instability, and karyotype variation

It is becoming increasingly evident that interspecific hybridization at the homoploid level or coupled with whole-genome duplication (i.e. allopolyploidization) has played a major role in biological evolution. However, the direct impacts of hybridization and allopolyploidization on genome structure and function, phenotype, and fitness remains to be fully understood. Synthetic hybrids and allopolyploids are trackable experimental systems that can be used to address this issue. In this study, we resynthesized a pair of reciprocal F1 hybrids and corresponding reciprocal allotetraploids using the two diploid progenitor species of bread wheat (Triticum aestivum, BBAADD), namely T. urartu (AA) and Aegilops tauschii (DD). By comparing phenotypes related to growth, development, and fitness, and by analysing genome expression in both hybrids and allotetraploids in relation to the parents, we found that the types and trends of karyotype variation in the immediately formed allotetraploids were correlated with both instability of meiosis and chromosome- and subgenome-biased expression. We determined clear advantages of allotetraploids over diploid F1 hybrids in several morphological traits including fitness that mirrored the tissue- and developmental stage-dependent subgenome-partitioning of the allotetraploids. The allotetraploids were meiotically unstable primarily due to homoeologous pairing that varied dramatically among the chromosomes. Nonetheless, the manifestation of organismal karyotype variation and the occurrence of meiotic irregularity were not concordant, suggesting a role of functional constraints probably imposed by subgenome- and chromosome-biased gene expression. Our results provide new insights into the direct impacts and consequences of hybridization and allopolyploidization that are relevant to evolution and likely to be informative for future crop improvement approaches using synthetic polyploids.

Post-hybridization introgression and natural selection promoted genomic divergence of Aegilops speltoides and the four S*-genome diploid species

Understanding how different driving forces have promoted biological divergence and speciation is one of the central issues in evolutionary biology. The Triticum/Aegilops species complex contains 13 diploid species belonging to the A-, B- and D-lineages and offers an ideal system to address the evolutionary dynamics of lineage fusion and splitting. Here, we sequenced the whole genomes of one S-genome species (Aegilops speltoides) of the B-lineage and four S*-genome diploid species (Aegilops bicornis, Aegilops longissima, Aegilops sharonensis and Aegilops searsii) of the D-lineage at the population level. We performed detailed comparisons of the five species and with the other four representative A-, B- and D-lineage species. Our estimates identified frequent genetic introgressions from A- and B-lineages to the D-lineage species. A remarkable observation is the contrasting distributions of putative introgressed loci by the A- and B-lineages along all the seven chromosomes to the extant D-lineage species. These genetic introgressions resulted in high levels of genetic divergence at centromeric regions between Ae. speltoides (B-lineage) and the other four S*-genome diploid species (D-lineage), while natural selection is a potential contributor to divergence among the four S*-genome species at telomeric regions. Our study provides a genome-wide view on how genetic introgression and natural selection acted together yet chromosome-regionally divided to promote genomic divergence among the five S- and S*-genome diploid species, which provides new and nuanced insights into the evolutionary history of the Triticum/Aegilops species complex.

Comparative genomic analysis of 5Mg chromosome of Aegilops geniculata and 5Uu chromosome of Aegilops umbellulata reveal genic diversity in the tertiary gene pool

Wheat is one of the most important cereal crops for the global food security. Due to its narrow genetic base, modern bread wheat cultivars face challenges from increasing abiotic and biotic stresses. Since genetic improvement is the most sustainable approach, finding novel genes and alleles is critical for enhancing the genetic diversity of wheat. The tertiary gene pool of wheat is considered a gold mine for genetic diversity as novel genes and alleles can be identified and transferred to wheat cultivars. Aegilops geniculata and Ae. umbellulata are the key members of the tertiary gene pool of wheat and harbor important genes against abiotic and biotic stresses. Homoeologous-group five chromosomes (5U^u and 5M^g) have been extensively studied from Ae. geniculata and Ae. umbellulata as they harbor several important genes including Lr57, Lr76, Yr40, Yr70, Sr53 and chromosomal pairing loci. In the present study, using chromosome DNA sequencing and RNAseq datasets, we performed comparative analysis to study homoeologous gene evolution in 5M^g, 5U^u, and group 5 wheat chromosomes. Our findings highlight the diversity of transcription factors and resistance genes, resulting from the differential expansion of the gene families. Both the chromosomes were found to be enriched with the "response to stimulus" category of genes providing resistance against biotic and abiotic stress. Phylogenetic study positioned the M genome closer to the D genome, with higher proximity to the A genome than the B genome. Over 4000 genes were impacted by SNPs on 5D, with 4-5% of those genes displaying non-disruptive variations that affect gene function.

Advances in the Mining of Disease Resistance Genes from Aegilops tauschii and the Utilization in Wheat

Aegilops tauschii is one of the malignant weeds that affect wheat production and is also the wild species ancestor of the D genome of hexaploid wheat (Triticum aestivum, AABBDD). It contains many disease resistance genes that have been lost in the long-term evolution of wheat and is an important genetic resource for the mining and utilization of wheat disease resistance genes. In recent years, the genome sequence of Aegilops tauschii has been preliminarily completed, which has laid a good foundation for the further exploration of wheat disease resistance genes in Aegilops tauschii. There are many studies on disease resistance genes in Aegilops tauschii; in order to provide better help for the disease resistance breeding of wheat, this paper analyzes and reviews the relationship between Aegilops tauschii and wheat, the research progress of Aegilops tauschii, the discovery of disease resistance genes from Aegilops tauschii, and the application of disease resistance genes from Aegilops tauschii to modern wheat breeding, providing a reference for the further exploration and utilization of Aegilops tauschii in wheat disease resistance breeding.

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

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

The impacts of allopolyploidization on Methyl-CpG-Binding Domain (MBD) gene family in Brassica napus

BACKGROUND

Polyploidization promotes species formation and is widespread in angiosperms. Genome changes dramatically bring opportunities and challenges to plants after polyploidy. Methyl-CpG-Binding Domain (MBD) proteins can recognize and bind to methylation sites and they play an important role in the physiological process related to methylation in animals and plants. However, research on the influence of the allopolyploidization process on the MBD gene family is still lacking, so it is necessary to conduct a comprehensive analysis.

RESULTS

In this study, twenty-two, ten and eleven MBD genes were identified in the genome of allotetraploid B. napus and its diploid ancestors, B. rapa and B. oleracea, respectively. Based on the clades of the MBD gene in Arabidopsis, rice and maize, we divided the new phylogenetic tree into 8 clades. Among them, the true MBD genes in Brassica existed in only 5 clades. Clade IV and Clade VI were unique in term of MBD genes in dicotyledons. Ka/Ks calculations showed that MBD genes underwent purifying selection in Brassica and may retain genes through sequence or functional differentiation early in evolution. In the process of allopolyploidization, the number of MBD gene introns increased, and the protein motifs changed. The MBD proteins had their own special motifs in each clade, and the MBD domains were only conserved in their clades. At the same time, the MBD genes were expressed in flower, leaf, silique, and stem tissues, and the expression levels of the different genes were significantly different, while the tissue specificity was not obvious. The allopolyploidization process may increase the number of cis-acting elements and activate the transposable elements. During allopolyploidization, the expression pattern of the MBD gene changes, which may be regulated by cis-acting elements and transposable elements. The number imbalance of cis-acting elements and transposable elements in A_n and C_n subgenomes may also lead to biased A_n subgenome expression of the MBD gene in B. napus.

CONCLUSIONS

In this study, by evaluating the number, structure, phylogeny and expression of the MBD gene in B. napus and its diploid ancestors, we increased the understanding of MBD genes in allopolyploids and provided a reference for future analysis of allopolyploidization.

Phylotranscriptomics Resolves the Phylogeny of Pooideae and Uncovers Factors for Their Adaptive Evolution

Adaptation to cool climates has occurred several times in different angiosperm groups. Among them, Pooideae, the largest grass subfamily with ∼3,900 species including wheat and barley, have successfully occupied many temperate regions and play a prominent role in temperate ecosystems. To investigate possible factors contributing to Pooideae adaptive evolution to cooling climates, we performed phylogenetic reconstruction using five gene sets (with 1,234 nuclear genes and their subsets) from 157 transcriptomes/genomes representing all 15 tribes and 24 of 26 subtribes. Our phylogeny supports the monophyly of all tribes (except Diarrheneae) and all subtribes with at least two species, with strongly supported resolution of their relationships. Molecular dating suggests that Pooideae originated in the late Cretaceous, with subsequent divergences under cooling conditions first among many tribes from the early middle to late Eocene and again among genera in the middle Miocene and later periods. We identified a cluster of gene duplications (CGD5) shared by the core Pooideae (with 80% Pooideae species) near the Eocene–Oligocene transition, coinciding with the transition from closed to open habitat and an upshift of diversification rate. Molecular evolutionary analyses homologs of CBF for cold resistance uncovered tandem duplications during the core Pooideae history, dramatically increasing their copy number and possibly promoting adaptation to cold habitats. Moreover, duplication of AP1/FUL-like genes before the Pooideae origin might have facilitated the regulation of the vernalization pathway under cold environments. These and other results provide new insights into factors that likely have contributed to the successful adaptation of Pooideae members to temperate regions.

Wheat speciation and adaptation: perspectives from reticulate evolution

Reticulate evolution through the interchanging of genetic components across organisms can impact significantly on the fitness and adaptation of species. Bread wheat (Triticum aestivum subsp. aestivum) is one of the most important crops in the world. Allopolyploid speciation, frequent hybridization, extensive introgression, and occasional horizontal gene transfer (HGT) have been shaping a typical paradigm of reticulate evolution in bread wheat and its wild relatives, which is likely to have a substantial influence on phenotypic traits and environmental adaptability of bread wheat. In this review, we outlined the evolutionary history of bread wheat and its wild relatives with a highlight on the interspecific hybridization events, demonstrating the reticulate relationship between species/subspecies in the genera Triticum and Aegilops. Furthermore, we discussed the genetic mechanisms and evolutionary significance underlying the introgression of bread wheat and its wild relatives. An in-depth understanding of the evolutionary process of Triticum species should be beneficial to future genetic study and breeding of bread wheat.

Nuclear-cytoplasmic balance: whole genome duplications induce elevated organellar genome copy number

The plant genome is partitioned across three distinct subcellular compartments: the nucleus, mitochondria, and plastids. Successful coordination of gene expression among these organellar genomes and the nuclear genome is critical for plant function and fitness. Whole genome duplication (WGD) events in the nucleus have played a major role in the diversification of land plants and are expected to perturb the relative copy number (stoichiometry) of nuclear, mitochondrial, and plastid genomes. Thus, elucidating the mechanisms whereby plant cells respond to the cytonuclear stoichiometric imbalance that follows WGDs represents an important yet underexplored question in understanding the evolutionary consequences of genome doubling. We used droplet digital PCR to investigate the relationship between nuclear and organellar genome copy numbers in allopolyploids and their diploid progenitors in both wheat and Arabidopsis. Polyploids exhibit elevated organellar genome copy numbers per cell, largely preserving the cytonuclear stoichiometry observed in diploids despite the change in nuclear genome copy number. To investigate the timescale over which cytonuclear stoichiometry may respond to WGD, we also estimated the organellar genome copy number in Arabidopsis synthetic autopolyploids and in a haploid-induced diploid line. We observed corresponding changes in organellar genome copy number in these laboratory-generated lines, indicating that at least some of the cellular response to cytonuclear stoichiometric imbalance is immediate following WGD. We conclude that increases in organellar genome copy numbers represent a common response to polyploidization, suggesting that maintenance of cytonuclear stoichiometry is an important component in establishing polyploid lineages.

Asymmetric epigenetic modification and homoeolog expression bias in the establishment and evolution of allopolyploid Brassica napus

This study explores how allopolyploidization reshapes the biased expression and asymmetric epigenetic modification of homoeologous gene pairs, and examines the regulation types and epigenetic basis of expression bias. We analyzed the gene expression and four epigenetic modifications (DNA methylation, H3K4me3, H3K27me3 and H3K27ac) of 29 976 homoeologous gene pairs in resynthesized, natural allopolyploid Brassica napus and an in silico 'hybrid'. We comprehensively elucidated the biased gene expression, asymmetric epigenetic modifications and the generational transmission characteristics of these homoeologous gene pairs in B. napus. We analyzed cis/trans effects and the epigenetic basis of homoeolog expression bias. There was a significant positive correlation between two active histone modifications and biased gene expression. We revealed that parental legacy was the dominant principle in the remodeling of homoeolog expression bias and asymmetric epigenetic modifications in B. napus, and further clarified that this depends on whether there were differences in the expression/epigenetic modifications of gene pairs in parents/progenitors. The maternal genome was dominant in the homoeolog expression bias of resynthesized B. napus, and this phenomenon was attenuated in natural B. napus. Furthermore, cis rather than trans effects were dominant when epigenetic modifications potentially affected biased expression of gene pairs in B. napus.

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

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

Pleiotropic function of the SQUAMOSA PROMOTER-BINDING PROTEIN-LIKE gene TaSPL14 in wheat plant architecture

MAIN CONCLUSION

The function of SQUAMOSA PROMOTER-BINDING PROTEIN-BOX gene TaSPL14 in wheat is similar to that of OsSPL14 in rice in regulating plant height, panicle length, spikelet number, and thousand-grain weight of wheat, but differs during tiller development. TaSPL14 may regulate spike development via ethylene-response gene EIN3-LIKE 1 (TaEIL1), ETHYLENE-RESPONSIVE TRANSCRIPTION FACTOR 2.11 (TaRAP2.11), and ETHYLENE-RESPONSIVE TRANSCRIPTION FACTOR 1 (TaERF1), but not DENSE AND ERECT PANICLE 1 (TaDEP1) in wheat. The SQUAMOSA PROMOTER-BINDING PROTEIN-LIKE gene OsSPL14 from rice is considered to be a major determinant of ideal plant architecture consisting of few unproductive tillers, more grains per spike, and high resistance of stems to lodging. However, the function of its orthologous gene, TaSPL14, in wheat is unknown. Here, we reported the functional similarities and differences between TaSPL14 and OsSPL14. Similar to OsSPL14 knock-outs in rice, wheat TaSPL14 knock-out plants exhibited decreased plant height, panicle length, spikelet number, and thousand-grain weight. In contrast to OsSPL14, however, TaSPL14 did not affect tiller number. Transcriptome analysis revealed that the expression of genes related to ethylene response was significantly decreased in young spikes of TaSPL14 knock-out lines as compared with wild type. TaSPL14 directly binds to the promoters of the ethylene-response genes TaEIL1, TaRAP2.11, and TaERF1, and promotes their expression, suggesting that TaSPL14 might regulate wheat spike development via the ethylene-response pathway. The elucidation of TaSPL14 will contribute to understanding of the molecular mechanisms that underlie wheat plant architecture.

Genome-wide sequence information reveals recurrent hybridization among diploid wheat wild relatives

Many conflicting hypotheses regarding the relationships among crops and wild species closely related to wheat (the genera Aegilops, Amblyopyrum, and Triticum) have been postulated. The contribution of hybridization to the evolution of these taxa is intensely discussed. To determine possible causes for this, and provide a phylogeny of the diploid taxa based on genome-wide sequence information, independent data were obtained from genotyping-by-sequencing and a target-enrichment experiment that returned 244 low-copy nuclear loci. The data were analyzed using Bayesian, likelihood and coalescent-based methods. D statistics were used to test if incomplete lineage sorting alone or together with hybridization is the source for incongruent gene trees. Here we present the phylogeny of all diploid species of the wheat wild relatives. We hypothesize that most of the wheat-group species were shaped by a primordial homoploid hybrid speciation event involving the ancestral Triticum and Am. muticum lineages to form all other species except Ae. speltoides. This hybridization event was followed by multiple introgressions affecting all taxa except Triticum. Mostly progenitors of the extant species were involved in these processes, while recent interspecific gene flow seems insignificant. The composite nature of many genomes of wheat-group taxa results in complicated patterns of diploid contributions when these lineages are involved in polyploid formation, which is, for example, the case for tetraploid and hexaploid wheats. Our analysis provides phylogenetic relationships and a testable hypothesis for the genome compositions in the basic evolutionary units within the wheat group of Triticeae.

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

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

Hybridization preceded radiation in diploid wheats

Evolutionary relationships among the Aegilops-Triticum relatives of cultivated wheats have been difficult to resolve owing to incomplete lineage sorting and reticulate evolution. Recent studies have suggested that the wheat D-genome lineage (progenitor of Ae. tauschii) originated through homoploid hybridization between the A-genome lineage (progenitor of Triticum s.str.) and the B-genome lineage (progenitor of Ae. speltoides). This scenario of reticulation has been debated, calling for adequate phylogenetic analyses based on comprehensive sampling. To reconstruct the evolution of Aegilops-Triticum diploids, we here combined high-throughput sequencing of 38 nuclear low-copy loci of multiple accessions of all 13 species with inferences of the species phylogeny using the full-parameterized MCMC_SEQ method. Phylogenies recovered a monophyletic Aegilops-Triticum lineage that began diversifying ~6.6 Ma ago and gave rise to four sublineages, i.e. the A- (2 species), B- (1 species), D- (9 species) and T- (Ae. mutica) genome lineage. Full-parameterized phylogenies as well as patterns of tree dilation and tree compression supported a hybrid origin of the D-genome lineage from A and B ~3.0-4.0 Ma ago, and did not indicate additional hybridization events. Conflicting ABBA-BABA tests suggestive of further reticulation were shown here to result from ancestral population structure rather than hybridization. This comprehensive and dated phylogeny of wheat relatives indicates that the origin of the hybrid D-genome was followed by intense diversification into the majority of extant diploid as well as allopolyploid wild wheats.

Cytonuclear Coevolution following Homoploid Hybrid Speciation in Aegilops tauschii

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

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.

Pervasive hybridizations in the history of wheat relatives

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

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

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

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

Background

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

Perspectives and Conclusion

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

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

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

Revisiting Pivotal-Differential Genome Evolution in Wheat

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

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

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.

Genetic Architecture of Anther Extrusion in Spring and Winter Wheat

Hybrid wheat breeding is gaining prominence worldwide because it ensures higher and more static yield than conventionally bred varieties. The cleistogamous floral architecture of wheat (Triticum aestivum L.) impedes anthers inside the floret, making it largely an inbreeder. For hybrid seed production, high anther extrusion is needed to promote cross pollination and to ensure a high level of pollen availability for the seed plant. This study, therefore, aimed at the genetic dissection of anther extrusion (AE) in panels of spring (SP), and winter wheat (WP) accessions by genome wide association studies (GWAS). We performed GWAS to identify the SNP markers potentially linked with AE in each panel separately. Phenotypic data were collected for 3 years for each panel. The average levels of Pearson's correlation (r) among all years and their best linear unbiased estimates (BLUEs) within both panels were high (r(SP) = 0.75, P < 0.0001;r(WP) = 0.72, P < 0.0001). Genotypic data (with minimum of 0.05 minor allele frequency applied) included 12,066 and 12,191 SNP markers for SP and WP, respectively. Both genotypes and environment influenced the magnitude of AE. In total, 23 significant (|log₁₀(P)| > 3.0) marker trait associations (MTAs) were detected (SP = 11; WP = 12). Anther extrusion behaved as a complex trait with significant markers having either favorable or unfavorable additive effects and imparting minor to moderate levels of phenotypic variance (R²(SP) = 9.75-14.24%; R² (WP) = 9.44-16.98%). All mapped significant markers as well as the markers within their significant linkage disequilibrium (r² ≥ 0.30) regions were blasted against wheat genome assembly (IWGSC1+popseq) to find the corresponding genes and their high confidence descriptions were retrieved. These genes and their orthologs in Hordeum vulgare, Brachypodium distachyon, Oryza sativa, and Sorghum bicolor revealed syntenic genomic regions potentially involved in flowering-related traits. Moreover, the expression data of these genes suggested potential candidates for AE. Our results suggest that the use of significant markers can help to introduce AE in high yielding varieties to increase cross fertilization rates and improve hybrid-seed production in wheat.

Transcriptome asymmetry in synthetic and natural allotetraploid wheats, revealed by RNA-sequencing

Allopolyploidization has occurred frequently within the Triticum-Aegilops complex which provides a suitable system to investigate how allopolyploidization shapes the expression patterns of duplicated homeologs. We have conducted transcriptome-profiling of leaves and young inflorescences in wild and domesticated tetraploid wheats, Triticum turgidum ssp. dicoccoides (BBAA) and ssp. durum (BBAA), an extracted tetraploid (BBAA), and a synthetic tetraploid (S(l) S(l) AA) wheat together with its diploid parents, Aegilops longissima (S(l) S(l) ) and Triticum urartu (AA). The two diploid species showed tissue-specific differences in genome-wide ortholog expression, which plays an important role in transcriptome shock-mediated homeolog expression rewiring and hence transcriptome asymmetry in the synthetic tetraploid. Further changes of homeolog expression apparently occurred in natural tetraploid wheats, which led to novel transcriptome asymmetry between the two subgenomes. In particular, our results showed that extremely biased homeolog expression can occur rapidly upon the allotetraploidzation and this trend is further enhanced in the course of domestication and evolution of polyploid wheats. Our results suggest that allopolyploidization is accompanied by distinct phases of homeolog expression changes, with parental legacy playing major roles in the immediate rewiring of homeolog expression upon allopolyploidization, while evolution and domestication under allotetraploidy drive further homeolog-expression changes toward re-established subgenome expression asymmetry.