Genome and secretome analysis of the hemibiotrophic fungal pathogen, Moniliophthora roreri, which causes frosty pod rot disease of cacao: mechanisms of the biotrophic and necrotrophic phases

On the origin of bird's nest fungi: Phylogenomic analyses of fungi in the Nidulariaceae (Agaricales, Basidiomycota)

Nidulariaceae, also known as bird's nest fungi, is an understudied group of mushroom-forming fungi. The common name is derived from their nest-like morphology. Bird's nest fungi are ubiquitous wood decomposers or saprobes on dung. Recent studies showed that species in the Nidulariaceae form a monophyletic group with five sub-clades. However, phylogenetic relationships among genera and placement of Nidulariaceae are still unclear. We present phylogenomic analyses of bird's nest fungi and related Agaricales fungi to gain insight into the evolution of Nidulariaceae. A species tree with 17 newly generated genomes of bird's nest fungi and representatives from all major clades of Agaricales was constructed using 1044 single-copy genes to explore the intergeneric relationships and pinpoint the placement of Nidulariaceae within Agaricales. We corroborated the hypothesis that bird's nest fungi are sister to Squamanitaceae, which includes mushroom-shaped fungi with a stipe and pileus that are saprobes and mycoparasites. Lastly, stochastic character mapping of discrete traits on phylogenies (SIMMAP) suggests that the ancestor of bird's nest fungi likely possessed an evanescent, globose peridium without strings attaching to the spore packets (funiculi). This analysis suggests that the funiculus was gained twice and that the persistent, cupulate peridium form was gained at least four times and lost once. However, alternative coding schemes and datasets with a wider array of Agaricales produced conflicting results during ancestral state reconstruction, indicating that there is some uncertainty in the number of peridium transitions and that taxon sampling may significantly alter ancestral state reconstructions. Overall, our results suggest that several key morphological characters of Nidulariaceae have been subject to homoplasy.

A critical review on bioaerosols-dispersal of crop pathogenic microorganisms and their impact on crop yield

Bioaerosols are potential sources of pathogenic microorganisms that can cause devastating outbreaks of global crop diseases. Various microorganisms, insects and viroids are known to cause severe crop diseases impeding global agro-economy. Such losses threaten global food security, as it is estimated that almost 821 million people are underfed due to global crisis in food production. It is estimated that global population would reach 10 billion by 2050. Hence, it is imperative to substantially increase global food production to about 60% more than the existing levels. To meet the increasing demand, it is essential to control crop diseases and increase yield. Better understanding of the dispersive nature of bioaerosols, seasonal variations, regional diversity and load would enable in formulating improved strategies to control disease severity, onset and spread. Further, insights on regional and global bioaerosol composition and dissemination would help in predicting and preventing endemic and epidemic outbreaks of crop diseases. Advanced knowledge of the factors influencing disease onset and progress, mechanism of pathogen attachment and penetration, dispersal of pathogens, life cycle and the mode of infection, aid the development and implementation of species-specific and region-specific preventive strategies to control crop diseases. Intriguingly, development of R gene-mediated resistant varieties has shown promising results in controlling crop diseases. Forthcoming studies on the development of an appropriately stacked R gene with a wide range of resistance to crop diseases would enable proper management and yield. The article reviews various aspects of pathogenic bioaerosols, pathogen invasion and infestation, crop diseases and yield.

A repertoire of candidate effector proteins of the fungus Ceratocystis cacaofunesta

The genus Ceratocystis includes many phytopathogenic fungi that affect different plant species. One of these is Ceratocystis cacaofunesta, which is pathogenic to the cocoa tree and causes Ceratocystis wilt, a lethal disease for the crop. However, little is known about how this pathogen interacts with its host. The knowledge and identification of possible genes encoding effector proteins are essential to understanding this pathosystem. The present work aimed to predict genes that code effector proteins of C. cacaofunesta from a comparative analysis of the genomes of five Ceratocystis species available in databases. We performed a new genome annotation through an in-silico analysis. We analyzed the secretome and effectorome of C. cacaofunesta using the characteristics of the peptides, such as the presence of signal peptide for secretion, absence of transmembrane domain, and richness of cysteine residues. We identified 160 candidate effector proteins in the C. cacaofunesta proteome that could be classified as cytoplasmic (102) or apoplastic (58). Of the total number of candidate effector proteins, 146 were expressed, presenting an average of 206.56 transcripts per million. Our database was created using a robust bioinformatics strategy, followed by manual curation, generating information on pathogenicity-related genes involved in plant interactions, including CAZymes, hydrolases, lyases, and oxidoreductases. Comparing proteins already characterized as effectors in Sordariomycetes species revealed five groups of protein sequences homologous to C. cacaofunesta. These data provide a valuable resource for studying the infection mechanisms of these pathogens in their hosts.

Clonal reproduction of Moniliophthora roreri and the emergence of unique lineages with distinct genomes during range expansion

The basidiomycete Moniliophthora roreri causes frosty pod rot of cacao (Theobroma cacao) in the western hemisphere. Moniliophthora roreri is considered asexual and haploid throughout its hemibiotrophic life cycle. To understand the processes driving genome modification, using long-read sequencing technology, we sequenced and assembled 5 high-quality M. roreri genomes out of a collection of 99 isolates collected throughout the pathogen's range. We obtained chromosome-scale assemblies composed of 11 scaffolds. We used short-read technology to sequence the genomes of 22 similarly chosen isolates. Alignments among the 5 reference assemblies revealed inversions, translocations, and duplications between and within scaffolds. Isolates at the front of the pathogens' expanding range tend to share lineage-specific structural variants, as confirmed by short-read sequencing. We identified, for the first time, 3 new mating type A locus alleles (5 in total) and 1 new potential mating type B locus allele (3 in total). Currently, only 2 mating type combinations, A1B1 and A2B2, are known to exist outside of Colombia. A systematic survey of the M. roreri transcriptome across 2 isolates identified an expanded candidate effector pool and provided evidence that effector candidate genes unique to the Moniliophthoras are preferentially expressed during the biotrophic phase of disease. Notably, M. roreri isolates in Costa Rica carry a chromosome segment duplication that has doubled the associated gene complement and includes secreted proteins and candidate effectors. Clonal reproduction of the haploid M. roreri genome has allowed lineages with unique genome structures and compositions to dominate as it expands its range, displaying a significant founder effect.

Genomic and Pathogenicity Mechanisms of the Main Theobroma cacao L. Eukaryotic Pathogens: A Systematic Review

A set of diseases caused by fungi and oomycetes are responsible for large losses in annual world cocoa production. Managing the impact caused by these diseases is very complex because a common solution has yet to be found for different pathogens. In this context, the systematic knowledge of Theobroma cacao L. pathogens' molecular characteristics may help researchers understand the possibilities and limitations of cocoa disease management strategies. This work systematically organized and summarized the main findings of omics studies of T. cacao eukaryotic pathogens, focusing on the plant-pathogen interaction and production dynamics. Using the PRISMA protocol and a semiautomated process, we selected papers from the Scopus and Web of Science databases and collected data from the selected papers. From the initial 3169 studies, 149 were selected. The first author's affiliations were mostly from two countries, Brazil (55%) and the USA (22%). The most frequent genera were Moniliophthora (105 studies), Phytophthora (59 studies) and Ceratocystis (13 studies). The systematic review database includes papers reporting the whole-genome sequence from six cocoa pathogens and evidence of some necrosis-inducing-like proteins, which are common in T. cacao pathogen genomes. This review contributes to the knowledge about T. cacao diseases, providing an integrated discussion of T. cacao pathogens' molecular characteristics, common mechanisms of pathogenicity and how this knowledge is produced worldwide.

Characterization of the microbiota dynamics associated with Moniliophthora roreri, causal agent of cocoa frosty pod rot disease, reveals new viral species

Introduction

Theobroma cacao, the cocoa tree, is a target for pathogens, such as fungi from the genera Phytophthora, Moniliophthora, Colletotrichum, Ceratocystis, among others. Some cacao pathogens are restricted to specific regions of the world, such as the Cacao swollen shoot virus (CSSV) in West African countries, while others are expanding geographically, such as Moniliophthora roreri in the Americas. M. roreri is one of the most threatening cacao pathogens since it directly attacks the cacao pods driving a significant reduction in production, and therefore economic losses. Despite its importance, the knowledge about the microenvironment of this pathogen and the cocoa pods is still poorly characterized.

Methods

Herein we performed RNA sequencing of spores in differential stages of culture in a medium supplemented with cacao pod extract and mycelium collected of the susceptible variety ICT 7121 naturally infected by the pathogen to evaluate the diversity and transcriptional activity of microorganisms associated with the in vitro sporulation of M. roreri.

Results

Our data revealed a great variety of fungi and bacteria associated with M. roreri, with an exceptional diversity of individuals from the genus Trichoderma sp. Interestingly, the dynamics of microorganisms from different kingdoms varied proportionally, suggesting they are somehow affected by M. roreri culture time. We also identified three sequences similar to viral genomes from the Narnaviridae family, posteriorly confirmed by phylogenetic analysis as members of the genus Narnavirus. Screening of M. roreri public datasets indicated the virus sequences circulating in samples from Ecuador, suggesting a wide spread of these elements. Of note, we did not identify traces of the viral sequences in the M. roreri genome or DNA sequencing, restricting the possibility of these sequences representing endogenized elements.

Discussion

To the best of our knowledge, this is the first report of viruses infecting the fungus of the genus Moniliophthora and only the third description of viruses that are able to parasite elements from the Marasmiaceae family.

TMT-based quantitative membrane proteomics identified PRRs potentially involved in the perception of MSP1 in rice leaves

Pathogen-associated molecular patterns (PAMPs) play a key role in triggering PAMPs triggered immunity (PTI) in plants. In the case of the rice-Magnaporthe oryzae pathosystem, fewer PAMPs and their pattern recognition receptors (PRRs) have been characterized. Recently, a M. oryzae snodprot1 homolog protein (MSP1) has been identified that functions as PAMP and triggering the PTI responses in rice. However, the molecular mechanism underlying MSP1-induced PTI is currently elusive. Therefore, we generated MSP1 overexpressed transgenic lines of rice, and a tandem mass tag (TMT)-based quantitative membrane proteomic analysis was employed to decipher the potential MSP1-induced signaling in rice using total cytosolic as well as membrane protein fractions. This approach led to the identification of 8033 proteins of which 1826 were differentially modulated in response to overexpression of MSP1 and/or exogenous jasmonic acid treatment. Of these, 20 plasma membrane-localized receptor-like kinases (RLKs) showed increased abundance in MSP1 overexpression lines. Moreover, activation of proteins related to the protein degradation and modification, calcium signaling, redox, and MAPK signaling was observed in transgenic lines expressing MSP1 in the apoplast. Taken together, our results identified potential PRR candidates involved in MSP1 recognition and suggested the overview mechanism of the MSP1-induced PTI signaling in rice leaves. SIGNIFICANCE: In plants, recognition of pathogen pathogen-derived molecules, such as PAMPs, by plant plant-derived PRRs has an essential role for in the activation of PTI against pathogen invasion. Typically, PAMPs are recognized by plasma membrane (PM) localized PRRs, however, identifying the PM-localized PRR proteins is challenging due to their low abundance. In this study, we performed an integrated membrane protein enrichment by microsomal membrane extraction (MME) method and subsequent TMT-labeling-based quantitative proteomic analysis using MSP1 overexpressed rice. Based on these results, we successfully identified various intracellular and membrane membrane-localized proteins that participated in the MSP1-induced immune response and characterized the potential PM-localized PRR candidates in rice.

Mode of action of potential biocontrol agents against Fusarium species and Cochliobolus sativus

The antagonistic activities of three potential biocontrol agents (BCAs), Clonostachys miodochialis, C. rosea, and Minimedusa polyspora, were tested in vitro against the target fungi Fusarium acuminatum, F. avenaceum, F. equiseti, F. graminearum, and Cochliobolus sativus. In vitro dual-culture assay revealed that the BCAs were able to inhibit the mycelial growth of most of the target fungi, with the most common mode of action being mycoparasitism. After 5 d of incubation, depending on the host-parasite interface, biotrophic mycoparasitism was observed via attachment structures on F. avenaceum, F. equiseti, and C. sativus. Haustoria and appressoria were formed by C. miodochialis and M. polyspora on F. avenaceum and C. sativus, respectively. Clamp-like structures were also produced by both BCAs, depending on the host fungus. Clonostachys rosea formed only contact points on its hosts. Coiling structure was observed with C. rosea and M. polyspora occasionally in pure culture but more abundantly in the dual-culture assay. After an additional 3 to 5 d of incubation, the BCAs started damaging their host cells. Asexual fruiting bodies of the BCAs spread necrotrophically on F. acuminatum and F. graminearum and began destroying their mycelia after the initial 5 d of incubation. Furthermore, mycelial dissolution of F. acuminatum by an excreted substance was observed at a distance before direct contact with C. miodochialis. The diffuse metabolite assay revealed that the highest inhibition of the three BCAs was on C. sativus and their least effectiveness was observed on F. graminearum. Overall, these results provide evidence that C. miodochialis, C. rosea, and M. polyspora are potential candidates for biological control of the tested target fungi. This constitutes the first report that these three BCAs are able to establish an initial biotrophic relationship followed by a necrotrophic lifestyle. In addition, C. miodochialis has not previously been reported as a BCA.

Cross-Kingdom Gene Coexpression Analysis Using a Stemphylium botryosum-Lens ervoides System Revealed Plasticity of Intercommunication Between the Pathogen Secretome and the Host Immune Systems

Necrotrophic pathogens are responsible for significant declines in crop yield and quality worldwide. During the infection process, a pathogen releases a series of secretory proteins to counteract the plant immune system, and this interaction of pathogen and host molecules determines whether the pathogen will successfully invade the host plant tissues. In this study, we adopted co-transcriptomic approaches to analyze the Lens ervoides-Stemphylium botryosum system, with a focus on 1,216 fungal genes coding for secretory proteins and 8,810 disease-responsive genes of the host 48, 96, and 144 h postinoculation, captured in two F₉ recombinant inbred lines (RILs) displaying contrasting disease responses. By constructing in planta gene coexpression networks (GCNs) for S. botryosum, we found that the pathogen tended to co-upregulate genes regulating cell wall degradation enzymes, effectors, oxidoreductases, and peptidases to a much higher degree in the susceptible host LR-66-577 than in the resistant RIL LR-66-637, indicating that the promotion of these digestive enzymes and toxins increased S. botryosum virulence. Construction of cross-kingdom GCNs between pathogen and plant for the two RILs revealed that the co-upregulation of these fungal digestive enzymes and toxins simultaneously promoted a series of defense responses such as redox change, expression of membrane-related genes and serine/threonine kinase, and stress and disease responses in the susceptible RIL which was not observed in the resistant RIL, indicating that these activities exacerbated susceptibility to S. botryosum.[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.

Adaptive evolution of Moniliophthora PR-1 proteins towards its pathogenic lifestyle

BACKGROUND

Plant pathogenesis related-1 (PR-1) proteins belong to the CAP superfamily and have been characterized as markers of induced defense against pathogens. Moniliophthora perniciosa and Moniliophthora roreri are hemibiotrophic fungi that respectively cause the witches' broom disease and frosty pod rot in Theobroma cacao. Interestingly, a large number of plant PR-1-like genes are present in the genomes of both species and many are up-regulated during the biotrophic interaction. In this study, we investigated the evolution of PR-1 proteins from 22 genomes of Moniliophthora isolates and 16 other Agaricales species, performing genomic investigation, phylogenetic reconstruction, positive selection search and gene expression analysis.

RESULTS

Phylogenetic analysis revealed conserved PR-1 genes (PR-1a, b, d, j), shared by many Agaricales saprotrophic species, that have diversified in new PR-1 genes putatively related to pathogenicity in Moniliophthora (PR-1f, g, h, i), as well as in recent specialization cases within M. perniciosa biotypes (PR-1c, k, l) and M. roreri (PR-1n). PR-1 families in Moniliophthora with higher evolutionary rates exhibit induced expression in the biotrophic interaction and positive selection clues, supporting the hypothesis that these proteins accumulated adaptive changes in response to host-pathogen arms race. Furthermore, although previous work showed that MpPR-1 can detoxify plant antifungal compounds in yeast, we found that in the presence of eugenol M. perniciosa differentially expresses only MpPR-1e, k, d, of which two are not linked to pathogenicity, suggesting that detoxification might not be the main function of most MpPR-1.

CONCLUSIONS

Based on analyses of genomic and expression data, we provided evidence that the evolution of PR-1 in Moniliophthora was adaptive and potentially related to the emergence of the parasitic lifestyle in this genus. Additionally, we also discuss how fungal PR-1 proteins could have adapted from basal conserved functions to possible roles in fungal pathogenesis.

Transcriptional profile of oil palm pathogen, Ganoderma boninense, reveals activation of lignin degradation machinery and possible evasion of host immune response

BACKGROUND

The white-rot fungi in the genus Ganoderma interact with both living and dead angiosperm tree hosts. Two Ganoderma species, a North American taxon, G. zonatum and an Asian taxon, G. boninense, have primarily been found associated with live palm hosts. During the host plant colonization process, a massive transcriptional reorganization helps the fungus evade the host immune response and utilize plant cell wall polysaccharides.

RESULTS

A publicly available transcriptome of G. boninense - oil palm interaction was surveyed to profile transcripts that were differentially expressed in planta. Ten percent of the G. boninense transcript loci had altered expression as it colonized oil palm plants one-month post inoculation. Carbohydrate active enzymes (CAZymes), particularly those with a role in lignin degradation, and auxiliary enzymes that facilitate lignin modification, like cytochrome P450s and haloacid dehalogenases, were up-regulated in planta. Several lineage specific proteins and secreted proteins that lack known functional domains were also up-regulated in planta, but their role in the interaction could not be established. A slowdown in G. boninense respiration during the interaction can be inferred from the down-regulation of proteins involved in electron transport chain and mitochondrial biogenesis. Additionally, pathogenicity related genes and chitin degradation machinery were down-regulated during the interaction indicating G. boninense may be evading detection by the host immune system.

CONCLUSIONS

This analysis offers an overview of the dynamic processes at play in G. boninense - oil palm interaction and provides a framework to investigate biology of Ganoderma fungi across plantations and landscape.

The Role of Fungi in the Cocoa Production Chain and the Challenge of Climate Change

Background: The role of fungi in cocoa crops is mainly associated with plant diseases and contamination of harvest with unwanted metabolites such as mycotoxins that can reach the final consumer. However, in recent years there has been interest in discovering other existing interactions in the environment that may be beneficial, such as antagonism, commensalism, and the production of specific enzymes, among others. Scope and approach: This review summarizes the different fungi species involved in cocoa production and the cocoa supply chain. In particular, it examines the presence of fungal species during cultivation, harvest, fermentation, drying, and storage, emphasizing the factors that possibly influence their prevalence in the different stages of production and the health risks associated with the production of mycotoxins in the light of recent literature. Key findings and conclusion: Fungi associated with the cocoa production chain have many different roles. They have evolved in a varied range of ecosystems in close association with plants and various habitats, affecting nearly all the cocoa chain steps. Reports of the isolation of 60 genera of fungi were found, of which only 19 were involved in several stages. Although endophytic fungi can help control some diseases caused by pathogenic fungi, climate change, with increased rain and temperatures, together with intensified exchanges, can favour most of these fungal infections, and the presence of highly aggressive new fungal genotypes increasing the concern of mycotoxin production. For this reason, mitigation strategies need to be determined to prevent the spread of disease-causing fungi and preserve beneficial ones.

Mycena genomes resolve the evolution of fungal bioluminescence

Mushroom-forming fungi in the order Agaricales represent an independent origin of bioluminescence in the tree of life; yet the diversity, evolutionary history, and timing of the origin of fungal luciferases remain elusive. We sequenced the genomes and transcriptomes of five bonnet mushroom species (Mycena spp.), a diverse lineage comprising the majority of bioluminescent fungi. Two species with haploid genome assemblies ∼150 Mb are among the largest in Agaricales, and we found that a variety of repeats between Mycena species were differentially mediated by DNA methylation. We show that bioluminescence evolved in the last common ancestor of mycenoid and the marasmioid clade of Agaricales and was maintained through at least 160 million years of evolution. Analyses of synteny across genomes of bioluminescent species resolved how the luciferase cluster was derived by duplication and translocation, frequently rearranged and lost in most Mycena species, but conserved in the Armillaria lineage. Luciferase cluster members were coexpressed across developmental stages, with the highest expression in fruiting body caps and stipes, suggesting fruiting-related adaptive functions. Our results contribute to understanding a de novo origin of bioluminescence and the corresponding gene cluster in a diverse group of enigmatic fungal species.

De novo sequencing, assembly and functional annotation of Armillaria borealis genome

BACKGROUND

Massive forest decline has been observed almost everywhere as a result of negative anthropogenic and climatic effects, which can interact with pests, fungi and other phytopathogens and aggravate their effects. Climatic changes can weaken trees and make fungi, such as Armillaria more destructive. Armillaria borealis (Marxm. & Korhonen) is a fungus from the Physalacriaceae family (Basidiomycota) widely distributed in Eurasia, including Siberia and the Far East. Species from this genus cause the root white rot disease that weakens and often kills woody plants. However, little is known about ecological behavior and genetics of A. borealis. According to field research data, A. borealis is less pathogenic than A. ostoyae, and its aggressive behavior is quite rare. Mainly A. borealis behaves as a secondary pathogen killing trees already weakened by other factors. However, changing environment might cause unpredictable effects in fungus behavior.

RESULTS

The de novo genome assembly and annotation were performed for the A. borealis species for the first time and presented in this study. The A. borealis genome assembly contained ~ 68 Mbp and was comparable with ~ 60 and ~ 79.5 Mbp for the A. ostoyae and A. mellea genomes, respectively. The N50 for contigs equaled 50,544 bp. Functional annotation analysis revealed 21,969 protein coding genes and provided data for further comparative analysis. Repetitive sequences were also identified. The main focus for further study and comparative analysis will be on the enzymes and regulatory factors associated with pathogenicity.

CONCLUSIONS

Pathogenic fungi such as Armillaria are currently one of the main problems in forest conservation. A comprehensive study of these species and their pathogenicity is of great importance and needs good genomic resources. The assembled genome of A. borealis presented in this study is of sufficiently good quality for further detailed comparative study on the composition of enzymes in other Armillaria species. There is also a fundamental problem with the identification and classification of species of the Armillaria genus, where the study of repetitive sequences in the genomes of basidiomycetes and their comparative analysis will help us identify more accurately taxonomy of these species and reveal their evolutionary relationships.

Draft genomic sequence of Armillaria gallica 012m: insights into its symbiotic relationship with Gastrodia elata

Armillaria species (Basidiomycota, Physalacriaceae) are well known as plant pathogens related to serious root rot disease on various trees in forests and plantations. Interestingly, some Armillaria species are essential symbionts of the rare Chinese medicinal herb Gastrodia elata, a rootless and leafless orchid used for over 2000 years. In this work, an 87.3-M draft genome of Armillaria gallica 012m strain, which was symbiotic with G. elata, was assembled. The genome includes approximately 23.6% repetitive sequences and encodes 26,261 predicted genes. In comparison with other four genomes of Armillaria, the following gene families related to pathogenicity/saprophytic phase, including cytochrome P450 monooxygenases, carbohydrate-active enzyme AA3, and hydrophobins, were significantly contracted in A. gallica 012m. These characteristics may be beneficial for G. elata to get less injuries. The genome-guided analysis of differential expression between rhizomorph (RH) and vegetative mycelium (VM) showed that a total of 2549 genes were differentially expressed, including 632 downregulated genes and 1917 upregulated genes. In the RH, most differentially expressed genes (DEGs) related to pathogenicity were significantly upregulated. To further elucidate gene function, Gene Ontology enrichment analysis showed that the upregulated DEGs significantly grouped into monooxygenase activity, hydrolase activity, glucosidase activity, extracellular region, fungal cell wall, response to xenobiotic stimulus, response to toxic substance, etc. These phenomena indicate that RH had better infection ability than VM. The infection ability of RH may be beneficial for G. elata to obtain nutrition, because the rhizomorph constantly infected the nutritional stems of G. elata and formed the hyphae that can be digested by G. elata. These results clarified the characteristics of A. gallica 012m and the reason why the strain 012m can establish a symbiotic relationship with G. elata in some extent from the perspective of genomics.

Genome Assembly and Pathway Analysis of Edible Mushroom Agrocybe cylindracea

Agrocybe cylindracea, an edible mushroom, is widely cultivated for its abundance of nutrients and flavor, and many of its metabolites are reported to have beneficial roles, such as medicinal effects on tumors and chronical illnesses. However, the lack of genomic information has hindered further molecular studies on this fungus. Here, we present a genome assembly of A. cylindracea together with comparative genomics and pathway analyses of Agaricales species. The draft, generated from both next-generation sequencing (NGS) and single-molecule real-time (SMRT) sequencing platforms to overcome high genetic heterozygosity, is composed of a 56.5 Mb sequence and 15,384 predicted genes. This mushroom possesses a complex reproductive system, including tetrapolar heterothallic and secondary homothallic mechanisms, and harbors several hydrolases and peptidases for gradual and effective degradation of various carbon sources. Our pathway analysis reveals complex processes involved in the biosynthesis of polysaccharides and other active substances, including B vitamins, unsaturated fatty acids, and N-acetylglucosamine. RNA-seq data show that A. cylindracea stipes tend to synthesize carbohydrate for carbon sequestration and energy storage, whereas pilei are more active in carbon utilization and unsaturated fatty acid biosynthesis. These results reflect diverse functions of the two anatomical structures of the fruiting body. Our comprehensive genomic and transcriptomic data, as well as preliminary comparative analyses, provide insights into the molecular details of the medicinal effects in terms of active compounds and nutrient components.

Transcriptome analysis of the edible mushroom Lentinula edodes in response to blue light

Lentinula edodes is one of the most popular edible mushrooms worldwide and contains important medicinal components such as lentinan, ergosterol, and eritadenine. Mushroom metabolism is regulated by the mycelia and fruit body using light; however, in mushrooms, the underlying molecular mechanisms controlling this process as well as light-induced gene expression remain unclear. Therefore, in this study, we compared morphological changes and gene expression in the fruit bodies of L. edodes cultivated under blue light and continuous darkness. Our results showed that blue light primarily induced pileus growth (diameter and thickness) compared to dark cultivation. Alternatively, stipe length development was promoted by dark cultivation. We also performed RNAseq on L. edodes under the blue light/dark cultivation conditions. A total of 12,051 genes were used for aligning the Illumina raw reads and 762 genes that showed fold change cut-offs of >|2| and significance p-values of <0.05 were selected under blue light condition. Among the genes which showed two-fold changed genes, 221 were upregulated and 541 were downregulated. In order to identify blue light induced candidate genes, differentially expressed genes (DEGs) were selected according to 4-fold changes and validated by RT-PCR. We identified 8 upregulated genes under blue light condition, such as DDR48-heat shock protein, Fasciclin-domain-containing protein and carbohydrate esterase family 4 protein, FAD NAD-binding domain-containing protein that are involved in morphological development of primordium and embryonic muscle development, cell adhesion and affect the structure of cellulosic and non-cellulosic cell walls of fruit body development, and photoreceptor of blue light signaling for fruit body and pigment development, respectively. This study provides valuable insights into the molecular mechanisms underlying the role of blue light in mushroom growth and development and can thus contribute to breeding programs to improve mushroom cultivation.

Plant Cells under Attack: Unconventional Endomembrane Trafficking during Plant Defense

Since plants lack specialized immune cells, each cell has to defend itself independently against a plethora of different pathogens. Therefore, successful plant defense strongly relies on precise and efficient regulation of intracellular processes in every single cell. Smooth trafficking within the plant endomembrane is a prerequisite for a diverse set of immune responses. Pathogen recognition, signaling into the nucleus, cell wall enforcement, secretion of antimicrobial proteins and compounds, as well as generation of reactive oxygen species, all heavily depend on vesicle transport. In contrast, pathogens have developed a variety of different means to manipulate vesicle trafficking to prevent detection or to inhibit specific plant responses. Intriguingly, the plant endomembrane system exhibits remarkable plasticity upon pathogen attack. Unconventional trafficking pathways such as the formation of endoplasmic reticulum (ER) bodies or fusion of the vacuole with the plasma membrane are initiated and enforced as the counteraction. Here, we review the recent findings on unconventional and defense-induced trafficking pathways as the plant´s measures in response to pathogen attack. In addition, we describe the endomembrane system manipulations by different pathogens, with a focus on tethering and fusion events during vesicle trafficking.

Time-Course RNAseq Reveals Exserohilum turcicum Effectors and Pathogenicity Determinants

Exserohilum turcicum (sexual stage Setosphaeria turcica) is the hemibiotrophic causal agent of northern leaf blight of maize and sorghum. This study aimed to identify the genes involved in host colonization during the biotrophic and necrotrophic phases of infection. It also aimed to identify race-specific differences in gene expression. RNAseq of maize seedlings inoculated with a race 13N or 23N E. turcicum isolate was conducted before inoculation and at 2, 5, 7, and 13 days post-inoculation (dpi). Biological replicates were pooled per time point for each race and sequenced. A bioinformatics pipeline was used to identify candidate effectors, and expression was validated for selected candidates. Fungal biomass was positively correlated with the percentages of E. turcicum reads mapped, which were low at early time points (2-7 dpi) with a significant increase at 13 dpi, indicating a lifestyle switch from biotrophy to necrotrophy between 7 and 13 dpi. AVRHt1 is the putative E. turcicum effector recognized by the maize resistance gene Ht1. Consistent with this, AVRHt1 was expressed in planta by race 23N, but transcripts were absent in race 13N. In addition, specific transposable elements were expressed in 23N only. Genes encoding the virulence-associated peptidases leupeptin-inhibiting protein 1 and fungalysin were expressed in planta. Transcriptional profiles of genes involved in secondary metabolite synthesis or cell wall degradation revealed the importance of these genes during late stages of infection (13 dpi). A total of 346 expressed candidate effectors were identified, including Ecp6 and proteins similar to the secreted in xylem (SIX) effectors common to formae speciales of Fusarium oxysporum, SIX13 and SIX5. Expression profiling of Ecp6 and SIX13-like indicated a peak in expression at 5 and 7 dpi compared to 2 and 13 dpi. Sequencing of SIX13-like from diverse isolates of E. turcicum revealed host-specific polymorphisms that were mostly non-synonymous, resulting in two groups of SIX13-like proteins that corresponded to the maize or sorghum origin of each isolate. This study suggests putative mechanisms whereby E. turcicum causes disease. Identification of the candidate effector SIX13-like is consistent with the infection mode of E. turcicum through the xylem of susceptible hosts.

Phyllospheric Microbiomes: Diversity, Ecological Significance, and Biotechnological Applications

The phyllosphere referred to the total aerial plant surfaces (above-ground portions), as habitat for microorganisms. Microorganisms establish compositionally complex communities on the leaf surface. The microbiome of phyllosphere is rich in diversity of bacteria, fungi, actinomycetes, cyanobacteria, and viruses. The diversity, dispersal, and community development on the leaf surface are based on the physiochemistry, environment, and also the immunity of the host plant. A colonization process is an important event where both the microbe and the host plant have been benefited. Microbes commonly established either epiphytic or endophytic mode of life cycle on phyllosphere environment, which helps the host plant and functional communication with the surrounding environment. To the scientific advancement, several molecular techniques like metagenomics and metaproteomics have been used to study and understand the physiology and functional relationship of microbes to the host and its environment. Based on the available information, this chapter describes the basic understanding of microbiome in leaf structure and physiology, microbial interactions, especially bacteria, fungi, and actinomycetes, and their adaptation in the phyllosphere environment. Further, the detailed information related to the importance of the microbiome in phyllosphere to the host plant and their environment has been analyzed. Besides, biopotentials of the phyllosphere microbiome have been reviewed.

Comparative genomic and transcriptomic analyses reveal different pathogenicity-related genes among three eucalyptus fungal pathogens

Ceratocystis fimbriata is an important plant pathogen known to cause Ceratocystis Wilt (CW), a prevalent fungal disease known to affect Eucalyptus spp. plantations in Brazil. To better understand the molecular mechanisms related to pathogenicity in eucalyptus, we generated a high-quality assembly and annotation of the Ce. fimbriata LPF1912 isolate (LPF1912) genome, as well as the first transcriptome of LPF1912 from 16 eucalyptus clones at three infection incubation periods (12, 18, and 24 h). The LPF1912 genome assembly contains 805 scaffolds, totaling 31.8 Mb, with 43% of the genome estimated to be coding sequence comprised of 7,390 protein-coding genes of which 626 (8.5%) were classified as secreted proteins, 120 ribosomal RNAs, and 532 transfer RNAs. Comparative genomic analysis among three eucalyptus fungal pathogens (Ce. fimbriata, Ce. eucalypticola, and Calonectria pseudoreteaudii), showed high similarity in the proteome (21.81%) and secretome (52.01%) of LPF1912 and Ce. eucalypticola. GO annotation of pathogenicity-related genes of LPF1912 and Ce. eucalypticola, revealed enrichment in cell wall degrading enzymes (CWDEs), and lipid/cutin metabolism for Ca. pseudoreteaudii. Additionally, a transcriptome analysis between resistant and susceptible eucalyptus clones to CW infection indicated that a majority (11) of LPF1912 differentially expressed genes had GO terms associated with enzymatic functions, such as the polygalacturonase gene family, confirming the crucial role of CWDEs for Ce. fimbriata pathogenicity. Finally, our genomic and transcriptomic analysis approach provides a better understanding of the mechanisms involved in Ce. fimbriata pathogenesis, as well as a framework for further studies.

Directing Trophic Divergence in Plant-Pathogen Interactions: Antagonistic Phytohormones With NO Doubt?

A fundamental process culminating in the mechanisms of plant-pathogen interactions is the regulation of trophic divergence into biotrophic, hemibiotrophic, and necrotrophic interactions. Plant hormones, of almost all types, play significant roles in this regulatory apparatus. In plant-pathogen interactions, two classical mechanisms underlying hormone-dependent trophic divergence are long recognized. While salicylic acid dominates in the execution of host defense response against biotrophic and early-stage hemibiotrophic pathogens, jasmonic acid, and ethylene are key players facilitating host defense response against necrotrophic and later-stage hemibiotrophic pathogens. Evidence increasingly suggests that trophic divergence appears to be modulated by more complex signaling networks. Acting antagonistically or agonistically, other hormones such as auxins, cytokinins, abscisic acid, gibberellins, brassinosteroids, and strigolactones, as well as nitric oxide, are emerging candidates in the regulation of trophic divergence. In this review, the latest advances in the dynamic regulation of trophic divergence are summarized, emphasizing common and contrasting hormonal and nitric oxide signaling strategies deployed in plant-pathogen interactions.

Insight into plant cell wall degradation and pathogenesis of Ganoderma boninense via comparative genome analysis

Background

G. boninense is a hemibiotrophic fungus that infects oil palms (Elaeis guineensis Jacq.) causing basal stem rot (BSR) disease and consequent massive economic losses to the oil palm industry. The pathogenicity of this white-rot fungus has been associated with cell wall degrading enzymes (CWDEs) released during saprophytic and necrotrophic stage of infection of the oil palm host. However, there is a lack of information available on the essentiality of CWDEs in wood-decaying process and pathogenesis of this oil palm pathogen especially at molecular and genome levels.

Methods

In this study, comparative genome analysis was carried out using the G. boninense NJ3 genome to identify and characterize carbohydrate-active enzyme (CAZymes) including CWDE in the fungal genome. Augustus pipeline was employed for gene identification in G. boninense NJ3 and the produced protein sequences were analyzed via dbCAN pipeline and PhiBase 4.5 database annotation for CAZymes and plant-host interaction (PHI) gene analysis, respectively. Comparison of CAZymes from G. boninense NJ3 was made against G. lucidum, a well-studied model Ganoderma sp. and five selected pathogenic fungi for CAZymes characterization. Functional annotation of PHI genes was carried out using Web Gene Ontology Annotation Plot (WEGO) and was used for selecting candidate PHI genes related to cell wall degradation of G. boninense NJ3.

Results

G. boninense was enriched with CAZymes and CWDEs in a similar fashion to G. lucidum that corroborate with the lignocellulolytic abilities of both closely-related fungal strains. The role of polysaccharide and cell wall degrading enzymes in the hemibiotrophic mode of infection of G. boninense was investigated by analyzing the fungal CAZymes with necrotrophic Armillaria solidipes, A. mellea, biotrophic Ustilago maydis, Melampsora larici-populina and hemibiotrophic Moniliophthora perniciosa. Profiles of the selected pathogenic fungi demonstrated that necrotizing pathogens including G. boninense NJ3 exhibited an extensive set of CAZymes as compared to the more CAZymes-limited biotrophic pathogens. Following PHI analysis, several candidate genes including polygalacturonase, endo β-1,3-xylanase, β-glucanase and laccase were identified as potential CWDEs that contribute to the plant host interaction and pathogenesis.

Discussion

This study employed bioinformatics tools for providing a greater understanding of the biological mechanisms underlying the production of CAZymes in G. boninense NJ3. Identification and profiling of the fungal polysaccharide- and lignocellulosic-degrading enzymes would further facilitate in elucidating the infection mechanisms through the production of CWDEs by G. boninense. Identification of CAZymes and CWDE-related PHI genes in G. boninense would serve as the basis for functional studies of genes associated with the fungal virulence and pathogenicity using systems biology and genetic engineering approaches.

Genome and transcriptome analysis of the latent pathogen Lasiodiplodia theobromae, an emerging threat to the cacao industry

Lasiodiplodia theobromae (Pat.) Griffon & Maubl., a member of the family Botryosphaeriaceae, is becoming a significant threat to crops and woody plants in many parts of the world, including the major cacao growing areas. While attempting to isolate Ceratobasidium theobromae, a causal agent of vascular streak dieback (VSD), from symptomatic cacao stems, 74% of isolated fungi were Lasiodiplodia spp. Sequence-based identification of 52 putative isolates of L. theobromae indicated that diverse species of Lasiodiplodia were associated with cacao in the studied areas, and the isolates showed variation in aggressiveness when assayed using cacao leaf discs. The present study reports a 43.75 Mb de novo assembled genome of an isolate of L. theobromae from cacao. Ab initio gene prediction generated 13 061 protein-coding genes, of which 2862 are unique to L. theobromae, when compared with other closely related Botryosphaeriaceae. Transcriptome analysis revealed that 11 860 predicted genes were transcriptionally active and 1255 were more highly expressed in planta compared with cultured mycelia. The predicted genes differentially expressed during infection were mainly those involved in carbohydrate, pectin, and lignin catabolism, cytochrome P450, necrosis-inducing proteins, and putative effectors. These findings significantly expand our knowledge of the genome of L. theobromae and the genes involved in virulence and pathogenicity.

Fungal evolution: major ecological adaptations and evolutionary transitions

Fungi are a highly diverse group of heterotrophic eukaryotes characterized by the absence of phagotrophy and the presence of a chitinous cell wall. While unicellular fungi are far from rare, part of the evolutionary success of the group resides in their ability to grow indefinitely as a cylindrical multinucleated cell (hypha). Armed with these morphological traits and with an extremely high metabolical diversity, fungi have conquered numerous ecological niches and have shaped a whole world of interactions with other living organisms. Herein we survey the main evolutionary and ecological processes that have guided fungal diversity. We will first review the ecology and evolution of the zoosporic lineages and the process of terrestrialization, as one of the major evolutionary transitions in this kingdom. Several plausible scenarios have been proposed for fungal terrestralization and we here propose a new scenario, which considers icy environments as a transitory niche between water and emerged land. We then focus on exploring the main ecological relationships of Fungi with other organisms (other fungi, protozoans, animals and plants), as well as the origin of adaptations to certain specialized ecological niches within the group (lichens, black fungi and yeasts). Throughout this review we use an evolutionary and comparative-genomics perspective to understand fungal ecological diversity. Finally, we highlight the importance of genome-enabled inferences to envision plausible narratives and scenarios for important transitions.

Chocolate Under Threat from Old and New Cacao Diseases

Theobroma cacao, the source of chocolate, is affected by destructive diseases wherever it is grown. Some diseases are endemic; however, as cacao was disseminated from the Amazon rain forest to new cultivation sites it encountered new pathogens. Two well-established diseases cause the greatest losses: black pod rot, caused by several species of Phytophthora, and witches' broom of cacao, caused by Moniliophthora perniciosa. Phytophthora megakarya causes the severest damage in the main cacao producing countries in West Africa, while P. palmivora causes significant losses globally. M. perniciosa is related to a sister basidiomycete species, M. roreri which causes frosty pod rot. These Moniliophthora species only occur in South and Central America, where they have significantly limited production since the beginnings of cacao cultivation. The basidiomycete Ceratobasidium theobromae causing vascular-streak dieback occurs only in South-East Asia and remains poorly understood. Cacao swollen shoot disease caused by Cacao swollen shoot virus is rapidly spreading in West Africa. This review presents contemporary research on the biology, taxonomy and genomics of what are often new-encounter pathogens, as well as the management of the diseases they cause.

Secretome Profiling Reveals Virulence-Associated Proteins of Fusarium proliferatum during Interaction with Banana Fruit

Secreted proteins are vital for the pathogenicity of many fungi through manipulating their hosts for efficient colonization. Fusarium proliferatum is a phytopathogenic fungus infecting many crops, vegetables, and fruit, including banana fruit. To access the proteins involved in pathogen–host interaction, we used label-free quantitative proteomics technology to comparatively analyze the secretomes of F. proliferatum cultured with and without banana peel in Czapek’s broth medium. By analyzing the secretomes of F. proliferatum, we have identified 105 proteins with 40 exclusively secreted and 65 increased in abundance in response to a banana peel. These proteins were involved in the promotion of invasion of banana fruit, and they were mainly categorized into virulence factors, cell wall degradation, metabolic process, response to stress, regulation, and another unknown biological process. The expressions of corresponding genes confirmed the existence of these secreted proteins in the banana peel. Furthermore, expression pattern suggested variable roles for these genes at different infection stages. This study expanded the current database of F. proliferatum secreted proteins which might be involved in the infection strategy of this fungus. Additionally, this study warranted the further attention of some secreted proteins that might initiate infection of F. proliferatum on banana fruit.

The biotrophy-associated secreted protein 4 (BAS4) participates in the transition of Magnaporthe oryzae from the biotrophic to the necrotrophic phase

The physiological and metabolic processes of host plants are manipulated and remodeled by phytopathogenic fungi during infection, revealed obvious signs of biotrophy of the hemibiotrophic pathogen. As we known that effector proteins play key roles in interaction of hemibiotrophic fungi and their host plants. BAS4 (biotrophy-associated secreted protein 4) is an EIHM (extrainvasive hyphal membrane) matrix protein that was highly expressed in infectious hyphae. In order to study whether BAS4 is involved in the transition of rice blast fungus from biotrophic to necrotrophic phase, The susceptible rice cultivar Lijiangxintuanheigu (LTH) that were pre-treated with prokaryotic expression product of BAS4 and then followed with inoculation of the blast strain, more serious blast disease symptom, more biomass such as sporulation and fungal relative growth, and lower expression level of pathogenicity-related genes appeared in lesion of the rice leaves than those of the PBS-pretreated-leaves followed with inoculation of the same blast strain, which demonstrating that BAS4 invitro changed rice defense system to facilitate infection of rice blast strain. And the susceptible rice cultivar (LTH) were inoculated withBAS4-overexpressed blast strain, we also found more serious blast disease symptom and more biomass also appeared in lesion of leaves inoculated with BAS4-overexpressed strain than those of leaves inoculated with the wild-type strain, and expression level of pathogenicity-related genes appeared lower in biotrophic phase and higher in necrotrophic phase of infection, indicating BAS4 maybe in vivo regulate defense system of rice to facilitate transition of biotrophic to necrotrophic phase. Our data demonstrates that BAS4 in vitro and in vivo participates in transition from the biotrophic to the necrotrophic phase of Magnaporthe oryzae.

Contemporaneous radiations of fungi and plants linked to symbiosis

Interactions between fungi and plants, including parasitism, mutualism, and saprotrophy, have been invoked as key to their respective macroevolutionary success. Here we evaluate the origins of plant-fungal symbioses and saprotrophy using a time-calibrated phylogenetic framework that reveals linked and drastic shifts in diversification rates of each kingdom. Fungal colonization of land was associated with at least two origins of terrestrial green algae and preceded embryophytes (as evidenced by losses of fungal flagellum, ca. 720 Ma), likely facilitating terrestriality through endomycorrhizal and possibly endophytic symbioses. The largest radiation of fungi (Leotiomyceta), the origin of arbuscular mycorrhizae, and the diversification of extant embryophytes occurred ca. 480 Ma. This was followed by the origin of extant lichens. Saprotrophic mushrooms diversified in the Late Paleozoic as forests of seed plants started to dominate the landscape. The subsequent diversification and explosive radiation of Agaricomycetes, and eventually of ectomycorrhizal mushrooms, were associated with the evolution of Pinaceae in the Mesozoic, and establishment of angiosperm-dominated biomes in the Cretaceous.

Plants and fungi interact widely and in diverse ways, from mutualism to parasitism and decomposition. Here, Lutzoni et al. analyse the timing of plant and fungal evolutionary radiations and identify four major periods in which plant-fungal interactions likely drove lineage diversification.

Suppression of Plant Immunity by Fungal Chitinase-like Effectors

Crop diseases caused by fungi constitute one of the most important problems in agriculture, posing a serious threat to food security [1]. To establish infection, phytopathogens interfere with plant immune responses [2, 3]. However, strategies to promote virulence employed by fungal pathogens, especially non-model organisms, remain elusive [4], mainly because fungi are more complex and difficult to study when compared to the better-characterized bacterial pathogens. Equally incomplete is our understanding of the birth of microbial virulence effectors. Here, we show that the cacao pathogen Moniliophthora perniciosa evolved an enzymatically inactive chitinase (MpChi) that functions as a putative pathogenicity factor. MpChi is among the most highly expressed fungal genes during the biotrophic interaction with cacao and encodes a chitinase with mutations that abolish its enzymatic activity. Despite the lack of chitinolytic activity, MpChi retains substrate binding specificity and prevents chitin-triggered immunity by sequestering immunogenic chitin fragments. Remarkably, its sister species M. roreri encodes a second non-orthologous catalytically impaired chitinase with equivalent function. Thus, a class of conserved enzymes independently evolved as putative virulence factors in these fungi. In addition to unveiling a strategy of host immune suppression by fungal pathogens, our results demonstrate that the neofunctionalization of enzymes may be an evolutionary pathway for the rise of new virulence factors in fungi. We anticipate that analogous strategies are likely employed by other pathogens.

Genome sequence and effectorome of Moniliophthora perniciosa and Moniliophthora roreri subpopulations

Background

The hemibiotrophic pathogens Moniliophthora perniciosa (witches’ broom disease) and Moniliophthora roreri (frosty pod rot disease) are among the most important pathogens of cacao. Moniliophthora perniciosa has a broad host range and infects a variety of meristematic tissues in cacao plants, whereas M. roreri infects only pods of Theobroma and Herrania genera. Comparative pathogenomics of these fungi is essential to understand Moniliophthora infection strategies, therefore the detection and in silico functional characterization of effector candidates are important steps to gain insight on their pathogenicity.

Results

Candidate secreted effector proteins repertoire were predicted using the genomes of five representative isolates of M. perniciosa subpopulations (three from cacao and two from solanaceous hosts), and one representative isolate of M. roreri from Peru. Many putative effectors candidates were identified in M. perniciosa: 157 and 134 in cacao isolates from Bahia, Brazil; 109 in cacao isolate from Ecuador, 92 and 80 in wild solanaceous isolates from Minas Gerais (Lobeira) and Bahia (Caiçara), Brazil; respectively. Moniliophthora roreri showed the highest number of effector candidates, a total of 243. A set of eight core effectors were shared among all Moniliophthora isolates, while others were shared either between the wild solanaceous isolates or among cacao isolates. Mostly, candidate effectors of M. perniciosa were shared among the isolates, whereas in M. roreri nearly 50% were exclusive to the specie. In addition, a large number of cell wall-degrading enzymes characteristic of hemibiotrophic fungi were found. From these, we highlighted the proteins involved in cell wall modification, an enzymatic arsenal that allows the plant pathogens to inhabit environments with oxidative stress, which promotes degradation of plant compounds and facilitates infection.

Conclusions

The present work reports six genomes and provides a database of the putative effectorome of Moniliophthora, a first step towards the understanding of the functional basis of fungal pathogenicity.

Electronic supplementary material

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

Moniliophthora roreri, causal agent of cacao frosty pod rot

Taxonomy: Moniliophthora roreri (Cif.) H.C. Evans et al. ; Phylum Basidiomycota; Class Agaricomycetes; Order Agaricales; Family Marasmiaceae; Genus Moniliophthora. Biology: Moniliophthora roreri attacks Theobroma and Herrania species causing frosty pod rot. Theobroma cacao (cacao) is the host of major economic concern. Moniliophthora roreri is a hemibiotroph with a long biotrophic phase (45-90 days). Spore masses, of apparent asexual origin, are produced on the pod surface after initiation of the necrotrophic phase. Spores are spread by wind, rain and human activity. Symptoms of the biotrophic phase can include necrotic flecks and, in some cases, pod malformation, but pods otherwise remain asymptomatic. Relationship to Moniliophthora perniciosa: Moniliophthora roreri and Moniliophthora perniciosa, causal agent of witches' broom disease of cacao, are closely related. Their genomes are similar, including many of the genes they carry which are considered to be important in the disease process. Moniliophthora perniciosa, also a hemibiotroph, has a typical basidiomycete lifestyle and morphology, forming clamp connections and producing mushrooms. Basidiospores infect meristematic tissues including flower cushions, stem tips and pods. Moniliophthora roreri does not form clamp connections or mushrooms and infects pods only. Both pathogens are limited to the Western Hemisphere and are a threat to cacao production around the world. Agronomic importance: Disease losses caused by frosty pod rot can reach 90% and result in field abandonment. Moniliophthora roreri remains in the invasive phase in the Western Hemisphere, not having reached Brazil, some islands within the Caribbean and a few specific regions within otherwise invaded countries.

DISEASE MANAGEMENT

The disease can be managed by a combination of cultural (for example, maintenance of tree height and removal of infected pods) and chemical methods. These methods benefit from regional application, but can be cost prohibitive. Breeding for disease resistance offers the greatest potential for frosty pod rot management and new tolerant materials are becoming available.

Morphological variants of Moniliophthora roreri on artificial media and the biotroph/necrotroph shift

Moniliophthora roreri (Mr) causes frosty pod rot of Theobroma cacao in a hemibiotrophic association. The Mr biotroph-like phase has not been studied in culture. Mr spores (isolates Co12, Co52, and B3) were germinated on high (V8) and low (BPMM) nutrients with different media hardness (0.5% to 3% agarose). Germination was high on V8 media. Hardness affected germination on BPMM. Most colonies on V8 were slow-growing, failing to sporulate. Colony morphology depended on the isolate. On BPMM, exaggerated mycelia formed of limited length with enlarged cells. On agarose, rapidly expanding sporulating necrotrophic colonies formed rarely. Co12 and B3 spores were germinated on V8 and BPMM with low melting point (LMP) agarose. Slow-growing colonies of B3 on BPMM were unstable on LMP agarose, often forming slow-growing/rapidly expanding hybrids. Slow-growing colonies are hypothesized to represent the biotrophic phase. One nucleus was common in Mr cells, other than spores. Binucleate cells were occasionally observed in aged cells of slow-growing mycelia. Co52 cells often had more than two nuclei per cell after germination. Mr mycelia cells typically carry a single nucleus, being considered haploid. Biotroph- and necrotroph-like mycelia displayed differential gene expression but results were inconsistent with published in vivo results and require further study.

Comparative Proteomic Analysis of Different Isolates of Fusarium oxysporum f.sp. lycopersici to Exploit the Differentially Expressed Proteins Responsible for Virulence on Tomato Plants

The vascular wilt of tomato caused by Fusarium oxysporum f.sp. lycopersici is an important soil borne pathogen causes severe yield loss. The molecular characterization and their interaction with its host is necessary to develop a protection strategy. 20 isolates of F. oxysporum f.sp. lycopersici (FOL) were isolated from wilt infected tomato plants across Tamil Nadu. They were subjected to cultural, morphological, molecular and virulence studies. The results revealed that all the isolates produced both micro and macro conidia with different size, number of cells. The colors of the culture and growth pattern were also varied. In addition, chlamydospores were observed terminally and intercalary. The PCR analysis with F. oxysporum species-specific primer significantly amplified an amplicon of 600 bp fragment in all the isolates. Based on the above characters and pathogenicity, isolate FOL-8 was considered as virulent and FOL-20 was considered as least virulent. Proteomics strategy was adopted to determine the virulence factors between the isolates of FOL-8 and FOL-20. The 2D analyses have showed the differential expression of 17 different proteins. Among them, three proteins were down regulated and 14 proteins were significantly up regulated in FOL-8 than FOL-20 isolate. Among the 17 proteins, 10 distinct spots were analyzed by MALDI-TOF. The functions of the analyzed proteins, suggested that they were involved in pathogenicity, symptom expression and disease development, sporulation, growth, and higher penetration rate on tomato root tissue. Overall, these experiments proves the role of proteome in pathogenicity of F. oxysporum f.sp. lycopersici in tomato and unravels the mechanism behinds the virulence of the pathogen in causing wilt disease.

Prediction of Cacao (Theobroma cacao) Resistance to Moniliophthora spp. Diseases via Genome-Wide Association Analysis and Genomic Selection

Cacao (Theobroma cacao) is a globally important crop, and its yield is severely restricted by disease. Two of the most damaging diseases, witches' broom disease (WBD) and frosty pod rot disease (FPRD), are caused by a pair of related fungi: Moniliophthora perniciosa and Moniliophthora roreri, respectively. Resistant cultivars are the most effective long-term strategy to address Moniliophthora diseases, but efficiently generating resistant and productive new cultivars will require robust methods for screening germplasm before field testing. Marker-assisted selection (MAS) and genomic selection (GS) provide two potential avenues for predicting the performance of new genotypes, potentially increasing the selection gain per unit time. To test the effectiveness of these two approaches, we performed a genome-wide association study (GWAS) and GS on three related populations of cacao in Ecuador genotyped with a 15K single nucleotide polymorphism (SNP) microarray for three measures of WBD infection (vegetative broom, cushion broom, and chirimoya pod), one of FPRD (monilia pod) and two productivity traits (total fresh weight of pods and % healthy pods produced). GWAS yielded several SNPs associated with disease resistance in each population, but none were significantly correlated with the same trait in other populations. Genomic selection, using one population as a training set to estimate the phenotypes of the remaining two (composed of different families), varied among traits, from a mean prediction accuracy of 0.46 (vegetative broom) to 0.15 (monilia pod), and varied between training populations. Simulations demonstrated that selecting seedlings using GWAS markers alone generates no improvement over selecting at random, but that GS improves the selection process significantly. Our results suggest that the GWAS markers discovered here are not sufficiently predictive across diverse germplasm to be useful for MAS, but that using all markers in a GS framework holds substantial promise in accelerating disease-resistance in cacao.

Biotrophy-necrotrophy switch in pathogen evoke differential response in resistant and susceptible sesame involving multiple signaling pathways at different phases

Infection stages of charcoal rot fungus Macrophomina phaseolina in sesame revealed for the first time a transition from biotrophy via BNS (biotrophy-to-necrotrophy switch) to necrotrophy as confirmed by transcriptional studies. Microscopy using normal and GFP-expressing pathogen showed typical constricted thick intercellular bitrophic hyphae which gave rise to thin intracellular necrotrophic hyphae during BNS and this stage was delayed in a resistant host. Results also show that as the pathogen switched its strategy of infection, the host tailored its defense strategy to meet the changing situation. Less ROS accumulation, upregulation of ROS signaling genes and higher antioxidant enzyme activities post BNS resulted in resistance. There was greater accumulation of secondary metabolites and upregulation of secondary metabolite-related genes after BNS. A total of twenty genes functioning in different aspects of plant defense that were monitored over a time course during the changing infection phases showed a coordinated response. Experiments using phytohormone priming and phytohormone inhibitors showed that resistance resulted from activation of JA-ET signaling pathway. Most importantly this defense response was more prompt in the resistant than the susceptible host indicating that a resistant host makes different choices from a susceptible host during infection which ultimately influences the severity of the disease.

Tree spacing impacts the individual incidence of Moniliophthora roreri disease in cacao agroforests

BACKGROUND

Using conventional pesticides in crop protection has raised serious environmental concerns and there is therefore a need for integrated pest management (IPM) methods. In this paper, we found that the spacing of trees can impact disease, which could result in a reduction in pesticide applications and may act as a potential IPM method. We studied Frosty Pod Rot (FPR) in 20 cacao agroforests in Costa Rica (Upala region).

RESULTS

Using a generalized linear mixed model, we analyzed the impact of the neighborhood composition and distance from a studied cacao individual on its individual FPR incidence. We found that the number of cacao tree neighbors in a radius of 3.7 m and the number of fruit trees in a radius of 4.3 m had a significant negative influence on the incidence of FPR on individual cacao trees. Moreover, cacao tree neighbors had the most significant local influence compared to the neighborhood of other taller categories such as fruit or forest trees.

CONCLUSION

The mechanisms involved are related to the barrier effect, due to the effectiveness of the cacao tree's architecture as an efficient barrier against FPR spore dispersal. This paper provides new insights into optimization of the spatial environment around each host as an original IPM method. © 2017 Society of Chemical Industry.

Lifestyle, gene gain and loss, and transcriptional remodeling cause divergence in the transcriptomes of Phytophthora infestans and Pythium ultimum during potato tuber colonization

BACKGROUND

How pathogen genomes evolve to support distinct lifestyles is not well-understood. The oomycete Phytophthora infestans, the potato blight agent, is a largely biotrophic pathogen that feeds from living host cells, which become necrotic only late in infection. The related oomycete Pythium ultimum grows saprophytically in soil and as a necrotroph in plants, causing massive tissue destruction. To learn what distinguishes their lifestyles, we compared their gene contents and expression patterns in media and a shared host, potato tuber.

RESULTS

Genes related to pathogenesis varied in temporal expression pattern, mRNA level, and family size between the species. A family's aggregate expression during infection was not proportional to size due to transcriptional remodeling and pseudogenization. Ph. infestans had more stage-specific genes, while Py. ultimum tended towards more constitutive expression. Ph. infestans expressed more genes encoding secreted cell wall-degrading enzymes, but other categories such as secreted proteases and ABC transporters had higher transcript levels in Py. ultimum. Species-specific genes were identified including new Pythium genes, perforins, which may disrupt plant membranes. Genome-wide ortholog analyses identified substantial diversified expression, which correlated with sequence divergence. Pseudogenization was associated with gene family expansion, especially in gene clusters.

CONCLUSION

This first large-scale analysis of transcriptional divergence within oomycetes revealed major shifts in genome composition and expression, including subfunctionalization within gene families. Biotrophy and necrotrophy seem determined by species-specific genes and the varied expression of shared pathogenicity factors, which may be useful targets for crop protection.

Different waves of effector genes with contrasted genomic location are expressed by Leptosphaeria maculans during cotyledon and stem colonization of oilseed rape

Leptosphaeria maculans, the causal agent of stem canker disease, colonizes oilseed rape (Brassica napus) in two stages: a short and early colonization stage corresponding to cotyledon or leaf colonization, and a late colonization stage during which the fungus colonizes systemically and symptomlessly the plant during several months before stem canker appears. To date, the determinants of the late colonization stage are poorly understood; L. maculans may either successfully escape plant defences, leading to stem canker development, or the plant may develop an 'adult-stage' resistance reducing canker incidence. To obtain an insight into these determinants, we performed an RNA-sequencing (RNA-seq) pilot project comparing fungal gene expression in infected cotyledons and in symptomless or necrotic stems. Despite the low fraction of fungal material in infected stems, sufficient fungal transcripts were detected and a large number of fungal genes were expressed, thus validating the feasibility of the approach. Our analysis showed that all avirulence genes previously identified are under-expressed during stem colonization compared with cotyledon colonization. A validation RNA-seq experiment was then performed to investigate the expression of candidate effector genes during systemic colonization. Three hundred and seven 'late' effector candidates, under-expressed in the early colonization stage and over-expressed in the infected stems, were identified. Finally, our analysis revealed a link between the regulation of expression of effectors and their genomic location: the 'late' effector candidates, putatively involved in systemic colonization, are located in gene-rich genomic regions, whereas the 'early' effector genes, over-expressed in the early colonization stage, are located in gene-poor regions of the genome.

Computational identification of Y-linked markers and genes in the grass carp genome by using a pool-and-sequence method

The molecular analysis of sex in vertebrates is important, as it has the potential to provide vital information for theoretical and applied research alike. Teleost fish are the ancient vertebrates that present a broad sex chromosome system but lack differentiated sex chromosomes in most species. Hence understanding the sex in fish would not only illuminate the sex determination evolution in vertebrates but also shed light on fish farming. In the present study, we used grass carp as a teleost fish model, studied the Y chromosome by using a pool-and-sequence strategy in combination with fragment-ratio method. In total, we identified five Y-linked scaffolds (totaling 347 Kb) and six Y-specific sequences that could be used as sex-specific markers, demonstrating the suitability of NGS-based re-sequencing of pooled DNAs for the identification of sex markers in fish. Moreover, 14 putative Y-linked genes were described for the first time. All the genes, except for un-y1, un-y2, and ubq-y, showed high similarity to their female homologs. RT-PCR revealed that ubq-y was only expressed in the male hypothalamus and pituitary. These findings provided an abundant resource for the Y chromosome of grass carp, and may help elucidate sex chromosome evolution in cyprinid fish.

Saprotrophic proteomes of biotypes of the witches' broom pathogen Moniliophthora perniciosa

Nine geographically diverse Moniliophthora perniciosa (witches' broom disease pathogen) isolates were cultured in vitro. They included six C-biotypes differing in virulence on cacao (Theobroma cacao), two S-biotypes (solanaceous hosts), and an L-biotype (liana hosts). Mycelial growth rates and morphologies differed considerably, but no characters were observed to correlate with virulence or biotype. In plant inoculations using basidiospores, one C-biotype caused symptoms on tomato (an S-biotype host), adding to evidence of limited host adaptation in these biotypes. Mycelial proteomes were analysed by two-dimensional gel electrophoresis (2-DE), and 619 gel spots were indexed on all replicate gels of at least one strain. Multivariate analysis of gel spots discriminated the L-biotype, but not the S-biotypes, from the remaining strains. The proteomic similarity of the S- and C-biotypes is consistent with their reported lack of phylogenetic distinction. Sequences from tandem mass spectrometry of tryptic peptides from major 2-DE spots were matched with Moniliophthora genome and transcript sequences on NCBI and WBD Transcriptome Atlas databases. Protein-spot identifications indicated that M. perniciosa saprotrophic mycelial proteomes expressed functions potentially connected with a 'virulence life-style', including peroxiredoxin, heat-shock proteins, nitrilase, formate dehydrogenase, a prominent complement of aldo-keto reductases, mannitol-1-phosphate dehydrogenase, and central metabolism enzymes with proposed pathogenesis functions.

In vitro secretomic analysis identifies putative pathogenicity-related proteins of Sporisorium scitamineum - The sugarcane smut fungus

Sporisorium scitamineum, the sugarcane smut pathogen, relies predominantly on its secretome to successfully colonise its host, in accordance with other related smut fungi. Considering the significance of deciphering its secretome, we have examined alterations in the in vitro secretome of S. scitamineum in response to synthetic and sugarcane meristem tissue-amended growth media, so as to identify host signal responsive secretory proteins. Secretory proteins that were differentially abundant and exclusively secreted in response to host extract media were identified by two-dimensional gel electrophoresis coupled with MALDI-TOF/TOF MS. Of the 16 differentially abundant and exclusively secreted proteins, nine proteins were identified. Among which, six were related to cell wall modification, morphogenesis, polysaccharide degradation, and carbohydrate metabolism. In planta gene expression profiling indicated that five in vitro secreted proteins were expressed in distinct patterns by S. scitamineum during different stages of infection with relatively higher expression at 1 day after inoculation, suggesting that these proteins could be aiding S. scitamineum at early time points in penetration and colonisation of sugarcane cells. The present study has provided insights into the alterations occurring in the secretome of S. scitamineum at in vitro conditions and has resulted in the identification of secretory proteins that are possibly associated with pathogenicity of the sugarcane smut fungus.

Global transcriptional analysis suggests Lasiodiplodia theobromae pathogenicity factors involved in modulation of grapevine defensive response

BACKGROUND

Lasiodiplodia theobromae is a fungus of the Botryosphaeriaceae that causes grapevine vascular disease, especially in regions with hot climates. Fungi in this group often remain latent within their host and become virulent under abiotic stress. Transcriptional regulation analysis of L. theobromae exposed to heat stress (HS) was first carried out in vitro in the presence of grapevine wood (GW) to identify potential pathogenicity genes that were later evaluated for in planta expression.

RESULTS

A total of 19,860 de novo assembled transcripts were obtained, forty-nine per cent of which showed homology to the Botryosphaeriaceae fungi, Neofusicoccum parvum or Macrophomina phaseolina. Three hundred ninety-nine have homology with genes involved in pathogenic processes and several belonged to expanded gene families in others fungal grapevine vascular pathogens. Gene expression analysis showed changes in fungal metabolism of phenolic compounds; where genes encoding for enzymes, with the ability to degrade salicylic acid (SA) and plant phenylpropanoid precursors, were up-regulated during in vitro HS response, in the presence of GW. These results suggest that the fungal L-tyrosine catabolism pathway could help the fungus to remove phenylpropanoid precursors thereby evading the host defense response. The in planta up-regulation of salicylate hydroxylase, intradiol ring cleavage dioxygenase and fumarylacetoacetase encoding genes, further supported this hypothesis. Those genes were even more up-regulated in HS-stressed plants, suggesting that fungus takes advantage of the increased phenylpropanoid precursors produced under stress. Pectate lyase was up-regulated while a putative amylase was down-regulated in planta, this could be associated with an intercellular growth strategy during the first stages of colonization.

CONCLUSIONS

L. theobromae transcriptome was established and validated. Its usefulness was demonstrated through the identification of genes expressed during the infection process. Our results support the hypothesis that heat stress facilitates fungal colonization, because of the fungus ability to use the phenylpropanoid precursors and SA, both compounds known to control host defense.

Geographic Differentiation and Population Genetic Structure of Moniliophthora roreri in the Principal Cocoa Production Areas in Colombia

Frosty pod rot (FPR) disease on cocoa, caused by Moniliophthora roreri, is one of the most devastating cocoa disease in the Western Hemisphere. In Colombia, the disease is particularly severe in the Magdalena Valley, which is considered the possible center of origin for the pathogen species. We analyzed the genetic diversity of isolates from the departments of Santander, Antioquia, Tolima, and Huila in Colombia using 23 simple-sequence repeats (SSR) markers. In total, 117 different multilocus genotypes were found among 120 isolates, each one representing a unique haplotype. High mutation rates in the SSR and gene flow can explain the high levels of diversity. Also, the observed and standardized indexes of association (I_A and řd) indicate that the populations of M. roreri are clonal. Furthermore, given the high haplotype diversity and the significant linkage disequilibrium observed, we hypothesize that M. roreri could be a primarily asexual species undergoing sporadic recombination or partial recombination through parasexuality. A Bayesian clustering analysis implemented by STRUCTURE showed that the most probable number of genetic groups in the data was three, confirming the geographical differentiation among isolates. Similar results were obtained by a discriminant analysis of principal components, a principal coordinate analysis, and a neighbor-joining tree from microsatellite loci base on Nei distance. Cacao genotypes and environmental variables did contribute to the genetic differentiation of the groups. We discuss how this information could be used to improve the management of FPR at the regional level.

Multi-genome analysis identifies functional and phylogenetic diversity of basidiomycete adenylate-forming reductases

Among the invaluable benefits of basidiomycete genomics is the dramatically enhanced insight into the potential capacity to biosynthesize natural products. This study focuses on adenylate-forming reductases, which is a group of natural product biosynthesis enzymes that resembles non-ribosomal peptide synthetases, yet serves to modify one substrate, rather than to condense two or more building blocks. Phylogenetically, these reductases fall in four classes. The phylogeny of Heterobasidion annosum (Russulales) and Serpula lacrymans (Boletales) adenylate-forming reductases was investigated. We identified a previously unrecognized phylogenetic branch within class III adenylate-forming reductases. Three representatives were heterologously produced and their substrate preferences determined in vitro: NPS9 and NPS11 of S. lacrymans preferred l-threonine and benzoic acid, respectively, while NPS10 of H. annosum accepted phenylpyruvic acid best. We also investigated two class IV adenylate-forming reductases of Coprinopsis cinerea, which each were active with l-alanine, l-valine, and l-serine as substrates. Our results show that adenylate-forming reductases are functionally more diverse than previously recognized. As none of the natural products known from the species investigated in this study includes the identified substrates of their respective reductases, our findings may help further explore the diversity of these basidiomycete secondary metabolomes.

OcculterCut: A Comprehensive Survey of AT-Rich Regions in Fungal Genomes

We present a novel method to measure the local GC-content bias in genomes and a survey of published fungal species. The method, enacted as "OcculterCut" (https://sourceforge.net/projects/occultercut, last accessed April 30, 2016), identified species containing distinct AT-rich regions. In most fungal taxa, AT-rich regions are a signature of repeat-induced point mutation (RIP), which targets repetitive DNA and decreases GC-content though the conversion of cytosine to thymine bases. RIP has in turn been identified as a driver of fungal genome evolution, as RIP mutations can also occur in single-copy genes neighboring repeat-rich regions. Over time RIP perpetuates "two speeds" of gene evolution in the GC-equilibrated and AT-rich regions of fungal genomes. In this study, genomes showing evidence of this process are found to be common, particularly among the Pezizomycotina. Further analysis highlighted differences in amino acid composition and putative functions of genes from these regions, supporting the hypothesis that these regions play an important role in fungal evolution. OcculterCut can also be used to identify genes undergoing RIP-assisted diversifying selection, such as small, secreted effector proteins that mediate host-microbe disease interactions.

The cacao pathogen Moniliophthora roreri (Marasmiaceae) possesses biallelic A and B mating loci but reproduces clonally

The cacao pathogen Moniliophthora roreri belongs to the mushroom-forming family Marasmiaceae, but it has never been observed to produce a fruiting body, which calls to question its capacity for sexual reproduction. In this study, we identified potential A (HD1 and HD2) and B (pheromone precursors and pheromone receptors) mating genes in M. roreri. A PCR-based method was subsequently devised to determine the mating type for a set of 47 isolates from across the geographic range of the fungus. We developed and generated an 11-marker microsatellite set and conducted association and linkage disequilibrium (standardized index of association, IA(s)) analyses. We also performed an ancestral reconstruction analysis to show that the ancestor of M. roreri is predicted to be heterothallic and tetrapolar, which together with sliding window analyses support that the A and B mating loci are likely unlinked and follow a tetrapolar organization within the genome. The A locus is composed of a pair of HD1 and HD2 genes, whereas the B locus consists of a paired pheromone precursor, Mr_Ph4, and receptor, STE3_Mr4. Two A and B alleles but only two mating types were identified. Association analyses divided isolates into two well-defined genetically distinct groups that correlate with their mating type; IA(s) values show high linkage disequilibrium as is expected in clonal reproduction. Interestingly, both mating types were found in South American isolates but only one mating type was found in Central American isolates, supporting a prior hypothesis of clonal dissemination throughout Central America after a single or very few introductions of the fungus from South America.