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.

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.

Genome of the early spider-orchid Ophrys sphegodes provides insights into sexual deception and pollinator adaptation.

Pollinator-driven evolution of floral traits is thought to be a major driver of angiosperm speciation and diversification. Ophrys orchids mimic female insects to lure male pollinators into pseudocopulation. This strategy, called sexual deception, is species-specific, thereby providing strong premating reproductive isolation. Identifying the genomic architecture underlying pollinator adaptation and speciation may shed light on the mechanisms of angiosperm diversification. Here, we report the 5.2 Gb chromosome-scale genome sequence of Ophrys sphegodes. We find evidence for transposable element expansion that preceded the radiation of the O. sphegodes group, and for gene duplication having contributed to the evolution of chemical mimicry. We report a highly differentiated genomic candidate region for pollinator-mediated evolution on chromosome 2. The Ophrys genome will prove useful for investigations into the repeated evolution of sexual deception, pollinator adaptation and the genomic architectures that facilitate evolutionary radiations.

A fatty acid anabolic pathway in specialized-cells sustains a remote signal that controls egg activation in Drosophila.

Egg activation, representing the critical oocyte-to-embryo transition, provokes meiosis completion, modification of the vitelline membrane to prevent polyspermy, and translation of maternally provided mRNAs. This transition is triggered by a calcium signal induced by spermatozoon fertilization in most animal species, but not in insects. In Drosophila melanogaster, mature oocytes remain arrested at metaphase-I of meiosis and the calcium-dependent activation occurs while the oocyte moves through the genital tract. Here, we discovered that the oenocytes of fruitfly females are required for egg activation. Oenocytes, cells specialized in lipid-metabolism, are located beneath the abdominal cuticle. In adult flies, they synthesize the fatty acids (FAs) that are the precursors of cuticular hydrocarbons (CHCs), including pheromones. The oenocyte-targeted knockdown of a set of FA-anabolic enzymes, involved in very-long-chain fatty acid (VLCFA) synthesis, leads to a defect in egg activation. Given that some but not all of the identified enzymes are required for CHC/pheromone biogenesis, this putative VLCFA-dependent remote control may rely on an as-yet unidentified CHC or may function in parallel to CHC biogenesis. Additionally, we discovered that the most posterior ventral oenocyte cluster is in close proximity to the uterus. Since oocytes dissected from females deficient in this FA-anabolic pathway can be activated in vitro, this regulatory loop likely operates upstream of the calcium trigger. To our knowledge, our findings provide the first evidence that a physiological extra-genital signal remotely controls egg activation. Moreover, our study highlights a potential metabolic link between pheromone-mediated partner recognition and egg activation.

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.

Growth-inhibiting effects of the unconventional plant APYRASE 7 of Arabidopsis thaliana influences the LRX/RALF/FER growth regulatory module.

Plant cell growth involves coordination of numerous processes and signaling cascades among the different cellular compartments to concomitantly enlarge the protoplast and the surrounding cell wall. The cell wall integrity-sensing process involves the extracellular LRX (LRR-Extensin) proteins that bind RALF (Rapid ALkalinization Factor) peptide hormones and, in vegetative tissues, interact with the transmembrane receptor kinase FERONIA (FER). This LRX/RALF/FER signaling module influences cell wall composition and regulates cell growth. The numerous proteins involved in or influenced by this module are beginning to be characterized. In a genetic screen, mutations in Apyrase 7 (APY7) were identified to suppress growth defects observed in lrx1 and fer mutants. APY7 encodes a Golgi-localized NTP-diphosphohydrolase, but opposed to other apyrases of Arabidopsis, APY7 revealed to be a negative regulator of cell growth. APY7 modulates the growth-inhibiting effect of RALF1, influences the cell wall architecture and -composition, and alters the pH of the extracellular matrix, all of which affect cell growth. Together, this study reveals a function of APY7 in cell wall formation and cell growth that is connected to growth processes influenced by the LRX/RALF/FER signaling module.

Mutagenesis of Wheat Powdery Mildew Reveals a Single Gene Controlling Both NLR and Tandem Kinase-Mediated Immunity.

Blumeria graminis f. sp. tritici (Bgt) is a globally important fungal wheat pathogen. Some wheat genotypes contain powdery mildew resistance (Pm) genes encoding immune receptors that recognize specific fungal-secreted effector proteins, defined as avirulence (Avr) factors. Identifying Avr factors is vital for understanding the mechanisms, functioning, and durability of wheat resistance. Here, we present AvrXpose, an approach to identify Avr genes in Bgt by generating gain-of-virulence mutants on Pm genes. We first identified six Bgt mutants with gain of virulence on Pm3b and Pm3c. They all had point mutations, deletions or insertions of transposable elements within the corresponding AvrPm3b2/c2 gene or its promoter region. We further selected six mutants on Pm3a, aiming to identify the yet unknown AvrPm3a3 recognized by Pm3a, in addition to the previously described AvrPm3a2/f2. Surprisingly, Pm3a virulence in the obtained mutants was always accompanied by an additional gain of virulence on the unrelated tandem kinase resistance gene WTK4. No virulence toward 11 additional R genes tested was observed, indicating that the gain of virulence was specific for Pm3a and WTK4. Several independently obtained Pm3a-WTK4 mutants have mutations in Bgt-646, a gene encoding a putative, nonsecreted ankyrin repeat-containing protein. Gene expression analysis suggests that Bgt-646 regulates a subset of effector genes. We conclude that Bgt-646 is a common factor required for avirulence on both a specific nucleotide-binding leucine-rich repeat and a WTK immune receptor. Our findings suggest that, beyond effectors, another type of pathogen protein can control the race-specific interaction between powdery mildew and wheat. [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.

An online database for einkorn wheat to aid in gene discovery and functional genomics studies.

Diploid A-genome wheat (einkorn wheat) presents a nutrition-rich option as an ancient grain crop and a resource for the improvement of bread wheat against abiotic and biotic stresses. Realizing the importance of this wheat species, reference-level assemblies of two einkorn wheat accessions were generated (wild and domesticated). This work reports an einkorn genome database that provides an interface to the cereals research community to perform comparative genomics, applied genetics and breeding research. It features queries for annotated genes, the use of a recent genome browser release, and the ability to search for sequence alignments using a modern BLAST interface. Other features include a comparison of reference einkorn assemblies with other wheat cultivars through genomic synteny visualization and an alignment visualization tool for BLAST results. Altogether, this resource will help wheat research and breeding. Database URL  https://wheat.pw.usda.gov/GG3/pangenome.

Wheat zinc finger protein TaZF interacts with both the powdery mildew AvrPm2 protein and the corresponding wheat Pm2a immune receptor.

Plant defense responses to pathogens are induced after direct or indirect perception of effector proteins or their activity on host proteins. In fungal-plant interactions, relatively little is known about whether, in addition to avirulence effectors and immune receptors, other proteins contribute to specific recognition. The nucleotide-binding leucine-rich repeat (NLR) immune receptor Pm2a in wheat recognizes the fungal powdery mildew effector AvrPm2. We found that the predicted wheat zinc finger TaZF interacts with both the fungal avirulence protein AvrPm2 and the wheat NLR Pm2a. We further demonstrated that the virulent AvrPm2-H2 variant does not interact with TaZF. TaZF silencing in wheat resulted in a reduction but not a loss of Pm2a-mediated powdery mildew resistance. Interaction studies showed that the leucine-rich repeat domain of Pm2a is the mediator of the interaction with TaZF. TaZF recruits both Pm2a and AvrPm2 from the cytosol to the nucleus, resulting in nuclear localization of Pm2a, TaZF, and AvrPm2 in wheat. We propose that TaZF acts as a facilitator of Pm2a-dependent AvrPm2 effector recognition. Our findings highlight the importance of identifying effector host targets for characterization of NLR-mediated effector recognition.

Integration of genetic and genomics resources in einkorn wheat enables precision mapping of important traits.

Einkorn wheat (Triticum monococcum) is an ancient grain crop and a close relative of the diploid progenitor (T. urartu) of polyploid wheat. It is the only diploid wheat species having both domesticated and wild forms and therefore provides an excellent system to identify domestication genes and genes for traits of interest to utilize in wheat improvement. Here, we leverage genomic advancements for einkorn wheat using an einkorn reference genome assembly combined with skim-sequencing of a large genetic population of 812 recombinant inbred lines (RILs) developed from a cross between a wild and a domesticated T. monococcum accession. We identify 15,919 crossover breakpoints delimited to a median and average interval of 114 Kbp and 219 Kbp, respectively. This high-resolution mapping resource enables us to perform fine-scale mapping of one qualitative (red coleoptile) and one quantitative (spikelet number per spike) trait, resulting in the identification of small physical intervals (400 Kb to 700 Kb) with a limited number of candidate genes. Furthermore, an important domestication locus for brittle rachis is also identified on chromosome 7A. This resource presents an exciting route to perform trait discovery in diploid wheat for agronomically important traits and their further deployment in einkorn as well as tetraploid pasta wheat and hexaploid bread wheat cultivars.

Einkorn genomics sheds light on history of the oldest domesticated wheat.

Einkorn (Triticum monococcum) was the first domesticated wheat species, and was central to the birth of agriculture and the Neolithic Revolution in the Fertile Crescent around 10,000 years ago1,2. Here we generate and analyse 5.2-Gb genome assemblies for wild and domesticated einkorn, including completely assembled centromeres. Einkorn centromeres are highly dynamic, showing evidence of ancient and recent centromere shifts caused by structural rearrangements. Whole-genome sequencing analysis of a diversity panel uncovered the population structure and evolutionary history of einkorn, revealing complex patterns of hybridizations and introgressions after the dispersal of domesticated einkorn from the Fertile Crescent. We also show that around 1% of the modern bread wheat (Triticum aestivum) A subgenome originates from einkorn. These resources and findings highlight the history of einkorn evolution and provide a basis to accelerate the genomics-assisted improvement of einkorn and bread wheat.

A survey of lineage-specific genes in Triticeae reveals de novo gene evolution from genomic raw material.

Diploid plant genomes typically contain ~35,000 genes, almost all belonging to highly conserved gene families. Only a small fraction are lineage-specific, which are found in only one or few closely related species. Little is known about how genes arise de novo in plant genomes and how often this occurs; however, they are believed to be important for plants diversification and adaptation. We developed a pipeline to identify lineage-specific genes in Triticeae, using newly available genome assemblies of wheat, barley, and rye. Applying a set of stringent criteria, we identified 5942 candidate Triticeae-specific genes (TSGs), of which 2337 were validated as protein-coding genes in wheat. Differential gene expression analyses revealed that stress-induced wheat TSGs are strongly enriched in putative secreted proteins. Some were previously described to be involved in Triticeae non-host resistance and cold response. Additionally, we show that 1079 TSGs have sequence homology to transposable elements (TEs), ~68% of them deriving from regulatory non-coding regions of Gypsy retrotransposons. Most importantly, we demonstrate that these TSGs are enriched in transmembrane domains and are among the most highly expressed wheat genes overall. To summarize, we conclude that de novo gene formation is relatively rare and that Triticeae probably possess ~779 lineage-specific genes per haploid genome. TSGs, which respond to pathogen and environmental stresses, may be interesting candidates for future targeted resistance breeding in Triticeae. Finally, we propose that non-coding regions of TEs might provide important genetic raw material for the functional innovation of TM domains and the evolution of novel secreted proteins.

Ancient variation of the AvrPm17 gene in powdery mildew limits the effectiveness of the introgressed rye Pm17 resistance gene in wheat.

Introgressions of chromosomal segments from related species into wheat are important sources of resistance against fungal diseases. The durability and effectiveness of introgressed resistance genes upon agricultural deployment is highly variable-a phenomenon that remains poorly understood, as the corresponding fungal avirulence genes are largely unknown. Until its breakdown, the Pm17 resistance gene introgressed from rye to wheat provided broad resistance against powdery mildew (Blumeria graminis). Here, we used quantitative trait locus (QTL) mapping to identify the corresponding wheat mildew avirulence effector AvrPm17. It is encoded by two paralogous genes that exhibit signatures of reoccurring gene conversion events and are members of a mildew sublineage specific effector cluster. Extensive haplovariant mining in wheat mildew and related sublineages identified several ancient virulent AvrPm17 variants that were present as standing genetic variation in wheat powdery mildew prior to the Pm17 introgression, thereby paving the way for the rapid breakdown of the Pm17 resistance. QTL mapping in mildew identified a second genetic component likely corresponding to an additional resistance gene present on the 1AL.1RS translocation carrying Pm17. This gene remained previously undetected due to suppressed recombination within the introgressed rye chromosomal segment. We conclude that the initial effectiveness of 1AL.1RS was based on simultaneous introgression of two genetically linked resistance genes. Our results demonstrate the relevance of pathogen-based genetic approaches to disentangling complex resistance loci in wheat. We propose that identification and monitoring of avirulence gene diversity in pathogen populations become an integral part of introgression breeding to ensure effective and durable resistance in wheat.

Global genomic analyses of wheat powdery mildew reveal association of pathogen spread with historical human migration and trade.

The fungus Blumeria graminis f. sp. tritici causes wheat powdery mildew disease. Here, we study its spread and evolution by analyzing a global sample of 172 mildew genomes. Our analyses show that B.g. tritici emerged in the Fertile Crescent during wheat domestication. After it spread throughout Eurasia, colonization brought it to America, where it hybridized with unknown grass mildew species. Recent trade brought USA strains to Japan, and European strains to China. In both places, they hybridized with local ancestral strains. Thus, although mildew spreads by wind regionally, our results indicate that humans drove its global spread throughout history and that mildew rapidly evolved through hybridization.

Transposable Element Populations Shed Light on the Evolutionary History of Wheat and the Complex Co-Evolution of Autonomous and Non-Autonomous Retrotransposons.

Wheat has one of the largest and most repetitive genomes among major crop plants, containing over 85% transposable elements (TEs). TEs populate genomes much in the way that individuals populate ecosystems, diversifying into different lineages, sub-families and sub-populations. The recent availability of high-quality, chromosome-scale genome sequences from ten wheat lines enables a detailed analysis how TEs evolved in allohexaploid wheat, its diploids progenitors, and in various chromosomal haplotype segments. LTR retrotransposon families evolved into distinct sub-populations and sub-families that were active in waves lasting several hundred thousand years. Furthermore, It is shown that different retrotransposon sub-families were active in the three wheat sub-genomes, making them useful markers to study and date polyploidization events and chromosomal rearrangements. Additionally, haplotype-specific TE sub-families are used to characterize chromosomal introgressions in different wheat lines. Additionally, populations of non-autonomous TEs co-evolved over millions of years with their autonomous partners, leading to complex systems with multiple types of autonomous, semi-autonomous and non-autonomous elements. Phylogenetic and TE population analyses revealed the relationships between non-autonomous elements and their mobilizing autonomous partners. TE population analysis provided insights into genome evolution of allohexaploid wheat and genetic diversity of species, and may have implication for future crop breeding.

Comparative Transcriptome Analysis of Wheat Lines in the Field Reveals Multiple Essential Biochemical Pathways Suppressed by Obligate Pathogens.

Mildew and rust are the most devastating cereal pathogens, and in wheat they can cause up to 50% yield loss every year. Wheat lines containing resistance genes are used to effectively control fungal diseases, but the molecular mechanisms underlying the interaction between wheat and its fungal pathogens are poorly understood. Here, we used RNA sequencing (RNA-Seq) to compare the transcriptomic landscape of susceptible and resistant wheat lines to identify genes and pathways that are targeted by obligate biotrophic fungal pathogens. The five lines differed in the expression of thousands of genes under infection as well as control conditions. Generally, mixed infection with powdery mildew and leaf rust resulted in downregulation of numerous genes in susceptible lines. Interestingly, transcriptomic comparison between the nearly isogenic lines Thatcher and Thatcher-Lr34 identified 753 genes that are uniquely downregulated in the susceptible line upon infection. Kyoto encyclopedia of genes and genomes (KEGG) enrichment analysis, revealed the suppression of six major biochemical pathways, namely nuclear transport, alternative splicing, DNA damage response, ubiquitin-mediated proteolysis, phosphoinositol signaling, and photosynthesis. We conclude that powdery mildew and leaf rust evade the wheat defense system by suppression of programmed cell death (PCD) and responses to cellular damage. Considering the broad range of the induced changes, we propose that the pathogen targets "master regulators" at critical steps in the respective pathways. Identification of these wheat genes targeted by the pathogen could inspire new directions for future wheat breeding.

A population-level invasion by transposable elements triggers genome expansion in a fungal pathogen.

Genome evolution is driven by the activity of transposable elements (TEs). The spread of TEs can have deleterious effects including the destabilization of genome integrity and expansions. However, the precise triggers of genome expansions remain poorly understood because genome size evolution is typically investigated only among deeply divergent lineages. Here, we use a large population genomics dataset of 284 individuals from populations across the globe of Zymoseptoria tritici, a major fungal wheat pathogen. We built a robust map of genome-wide TE insertions and deletions to track a total of 2456 polymorphic loci within the species. We show that purifying selection substantially depressed TE frequencies in most populations, but some rare TEs have recently risen in frequency and likely confer benefits. We found that specific TE families have undergone a substantial genome-wide expansion from the pathogen's center of origin to more recently founded populations. The most dramatic increase in TE insertions occurred between a pair of North American populations collected in the same field at an interval of 25 years. We find that both genome-wide counts of TE insertions and genome size have increased with colonization bottlenecks. Hence, the demographic history likely played a major role in shaping genome evolution within the species. We show that both the activation of specific TEs and relaxed purifying selection underpin this incipient expansion of the genome. Our study establishes a model to recapitulate TE-driven genome evolution over deeper evolutionary timescales.

Defects in Cell Wall Differentiation of the Arabidopsis Mutant rol1-2 Is Dependent on Cyclin-Dependent Kinase CDK8.

Plant cells are encapsulated by cell walls whose properties largely determine cell growth. We have previously identified the rol1-2 mutant, which shows defects in seedling root and shoot development. rol1-2 is affected in the Rhamnose synthase 1 (RHM1) and shows alterations in the structures of Rhamnogalacturonan I (RG I) and RG II, two rhamnose-containing pectins. The data presented here shows that root tissue of the rol1-2 mutant fails to properly differentiate the cell wall in cell corners and accumulates excessive amounts of callose, both of which likely alter the physical properties of cells. A surr (suppressor of the rol1-2 root developmental defect) mutant was identified that alleviates the cell growth defects in rol1-2. The cell wall differentiation defect is re-established in the rol1-2 surr mutant and callose accumulation is reduced compared to rol1-2. The surr mutation is an allele of the cyclin-dependent kinase 8 (CDK8), which encodes a component of the mediator complex that influences processes central to plant growth and development. Together, the identification of the surr mutant suggests that changes in cell wall composition and turnover in the rol1-2 mutant have a significant impact on cell growth and reveals a function of CDK8 in cell wall architecture and composition.

Chromosome-scale genome assembly provides insights into rye biology, evolution and agronomic potential.

Rye (Secale cereale L.) is an exceptionally climate-resilient cereal crop, used extensively to produce improved wheat varieties via introgressive hybridization and possessing the entire repertoire of genes necessary to enable hybrid breeding. Rye is allogamous and only recently domesticated, thus giving cultivated ryes access to a diverse and exploitable wild gene pool. To further enhance the agronomic potential of rye, we produced a chromosome-scale annotated assembly of the 7.9-gigabase rye genome and extensively validated its quality by using a suite of molecular genetic resources. We demonstrate applications of this resource with a broad range of investigations. We present findings on cultivated rye's incomplete genetic isolation from wild relatives, mechanisms of genome structural evolution, pathogen resistance, low-temperature tolerance, fertility control systems for hybrid breeding and the yield benefits of rye-wheat introgressions.