Multiple wheat genomes reveal global variation in modern breeding

Expression interplay of genes coding for calcium-binding proteins and transcription factors during the osmotic phase provides insights on salt stress response mechanisms in bread wheat

Bread wheat is an important crop for the human diet, but the increasing soil salinization is reducing the yield. The Ca2+ signaling events at the early stages of the osmotic phase of salt stress are crucial for the acclimation response of the plants through the performance of calcium-sensing proteins, which activate or repress transcription factors (TFs) that affect the expression of downstream genes. Physiological, genetic mapping, and transcriptomics studies performed with the contrasting genotypes Syn86 (synthetic, salt-susceptible) and Zentos (elite cultivar, salt-tolerant) were integrated to gain a comprehensive understanding of the salt stress response. The MACE (Massive Analysis of cDNA 3'-Ends) based transcriptome analysis until 4 h after stress exposure revealed among the salt-responsive genes, the over-representation of genes coding for calcium-binding proteins. The functional and structural diversity within this category was studied and linked with the expression levels during the osmotic phase in the contrasting genotypes. The non-EF-hand category from calcium-binding proteins was found to be enriched for the susceptibility response. On the other side, the tolerant genotype was characterized by a faster and higher up-regulation of genes coding for proteins with EF-hand domain, such as RBOHD orthologs, and TF members. This study suggests that the interplay of calcium-binding proteins, WRKY, and AP2/ERF TF families in signaling pathways at the start of the osmotic phase can affect the expression of downstream genes. The identification of SNPs in promoter sequences and 3' -UTR regions provides insights into the molecular mechanisms controlling the differential expression of these genes through differential transcription factor binding affinity or altered mRNA stability.

A stepwise guide for pangenome development in crop plants: an alfalfa (Medicago sativa) case study

BACKGROUND

The concept of pangenomics and the importance of structural variants is gaining recognition within the plant genomics community. Due to advancements in sequencing and computational technology, it has become feasible to sequence the entire genome of numerous individuals of a single species at a reasonable cost. Pangenomes have been constructed for many major diploid crops, including rice, maize, soybean, sorghum, pearl millet, peas, sunflower, grapes, and mustards. However, pangenomes for polyploid species are relatively scarce and are available in only few crops including wheat, cotton, rapeseed, and potatoes.

MAIN BODY

In this review, we explore the various methods used in crop pangenome development, discussing the challenges and implications of these techniques based on insights from published pangenome studies. We offer a systematic guide and discuss the tools available for constructing a pangenome and conducting downstream analyses. Alfalfa, a highly heterozygous, cross pollinated and autotetraploid forage crop species, is used as an example to discuss the concerns and challenges offered by polyploid crop species. We conducted a comparative analysis using linear and graph-based methods by constructing an alfalfa graph pangenome using three publicly available genome assemblies. To illustrate the intricacies captured by pangenome graphs for a complex crop genome, we used five different gene sequences and aligned them against the three graph-based pangenomes. The comparison of the three graph pangenome methods reveals notable variations in the genomic variation captured by each pipeline.

CONCLUSION

Pangenome resources are proving invaluable by offering insights into core and dispensable genes, novel gene discovery, and genome-wide patterns of variation. Developing user-friendly online portals for linear pangenome visualization has made these resources accessible to the broader scientific and breeding community. However, challenges remain with graph-based pangenomes including compatibility with other tools, extraction of sequence for regions of interest, and visualization of genetic variation captured in pangenome graphs. These issues necessitate further refinement of tools and pipelines to effectively address the complexities of polyploid, highly heterozygous, and cross-pollinated species.

Carbohydrate-active enzymes involved in rice cell wall metabolism

Plant cell walls are complex, multifunctional structures, built up of polysaccharides and proteins. The configuration and abundance of cell wall constituents determine cellular elongation and plant growth. The emphasis of this review is on rice, a staple crop with economic importance, serving as model for grasses/cereals. Recent advancements have contributed to a better understanding of the grass/cereal cell wall. This review brings together current knowledge of the organization and metabolism of the rice cell wall, and addresses gaps in the information regarding the cell wall and enzymes involved. Several cell wall fractions, including cellulose, mixed-linkage glucans, and glucuronoarabinoxylans, are well understood in rice and other grasses/grains. Conversely, there are still open questions and missing links in relation to xyloglucans, glucomannans, pectin, lignin, and arabinogalactan proteins. There is still a large and untapped potential to identify carbohydrate-active enzymes (CAZymes), to characterize their activity, and to elucidate their involvement in the metabolism of the mentioned cell wall fractions. This review highlights the involvement of carbohydrate-active enzymes in rice cell wall metabolism, providing an update of current understanding with the aim of demarcating research areas with potential for further investigations.

Natural alleles of LEAFY and WAPO1 interact to regulate spikelet number per spike in wheat

KEY MESSAGE

Specific combinations of LFY and WAPO1 natural alleles maximize spikelet number per spike in wheat. Spikelet number per spike (SNS) is an important yield component in wheat that determines the maximum number of grains that can be formed in a wheat spike. In wheat, loss-of-function mutations in LEAFY (LFY) or its interacting protein WHEAT ORTHOLOG OF APO1 (WAPO1) significantly reduce SNS by reducing the rate of formation of spikelet meristems. In previous studies, we identified a natural amino acid change in WAPO1 (C47F) that significantly increases SNS in hexaploid wheat. In this study, we searched for natural variants in LFY that were associated with differences in SNS and detected significant effects in the LFY-B region in a nested association mapping population. We generated a large mapping population and confirmed that the LFY-B polymorphism R80S is linked with the differences in SNS, suggesting that LFY-B is the likely causal gene. A haplotype analysis revealed two amino acid changes P34L and R80S, which were both enriched during wheat domestication and breeding suggesting positive selection. We also explored the interactions between the LFY and WAPO1 natural variants for SNS using biparental populations and identified significant interaction, in which the positive effect of the 80S and 34L alleles from LFY-B was only detected in the WAPO-A1 47F background but not in the 47C background. Based on these results, we propose that the allele combination WAPO-A1-47F/LFY-B 34L 80S can be used in wheat breeding programs to maximize SNS and increase grain yield potential in wheat.

Diverse epistatic effects in barley-powdery mildew interactions localize to host chromosome hotspots

Barley Mildew locus a (Mla) encodes a multi-allelic series of nucleotide-binding leucine-rich repeat (NLR) receptors that specify recognition to diverse cereal diseases. We exploited time-course transcriptome dynamics of barley and derived immune mutants infected with the powdery mildew fungus, Blumeria hordei (Bh), to infer gene effects governed by Mla6 and two other loci significant to disease development, Blufensin1 (Bln1), and Required for Mla6 resistance3 (rar3 = Sgt1 ΔKL308-309 ). Interactions of Mla6 and Bln1 resulted in diverse epistatic effects on the Bh-induced barley transcriptome, differential immunity to Pseudomonas syringae expressing the effector protease AvrPphB, and reaction to Bh. From a total of 468 barley NLRs, 115 were grouped under different gene effect models; genes classified under these models localized to host chromosome hotspots. The corresponding Bh infection transcriptome was classified into nine co-expressed modules, linking differential expression with pathogen structures, signifying that disease is regulated by an inter-organismal network that diversifies the response.

k-mer genome-wide association study for anthracnose and BCMV resistance in a Phaseolus vulgaris Andean Diversity Panel

Access to broad genomic resources and closely linked marker-trait associations for common beans (Phaseolus vulgaris L.) can facilitate development of improved varieties with increased yield, improved market quality traits, and enhanced disease resistance. The emergence of virulent races of anthracnose (caused by Colletotrichum lindemuthianum) and bean common mosaic virus (BCMV) highlight the need for improved methods to identify and incorporate pan-genomic variation in breeding for disease resistance. We sequenced the P. vulgaris Andean Diversity Panel (ADP) and performed a genome-wide association study (GWAS) to identify associations for resistance to BCMV and eight races of anthracnose. Historical single nucleotide polymorphism (SNP)-chip and phenotypic data enabled a three-way comparison between SNP-chip, reference-based whole genome shotgun sequence (WGS)-SNP, and reference-free k-mer (short nucleotide subsequence) GWAS. Across all traits, there was excellent concordance between SNP-chip, WGS-SNP, and k-mer GWAS results-albeit at a much higher marker resolution for the WGS data sets. Significant k-mer haplotype variation revealed selection of the linked I-gene and Co-u traits in North American breeding lines and cultivars. Due to structural variation, only 9.1 to 47.3% of the significantly associated k-mers could be mapped to the reference genome. Thus, to determine the genetic context of cis-associated k-mers, we generated draft whole genome assemblies of four ADP accessions and identified an expanded local repertoire of disease resistance genes associated with resistance to anthracnose and BCMV. With access to variant data in the context of a pan-genome, high resolution mapping of agronomic traits for common bean is now feasible.

Trans-crop applications of atypical R genes for multipathogen resistance

Genetic resistance to plant diseases is essential for global food security. Significant progress has been achieved for plant disease-resistance (R) genes comprising nucleotide-binding domain, leucine-rich repeat-containing receptors (NLRs), and membrane-localized receptor-like kinases or proteins (RLKs/RLPs), which we refer to as typical R genes. However, there is a knowledge gap in how non-receptor-type or atypical R genes contribute to plant immunity. Here, we summarize resources and technologies facilitating the study of atypical R genes, examine diverse atypical R proteins for broad-spectrum resistance, and outline potential approaches for trans-crop applications of atypical R genes. Studies of atypical R genes are important for a holistic understanding of plant immunity and the development of novel strategies in disease control and crop improvement.

Analysis of a global wheat panel reveals a highly diverse introgression landscape and provides evidence for inter-homoeologue chromosomal recombination

KEY MESSAGE

This study highlights the agronomic potential of rare introgressions, as demonstrated by a major QTL for powdery mildew resistance on chromosome 7D. It further shows evidence for inter-homoeologue recombination in wheat. Agriculturally important genes are often introgressed into crops from closely related donor species or landraces. The gene pool of hexaploid bread wheat (Triticum aestivum) is known to contain numerous such "alien" introgressions. Recently established high-quality reference genome sequences allow prediction of the size, frequency and identity of introgressed chromosome regions. Here, we characterise chromosomal introgressions in bread wheat using exome capture data from the WHEALBI collection. We identified 24,981 putative introgression segments of at least 2 Mb across 434 wheat accessions. Detailed study of the most frequent introgressions identified T. timopheevii or its close relatives as a frequent donor species. Importantly, 118 introgressions of at least 10 Mb were exclusive to single wheat accessions, revealing that large populations need to be studied to assess the total diversity of the wheat pangenome. In one case, a 14 Mb introgression in chromosome 7D, exclusive to cultivar Pamukale, was shown by QTL mapping to harbour a recessive powdery mildew resistance gene. We identified multiple events where distal chromosomal segments of one subgenome were duplicated in the genome and replaced the homoeologous segment in another subgenome. We propose that these examples are the results of inter-homoeologue recombination. Our study produced an extensive catalogue of the wheat introgression landscape, providing a resource for wheat breeding. Of note, the finding that the wheat gene pool contains numerous rare, but potentially important introgressions and chromosomal rearrangements has implications for future breeding.

Tracing post-domestication historical events and screening pre-breeding germplasm from large gene pools in wheat in the absence of phenotype data

Wheat, particularly common wheat (Triticum aestivum L.), is a major crop accounting for 25% of the world cereal production and thriving in diverse ecogeographic regions. Its adaptation to diverse environments arises from its three distinct genomes adapted to different environments and post-domestication anthropogenic interventions. In search of key genomic regions revealing historic events and breeding significance to common wheat, we performed genome scan and genome-environment association (GEA) analyses using high-marker density genotype datasets. Whole-genome scans revealed highly differentiated regions on chromosomes 2A, 3B, and 4A. In-depth analyses corroborated our previous prediction of the 4A differentiated region signifying the separation between Spelt/Macha and other wheat types. Individual chromosome scans captured key introgressions, including one from T. timopheevii and one from Thinopyrum ponticum on 2B and 3D, respectively, as well as known genes such as Vrn-A1 on 5A. GEA highlighted loci linked to latitude-induced environmental variations, influencing traits such as photoperiodism and responses to abiotic stress. Variation at the Vrn-A1 locus on 5A assigned accessions to two haplotypes (6% and 94%). Further analysis on Vrn-A1 coding gene revealed four subgroups of the major haplotype, while the minor haplotype remained undifferentiated. Analyses at differentiated loci mostly dichotomized the population, illustrating the possibility of isolating pre-breeding materials with desirable traits from large gene pools in the absence of phenotype data. Given the current availability of broad genetic data, the genome-scan-GEA hybrid can be an efficient and cost-effective approach for pinpointing environmentally resilient pre-breeding germplasm from vast gene pools, including gene banks regardless of trait characterization.

On the evolution and genetic diversity of the bread wheat D genome

Bread wheat (Triticum aestivum) became a globally dominant crop after incorporating the D genome from the donor species Aegilops tauschii, but the evolutionary history that shaped the D genome during this process remains to be clarified. Here, we propose a renewed evolutionary model linking Ae. tauschii and the hexaploid wheat D genome by constructing an ancestral haplotype map covering 762 Ae. tauschii and hexaploid wheat accessions. We dissected the evolutionary trajectories of Ae. tauschii lineages and reported a few independent intermediate accessions, demonstrating that low-frequency inter-sublineage gene flow had enriched the diversity of Ae. tauschii. We discovered that the D genome of hexaploid wheat was inherited from a unified ancestral template, but with a mosaic composition that was highly mixed and derived mainly from three Ae. tauschii L2 sublineages located in the Caspian coastal region. This result suggests that early agricultural activities facilitated innovations in D-genome composition and finalized the success of hexaploidization. We found that the majority (51.4%) of genetic diversity was attributed to novel mutations absent in Ae. tauschii, and we identified large Ae. tauschii introgressions from various lineages, which expanded the diversity of the wheat D genome and introduced beneficial alleles. This work sheds light on the process of wheat hexaploidization and highlights the evolutionary significance of the multi-layered genetic diversity of the bread wheat D genome.

Wheat improvement through advances in single nucleotide polymorphism (SNP) detection and genotyping with a special emphasis on rust resistance

KEY MESSAGE

Single nucleotide polymorphism (SNP) markers in wheat and their prospects in breeding with special reference to rust resistance. Single nucleotide polymorphism (SNP)-based markers are increasingly gaining momentum for screening and utilizing vital agronomic traits in wheat. To date, more than 260 million SNPs have been detected in modern cultivars and landraces of wheat. This rapid SNP discovery was made possible through the release of near-complete reference and pan-genome assemblies of wheat and its wild relatives, coupled with whole genome sequencing (WGS) of thousands of wheat accessions. Further, genotyping customized SNP sites were facilitated by a series of arrays (9 to 820Ks), a cost effective substitute WGS. Lately, germplasm-specific SNP arrays have been introduced to characterize novel traits and detect closely linked SNPs for marker-assisted breeding. Subsequently, the kompetitive allele-specific PCR (KASP) assay was introduced for rapid and large-scale screening of specific SNP markers. Moreover, with the advances and reduction in sequencing costs, ample opportunities arise for generating SNPs artificially through mutations and in combination with next-generation sequencing and comparative genomic analyses. In this review, we provide historical developments and prospects of SNP markers in wheat breeding with special reference to rust resistance where over 50 genetic loci have been characterized through SNP markers. Rust resistance is one of the most essential traits for wheat breeding as new strains of the Puccinia fungus, responsible for rust diseases, evolve frequently and globally.

Hi-C techniques: from genome assemblies to transcription regulation

The invention of chromosome conformation capture (3C) techniques, in particular the key method Hi-C providing genome-wide information about chromatin contacts, revolutionized the way we study the three-dimensional organization of the nuclear genome and how it affects transcription, replication, and DNA repair. Because the frequency of chromatin contacts between pairs of genomic segments predictably relates to the distance in the linear genome, the information obtained by Hi-C has also proved useful for scaffolding genomic sequences. Here, we review recent improvements in experimental procedures of Hi-C and its various derivatives, such as Micro-C, HiChIP, and Capture Hi-C. We assess the advantages and limitations of the techniques, and present examples of their use in recent plant studies. We also report on progress in the development of computational tools used in assembling genome sequences.

Genome architecture of the allotetraploid wild grass Aegilops ventricosa reveals its evolutionary history and contributions to wheat improvement

The allotetraploid wild grass Aegilops ventricosa (2n = 4x = 28, genome DvDvNvNv) has been recognized as an important germplasm resource for wheat improvement owing to its ability to tolerate biotic stresses. In particular, the 2NvS segment from Ae. ventricosa, as a stable and effective resistance source, has contributed greatly to wheat improvement. The 2NvS/2AS translocation is a prevalent chromosomal translocation between common wheat and wild relatives, ranking just behind the 1B/1R translocation in importance for modern wheat breeding. Here, we assembled a high-quality chromosome-level reference genome of Ae. ventricosa RM271 with a total length of 8.67 Gb. Phylogenomic analyses revealed that the progenitor of the Dv subgenome of Ae. ventricosa is Ae. tauschii ssp. tauschii (genome DD); by contrast, the progenitor of the D subgenome of bread wheat (Triticum aestivum L.) is Ae. tauschii ssp. strangulata (genome DD). The oldest polyploidization time of Ae. ventricosa occurred ∼0.7 mya. The Dv subgenome of Ae. ventricosa is less conserved than the D subgenome of bread wheat. Construction of a graph-based pangenome of 2AS/6NvL (originally known as 2NvS) segments from Ae. ventricosa and other genomes in the Triticeae enabled us to identify candidate resistance genes sourced from Ae. ventricosa. We identified 12 nonredundant introgressed segments from the Dv and Nv subgenomes using a large winter wheat collection representing the full diversity of the European wheat genetic pool, and 29.40% of European wheat varieties inherit at least one of these segments. The high-quality RM271 reference genome will provide a basis for cloning key genes, including the Yr17-Lr37-Sr38-Cre5 resistance gene cluster in Ae. ventricosa, and facilitate the full use of elite wild genetic resources to accelerate wheat improvement.

Expression patterns of candidate genes for the Lr46/Yr29 "slow rust" locus in common wheat (Triticum aestivum L.) and associated miRNAs inform of the gene conferring the Puccinia triticina resistance trait

Leaf rust caused by Puccinia triticina (Pt) is one of the most impactful diseases causing substantial losses in common wheat (Triticum aestivum L.) crops. In adult plants resistant to Pt, a horizontal adult plant resistance (APR) is observed: APR protects the plant against multiple pathogen races and is distinguished by durable persistence under production conditions. The Lr46/Yr29 locus was mapped to chromosome 1B of common wheat genome, but the identity of the underlying gene has not been demonstrated although several candidate genes have been proposed. This study aimed to analyze the expression of nine candidate genes located at the Lr46/Yr29 locus and their four complementary miRNAs (tae-miR5384-3p, tae-miR9780, tae-miR9775, and tae-miR164), in response to Pt infection. The plant materials tested included five reference cultivars in which the molecular marker csLV46G22 associated with the Lr46/Yr29-based Pt resistance was identified, as well as one susceptible control cultivar. Biotic stress was induced in adult plants by inoculation with fungal spores under controlled conditions. Plant material was sampled before and at 6, 12, 24, 48 hours post inoculation (hpi). Differences in expression of candidate genes at the Lr46/Yr29 locus were analyzed by qRT-PCR and showed that the expression of the genes varied at the analyzed time points. The highest expression of Lr46/Yr29 candidate genes (Lr46-Glu1, Lr46-Glu2, Lr46-Glu3, Lr46-RLK1, Lr46-RLK2, Lr46-RLK3, Lr46-RLK4, Lr46-Snex, and Lr46-WRKY) occurred at 12 and 24 hpi and such expression profiles were obtained only for one candidate gene among the nine genes analyzed (Lr46-Glu2), indicating that it may be a contributing factor in the resistance response to Pt infection.

Facts, uncertainties, and opportunities in wheat molecular improvement

The year 2020 was a landmark year for wheat. The wheat HB4 event harboring a drought-resistant gene from sunflowers, received regulatory approval and was grown commercially in Argentina, with approval for food and feed in other countries. This, indeed, is many years after the adoption of genetic modifications in other crops. The lack of consumer acceptance and resulting trade barriers halted the commercialization of the earliest events and had a chilling effect on, especially, private Research & Development (R&D) investments. As regulations for modern breeding technologies such as genome-edited cultivars are being discussed and/or adopted across the globe, we would like to propose a framework to ensure that wheat is not left behind a second time as the potential benefits far outweigh the perceived risks. In this paper, after a review of the technical challenges wheat faces with the generation of trans- and cis-genic wheat varieties, we discuss some of the factors that could help demystify the risk/reward equation and thereby the consumer's reluctance or acceptance of these techniques for future wheat improvement. The advent of next-generation sequencing is shedding light on natural gene transfer between species and the number of perturbations other accepted techniques like mutagenesis create. The transition from classic breeding techniques and embracing transgenic, cisgenic, and genome editing approaches feels inevitable for wheat improvement if we are to develop climate-resilient wheat varieties to feed a growing world population.

Mapping rust resistance in European winter wheat: many QTLs for yellow rust resistance, but only a few well characterized genes for stem rust resistance

KEY MESSAGE

Stem rust resistance was mainly based on a few, already known resistance genes; for yellow rust resistance there was a combination of designated genes and minor QTLs. Yellow rust (YR) caused by Puccinia striiformis f. sp. tritici (Pst) and stem rust (SR) caused by Puccinia graminis f. sp. tritici (Pgt) are among the most damaging wheat diseases. Although, yellow rust has occurred regularly in Europe since the advent of the Warrior race in 2011, damaging stem rust epidemics are still unusual. We analyzed the resistance of seven segregating populations at the adult growth stage with the parents being selected for YR and SR resistances across three to six environments (location-year combinations) following inoculation with defined Pst and Pgt races. In total, 600 progenies were phenotyped and 563 were genotyped with a 25k SNP array. For SR resistance, three major resistance genes (Sr24, Sr31, Sr38/Yr17) were detected in different combinations. Additional QTLs provided much smaller effects except for a gene on chromosome 4B that explained much of the genetic variance. For YR resistance, ten loci with highly varying percentages of explained genetic variance (pG, 6-99%) were mapped. Our results imply that introgression of new SR resistances will be necessary for breeding future rust resistant cultivars, whereas YR resistance can be achieved by genomic selection of many of the detected QTLs.

The era of panomics-driven gene discovery in plants

Panomics is an approach to integrate multiple 'omics' datasets, generated using different individuals or natural variations. Considering their diverse phenotypic spectrum, the phenome is inherently associated with panomics-based science, which is further combined with genomics, transcriptomics, metabolomics, and other omics techniques, either independently or collectively. Panomics has been accelerated through recent technological advancements in the field of genomics that enable the detection of population-wide structural variations (SVs) and hence offer unprecedented insights into the genetic variations contributing to important agronomic traits. The present review provides the recent trends of panomics-driven gene discovery toward various traits related to plant development, stress tolerance, accumulation of specialized metabolites, and domestication/dedomestication. In addition, the success stories are highlighted in the broader context of enhancing crop productivity.

Engineering a One Health Super Wheat

Wheat is the predominant crop worldwide, contributing approximately 20% of protein and calories to the human diet. However, the yield potential of wheat faces limitations due to pests, diseases, and abiotic stresses. Although conventional breeding has improved desirable traits, the use of modern transgenesis technologies has been limited in wheat in comparison to other crops such as maize and soybean. Recent advances in wheat gene cloning and transformation technology now enable the development of a super wheat consistent with the One Health goals of sustainability, food security, and environmental stewardship. This variety combines traits to enhance pest and disease resistance, elevate grain nutritional value, and improve resilience to climate change. In this review, we explore ways to leverage current technologies to combine and transform useful traits into wheat. We also address the requirements of breeders and legal considerations such as patents and regulatory issues.

The prominent multiplication of Japanese soil-borne wheat mosaic virus co-infected with barley yellow mosaic virus in barley

Various members of the viral genera Furovirus and Bymovirus are damaging pathogens of a range of crop species. Infection of the soil-borne plasmodiophorid Polymyxa graminis transmits both Japanese soil-borne wheat mosaic virus (JSBWMV) and the barley yellow mosaic virus (BaYMV) to barley, but their interaction during an episode of their co-infection has not been characterized to date. Here, we present an analysis of the titer of JSBWMV and BaYMV in plants of winter barley growing over a five-month period from late fall until mid-spring. Although JSBWMV was detectable in the plants' roots four weeks earlier than BaYMV, the translocation of both viruses from the root to the leaves occurred nearly simultaneously. Both viruses were co-localized in the roots, leaf sheathes, and leaf blades; however, in some stripes of leaf veins where infection by JSBWMV was prominent, BaYMV was not detectable. A substantial titer of both viruses persisted until early spring, after which JSBWMV became more prominent, being in a range of 10 to 100 times abundant of BaYMV. However, JSBWMV was only able to infect a single wheat accession (cv. Norin 61), whereas all of the wheat entries assayed appeared to be immune to BaYMV infection. Overall, our findings highlight the importance of resistance mechanisms against soil-borne viruses in cereal crops, expanding our understanding of plant-virus interactions and potentially informing strategies for crop protection against viral pathogens.

Origin and evolution of the bread wheat D genome

Bread wheat (Triticum aestivum) is a globally dominant crop and major source of calories and proteins for the human diet. Compared with its wild ancestors, modern bread wheat shows lower genetic diversity, caused by polyploidisation, domestication and breeding bottlenecks1,2. Wild wheat relatives represent genetic reservoirs, and harbour diversity and beneficial alleles that have not been incorporated into bread wheat. Here we establish and analyse extensive genome resources for Tausch's goatgrass (Aegilops tauschii), the donor of the bread wheat D genome. Our analysis of 46 Ae. tauschii genomes enabled us to clone a disease resistance gene and perform haplotype analysis across a complex disease resistance locus, allowing us to discern alleles from paralogous gene copies. We also reveal the complex genetic composition and history of the bread wheat D genome, which involves contributions from genetically and geographically discrete Ae. tauschii subpopulations. Together, our results reveal the complex history of the bread wheat D genome and demonstrate the potential of wild relatives in crop improvement.

Genomic insights into the modifications of spike morphology traits during wheat breeding

Over the past century, environmental changes have significantly impacted wheat spike morphology, crucial for adaptation and grain yield. However, the changes in wheat spike modifications during this period remain largely unknown. This study examines 16 spike morphology traits in 830 accessions released from 1900 to 2020. It finds that spike weight, grain number per spike (GN), and thousand kernel weight have significantly increased, while spike length has no significant change. The increase in fertile spikelets is due to fewer degenerated spikelets, resulting in a higher GN. Genome-wide association studies identified 49,994 significant SNPs, grouped into 293 genomic regions. The accumulation of favorable alleles in these genomic regions indicates the genetic basis for modification in spike morphology traits. Genetic network analysis of these genomic regions reveals the genetic basis for phenotypic correlations among spike morphology traits. The haplotypes of the identified genomic regions display obvious geographical differentiation in global accessions and environmental adaptation over the past 120 years. In summary, we reveal the genetic basis of adaptive evolution and the interactions of spike morphology, offering valuable resources for the genetic improvement of spike morphology to enhance environmental adaptation.

Genome-wide characterization and expression analysis of the CINNAMYL ALCOHOL DEHYDROGENASE gene family in Triticum aestivum

BACKGROUND

CINNAMYL ALCOHOL DEHYDROGENASE (CAD) catalyzes the NADPH-dependent reduction of cinnamaldehydes into cinnamyl alcohols and is a key enzyme found at the final step of the monolignol pathway. Cinnamyl alcohols and their conjugates are subsequently polymerized in the secondary cell wall to form lignin. CAD genes are typically encoded by multi-gene families and thus traditionally organized into general classifications of functional relevance.

RESULTS

In silico analysis of the hexaploid Triticum aestivum genome revealed 47 high confidence TaCAD copies, of which three were determined to be the most significant isoforms (class I) considered bone fide CADs. Class I CADs were expressed throughout development both in RNAseq data sets as well as via qRT-PCR analysis. Of the 37 class II TaCADs identified, two groups were observed to be significantly co-expressed with class I TaCADs in developing tissue and under chitin elicitation in RNAseq data sets. These co-expressed class II TaCADs were also found to be phylogenetically unrelated to a separate clade of class II TaCADs previously reported to be an influential resistance factor to pathogenic fungal infection. Lastly, two groups were phylogenetically identified as class III TaCADs, which possess distinct conserved gene structures. However, the lack of data supporting their catalytic activity for cinnamaldehydes and their bereft transcriptional presence in lignifying tissues challenges their designation and function as CADs.

CONCLUSIONS

Taken together, our comprehensive transcriptomic analyses suggest that TaCAD genes contribute to overlapping but nonredundant functions during T. aestivum growth and development across a wide variety of agroecosystems and provide tolerance to various stressors.

Wheat genomics: genomes, pangenomes, and beyond

There is an urgent need to improve wheat for upcoming challenges, including biotic and abiotic stresses. Sustainable wheat improvement requires the introduction of new genes and alleles in high-yielding wheat cultivars. Using new approaches, tools, and technologies to identify and introduce new genes in wheat cultivars is critical. High-quality genomes, transcriptomes, and pangenomes provide essential resources and tools to examine wheat closely to identify and manipulate new and targeted genes and alleles. Wheat genomics has improved excellently in the past 5 years, generating multiple genomes, pangenomes, and transcriptomes. Leveraging these resources allows us to accelerate our crop improvement pipelines. This review summarizes the progress made in wheat genomics and trait discovery in the past 5 years.

Fine mapping of stem rust resistance derived from soft red winter wheat cultivar AGS2000 to an NLR gene cluster on chromosome 6D

The Puccinia graminis f. sp. tritici (Pgt) Ug99-emerging virulent races present a major challenge to global wheat production. To meet present and future needs, new sources of resistance must be found. Identification of markers that allow tracking of resistance genes is needed for deployment strategies to combat highly virulent pathogen races. Field evaluation of a DH population located a QTL for stem rust (Sr) resistance, QSr.nc-6D from the breeding line MD01W28-08-11 to the distal region of chromosome arm 6DS where Sr resistance genes Sr42, SrCad, and SrTmp have been identified. A locus for seedling resistance to Pgt race TTKSK was identified in a DH population and an RIL population derived from the cross AGS2000 × LA95135. The resistant cultivar AGS2000 is in the pedigree of MD01W28-08-11 and our results suggest that it is the source of Sr resistance in this breeding line. We exploited published markers and exome capture data to enrich marker density in a 10 Mb region flanking QSr.nc-6D. Our fine mapping in heterozygous inbred families identified three markers co-segregating with resistance and delimited QSr.nc-6D to a 1.3 Mb region. We further exploited information from other genome assemblies and identified collinear regions of 6DS harboring clusters of NLR genes. Evaluation of KASP assays corresponding to our co-segregating SNP suggests that they can be used to track this Sr resistance in breeding programs. However, our results also underscore the challenges posed in identifying genes underlying resistance in such complex regions in the absence of genome sequence from the resistant genotypes.

Identification of stripe rust resistance gene YrBDT in Chinese landrace wheat Baidatou using BSE-seq and BSR-seq

KEY MESSAGE

A new stripe rust resistance gene YrBDT in Chinese landrace wheat Baidatou was mapped to a 943.6-kb interval on chromosome arm 6DS and co-segregated with a marker CAPS3 developed from candidate gene TraesCS6D03G0027300. Stripe rust caused by Puccinia striiformis f. sp. tritici (Pst) is a devastating foliar disease of wheat. Chinese landrace wheat Baidatou has shown high resistance to a broad spectrum of Pst races at both the seedling and adult-plant stages for decades in the Longnan region of Gansu province, a hot spot for stripe rust epidemics. Here, we report fine mapping and candidate gene analysis of stripe rust resistance gene YrBDT in Baidatou. Analysis of F1, F2 plants and F2:3 lines indicated that resistance in Baidatou to Pst race CYR31 was conferred by a single dominant gene, temporarily designated YrBDT. Bulked segregant exome capture sequencing (BSE-seq) analysis revealed 61 high-confidence polymorphic SNPs concentrated in a 5.4-Mb interval at the distal of chromosome arm 6DS. Several SNPs and InDels were also identified by genome mining of DNA sampled from the parents and contrasting bulks. The YrBDT locus was mapped to a 943.6-kb (4,658,322-5,601,880 bp) genomic region spanned by markers STS2 and STS3 based on IWGSC RefSeq v2.1, including five putative disease resistance genes. There was high collinearity of the target interval among Chinese Spring RefSeq v2.1, Ae. tauschii AL8/78 and Fielder genomes. The expression level of TraesCS6D03G0027300 showed significant association with Pst infection, and a gene-specific marker CAPS3 developed from TraesCS6D03G0027300 co-segregated with YrBDT suggesting this gene as a candidate of YrBDT. The resistance gene and flanking markers can be used in marker-assisted selection for improvement of stripe rust resistance.

Evolution of Einkorn wheat centromeres is driven by the mutualistic interplay of two LTR retrotransposons

BACKGROUND

Centromere function is highly conserved across eukaryotes, but the underlying centromeric DNA sequences vary dramatically between species. Centromeres often contain a high proportion of repetitive DNA, such as tandem repeats and/or transposable elements (TEs). Einkorn wheat centromeres lack tandem repeat arrays and are instead composed mostly of the two long terminal repeat (LTR) retrotransposon families RLG_Cereba and RLG_Quinta which specifically insert in centromeres. However, it is poorly understood how these two TE families relate to each other and if and how they contribute to centromere function and evolution.

RESULTS

Based on conservation of diagnostic motifs (LTRs, integrase and primer binding site and polypurine-tract), we propose that RLG_Cereba and RLG_Quinta are a pair of autonomous and non-autonomous partners, in which the autonomous RLG_Cereba contributes all the proteins required for transposition, while the non-autonomous RLG_Quinta contributes GAG protein. Phylogenetic analysis of predicted GAG proteins showed that the RLG_Cereba lineage was present for at least 100 million years in monocotyledon plants. In contrast, RLG_Quinta evolved from RLG_Cereba between 28 and 35 million years ago in the common ancestor of oat and wheat. Interestingly, the integrase of RLG_Cereba is fused to a so-called CR-domain, which is hypothesized to guide the integrase to the functional centromere. Indeed, ChIP-seq data and TE population analysis show only the youngest subfamilies of RLG_Cereba and RLG_Quinta are found in the active centromeres. Importantly, the LTRs of RLG_Quinta and RLG_Cereba are strongly associated with the presence of the centromere-specific CENH3 histone variant. We hypothesize that the LTRs of RLG_Cereba and RLG_Quinta contribute to wheat centromere integrity by phasing and/or placing CENH3 nucleosomes, thus favoring their persistence in the competitive centromere-niche.

CONCLUSION

Our data show that RLG_Cereba cross-mobilizes the non-autonomous RLG_Quinta retrotransposons. New copies of both families are specifically integrated into functional centromeres presumably through direct binding of the integrase CR domain to CENH3 histone variants. The LTRs of newly inserted RLG_Cereba and RLG_Quinta elements, in turn, recruit and/or phase new CENH3 deposition. This mutualistic interplay between the two TE families and the plant host dynamically maintains wheat centromeres.

PanVA: Pangenomic Variant Analysis

Genomics researchers increasingly use multiple reference genomes to comprehensively explore genetic variants underlying differences in detectable characteristics between organisms. Pangenomes allow for an efficient data representation of multiple related genomes and their associated metadata. However, current visual analysis approaches for exploring these complex genotype-phenotype relationships are often based on single reference approaches or lack adequate support for interpreting the variants in the genomic context with heterogeneous (meta)data. This design study introduces PanVA, a visual analytics design for pangenomic variant analysis developed with the active participation of genomics researchers. The design uniquely combines tailored visual representations with interactions such as sorting, grouping, and aggregation, allowing users to navigate and explore different perspectives on complex genotype-phenotype relations. Through evaluation in the context of plants and pathogen research, we show that PanVA helps researchers explore variants in genes and generate hypotheses about their role in phenotypic variation.

Plant pangenomes for crop improvement, biodiversity and evolution

Plant genome sequences catalogue genes and the genetic elements that regulate their expression. Such inventories further research aims as diverse as mapping the molecular basis of trait diversity in domesticated plants or inquiries into the origin of evolutionary innovations in flowering plants millions of years ago. The transformative technological progress of DNA sequencing in the past two decades has enabled researchers to sequence ever more genomes with greater ease. Pangenomes - complete sequences of multiple individuals of a species or higher taxonomic unit - have now entered the geneticists' toolkit. The genomes of crop plants and their wild relatives are being studied with translational applications in breeding in mind. But pangenomes are applicable also in ecological and evolutionary studies, as they help classify and monitor biodiversity across the tree of life, deepen our understanding of how plant species diverged and show how plants adapt to changing environments or new selection pressures exerted by human beings.

Harnessing landrace diversity empowers wheat breeding

Harnessing genetic diversity in major staple crops through the development of new breeding capabilities is essential to ensure food security1. Here we examined the genetic and phenotypic diversity of the A. E. Watkins landrace collection2 of bread wheat (Triticum aestivum), a major global cereal, by whole-genome re-sequencing of 827 Watkins landraces and 208 modern cultivars and in-depth field evaluation spanning a decade. We found that modern cultivars are derived from two of the seven ancestral groups of wheat and maintain very long-range haplotype integrity. The remaining five groups represent untapped genetic sources, providing access to landrace-specific alleles and haplotypes for breeding. Linkage disequilibrium-based haplotypes and association genetics analyses link Watkins genomes to the thousands of identified high-resolution quantitative trait loci and significant marker-trait associations. Using these structured germplasm, genotyping and informatics resources, we revealed many Watkins-unique beneficial haplotypes that can confer superior traits in modern wheat. Furthermore, we assessed the phenotypic effects of 44,338 Watkins-unique haplotypes, introgressed from 143 prioritized quantitative trait loci in the context of modern cultivars, bridging the gap between landrace diversity and current breeding. This study establishes a framework for systematically utilizing genetic diversity in crop improvement to achieve sustainable food security.

How Important Are Genetic Diversity and Cultivar Uniformity in Wheat? The Case of Gliadins

Improvements in self-pollinated crops rely on crosses between different genotypes. It has been suggested that the repeated use of "the best" genotypes may lead to the restriction of the genetic diversity of the crop. In wheat, the analysis of gliadin (storage protein) polymorphism has provided evidence that genetic diversity was high and stable throughout the 20th century. Moreover, a worldwide analysis of gliadin polymorphism shows that genetic diversity is structured spatially across countries and their regions. Therefore, the analysis of gliadin genotypes in a given grain sample can provide reliable information about the origin of grains in this sample. An unexpected finding is that many registered common wheat cultivars are genetically non-uniform and composed of authentic biotypes (genotypically related lines originated from the initial cross) in spite of current crop-registration rules that include a strict demand for each new cultivar to be genetically uniform (DUS rules). In summary, the results suggest that each cultivar is the fruit of joint effects of a breeder and of a region's environmental factors. We believe this finding will not be restricted to wheat and suggest there may be a need to re-evaluate relevant rules of cultivar registration for crop species in general.

Genome-wide atlas of rust resistance loci in wheat

Rust diseases, including leaf rust, stripe/yellow rust, and stem rust, significantly impact wheat (Triticum aestivum L.) yields, causing substantial economic losses every year. Breeding and deployment of cultivars with genetic resistance is the most effective and sustainable approach to control these diseases. The genetic toolkit for wheat breeders to select for rust resistance has rapidly expanded with a multitude of genetic loci identified using the latest advances in genomics, mapping and cloning strategies. The goal of this review was to establish a wheat genome atlas that provides a comprehensive summary of reported loci associated with rust resistance. Our atlas provides a summary of mapped quantitative trait loci (QTL) and characterised genes for the three rusts from 170 publications over the past two decades. A total of 920 QTL or resistance genes were positioned across the 21 chromosomes of wheat based on the latest wheat reference genome (IWGSC RefSeq v2.1). Interestingly, 26 genomic regions contained multiple rust loci suggesting they could have pleiotropic effects on two or more rust diseases. We discuss a range of strategies to exploit this wealth of genetic information to efficiently utilise sources of resistance, including genomic information to stack desirable and multiple QTL to develop wheat cultivars with enhanced resistance to rust disease.

Methyltransferase TaSAMT1 mediates wheat freezing tolerance by integrating brassinosteroid and salicylic acid signaling

Cold injury is a major environmental stress affecting the growth and yield of crops. Brassinosteroids (BRs) and salicylic acid (SA) play important roles in plant cold tolerance. However, whether or how BR signaling interacts with the SA signaling pathway in response to cold stress is still unknown. Here, we identified an SA methyltransferase, TaSAMT1 that converts SA to methyl SA (MeSA) and confers freezing tolerance in wheat (Triticum aestivum). TaSAMT1 overexpression greatly enhanced wheat freezing tolerance, with plants accumulating more MeSA and less SA, whereas Tasamt1 knockout lines were sensitive to freezing stress and accumulated less MeSA and more SA. Spraying plants with MeSA conferred freezing tolerance to Tasamt1 mutants, but SA did not. We revealed that BRASSINAZOLE-RESISTANT 1 (TaBZR1) directly binds to the TaSAMT1 promoter and induces its transcription. Moreover, TaBZR1 interacts with the histone acetyltransferase TaHAG1, which potentiates TaSAMT1 expression via increased histone acetylation and modulates the SA pathway during freezing stress. Additionally, overexpression of TaBZR1 or TaHAG1 altered TaSAMT1 expression and improved freezing tolerance. Our results demonstrate a key regulatory node that connects the BR and SA pathways in the plant cold stress response. The regulatory factors or genes identified could be effective targets for the genetic improvement of freezing tolerance in crops.

High-Copy Transposons from a Pathogen Give Rise to a Conserved sRNA Family with a Novel Host Immunity Target

Small RNAs (sRNAs) are involved in gene silencing in multiple ways, including through cross-kingdom transfers from parasites to their hosts. Little is known about the evolutionary mechanisms enabling eukaryotic microbes to evolve functional mimics of host small regulatory RNAs. Here, we describe the identification and functional characterization of SINE_sRNA1, an sRNA family derived from highly abundant short interspersed nuclear element (SINE) retrotransposons in the genome of the wheat powdery mildew pathogen. SINE_sRNA1 is encoded by a sequence motif that is conserved in multiple SINE families and corresponds to a functional plant microRNA (miRNA) mimic targeting Tae_AP1, a wheat gene encoding an aspartic protease only found in monocots. Tae_AP1 has a novel function enhancing both pattern-triggered immunity (PTI) and effector-triggered immunity (ETI), thereby contributing to the cross activation of plant defenses. We conclude that SINE_sRNA1 and Tae_AP1 are functional innovations, suggesting the contribution of transposons to the evolutionary arms race between a parasite and its host. [Formula: see text] Copyright © 2024 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

A telomere-to-telomere gap-free reference genome assembly of avocado provides useful resources for identifying genes related to fatty acid biosynthesis and disease resistance

Avocado (Persea americana Mill.) is an economically valuable plant because of the high fatty acid content and unique flavor of its fruits. Its fatty acid content, especially the relatively high unsaturated fatty acid content, provides significant health benefits. We herein present a telomere-to-telomere gapless genome assembly (841.6 Mb) of West Indian avocado. The genome contains 40 629 predicted protein-coding genes. Repeat sequences account for 57.9% of the genome. Notably, all telomeres, centromeres, and a nucleolar organizing region are included in this genome. Fragments from these three regions were observed via fluorescence in situ hybridization. We identified 376 potential disease resistance-related nucleotide-binding leucine-rich repeat genes. These genes, which are typically clustered on chromosomes, may be derived from gene duplication events. Five NLR genes (Pa11g0262, Pa02g4855, Pa07g3139, Pa07g0383, and Pa02g3196) were highly expressed in leaves, stems, and fruits, indicating they may be involved in avocado disease responses in multiple tissues. We also identified 128 genes associated with fatty acid biosynthesis and analyzed their expression patterns in leaves, stems, and fruits. Pa02g0113, which encodes one of 11 stearoyl-acyl carrier protein desaturases mediating C18 unsaturated fatty acid synthesis, was more highly expressed in the leaves than in the stems and fruits. These findings provide valuable insights that enhance our understanding of fatty acid biosynthesis in avocado.

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.

Tagging large CNV blocks in wheat boosts digitalization of germplasm resources by ultra-low-coverage sequencing

BACKGROUND

The massive structural variations and frequent introgression highly contribute to the genetic diversity of wheat, while the huge and complex genome of polyploid wheat hinders efficient genotyping of abundant varieties towards accurate identification, management, and exploitation of germplasm resources.

RESULTS

We develop a novel workflow that identifies 1240 high-quality large copy number variation blocks (CNVb) in wheat at the pan-genome level, demonstrating that CNVb can serve as an ideal DNA fingerprinting marker for discriminating massive varieties, with the accuracy validated by PCR assay. We then construct a digitalized genotyping CNVb map across 1599 global wheat accessions. Key CNVb markers are linked with trait-associated introgressions, such as the 1RS·1BL translocation and 2NvS translocation, and the beneficial alleles, such as the end-use quality allele Glu-D1d (Dx5 + Dy10) and the semi-dwarf r-e-z allele. Furthermore, we demonstrate that these tagged CNVb markers promote a stable and cost-effective strategy for evaluating wheat germplasm resources with ultra-low-coverage sequencing data, competing with SNP array for applications such as evaluating new varieties, efficient management of collections in gene banks, and describing wheat germplasm resources in a digitalized manner. We also develop a user-friendly interactive platform, WheatCNVb ( http://wheat.cau.edu.cn/WheatCNVb/ ), for exploring the CNVb profiles over ever-increasing wheat accessions, and also propose a QR-code-like representation of individual digital CNVb fingerprint. This platform also allows uploading new CNVb profiles for comparison with stored varieties.

CONCLUSIONS

The CNVb-based approach provides a low-cost and high-throughput genotyping strategy for enabling digitalized wheat germplasm management and modern breeding with precise and practical decision-making.

Haplotype-resolved genomes of octoploid species in Phyllanthaceae family reveal a critical role for polyploidization and hybridization in speciation

The Phyllanthaceae family comprises a diverse range of plants with medicinal, edible, and ornamental value, extensively cultivated worldwide. Polyploid species commonly occur in Phyllanthaceae. Due to the rather complex genomes and evolutionary histories, their speciation process has been still lacking in research. In this study, we generated chromosome-scale haplotype-resolved genomes of two octoploid species (Phyllanthus emblica and Sauropus spatulifolius) in Phyllanthaceae family. Combined with our previously reported one tetraploid (Sauropus androgynus) and one diploid species (Phyllanthus cochinchinensis) from the same family, we explored their speciation history. The three polyploid species were all identified as allopolyploids with subgenome A/B. Each of their two distinct subgenome groups from various species was uncovered to independently share a common diploid ancestor (Ancestor-AA and Ancestor-BB). Via different evolutionary routes, comprising various scenarios of bifurcating divergence, allopolyploidization (hybrid polyploidization), and autopolyploidization, they finally evolved to the current tetraploid S. androgynus, and octoploid S. spatulifolius and P. emblica, respectively. We further discuss the variations in copy number of alleles and the potential impacts within the two octoploids. In addition, we also investigated the fluctuation of metabolites with medical values and identified the key factor in its biosynthesis process in octoploids species. Our study reconstructed the evolutionary history of these Phyllanthaceae species, highlighting the critical roles of polyploidization and hybridization in their speciation processes. The high-quality genomes of the two octoploid species provide valuable genomic resources for further research of evolution and functional genomics.

MicroRNAs potentially targeting DDR-related genes are differentially expressed upon exposure to γ-rays during seed germination in wheat

DNA damage response (DDR), a complex network of cellular pathways that cooperate to sense and repair DNA lesions, is regulated by several mechanisms, including microRNAs. As small, single-stranded RNA molecules, miRNAs post-transcriptionally regulate their target genes by mRNA cleavage or translation inhibition. Knowledge regarding miRNAs influence on DDR-associated genes is still scanty in plants. In this work, an in silico analysis was performed to identify putative miRNAs that could target DDR sensors, signal transducers and effector genes in wheat. Selected putative miRNA-gene pairs were tested in an experimental system where seeds from two wheat mutant lines were irradiated with 50 Gy and 300 Gy gamma(γ)-rays. To evaluate the effect of the treatments on wheat germination, phenotypic and molecular (DNA damage, ROS accumulation, gene/miRNA expression profile) analyses have been carried out. The results showed that in dry seeds ROS accumulated immediately after irradiation and decayed soon after while the negative impact on seedling growth was supported by enhanced accumulation of DNA damage. When a qRT-PCR analysis was performed, the selected miRNAs and DDR-related genes were differentially modulated by the γ-rays treatments in a dose-, time- and genotype-dependent manner. A significant negative correlation was observed between the expression of tae-miR5086 and the RAD50 gene, involved in double-strand break sensing and homologous recombination repair, one of the main processes that repairs DNA breaks induced by γ-rays. The results hereby reported can be relevant for wheat breeding programs and screening of the radiation response and tolerance of novel wheat varieties.

A k-mer-based pangenome approach for cataloging seed-storage-protein genes in wheat to facilitate genotype-to-phenotype prediction and improvement of end-use quality

Wheat is a staple food for more than 35% of the world's population, with wheat flour used to make hundreds of baked goods. Superior end-use quality is a major breeding target; however, improving it is especially time-consuming and expensive. Furthermore, genes encoding seed-storage proteins (SSPs) form multi-gene families and are repetitive, with gaps commonplace in several genome assemblies. To overcome these barriers and efficiently identify superior wheat SSP alleles, we developed "PanSK" (Pan-SSP k-mer) for genotype-to-phenotype prediction based on an SSP-based pangenome resource. PanSK uses 29-mer sequences that represent each SSP gene at the pangenomic level to reveal untapped diversity across landraces and modern cultivars. Genome-wide association studies with k-mers identified 23 SSP genes associated with end-use quality that represent novel targets for improvement. We evaluated the effect of rye secalin genes on end-use quality and found that removal of ω-secalins from 1BL/1RS wheat translocation lines is associated with enhanced end-use quality. Finally, using machine-learning-based prediction inspired by PanSK, we predicted the quality phenotypes with high accuracy from genotypes alone. This study provides an effective approach for genome design based on SSP genes, enabling the breeding of wheat varieties with superior processing capabilities and improved end-use quality.

New Observations of the Effects of the Cytoplasm of Aegilops kotschyi Boiss. in Bread Wheat Triticum aestivum L

The cytoplasm of Aegilops kotschyi is known for the induction of male sterility and haploidy in wheat. Both systems originally appeared rather simple, but manipulation of the standard chromosome constitution of the nuclear genome revealed additional interactions. This study shows that while there is little or no allelic variation at the main fertility restorer locus Rfmulti on chromosome arm 1BS, additional genes may also be involved in the nuclear-mitochondrial genome interactions, affecting not only male fertility but also the growth rate, from pollen competition for fertilization and early endosperm divisions all the way to seed size and plant maturity. Some of these effects appear to be of a sporophytic nature; others are gametophytic. Induction of parthenogenesis by a rye inducer in conjunction with the Ae. kotschyi cytoplasm is well known. However, here we show that the cytoplasmic-nuclear interactions affect all aspects of double fertilization: producing maternal haploids from unfertilized eggs, diploids from fertilized eggs or synergids, embryo-less kernels, and fertilized eggs without fertilization of the double nucleus in the embryo sack. It is unclear how frequent the inducers of parthenogenesis are, as variation, if any, is obscured by suppressors present in the wheat genome. Genetic dissection of a single wheat accession revealed five distinct loci affecting the rate of maternal haploid production: four acting as suppressors and one as an enhancer. Only when the suppressing haplotypes are confirmed may it be possible to the identify genetic variation of haploidy inducers, map their position(s), and determine their nature and the mode of action.

Impact of whole-genome duplications on structural variant evolution in Cochlearia

Polyploidy, the result of whole-genome duplication (WGD), is a major driver of eukaryote evolution. Yet WGDs are hugely disruptive mutations, and we still lack a clear understanding of their fitness consequences. Here, we study whether WGDs result in greater diversity of genomic structural variants (SVs) and how they influence evolutionary dynamics in a plant genus, Cochlearia (Brassicaceae). By using long-read sequencing and a graph-based pangenome, we find both negative and positive interactions between WGDs and SVs. Masking of recessive mutations due to WGDs leads to a progressive accumulation of deleterious SVs across four ploidal levels (from diploids to octoploids), likely reducing the adaptive potential of polyploid populations. However, we also discover putative benefits arising from SV accumulation, as more ploidy-specific SVs harbor signals of local adaptation in polyploids than in diploids. Together, our results suggest that SVs play diverse and contrasting roles in the evolutionary trajectories of young polyploids.

Chromosome-level assembly of the synthetic hexaploid wheat-derived cultivar Chuanmai 104

Synthetic hexaploid wheats (SHWs) are effective genetic resources for transferring agronomically important genes from wild relatives to common wheat (Triticum aestivum L.). Dozens of reference-quality pseudomolecule assemblies of hexaploid wheat have been generated, but none is reported for SHW-derived cultivars. Here, we generated a chromosome-scale assembly for the SHW-derived cultivar 'Chuanmai 104' based on PacBio HiFi reads and chromosome conformation capture sequencing. The total assembly size was 14.81 Gb with a contig N50 length of 58.25 Mb. A BUSCO analysis yielded a completeness score of 99.30%. In total, repetitive elements comprised 81.36% of the genome and 122,554 high-confidence protein-coding gene models were predicted. In summary, the first chromosome-level assembly for a SHW-derived cultivar presents a promising outlook for the study and utilization of SHWs in wheat improvement, which is essential to meet the global food demand.

Cicer super-pangenome provides insights into species evolution and agronomic trait loci for crop improvement in chickpea

Chickpea (Cicer arietinum L.)-an important legume crop cultivated in arid and semiarid regions-has limited genetic diversity. Efforts are being undertaken to broaden its diversity by utilizing its wild relatives, which remain largely unexplored. Here, we present the Cicer super-pangenome based on the de novo genome assemblies of eight annual Cicer wild species. We identified 24,827 gene families, including 14,748 core, 2,958 softcore, 6,212 dispensable and 909 species-specific gene families. The dispensable genome was enriched for genes related to key agronomic traits. Structural variations between cultivated and wild genomes were used to construct a graph-based genome, revealing variations in genes affecting traits such as flowering time, vernalization and disease resistance. These variations will facilitate the transfer of valuable traits from wild Cicer species into elite chickpea varieties through marker-assisted selection or gene-editing. This study offers valuable insights into the genetic diversity and potential avenues for crop improvement in chickpea.

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.

The wheat powdery mildew resistance gene Pm4 also confers resistance to wheat blast

Wheat blast, caused by the fungus Magnaporthe oryzae, threatens global cereal production since its emergence in Brazil in 1985 and recently spread to Bangladesh and Zambia. Here we demonstrate that the AVR-Rmg8 effector, common in wheat-infecting isolates, is recognized by the gene Pm4, previously shown to confer resistance to specific races of Blumeria graminis f. sp. tritici, the cause of powdery mildew of wheat. We show that Pm4 alleles differ in their recognition of different AVR-Rmg8 alleles, and some confer resistance only in seedling leaves but not spikes, making it important to select for those alleles that function in both tissues. This study has identified a gene recognizing an important virulence factor present in wheat blast isolates in Bangladesh and Zambia and represents an important first step towards developing durably resistant wheat cultivars for these regions.

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.

Genomic insights into glume pubescence in durum wheat: GWAS and haplotype analysis implicates TdELD1-1A as a candidate gene

Glume pubescence is an important morphological trait for the characterization of wheat cultivars. It shows tolerance to biotic and abiotic stresses to some extent. Hg1 (formerly named Hg) locus on chromosome 1AS controls glume pubescence in wheat. Its genetic analysis, fine-mapping and candidate gene analysis have been widely studied recently, however, the cloning of Hg1 has not yet been reported. Here, we conducted a GWAS between a dense panel of 171,103 SNPs and glume pubescence (Gp) in a durum wheat population of 145 lines, and further analyzed the candidate genes of Hg1 combined with the gene expression, functional annotation, and haplotype analysis. As a results, TRITD0Uv1G104670 (TdELD1-1A), encoding glycosyltransferase-like ELD1/KOBITO 1, was detected as the most promising candidate gene of Hg1 for glume pubescence in durum wheat. Our findings not only contribute to a deeper understanding of its cloning and functional validation but also underscore the significance of accurate genome sequences and annotations. Additionally, our study highlights the relevance of unanchored sequences in chrUn and the application of bioinformatics analysis for gene discovery in durum wheat.