Variability in Chromosome 1 of select Moroccan P. teres f. teres isolates enables isolates to overcome 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.

Beet Soil-Borne Virus Is a Helper Virus for the Novel Beta vulgaris Satellite Virus 1A.

Sugar beet (Beta vulgaris) is grown in temperate regions around the world as a source of sucrose used for natural sweetening. Sugar beet is susceptible to a number of viral diseases, but identification of the causal agent(s) under field conditions is often difficult due to mixtures of viruses that may be responsible for disease symptoms. In this study, the application of RNAseq to RNA extracted from diseased sugar beet roots obtained from the field and from greenhouse-reared plants grown in soil infested with the virus disease rhizomania (causal agent beet necrotic yellow vein virus; BNYVV) yielded genome-length sequences from BNYVV, as well as beet soil-borne virus (BSBV). The nucleotide identities of the derived consensus sequence of BSBV RNAs ranged from 99.4 to 96.7% (RNA1), 99.3 to 95.3% (RNA2), and 98.3 to 95.9% (RNA3) compared with published BSBV sequences. Based on the BSBV genome consensus sequence, clones of the genomic RNAs 1, 2, and 3 were obtained to produce RNA copies of the genome through in vitro transcription. Capped RNA produced from the clones was infectious when inoculated into leaves of Chenopodium quinoa and B. vulgaris, and extracts from transcript-infected C. quinoa leaves could infect sugar beet seedling roots through a vortex inoculation method. Subsequent exposure of these infected sugar beet seedling roots to aviruliferous Polymyxa betae, the protist vector of both BNYVV and BSBV, confirmed that BSBV derived from the infectious clones could be transmitted by the vector. Co-inoculation of BSBV synthetic transcripts with transcripts of a cloned putative satellite virus designated Beta vulgaris satellite virus 1A (BvSat1A) resulted in the production of lesions on leaves of C. quinoa similar to those produced by inoculation with BSBV alone. Nevertheless, accumulation of genomic RNA and the encoded protein of the satellite virus in co-inoculated leaves was readily detected on Northern and Western blots, respectively, whereas no accumulation of satellite virus products occurred when satellite virus RNA was inoculated alone. The predicted sequence of the detected protein encoded by BvSat1A bears hallmarks of coat proteins of other satellite viruses, and virions of a size consistent with a satellite virus were observed in samples testing positive for the virus. The results demonstrate that BSBV is a helper virus for the novel satellite virus BvSat1A.

Emergence and spread of the barley net blotch pathogen coincided with crop domestication and cultivation history.

Fungal pathogens cause devastating disease in crops. Understanding the evolutionary origin of pathogens is essential to the prediction of future disease emergence and the potential of pathogens to disperse. The fungus Pyrenophora teres f. teres causes net form net blotch (NFNB), an economically significant disease of barley. In this study, we have used 104 P. teres f. teres genomes from four continents to explore the population structure and demographic history of the fungal pathogen. We showed that P. teres f. teres is structured into populations that tend to be geographically restricted to different regions. Using Multiple Sequentially Markovian Coalescent and machine learning approaches we demonstrated that the demographic history of the pathogen correlates with the history of barley, highlighting the importance of human migration and trade in spreading the pathogen. Exploring signatures of natural selection, we identified several population-specific selective sweeps that colocalized with genomic regions enriched in putative virulence genes, and loci previously identified as determinants of virulence specificities by quantitative trait locus analyses. This reflects rapid adaptation to local hosts and environmental conditions of P. teres f. teres as it spread with barley. Our research highlights how human activities can contribute to the spread of pathogens that significantly impact the productivity of field crops.

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.

First report of Pepper yellow dwarf strain of Beet curly top virus and Spinach curly top Arizona virus in red table beet in Idaho, United States.

In the fall 2021, red table beet plants (Beta vulgaris L. cv 'Eagle') exhibiting stunted growth with shorter petioles were observed at an incidence of 10 to 15 percent in a production field in Payette County, Idaho, United States. In addition to stunting, beet leaves displayed yellowing and mild curling and crumpling, and the roots exhibited hairy root symptoms (sFig.1). To identify potential causal viruses, total RNA was isolated from the leaf and root tissue using RNeasy Plant Mini Kit (Qiagen, Valencia, CA) and subjected to high-throughput sequencing (HTS). Two libraries were prepared, one for the leaf sample and another for the root sample using a ribo-minus TruSeq Stranded Total RNA Library Prep kit (Illumina, San Diego, CA). HTS was performed with 150 bp paired-end sequencing on a NovaSeq 6000 (Novogene, Sacramento, CA). Following adapter trimming and removal of host transcripts, 5.9 and 16.2 million reads were obtained from the leaf and root samples, respectively. These reads were de novo assembled using the SPAdes assembler (Bankevitch et al., 2012; Prjibelski et al., 2020). The assembled leaf sample contigs were aligned to the NCBI non-redundant database to identify contigs matching known viruses. A single contig of 2845 nts that shared 96% coverage and 95.6% sequence identity to the pepper yellow dwarf strain of beet curly top virus (BCTV-PeYD, EU921828; Varsani et al., 2014), and 98% coverage and 98.39% identity with an isolate of BCTV-PeYD (KX529650) from Mexico, was identified in the leaf sample (GenBank Accession OP477336). To validate the HTS detection of BCTV-PeYD, total DNA was isolated from the leaf sample and a 454 bp fragment of the C1 gene (replication-associate protein) was PCR amplified and Sanger sequencing of the amplicon revealed 99.7% identity to the HTS assembled BCTV-PeYD sequence. In addition to the PeYD strain of BCTV, the Worland strain of BCTV (BCTV-Wor) was detected as a single 2930 nt contig with 100% coverage and 97.3% identity to the BCTV-Wor isolate CTS14-015 (KX867045) known to infect sugar beet in Idaho. Of note, there are 11 strains of BCTV and among those, the BCTV-Wor strain induces mild symptoms in sugar beet (Strausbaugh et al., 2017), whereas BCTV-PeYD was found only in pepper from New Mexico. Further, two contigs of 2201 nts and 523 nts were assembled generating a nearly complete genome of spinach curly top Arizona virus (SpCTAV) in the leaf sample with 99% coverage and 99.3% identity (GenBank Accession OQ703946) to the reference genome of SpCTAV (HQ443515; Hernandez-Zepeda et al., 2013). To validate the HTS results, total DNA was isolated from the leaf tissue and PCR amplified a 442 bp fragment that overlaps the V1, V2, and V3 ORFs and its sequence revealed 100% identity with the HTS assembled SpCTAV. The roots sample also showed HTS reads corresponding to BCTV-PeYD and SpCTAV. In addition, beet necrotic yellow vein virus (BNYVV) was detected in the root sample with 30% coverage, but no sequence reads matching to BNYVV was detected in the leaf sample. BNYVV is known to infect sugar beet causing rhizomania (Tamada et al., 1973; Schirmer et al., 2005). To further confirm the BNYVV HTS results, total RNA was extracted separately from the root and leaf tissue, and RT-PCR was performed with primers that were designed to amplify portions of BNYVV RNAs (Weiland et al., 2020). RT-PCR analysis generated the appropriate amplicons with expected sequences corresponding to the RNA-1, RNA-2, RNA-3, and RNA-4 of BNYVV as determined by Sanger sequencing implying BNYVV the causal agent of hairy root symptoms. Similar to observations seen for BNYVV infection in conventional sugar beet varieties, no amplification was detected for BNYVV in the RNA extracted from leaf tissue, indicating that the RT-PCR results are consistent with the HTS analysis. This is the first report of BCTV-PeYD and SpCTAV observed naturally infecting red table beet in Idaho suggesting the geographical expansion of these viruses. The co-existence of BCTV-PeYD and SpCTAV with limited host range needs to be investigated to determine the actual cause of the observed foliar symptoms. This report provides the basis for further research to understand the pathogenic nature of these viruses and their potential threat to red table beet and sugar beet production in Idaho.

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.

Resistance to QoI and DMI Fungicides Does Not Reduce Virulence of C. beticola Isolates in North Central United States.

Cercospora leaf spot (CLS) is a destructive disease limiting sugar beet production and is managed using resistant cultivars, crop rotation, and timely applications of effective fungicides. Since 2016, its causal agent, Cercospora beticola, has been reported to be resistant to quinone outside inhibitors (QoIs) and to have reduced sensitive to demethylation inhibitors (DMIs) in sugar beet growing areas in North Dakota and Minnesota. Isolates of C. beticola resistant to QoIs, DMIs, and both QoIs and DMIs were collected from fields in Foxhome, Minnesota, in 2017. Fitness of these resistant isolates was compared with that of QoI- and DMI-sensitive isolates in laboratory and greenhouse studies. In the lab, mycelial growth, spore production, and spore germination were measured. The results showed that resistant isolates had significantly less mycelial growth and spore production than sensitive isolates, while no significant difference in spore germination was detected. In the greenhouse, six leaf-stage sugar beets were inoculated with a spore suspension made from each resistant group and incubated in separate humidity chambers. CLS disease severity was evaluated visually at 7, 14, and 21 days after inoculation (DAI), and the areas under disease progress curve (AUDPC) were calculated. Resistant isolates had significantly smaller AUDPC but still caused as high disease severity as the sensitive ones at 21 DAI. Although QoI- and/or DMI-resistant isolates had a relatively slower disease development, they still caused high disease severity and need to be factored in disease management practices.

First report of Tomato bushy stunt virus naturally infecting sugar beet in the United States.

Sugar beet (Beta vulgaris L.) is an important crop grown for its sucrose content used in sugar production around the world. Tomato bushy stunt virus (TBSV) is an RNA virus that belongs to the Tombusvirus genus of the family Tombusviridae (Hearne et al., 1990). The virus was first isolated from tomato, and it is known to infect a wide range of plants (Smith, 1935; Martelli et al., 1988; Hafez et al., 2010). In 1980, a natural infection of TBSV was reported in sugar beet leaves with chlorotic and necrotic ring spots and line pattern symptoms based on serological affinity to TBSV anti-sera in Czechoslovakia (Novak and Lanzova, 1980). In March 2021, sugarbeet plants showing stunted and bushy growth with yellowing and necrotic leaves were observed in a production field in the Imperial Valley of California. Harvested roots exhibited stunted and abnormal growth compared to roots from healthy plants (sFig. 1A). These symptoms prompted a screen for potential infection by TBSV. Root-tissue harvested from the symptomatic sugar beet was initially screened using a TBSV double antibody sandwich enzyme-linked immunosorbent assay (DAS-ELISA; Agdia, Inc., Elkhart, IN), which reacted positive for TBSV. To obtain the full-length sequence of TBSV and potentially other viruses in the sample, total RNA isolated using the RNeasy Plant Mini Kit (Qiagen, Valencia, CA) from the root-tissue was subjected to high-throughput sequencing (HTS). Libraries were prepared using the TruSeq Stranded Total RNA Library Prep kit (Illumina, San Diego, CA) and sequenced using Illumina NovoSeq 6000 paired-end platform (Novogene, Sacramento, CA). A total of 52 million reads were obtained after removing the adapters and reads mapping to the host genome. These high-quality reads were de novo assembled into 75,891 contigs that are larger than 500 base pairs using the SPAdes assembler (Bankevitch et al., 2012; Prjibelski et al., 2020). The resulting contigs were searched for matching sequences to known viruses using the NCBI non-redundant database. A single contig of 4770 nts representing the full-length genome of TBSV was generated (Accession number OP477335), which showed 100% coverage to previously reported TBSV isolates 'statice' (AJ249740.1) and 'nipplefruit' (AY579432.1) with 92.19% and 91.25% nucleotide sequence identities, respectively, and thus confirming the presence of TBSV in sugar beet root-tissue. However, it showed 74% coverage with only 87% nucleotide identity to a previously reported Lettuce necrotic stunt virus (LNSV) from sugar beet, a tombusvirus that was re-classified as Moroccan pepper virus (MPV) due to high degree (>97%) of sequence identity (Obermeier et al., 2001; Wintermantel and Anchieta, 2012; Wintermantel and Hladky, 2013). The coat protein is conserved within species in tombusvirus, and it plays a significant role by providing serological relationships to tombusvirus taxonomy. The coat protein of TBSV-isolate of this study shared 98.45% and 96.91% identities at amino acid level with TBSV 'nipplefruit' (AY579432.1) and TBSV 'statice' (AJ249740.1) isolates, respectively. In contrast, it showed only 61.56% identity with the coat protein of MPV as shown in the phylogenetic tree indicating that the TBSV-isolate reported here is different from MPV (sFig. 2). To confirm the presence of TBSV, reverse-transcription (RT)-PCR was performed using the total RNA isolated from the root-tissue with primers (VR306: 5'-CGCTCACGAGCCCAGCATCCTTGA-3' and VR297: 5'-ACACCGCCACAGGAGCCATGATTG-3') designed based on the HTS data to amplify a portion of the TBSV genome. Sequencing of the RT-PCR product confirmed the presence of TBSV sequence with 99.1% identity to the TBSV-isolate identified in this study. Further, mechanical inoculation of total RNA isolated from the symptomatic sugar beet roots produced local lesions and systemic necrosis symptoms on the leaves of Chenopodium quinoa (sFig. 1B). Sequencing of the amplicon obtained using RT-PCR with primers VR306 and VR297 confirmed the presence of TBSV in C. quinoa. In addition to TBSV, several viral contigs representing Beet necrotic yellow vein virus were identified in the root-tissue indicating mixed infection in the field. To our knowledge, this is the first report that documents the occurrence of TBSV in sugar beet in the United States. Since TBSV is a soil-borne virus, our findings indicate the need for further studies focused on the frequency and coexistence of the TBSV with BNYVV in sugar beet production fields to understand the disease complexity resulting from potential mixed infections.

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.

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.

A genome-wide genetic linkage map and reference quality genome sequence for a new race in the wheat pathogen Pyrenophora tritici-repentis.

Pyrenophora tritici-repentis is an ascomycete fungus that causes tan spot of wheat. The disease has a worldwide distribution and can cause significant yield and quality losses in wheat production. The fungal pathogen is homothallic in nature, which means it can undergo sexual reproduction by selfing to produce pseudothecia on wheat stubble for seasonal survival. Since homothallism precludes the development of bi-parental fungal populations, no genetic linkage map has been developed for P. tritici-repentis for mapping and map-based cloning of fungal virulence genes. In this work, we created two heterothallic strains by deleting one of the mating type genes in each of two parental isolates 86-124 (race 2) and AR CrossB10 (a new race) and developed a bi-parental fungal population between them. The draft genome sequences of the two parental isolates were aligned to the Pt-1C-BFP reference sequence to mine single nucleotide polymorphisms (SNPs). A total of 225 SNP markers were developed for genotyping the entire population. Additionally, 75 simple sequence repeat, and two gene markers were also developed and used in the genotyping. The resulting linkage map consisted of 13 linkage groups spanning 5,075.83 cM in genetic distance. Because the parental isolate AR CrossB10 is a new race and produces Ptr ToxC, it was sequenced using long-read sequencing platforms and de novo assembled into contigs. The majority of the contigs were further anchored into chromosomes with the aid of the linkage maps. The whole genome comparison of AR CrossB10 to the reference genome of M4 revealed a few chromosomal rearrangements. The genetic linkage map and the new AR CrossB10 genome sequence are valuable tools for gene cloning in P. tritici-repentis.

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.

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.

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.

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.

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.

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.

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.

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.