Assembly and evaluation of a confocal microscopy image analysis pipeline useful in revealing the secrets of plant-fungal interactions.

The ability of laser scanning confocal microscopy to generate high-contrast 2D and 3D images has become essential in studying plant-fungal interactions. Techniques such as visualization of native fluorescence, fluorescent protein tagging of microbes, GFP/RFP-fusion proteins, and fluorescent labelling of plant and fungal proteins have been widely used to aid in these investigations. Use of fluorescent proteins has several pitfalls including variability of expression in planta and the requirement of gene transformation. Here we used the unlabeled pathogens Parastagonospora nodorum, Pyrenophora teres f. teres, and Cercospora beticola infecting wheat, barley, and sugar beet respectively, to show the utility of a staining and imaging pipeline that uses propidium iodide (PI), which stains RNA and DNA, and wheat germ agglutinin labeled with fluorescein isothiocyanate (WGA-FITC), which stains chitin, to visualize fungal colonization of plants. This pipeline relies on the use of KOH to remove the cutin layer of the leaf, increasing its permeability, allowing the different stains to penetrate and effectively bind to their targets, resulting in a consistent visualization of cellular structures. To expand the utility of this pipeline, we used the staining techniques in conjunction with machine learning to analyze fungal biomass through volume analysis, as well as quantifying nuclear breakdown, an early indicator of programmed cell death (PCD). This pipeline is simple to use, robust, consistent across host and fungal species and can be applied to most plant-fungal interactions. Therefore, this pipeline can be used to characterize model systems as well as non-model interactions where transformation is not routine.

Variability in Chromosome 1 of Select Moroccan Pyrenophora teres f. teres Isolates Overcomes a Highly Effective Barley Chromosome 6H Source of Resistance.

Barley net form net blotch (NFNB) is a destructive foliar disease caused by Pyrenophora teres f. teres. Barley line CIho5791, which harbors the broadly effective chromosome 6H resistance gene Rpt5, displays dominant resistance to P. teres f. teres. To genetically characterize P. teres f. teres avirulence/virulence on the barley line CIho5791, we generated a P. teres f. teres mapping population using a cross between the Moroccan CIho5791-virulent isolate MorSM40-3 and the avirulent reference isolate 0-1. Full genome sequences were generated for 103 progenies. Saturated chromosome-level genetic maps were generated, and quantitative trait locus (QTL) mapping identified two major QTL associated with P. teres f. teres avirulence/virulence on CIho5791. The most significant QTL mapped to chromosome (Ch) 1, where the virulent allele was contributed by MorSM40-3. A second QTL mapped to Ch8; however, this virulent allele was contributed by the avirulent parent 0-1. The Ch1 and Ch8 loci accounted for 27 and 15% of the disease variation, respectively, and the avirulent allele at the Ch1 locus was epistatic over the virulent allele at the Ch8 locus. As a validation, we used a natural P. teres f. teres population in a genome-wide association study that identified the same Ch1 and Ch8 loci. We then generated a new reference quality genome assembly of parental isolate MorSM40-3 with annotation supported by deep transcriptome sequencing of infection time points. The annotation identified candidate genes predicted to encode small, secreted proteins, one or more of which are likely responsible for overcoming the CIho5791 resistance. [Formula: see text] The author(s) have dedicated the work to the public domain under the Creative Commons CC0 "No Rights Reserved" license by waiving all of his or her rights to the work worldwide under copyright law, including all related and neighboring rights, to the extent allowed by law, 2024.

A Moroccan Pyrenophora teres f. teres Population Defeats Rpt5, the Broadly Effective Resistance on Barley Chromosome 6H.

Net form net blotch (NFNB), caused by Pyrenophora teres f. teres, is an important barley disease. The centromeric region of barley chromosome 6H has often been associated with resistance or susceptibility to NFNB, including the broadly effective dominant resistance gene Rpt5 derived from barley line CIho 5791. We characterized a population of Moroccan P. teres f. teres isolates that had overcome Rpt5 resistance and identified quantitative trait loci (QTL) that were effective against these isolates. Eight Moroccan P. teres f. teres isolates were phenotyped on barley lines CIho 5791 and Tifang. Six isolates were virulent on CIho 5791, and two were avirulent. A CIho 5791 × Tifang recombinant inbred line (RIL) population was phenotyped with all eight isolates and confirmed the defeat of the 6H resistance locus formerly mapped as Rpt5 in barley line CI9819. A major QTL on chromosome 3H with the resistance allele derived from Tifang, as well as minor QTL, was identified and provided resistance against these isolates. F2 segregation ratios supported dominant inheritance for both the 3H and 6H resistance. Furthermore, inoculation of progeny isolates derived from a cross of P. teres f. teres isolates 0-1 (virulent on Tifang/avirulent on CIho 5791) and MorSM 40-3 (avirulent on Tifang/virulent on CIho 5791) onto the RIL and F2 populations determined that recombination between isolates can generate novel genotypes that overcome both resistance genes. Markers linked to the QTL identified in this study can be used to incorporate both resistance loci into elite barley cultivars for durable resistance.

A triple threat: the Parastagonospora nodorum SnTox267 effector exploits three distinct host genetic factors to cause disease in wheat.

Parastagonospora nodorum is a fungal pathogen of wheat. As a necrotrophic specialist, it deploys effector proteins that target dominant host susceptibility genes to elicit programmed cell death (PCD). Here we identify and functionally validate the effector targeting the host susceptibility genes Snn2, Snn6 and Snn7. We utilized whole-genome sequencing, association mapping, gene-disrupted mutants, gain-of-function transformants, virulence assays, bioinformatics and quantitative PCR to characterize these interactions. A single proteinaceous effector, SnTox267, targeted Snn2, Snn6 and Snn7 to trigger PCD. Snn2 and Snn6 functioned cooperatively to trigger PCD in a light-dependent pathway, whereas Snn7-mediated PCD functioned in a light-independent pathway. Isolates harboring 20 SnTox267 protein isoforms quantitatively varied in virulence. The diversity and distribution of isoforms varied between populations, indicating adaptation to local selection pressures. SnTox267 deletion resulted in the upregulation of effector genes SnToxA, SnTox1 and SnTox3. We validated a novel effector operating in an inverse-gene-for-gene manner to target three genetically distinct host susceptibility genes and elicit PCD. The discovery of the complementary gene action of Snn2 and Snn6 indicates their potential function in a guard or decoy model. Additionally, differences in light dependency in the elicited pathways and upregulation of unlinked effectors sheds new light onto a complex fungal necrotroph-host interaction.

The Parastagonospora nodorum necrotrophic effector SnTox5 targets the wheat gene Snn5 and facilitates entry into the leaf mesophyll.

Parastagonospora nodorum is an economically important necrotrophic fungal pathogen of wheat. Parastagonospora nodorum secretes necrotrophic effectors that target wheat susceptibility genes to induce programmed cell death (PCD). In this study, we cloned and functionally validated SnTox5 and characterized its role in pathogenesis. We used whole genome sequencing, genome-wide association study (GWAS) mapping, CRISPR-Cas9-based gene disruption, gain-of-function transformation, quantitative trait locus (QTL) analysis, haplotype and isoform analysis, protein modeling, quantitative PCR, and laser confocal microscopy to validate SnTox5 and functionally characterize SnTox5. SnTox5 is a mature 16.26 kDa protein with high structural similarity to SnTox3. Wild-type and mutant P. nodorum strains and wheat genotypes of SnTox5 and Snn5, respectively, were used to show that SnTox5 not only targets Snn5 to induce PCD but also facilitates the colonization of the mesophyll layer even in the absence of Snn5. Here we show that SnTox5 facilitates the efficient colonization of the mesophyll tissue and elicits PCD specific to host lines carrying Snn5. The homology to SnTox3 and the ability of SnTox5 to facilitate the colonizing of the mesophyll also suggest a role in the suppression of host defense before PCD induction.

Host and pathogen genetics reveal an inverse gene-for-gene association in the P. teres f. maculata-barley pathosystem.

KEY MESSAGE: Pathogen and host genetics were used to uncover an inverse gene-for-gene interaction where virulence genes from the pathogen Pyrenophora teres f. maculata target barley susceptibility genes, resulting in disease. Although models have been proposed to broadly explain how plants and pathogens interact and coevolve, each interaction evolves independently, resulting in various scenarios of host manipulation and plant defense. Spot form net blotch is a foliar disease of barley caused by Pyrenophora teres f. maculata. We developed a barley population (Hockett × PI 67381) segregating for resistance to a diverse set of P. teres f. maculata isolates. Quantitative trait locus analysis identified major loci on barley chromosomes (Chr) 2H and 7H associated with resistance/susceptibility. Subsequently, we used avirulent and virulent P. teres f. maculata isolates to develop a pathogen population, identifying two major virulence loci located on Chr1 and Chr2. To further characterize this host-pathogen interaction, progeny from the pathogen population harboring virulence alleles at either the Chr1 or Chr2 locus was phenotyped on the Hockett × PI 67381 population. Progeny harboring only the Chr1 virulence allele lost the barley Chr7H association but maintained the 2H association. Conversely, isolates harboring only the Chr2 virulence allele lost the barley Chr2H association but maintained the 7H association. Hockett × PI 67381 F2 individuals showed susceptible/resistant ratios not significantly different than 15:1 and results from F2 inoculations using the single virulence genotypes were not significantly different from a 3:1 (S:R) ratio, indicating two dominant susceptibility genes. Collectively, this work shows that P. teres f. maculata virulence alleles at the Chr1 and Chr2 loci are targeting the barley 2H and 7H susceptibility alleles in an inverse gene-for-gene manner to facilitate colonization.

Four Reference Quality Genome Assemblies of Pyrenophora teres f. maculata: A Resource for Studying the Barley Spot Form Net Blotch Interaction.

Pyrenophora teres is the causal agent of net blotch, the most devastating foliar disease of barley. In nature, net blotch is seen in two forms, net form net blotch, caused by P. teres f. teres, and spot form net blotch, caused by P. teres f. maculata. To date, 11 P. teres f. teres genomes have been sequenced and deposited in publicly available repositories, but only one P. teres f. maculata genome has been publicly deposited. Here, we present four additional reference-quality full-genome sequences of P. teres f. maculata isolates with good geographical and phenotypic diversity, with accompanying RNA sequencing-based genome annotations. These additional P. teres f. maculata genomes will aid in the understanding of the genomic complexities of this important barley pathogen.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

A Comparative Genomic Analysis of the Barley Pathogen Pyrenophora teres f. teres Identifies Subtelomeric Regions as Drivers of Virulence.

Pyrenophora teres f. teres causes net form net blotch of barley and is an economically important pathogen throughout the world. However, P. teres f. teres is lacking in the genomic resources necessary to characterize the mechanisms of virulence. Recently a high-quality reference genome was generated for P. teres f. teres isolate 0-1. Here, we present the reference quality sequence and annotation of four new isolates and we use the five available P. teres f. teres genomes for an in-depth comparison, resulting in the generation of hypotheses pertaining to the potential mechanisms and evolution of virulence. Comparative analyses were performed between all five P. teres f. teres genomes, examining genomic organization, structural variations, and core and accessory genomic content, specifically focusing on the genomic characterization of known virulence loci and the localization of genes predicted to encode secreted and effector proteins. We showed that 14 of 15 currently published virulence quantitative trait loci (QTL) span accessory genomic regions, consistent with these accessory regions being important drivers of host adaptation. Additionally, these accessory genomic regions were frequently found in subtelomeric regions of chromosomes, with 10 of the 14 accessory region QTL localizing to subtelomeric regions. Comparative analysis of the subtelomeric regions of P. teres f. teres chromosomes revealed translocation events in which homology was detected between nonhomologous chromosomes at a significantly higher rate than the rest of the genome. These results indicate that the subtelomeric accessory genomic compartments not only harbor most of the known virulence loci but, also, that these regions have the capacity to rapidly evolve.

Genetic analysis of virulence in the Pyrenophora teres f. teres population BB25×FGOH04Ptt-21.

Pyrenophora teres f. teres is the causal agent of net form net blotch (NFNB) of barley. In order to map the genetics of avirulence/virulence in P. teres f. teres, a fungal population was developed using P. teres f. teres isolates BB25 (Denmark) and FGOH04Ptt-21 (North Dakota, USA) due to these two isolates differing in virulence on several common barley lines. 109 progeny isolates were obtained from the BB25 by FGOH04Ptt-21 cross that were then used for NFNB disease evaluation across eight barley lines, four of which have been used commonly as NFNB differential lines as well as four cultivars commonly used in barley production in the Northern Great Plains. BB25 was virulent on one of the barley lines and avirulent on seven of the barley lines whereas, FGOH04Ptt-21 was virulent on all eight barley lines evaluated. Genetic maps were generated with single nucleotide polymorphism (SNP) markers obtained using a restriction associated DNA genotyping by sequencing (RAD-GBS) approach. Sixteen linkage groups were formed and were used to identify quantitative trait loci (QTL) associated with avirulence/virulence. Nine unique QTL were identified on eight linkage groups out of which three QTL had major effects (R2≥45%) while the remaining six QTL were relatively minor (R2<20%). One or two major effect loci were identified for the lines commonly used as differentials. Conversely, variation in virulence on the local barley cultivars was mostly associated with small effect loci that contributed quantitatively to disease.

Genetic drift, historic migration, and limited gene flow contributing to the subpopulation divergence in wild sea beet (Beta vulgaris ssp. maritima (L.) Arcang).

Cultivated beet (Beta vulgaris L. ssp. vulgaris) originated from sea beet (B. vulgaris ssp. maritima (L.) Arcang), a wild beet species widely distributed along the coasts of the Mediterranean Sea and Atlantic Ocean, as well as northern Africa. Understanding the evolution of sea beet will facilitate its efficient use in sugarbeet improvement. We used SNPs (single nucleotide polymorphisms) covering the whole genome to analyze 599 sea beet accessions collected from the north Atlantic Ocean and Mediterranean Sea coasts. All B. maritima accessions can be grouped into eight clusters with each corresponding to a specific geographic region. Clusters 2, 3 and 4 with accessions mainly collected from Mediterranean coasts are genetically close to each other as well as to Cluster 6 that contained mainly cultivated beet. Other clusters were relatively distinct from cultivated beets with Clusters 1 and 5 containing accessions from north Atlantic Ocean coasts, Clusters 7 and Cluster 8 mainly have accessions from northern Egypt and southern Europe, and northwest Morocco, respectively. Distribution of B. maritima subpopulations aligns well with the direction of marine currents that was considered a main dynamic force in spreading B. maritima during evolution. Estimation of genetic diversity indices supported the formation of B. maritima subpopulations due to local genetic drift, historic migration, and limited gene flow. Our results indicated that B. maritima originated from southern Europe and then spread to other regions through marine currents to form subpopulations. This research provides vital information for conserving, collecting, and utilizing wild sea beet to sustain sugarbeet improvement.

Meta-transcriptomic analysis reveals the geographical expansion of known sugarbeet-infecting viruses and the occurrence of a novel virus in sugarbeet in the United States.

In this study, meta-transcriptome sequencing was conducted on a total of 18 sugarbeet (Beta vulgaris L. subsp. vulgaris) sample libraries to profile the virome of field-grown sugarbeet to identify the occurrence and distribution of known and potentially new viruses from five different states in the United States. Sugarbeet roots with symptoms resembling rhizomania caused by beet necrotic yellow vein virus (BNYVV), or leaves exhibiting leaf-curling, yellowing to browning, or green mosaic were collected from the sugarbeet growing areas of California, Colorado, Idaho, Minnesota, and North Dakota. In silico analysis of de novo assembled contigs revealed the presence of nearly full-length genomes of BNYVV, beet soil-borne virus (BSBV), and beet soil-borne mosaic virus (BSBMV), which represent known sugarbeet-infecting viruses. Among those, BNYVV was widespread across the locations, whereas BSBV was prevalent in Minnesota and Idaho, and BSBMV was only detected in Minnesota. In addition, two recently reported Beta vulgaris satellite virus isoforms (BvSatV-1A and BvSatV-1B) were detected in new locations, indicating the geographical expansion of this known virus. Besides these known sugarbeet-infecting viruses, the bioinformatic analysis identified the widespread occurrence of a new uncharacterized Erysiphe necator-associated abispo virus (En_abispoV), a fungus-related virus that was identified in all 14 libraries. En_abispoV contains two RNA components, and nearly complete sequences of both RNA1 and RNA2 were obtained from RNASeq and were further confirmed by primer-walking RT-PCR and Sanger sequencing. Phylogenetic comparison of En_abispoV isolates obtained in this study showed varying levels of genetic diversity within RNA1 and RNA2 compared to previously reported isolates. The undertaken meta-transcriptomic approach revealed the widespread nature of coexisting viruses associated with field-grown sugarbeet exhibiting virus disease-like symptoms in the United States.

Hybrid inferiority and genetic incompatibilities drive divergence of fungal pathogens infecting the same host.

Agro-ecosystems provide environments that are conducive for rapid evolution and dispersal of plant pathogens. Previous studies have demonstrated that hybridization of crop pathogens can give rise to new lineages with altered virulence profiles. Currently, little is known about either the genetics of fungal pathogen hybridization or the mechanisms that may prevent hybridization between related species. The fungus Pyrenophora teres is a global pathogen of barley. The pathogenic fungus P. teres exists as two distinct lineages P. teres f. teres and P. teres f. maculata (Ptt and Ptm, respectively), which both infect barley but produce very distinct lesions and rarely interbreed. Interestingly, Ptt and Ptm can, by experimental mating, produce viable progenies. Here, we addressed the underlying genetics of reproductive barriers of P. teres. We hypothesize that Ptt and Ptm diverged in the past, possibly by adapting to distinct hosts, and only more recently colonized the same host in agricultural fields. Using experimental mating and in planta phenotyping in barley cultivars susceptible to both P. teres forms, we demonstrate that hybrids produce mixed infection phenotypes but overall show inferior pathogenic fitness relative to the pure parents. Based on analyses of 104 hybrid genomes, we identify signatures of negative epistasis between parental alleles at distinct loci (Dobzhansky-Müller incompatibilities). Most DMI regions are not involved in virulence but certain genes are predicted or known to play a role in virulence. These results potentially suggest that divergent niche adaptation-albeit in the same host plant-contributes to speciation in P. teres.

Research advances in the Pyrenophora teres-barley interaction.

Pyrenophora teres f. teres and P. teres f. maculata are significant pathogens that cause net blotch of barley. An increased number of loci involved in P. teres resistance or susceptibility responses of barley as well as interacting P. teres virulence effector loci have recently been identified through biparental and association mapping studies of both the pathogen and host. Characterization of the resistance/susceptibility loci in the host and the interacting effector loci in the pathogen will provide a path for targeted gene validation for better-informed release of resistant barley cultivars. This review assembles concise consensus maps for all loci published for both the host and pathogen, providing a useful resource for the community to be used in pathogen characterization and barley breeding for resistance to both forms of P. teres.

Reference Assembly and Annotation of the Pyrenophora teres f. teres Isolate 0-1.

Pyrenophora teres f. teres, the causal agent of net form net blotch (NFNB) of barley, is a destructive pathogen in barley-growing regions throughout the world. Typical yield losses due to NFNB range from 10 to 40%; however, complete loss has been observed on highly susceptible barley lines where environmental conditions favor the pathogen. Currently, genomic resources for this economically important pathogen are limited to a fragmented draft genome assembly and annotation, with limited RNA support of the P. teres f. teres isolate 0-1. This research presents an updated 0-1 reference assembly facilitated by long-read sequencing and scaffolding with the assistance of genetic linkage maps. Additionally, genome annotation was mediated by RNAseq analysis using three infection time points and a pure culture sample, resulting in 11,541 high-confidence gene models. The 0-1 genome assembly and annotation presented here now contains the majority of the repetitive content of the genome. Analysis of the 0-1 genome revealed classic characteristics of a "two-speed" genome, being compartmentalized into GC-equilibrated and AT-rich compartments. The assembly of repetitive AT-rich regions will be important for future investigation of genes known as effectors, which often reside in close proximity to repetitive regions. These effectors are responsible for manipulation of the host defense during infection. This updated P. teres f. teres isolate 0-1 reference genome assembly and annotation provides a robust resource for the examination of the barley-P. teres f. teres host-pathogen coevolution.

Reference Quality Genome Assemblies of Three Parastagonospora nodorum Isolates Differing in Virulence on Wheat.

Parastagonospora nodorum, the causal agent of Septoria nodorum blotch in wheat, has emerged as a model necrotrophic fungal organism for the study of host-microbe interactions. To date, three necrotrophic effectors have been identified and characterized from this pathogen, including SnToxA, SnTox1, and SnTox3. Necrotrophic effector identification was greatly aided by the development of a draft genome of Australian isolate SN15 via Sanger sequencing, yet it remained largely fragmented. This research presents the development of nearly finished genomes of P. nodorum isolates Sn4, Sn2000, and Sn79-1087 using long-read sequencing technology. RNAseq analysis of isolate Sn4, consisting of eight time points covering various developmental and infection stages, mediated the annotation of 13,379 genes. Analysis of these genomes revealed large-scale polymorphism between the three isolates, including the complete absence of contig 23 from isolate Sn79-1087, and a region of genome expansion on contig 10 in isolates Sn4 and Sn2000. Additionally, these genomes exhibit the hallmark characteristics of a "two-speed" genome, being partitioned into two distinct GC-equilibrated and AT-rich compartments. Interestingly, isolate Sn79-1087 contains a lower proportion of AT-rich segments, indicating a potential lack of evolutionary hotspots. These newly sequenced genomes, consisting of telomere-to-telomere assemblies of nearly all 23 P. nodorum chromosomes, provide a robust foundation for the further examination of effector biology and genome evolution.

Transposable Element Genomic Fissuring in Pyrenophora teres Is Associated With Genome Expansion and Dynamics of Host-Pathogen Genetic Interactions.

Pyrenophora teres, P. teres f. teres (PTT) and P. teres f. maculata (PTM) cause significant diseases in barley, but little is known about the large-scale genomic differences that may distinguish the two forms. Comprehensive genome assemblies were constructed from long DNA reads, optical and genetic maps. As repeat masking in fungal genomes influences the final gene annotations, an accurate and reproducible pipeline was developed to ensure comparability between isolates. The genomes of the two forms are highly collinear, each composed of 12 chromosomes. Genome evolution in P. teres is characterized by genome fissuring through the insertion and expansion of transposable elements (TEs), a process that isolates blocks of genic sequence. The phenomenon is particularly pronounced in PTT, which has a larger, more repetitive genome than PTM and more recent transposon activity measured by the frequency and size of genome fissures. PTT has a longer cultivated host association and, notably, a greater range of host-pathogen genetic interactions compared to other Pyrenophora spp., a property which associates better with genome size than pathogen lifestyle. The two forms possess similar complements of TE families with Tc1/Mariner and LINE-like Tad-1 elements more abundant in PTT. Tad-1 was only detectable as vestigial fragments in PTM and, within the forms, differences in genome sizes and the presence and absence of several TE families indicated recent lineage invasions. Gene differences between P. teres forms are mainly associated with gene-sparse regions near or within TE-rich regions, with many genes possessing characteristics of fungal effectors. Instances of gene interruption by transposons resulting in pseudogenization were detected in PTT. In addition, both forms have a large complement of secondary metabolite gene clusters indicating significant capacity to produce an array of different molecules. This study provides genomic resources for functional genetics to help dissect factors underlying the host-pathogen interactions.

Characterizing the Pyrenophora teres f. maculata-Barley Interaction Using Pathogen Genetics.

Pyrenophora teres f. maculata is the cause of the foliar disease spot form net blotch (SFNB) on barley. To evaluate pathogen genetics underlying the P. teres f. maculata-barley interaction, we developed a 105-progeny population by crossing two globally diverse isolates, one from North Dakota and the other from Western Australia. Progeny were phenotyped on a set of four barley genotypes showing a differential reaction to the parental isolates, then genotyped using a restriction site-associated-genotype-by-sequencing (RAD-GBS) approach. Genetic maps were developed for use in quantitative trait locus (QTL) analysis to identify virulence-associated QTL. Six QTL were identified on five different linkage groups and individually accounted for 20-37% of the disease variation, with the number of significant QTL ranging from two to four for the barley genotypes evaluated. The data presented demonstrate the complexity of virulence involved in the P. teres f. maculata-barley pathosystem and begins to lay the foundation for understanding this important interaction.