Identification of a locus conferring dominant susceptibility to Pyrenophora tritici-repentis in barley

The fungus Pyrenophora tritici-repentis (Ptr) causes tan spot, a destructive foliar disease of wheat worldwide. The pathogen produces several necrotrophic effectors, which induce necrosis or chlorosis on susceptible wheat lines. Multiple races of Ptr have been identified, based on their ability to produce one or more of these effectors. Ptr has a wide host range of cereal and non-cereal grasses, but is known to cause damage only on wheat. Previously, we showed that Ptr can interact specifically with cultivated barley (Hordeum vulgare ssp. vulgare), and that the necrotrophic effector Ptr ToxB induces mild chlorosis in a highly selective manner when infiltrated into certain barley genotypes. In the present study, a barley doubled-haploid (DH) population was evaluated for reaction to Ptr race 5, a Ptr ToxB-producer. Then a comprehensive genetic map composed of 381 single nucleotide polymorphism (SNP) markers was used to map the locus conditioning this chlorosis. The F1 seedlings, and 92 DH lines derived from a cross between the resistant Japanese malting barley cultivar Haruna Nijo and the susceptible wild barley (H. vulgare ssp. spontaneum) OUH602 were inoculated with a conidial suspension of Ptr race 5 isolate at the two-leaf stage. The seedlings were monitored daily for symptoms and assessed for chlorosis development on the second leaf, 6 days after inoculation. All tested F1 seedlings exhibited chlorosis symptoms similar to the susceptible parent, and the DH lines segregated 1:1 for susceptible:resistant phenotypes, indicating the involvement of a single locus. Marker-trait linkage analysis based on interval mapping identified a single locus on the distal region of the short arm of chromosome 2H. We designate this locus Susceptibility to P. tritici-repentis1 (Spr1). The region encompassing this locus has 99 high confidence gene models, including membrane receptor-like kinases (RLKs), intracellular nucleotide-binding, leucine-rich repeat receptors (NLRs), and ankyrin-repeat proteins (ANKs). This shows the involvement of a dominant locus conferring susceptibility to Ptr in barley. Further work using high-resolution mapping and transgenic complementation will be required to identify the underlying gene

The pangenome of the wheat pathogen <i>Pyrenophora tritici-repentis</i> reveals novel transposons associated with necrotrophic effectors <i>ToxA</i> and <i>ToxB</i>

BACKGROUND: In fungal plant pathogens, genome rearrangements followed by selection pressure for adaptive traits have facilitated the co-evolutionary arms race between hosts and their pathogens. Pyrenophora tritici-repentis (Ptr) has emerged recently as a foliar pathogen of wheat worldwide and its populations consist of isolates that vary in their ability to produce combinations of different necrotrophic effectors. These effectors play vital roles in disease development. Here, we sequenced the genomes of a global collection (40 isolates) of Ptr to gain insights into its gene content and genome rearrangements. RESULTS: A comparative genome analysis revealed an open pangenome, with an abundance of accessory genes (~ 57%) reflecting Ptr’s adaptability. A clear distinction between pathogenic and non-pathogenic genomes was observed in size, gene content, and phylogenetic relatedness. Chromosomal rearrangements and structural organization, specifically around effector coding genes, were detailed using long-read assemblies (PacBio RS II) generated in this work in addition to previously assembled genomes. We also discovered the involvement of large mobile elements associated with Ptr’s effectors: ToxA, the gene encoding for the necrosis effector, was found as a single copy within a 143-kb ‘Starship’ transposon (dubbed ‘Horizon’) with a clearly defined target site and target site duplications. ‘Horizon’ was located on different chromosomes in different isolates, indicating mobility, and the previously described ToxhAT transposon (responsible for horizontal transfer of ToxA) was nested within this newly identified Starship. Additionally, ToxB, the gene encoding the chlorosis effector, was clustered as three copies on a 294-kb element, which is likely a different putative ‘Starship’ (dubbed ‘Icarus’) in a ToxB-producing isolate. ToxB and its putative transposon were missing from the ToxB non-coding reference isolate, but the homolog toxb and ‘Icarus’ were both present in a different non-coding isolate. This suggests that ToxB may have been mobile at some point during the evolution of the Ptr genome which is contradictory to the current assumption of ToxB vertical inheritance. Finally, the genome architecture of Ptr was defined as ‘one-compartment’ based on calculated gene distances and evolutionary rates. CONCLUSIONS: These findings together reflect on the highly plastic nature of the Ptr genome which has likely helped to drive its worldwide adaptation and has illuminated the involvement of giant transposons in facilitating the evolution of virulence in Ptr. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1186/s12915-022-01433-w

Genomic Prediction Accuracy of Stripe Rust in Six Spring Wheat Populations by Modeling Genotype by Environment Interaction

Some previous studies have assessed the predictive ability of genome-wide selection on stripe (yellow) rust resistance in wheat, but the effect of genotype by environment interaction (GEI) in prediction accuracies has not been well studied in diverse genetic backgrounds. Here, we compared the predictive ability of a model based on phenotypic data only (M1), the main effect of phenotype and molecular markers (M2), and a model that incorporated GEI (M3) using three cross-validations (CV1, CV2, and CV0) scenarios of interest to breeders in six spring wheat populations. Each population was evaluated at three to eight field nurseries and genotyped with either the DArTseq technology or the wheat 90K single nucleotide polymorphism arrays, of which a subset of 1,058- 23,795 polymorphic markers were used for the analyses. In the CV1 scenario, the mean prediction accuracies of the M1, M2, and M3 models across the six populations varied from 0.11 to 0.07, from 0.22 to 0.49, and from 0.19 to 0.48, respectively. Mean accuracies obtained using the M3 model in the CV1 scenario were significantly greater than the M2 model in two populations, the same in three populations, and smaller in one population. In both the CV2 and CV0 scenarios, the mean prediction accuracies of the three models varied from 0.53 to 0.84 and were not significantly different in all populations, except the Attila/CDC Go in the CV2, where the M3 model gave greater accuracy than both the M1 and M2 models. Overall, the M3 model increased prediction accuracies in some populations by up to 12.4% and decreased accuracy in others by up to 17.4%, demonstrating inconsistent results among genetic backgrounds that require considering each population separately. This is the first comprehensive genome-wide prediction study that investigated details of the effect of GEI on stripe rust resistance across diverse spring wheat populations

Identification of Disease Resistance Parents and Genome-Wide Association Mapping of Resistance in Spring Wheat

The likelihood of success in developing modern cultivars depend on multiple factors, including the identification of suitable parents to initiate new crosses, and characterizations of genomic regions associated with target traits. The objectives of the present study were to (a) determine the best economic weights of four major wheat diseases (leaf spot, common bunt, leaf rust, and stripe rust) and grain yield for multi-trait restrictive linear phenotypic selection index (RLPSI), (b) select the top 10% cultivars and lines (hereafter referred as genotypes) with better resistance to combinations of the four diseases and acceptable grain yield as potential parents, and (c) map genomic regions associated with resistance to each disease using genome-wide association study (GWAS). A diversity panel of 196 spring wheat genotypes was evaluated for their reaction to stripe rust at eight environments, leaf rust at four environments, leaf spot at three environments, common bunt at two environments, and grain yield at five environments. The panel was genotyped with the Wheat 90K SNP array and a few KASP SNPs of which we used 23,342 markers for statistical analyses. The RLPSI analysis performed by restricting the expected genetic gain for yield displayed significant (p \u3c 0.05) differences among the 3125 economic weights. Using the best four economic weights, a subset of 22 of the 196 genotypes were selected as potential parents with resistance to the four diseases and acceptable grain yield. GWAS identified 37 genomic regions, which included 12 for common bunt, 13 for leaf rust, 5 for stripe rust, and 7 for leaf spot. Each genomic region explained from 6.6 to 16.9% and together accounted for 39.4% of the stripe rust, 49.1% of the leaf spot, 94.0% of the leaf rust, and 97.9% of the common bunt phenotypic variance combined across all environments. Results from this study provide valuable information for wheat breeders selecting parental combinations for new crosses to develop improved germplasm with enhanced resistance to the four diseases as well as the physical positions of genomic regions that confer resistance, which facilitates direct comparisons for independent mapping studies in the future

Evaluating the importance of the tan spot ToxA–Tsn1 interaction in Australian wheat varieties

The necrotrophic fungal pathogen Pyrenophora tritici-repentis (Ptr) causes the major wheat disease tan spot, and produces multiple necrotrophic effectors that contribute to virulence. The proteinaceous effector ToxA induces necrosis in wheat genotypes possessing the Tsn1 gene, although the importance of the ToxA–Tsn1 interaction itself in varietal disease development has not been well studied. Here, 40 Australian spring wheat varieties were assessed for ToxA sensitivity and disease response to a race 1 wildtype Ptr isolate and ToxA-deleted strain at both seedling and tillering growth stages. ToxA sensitivity was generally associated with disease susceptibility, but did not always predict spreading necrotic symptoms. Whilst the majority of Tsn1 varieties exhibited lower disease scores following toxa mutant infection, several exhibited no distinct differences between wildtype and toxa symptoms. This implies that ToxA is not the major determinant in tan spot disease development in some host backgrounds and indicates the presence of additional effectors. Unexpectedly, several tsn1 varieties exhibited a reduction in disease severity following toxa mutant inoculation, which may suggest an indirect role for ToxA in pathogen fitness. Additionally, increased chlorosis was observed following toxa mutant infection in three varieties, and further work is required to determine whether this is likely to be due to ToxA epistasis of ToxC symptoms. Taken together, these observations demonstrate that Ptr interacts with the host in a complex and intricate manner, leading to a variety of disease reactions that are dependent or independent of the ToxA–Tsn1 interaction

Comparative genomics of a plant-pathogenic fungus, pyrenophora tritici-repentis, reveals transduplication and the impact of repeat elements on pathogenicity and population divergence

Pyrenophora tritici-repentis is a necrotrophic fungus causal to the disease tan spot of wheat, whose contribution to crop loss has increased significantly during the last few decades. Pathogenicity by this fungus is attributed to the production of host-selective toxins (HST), which are recognized by their host in a genotype-specific manner. To better understand the mechanisms that have led to the increase in disease incidence related to this pathogen, we sequenced the genomes of three P. tritici-repentis isolates. A pathogenic isolate that produces two known HSTs was used to assemble a reference nuclear genome of approximately 40 Mb composed of 11 chromosomes that encode 12,141 predicted genes. Comparison of the reference genome with those of a pathogenic isolate that produces a third HST, and a nonpathogenic isolate, showed the nonpathogen genome to be more diverged than those of the two pathogens. Examination of gene-coding regions has provided candidate pathogen-specific proteins and revealed gene families that may play a role in a necrotrophic lifestyle. Analysis of transposable elements suggests that their presence in the genome of pathogenic isolates contributes to the creation of novel genes, effector diversification, possible horizontal gene transfer events, identified copy number variation, and the first example of transduplication by DNA transposable elements in fungi.Overall, comparative analysis of these genomes provides evidence that pathogenicity in this species arose through an influx of transposable elements, which created a genetically flexible landscape that can easily respond to environmental changes

Causes and Consequences of Genome Expansion in Fungi

Fungi display a large diversity in genome size and complexity, variation that is often considered to be adaptive. But because nonadaptive processes can also have important consequences on the features of genomes, we investigated the relationship of genetic drift and genome size in the phylum Ascomycota using multiple indicators of genetic drift. We detected a complex relationship between genetic drift and genome size in fungi: genetic drift is associated with genome expansion on broad evolutionary timescales, as hypothesized for other eukaryotes; but within subphyla over smaller timescales, the opposite trend is observed. Moreover, fungi and bacteria display similar patterns of genome degradation that are associated with initial effects of genetic drift. We conclude that changes in genome size within Ascomycota have occurred using two different routes: large-scale genome expansions are catalyzed by increasing drift as predicted by the mutation-hazard model of genome evolution and small-scale modifications in genome size are independent of drift

Multiple Translocation of the AVR-Pita Effector Gene among Chromosomes of the Rice Blast Fungus Magnaporthe oryzae and Related Species

Magnaporthe oryzae is the causal agent of rice blast disease, a devastating problem worldwide. This fungus has caused breakdown of resistance conferred by newly developed commercial cultivars. To address how the rice blast fungus adapts itself to new resistance genes so quickly, we examined chromosomal locations of AVR-Pita, a subtelomeric gene family corresponding to the Pita resistance gene, in various isolates of M. oryzae (including wheat and millet pathogens) and its related species. We found that AVR-Pita (AVR-Pita1 and AVR-Pita2) is highly variable in its genome location, occurring in chromosomes 1, 3, 4, 5, 6, 7, and supernumerary chromosomes, particularly in rice-infecting isolates. When expressed in M. oryzae, most of the AVR-Pita homologs could elicit Pita-mediated resistance, even those from non-rice isolates. AVR-Pita was flanked by a retrotransposon, which presumably contributed to its multiple translocation across the genome. On the other hand, family member AVR-Pita3, which lacks avirulence activity, was stably located on chromosome 7 in a vast majority of isolates. These results suggest that the diversification in genome location of AVR-Pita in the rice isolates is a consequence of recognition by Pita in rice. We propose a model that the multiple translocation of AVR-Pita may be associated with its frequent loss and recovery mediated by its transfer among individuals in asexual populations. This model implies that the high mobility of AVR-Pita is a key mechanism accounting for the rapid adaptation toward Pita. Dynamic adaptation of some fungal plant pathogens may be achieved by deletion and recovery of avirulence genes using a population as a unit of adaptation

Influence of nitrogen sources on growth and mycotoxin production by isolates of Pyrenophora tritici-repentis from wheat

The fungus Pyrenophora tritici-repentis (Died.) Drechs. infects the leaves and kernels of wheat, causing tan spot and red smudge, respectively. Isolates of P. tritici-repentis have been reported to be both phytopathogenic and mycotoxigenic. This research investigates the influence of nitrogen sources on growth and production of mycotoxins by eight different isolates of P. tritici-repentis. A synthetic agar medium (SAM) was used with different nitrogen sources, both inorganic [(NH4Cl, NH4NO3 and (NH4)2SO4)] and organic (l-alanine, l-histidine, and l-lysine), at a concentration of 37.5 mmol L− 1. Individual isolates exhibited different growth rates that varied according to the nitrogen source added to the medium. The choice of nitrogen source also had a major effect on production of the mycotoxins emodin, catenarin and islandicin. The highest concentrations of emodin, 54.40 ± 4.46 μg g− 1, 43.07 ± 23.39 μg g− 1 and 28.91 ± 4.64 μg g− 1 of growth medium, were produced on the complex medium (V8-potato dextrose agar) by the isolates Alg-H2, 331-2 and TS93-71B, respectively. A relatively high concentration of emodin also was produced by isolates Az35-5 (28.29 ± 4.71 μg g− 1 of medium) and TS93-71B (27.03 ± 4.09 μg g− 1 of medium) on synthetic medium supplemented with l-alanine. The highest concentrations of catenarin (174.54 ± 14.46 μg g− 1 and 104.87 ± 6.13 μg g− 1 of medium) were recorded for isolates TS93-71B and Alg-H2 on synthetic medium supplemented with l-alanine and NH4Cl, respectively. The highest concentration of islandicin (4.64 ± 0.36 μg g− 1 medium) was observed for isolate 331-2 in the presence of l-lysine. There was not a close relationship between mycelial growth and mycotoxin production by the fungal isolates. This is the first report on the influence of nitrogen sources on the production of mycotoxins by P. tritici-repentis