Comparative genome structure, secondary metabolite, and effector coding capacity across Cochliobolus pathogens

A barley MLA immune receptor is activated by a fungal nonribosomal peptide effector for disease susceptibility

The barley Mla locus contains functionally diversified genes that encode intracellular nucleotide-binding leucine-rich repeat receptors (NLRs) and confer strain-specific immunity to biotrophic and hemibiotrophic fungal pathogens. In this study, we isolated a barley gene Scs6, which is an allelic variant of Mla genes but confers susceptibility to the isolate ND90Pr (BsND90Pr) of the necrotrophic fungus Bipolaris sorokiniana. We generated Scs6 transgenic barley lines and showed that Scs6 is sufficient to confer susceptibility to BsND90Pr in barley genotypes naturally lacking the receptor. The Scs6-encoded NLR (SCS6) is activated by a nonribosomal peptide (NRP) effector produced by BsND90Pr to induce cell death in barley and Nicotiana benthamiana. Domain swaps between MLAs and SCS6 reveal that the SCS6 leucine-rich repeat domain is a specificity determinant for receptor activation by the NRP effector. Scs6 is maintained in both wild and domesticated barley populations. Our phylogenetic analysis suggests that Scs6 is a Hordeum-specific innovation. We infer that SCS6 is a bona fide immune receptor that is likely directly activated by the nonribosomal peptide effector of BsND90Pr for disease susceptibility in barley. Our study provides a stepping stone for the future development of synthetic NLR receptors in crops that are less vulnerable to modification by necrotrophic pathogens.

Regional diversity and leaf microbiome interactions of the fungal maize pathogen Exserohilum turcicum in Switzerland: A metagenomic analysis

The spread and adaptation of fungal plant pathogens in agroecosystems are facilitated by environmental homogeneity. Metagenomic sequencing of infected tissues allowed us to monitor eco-evolutionary dynamics and interactions between host, pathogen and plant microbiome. Exserohilum turcicum, the causal agent of northern corn leaf blight (NCLB) in maize, is distributed in multiple clonal lineages throughout Europe. To characterize regional pathogen diversity, we conducted metagenomic DNA sequencing on 241 infected leaf samples from the highly susceptible Swiss maize landrace Rheintaler Ribelmais, collected over 3 years (2016-2018) from an average of 14 agricultural farms within the Swiss Rhine Valley. All major European clonal lineages of E. turcicum were identified. Lineages differ by their mating types which indicates potential for sexual recombination and rapid evolution of new pathogen strains, although we found no evidence of recent recombination. The associated eukaryotic and prokaryotic leaf microbiome exhibited variation in taxonomic diversity between years and locations and is likely influenced by local weather conditions. A network analysis revealed distinct clusters of eukaryotic and prokaryotic taxa that correlates with the frequency of E. turcicum sequencing reads, suggesting causal interactions. Notably, the yeast genus Metschnikowia exhibited a strongly negative association with E. turcicum, supporting its known potential as biological control agent against fungal pathogens. Our findings show that metagenomic sequencing is a useful tool for analysing the role of environmental factors and potential pathogen-microbiome interactions in shaping pathogen dynamics and evolution, suggesting their potential for effective pathogen management strategies.

Recent advances in the population biology and management of maize foliar fungal pathogens Exserohilum turcicum, Cercospora zeina and Bipolaris maydis in Africa

Maize is the most widely cultivated and major security crop in sub-Saharan Africa. Three foliar diseases threaten maize production on the continent, namely northern leaf blight, gray leaf spot, and southern corn leaf blight. These are caused by the fungi Exserohilum turcicum, Cercospora zeina, and Bipolaris maydis, respectively. Yield losses of more than 10% can occur if these pathogens are diagnosed inaccurately or managed ineffectively. Here, we review recent advances in understanding the population biology and management of the three pathogens, which are present in Africa and thrive under similar environmental conditions during a single growing season. To effectively manage these pathogens, there is an increasing adoption of breeding for resistance at the small-scale level combined with cultural practices. Fungicide usage in African cropping systems is limited due to high costs and avoidance of chemical control. Currently, there is limited knowledge available on the population biology and genetics of these pathogens in Africa. The evolutionary potential of these pathogens to overcome host resistance has not been fully established. There is a need to conduct large-scale sampling of isolates to study their diversity and trace their migration patterns across the continent.

Complete genome of the Medicago anthracnose fungus, Colletotrichum destructivum, reveals a mini-chromosome-like region within a core chromosome

Colletotrichum destructivum (Cd) is a phytopathogenic fungus causing significant economic losses on forage legume crops (Medicago and Trifolium species) worldwide. To gain insights into the genetic basis of fungal virulence and host specificity, we sequenced the genome of an isolate from Medicago sativa using long-read (PacBio) technology. The resulting genome assembly has a total length of 51.7 Mb and comprises ten core chromosomes and two accessory chromosomes, all of which were sequenced from telomere to telomere. A total of 15, 631 gene models were predicted, including genes encoding potentially pathogenicity-related proteins such as candidate-secreted effectors (484), secondary metabolism key enzymes (110) and carbohydrate-active enzymes (619). Synteny analysis revealed extensive structural rearrangements in the genome of Cd relative to the closely related Brassicaceae pathogen, Colletotrichum higginsianum. In addition, a 1.2 Mb species-specific region was detected within the largest core chromosome of Cd that has all the characteristics of fungal accessory chromosomes (transposon-rich, gene-poor, distinct codon usage), providing evidence for exchange between these two genomic compartments. This region was also unique in having undergone extensive intra-chromosomal segmental duplications. Our findings provide insights into the evolution of accessory regions and possible mechanisms for generating genetic diversity in this asexual fungal pathogen.

Molecular Characterization and Expression Analysis of a Gene Encoding 3-Hydroxy-3-Methylglutaryl-CoA Reductase (HMGR) from Bipolaris eleusines, an Ophiobolin A-Producing Fungus

Ophibolin A, a fungal sesterterpene, exerts a pivotal influence in a diverse array of biological processes, encompassing herbicidal, bactericidal, fungicidal, and cytotoxic activities. Sixty genes associated with sesterterpene compound biosynthesis were obtained from Bipolaris eleusines via transcriptome sequencing, and those closely linked to ophiobolin A biosynthesis were subsequently filtered. A gene encoding 3-hydroxy-3-methylglutaryl-CoA reductase (HMGR) that catalyzes the first committed step of ophiobolin biosynthesis in the mevalonic acid (MVA) pathway was isolated and characterized using RACE (Rapid Amplification of cDNA Ends) technology from ophiobolin A-producing fungus, B. eleusines. The full-length cDNA of the B. eleusines HMGR gene (BeHMGR) was 3906 bp and contained a 3474 bp open reading frame (ORF) encoding 1157 amino acids. Sequence analysis revealed that deduced BeHMGR had high homology to the known HMGRs from Pyrenophora tritici-repentis and Leptosphaeria maculans. It had a calculated molecular mass of about 124.65 kDa and an isoelectric point (pI) of 6.90. It contained two putative HMG-CoA-binding motifs and two NADP(H)-binding motifs. Induced expression analysis of the BeHMGR gene by methyl jasmonate treatment using quantitative fluorescence PCR showed that it significantly elevated after 3 h of methyl jasmonate treatment, peaked at 6 h, and then gradually decreased. This demonstrates that BeHMGR gene expression is induced by methyl jasmonate.

Interactive transcriptome analyses of Northern Wild Rice (Zizania palustris L.) and Bipolaris oryzae show convoluted communications during the early stages of fungal brown spot development

Fungal diseases, caused mainly by Bipolaris spp., are past and current threats to Northern Wild Rice (NWR) grain production and germplasm preservation in both natural and cultivated settings. Genetic resistance against the pathogen is scarce. Toward expanding our understanding of the global gene communications of NWR and Bipolaris oryzae interaction, we designed an RNA sequencing study encompassing the first 12 h and 48 h of their encounter. NWR activated numerous plant recognition receptors after pathogen infection, followed by active transcriptional reprogramming of signaling mechanisms driven by Ca2+ and its sensors, mitogen-activated protein kinase cascades, activation of an oxidative burst, and phytohormone signaling-bound mechanisms. Several transcription factors associated with plant defense were found to be expressed. Importantly, evidence of diterpenoid phytoalexins, especially phytocassane biosynthesis, among expression of other defense genes was found. In B. oryzae, predicted genes associated with pathogenicity including secreted effectors that could target plant defense mechanisms were expressed. This study uncovered the early molecular communication between the NWR-B. oryzae pathosystem, which could guide selection for allele-specific genes to boost NWR defenses, and overall aid in the development of more efficient selection methods in NWR breeding through the use of the most virulent fungal isolates.

Transcriptome analysis of Bipolaris sorokiniana - Hordeum vulgare provides insights into mechanisms of host-pathogen interaction

Spot blotch disease incited by Bipolaris sorokiniana severely affects the cultivation of barley. The resistance to B. sorokiniana is quantitative in nature and its interaction with the host is highly complex which necessitates in-depth molecular analysis. Thus, the study aimed to conduct the transcriptome analysis to decipher the mechanisms and pathways involved in interactions between barley and B. sorokiniana in both the resistant (EC0328964) and susceptible (EC0578292) genotypes using the RNA Seq approach. In the resistant genotype, 6,283 genes of Hordeum vulgare were differentially expressed out of which 5,567 genes were upregulated and 716 genes were downregulated. 1,158 genes of Hordeum vulgare were differentially expressed in the susceptible genotype, out of which 654 genes were upregulated and 504 genes were downregulated. Several defense-related genes like resistant gene analogs (RGAs), disease resistance protein RPM1, pathogenesis-related protein PRB1-2-like, pathogenesis-related protein 1, thaumatin-like protein PWIR2 and defensin Tm-AMP-D1.2 were highly expressed exclusively in resistant genotype only. The pathways involved in the metabolism and biosynthesis of secondary metabolites were the most prominently represented pathways in both the resistant and susceptible genotypes. However, pathways involved in MAPK signaling, plant-pathogen interaction, and plant hormone signal transduction were highly enriched in resistant genotype. Further, a higher number of pathogenicity genes of B. sorokiniana was found in response to the susceptible genotype. The pathways encoding for metabolism, biosynthesis of secondary metabolites, ABC transporters, and ubiquitin-mediated proteolysis were highly expressed in susceptible genotype in response to the pathogen. 14 and 11 genes of B. sorokiniana were identified as candidate effectors from susceptible and resistant host backgrounds, respectively. This investigation will offer valuable insights in unraveling the complex mechanisms involved in barley- B. sorokiniana interaction.

Genome Mining Reveals the Biosynthesis of Sativene and Its Oxidative Conversion to seco-Sativene

Sativene (1) and seco-sativene are an important family of fungal sesquiterpenoids that feature unique tricyclo[4.4.0.01,7]decane and bicyclo[3.2.1]octane skeletons, respectively. Herein, we identify a three-enzyme cassette: SatA cyclizes farnesyl diphosphate (FPP) to form compound 1; CYP450 SatB catalyzes C14-C15 dihydroxylations and subsequent bond cleavage; and reductase SatC regioselectively reduces C14 aldehyde and mediates hemiacetal ring closure to generate prehelminthosporol (2). Our findings clarify the synthetic step of sativene and its oxidative transformation processes into seco-sativene.

Infection-specific transcriptional patterns of the maize pathogen Cochliobolus heterostrophus unravel genes involved in asexual development and virulence

Southern corn leaf blight (SCLB) caused by Cochliobolus heterostrophus is a destructive disease that threatens global maize (Zea mays) production. Despite many studies being conducted, very little is known about molecular processes employed by the pathogen during infection. There is a need to understand the fungal arms strategy and identify novel functional genes as targets for fungicide development. Transcriptome analysis based on RNA sequencing was carried out across conidia germination and host infection by C. heterostrophus. The present study revealed major changes in C. heterostrophus gene expression during host infection. Several differentially expressed genes (DEGs) induced during C. heterostrophus infection could be involved in the biosynthesis of secondary metabolites, peroxisome, energy metabolism, amino acid degradation and oxidative phosphorylation. In addition, histone acetyltransferase, secreted proteins, peroxisomal proteins, NADPH oxidase and transcription factors were selected for further functional validation. Here, we demonstrated that histone acetyltransferases (Hat2 and Rtt109), secreted proteins (Cel61A and Mep1), peroxisomal proteins (Pex11A and Pex14), NADPH oxidases (NoxA, NoxD and NoxR) and transcription factors (Crz1 and MtfA) play essential roles in C. heterostrophus conidiation, stress adaption and virulence. Taken together, our study revealed major changes in gene expression associated with C. heterostrophus infection and identified a diverse repertoire of genes critical for successful infection.

Population genomic analyses suggest recent dispersal events of the pathogen Cercospora zeina into East and Southern African maize cropping systems

A serious factor hampering global maize production is gray leaf spot disease. Cercospora zeina is one of the causative pathogens, but population genomics analysis of C. zeina is lacking. We conducted whole-genome Illumina sequencing of a representative set of 30 C. zeina isolates from Kenya and Uganda (East Africa) and Zambia, Zimbabwe, and South Africa (Southern Africa). Selection of the diverse set was based on microsatellite data from a larger collection of the pathogen. Pangenome analysis of the C. zeina isolates was done by (1) de novo assembly of the reads with SPAdes, (2) annotation with BRAKER, and (3) protein clustering with OrthoFinder. A published long-read assembly of C. zeina (CMW25467) from Zambia was included and annotated using the same pipeline. This analysis revealed 790 non-shared accessory and 10,677 shared core orthogroups (genes) between the 31 isolates. Accessory gene content was largely shared between isolates from all countries, with a few genes unique to populations from Southern Africa (32) or East Africa (6). There was a significantly higher proportion of effector genes in the accessory secretome (44%) compared to the core secretome (24%). PCA, ADMIXTURE, and phylogenetic analysis using a neighbor-net network indicated a population structure with a geographical subdivision between the East African isolates and the Southern African isolates, although gene flow was also evident. The small pangenome and partial population differentiation indicated recent dispersal of C. zeina into Africa, possibly from 2 regional founder populations, followed by recurrent gene flow owing to widespread maize production across sub-Saharan Africa.

Population structure and mixed reproductive strategies in Bipolaris maydis from single and multiple corn cultivars in Fujian Province, China

Bipolaris maydis is the pathogenic microorganism of southern corn leaf blight, a persistent biotic constraint responsible for substantial yield losses of corn worldwide. In the present study, 96 isolates from six representative fields growing single and multiple sweet corn cultivars in Pingnan, Fuqing, and Jian'ou in Fujian Province, which are characterized by different geographical characteristics and cropping patterns, were genetically analyzed using inter-simple sequence repeat (ISSR) markers to assess the impact of geographical origins and corn cultivars on B. maydis population differentiation. B. maydis isolates originated from diverse regions possessed higher genetic variety than those from single and multiple sweet corn cultivars. Phylogenetic analysis showed that the isolates from single and multiple sweet corn cultivars were randomly grouped into different clusters, with those from the same location tending to form clusters. A greater genetic differentiation among different geographical populations than between those from single and multiple sweet corn cultivars was observed by pairwise comparison. Hierarchical analysis indicated that among-population variation was higher when comparatively analyzed B. maydis populations from different locations than in those from single and multiple sweet corn cultivars. In conclusion, these results suggest that geographical origin acts a more considerable role in genetic differentiation of B. maydis than corn cultivar. Two divided genetic clusters were detected in the B. maydis populations from single and multiple sweet corn cultivars at the three locations in Fujian Province, with major genetic variation being derived within populations. The high haplotypic diversity and expected mating type ratio of 1:1 in combination with significant linkage disequilibrium suggested that a mixed reproductive strategy occurs in the B. maydis population in Fujian Province. This study will enrich the information on the role that geographical origins and corn cultivars play in the population structure of the pathogen as well as the reproductive strategies in B. maydis population in Fujian Province.

Prevalence and Importance of the Necrotrophic Effector Gene ToxA in Bipolaris sorokiniana Populations Collected from Spring Wheat and Barley

Bipolaris sorokiniana is a necrotrophic fungal pathogen that causes foliar and root diseases on wheat and barley. These diseases are common in all wheat- and barley-growing regions, with more severe outbreaks occurring under warm and humid conditions. B. sorokiniana can also infect a wide range of grass species in the family Poaceae and secrete ToxA, an important necrotrophic effector also identified other wheat leaf spotting pathogens. In this study, the prevalence and virulence role of ToxA were investigated in a collection of 278 B. sorokiniana isolates collected from spring wheat and barley in the Upper Midwest of the United States or other places, including 169 from wheat leaves, 75 from wheat roots, 30 from barley leaves, and 4 from wild quack grass leaves. ToxA was present in the isolates from wheat leaves, wheat roots, and wild grass leaves but was absent from isolates collected from barley leaves. Prevalence of ToxA in wheat leaf isolates (34.3%) was much higher than that in wheat root isolates (16%). Sequencing analysis revealed the presence of two haplotypes, with the majority being BsH2. All ToxA+ isolates produced the functional effector in liquid cultures. Pathogenicity assays revealed that ToxA+ isolates caused significantly more disease on spring wheat lines harboring Tsn1 than their tsn1 mutants, suggesting that the ToxA-Tsn1 interaction plays an important role in spot blotch development. This work confirms the importance of ToxA in B. sorokiniana populations infecting wheat and, thus, the need to eliminate Tsn1 from spring wheat cultivars to reduce susceptibility to spot blotch.

High-Quality Nuclear Genome and Mitogenome of Bipolaris sorokiniana LK93, a Devastating Pathogen Causing Wheat Root Rot

Bipolaris sorokiniana, one of the most devastating hemibiotrophic fungal pathogens, causes root rot, crown rot, leaf blotching, and black embryos of gramineous crops worldwide, posing a serious threat to global food security. However, the host-pathogen interaction mechanism between B. sorokiniana and wheat remains poorly understood. To facilitate related studies, we sequenced and assembled the genome of B. sorokiniana LK93. Nanopore long reads and next generation sequencing short reads were applied in the genome assembly, and the final 36.4-Mb genome assembly contains 16 contigs with the contig N50 of 2.3 Mb. Subsequently, we annotated 11,811 protein-coding genes. Of these, 10,620 were functional genes, 258 of which were identified as secretory proteins, including 211 predicted effectors. Additionally, the 111,581-bp mitogenome of LK93 was assembled and annotated. The LK93 genomes presented in this study will facilitate research in the B. sorokiniana-wheat pathosystem for better control of crop diseases. [Formula: see text] Copyright © 2023 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Unraveling the role of effector proteins in Bipolaris oryzae infecting North East Indian rice cultivars through time-course transcriptomics analysis

Bipolaris oryzae, causing brown spot disease in rice, is one of the neglected diseases reducing rice productivity. Limited knowledge is available on the genetics of host-pathogen interaction. Here, we used time-course transcriptome sequencing to elucidate the differential transcriptional responses of the pathogen genes in two contradictory infection-responsive rice hosts. Evaluation of transcriptome data showed similar regulation of fungal genes within susceptible (1733) and resistant (1846) hosts at an early stage however, in the later stage, the number was significantly higher in susceptible (2877) compared to resistant (1955) hosts. GO enrichment terms for upregulated genes showed a similar pattern in both the hosts at an early stage, but in the later stage terms related to degradation of carbohydrates, carbohydrate transport, and pathogenesis are enriched extensively within the susceptible host. Likewise, similar expression responses were observed with the secretory and effector proteins. Plant pathogenic homologs genes such as those involved in appressorium and conidia formation, host cell wall degradative enzymes, etc. were reported to be highly upregulated within the susceptible host. This study predicts the successful establishment of B. oryzae BO1 in both the host surfaces at an early stage, while disease progression only occurs in the susceptible host in later stage.

T-Toxin Virulence Genes: Unconnected Dots in a Sea of Repeats

In 1970, the Southern Corn Leaf Blight epidemic ravaged U.S. fields to great economic loss. The outbreak was caused by never-before-seen, supervirulent, Race T of the fungus Cochliobolus heterostrophus. The functional difference between Race T and O, the previously known, far less aggressive strain, is production of T-toxin, a host-selective polyketide. Supervirulence is associated with ~1 Mb of Race T-specific DNA; only a fraction encodes T-toxin biosynthetic genes (Tox1). Tox1 is genetically and physically complex, with unlinked loci (Tox1A, Tox1B) genetically inseparable from breakpoints of a Race O reciprocal translocation that generated hybrid Race T chromosomes. Previously, we identified 10 genes for T-toxin biosynthesis. Unfortunately, high-depth, short-read sequencing placed these genes on four small, unconnected scaffolds surrounded by repeated A+T rich sequence, concealing context. To sort out Tox1 topology and pinpoint the hypothetical Race O translocation breakpoints corresponding to Race T-specific insertions, we undertook PacBio long-read sequencing which revealed Tox1 gene arrangement and the breakpoints. Six Tox1A genes are arranged as three small islands in a Race T-specific sea (~634 kb) of repeats. Four Tox1B genes are linked, on a large loop of Race T-specific DNA (~210 kb). The race O breakpoints are short sequences of race O-specific DNA; corresponding positions in race T are large insertions of race T-specific, A+T rich DNA, often with similarity to transposable (predominantly Gypsy) elements. Nearby, are 'Voyager Starship' elements and DUF proteins. These elements may have facilitated Tox1 integration into progenitor Race O and promoted large scale recombination resulting in race T. IMPORTANCE In 1970 a corn disease epidemic ravaged fields in the United States to great economic loss. The outbreak was caused by a never-before seen, supervirulent strain of the fungal pathogen Cochliobolus heterostrophus. This was a plant disease epidemic, however, the current COVID-19 pandemic of humans is a stark reminder that novel, highly virulent, pathogens evolve with devastating consequences, no matter what the host-animal, plant, or other organism. Long read DNA sequencing technology allowed in depth structural comparisons between the sole, previously known, much less aggressive, version of the pathogen and the supervirulent version and revealed, in meticulous detail, the structure of the unique virulence-causing DNA. These data are foundational for future analysis of mechanisms of DNA acquisition from a foreign source.

The Chromosome-Scale Genomes of Exserohilum rostratum and Bipolaris zeicola Pathogenic Fungi Causing Rice Spikelet Rot Disease

Rice spikelet rot disease occurs mainly in the late stages of rice growth. Pathogenicity and biological characteristics of the pathogenic fungus and the infestation site have been the primary focus of research on the disease. To learn more about the disease, we performed whole-genome sequencing of Exserohilum rostratum and Bipolaris zeicola for predicting potentially pathogenic genes. The fungus B. zeicola was only recently identified in rice.We obtained 16 and 15 scaffolds down to the chromosome level for E. rostratum LWI and B. zeicola LWII, respectively. The length of LWI strain was approximately 34.05 Mb, and the G + C content of the whole genome was 50.56%. The length of the LWII strain was approximately 32.21 Mb, and the G + C content of the whole genome was 50.66%. After the prediction and annotation of E. rostratum LWI and B. zeicola LWII, we predicted that the LWI strain and LWII strain contain 8 and 13 potential pathogenic genes, respectively, which may be related to rice infection. These results improve our understanding of the genomes of E. rostratum and B. zeicola and update the genomic databases of these two species. It benefits subsequent studies on the mechanisms of E. rostratum and B. zeicola interactions with rice and helps to develop efficient control measures against rice spikelet rot disease.

Whole Genome Sequencing and Comparative Genomics of Indian Isolates of Wheat Spot Blotch Pathogen Bipolaris sorokiniana Reveals Expansion of Pathogenicity Gene Clusters

Spot blotch is a highly destructive disease in wheat caused by the fungal pathogen Bipolaris sorokiniana (teleomorph, Cochliobolus sativus). It is prevalent in warm and humid areas, including Africa, Asia, Latin America, and the USA. In the present study, twelve isolates of B. sorokiniana were collected from wheat fields in three different geographical locations in India. The pathogenicity of seven sporulating isolates was assessed on 'DDK 1025', a spot blotch-susceptible wheat variety under greenhouse conditions. The isolate 'D2' illustrated the highest virulence, followed by 'SI' and 'BS52'. These three isolates were sequenced using the Illumina HiSeq1000 platform. The estimated genome sizes of the isolates BS52, D2, and SI were 35.19 MB, 39.32 MB, and 32.76 MB, with GC contents of 48.48%, 50.43%, and 49.42%, respectively. The numbers of pathogenicity genes identified in BS52, D2, and SI isolates were 2015, 2476, and 2018, respectively. Notably, the isolate D2 exhibited a relatively larger genome with expanded arsenals of Biosynthetic Gene Clusters (BGCs), CAZymes, secretome, and pathogenicity genes, which could have contributed to its higher virulence among the tested isolates. This study provides the first comparative genome analysis of the Indian isolates of B. sorokiniana using whole genome sequencing.

Genetic diversity and structure of Bipolaris oryzae and Exserohilum rostratum populations causing brown spot of rice in Burkina Faso based on genotyping-by-sequencing

In recent years, Brown spot disease of rice (BSR) has been observed on leaves and seeds of rice in all rice-growing areas of Burkina Faso. Bipolaris oryzae and Exserohilum rostratum are the main fungal species isolated from BSR infected tissues and they are frequently observed in the same field. However, we are lacking information on the genetic diversity and population structure of these fungi in Burkina Faso. The mode of reproduction is also unknown. The genetic diversity of isolates of B. oryzae (n=61) and E. rostratum (n=151), collected from major rice-growing areas of Burkina Faso, was estimated using genotyping-by-sequencing (GBS). The mean values for nucleotide diversity (π) were 1.9 x10^-4 for B. oryzae and 4.8 x10^-4 for E. rostratum. There is no genetic differentiation between the geographical populations of each species. The analysis of molecular variance revealed that 89% and 94% of the genetic variances were within the populations of B. oryzae and E. rostratum, respectively. For each species, four genetic clusters were identified by two clustering methods (DAPC and sNMF). The distribution of these genetic groups was independent of the geographical origin of the isolates. Evidence of recombination was detected in the populations of B. oryzae and E. rostratum. For B. oryzae balanced mating type ratios were supporting sexual reproduction. For E. rostratum overrepresentation of MAT1-2 isolates (79%) suggested a predominant asexual reproduction. This study provides important information on the biology and genetics of the two major fungi causing brown spot disease of rice in Burkina Faso.

Comparative genomic analysis reveals cellulase plays an important role in the pathogenicity of Setosphaeria turcica f. sp. zeae

Setosphaeria turcica f. sp. zeae and S. turcica f. sp. sorghi, the two formae speciales of S. turcica, cause northern leaf blight disease of corn and sorghum, respectively, and often cause serious economic losses. They have obvious physiological differentiation and show complete host specificity. Host specificity is often closely related to pathogen virulence factors, including secreted protein effectors and secondary metabolites. Genomic sequencing can provide more information for understanding the virulence mechanisms of pathogens. However, the complete genomic sequence of S. turcica f. sp. sorghi has not yet been reported, and no comparative genomic information is available for the two formae speciales. In this study, S. turcica f. sp. zeae was predicted to have fewer secreted proteins, pathogen-host interaction (PHI) genes and carbohydrate-active enzymes (CAZys) than S. turcica f. sp. sorghi. Fifteen and 20 polyketide synthase (PKS) genes were identified in S. turcica f. sp. zeae and S. turcica f. sp. sorghi, respectively, which maintained high homology. There were eight functionally annotated effector protein-encoding genes specifically in S. turcica f. sp. zeae, among which the encoding gene StCEL2 of endo-1, 4-β-D-glucanase, an important component of cellulase, was significantly up-regulated during the interaction process. Finally, gluconolactone inhibited cellulase activity and decreased infection rate and pathogenicity, which indicates that cellulase is essential for maintaining virulence. These findings demonstrate that cellulase plays an important role in the pathogenicity of S. turcica f. sp. zeae. Our results also provide a theoretical basis for future research on the molecular mechanisms underlying the pathogenicity of the two formae speciales and for identifying any associated genes.

The Secreted Ribonuclease SRE1 Contributes to Setosphaeria turcica Virulence and Activates Plant Immunity

During the plant infection process, pathogens can secrete several effectors. Some of the effectors are well-known for their roles in regulating plant immunity and promoting successful pathogen colonization. However, there are few studies on the ribonuclease (RNase) effectors secreted by fungi. In the present study, we discovered a secretable RNase (SRE1) in the secretome of Setosphaeria turcica that was significantly upregulated during the early stages of S. turcica infection in maize. Knockdown of SRE1 significantly reduced the virulence of S. turcica. SRE1 can induce cell death in maize and Nicotiana benthamiana. However, unlike the conventional hypersensitive response (HR) caused by other effectors, SRE1 is not dependent on its signal peptide (SP) or plant receptor kinases (such as BAK1 and SOBIR1). SRE1-induced cell death depends upon its enzymatic activity and the N-terminal β-hairpin structure. SRE1 relies on its N-terminal β-hairpin structure to enter cells, and then degrades plant's RNA through its catalytic activity causing cytotoxic effects. Additionally, SRE1 enhances N. benthamiana's resistance to pathogenic fungi and oomycetes. In summary, SRE1 promotes the virulence of S. turcica, inducing plant cell death and activating plant immune responses.

The Interpretation of the Role of a Polyketide Synthase ClPKS18 in the Pathogenicity of Curvularia lunata

Plant pathogenic fungus Curvularia lunata (Wakker) Boedijn causes leaf spot diseases in several plants such as Oryza sativa, Sorghum bicolor (L.) Moench, and Capsicum frutescens. It has been spread worldwide, specifically in maize-growing regions. The polyketide synthase (PKS) plays a significant role in secondary metabolite production and its effect on virulence. The Clpks18 of C. lunata strongly correlated with its pathogenicity. The role of Clpks18 gene on the pathogenic activity of C. lunata remains unclear. Hence, in this study, we analyzed the importance of Clpks18 gene on the hyphae and conidial melanization and on the sporulation and hyphal growth. The deletion of Clpks18 gene reduced the production of methyl 5-(hydroxymethyl)furan-2-carboxylate toxin. The virulence of ΔClpks18 mutant was significantly reduced compared with the wild type. The metabolomics data revealed that (R)-(-)-mellein was a vital factor in the virulence of C. lunata. The (R)-(-)-mellein and the toxin produced by C. lunata were detected in the maize leaves during its infestation. In addition, the metabolomic analysis showed that the Clpks18 gene influences glycerolipid, non-ribosomal peptide biosynthesis, and its metabolism. This study demonstrates that the Clpks18 gene is important for the pathogenicity of C. lunata by influencing the complex metabolic network.

Genome-Wide Associations with Resistance to Bipolaris Leaf Spot (Bipolaris oryzae (Breda de Haan) Shoemaker) in a Northern Switchgrass Population (Panicum virgatum L.)

Switchgrass (Panicum virgatum L.), a northern native perennial grass, suffers from yield reduction from Bipolaris leaf spot caused by Bipolaris oryzae (Breda de Haan) Shoemaker. This study aimed to determine the resistant populations via multiple phenotyping approaches and identify potential resistance genes from genome-wide association studies (GWAS) in the switchgrass northern association panel. The disease resistance was evaluated from both natural (field evaluations in Ithaca, New York and Phillipsburg, Philadelphia) and artificial inoculations (detached leaf and leaf disk assays). The most resistant populations based on a combination of three phenotyping approaches—detached leaf, leaf disk, and mean from two locations—were ‘SW788’, ‘SW806’, ‘SW802’, ‘SW793’, ‘SW781’, ‘SW797’, ‘SW798’, ‘SW803’, ‘SW795’, ‘SW805’. The GWAS from the association panel showed 27 significant SNPs on 12 chromosomes: 1K, 2K, 2N, 3K, 3N, 4N, 5K, 5N, 6N, 7K, 7N, and 9N. These markers accumulatively explained the phenotypic variance of the resistance ranging from 3.28 to 26.52%. Within linkage disequilibrium of 20 kb, these SNP markers linked with the potential resistance genes included the genes encoding for NBS-LRR, PPR, cell-wall related proteins, homeostatic proteins, anti-apoptotic proteins, and ABC transporter.

A novel effector, CsSp1, from Bipolaris sorokiniana, is essential for colonization in wheat and is also involved in triggering host immunity

The hemibiotrophic pathogen Bipolaris sorokiniana causes root rot, leaf blotching, and black embryos in wheat and barley worldwide, resulting in significant yield and quality reductions. However, the mechanism underlying the host-pathogen interactions between B. sorokiniana and wheat or barley remains unknown. The B. sorokiniana genome encodes a large number of uncharacterized putative effector proteins. In this study, we identified a putative secreted protein, CsSp1, with a classic N-terminal signal peptide, that is induced during early infection. A split-marker approach was used to knock out CsSP1 in the Lankao 9-3 strain. Compared with the wild type, the deletion mutant ∆Cssp1 displayed less radial growth on potato dextrose agar plates and produced fewer spores, and complementary transformation completely restored the phenotype of the deletion mutant to that of the wild type. The pathogenicity of the deletion mutant in wheat was attenuated even though appressoria still penetrated the host. Additionally, the infectious hyphae in the deletion mutant became swollen and exhibited reduced growth in plant cells. The signal peptide of CsSp1 was functionally verified through a yeast YTK12 secretion system. Transient expression of CsSp1 in Nicotiana benthamiana inhibited lesion formation caused by Phytophthora capsici. Moreover, CsSp1 localized in the nucleus and cytoplasm of plant cells. In B. sorokiniana-infected wheat leaves, the salicylic acid-regulated genes TaPAL, TaPR1, and TaPR2 were down-regulated in the ∆Cssp1 strain compared with the wild-type strain under the same conditions. Therefore, CsSp1 is a virulence effector and is involved in triggering host immunity.

Comparative genome analysis of plant ascomycete fungal pathogens with different lifestyles reveals distinctive virulence strategies

BACKGROUND

Pathogens have evolved diverse lifestyles and adopted pivotal new roles in both natural ecosystems and human environments. However, the molecular mechanisms underlying their adaptation to new lifestyles are obscure. Comparative genomics was adopted to determine distinct strategies of plant ascomycete fungal pathogens with different lifestyles and to elucidate their distinctive virulence strategies.

RESULTS

We found that plant ascomycete biotrophs exhibited lower gene gain and loss events and loss of CAZyme-encoding genes involved in plant cell wall degradation and biosynthesis gene clusters for the production of secondary metabolites in the genome. Comparison with the candidate effectome detected distinctive variations between plant biotrophic pathogens and other groups (including human, necrotrophic and hemibiotrophic pathogens). The results revealed the biotroph-specific and lifestyle-conserved candidate effector families. These data have been configured in web-based genome browser applications for public display ( http://lab.malab.cn/soft/PFPG ). This resource allows researchers to profile the genome, proteome, secretome and effectome of plant fungal pathogens.

CONCLUSIONS

Our findings demonstrated different genome evolution strategies of plant fungal pathogens with different lifestyles and explored their lifestyle-conserved and specific candidate effectors. It will provide a new basis for discovering the novel effectors and their pathogenic mechanisms.

Novel factors contributing to fungal pathogenicity at early stages of Setosphaeria turcica infection

The fungal pathogen Setosphaeria turcica causes leaf blight on maize, which leads to considerable crop losses. However, how S. turcica establishes sustained systemic infection is largely unknown. Here, we report several novel factors contributing to S. turcica pathogenicity, identified using a genomic and transcriptional screen at different stages of S. turcica appressorium development. We identified two cytoskeleton regulators, SLM1 and SLM2, that are crucial for hypha and appressorium development. The SLM1 and SLM2 transcripts accumulated during germling stage but their levels were notably reduced at the appressorium stage. Deletion of SLM2 dramatically affected cell morphology, penetration ability, and pathogenicity. We also identified three different types of S. turcica glycosyl hydrolases that are critical for plant cell wall degradation. Their transcripts accumulated during the appressorium infection stage induced by cellophane and maize leaf. Most importantly, we characterized a novel and specific S. turcica effector, appressorium-coupled effector 1 (StACE1), whose expression is coupled to appressorium formation in S. turcica. This protein is required for maize infection and induces cell death on expression in Nicotiana benthamiana. These observations suggest that the phytopathogen S. turcica is primed in advance with multiple strategies for maize infection, which are coupled to appressorium formation at the early infection stages.

Aspergillus sp. A31 and Curvularia geniculata P1 mitigate mercury toxicity to Oryza sativa L

Aspergillus sp. A31 and Curvularia geniculata P1 are endophytes that colonize the roots of Aeschynomene fluminensis Vell. and Polygonum acuminatum Kunth. in humid environments contaminated with mercury. The two strains mitigated mercury toxicity and promoted Oryza sativa L growth. C. geniculata P1 stood out for increasing the host biomass by fourfold and reducing the negative effects of the metal on photosynthesis. Assembling and annotation of Aspergillus sp. A31 and C. geniculata P1 genomes resulted in 28.60 Mb (CG% 53.1; 10,312 coding DNA sequences) and 32.92 Mb (CG% 50.72; 8,692 coding DNA sequences), respectively. Twelve and 27 genomes of Curvularia/Bipolaris and Aspergillus were selected for phylogenomic analyzes, respectively. Phylogenetic analysis inferred the separation of species from the genus Curvularia and Bipolaris into different clades, and the separation of species from the genus Aspergillus into three clades; the species were distinguished by occupied niche. The genomes had essential gene clusters for the adaptation of microorganisms to high metal concentrations, such as proteins of the phytoquelatin-metal complex (GO: 0090423), metal ion binders (GO: 0046872), ABC transporters (GO: 0042626), ATPase transporters (GO: 0016887), and genes related to response to reactive oxygen species (GO: 0000302) and oxidative stress (GO: 0006979). The results reported here help to understand the unique regulatory mechanisms of mercury tolerance and plant development.

FunOrder: A robust and semi-automated method for the identification of essential biosynthetic genes through computational molecular co-evolution

Secondary metabolites (SMs) are a vast group of compounds with different structures and properties that have been utilized as drugs, food additives, dyes, and as monomers for novel plastics. In many cases, the biosynthesis of SMs is catalysed by enzymes whose corresponding genes are co-localized in the genome in biosynthetic gene clusters (BGCs). Notably, BGCs may contain so-called gap genes, that are not involved in the biosynthesis of the SM. Current genome mining tools can identify BGCs, but they have problems with distinguishing essential genes from gap genes. This can and must be done by expensive, laborious, and time-consuming comparative genomic approaches or transcriptome analyses. In this study, we developed a method that allows semi-automated identification of essential genes in a BGC based on co-evolution analysis. To this end, the protein sequences of a BGC are blasted against a suitable proteome database. For each protein, a phylogenetic tree is created. The trees are compared by treeKO to detect co-evolution. The results of this comparison are visualized in different output formats, which are compared visually. Our results suggest that co-evolution is commonly occurring within BGCs, albeit not all, and that especially those genes that encode for enzymes of the biosynthetic pathway are co-evolutionary linked and can be identified with FunOrder. In light of the growing number of genomic data available, this will contribute to the studies of BGCs in native hosts and facilitate heterologous expression in other organisms with the aim of the discovery of novel SMs.

The discovery and description of novel fungal secondary metabolites promises novel antibiotics, pharmaceuticals, and other useful compounds. A way to identify novel secondary metabolites is to express the corresponding genes in a suitable expression host. Consequently, a detailed knowledge or an accurate prediction of these genes is necessary. In fungi, the genes are co-localized in so-called biosynthetic gene clusters. Notably, the clusters may also contain genes that are not necessary for the biosynthesis of the secondary metabolites, so-called gap genes. We developed a method to detect co-evolved genes within the clusters and demonstrated that essential genes are co-evolving and can thus be differentiated from the gap genes. This adds an additional layer of information, which can support researchers with their decisions on which genes to study and express for the discovery of novel secondary metabolites.

The Conserved Colletotrichum spp. Effector Candidate CEC3 Induces Nuclear Expansion and Cell Death in Plants

Plant pathogens secrete proteins, known as effectors, that promote infection by manipulating host cells. Members of the phytopathogenic fungal genus Colletotrichum collectively have a broad host range and generally adopt a hemibiotrophic lifestyle that includes an initial biotrophic phase and a later necrotrophic phase. We hypothesized that Colletotrichum fungi use a set of conserved effectors during infection to support the two phases of their hemibiotrophic lifestyle. This study aimed to examine this hypothesis by identifying and characterizing conserved effectors among Colletotrichum fungi. Comparative genomic analyses using genomes of ascomycete fungi with different lifestyles identified seven effector candidates that are conserved across the genus Colletotrichum. Transient expression assays showed that one of these putative conserved effectors, CEC3, induces nuclear expansion and cell death in Nicotiana benthamiana, suggesting that CEC3 is involved in promoting host cell death during infection. Nuclear expansion and cell death induction were commonly observed in CEC3 homologs from four different Colletotrichum species that vary in host specificity. Thus, CEC3 proteins could represent a novel class of core effectors with functional conservation in the genus Colletotrichum.

Identification of non-ribosomal peptide synthetase in Ganoderma boninense Pat. that was expressed during the interaction with oil palm

Basal stem rot (BSR) of oil palm is a disastrous disease caused by a white-rot fungus Ganoderma boninense Pat. Non-ribosomal peptides (NRPs) synthesized by non-ribosomal peptide synthetases (NRPSs) are a group of secondary metabolites that act as fungal virulent factors during pathogenesis in the host. In this study, we aimed to isolate NRPS gene of G. boninense strain UPMGB001 and investigate the role of this gene during G. boninense-oil palm interaction. The isolated NRPS DNA fragment of 8322 bp was used to predict the putative peptide sequence of different domains and showed similarity with G. sinense (85%) at conserved motifs of three main NRPS domains. Phylogenetic analysis of NRPS peptide sequences demonstrated that NRPS of G. boninense belongs to the type VI siderophore family. The roots of 6-month-old oil palm seedlings were artificially inoculated for studying NRPS gene expression and disease severity in the greenhouse. The correlation between high disease severity (50%) and high expression (67-fold) of G. boninense NRPS gene at 4 months after inoculation and above indicated that this gene played a significant role in the advancement of BSR disease. Overall, these findings increase our knowledge on the gene structure of NRPS in G. boninense and its involvement in BSR pathogenesis as an effector gene.

Increased abundance of secreted hydrolytic enzymes and secondary metabolite gene clusters define the genomes of latent plant pathogens in the Botryosphaeriaceae

BACKGROUND

The Botryosphaeriaceae are important plant pathogens, but also have the ability to establish asymptomatic infections that persist for extended periods in a latent state. In this study, we used comparative genome analyses to shed light on the genetic basis of the interactions of these fungi with their plant hosts. For this purpose, we characterised secreted hydrolytic enzymes, secondary metabolite biosynthetic gene clusters and general trends in genomic architecture using all available Botryosphaeriaceae genomes, and selected Dothideomycetes genomes.

RESULTS

The Botryosphaeriaceae genomes were rich in carbohydrate-active enzymes (CAZymes), proteases, lipases and secondary metabolic biosynthetic gene clusters (BGCs) compared to other Dothideomycete genomes. The genomes of Botryosphaeria, Macrophomina, Lasiodiplodia and Neofusicoccum, in particular, had gene expansions of the major constituents of the secretome, notably CAZymes involved in plant cell wall degradation. The Botryosphaeriaceae genomes were shown to have moderate to high GC contents and most had low levels of repetitive DNA. The genomes were not compartmentalized based on gene and repeat densities, but genes of secreted enzymes were slightly more abundant in gene-sparse regions.

CONCLUSION

The abundance of secreted hydrolytic enzymes and secondary metabolite BGCs in the genomes of Botryosphaeria, Macrophomina, Lasiodiplodia, and Neofusicoccum were similar to those in necrotrophic plant pathogens and some endophytes of woody plants. The results provide a foundation for comparative genomic analyses and hypotheses to explore the mechanisms underlying Botryosphaeriaceae host-plant interactions.

Reference genome assembly for Australian Ascochyta lentis isolate Al4

Ascochyta lentis causes ascochyta blight in lentil (Lens culinaris Medik.) and yield loss can be as high as 50%. With careful agronomic management practices, fungicide use, and advances in breeding resistant lentil varieties, disease severity and impact to farmers have been largely controlled. However, evidence from major lentil producing countries, Canada and Australia, suggests that A. lentis isolates can change their virulence profile and level of aggressiveness over time and under different selection pressures. In this paper, we describe the first genome assembly for A. lentis for the Australian isolate Al4, through the integration of data from Illumina and PacBio SMRT sequencing. The Al4 reference genome assembly is almost 42 Mb in size and encodes 11,638 predicted genes. The Al4 genome comprises 21 full-length and gapless chromosomal contigs and two partial chromosome contigs each with one telomere. We predicted 31 secondary metabolite clusters, and 38 putative protein effectors, many of which were classified as having an unknown function. Comparison of A. lentis genome features with the recently published reference assembly for closely related A. rabiei show that genome synteny between these species is highly conserved. However, there are several translocations and inversions of genome sequence. The location of secondary metabolite clusters near transposable element and repeat-rich genomic regions was common for A. lentis as has been reported for other fungal plant pathogens.

Genome sequence of Monilinia vaccinii-corymbosi sheds light on mummy berry disease infection of blueberry and mating type

Mummy berry disease, caused by the fungal pathogen Monilinia vaccinii-corymbosi (Mvc), is one of the most economically important diseases of blueberries in North America. Mvc is capable of inducing two separate blighting stages during its life cycle. Infected fruits are rendered mummified and unmarketable. Genomic data for this pathogen is lacking, but could be useful in understanding the reproductive biology of Mvc and the mechanisms it deploys to facilitate host infection. In this study, PacBio sequencing and Hi-C interaction data were utilized to create a chromosome-scale reference genome for Mvc. The genome comprises nine chromosomes with a total length of 30 Mb, an N50 length of 4.06 Mb, and an average 413X sequence coverage. A total of 9399 gene models were predicted and annotated, and BUSCO analysis revealed that 98% of 1,438 searched conserved eukaryotic genes were present in the predicted gene set. Potential effectors were identified, and the mating-type (MAT) locus was characterized. Biotrophic effectors allow the pathogen to avoid recognition by the host plant and evade or mitigate host defense responses during the early stages of fruit infection. Following locule colonization, necrotizing effectors promote the mummification of host tissues. Potential biotrophic effectors utilized by Mvc include chorismate mutase for reducing host salicylate and necrotrophic effectors include necrosis-inducing proteins and hydrolytic enzymes for macerating host tissue. The MAT locus sequences indicate the potential for homothallism in the reference genome, but a deletion allele of the MAT locus, characterized in a second isolate, indicates heterothallism. Further research is needed to verify the roles of individual effectors in virulence and to determine the role of the MAT locus in outcrossing and population genotypic diversity.

The Role of Clt1-Regulated Xylan Metabolism in Melanin and Toxin Formation for the Pathogenicity of Curvularia lunata in Maize

We previously reported that the BTB (brica-brac, tramtrack, and broad) domain-containing protein Clt1 regulates melanin and toxin synthesis, conidiation, and pathogenicity in Curvularia lunata, but the interacting proteins and regulative mechanism of Clt1 are unclear. In this research, we identified two proteins, which respectively correspond to xylanase (Clxyn24) and acetyl xylan esterase (Claxe43) from C. lunata, that were regulated by Clt1. Yeast two-hybrid (Y2H) and bimolecular fluorescence complementation assays were conducted to verify the interaction of Clt1 with full-length Clxyn24 and Claxe43. Furthermore, the Y2H assay revealed that Clt1 physically interacted with Clxyn24 and Claxe43 through its BTB domain to degrade xylan, which was used as a carbon source for C. lunata growth. The utilization of xylan provides acetyl-CoA for the synthesis of melanin and toxin as well as energy and other intermediate metabolites for conidiation. Furthermore, transcriptome analysis revealed that PKS18 and its 13 flanking genes found clustered in a region spanning 57.89 kb on scaffold 9 of the C. lunata CX-3 genome were down-regulated in toxin production-deficient mutant T806, and this cluster is possibly responsible for toxin biosynthesis of C. lunata.[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.

Distribution of methionine sulfoxide reductases in fungi and conservation of the free-methionine-R-sulfoxide reductase in multicellular eukaryotes

Methionine, either as a free amino acid or included in proteins, can be oxidized into methionine sulfoxide (MetO), which exists as R and S diastereomers. Almost all characterized organisms possess thiol-oxidoreductases named methionine sulfoxide reductase (Msr) enzymes to reduce MetO back to Met. MsrA and MsrB reduce the S and R diastereomers of MetO, respectively, with strict stereospecificity and are found in almost all organisms. Another type of thiol-oxidoreductase, the free-methionine-R-sulfoxide reductase (fRMsr), identified so far in prokaryotes and a few unicellular eukaryotes, reduces the R MetO diastereomer of the free amino acid. Moreover, some bacteria possess molybdenum-containing enzymes that reduce MetO, either in the free or protein-bound forms. All these Msrs play important roles in the protection of organisms against oxidative stress. Fungi are heterotrophic eukaryotes that colonize all niches on Earth and play fundamental functions, in organic matter recycling, as symbionts, or as pathogens of numerous organisms. However, our knowledge on fungal Msrs is still limited. Here, we performed a survey of msr genes in almost 700 genomes across the fungal kingdom. We show that most fungi possess one gene coding for each type of methionine sulfoxide reductase: MsrA, MsrB, and fRMsr. However, several fungi living in anaerobic environments or as obligate intracellular parasites were devoid of msr genes. Sequence inspection and phylogenetic analyses allowed us to identify non-canonical sequences with potentially novel enzymatic properties. Finaly, we identified several ocurences of msr horizontal gene transfer from bacteria to fungi.

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.

Design of a diagnostic system based on molecular markers derived from the ascomycetes pan-genome analysis: The case of Fusarium dieback disease

A key factor to take actions against phytosanitary problems is the accurate and rapid detection of the causal agent. Here, we develop a molecular diagnostics system based on comparative genomics to easily identify fusariosis and specific pathogenic species as the Fusarium kuroshium, the symbiont of the ambrosia beetle Euwallaceae kuroshio Gomez and Hulcr which is responsible for Fusarium dieback disease in San Diego CA, USA. We performed a pan-genome analysis using sixty-three ascomycetes fungi species including phytopathogens and fungi associated with the ambrosia beetles. Pan-genome analysis revealed that 2,631 orthologue genes are only shared by Fusarium spp., and on average 3,941 (SD ± 1,418.6) are species-specific genes. These genes were used for PCR primer design and tested on DNA isolated from i) different strains of ascomycete species, ii) artificially infected avocado stems and iii) plant tissue of field-collected samples presumably infected. Our results let us propose a useful set of primers to either identify any species from Fusarium genus or, in a specific manner, species such as F. kuroshium, F. oxysporum, and F. graminearum. The results suggest that the molecular strategy employed in this study can be expanded to design primers against different types of pathogens responsible for provoking critical plant diseases.

Developmental Roles of the Hog1 Protein Phosphatases of the Maize Pathogen Cochliobolus heterostrophus

Protein phosphorylation cascades are universal in cell signaling. While kinome diversity allows specific phosphorylation events, relatively few phosphatases dephosphorylate key signaling proteins. Fungal mitogen activated protein kinases (MAPK), in contrast to their mammalian counterparts, often show detectable basal phosphorylation levels. Dephosphorylation, therefore, could act as a signal. In Cochliobolus heterostrophus, the Dothideomycete causing Southern corn leaf blight, ferulic acid (FA)—an abundant phenolic found in plant host cell walls—acts as a signal to rapidly dephosphorylate the stress-activated MAP kinase Hog1 (High Osmolarity Glycerol 1). In order to identify the protein phosphatases responsible, we constructed mutants in Hog1 phosphatases predicted from the genome by homology to yeast and other species. We found that Cochliobolus heterostrophus mutants lacking PtcB, a member of the PP2C family, exhibited altered growth, sporulation, and attenuated dephosphorylation in response to FA. The loss of the dual-specificity phosphatase CDC14 led to slow growth, decreased virulence, and attenuated dephosphorylation. Mutants in two predicted tyrosine phosphatase genes PTP1 and PTP2 showed normal development and virulence. Our results suggest that a network of phosphatases modulate Hog1’s dual phosphorylation levels. The mutants we constructed in this work provide a starting point to further unravel the signaling hierarchy by which exposure to FA leads to stress responses in the pathogen.

Validation of an in vitro system to trigger changes in the gene expression of effectors of Sclerotinia sclerotiorum

AIMS

Sclerotinia sclerotiorum, the causal agent of white mold, can infect several host species, including economically important crops. In this study, we propose and validate a new in vitro system able to mimic the conditions of interaction with the host and promote the induction of S. sclerotiorum effectors.

METHODS AND RESULTS

For culture media production, we selected three plant species, common bean (Phaseolus vulgaris L, cv. Requinte.), maize (Zea mays, cv. BRS1030) and beggarticks (Bidens pilosa). To validate this system as an in vitro inducer of effectors, the qRT-PCR technique was used to investigate the expression profile of some S. sclerotiorum effector genes in each growth medium at different times after inoculation.

CONCLUSION

The results obtained in this study provide a validation of a new method to study S. sclerotiorum during mimetic interaction with different hosts. Although leaf extract does not fully represent the plant environment, the presence of plant components in the culture medium seems to induce effector genes, mimicking in planta conditions. The use of MEVM is simpler than in planta growth, bypasses problems such as the amount of mycelium produced, as well as contamination of host cells during transcriptomic and proteomic analyses.

SIGNIFICANCE AND IMPACT OF THE STUDY

We have devised MEVM media as a model mimicking the interaction of S. sclerotiorum and its hosts and used it to evaluate in vitro expression of effectors normally expressed only in planta.

Sphinganine-Analog Mycotoxins (SAMs): Chemical Structures, Bioactivities, and Genetic Controls

Sphinganine-analog mycotoxins (SAMs) including fumonisins and A. alternata f. sp. Lycopersici (AAL) toxins are a group of related mycotoxins produced by plant pathogenic fungi in the Fusarium genus and in Alternaria alternata f. sp. Lycopersici, respectively. SAMs have shown diverse cytotoxicity and phytotoxicity, causing adverse impacts on plants, animals, and humans, and are a destructive force to crop production worldwide. This review summarizes the structural diversity of SAMs and encapsulates the relationships between their structures and biological activities. The toxicity of SAMs on plants and animals is mainly attributed to their inhibitory activity against the ceramide biosynthesis enzyme, influencing the sphingolipid metabolism and causing programmed cell death. We also reviewed the detoxification methods against SAMs and how plants develop resistance to SAMs. Genetic and evolutionary analyses revealed that the FUM (fumonisins biosynthetic) gene cluster was responsible for fumonisin biosynthesis in Fusarium spp. Sequence comparisons among species within the genus Fusarium suggested that mutations and multiple horizontal gene transfers involving the FUM gene cluster were responsible for the interspecific difference in fumonisin synthesis. We finish by describing methods for monitoring and quantifying SAMs in food and agricultural products.

Victorin, the host-selective cyclic peptide toxin from the oat pathogen Cochliobolus victoriae, is ribosomally encoded

The necrotrophic fungal pathogen Cochliobolus victoriae produces victorin, a host-selective toxin (HST) essential for pathogenicity to certain oat cultivars with resistance against crown rust. Victorin is a mixture of highly modified heterodetic cyclic hexapeptides, previously assumed to be synthesized by a nonribosomal peptide synthetase. Herein, we demonstrate that victorin is a member of the ribosomally synthesized and posttranslationally modified peptide (RiPP) family of natural products. Analysis of a newly generated long-read assembly of the C. victoriae genome revealed three copies of precursor peptide genes (vicA1-3) with variable numbers of "GLKLAF" core peptide repeats corresponding to the victorin peptide backbone. vicA1-3 are located in repeat-rich gene-sparse regions of the genome and are loosely clustered with putative victorin biosynthetic genes, which are supported by the discovery of compact gene clusters harboring corresponding homologs in two distantly related plant-associated Sordariomycete fungi. Deletion of at least one copy of vicA resulted in strongly diminished victorin production. Deletion of a gene encoding a DUF3328 protein (VicYb) abolished the production altogether, supporting its predicted role in oxidative cyclization of the core peptide. In addition, we uncovered a copper amine oxidase (CAO) encoded by vicK, in which its deletion led to the accumulation of new glycine-containing victorin derivatives. The role of VicK in oxidative deamination of the N-terminal glycyl moiety of the hexapeptides to the active glyoxylate forms was confirmed in vitro. This study finally unraveled the genetic and molecular bases for biosynthesis of one of the first discovered HSTs and expanded our understanding of underexplored fungal RiPPs.

Phylogenomic Analyses of Non-Dikarya Fungi Supports Horizontal Gene Transfer Driving Diversification of Secondary Metabolism in the Amphibian Gastrointestinal Symbiont, Basidiobolus

Research into secondary metabolism (SM) production by fungi has resulted in the discovery of diverse, biologically active compounds with significant medicinal applications. The fungi rich in SM production are taxonomically concentrated in the subkingdom Dikarya, which comprises the phyla Ascomycota and Basidiomycota. Here, we explore the potential for SM production in Mucoromycota and Zoopagomycota, two phyla of nonflagellated fungi that are not members of Dikarya, by predicting and identifying core genes and gene clusters involved in SM. The majority of non-Dikarya have few genes and gene clusters involved in SM production except for the amphibian gut symbionts in the genus Basidiobolus. Basidiobolus genomes exhibit an enrichment of SM genes involved in siderophore, surfactin-like, and terpene cyclase production, all these with evidence of constitutive gene expression. Gene expression and chemical assays also confirm that Basidiobolus has significant siderophore activity. The expansion of SMs in Basidiobolus are partially due to horizontal gene transfer from bacteria, likely as a consequence of its ecology as an amphibian gut endosymbiont.

Host-specificity factors in plant pathogenic fungi

Fortunately, no fungus can cause disease on all plant species, and although some plant-pathogenic fungi have quite a broad host range, most are highly limited in the range of plant species or even cultivars that they cause disease in. The mechanisms of host specificity have been extensively studied in many plant-pathogenic fungi, especially in fungal pathogens causing disease on economically important crops. Specifically, genes involved in host specificity have been identified during the last few decades. In this overview, we describe and discuss these host-specificity genes. These genes encode avirulence (Avr) proteins, proteinaceous host-specific toxins or secondary metabolites. We discuss the genomic context of these genes, their expression, polymorphism, horizontal transfer and involvement in pathogenesis.

The rice/maize pathogen Cochliobolus spp. infect and reproduce on Arabidopsis revealing differences in defensive phytohormone function between monocots and dicots

The fungal genus Cochliobolus describes necrotrophic pathogens that give rise to significant losses on rice, wheat, and maize. Revealing plant mechanisms of non-host resistance (NHR) against Cochliobolus will help to uncover strategies that can be exploited in engineered cereals. Therefore, we developed a heterogeneous pathosystem and studied the ability of Cochliobolus to infect dicotyledons. We report here that C. miyabeanus and C. heterostrophus infect Arabidopsis accessions and produce functional conidia, thereby demonstrating the ability to accept Brassica spp. as host plants. Some ecotypes exhibited a high susceptibility, whereas others hindered the necrotrophic disease progression of the Cochliobolus strains. Natural variation in NHR among the tested Arabidopsis accessions can advance the identification of genetic loci that prime the plant's defence repertoire. We found that applied phytotoxin-containing conidial fluid extracts of C. miyabeanus caused necrotic lesions on rice leaves but provoked only minor irritations on Arabidopsis. This result implies that C. miyabeanus phytotoxins are insufficiently adapted to promote dicot colonization, which corresponds to a retarded infection progression. Previous studies on rice demonstrated that ethylene (ET) promotes C. miyabeanus infection, whereas salicylic acid (SA) and jasmonic acid (JA) exert a minor function. However, in Arabidopsis, we revealed that the genetic disruption of the ET and JA signalling pathways compromises basal resistance against Cochliobolus, whereas SA biosynthesis mutants showed a reduced susceptibility. Our results refer to the synergistic action of ET/JA and indicate distinct defence systems between Arabidopsis and rice to confine Cochliobolus propagation. Moreover, this heterogeneous pathosystem may help to reveal mechanisms of NHR and associated defensive genes against Cochliobolus infection.

Biosynthetic gene clusters and the evolution of fungal chemodiversity

Covering: up to 2019Fungi produce a remarkable diversity of secondary metabolites: small, bioactive molecules not required for growth but which are essential to their ecological interactions with other organisms. Genes that participate in the same secondary metabolic pathway typically reside next to each other in fungal genomes and form biosynthetic gene clusters (BGCs). By synthesizing state-of-the-art knowledge on the evolution of BGCs in fungi, we propose that fungal chemodiversity stems from three molecular evolutionary processes involving BGCs: functional divergence, horizontal transfer, and de novo assembly. We provide examples of how these processes have contributed to the generation of fungal chemodiversity, discuss their relative importance, and outline major, outstanding questions in the field.