Genome characterization of Pyrenophora tritici-repentis isolates reveals high plasticity and independent chromosomal location of ToxA and ToxB

Identification of a novel genetic locus conferring virulence in the wheat tan spot pathogen Pyrenophora tritici-repentis

The ascomycete Pyrenophora tritici-repentis (Ptr) is the causal agent of tan spot, a common and economically important wheat disease worldwide. Three necrotrophic effectors (NEs), known as Ptr ToxA, Ptr ToxB, and Ptr ToxC, have been identified from the fungal pathogen as major virulence factors. The race 2 isolate 86-124 which produces Ptr ToxA is capable of causing disease on wheat lines that is insensitive to Ptr ToxA, suggesting the use of additional NEs during the infection. To identify new NE gene(s) from 86-124, we developed a biparental fungal population from a cross between this isolate and the race 5 isolate DW5 using genetically modified heterothallic strains. The fungal population was genotyped with SNP and SSR markers as well as the ToxA gene, the mating type genes, and six ToxB loci. Each progeny was phenotyped onto the hard red spring wheat line CDC-Osler, which is insensitive to both Ptr ToxA and Ptr ToxB, but is highly susceptible to 86-124. The constructed genetic map consisted of 11 linkage groups that corresponded to the 11 chromosomes (chr) of the Ptr reference genome. ToxA and mating type genes mapped to the expected positions. Five of the six ToxB copies were tightly linked, residing at the distal end of chr 11, while the sixth copy was localized to the distal end of chr 5. Composite interval mapping revealed a major QTL on the distal end of chr 2 conferring virulence toward CDC-Osler by 86-124. This locus was designated as VirOsler1. Genomic sequence alignment at the locus showed a region of approximately 900 kb at the end of chr 2 absent in DW5. The identification of VirOsler1 locus provides clear evidence that the wheat tan spot pathogen uses additional virulence factors that interact with unidentified host factors for disease susceptibility.

What are the 100 most cited fungal genera?

The global diversity of fungi has been estimated between 2 to 11 million species, of which only about 155 000 have been named. Most fungi are invisible to the unaided eye, but they represent a major component of biodiversity on our planet, and play essential ecological roles, supporting life as we know it. Although approximately 20 000 fungal genera are presently recognised, the ecology of most remains undetermined. Despite all this diversity, the mycological community actively researches some fungal genera more commonly than others. This poses an interesting question: why have some fungal genera impacted mycology and related fields more than others? To address this issue, we conducted a bibliometric analysis to identify the top 100 most cited fungal genera. A thorough database search of the Web of Science, Google Scholar, and PubMed was performed to establish which genera are most cited. The most cited 10 genera are Saccharomyces, Candida, Aspergillus, Fusarium, Penicillium, Trichoderma, Botrytis, Pichia, Cryptococcus and Alternaria. Case studies are presented for the 100 most cited genera with general background, notes on their ecology and economic significance and important research advances. This paper provides a historic overview of scientific research of these genera and the prospect for further research. Citation: Bhunjun CS, Chen YJ, Phukhamsakda C, Boekhout T, Groenewald JZ, McKenzie EHC, Francisco EC, Frisvad JC, Groenewald M, Hurdeal VG, Luangsa-ard J, Perrone G, Visagie CM, Bai FY, Błaszkowski J, Braun U, de Souza FA, de Queiroz MB, Dutta AK, Gonkhom D, Goto BT, Guarnaccia V, Hagen F, Houbraken J, Lachance MA, Li JJ, Luo KY, Magurno F, Mongkolsamrit S, Robert V, Roy N, Tibpromma S, Wanasinghe DN, Wang DQ, Wei DP, Zhao CL, Aiphuk W, Ajayi-Oyetunde O, Arantes TD, Araujo JC, Begerow D, Bakhshi M, Barbosa RN, Behrens FH, Bensch K, Bezerra JDP, Bilański P, Bradley CA, Bubner B, Burgess TI, Buyck B, Čadež N, Cai L, Calaça FJS, Campbell LJ, Chaverri P, Chen YY, Chethana KWT, Coetzee B, Costa MM, Chen Q, Custódio FA, Dai YC, Damm U, de Azevedo Santiago ALCM, De Miccolis Angelini RM, Dijksterhuis J, Dissanayake AJ, Doilom M, Dong W, Alvarez-Duarte E, Fischer M, Gajanayake AJ, Gené J, Gomdola D, Gomes AAM, Hausner G, He MQ, Hou L, Iturrieta-González I, Jami F, Jankowiak R, Jayawardena RS, Kandemir H, Kiss L, Kobmoo N, Kowalski T, Landi L, Lin CG, Liu JK, Liu XB, Loizides M, Luangharn T, Maharachchikumbura SSN, Makhathini Mkhwanazi GJ, Manawasinghe IS, Marin-Felix Y, McTaggart AR, Moreau PA, Morozova OV, Mostert L, Osiewacz HD, Pem D, Phookamsak R, Pollastro S, Pordel A, Poyntner C, Phillips AJL, Phonemany M, Promputtha I, Rathnayaka AR, Rodrigues AM, Romanazzi G, Rothmann L, Salgado-Salazar C, Sandoval-Denis M, Saupe SJ, Scholler M, Scott P, Shivas RG, Silar P, Souza-Motta CM, Silva-Filho AGS, Spies CFJ, Stchigel AM, Sterflinger K, Summerbell RC, Svetasheva TY, Takamatsu S, Theelen B, Theodoro RC, Thines M, Thongklang N, Torres R, Turchetti B, van den Brule T, Wang XW, Wartchow F, Welti S, Wijesinghe SN, Wu F, Xu R, Yang ZL, Yilmaz N, Yurkov A, Zhao L, Zhao RL, Zhou N, Hyde KD, Crous PW (2024). What are the 100 most cited fungal genera? Studies in Mycology 108: 1-411. doi: 10.3114/sim.2024.108.01.

Evolution of the Toxb Gene in Pyrenophora tritici-repentis and Related Species

Pyrenophora tritici-repentis (tan spot) is a destructive foliar pathogen of wheat with global impact. This ascomycete fungus possesses a highly plastic open pangenome shaped by the gain and loss of effector genes. This study investigated the allelic variations in the chlorosis-encoding gene ToxB across 422 isolates representing all identified pathotypes and worldwide origins. To gain better insights into ToxB evolution, we examined its presence and variability in other Pyrenophora spp. A ToxB haplotype network was constructed, revealing the evolutionary relationships of this gene (20 haplotypes) across four Pyrenophora species. Notably, toxb, the homolog of ToxB, was detected for the first time in the barley pathogen Pyrenophora teres. The ToxB/toxb genes display evidence of selection that is characterized by loss of function, duplication, and diverse mutations. Within the ToxB/toxb open reading frame, 72 mutations were identified, including 14 synonymous, 55 nonsynonymous, and 3 indel mutations. Remarkably, a, ∼5.6-kb Copia-like retrotransposon, named Copia-1_Ptr, was found inserted in the toxb gene of a race 3 isolate. This insert disrupted the ToxB gene's function, a first case of effector gene disruption by a transposable element in P. tritici-repentis. Additionally, a microsatellite with 25 nucleotide repeats (0 to 10) in the upstream region of ToxB suggested a potential mechanism influencing ToxB expression and regulation. Exploring ToxB-like protein distribution in other ascomycetes revealed the presence of ToxB-like proteins in 19 additional species, including the Leotiomycetes class for the first time. The presence/absence pattern of ToxB-like proteins defied species relatedness compared with a phylogenetic tree, suggesting a past horizontal gene transfer event during the evolution of the ToxB gene. [Formula: see text] Copyright © 2024 His Majesty the King in Right of Canada, as represented by the Minister of Agriculture and Agri-Food. This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Starships are active eukaryotic transposable elements mobilized by a new family of tyrosine recombinases

Transposable elements in eukaryotic organisms have historically been considered "selfish," at best conferring indirect benefits to their host organisms. The Starships are a recently discovered feature in fungal genomes that are, in some cases, predicted to confer beneficial traits to their hosts and also have hallmarks of being transposable elements. Here, we provide experimental evidence that Starships are indeed autonomous transposons, using the model Paecilomyces variotii, and identify the HhpA "Captain" tyrosine recombinase as essential for their mobilization into genomic sites with a specific target site consensus sequence. Furthermore, we identify multiple recent horizontal gene transfers of Starships, implying that they jump between species. Fungal genomes have mechanisms to defend against mobile elements, which are frequently detrimental to the host. We discover that Starships are also vulnerable to repeat-induced point mutation defense, thereby having implications on the evolutionary stability of such elements.

Important wheat diseases in the US and their management in the 21st century

Wheat is a crop of historical significance, as it marks the turning point of human civilization 10,000 years ago with its domestication. Due to the rapid increase in population, wheat production needs to be increased by 50% by 2050 and this growth will be mainly based on yield increases, as there is strong competition for scarce productive arable land from other sectors. This increasing demand can be further achieved using sustainable approaches including integrated disease pest management, adaption to warmer climates, less use of water resources and increased frequency of abiotic stress tolerances. Out of 200 diseases of wheat, 50 cause economic losses and are widely distributed. Each year, about 20% of wheat is lost due to diseases. Some major wheat diseases are rusts, smut, tan spot, spot blotch, fusarium head blight, common root rot, septoria blotch, powdery mildew, blast, and several viral, nematode, and bacterial diseases. These diseases badly impact the yield and cause mortality of the plants. This review focuses on important diseases of the wheat present in the United States, with comprehensive information of causal organism, economic damage, symptoms and host range, favorable conditions, and disease management strategies. Furthermore, major genetic and breeding efforts to control and manage these diseases are discussed. A detailed description of all the QTLs, genes reported and cloned for these diseases are provided in this review. This study will be of utmost importance to wheat breeding programs throughout the world to breed for resistance under changing environmental conditions.

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

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.

A Conserved Hypothetical Gene Is Required but Not Sufficient for Ptr ToxC Production in Pyrenophora tritici-repentis

The fungus Pyrenophora tritici-repentis causes tan spot, an important foliar disease of wheat worldwide. The fungal pathogen produces three necrotrophic effectors, namely Ptr ToxA, Ptr ToxB, and Ptr ToxC to induce necrosis or chlorosis in wheat. Both Ptr ToxA and Ptr ToxB are proteins, and their encoding genes have been cloned. Ptr ToxC was characterized as a low-molecular weight molecule 20 years ago but the one or more genes controlling its production in P. tritici-repentis are unknown. Here, we report the genetic mapping, molecular cloning, and functional analysis of a fungal gene that is required for Ptr ToxC production. The genetic locus controlling the production of Ptr ToxC, termed ToxC, was mapped to a subtelomeric region using segregating biparental populations, genome sequencing, and association analysis. Additional marker analysis further delimited ToxC to a 173-kb region. The predicted genes in the region were examined for presence/absence polymorphism in different races and isolates leading to the identification of a single candidate gene. Functional validation showed that this gene was required but not sufficient for Ptr ToxC production, thus it is designated as ToxC1. ToxC1 encoded a conserved hypothetical protein likely located on the vacuole membrane. The gene was highly expressed during infection, and only one haplotype was identified among 120 isolates sequenced. Our work suggests that Ptr ToxC is not a protein and is likely produced through a cascade of biosynthetic pathway. The identification of ToxC1 is a major step toward revealing the Ptr ToxC biosynthetic pathway and studying its molecular interactions with host factors.[Formula: see text] Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Fungal Pathogens in Grasslands

Grasslands are major primary producers and function as major components of important watersheds. Although a concise definition of grasslands cannot be given using a physiognomic or structural approach, grasslands can be described as vegetation communities experiencing periodical droughts and with canopies dominated by grasses and grass-like plants. Grasslands have a cosmopolitan distribution except for the Antarctic region. Fungal interactions with grasses can be pathogenic or symbiotic. Herbivorous mammals, insects, other grassland animals, and fungal pathogens are known to play important roles in maintaining the biomass and biodiversity of grasslands. Although most pathogenicity studies on the members of Poaceae have been focused on economically important crops, the plant-fungal pathogenic interactions involved can extend to the full range of ecological circumstances that exist in nature. Hence, it is important to delineate the fungal pathogen communities and their interactions in man-made monoculture systems and highly diverse natural ecosystems. A better understanding of the key fungal players can be achieved by combining modern techniques such as next-generation sequencing (NGS) together with studies involving classic phytopathology, taxonomy, and phylogeny. It is of utmost importance to develop experimental designs that account for the ecological complexity of the relationships between grasses and fungi, both above and below ground. In grasslands, loss in species diversity increases interactions such as herbivory, mutualism, predation or infectious disease transmission. Host species density and the presence of heterospecific host species, also affect the disease dynamics in grasslands. Many studies have shown that lower species diversity increases the severity as well as the transmission rate of fungal diseases. Moreover, communities that were once highly diverse but have experienced decreased species richness and dominancy have also shown higher pathogenicity load due to the relaxed competition, although this effect is lower in natural communities. This review addresses the taxonomy, phylogeny, and ecology of grassland fungal pathogens and their interactions in grassland ecosystems.

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.

PacBio genome sequencing reveals new insights into the genomic organisation of the multi-copy ToxB gene of the wheat fungal pathogen Pyrenophora tritici-repentis

Background

Necrotrophic effector proteins secreted by fungal pathogens are important virulence factors that mediate the development of disease in wheat. Pyrenophora tritici-repentis (Ptr), the causal agent of wheat tan spot, has a race structure dependent on the combination of effectors. In Ptr, ToxA and ToxB are known proteinaceous effectors responsible for necrosis and chlorosis respectively. While Ptr ToxA is encoded by the single gene ToxA, ToxB has multiple loci in the Ptr genome, which is postulated to be directly related to the level of ToxB production and leaf chlorosis. Although previous analysis has indicated that the majority of the ToxB loci lie on a single chromosome, the exact number and chromosomal locations for all the ToxB loci have not been fully identified.

Results

In this study, we have sequenced the genome of a race 5 ToxB-producing isolate (DW5), using PacBio long read technology, and found that ToxB duplications are nested in the complex subtelomeric chromosomal regions. A total of ten identical ToxB gene copies were identified and based on flanking sequence identity, nine loci appeared associated with chromosome 11 and a single copy with chromosome 5. Chromosome 11 multiple ToxB gene loci were separated by large sequence regions between 31 and 66 kb within larger segmental duplications in an alternating pattern related to loci strand, and flanked by transposable elements.

Conclusion

This work provides for the first time the full accompaniment of ToxB loci and surrounding regions, and identifies the organization and distribution of ten ToxB loci to subtelomeric regions. To our knowledge, this is the first report of an interwoven strand-related duplication pattern event. This study further highlights the importance of resolving the highly complex distal chromosomal regions, that remain difficult to assemble, and can harbour important effectors and virulence factors.

Supplementary information

Supplementary information accompanies this paper at 10.1186/s12864-020-07029-4.

Pyrenophora tritici-repentis in Tunisia: Race Structure and Effector Genes

Tan spot is a destructive foliar wheat disease worldwide and caused by the ascomycete fungus Pyrenophora tritici-repentis (Ptr); it has become more frequent in Tunisia over the last decade. In this study, the virulence of 73 single-spore isolates, collected from durum and bread wheat fields during 2017-2018 growing season, was evaluated on four differential wheat genotypes. This was followed by polymerase chain reaction tests with specific primers for the effector genes ToxA, ToxB, and toxb (ToxB-homolog). Sequence analysis to validate the identity of the amplified genes was followed, and ToxA amplicons from a subset of 22 isolates were analyzed to determine its haplotype identity. Ptr isolates from Tunisia were grouped in races 2, 4, 5, and 7, and 44% of the tested isolates did not fit under any known race, and were denoted here as atypical. These atypical isolates induced the same symptoms as race 7 isolates, extensive necrosis, and chlorosis on susceptible genotypes, but lacked the ToxA gene. ToxA is the only identified necrosis-inducing effector in Ptr, and was amplified in 51% of tested isolates, and shared identical sequence to previously identified haplotype (H15). ToxB and its homolog toxb were present in 97% and 93% of tested isolates, respectively. Ptr in Tunisia lacked Ptr ToxC activity, and none of the tested isolates induced the specific symptoms of that effector. Race 7 and the atypical isolates dominated the Tunisian Ptr population, while races 2, 4, and 5 were found at low percentages. In conclusion, ToxB and its homolog were the most dominant genes in Ptr from Tunisia, and the majority of the isolates induced necrosis and chlorosis on Ptr ToxA and Ptr ToxB susceptible wheat genotypes. However, only about half of that necrosis can be attributed to ToxA presence, this result necessitates further research to investigate the prevalence of additional necrotic effector(s). Terminology: in this paper, Pyrenophora tritici-repentis abbreviated as Ptr, the effectors are referred to by Ptr ToxA, Ptr ToxB and Ptr ToxC, and the genes coding for them are written in italic as ToxA, ToxB, and ToxC, respectively.

A new PacBio genome sequence of an Australian Pyrenophora tritici-repentis race 1 isolate

OBJECTIVES

The necrotrophic fungal pathogen Pyrenophora tritici-repentis (Ptr) is the causal agent of tan spot a major disease of wheat. We have generated a new genome resource for an Australian Ptr race 1 isolate V1 to support comparative 'omics analyses. In particular, the V1 PacBio Biosciences long-read sequence assembly was generated to confirm the stability of large-scale genome rearrangements of the Australian race 1 isolate M4 when compared to the North American race 1 isolate Pt-1C-BFP.

RESULTS

Over 1.3 million reads were sequenced by PacBio Sequel small-molecule real-time sequencing (SRMT) cell to yield 11.4 Gb for the genome assembly of V1 (285X coverage), with median and maximum read lengths of 8959 bp and 72,292 bp respectively. The V1 genome was assembled into 33 contiguous sequences with a of total length 40.4 Mb and GC content of 50.44%. A total of 14,050 protein coding genes were predicted and annotated for V1. Of these 11,519 genes were orthologous to both Pt-1C-BFP and M4. Whole genome alignment of the Australian long-read assemblies (V1 to M4) confirmed previously identified large-scale genome rearrangements between M4 and Pt-1C-BFP and presented small scale variations, which included a sequence break within a race-specific region for ToxA, a well-known necrotrophic effector gene.

Transposon-Mediated Horizontal Transfer of the Host-Specific Virulence Protein ToxA between Three Fungal Wheat Pathogens

This work dissects the tripartite horizontal transfer of ToxA, a gene that has a direct negative impact on global wheat yields. Defining the extent of horizontally transferred DNA is important because it can provide clues to the mechanisms that facilitate HGT. Our analysis of ToxA and its surrounding 14 kb suggests that this gene was horizontally transferred in two independent events, with one event likely facilitated by a type II DNA transposon. These horizontal transfer events are now in various processes of decay in each species due to the repeated insertion of new transposons and subsequent rounds of targeted mutation by a fungal genome defense mechanism known as repeat induced point mutation. This work highlights the role that HGT plays in the evolution of host adaptation in eukaryotic pathogens. It also increases the growing body of evidence indicating that transposons facilitate adaptive HGT events between fungi present in similar environments and hosts.

Most known examples of horizontal gene transfer (HGT) between eukaryotes are ancient. These events are identified primarily using phylogenetic methods on coding regions alone. Only rarely are there examples of HGT where noncoding DNA is also reported. The gene encoding the wheat virulence protein ToxA and the surrounding 14 kb is one of these rare examples. ToxA has been horizontally transferred between three fungal wheat pathogens (Parastagonospora nodorum, Pyrenophora tritici-repentis, and Bipolaris sorokiniana) as part of a conserved ∼14 kb element which contains coding and noncoding regions. Here we used long-read sequencing to define the extent of HGT between these three fungal species. Construction of near-chromosomal-level assemblies enabled identification of terminal inverted repeats on either end of the 14 kb region, typical of a type II DNA transposon. This is the first description of ToxA with complete transposon features, which we call ToxhAT. In all three species, ToxhAT resides in a large (140-to-250 kb) transposon-rich genomic island which is absent in isolates that do not carry the gene (annotated here as toxa⁻). We demonstrate that the horizontal transfer of ToxhAT between P. tritici-repentis and P. nodorum occurred as part of a large (∼80 kb) HGT which is now undergoing extensive decay. In B. sorokiniana, in contrast, ToxhAT and its resident genomic island are mobile within the genome. Together, these data provide insight into the noncoding regions that facilitate HGT between eukaryotes and into the genomic processes which mask the extent of HGT between these species.

Comparative genomics of the wheat fungal pathogen Pyrenophora tritici-repentis reveals chromosomal variations and genome plasticity

Background

Pyrenophora tritici-repentis (Ptr) is a necrotrophic fungal pathogen that causes the major wheat disease, tan spot. We set out to provide essential genomics-based resources in order to better understand the pathogenicity mechanisms of this important pathogen.

Results

Here, we present eight new Ptr isolate genomes, assembled and annotated; representing races 1, 2 and 5, and a new race. We report a high quality Ptr reference genome, sequenced by PacBio technology with Illumina paired-end data support and optical mapping. An estimated 98% of the genome coverage was mapped to 10 chromosomal groups, using a two-enzyme hybrid approach. The final reference genome was 40.9 Mb and contained a total of 13,797 annotated genes, supported by transcriptomic and proteogenomics data sets.

Conclusions

Whole genome comparative analysis revealed major chromosomal segmental rearrangements and fusions, highlighting intraspecific genome plasticity in this species. Furthermore, the Ptr race classification was not supported at the whole genome level, as phylogenetic analysis did not cluster the ToxA producing isolates. This expansion of available Ptr genomics resources will directly facilitate research aimed at controlling tan spot disease.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4680-3) contains supplementary material, which is available to authorized users.

Fungal Plant Pathogenesis Mediated by Effectors

The interactions between fungi and plants encompass a spectrum of ecologies ranging from saprotrophy (growth on dead plant material) through pathogenesis (growth of the fungus accompanied by disease on the plant) to symbiosis (growth of the fungus with growth enhancement of the plant). We consider pathogenesis in this article and the key roles played by a range of pathogen-encoded molecules that have collectively become known as effectors.

Phylogenomic evolutionary surveys of subtilase superfamily genes in fungi

Subtilases belong to a superfamily of serine proteases which are ubiquitous in fungi and are suspected to have developed distinct functional properties to help fungi adapt to different ecological niches. In this study, we conducted a large-scale phylogenomic survey of subtilase protease genes in 83 whole genome sequenced fungal species in order to identify the evolutionary patterns and subsequent functional divergences of different subtilase families among the main lineages of the fungal kingdom. Our comparative genomic analyses of the subtilase superfamily indicated that extensive gene duplications, losses and functional diversifications have occurred in fungi, and that the four families of subtilase enzymes in fungi, including proteinase K-like, Pyrolisin, kexin and S53, have distinct evolutionary histories which may have facilitated the adaptation of fungi to a broad array of life strategies. Our study provides new insights into the evolution of the subtilase superfamily in fungi and expands our understanding of the evolution of fungi with different lifestyles.

De novo genome assembly of the fungal plant pathogen Pyrenophora semeniperda

Pyrenophora semeniperda (anamorph Drechslera campulata) is a necrotrophic fungal seed pathogen that has a wide host range within the Poaceae. One of its hosts is cheatgrass (Bromus tectorum), a species exotic to the United States that has invaded natural ecosystems of the Intermountain West. As a natural pathogen of cheatgrass, P. semeniperda has potential as a biocontrol agent due to its effectiveness at killing seeds within the seed bank; however, few genetic resources exist for the fungus. Here, the genome of P. semeniperda isolate assembled from sequence reads of 454 pyrosequencing is presented. The total assembly is 32.5 Mb and includes 11,453 gene models encoding putative proteins larger than 24 amino acids. The models represent a variety of putative genes that are involved in pathogenic pathways typically found in necrotrophic fungi. In addition, extensive rearrangements, including inter- and intrachromosomal rearrangements, were found when the P. semeniperda genome was compared to P. tritici-repentis, a related fungal species.

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.

Karyotype polymorphism and chromosomal rearrangement in populations of the phytopathogenic fungus, Ascochyta rabiei

The fungus Ascochyta rabiei is the causal agent of Ascochyta blight of chickpea and the most serious threat to chickpea production. Little is currently known about the genome size or organization of A. rabiei. Given recent genome sequencing efforts, characterization of the genome at a population scale will provide a framework for genome interpretation and direction of future resequencing efforts. Electrophoretic karyotype profiles of 112 isolates from 21 countries revealed 12-16 chromosomes between 0.9 Mb and 4.6 Mb with an estimated genome size of 23 Mb-34 Mb. Three general karyotype profiles A, B, and C were defined by the arrangement of the largest chromosomes. Approximately one-third of isolates (group A) possessed a chromosome larger than 4.0 Mb that was absent from group B and C isolates. The ribosomal RNA gene (rDNA) cluster was assigned to the largest chromosome in all except four isolates (group C) whose rDNA cluster was located on the second largest chromosome (3.2 Mb). Analysis of progeny from an in vitro sexual cross between two group B isolates revealed one of 16 progeny with an rDNA-encoding chromosome larger than 4.0 Mb similar to group A isolates, even though a chromosome of this size was not present in either parent. No expansion of the rDNA cluster was detected in the progeny, indicating the increase in chromosome size was not due to an expansion in number of rDNA repeats. The karyotype of A. rabiei is relatively conserved when compared with published examples of asexual ascomycetes, but labile with the potential for large scale chromosomal rearrangements during meiosis. The results of this study will allow for the targeted sequencing of specific isolates to determine the molecular mechanisms of karyotype variation within this species.

RNA-mediated gene silencing of ToxB in Pyrenophora tritici-repentis

The fungus Pyrenophora tritici-repentis causes tan spot, a wheat leaf disease of worldwide importance. The pathogen produces three host-selective toxins, including Ptr ToxB, which causes chlorophyll degradation and foliar chlorosis on toxin-sensitive wheat genotypes. The ToxB gene, which codes for Ptr ToxB, was silenced in a wild-type race 5 isolate of the fungus thorough a sense- and antisense-mediated silencing mechanism. Toxin production by the silenced strains was evaluated in culture filtrates of the fungus via Western blotting analysis, and plant bioassays were conducted to test the virulence of the transformants in planta. The chlorosis-inducing ability of the silenced strains was correlated with the quantity of Ptr ToxB, and transformants in which toxin production was strongly decreased also caused very little disease on toxin-sensitive wheat genotypes. Cytological analysis of the infection process revealed that, in addition to a reduced capacity to induce chlorosis, the silenced strains with the greatest decrease in the levels of Ptr ToxB produced significantly fewer appressoria than the wild-type isolate, 12 and 24 h after inoculation onto wheat leaves. The results provide strong support for the suggestion that the amount of Ptr ToxB protein produced by fungal isolates plays a significant role in the quantitative variation in the virulence of P. tritici-repentis.

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.

Pyrenophora bromi, causal agent of brownspot of bromegrass, expresses a gene encoding a protein with homology and similar activity to Ptr ToxB, a host-selective toxin of wheat

Ptr ToxB, encoded by ToxB, is one of multiple host-selective toxins (HST) produced by the wheat pathogen Pyrenophora tritici-repentis. Homologs of ToxB are found in several ascomycetes, including sister species Pyrenophora bromi, causal agent of brownspot of bromegrass. Due to the close evolutionary relatedness of P. tritici-repentis and P. bromi and that of their grass hosts, we hypothesized that homologs of ToxB in P. bromi may act as HST in the disease interaction between P. bromi and bromegrass. A representative set of transcriptionally active P. bromi ToxB genes were heterologously expressed in Pichia pastoris and the resultant proteins tested for their ability to act as HST on bromegrass. The tested Pyrenophora bromi ToxB (Pb ToxB) proteins were not toxic to bromegrass; thus, Pb ToxB does not appear to function as an HST in the P. bromi-bromegrass interaction. Instead, we revealed that the Pb ToxB proteins can be toxic to Ptr ToxB-sensitive wheat, at levels similar to Ptr ToxB, and the corresponding P. bromi ToxB genes are expressed in P. bromi-inoculated wheat. Our data suggest that P. bromi possesses the potential to become a wheat pathogen and highlights the importance of investigating the interaction between P. bromi and wheat.