Analysis of the mating-type loci of co-occurring and phylogenetically related species of Ascochyta and Phoma

Genome and Transcriptome Analysis of Ascochyta pisi Provides Insights into the Pathogenesis of Ascochyta Blight of Pea

Ascochyta blight caused by Ascochyta pisi is a major constraint to pea (Pisum sativum L.) production worldwide. Deciphering the pathogenic mechanism of A. pisi on peas will help in breeding resistant pea varieties and developing effective approaches for disease management. However, little is known about the genomic features and pathogenic factors of A. pisi. In this study, we first report that A. pisi is one of the causal agents of ascochyta blight disease of pea in China. The genome of the representative isolate A. pisi HNA23 was sequenced using PacBio and Illumina sequencing technologies. The HNA23 genome assembly is almost 41.5 Mb in size and harbors 10,796 putative protein-encoding genes. We predicted 555 carbohydrate-active enzymes (CAZymes), 1,008 secreted proteins, 74 small secreted cysteine-rich proteins (SSCPs), and 26 secondary metabolite biosynthetic gene clusters (SMGCs). A comparison of A. pisi genome features with the features of 6 other available genomes of Ascochyta species showed that CAZymes, the secretome, and SMGCs of this genus are considerably conserved. Importantly, the transcriptomes of HNA23 during infection of peas at three stages were further analyzed. We found that 245 CAZymes and 29 SSCPs were upregulated at all three tested infection stages. SMGCs were also trigged, but most of them were induced at only one stage of infection. Together, our results provide important genomic information on Ascochyta spp. and offer insights into the pathogenesis of A. pisi. IMPORTANCE Ascochyta blight is a major disease of legumes worldwide. Ascochyta pisi and other Ascochyta species have been identified as pathogens of ascochyta blight. Here, we first report that A. pisi causes ascochyta blight of pea in China, and we report the high-quality, fully annotated genome of A. pisi. Comparative genome analysis was performed to elucidate the differences and similarities among 7 Ascochyta species. We predict abundant CAZymes (569 per species), secreted proteins (851 per species), and prolific secondary metabolite gene clusters (29 per species) in these species. We identified a set of genes that may be responsible for fungal virulence based on transcriptomes in planta, including CAZymes, SSCPs, and secondary metabolites. The findings from the comparative genome analysis highlight the genetic diversity and help in understanding the evolutionary relationship of Ascochyta species. In planta transcriptome analysis provides reliable information for further investigation of the mechanism of the interaction between Ascochyta spp. and legumes.

De Novo Long-Read Whole-Genome Assemblies and the Comparative Pan-Genome Analysis of Ascochyta Blight Pathogens Affecting Field Pea

Ascochyta Blight (AB) is a major disease of many cool-season legumes globally. In field pea, three fungal pathogens have been identified to be responsible for this disease in Australia, namely Peyronellaea pinodes, Peyronellaea pinodella and Phoma koolunga. Limited genomic resources for these pathogens have been generated, which has hampered the implementation of effective management strategies and breeding for resistant cultivars. Using Oxford Nanopore long-read sequencing, we report the first high-quality, fully annotated, near-chromosome-level nuclear and mitochondrial genome assemblies for 18 isolates from the Australian AB complex. Comparative genome analysis was performed to elucidate the differences and similarities between species and isolates using phylogenetic relationships and functional diversity. Our data indicated that P. pinodella and P. koolunga are heterothallic, while P. pinodes is homothallic. More homology and orthologous gene clusters are shared between P. pinodes and P. pinodella compared to P. koolunga. The analysis of the repetitive DNA content showed differences in the transposable repeat composition in the genomes and their expression in the transcriptomes. Significant repeat expansion in P. koolunga's genome was seen, with strong repeat-induced point mutation (RIP) activity being evident. Phylogenetic analysis revealed that genetic diversity can be exploited for species marker development. This study provided the much-needed genetic resources and characterization of the AB species to further drive research in key areas such as disease epidemiology and host-pathogen interactions.

The role of internal transcribed spacer 2 secondary structures in classifying mycoparasitic Ampelomyces

Many fungi require specific growth conditions before they can be identified. Direct environmental DNA sequencing is advantageous, although for some taxa, specific primers need to be used for successful amplification of molecular markers. The internal transcribed spacer region is the preferred DNA barcode for fungi. However, inter- and intra-specific distances in ITS sequences highly vary among some fungal groups; consequently, it is not a solely reliable tool for species delineation. Ampelomyces, mycoparasites of the fungal phytopathogen order Erysiphales, can have ITS genetic differences up to 15%; this may lead to misidentification with other closely related unknown fungi. Indeed, Ampelomyces were initially misidentified as other pycnidial mycoparasites, but subsequent research showed that they differ in pycnidia morphology and culture characteristics. We investigated whether the ITS2 nucleotide content and secondary structure was different between Ampelomyces ITS2 sequences and those unrelated to this genus. To this end, we retrieved all ITS sequences referred to as Ampelomyces from the GenBank database. This analysis revealed that fungal ITS environmental DNA sequences are still being deposited in the database under the name Ampelomyces, but they do not belong to this genus. We also detected variations in the conserved hybridization model of the ITS2 proximal 5.8S and 28S stem from two Ampelomyces strains. Moreover, we suggested for the first time that pseudogenes form in the ITS region of this mycoparasite. A phylogenetic analysis based on ITS2 sequences-structures grouped the environmental sequences of putative Ampelomyces into a different clade from the Ampelomyces-containing clades. Indeed, when conducting ITS2 analysis, resolution of genetic distances between Ampelomyces and those putative Ampelomyces improved. Each clade represented a distinct consensus ITS2 S2, which suggested that different pre-ribosomal RNA (pre-rRNA) processes occur across different lineages. This study recommends the use of ITS2 S2s as an important tool to analyse environmental sequencing and unveiling the underlying evolutionary processes.

Cryptic diversity found in Didymellaceae from Australian native legumes

Ascochytakoolunga (Didymellaceae, Pleosporales) was first described in 2009 (as Phomakoolunga) and identified as the causal agent of Ascochyta blight of Pisumsativum (field pea) in South Australia. Since then A.koolunga has not been reported anywhere else in the world, and its origins and occurrence on other legume (Fabaceae) species remains unknown. Blight and leaf spot diseases of Australian native, pasture and naturalised legumes were studied to investigate a possible native origin of A.koolunga.

Ascochytakoolunga was not detected on native, naturalised or pasture legumes that had leaf spot symptoms, in any of the studied regions in southern Australia, and only one isolate was recovered from P.sativum. However, we isolated five novel species in the Didymellaceae from leaf spots of Australian native legumes from commercial field pea regions throughout southern Australia. The novel species were classified on the basis of morphology and phylogenetic analyses of the internal transcribed spacer region and part of the RNA polymerase II subunit B gene region. Three of these species, Nothophomagarlbiwalawardasp. nov., Nothophomanaiawusp. nov. and Nothophomangayawangsp. nov., were isolated from Sennaartemisioides. The other species described here are Epicoccumdjirangnandirisp. nov. from Swainsonagalegifolia and Neodidymelliopsistinkyukukusp. nov. from Hardenbergiaviolacea. In addition, we report three new host-pathogen associations in Australia, namely Didymellapinodes on S.artemisioides and Viciacracca, and D.lethalis on Lathyrustingitanus. This is also the first report of Didymellaprosopidis in Australia.

Repeated evolution of self-compatibility for reproductive assurance

Sexual reproduction in eukaryotes requires the fusion of two compatible gametes of opposite sexes or mating types. To meet the challenge of finding a mating partner with compatible gametes, evolutionary mechanisms such as hermaphroditism and self-fertilization have repeatedly evolved. Here, by combining the insights from comparative genomics, computer simulations and experimental evolution in fission yeast, we shed light on the conditions promoting separate mating types or self-compatibility by mating-type switching. Analogous to multiple independent transitions between switchers and non-switchers in natural populations mediated by structural genomic changes, novel switching genotypes readily evolved under selection in the experimental populations. Detailed fitness measurements accompanied by computer simulations show the benefits and costs of switching during sexual and asexual reproduction, governing the occurrence of both strategies in nature. Our findings illuminate the trade-off between the benefits of reproductive assurance and its fitness costs under benign conditions facilitating the evolution of self-compatibility.

Mating-type switching enables self-compatible reproduction in fungi, but switching ability is variable even within species. Here, the authors find de novo evolution of switching genotypes in experimentally evolved fission yeast populations and show a trade-off between mating success and growth.

Mixed-reproductive strategies, competitive mating-type distribution and life cycle of fourteen black morel species

Morchella species are well known world-round as popular and prized edible fungi due to their unique culinary flavor. Recently, several species have been successfully cultivated in China. However, their reproductive modes are still unknown, and their basic biology needs to be elucidated. Here, we use the morel genome information to investigate mating systems and life cycles of fourteen black morel species. Mating type-specific primers were developed to screen and genotype ascospores, hymenia and stipes from 223 ascocarps of the 14 species from Asia and Europe. Our data indicated that they are all heterothallic and their life cycles are predominantly haploid, but sterile haploid fruiting also exists. Ascospores in all species are mostly haploid, homokaryotic, and multinuclear, whereas aborted ascospores without any nuclei were also detected. Interestingly, we monitored divergent spatial distribution of both mating types in natural morel populations and cultivated sites, where the fertile tissue of fruiting bodies usually harbored both mating types, whereas sterile tissue of wild morels constantly had one MAT allele, while the sterile tissue of cultivated strains always exhibited both MAT alleles. Furthermore, MAT1-1-1 was detected significantly more commonly than MAT1-2-1 in natural populations, which strongly suggested a competitive advantage for MAT1-1 strains.

Mating-Type Gene Structure and Spatial Distribution of Didymella tanaceti in Pyrethrum Fields

Tan spot of pyrethrum (Tanacetum cinerariifolium) is caused by the ascomycete Didymella tanaceti. To assess the evolutionary role of ascospores in the assumed asexual species, the structure and arrangement of mating-type (MAT) genes were examined. A single MAT1-1 or MAT1-2 idiomorph was identified in all isolates examined, indicating that the species is heterothallic. The idiomorphs were flanked upstream and downstream by regions encoding pyridoxamine phosphate oxidase-like and DNA lyase-like proteins, respectively. A multiplex MAT-specific polymerase chain reaction assay was developed and used to genotype 325 isolates collected within two transects in each of four fields in Tasmania, Australia. The ratio of isolates of each mating-type in each transect was consistent with a 1:1 ratio. The spatial distribution of the isolates of the two mating-types within each transect was random for all except one transect for MAT1-1 isolates, indicating that clonal patterns of each mating-type were absent. However, evidence of a reduced selection pressure on MAT1-1 isolates was observed, with a second haplotype of the MAT1-1-1 gene identified in 4.4% of MAT1-1 isolates. In vitro crosses between isolates with opposite mating-types failed to produce ascospores. Although the sexual morph could not be induced, the occurrence of both mating-types in equal frequencies suggested that a cryptic sexual mode of reproduction may occur within field populations.

Development of a loop-mediated isothermal amplification method for the rapid diagnosis of Ascochyta rabiei L. in chickpeas

Ascochyta blight (AB) is a devastating fungal disease of chickpeas that has spread to nearly all of the chickpea cultivating regions of the world. The rapid diagnosis of Ascochyta rabiei L. (A. rabiei), the pathogen that causes AB, plays an important role in A. rabiei epidemic tracking and AB management. In this study, a group of loop-mediated isothermal amplification (LAMP) primers was designed to detect A. rabiei in chickpea plants and seeds via a LAMP method and a conventional PCR method based on an internal transcribed spacer (ITS) sequence analysis of A. rabiei. Compared with the conventional PCR method, the LAMP method not only exhibited greater sensitivity and specificity in the detection of A. rabiei but also used simpler equipment and required less operational time. The minimum detectable concentration of the A. rabiei genomic DNA solution with the LAMP method was 6.01 × 10(-6 )ng/μl, which was 100 times lower than that of the conventional PCR method with the same outer primers. The greatest advantage of the LAMP method is that results can be observed via the visualization of color changes in SYBR Green I dye with the naked eye, and it does not require expensive instruments, also with less time consumption.

Characterization and distribution of mating-type genes of the turfgrass pathogen Sclerotinia homoeocarpa on a global scale

Sclerotinia homoeocarpa F.T. Bennett is a filamentous member of Ascomycota that causes dollar spot, the most economically important disease of turfgrass worldwide. We sequenced and characterized the mating-type (MAT) locus of four recently-collected contemporary strains causing dollar spot, four historical type strains used to describe the fungus, and three species of Rutstroemiaceae. Moreover, we developed a multiplex PCR assay to screen 1019 contemporary isolates for mating-type. The organization of the MAT loci of all strains examined could be classified into one of four categories: (1) putatively heterothallic, as exemplified by all contemporary strains and three of four historical type strains; (2) putatively heterothallic with a deleted putative gene in the MAT1-2 idiomorph, as detected in strains from two recently-collected populations in the United Kingdom that show more similarity to historical strains; (3) putatively homothallic with close physical linkage between MAT1-1-1 and MAT1-2-1, as found in one historical type strain of S. homoeocarpa and two strains of Rutstroemia cuniculi; and (4) an unresolved but apparently homothallic organization in which strains contained both MAT1-1-1 and MAT1-2-1 but linkage between these genes and between the two flanking genes could not be confirmed, as identified in R. paludosa and Poculum henningsianum. In contemporary S. homoeocarpa populations there was no significant difference in the frequency of the two mating types in clone-corrected samples when analyzed on regional and local scales, suggesting sex may be possible in this pathogen. However, two isolates from Italy and twenty from California were heterokaryotic for both complete heterothallic MAT idiomorphs. Results from this study contribute to knowledge about mating systems in filamentous fungi and enhance our understanding of the evolution and biology of an important plant pathogen.

Characterization of mating type genes supports the hypothesis that Stagonosporopsis chrysanthemi is homothallic and provides evidence that Stagonosporopsis tanaceti is heterothallic

To understand the organization of the mating type locus of Stagonosporopsis tanaceti and Stagonosporopsis chrysanthemi, and its potential role in the epidemiology of ray blight of pyrethrum and chrysanthemum, respectively, the mating type (MAT) locus of these species was cloned and characterized using PCR-based techniques. The complete MAT locus of each species was cloned and annotated including complete and/or partial hypothetical genes flanking the idiomorphs. Analysis of the MAT locus organization indicated that S. chrysanthemi is likely homothallic with both MAT1-2-1 and MAT1-1-1 co-located within the idiomorph, and this was supported by production of the teleomorph in cultures of single-conidial-derived isolates. Sequencing of the MAT locus and flanking genes of S. tanaceti demonstrated that only a single MAT gene, MAT1-1-1, was located within this idiomorph and suggesting that S. tanaceti is heterothallic. MAT-specific PCR primers were developed and used to determine mating type of isolates sampled from diseased pyrethrum fields in Australia. These results indicated that only one mating type of S. tanaceti was present in Tasmania, Australia. The absence of a second mating type suggests that this species does not reproduce sexually in Tasmania, Australia and that ascospores are unlikely to be a source of inoculum for ray blight of pyrethrum. The MAT-specific PCR assay will be a valuable tool to distinguish mating types present among isolates of S. tanaceti, to monitor populations of S. tanaceti for the introduction of a second mating type and to differentiate S. tanaceti from S. chrysanthemi.

Phoma-like fungi on soybeans

Numerous coelomycetous fungi classified in Ascochyta, Phoma and Phyllosticta, and lately established and/or re-classified genera and species, namely Boeremia and Peyronellaea have been recorded from spots on leaves and pods of soybeans. These rarely observed pathogens are cosmopolitan, ubiquitous species on diseased and dead plant materials, and define frequently as weak or opportunistic parasites. Based on the Genealogical Concordance Phylogenetic Species Recognition, the authors summarize the re-evaluation of the taxonomic status of Phoma sojicola (syn. Ascochyta sojicola) and Phyllosticta sojicola. Inspite of the former delimitation of Ph. sojicola based on small differences in morphological features, it has proved to be identical to Peyronellaea pinodella (syn. Phoma pinodella). Similarly, it was also confirmed that Ph. sojicola was identical to Boeremia exigua var. exigua (syn. Phoma exigua var. exigua). The authors and co-workers contributed to the identification of Phoma-like fungi by combined conventional and molecular methods. Protein-encoding genes (TEF1 and β-tubulin) were successfully applied within the Phoma genus to infer phylogenetic relationships.

The heterothallic sugarbeet pathogen Cercospora beticola contains exon fragments of both MAT genes that are homogenized by concerted evolution

Dothideomycetes is one of the most ecologically diverse and economically important classes of fungi. Sexual reproduction in this group is governed by mating type (MAT) genes at the MAT1 locus. Self-sterile (heterothallic) species contain one of two genes at MAT1 (MAT1-1-1 or MAT1-2-1) and only isolates of opposite mating type are sexually compatible. In contrast, self-fertile (homothallic) species contain both MAT genes at MAT1. Knowledge of the reproductive capacities of plant pathogens are of particular interest because recombining populations tend to be more difficult to manage in agricultural settings. In this study, we sequenced MAT1 in the heterothallic Dothideomycete fungus Cercospora beticola to gain insight into the reproductive capabilities of this important plant pathogen. In addition to the expected MAT gene at MAT1, each isolate contained fragments of both MAT1-1-1 and MAT1-2-1 at ostensibly random loci across the genome. When MAT fragments from each locus were manually assembled, they reconstituted MAT1-1-1 and MAT1-2-1 exons with high identity, suggesting a retroposition event occurred in a homothallic ancestor in which both MAT genes were fused. The genome sequences of related taxa revealed that MAT gene fragment pattern of Cercospora zeae-maydis was analogous to C. beticola. In contrast, the genome of more distantly related Mycosphaerella graminicola did not contain MAT fragments. Although fragments occurred in syntenic regions of the C. beticola and C. zeae-maydis genomes, each MAT fragment was more closely related to the intact MAT gene of the same species. Taken together, these data suggest MAT genes fragmented after divergence of M. graminicola from the remaining taxa, and concerted evolution functioned to homogenize MAT fragments and MAT genes in each species.

Unequal recombination and evolution of the mating-type (MAT) loci in the pathogenic fungus Grosmannia clavigera and relatives

Sexual reproduction in fungi is regulated by the mating-type (MAT) locus where recombination is suppressed. We investigated the evolution of MAT loci in eight fungal species belonging to Grosmannia and Ophiostoma (Sordariomycetes, Ascomycota) that include conifer pathogens and beetle symbionts. The MAT1-2 idiomorph/allele was identified from the assembled and annotated Grosmannia clavigera genome, and the MAT locus is flanked by genes coding for cytoskeleton protein (SLA) and DNA lyase. The synteny of these genes is conserved and consistent with other members in Ascomycota. Using sequences from SLA and flanking regions, we characterized the MAT1-1 idiomorph from other isolates of G. clavigera and performed dotplot analysis between the two idiomorphs. Unexpectedly, the MAT1-2 idiomorph contains a truncated MAT1-1-1 gene upstream of the MAT1-2-1 gene that bears the high-mobility-group domain. The nucleotide and amino acid sequence of the truncated MAT1-1-1 gene is similar to its homologous copy in the MAT1-1 idiomorph in the opposite mating-type isolate, except that positive selection is acting on the truncated gene and the alpha(α)-box that encodes the transcription factor has been deleted. The MAT idiomorphs sharing identical gene organization were present in seven additional species in the Ophiostomatales, suggesting that the presence of truncated MAT1-1-1 gene is a general pattern in this order. We propose that an ancient unequal recombination event resulted in the ancestral MAT1-1-1 gene integrated into the MAT1-2 idiomorph and surviving as the truncated MAT1-1-1 genes. The α-box domain of MAT1-1-1 gene, located at the same MAT locus adjacent to the MAT1-2-1 gene, could have been removed by deletion after recombination due to mating signal interference. Our data confirmed a 1:1 MAT/sex ratio in two pathogen populations, and showed that all members of the Ophiostomatales studied here including those that were previously deemed asexual have the potential to reproduce sexually. This ability can potentially increase genetic variability and can enhance fitness in new, ecological niches.

The Sclerotinia sclerotiorum mating type locus (MAT) contains a 3.6-kb region that is inverted in every meiotic generation

Sclerotinia sclerotiorum is a fungal plant pathogen and the causal agent of lettuce drop, an economically important disease of California lettuce. The structure of the S. sclerotiorum mating type locus MAT has previously been reported and consists of two idiomorphs that are fused end-to-end as in other homothallics. We investigated the diversity of S. sclerotiorum MAT using a total of 283 isolates from multiple hosts and locations, and identified a novel MAT allele that differed by a 3.6-kb inversion and was designated Inv+, as opposed to the previously known S. sclerotiorum MAT that lacked the inversion and was Inv-. The inversion affected three of the four MAT genes: MAT1-2-1 and MAT1-2-4 were inverted and MAT1-1-1 was truncated at the 3'-end. Expression of MAT genes differed between Inv+ and Inv- isolates. In Inv+ isolates, only one of the three MAT1-2-1 transcript variants of Inv- isolates was detected, and the alpha1 domain of Inv+ MAT1-1-1 transcripts was truncated. Both Inv- and Inv+ isolates were self-fertile, and the inversion segregated in a 1∶1 ratio regardless of whether the parent was Inv- or Inv+. This suggested the involvement of a highly regulated process in maintaining equal proportions of Inv- and Inv+, likely associated with the sexual state. The MAT inversion region, defined as the 3.6-kb MAT inversion in Inv+ isolates and the homologous region of Inv- isolates, was flanked by a 250-bp inverted repeat on either side. The 250-bp inverted repeat was a partial MAT1-1-1 that through mediation of loop formation and crossing over, may be involved in the inversion process. Inv+ isolates were widespread, and in California and Nebraska constituted half of the isolates examined. We speculate that a similar inversion region may be involved in mating type switching in the filamentous ascomycetes Chromocrea spinulosa, Sclerotinia trifoliorum and in certain Ceratocystis species.

Evaluation of the potential for sexual reproduction in field populations of Cercospora beticola from USA

Cercospora leaf spot, caused by the hemibiotrophic fungal pathogen Cercospora beticola, is the most economically damaging foliar disease of sugarbeet worldwide. Although most C. beticola populations display characteristics reminiscent of sexual recombination, no teleomorph has been described. To assess whether populations in northern United States have characteristics consistent with sexual reproduction, 1024 isolates collected over a 3-y period were analyzed for frequency and distribution of mating type genes. After clone correction, an approximately equal distribution of mating types was found for each sampling year. Mating type frequency was also assessed in individual lesions. Lesions always consisted of isolates with a single mating type and microsatellite haplotype, but both mating types and up to five microsatellite haplotypes could be found on an individual leaf. The MAT1-1-1 and MAT1-2-1 genes were sequenced from 28 MAT1-1 and 28 MAT1-2 isolates, respectively. Three MAT1-1-1 nucleotide haplotypes were identified that encoded a single amino acid sequence. For MAT1-2-1, five nucleotide haplotypes were identified that encoded four protein variants. MAT1-1-1 and MAT1-2-1 gene expression analyses were conducted on plants inoculated with either or both mating types. MAT1-1-1 expression remained low, but MAT1-2-1 spiked during late stages of colonization. A segment of the MAT1-2-1 coding sequence was also found in MAT1-1 isolates. Taken together, these results suggest that C. beticola has the potential for sexual reproduction.