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

Potential of tetraploid wheats in plant breeding: A review

Domestication syndrome, selection pressure, and modern plant breeding programs have reduced the genetic diversity of the wheat germplasm. For the genetic gains of breeding programs to be sustainable, plant breeders require a diverse gene pool to select genes for resistance to biotic stress factors, tolerance to abiotic stress factors, and improved quality and yield components. Thus, old landraces, subspecies and wild ancestors are rich sources of genetic diversity that have not yet been fully exploited, and it is possible to utilize this diversity. Compared with durum wheat, tetraploid wheat subspecies have retained much greater genetic diversity despite genetic drift and various environmental influences, and the identification and utilization of this diversity can make a greater contribution to the genetic enrichment of wheat. In addition, using the pre-breeding method, the valuable left-behind alleles in the wheat gene pool can be re-introduced through hybridization and introgressive gene flow to create a sustainable opportunity for the genetic gain of wheat. This review provides some insights about the potential of tetraploid wheats in plant breeding and the genetic gains made by them in plant breeding across past decades, and gathers the known functional information on genes/QTLs, metabolites, traits and their direct involvement in wheat resistance/tolerance to biotic/abiotic stresses.

A kinase fusion protein from Aegilops longissima confers resistance to wheat powdery mildew

Many disease resistance genes have been introgressed into wheat from its wild relatives. However, reduced recombination within the introgressed segments hinders the cloning of the introgressed genes. Here, we have cloned the powdery mildew resistance gene Pm13, which is introgressed into wheat from Aegilops longissima, using a method that combines physical mapping with radiation-induced chromosomal aberrations and transcriptome sequencing analysis of ethyl methanesulfonate (EMS)-induced loss-of-function mutants. Pm13 encodes a kinase fusion protein, designated MLKL-K, with an N-terminal domain of mixed lineage kinase domain-like protein (MLKL_NTD domain) and a C-terminal serine/threonine kinase domain bridged by a brace. The resistance function of Pm13 is validated through transient and stable transgenic complementation assays. Transient over-expression analyses in Nicotiana benthamiana leaves and wheat protoplasts reveal that the fragment Brace-Kinase122-476 of MLKL-K is capable of inducing cell death, which is dependent on a functional kinase domain and the three α-helices in the brace region close to the N-terminus of the kinase domain.

Wild grass-derived alleles represent a genetic architecture for the resilience of modern common wheat to stresses

This review explores the integration of wild grass-derived alleles into modern bread wheat breeding to tackle the challenges of climate change and increasing food demand. With a focus on synthetic hexaploid wheat, this review highlights the potential of genetic variability in wheat wild relatives, particularly Aegilops tauschii, for improving resilience to multifactorial stresses like drought, heat, and salinity. The evolutionary journey of wheat (Triticum spp.) from diploid to hexaploid species is examined, revealing significant genetic contributions from wild grasses. We also emphasize the importance of understanding incomplete lineage sorting in the genomic evolution of wheat. Grasping this information is crucial as it can guide breeders in selecting the appropriate alleles from the gene pool of wild relatives to incorporate into modern wheat varieties. This approach improves the precision of phylogenetic relationships and increases the overall effectiveness of breeding strategies. This review also addresses the challenges in utilizing the wheat wild genetic resources, such as the linkage drag and cross-compatibility issues. Finally, we culminate the review with future perspectives, advocating for a combined approach of high-throughput phenotyping tools and advanced genomic techniques to comprehensively understand the genetic and regulatory architectures of wheat under stress conditions, paving the way for more precise and efficient breeding strategies.

Integrative gene duplication and genome-wide analysis as an approach to facilitate wheat reverse genetics: An example in the TaCIPK family

INTRODUCTION

Reverse genetic studies conducted in the plant with a complex or polyploidy genome enriched with large gene families (like wheat) often meet challenges in identifying the key candidate genes related to important traits and prioritizing the genes for functional experiments.

OBJECTIVE

To overcome the above-mentioned challenges of reverse genetics, this work aims to establish an efficient multi-species strategy for genome-wide gene identification and prioritization of the key candidate genes.

METHODS

We established the integrative gene duplication and genome-wide analysis (iGG analysis) as a strategy for pinpointing key candidate genes deserving functional research. The iGG captures the evolution, and the expansion/contraction of large gene families across phylogeny-related species and integrates spatial-temporal expression information for gene function inference. Transgenic approaches were also employed to functional validation.

RESULTS

As a proof-of-concept for the iGG analysis, we took the wheat calcineurin B-like protein-interacting protein kinases (CIPKs) family as an example. We identified CIPKs from seven monocot species, established the orthologous relationship of CIPKs between rice and wheat, and characterized Triticeae-specific CIPK duplicates (e.g., CIPK4 and CIPK17). Integrated with our analysis of CBLs and CBL-CIPK interaction, we revealed that divergent expressions of TaCBLs and TaCIPKs could play an important role in keeping the stoichiometric balance of CBL-CIPK. Furthermore, we validated the function of TaCIPK17-A2 in the regulation of drought tolerance by using transgenic approaches. Overexpression of TaCIPK17 enhanced antioxidant capacity and improved drought tolerance in wheat.

CONCLUSION

The iGG analysis leverages evolutionary and comparative genomics of crops with large genomes to rapidly highlight the duplicated genes potentially associated with speciation, domestication and/or particular traits that deserve reverse-genetic functional studies. Through the identification of Triticeae-specific TaCIPK17 duplicates and functional validation, we demonstrated the effectiveness of the iGG analysis and provided a new target gene for improving drought tolerance in wheat.

Evolutionary genomics of two diploid goat grass species belonging to the section Sitopsis of Aegilops, Aegilops longissima, and Aegilops sharonensis

Aegilops longissima and Ae. sharonensis, being classified into the Sitopsis section of genus Aegilops, are distinct species both taxonomically and ecologically. Nevertheless, earlier observations indicate that the two species are not reproductively isolated to full extent and can inter-bred upon secondary contact. However, the genomic underpinnings of the morpho-ecological differentiation between the two foci species remained unexplored. Here, we resequenced 31 representative accessions of the two species and conducted in-depth comparative genomic analyses. We demonstrate recurrent and ongoing natural hybridizations between Ae. longissima and Ae. sharonensis, and depict features of genome composition of the resultant hybrids at both individual and population levels. We also delineate genomic regions and candidate genes potentially underpinning the differential morphological and edaphic adaptations of the two species. Intriguingly, a binary morphology was observed in the hybrids, suggesting existence of highly diverged genomic regions that remain uneroded by the admixtures. Together, our results provide new insights into the molding effects of interspecific hybridization on genome composition and mechanisms preventing merge of the two species.

Pm57 from Aegilops searsii encodes a tandem kinase protein and confers wheat powdery mildew resistance

Powdery mildew is a devastating disease that affects wheat yield and quality. Wheat wild relatives represent valuable sources of disease resistance genes. Cloning and characterization of these genes will facilitate their incorporation into wheat breeding programs. Here, we report the cloning of Pm57, a wheat powdery mildew resistance gene from Aegilops searsii. It encodes a tandem kinase protein with putative kinase-pseudokinase domains followed by a von Willebrand factor A domain (WTK-vWA), being ortholog of Lr9 that mediates wheat leaf rust resistance. The resistance function of Pm57 is validated via independent mutants, gene silencing, and transgenic assays. Stable Pm57 transgenic wheat lines and introgression lines exhibit high levels of all-stage resistance to diverse isolates of the Bgt fungus, and no negative impacts on agronomic parameters are observed in our experimental set-up. Our findings highlight the emerging role of kinase fusion proteins in plant disease resistance and provide a valuable gene for wheat breeding.

Spontaneous chromosome instability and tissue culture-induced karyotypic alteration in wheat-Thinopyrum intermedium alien addition lines

MAIN CONCLUSION

Wheat lines harboring wild-relative chromosomes can be karyotypically unstable during long-term maintenance. Tissue culture exacerbates chromosomal instability but appears inefficient to induce somatic homoeologous exchange between alien and wheat chromosomes. We assessed if long-term refrigerator storage with regular renewal via self-fertilization, a widely used practice for crop germplasm maintenance, would ensure genetic fidelity of alien addition lines, and explored the possibility of inducing somatic homoeologues exchange by tissue culture. We cytogenetically characterized sampled stock seeds of originally confirmed 12 distinct wheat-Thinopyrum intermedium alien addition lines (dubbed TAI lines), and subjected immature embryos of the TAI lines to tissue culture. We find eight of the 12 TAI lines were karyotypically departed from their original identity as bona fide disomic alien addition lines due to extensive loss of whole-chromosomes of both Th. intermedium and wheat origins during the ca. 3-decade storage. Rampant numerical chromosome variations (NCVs) involving both alien and wheat chromosomes were detected in regenerated plants of all 12 studied TAI lines, but at variable rates among the wheat sub-genomes and chromosomes. Compared with NCVs, structural chromosome variations (SCVs) occurred at substantially lower rates, and no SCV involving the added alien chromosomes was observed. The NCVs manifested only moderate effects on phenotypes of the regenerated plants under field conditions.

Variation in TaSPL6-D confers salinity tolerance in bread wheat by activating TaHKT1;5-D while preserving yield-related traits

Na+ exclusion from above-ground tissues via the Na+-selective transporter HKT1;5 is a major salt-tolerance mechanism in crops. Using the expression genome-wide association study and yeast-one-hybrid screening, we identified TaSPL6-D, a transcriptional suppressor of TaHKT1;5-D in bread wheat. SPL6 also targeted HKT1;5 in rice and Brachypodium. A 47-bp insertion in the first exon of TaSPL6-D resulted in a truncated peptide, TaSPL6-DIn, disrupting TaHKT1;5-D repression exhibited by TaSPL6-DDel. Overexpressing TaSPL6-DDel, but not TaSPL6-DIn, led to inhibited TaHKT1;5-D expression and increased salt sensitivity. Knockout of TaSPL6-DDel in two wheat genotypes enhanced salinity tolerance, which was attenuated by a further TaHKT1;5-D knockdown. Spike development was preserved in Taspl6-dd mutants but not in Taspl6-aabbdd mutants. TaSPL6-DIn was mainly present in landraces, and molecular-assisted introduction of TaSPL6-DIn from a landrace into a leading wheat cultivar successfully improved yield on saline soils. The SPL6-HKT1;5 module offers a target for the molecular breeding of salt-tolerant crops.

Evolution of wheat blast resistance gene Rmg8 accompanied by differentiation of variants recognizing the powdery mildew fungus

Wheat blast, a devastating disease having spread recently from South America to Asia and Africa, is caused by Pyricularia oryzae (synonym of Magnaporthe oryzae) pathotype Triticum, which first emerged in Brazil in 1985. Rmg8 and Rmg7, genes for resistance to wheat blast found in common wheat and tetraploid wheat, respectively, recognize the same avirulence gene, AVR-Rmg8. Here we show that an ancestral resistance gene, which had obtained an ability to recognize AVR-Rmg8 before the differentiation of Triticum and Aegilops, has expanded its target pathogens. Molecular cloning revealed that Rmg7 was an allele of Pm4, a gene for resistance to wheat powdery mildew on 2AL, whereas Rmg8 was its homoeologue on 2BL ineffective against wheat powdery mildew. Rmg8 variants with the ability to recognize AVR-Rmg8 were distributed not only in Triticum spp. but also in Aegilops speltoides, Aegilops umbellulata and Aegilops comosa. This result suggests that the origin of resistance gene(s) recognizing AVR-Rmg8 dates back to the time before differentiation of A, B, S, U and M genomes, that is, ~5 Myr before the emergence of its current target, the wheat blast fungus. Phylogenetic analyses suggested that, in the evolutionary process thereafter, some of their variants gained the ability to recognize the wheat powdery mildew fungus and evolved into genes controlling dual resistance to wheat powdery mildew and wheat blast.

Unveiling the distinctive traits of functional rye centromeres: minisatellites, retrotransposons, and R-loop formation

Centromeres play a vital role in cellular division by facilitating kinetochore assembly and spindle attachments. Despite their conserved functionality, centromeric DNA sequences exhibit rapid evolution, presenting diverse sizes and compositions across species. The functional significance of rye centromeric DNA sequences, particularly in centromere identity, remains unclear. In this study, we comprehensively characterized the sequence composition and organization of rye centromeres. Our findings revealed that these centromeres are primarily composed of long terminal repeat retrotransposons (LTR-RTs) and interspersed minisatellites. We systematically classified LTR-RTs into five categories, highlighting the prevalence of younger CRS1, CRS2, and CRS3 of CRSs (centromeric retrotransposons of Secale cereale) were primarily located in the core centromeres and exhibited a higher association with CENH3 nucleosomes. The minisatellites, mainly derived from retrotransposons, along with CRSs, played a pivotal role in establishing functional centromeres in rye. Additionally, we observed the formation of R-loops at specific regions of CRS1, CRS2, and CRS3, with both rye pericentromeres and centromeres exhibiting enrichment in R-loops. Notably, these R-loops selectively formed at binding regions of the CENH3 nucleosome in rye centromeres, suggesting a potential role in mediating the precise loading of CENH3 to centromeres and contributing to centromere specification. Our work provides insights into the DNA sequence composition, distribution, and potential function of R-loops in rye centromeres. This knowledge contributes valuable information to understanding the genetics and epigenetics of rye centromeres, offering implications for the development of synthetic centromeres in future plant modifications and beyond.

Future of durum wheat research and breeding: Insights from early career researchers

Durum wheat (Triticum turgidum ssp. durum) is globally cultivated for pasta, couscous, and bulgur production. With the changing climate and growing world population, the need to significantly increase durum production to meet the anticipated demand is paramount. This review summarizes recent advancements in durum research, encompassing the exploitation of existing and novel genetic diversity, exploration of potential new diversity sources, breeding for climate-resilient varieties, enhancements in production and management practices, and the utilization of modern technologies in breeding and cultivar development. In comparison to bread wheat (T. aestivum), the durum wheat community and production area are considerably smaller, often comprising many small-family farmers, notably in African and Asian countries. Public breeding programs such as the International Maize and Wheat Improvement Center (CIMMYT) and the International Center for Agricultural Research in the Dry Areas (ICARDA) play a pivotal role in providing new and adapted cultivars for these small-scale growers. We spotlight the contributions of these and others in this review. Additionally, we offer our recommendations on key areas for the durum research community to explore in addressing the challenges posed by climate change while striving to enhance durum production and sustainability. As part of the Wheat Initiative, the Expert Working Group on Durum Wheat Genomics and Breeding recognizes the significance of collaborative efforts in advancing toward a shared objective. We hope the insights presented in this review stimulate future research and deliberations on the trajectory for durum wheat genomics and breeding.

High-quality genomes of Bombax ceiba and Ceiba pentandra provide insights into the evolution of Malvaceae species and differences in their natural fiber development

Members of the Malvaceae family, including Corchorus spp., Gossypium spp., Bombax spp., and Ceiba spp., are important sources of natural fibers. In the past decade, the genomes of several Malvaceae species have been assembled; however, the evolutionary history of Malvaceae species and the differences in their fiber development remain to be clarified. Here, we report the genome assembly and annotation of two natural fiber plants from the Malvaceae, Bombax ceiba and Ceiba pentandra, whose assembled genome sizes are 783.56 Mb and 1575.47 Mb, respectively. Comparative analysis revealed that whole-genome duplication and Gypsy long terminal repeat retroelements have been the major causes of differences in chromosome number (2n = 14 to 2n = 96) and genome size (234 Mb to 2676 Mb) among Malvaceae species. We also used comparative genomic analyses to reconstruct the ancestral Malvaceae karyotype with 11 proto-chromosomes, providing new insights into the evolutionary trajectories of Malvaceae species. MYB-MIXTA-like 3 is relatively conserved among the Malvaceae and functions in fiber cell-fate determination in the epidermis. It appears to perform this function in any tissue where it is expressed, i.e. in fibers on the endocarp of B. ceiba and in ovule fibers of cotton. We identified a structural variation in a cellulose synthase gene and a higher copy number of cellulose synthase-like genes as possible causes of the finer, less spinnable, weaker fibers of B. ceiba. Our study provides two high-quality genomes of natural fiber plants and offers insights into the evolution of Malvaceae species and differences in their natural fiber formation and development through multi-omics analysis.

Mapping and characterization of rust resistance genes Lr53 and Yr35 introgressed from Aegilops species

KEY MESSAGE

The rust resistance genes Lr53 and Yr35 were introgressed into bread wheat from Aegilops longissima or Aegilops sharonensis or their S-genome containing species and mapped to the telomeric region of chromosome arm 6BS. Wheat leaf and stripe rusts are damaging fungal diseases of wheat worldwide. Breeding for resistance is a sustainable approach to control these two foliar diseases. In this study, we used SNP analysis, sequence comparisons, and cytogenetic assays to determine that the chromosomal segment carrying Lr53 and Yr35 was originated from Ae.longissima or Ae. sharonensis or their derived species. In seedling tests, Lr53 conferred strong resistance against all five Chinese Pt races tested, and Yr35 showed effectiveness against Pst race CYR34 but susceptibility to race CYR32. Using a large population (3892 recombinant gametes) derived from plants homozygous for the ph1b mutation obtained from the cross 98M71 × CSph1b, both Lr53 and Yr35 were successfully mapped to a 6.03-Mb telomeric region of chromosome arm 6BS in the Chinese Spring reference genome v1.1. Co-segregation between Lr53 and Yr35 was observed within this large mapping population. Within the candidate region, several nucleotide-binding leucine-rich repeat genes and protein kinases were identified as candidate genes. Marker pku6B3127 was completely linked to both genes and accurately predicted the absence or presence of alien segment harboring Lr53 and Yr35 in 87 tetraploid and 149 hexaploid wheat genotypes tested. We developed a line with a smaller alien segment (< 6.03 Mb) to reduce any potential linkage drag and demonstrated that it conferred resistance levels similar to those of the original donor parent 98M71. The newly developed introgression line and closely linked PCR markers will accelerate the deployment of Lr53 and Yr35 in wheat breeding programs.

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.

Phylogenomic profiles of whole-genome duplications in Poaceae and landscape of differential duplicate retention and losses among major Poaceae lineages

Poaceae members shared a whole-genome duplication called rho. However, little is known about the evolutionary pattern of the rho-derived duplicates among Poaceae lineages and implications in adaptive evolution. Here we present phylogenomic/phylotranscriptomic analyses of 363 grasses covering all 12 subfamilies and report nine previously unknown whole-genome duplications. Furthermore, duplications from a single whole-genome duplication were mapped to multiple nodes on the species phylogeny; a whole-genome duplication was likely shared by woody bamboos with possible gene flow from herbaceous bamboos; and recent paralogues of a tetraploid Oryza are implicated in tolerance of seawater submergence. Moreover, rho duplicates showing differential retention among subfamilies include those with functions in environmental adaptations or morphogenesis, including ACOT for aquatic environments (Oryzoideae), CK2β for cold responses (Pooideae), SPIRAL1 for rapid cell elongation (Bambusoideae), and PAI1 for drought/cold responses (Panicoideae). This study presents a Poaceae whole-genome duplication profile with evidence for multiple evolutionary mechanisms that contribute to gene retention and losses.

Exploring the patterns of evolution: Core thoughts and focus on the saltational model

The Modern Synthesis, a pillar in biological thought, united Darwin's species origin concepts with Mendel's laws of character heredity, providing a comprehensive understanding of evolution within species. Highlighting phenotypic variation and natural selection, it elucidated the environment's role as a selective force, shaping populations over time. This framework integrated additional mechanisms, including genetic drift, random mutations, and gene flow, predicting their cumulative effects on microevolution and the emergence of new species. Beyond the Modern Synthesis, the Extended Evolutionary Synthesis expands perspectives by recognizing the role of developmental plasticity, non-genetic inheritance, and epigenetics. We suggest that these aspects coexist in the plant evolutionary process; in this context, we focus on the saltational model, emphasizing how saltation events, such as dichotomous saltation, chromosomal mutations, epigenetic phenomena, and polyploidy, contribute to rapid evolutionary changes. The saltational model proposes that certain evolutionary changes, such as the rise of new species, may result suddenly from single macromutations rather than from gradual changes in DNA sequences and allele frequencies within a species over time. These events, observed in domesticated and wild higher plants, provide well-defined mechanistic bases, revealing their profound impact on plant diversity and rapid evolutionary events. Notably, next-generation sequencing exposes the likely crucial role of allopolyploidy and autopolyploidy (saltational events) in generating new plant species, each characterized by distinct chromosomal complements. In conclusion, through this review, we offer a thorough exploration of the ongoing dissertation on the saltational model, elucidating its implications for our understanding of plant evolutionary processes and paving the way for continued research in this intriguing field.

Reshaped DNA methylation cooperating with homoeolog-divergent expression promotes improved root traits in synthesized tetraploid wheat

Polyploidization is a major event driving plant evolution and domestication. However, how reshaped epigenetic modifications coordinate gene transcription to generate phenotypic variations during wheat polyploidization is currently elusive. Here, we profiled transcriptomes and DNA methylomes of two diploid wheat accessions (SlSl and AA) and their synthetic allotetraploid wheat line (SlSlAA), which displayed elongated root hair and improved root capability for nitrate uptake and assimilation after tetraploidization. Globally decreased DNA methylation levels with a reduced difference between subgenomes were observed in the roots of SlSlAA. DNA methylation changes in first exon showed strong connections with altered transcription during tetraploidization. Homoeolog-specific transcription was associated with biased DNA methylation as shaped by homoeologous sequence variation. The hypomethylated promoters showed significantly enriched binding sites for MYB, which may affect gene transcription in response to root hair growth. Two master regulators in root hair elongation pathway, AlCPC and TuRSL4, exhibited upregulated transcription levels accompanied by hypomethylation in promoter, which may contribute to the elongated root hair. The upregulated nitrate transporter genes, including NPFs and NRTs, also are significantly associated with hypomethylation, indicating an epigenetic-incorporated regulation manner in improving nitrogen use efficiency. Collectively, these results provided new insights into epigenetic changes in response to crop polyploidization and underscored the importance of epigenetic regulation in improving crop traits.

Wheat powdery mildew resistance gene Pm13 encodes a mixed lineage kinase domain-like protein

Wheat powdery mildew is one of the most destructive diseases threatening global wheat production. The wild relatives of wheat constitute rich sources of diversity for powdery mildew resistance. Here, we report the map-based cloning of the powdery mildew resistance gene Pm13 from the wild wheat species Aegilops longissima. Pm13 encodes a mixed lineage kinase domain-like (MLKL) protein that contains an N-terminal-domain of MLKL (MLKL_NTD) domain in its N-terminus and a C-terminal serine/threonine kinase (STK) domain. The resistance function of Pm13 is validated by mutagenesis, gene silencing, transgenic assay, and allelic association analyses. The development of introgression lines with significantly reduced chromosome segments of Ae. longissima encompassing Pm13 enables widespread deployment of this gene into wheat cultivars. The cloning of Pm13 may provide valuable insights into the molecular mechanisms underlying Pm13-mediated powdery mildew resistance and highlight the important roles of kinase fusion proteins (KFPs) in wheat immunity.

Evolutionary trajectory of organelle-derived nuclear DNAs in the Triticum/Aegilops complex species

Organelle-derived nuclear DNAs, nuclear plastid DNAs (NUPTs), and nuclear mitochondrial DNAs (NUMTs) have been identified in plants. Most, if not all, genes residing in NUPTs/NUMTs (NUPGs/NUMGs) are known to be inactivated and pseudogenized. However, the role of epigenetic control in silencing NUPGs/NUMGs and the dynamic evolution of NUPTs/NUMTs with respect to organismal phylogeny remain barely explored. Based on the available nuclear and organellar genomic resources of wheat (genus Triticum) and goat grass (genus Aegilops) within Triticum/Aegilops complex species, we investigated the evolutionary fates of NUPTs/NUMTs in terms of their epigenetic silencing and their dynamic occurrence rates in the nuclear diploid genomes and allopolyploid subgenomes. NUPTs and NUMTs possessed similar genomic atlas, including (i) predominantly located in intergenic regions and preferential integration to gene regulation regions and (ii) generating sequence variations in the nuclear genome. Unlike nuclear indigenous genes, the alien NUPGs/NUMGs were associated with repressive epigenetic signals, namely high levels of DNA methylation and low levels of active histone modifications. Phylogenomic analyses suggested that the species-specific and gradual accumulation of NUPTs/NUMTs accompanied the speciation processes. Moreover, based on further pan-genomic analyses, we found significant subgenomic asymmetry in the NUPT/NUMT occurrence, which accumulated during allopolyploid wheat evolution. Our findings provide insight into the dynamic evolutionary fates of organelle-derived nuclear DNA in plants.

Dissection of a rapidly evolving wheat resistance gene cluster by long-read genome sequencing accelerated the cloning of Pm69

Gene cloning in repeat-rich polyploid genomes remains challenging. Here, we describe a strategy for overcoming major bottlenecks in cloning of the powdery mildew resistance gene (R-gene) Pm69 derived from tetraploid wild emmer wheat. A conventional positional cloning approach was not effective owing to suppressed recombination. Chromosome sorting was compromised by insufficient purity. A Pm69 physical map, constructed by assembling Oxford Nanopore Technology (ONT) long-read genome sequences, revealed a rapidly evolving nucleotide-binding leucine-rich repeat (NLR) R-gene cluster with structural variations. A single candidate NLR was identified by anchoring RNA sequencing reads from susceptible mutants to ONT contigs and was validated by virus-induced gene silencing. Pm69 is likely a newly evolved NLR and was discovered in only one location across the wild emmer wheat distribution range in Israel. Pm69 was successfully introgressed into cultivated wheat, and a diagnostic molecular marker was used to accelerate its deployment and pyramiding with other R-genes.

Cytological analysis of the diploid-like inheritance of newly synthesized allotetraploid wheat

Polyploidization is a process which is related to species hybridization and whole genome duplication. It is widespread among angiosperm evolution and is essential for speciation and diversification. Allopolyploidization is mainly derived from interspecific hybridization and is believed to pose chromosome imbalances and genome instability caused by meiotic irregularity. However, the self-compatible allopolyploid in wild nature is cytogenetically and genetically stable. Whether this stabilization form was achieved in initial generation or a consequence of long term of evolution was largely unknown. Here, we synthesized a series of nascent allotetraploid wheat derived from three diploid genomes of A, S*, and D. The chromosome numbers of the majority of the progeny derived from these newly formed allotetraploid wheat plants were found to be relatively consistent, with each genome containing 14 chromosomes. In meiosis, bivalent was the majority of the chromosome configuration in metaphase I which supports the stable chromosome number inheritance in the nascent allotetraploid. These findings suggest that diploidization occurred in the newly formed synthetic allotetraploid wheat. However, we still detected aneuploids in a proportion of newly formed allotetraploid wheat, and meiosis of these materials present more irregular chromosome behavior than the euploid. We found that centromere pairing and centromere clustering in meiosis was affected in the aneuploids, which suggest that aneuploidy may trigger the irregular interactions of centromere in early meiosis which may take participate in promoting meiosis stabilization in newly formed allotetraploid wheat.

Effects of Aegilops longissima chromosome 1Sl on wheat bread-making quality in two types of translocation lines

KEY MESSAGE

Two wheat-Ae. longissima translocation chromosomes (1BS·1SlL and 1SlS·1BL) were transferred into three commercial wheat varieties, and the new advanced lines showed improved bread-making quality compared to their recurrent parents. Aegilops longissima chromosome 1Sl encodes specific types of gluten subunits that may positively affect wheat bread-making quality. The most effective method of introducing 1Sl chromosomal fragments containing the target genes into wheat is chromosome translocation. Here, a wheat-Ae. longissima 1BS·1SlL translocation line was developed using molecular marker-assisted chromosome engineering. Two types of translocation chromosomes developed in a previous study, 1BS·1SlL and 1SlS·1BL, were introduced into three commercial wheat varieties (Ningchun4, Ningchun50, and Westonia) via backcrossing with marker-assisted selection. Advanced translocation lines were confirmed through chromosome in situ hybridization and genotyping by target sequencing using the wheat 40 K system. Bread-making quality was found to be improved in the two types of advanced translocation lines compared to the corresponding recurrent parents. Furthermore, 1SlS·1BL translocation lines displayed better bread-making quality than 1BS·1SlL translocation lines in each genetic background. Further analysis revealed that high molecular weight glutenin subunit (HMW-GS) contents and expression levels of genes encoding low molecular weight glutenin subunits (LMW-GSs) were increased in 1SlS·1BL translocation lines. Gliadin and gluten-related transcription factors were also upregulated in the grains of the two types of advanced translocation lines compared to the recurrent parents. This study clarifies the impacts of specific glutenin subunits on bread-making quality and provides novel germplasm resources for further improvement of wheat quality through molecular breeding.

Diploid species phylogeny and evolutionary reticulation indicate early radiation of Ephedra in the Tethys coast

Reconstructing a robust species phylogeny and disentangling the evolutionary and biogeographic history of the gymnosperm genus Ephedra, which has a large genome and rich polyploids, remain a big challenge. Here we reconstructed a transcriptome-based phylogeny of 19 diploid Ephedra species, and explored evolutionary reticulations in this genus represented by 50 diploid and polyploid species, using four low-copy nuclear and nine plastid genes. The diploid species phylogeny indicates that the Mediterranean species diverged first, and the remaining species split into three clades, including the American species (Clade A), E. rhytidosperma, and all other Asian species (Clade B). The single-gene trees placed E. rhytidosperma sister to Clade A, Clade B, or Clades A + B in similar proportions, suggesting that radiation and gene flow likely occurred in the early evolution of Ephedra. In addition, reticulate evolution occurred not only among the deep nodes, but also in the recently evolved South American species, which further caused difficulty in phylogenetic reconstruction. Moreover, we found that allopolyploid speciation was pervasive in Ephedra. Our study also suggests that Ephedra very likely originated in the Tethys coast during the late Cretaceous, and the South American Ephedra species have a single origin by dispersal from Mexico or North America.

Subgenome phasing for complex allopolyploidy: case-based benchmarking and recommendations

Accurate subgenome phasing is crucial for understanding the origin, evolution and adaptive potential of polyploid genomes. SubPhaser and WGDI software are two common methodologies for subgenome phasing in allopolyploids, particularly in scenarios lacking known diploid progenitors. Triggered by a recent debate over the subgenomic origins of the cultivated octoploid strawberry, we examined four well-documented complex allopolyploidy cases as benchmarks, to evaluate and compare the accuracy of the two software. Our analysis demonstrates that the subgenomic structure phased by both software is in line with prior research, effectively tracing complex allopolyploid evolutionary trajectories despite the limitations of each software. Furthermore, using these validated methodologies, we revisited the controversial issue regarding the progenitors of the octoploid strawberry. The results of both methodologies reaffirm Fragaria vesca and Fragaria iinumae as progenitors of the octoploid strawberry. Finally, we propose recommendations for enhancing the accuracy of subgenome phasing in future studies, recognizing the potential of integrated tools for advanced complex allopolyploidy research and offering a new roadmap for robust subgenome-based phylogenetic analysis.

Pangenome of water caltrop reveals structural variations and asymmetric subgenome divergence after allopolyploidization

Water caltrop (Trapa spp., Lythraceae) is a traditional but currently underutilized non-cereal crop. Here, we generated chromosome-level genome assemblies for the two diploid progenitors of allotetraploid Trapa. natans (4x, AABB), i.e., diploid T. natans (2x, AA) and Trapa incisa (2x, BB). In conjunction with four published (sub)genomes of Trapa, we used gene-based and graph-based pangenomic approaches and a pangenomic transposable element (TE) library to develop Trapa genomic resources. The pangenome displayed substantial gene-content variation with dispensable and private gene clusters occupying a large proportion (51.95%) of the total cluster sets in the six (sub)genomes. Genotyping of presence-absence variation (PAVs) identified 40 453 PAVs associated with 2570 genes specific to A- or B-lineages, of which 1428 were differentially expressed, and were enriched in organ development process, organic substance metabolic process and response to stimulus. Comparative genome analyses showed that the allotetraploid T. natans underwent asymmetric subgenome divergence, with the B-subgenome being more dominant than the A-subgenome. Multiple factors, including PAVs, asymmetrical amplification of TEs, homeologous exchanges (HEs), and homeolog expression divergence, together affected genome evolution after polyploidization. Overall, this study sheds lights on the genome architecture and evolution of Trapa, and facilitates its functional genomic studies and breeding program.

Deciphering the evolution and complexity of wheat germplasm from a genomic perspective

Bread wheat provides an essential fraction of the daily calorific intake for humanity. Due to its huge and complex genome, progress in studying on the wheat genome is substantially trailed behind those of the other two major crops, rice and maize, for at least a decade. With rapid advances in genome assembling and reduced cost of high-throughput sequencing, emerging de novo genome assemblies of wheat and whole-genome sequencing data are leading to a paradigm shift in wheat research. Here, we review recent progress in dissecting the complex genome and germplasm evolution of wheat since the release of the first high-quality wheat genome. New insights have been gained in the evolution of wheat germplasm during domestication and modern breeding progress, genomic variations at multiple scales contributing to the diversity of wheat germplasm, and complex transcriptional and epigenetic regulations of functional genes in polyploid wheat. Genomics databases and bioinformatics tools meeting the urgent needs of wheat genomics research are also summarized. The ever-increasing omics data, along with advanced tools and well-structured databases, are expected to accelerate deciphering the germplasm and gene resources in wheat for future breeding advances.

Discovery of isoflavone phytoalexins in wheat reveals an alternative route to isoflavonoid biosynthesis

Isoflavones are a group of phenolic compounds mostly restricted to plants of the legume family, where they mediate important interactions with plant-associated microbes, including in defense from pathogens and in nodulation. Their well-studied health promoting attributes have made them a prime target for metabolic engineering, both for bioproduction of isoflavones as high-value molecules, and in biofortification of food crops. A key gene in their biosynthesis, isoflavone synthase, was identified in legumes over two decades ago, but little is known about formation of isoflavones outside of this family. Here we identify a specialized wheat-specific isoflavone synthase, TaCYP71F53, which catalyzes a different reaction from the leguminous isoflavone synthases, thus revealing an alternative path to isoflavonoid biosynthesis and providing a non-transgenic route for engineering isoflavone production in wheat. TaCYP71F53 forms part of a biosynthetic gene cluster that produces a naringenin-derived O-methylated isoflavone, 5-hydroxy-2',4',7-trimethoxyisoflavone, triticein. Pathogen-induced production and in vitro antimicrobial activity of triticein suggest a defense-related role for this molecule in wheat. Genomic and metabolic analyses of wheat ancestral grasses further show that the triticein gene cluster was introduced into domesticated emmer wheat through natural hybridization ~9000 years ago, and encodes a pathogen-responsive metabolic pathway that is conserved in modern bread wheat varieties.

Genomic characterization and gene bank curation of Aegilops: the wild relatives of wheat

Genetic diversity found in crop wild relatives is critical to preserve and utilize for crop improvement to achieve sustainable food production amid climate change and increased demand. We genetically characterized a large collection of 1,041 Aegilops accessions distributed among 23 different species using more than 45K single nucleotide polymorphisms identified by genotyping-by-sequencing. The Wheat Genetics Resource Center (WGRC) Aegilops germplasm collection was curated through the identification of misclassified and redundant accessions. There were 49 misclassified and 28 sets of redundant accessions within the four diploid species. The curated germplasm sets now have improved utility for genetic studies and wheat improvement. We constructed a phylogenetic tree and principal component analysis cluster for all Aegilops species together, giving one of the most comprehensive views of Aegilops. The Sitopsis section and the U genome Aegilops clade were further scrutinized with in-depth population analysis. The genetic relatedness among the pair of Aegilops species provided strong evidence for the species evolution, speciation, and diversification. We inferred genome symbols for two species Ae. neglecta and Ae. columnaris based on the sequence read mapping and the presence of segregating loci on the pertinent genomes as well as genetic clustering. The high genetic diversity observed among Aegilops species indicated that the genus could play an even greater role in providing the critical need for untapped genetic diversity for future wheat breeding and improvement. To fully characterize these Aegilops species, there is an urgent need to generate reference assemblies for these wild wheats, especially for the polyploid Aegilops.

Crop adaptation to climate change: An evolutionary perspective

The disciplines of evolutionary biology and plant and animal breeding have been intertwined throughout their development, with responses to artificial selection yielding insights into the action of natural selection and evolutionary biology providing statistical and conceptual guidance for modern breeding. Here we offer an evolutionary perspective on a grand challenge of the 21st century: feeding humanity in the face of climate change. We first highlight promising strategies currently under way to adapt crops to current and future climate change. These include methods to match crop varieties with current and predicted environments and to optimize breeding goals, management practices, and crop microbiomes to enhance yield and sustainable production. We also describe the promise of crop wild relatives and recent technological innovations such as speed breeding, genomic selection, and genome editing for improving environmental resilience of existing crop varieties or for developing new crops. Next, we discuss how methods and theory from evolutionary biology can enhance these existing strategies and suggest novel approaches. We focus initially on methods for reconstructing the evolutionary history of crops and their pests and symbionts, because such historical information provides an overall framework for crop-improvement efforts. We then describe how evolutionary approaches can be used to detect and mitigate the accumulation of deleterious mutations in crop genomes, identify alleles and mutations that underlie adaptation (and maladaptation) to agricultural environments, mitigate evolutionary trade-offs, and improve critical proteins. Continuing feedback between the evolution and crop biology communities will ensure optimal design of strategies for adapting crops to climate change.

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.

New insights on the evolution of nucleolar dominance in newly resynthesized hexaploid wheat Triticum zhukovskyi

Nucleolar dominance (ND) is a widespread epigenetic phenomenon in hybridizations where nucleolus transcription fails at the nucleolus organizer region (NOR). However, the dynamics of NORs during the formation of Triticum zhukovskyi (GGAu Au Am Am ), another evolutionary branch of allohexaploid wheat, remains poorly understood. Here, we elucidated genetic and epigenetic changes occurring at the NOR loci within the Am , G, and D subgenomes during allopolyploidization by synthesizing hexaploid wheat GGAu Au Am Am and GGAu Au DD. In T. zhukovskyi, Au genome NORs from T. timopheevii (GGAu Au ) were lost, while the second incoming NORs from T. monococcum (Am Am ) were retained. Analysis of the synthesized T. zhukovskyi revealed that rRNA genes from the Am genome were silenced in F1 hybrids (GAu Am ) and remained inactive after genome doubling and subsequent self-pollinations. We observed increased DNA methylation accompanying the inactivation of NORs in the Am genome and found that silencing of NORs in the S1 generation could be reversed by a cytidine methylase inhibitor. Our findings provide insights into the ND process during the evolutionary period of T. zhukovskyi and highlight that inactive rDNA units may serve as a 'first reserve' in the form of R-loops, contributing to the successful evolution of T. zhukovskyi.

All families of transposable elements were active in the recent wheat genome evolution and polyploidy had no impact on their activity

Bread wheat (Triticum aestivum L.) is a major crop and its genome is one of the largest ever assembled at reference-quality level. It is 15 Gb, hexaploid, with 85% of transposable elements (TEs). Wheat genetic diversity was mainly focused on genes and little is known about the extent of genomic variability affecting TEs, transposition rate, and the impact of polyploidy. Multiple chromosome-scale assemblies are now available for bread wheat and for its tetraploid and diploid wild relatives. In this study, we computed base pair-resolved, gene-anchored, whole genome alignments of A, B, and D lineages at different ploidy levels in order to estimate the variability that affects the TE space. We used assembled genomes of 13 T. aestivum cultivars (6x = AABBDD) and a single genome for Triticum durum (4x = AABB), Triticum dicoccoides (4x = AABB), Triticum urartu (2x = AA), and Aegilops tauschii (2x = DD). We show that 5%-34% of the TE fraction is variable, depending on the species divergence. Between 400 and 13,000 novel TE insertions per subgenome were detected. We found lineage-specific insertions for nearly all TE families in di-, tetra-, and hexaploids. No burst of transposition was observed and polyploidization did not trigger any boost of transposition. This study challenges the prevailing idea of wheat TE dynamics and is more in agreement with an equilibrium model of evolution.

Pangenome-based trajectories of intracellular gene transfers in Poaceae unveil high cumulation in Triticeae

Intracellular gene transfers (IGTs) between the nucleus and organelles, including plastids and mitochondria, constantly reshape the nuclear genome during evolution. Despite the substantial contribution of IGTs to genome variation, the dynamic trajectories of IGTs at the pangenomic level remain elusive. Here, we developed an approach, IGTminer, that maps the evolutionary trajectories of IGTs using collinearity and gene reannotation across multiple genome assemblies. We applied IGTminer to create a nuclear organellar gene (NOG) map across 67 genomes covering 15 Poaceae species, including important crops. The resulting NOGs were verified by experiments and sequencing data sets. Our analysis revealed that most NOGs were recently transferred and lineage specific and that Triticeae species tended to have more NOGs than other Poaceae species. Wheat (Triticum aestivum) had a higher retention rate of NOGs than maize (Zea mays) and rice (Oryza sativa), and the retained NOGs were likely involved in photosynthesis and translation pathways. Large numbers of NOG clusters were aggregated in hexaploid wheat during 2 rounds of polyploidization, contributing to the genetic diversity among modern wheat accessions. We implemented an interactive web server to facilitate the exploration of NOGs in Poaceae. In summary, this study provides resources and insights into the roles of IGTs in shaping interspecies and intraspecies genome variation and driving plant genome evolution.

Comparative Analysis Reveals Different Evolutionary Fates and Biological Functions in Wheat Duplicated Genes (Triticum aestivum L.)

Wheat (Triticum aestivum L.) is a staple food crop that provides 20% of total human calorie consumption. Gene duplication has been considered to play an important role in evolution by providing new genetic resources. However, the evolutionary fates and biological functions of the duplicated genes in wheat remain to be elucidated. In this study, the resulting data showed that the duplicated genes evolved faster with shorter gene lengths, higher codon usage bias, lower expression levels, and higher tissue specificity when compared to non-duplicated genes. Our analysis further revealed functions of duplicated genes in various biological processes with significant enrichment to environmental stresses. In addition, duplicated genes derived from dispersed, proximal, tandem, transposed, and whole-genome duplication differed in abundance, evolutionary rate, gene compactness, expression pattern, and genetic diversity. Tandem and proximal duplicates experienced stronger selective pressure and showed a more compact gene structure with diverse expression profiles than other duplication modes. Moreover, genes derived from different duplication modes showed an asymmetrical evolutionary pattern for wheat A, B, and D subgenomes. Several candidate duplication hotspots associated with wheat domestication or polyploidization were characterized as potential targets for wheat molecular breeding. Our comprehensive analysis revealed the evolutionary trajectory of duplicated genes and laid the foundation for future functional studies on wheat.

Fighting wheat powdery mildew: from genes to fields

KEY MESSAGE

Host resistance conferred by Pm genes provides an effective strategy to control powdery mildew. The study of Pm genes helps modern breeding develop toward more intelligent and customized. Powdery mildew of wheat is one of the most destructive diseases seriously threatening the crop yield and quality worldwide. The genetic research on powdery mildew (Pm) resistance has entered a new era. Many Pm genes from wheat and its wild and domesticated relatives have been mined and cloned. Meanwhile, modern breeding strategies based on high-throughput sequencing and genome editing are emerging and developing toward more intelligent and customized. This review highlights mining and cloning of Pm genes, molecular mechanism studies on the resistance and avirulence genes, and prospects for genomic-assisted breeding for powdery mildew resistance in wheat.

Identification of acyl-CoA-binding protein gene in Triticeae species reveals that TaACBP4A-1 and TaACBP4A-2 positively regulate powdery mildew resistance in wheat

Plant acyl-CoA-binding proteins (ACBPs), which contain the conserved ACB domain, participate in multiple biological processes, however, there are few reports on wheat ACBPs. In this study, the ACBP genes from nine different species were identified comprehensively. The expression patterns of TaACBP genes in multiple tissues and under various biotic stresses were determined by qRT-PCR. The function of selected TaACBP genes was studied by virus-induced gene silencing. A total of 67 ACBPs were identified from five monocotyledonous and four dicotyledonous species and divided into four classes. Tandem duplication analysis of the ACBPs suggested that tandem duplication events occurred in Triticum dicoccoides, but there was no tandem duplication event in wheat ACBP genes. Evolutionary analysis suggested that the TdACBPs may have experienced gene introgression during tetraploid evolution, while TaACBP gene loss events occurred during hexaploid wheat evolution. The expression pattern showed that all the TaACBP genes were expressed, and most of them were responsive to induction by Blumeria graminis f. sp. tritici or Fusarium graminearum. Silencing of TaACBP4A-1 and TaACBP4A-2 increased powdery mildew susceptibility in the common wheat BainongAK58. Furthermore, TaACBP4A-1, which belonged to class III, physically interacted with autophagy-related ubiquitin-like protein TaATG8g in yeast cells. This study provided a valuable reference for further investigations into the functional and molecular mechanisms of the ACBP gene family.

Genome-wide association mapping of rust resistance in Aegilops longissima

The rust diseases, including leaf rust caused by Puccinia triticina (Pt), stem rust caused by P. graminis f. sp. tritici (Pgt), and stripe rust caused by P. striiformis f. sp. tritici (Pst), are major limiting factors in wheat production worldwide. Identification of novel sources of rust resistance genes is key to developing cultivars resistant to rapidly evolving pathogen populations. Aegilops longissima is a diploid wild grass native to the Levant and closely related to the modern bread wheat D subgenome. To explore resistance genes in the species, we evaluated a large panel of Ae. longissima for resistance to several races of Pt, Pgt, and Pst, and conducted a genome-wide association study (GWAS) to map rust resistance loci in the species. A panel of 404 Ae. longissima accessions, mostly collected from Israel, were screened for seedling-stage resistance to four races of Pt, four races of Pgt, and three races of Pst. Out of the 404 accessions screened, two were found that were resistant to all 11 races of the three rust pathogens screened. The percentage of all accessions screened that were resistant to a given rust pathogen race ranged from 18.5% to 99.7%. Genotyping-by-sequencing (GBS) was performed on 381 accessions of the Ae. longissima panel, wherein 125,343 single nucleotide polymorphisms (SNPs) were obtained after alignment to the Ae. longissima reference genome assembly and quality control filtering. Genetic diversity analysis revealed the presence of two distinct subpopulations, which followed a geographic pattern of a northern and a southern subpopulation. Association mapping was performed in the genotyped portion of the collection (n = 381) and in each subpopulation (n = 204 and 174) independently via a single-locus mixed-linear model, and two multi-locus models, FarmCPU, and BLINK. A large number (195) of markers were significantly associated with resistance to at least one of 10 rust pathogen races evaluated, nine of which are key candidate markers for further investigation due to their detection via multiple models and/or their association with resistance to more than one pathogen race. The novel resistance loci identified will provide additional diversity available for use in wheat breeding.

Sequencing trait-associated mutations to clone wheat rust-resistance gene YrNAM

Stripe (yellow) rust, caused by Puccinia striiformis f. sp. tritici (Pst), can significantly affect wheat production. Cloning resistance genes is critical for efficient and effective breeding of stripe rust resistant wheat cultivars. One resistance gene (Yr10CG) underlying the Pst resistance locus Yr10 has been cloned. However, following haplotype and linkage analyses indicate the presence of additional Pst resistance gene(s) underlying/near Yr10 locus. Here, we report the cloning of the Pst resistance gene YrNAM in this region using the method of sequencing trait-associated mutations (STAM). YrNAM encodes a non-canonical resistance protein with a NAM domain and a ZnF-BED domain. We show that both domains are required for resistance. Transgenic wheat harboring YrNAM gene driven by its endogenous promoter confers resistance to stripe rust races CYR32 and CYR33. YrNAM is an ancient gene and present in wild wheat species Aegilops longissima and Ae. sharonensis; however, it is absent in most wheat cultivars, which indicates its breeding value.

Genome-wide characterization of the VQ genes in Triticeae and their functionalization driven by polyploidization and gene duplication events in wheat

Valine-glutamine motif-containing (VQ) proteins are transcriptional cofactors widely involved in plant growth, development, and response to various stresses. Although the VQ family has been genome-wide identified in some species, but the knowledge regarding duplication-driven functionalization of VQ genes among evolutionarily related species is still lacking. Here, 952 VQ genes have been identified from 16 species, emphasizing seven Triticeae species including the bread wheat. Comprehensive phylogenetic and syntenic analyses allow us to establish the orthologous relationship of VQ genes from rice (Oryza sativa) to bread wheat (Triticum aestivum). The evolutionary analysis revealed that whole-genome duplication (WGD) drives the expansion of OsVQs, while TaVQs expansion is associated with a recent burst of gene duplication (RBGD). We also analyzed the motif composition and molecular properties of TaVQ proteins, enriched biological functions, and expression patterns of TaVQs. We demonstrate that WGD-derived TaVQs have become divergent in both protein motif composition and expression pattern, while RBGD-derived TaVQs tend to adopt specific expression patterns, suggesting their functionalization in certain biological processes or in response to specific stresses. Furthermore, some RBGD-derived TaVQs are found to be associated with salt tolerance. Several of the identified salt-related TaVQ proteins were located in the cytoplasm and nucleus and their salt-responsive expression patterns were validated by qPCR analysis. Yeast-based functional experiments confirmed that TaVQ27 may be a new regulator to salt response and regulation. Overall, this study lays the foundation for further functional validation of VQ family members within the Triticeae species.

New insights into polyploid evolution and dynamic nature of Ludwigia section Isnardia (Onagraceae)

BACKGROUND

While polyploids are common in plants, the evolutionary history and natural dynamics of most polyploid groups are still unclear. Owing to plentiful earlier systematic studies, Ludwigia sect. Isnardia (comprising 22 wetland taxa) is an ideal allopolyploid complex to investigate polyploid evolution and natural dynamics within and among taxa. With a considerable sampling, we concentrated on revisiting earlier phylogenies of Isnardia, reevaluating the earlier estimated age of the most recent common ancestor (TMRCA), exploring the correlation between infraspecific genetic diversity and ploidy levels, and inspecting interspecific gene flows among taxa.

RESULTS

Phylogenetic trees and network concurred with earlier phylogenies and hypothesized genomes by incorporating 192 atpB-rbcL and ITS sequences representing 91% of Isnardia taxa. Moreover, we detected three multi-origin taxa. Our findings on L. repens and L. sphaerocarpa were consistent with earlier studies; L. arcuata was reported as a multi-origin taxon here, and an additional evolutionary scenario of L. sphaerocarpa was uncovered, both for the first time. Furthermore, estimated Isnardia TMRCA ages based on our data (5.9 or 8.9 million years ago) are in accordance with earlier estimates, although younger than fossil dates (Middle Miocene). Surprisingly, infraspecific genetic variations of Isnardia taxa did not increase with ploidy levels as anticipated from many other polyploid groups. In addition, the exuberant, low, and asymmetrical gene flows among Isnardia taxa indicated that the reproductive barriers may be weakened owing to allopolyploidization, which has rarely been reported.

CONCLUSIONS

The present research gives new perceptions of the reticulate evolution and dynamic nature of Isnardia and points to gaps in current knowledge about allopolyploid evolution.

An unusual tandem kinase fusion protein confers leaf rust resistance in wheat

The introgression of chromosome segments from wild relatives is an established strategy to enrich crop germplasm with disease-resistance genes¹. Here we use mutagenesis and transcriptome sequencing to clone the leaf rust resistance gene Lr9, which was introduced into bread wheat from the wild grass species Aegilops umbellulata². We established that Lr9 encodes an unusual tandem kinase fusion protein. Long-read sequencing of a wheat Lr9 introgression line and the putative Ae. umbellulata Lr9 donor enabled us to assemble the ~28.4-Mb Lr9 translocation and to identify the translocation breakpoint. We likewise cloned Lr58, which was reportedly introgressed from Aegilops triuncialis³, but has an identical coding sequence compared to Lr9. Cytogenetic and haplotype analyses corroborate that the two genes originate from the same translocation event. Our work sheds light on the emerging role of kinase fusion proteins in wheat disease resistance, expanding the repertoire of disease-resistance genes for breeding.

Differential alternative polyadenylation of homoeologous genes of allohexaploid wheat ABD subgenomes during drought stress response

Because allohexaploid wheat genome contains ABD subgenomes, how the expression of homoeologous genes is coordinated remains largely unknown, particularly at the co-transcriptional level. Alternative polyadenylation (APA) is an important part of co-transcriptional regulation, which is crucial in developmental processes and stress responses. Drought stress is a major threat to the stable yield of wheat. Focusing on APA, we used poly(A) tag sequencing to track poly(A) site dynamics in wheat under drought stress. The results showed that drought stress led to extensive APA involving 37-47% of differentially expressed genes in wheat. Significant poly(A) site switching was found in stress-responsive genes. Interestingly, homoeologous genes exhibit unequal numbers of poly(A) sites, divergent APA patterns with tissue specificity and time-course dynamics, and distinct 3'-UTR length changes. Moreover, differentially expressed transcripts in leaves and roots used different poly(A) signals, the up- and downregulated isoforms had distinct preferences for non-canonical poly(A) sites. Genes that encode key polyadenylation factors showed differential expression patterns under drought stress. In summary, poly(A) signals and the changes in core poly(A) factors may widely affect the selection of poly(A) sites and gene expression levels during the response to drought stress, and divergent APA patterns among homoeologous genes add extensive plasticity to this responsive network. These results not only reveal the significant role of APA in drought stress response, but also provide a fresh perspective on how homoeologous genes contribute to adaptability through transcriptome diversity. In addition, this work provides information about the ends of transcripts for a better annotation of the wheat genome.

The Influence of Edaphic Factors on DNA Damage and Repair in Wild Wheat Triticum dicoccoides Körn. (Poaceae, Triticeae)

A complex DNA repair network maintains genome integrity and genetic stability. In this study, the influence of edaphic factors on DNA damage and repair in wild wheat Triticum dicoccoides was addressed. Plants inhabiting two abutting microsites with dry terra rossa and humid basalt soils were studied. The relative expression level of seven genes involved in DNA repair pathways-RAD51, BRCA1, LigIV, KU70, MLH1, MSH2, and MRE11-was assessed using quantitative real-time PCR (qPCR). Immunolocalization of RAD51, LigIV, γH2AX, RNA Polymerase II, and DNA-RNA hybrid [S9.6] (R-loops) in somatic interphase nuclei and metaphase chromosomes was carried out in parallel. The results showed a lower expression level of genes involved in DNA repair and a higher number of DNA double-strand breaks (DSBs) in interphase nuclei in plants growing in terra rossa soil compared with plants in basalt soil. Further, the number of DSBs and R-loops in metaphase chromosomes was also greater in plants growing on terra rossa soil. Finally, RAD51 and LigIV foci on chromosomes indicate ongoing DSB repair during the M-phase via the Homologous Recombination and Non-Homologous End Joining pathways. Together, these results show the impact of edaphic factors on DNA damage and repair in the wheat genome adapted to contrasting environments.

CRISPR/Cas9-targeted mutagenesis of TaDCL4, TaDCL5 and TaRDR6 induces male sterility in common wheat

Phased, small interfering RNAs (phasiRNAs) are important for plant anther development, especially for male sterility. PhasiRNA biogenesis is dependent on genes like RNA polymerase 6 (RDR6), DICER-LIKE 4 (DCL4), or DCL5 to produce 21- or 24 nucleotide (nt) double-strand small RNAs. Here, we generated mutants of DCL4, DCL5 and RDR6 using CRISPR/Cas9 system and studied their effects on plant reproductive development and phasiRNA production in wheat. We found that RDR6 mutation caused sever consequence throughout plant development starting from seed germination and the dcl4 mutants grew weaker with thorough male sterility, while dcl5 plants developed normally but exhibited male sterility. Correspondingly, DCL4 and DCL5, respectively, specified 21- and 24-nt phasiRNA biogenesis, while RDR6 contributed to both. Also, the three key genes evolved differently in wheat, with TaDCL5-A/B becoming non-functioning and TaRDR6-A being lost after polyploidization. Furthermore, we found that PHAS genes (phasiRNA precursors) identified via phasiRNAs diverged rapidly among sub-genomes of polyploid wheat. Despite no similarity being found among phasiRNAs of grasses, their targets were enriched for similar biological functions. In light of the important roles of phasiRNA pathways in gametophyte development, genetic dissection of the function of key genes may help generate male sterile lines suitable for hybrid wheat breeding.

Orthologous genes Pm12 and Pm21 from two wild relatives of wheat show evolutionary conservation but divergent powdery mildew resistance

Wheat powdery mildew, caused by Blumeria graminis f. sp. tritici (Bgt), is a devastating disease that threatens wheat production worldwide. Pm12, which originated from Aegilops speltoides, a wild relative of wheat, confers strong resistance to powdery mildew and therefore has potential use in wheat breeding. Using susceptible mutants induced by gamma irradiation, we physically mapped and isolated Pm12 and showed it to be orthologous to Pm21 from Dasypyrum villosum, also a wild relative of wheat. The resistance function of Pm12 was validated via ethyl methanesulfonate mutagenesis, virus-induced gene silencing, and stable genetic transformation. Evolutionary analysis indicates that the Pm12/Pm21 loci in wheat species are relatively conserved but dynamic. Here, we demonstrated that the two orthologous genes, Pm12 and Pm21, possess differential resistance against the same set of Bgt isolates. Overexpression of the coiled-coil domains of both PM12 and PM21 induces cell death in Nicotiana benthamiana leaves. However, their full-length forms display different cell death-inducing activities caused by their distinct intramolecular interactions. Cloning of Pm12 will facilitate its application in wheat breeding programs. This study also gives new insight into two orthologous resistance genes, Pm12 and Pm21, which show different race specificities and intramolecular interaction patterns.

Genome sequencing of Sitopsis species provides insights into their contribution to the B subgenome of bread wheat

Wheat (Triticum aestivum, BBAADD) is an allohexaploid species that originated from two polyploidization events. The progenitors of the A and D subgenomes have been identified as Triticum urartu and Aegilops tauschii, respectively. Current research suggests that Aegilops speltoides is the closest but not the direct ancestor of the B subgenome. However, whether Ae. speltoides has contributed genomically to the wheat B subgenome and which chromosome regions are conserved between Ae. speltoides and the B subgenome remain unclear. Here, we assembled a high-quality reference genome for Ae. speltoides, resequenced 53 accessions from seven species (Aegilops bicornis, Aegilops longissima, Aegilops searsii, Aegilops sharonensis, Ae. speltoides, Aegilops mutica [syn. Amblyopyrum muticum], and Triticum dicoccoides) and revealed their genomic contributions to the wheat B subgenome. Our results showed that centromeric regions were particularly conserved between Aegilops and Triticum and revealed 0.17 Gb of conserved blocks between Ae. speltoides and the B subgenome. We classified five groups of conserved and non-conserved genes between Aegilops and Triticum, revealing their biological characteristics, differentiation in gene expression patterns, and collinear relationships between Ae. speltoides and the wheat B subgenome. We also identified gene families that expanded in Ae. speltoides during its evolution and 789 genes specific to Ae. speltoides. These genes can serve as genetic resources for improvement of adaptability to biotic and abiotic stress. The newly constructed reference genome and large-scale resequencing data for Sitopsis species will provide a valuable genomic resource for wheat genetic improvement and genomic studies.

A chromosome-scale genome assembly of Dasypyrum villosum provides insights into its application as a broad-spectrum disease resistance resource for wheat improvement

Dasypyrum villosum is one of the most valuable gene resources in wheat improvement, especially for disease resistance. The mining of favorable genes from D. villosum is frustrated by the lack of a whole genome sequence. In this study, we generated a doubled-haploid line, 91C43^DH, using microspore culture and obtained a 4.05-GB high-quality, chromosome-scale genome assembly for D. villosum. The assembly contains39 727 high-confidence genes, and 85.31% of the sequences are repetitive. Two reciprocal translocation events were detected, and 7VS-4VL is a unique translocation in D. villosum. The prolamin seed storage protein-coding genes were found to be duplicated; in particular, the genes encoding low-molecular-weight glutenin at the Glu-V3 locus were significantly expanded. RNA sequencing (RNA-seq) analysis indicated that, after Blumeria graminearum f.sp tritici (Bgt) inoculation, there were more upregulated genes involved in the pattern-triggered immunity and effector-triggered immunity defense pathways in D. villosum than in Triticum urartu. MNase hypersensitive sequencing (MH-seq) identified two Bgt-inducible MH sites (MHSs), one in the promoter and one in the 3' terminal region of the powdery mildew resistance (Pm) gene NLR1-V. Each site had two subpeaks and they were termed MHS1 (MHS1.1/1.2) and MHS2 (MHS2.1/2.2). Bgt-inducible MHS2.2 was uniquely present in D. villosum, and MHS1.1 was more inducible in D. villosum than in wheat, suggesting that MHSs may be critical for regulation of NLR1-V expression and plant defense. In summary, this study provides a valuable genome resource for functional genomics studies and wheat-D. villosum introgression breeding. The identified regulatory mechanisms may also be exploited to develop new strategies for enhancing Pm resistance by optimizing gene expression in wheat.