Transcriptomic analysis identifies candidate genes for Aphanomyces root rot disease resistance in pea

BACKGROUND

Aphanomyces euteiches is a soil-borne oomycete that causes root rot in pea and other legume species. Symptoms of Aphanomyces root rot (ARR) include root discoloration and wilting, leading to significant yield losses in pea production. Resistance to ARR is known to be polygenic but the roles of single genes in the pea immune response are still poorly understood. This study uses transcriptomics to elucidate the immune response of two pea genotypes varying in their levels of resistance to A. euteiches.

RESULTS

In this study, we inoculated roots of the pea (P. sativum L.) genotypes 'Linnea' (susceptible) and 'PI180693' (resistant) with two different A. euteiches strains varying in levels of virulence. The roots were harvested at 6 h post-inoculation (hpi), 20 hpi and 48 hpi, followed by differential gene expression analysis. Our results showed a time- and genotype-dependent immune response towards A. euteiches infection, involving several WRKY and MYB-like transcription factors, along with genes associated with jasmonic acid (JA) and abscisic acid (ABA) signaling. By cross-referencing with genes segregating with partial resistance to ARR, we identified 39 candidate disease resistance genes at the later stage of infection. Among the genes solely upregulated in the resistant genotype 'PI180693', Psat7g091800.1 was polymorphic between the pea genotypes and encoded a Leucine-rich repeat receptor-like kinase reminiscent of the Arabidopsis thaliana FLAGELLIN-SENSITIVE 2 receptor.

CONCLUSIONS

This study provides new insights into the gene expression dynamics controlling the immune response of resistant and susceptible pea genotypes to A. euteiches infection. We present a set of 39 candidate disease resistance genes for ARR in pea, including the putative immune receptor Psat7g091800.1, for future functional validation.

Verticillium longisporum phospholipase VlsPLA2 is a virulence factor that targets host nuclei and modulates plant immunity

Phospholipase A2 (PLA2 ) is a lipolytic enzyme that hydrolyses phospholipids in the cell membrane. In the present study, we investigated the role of secreted PLA2 (VlsPLA2 ) in Verticillium longisporum, a fungal phytopathogen that mostly infects plants belonging to the Brassicaceae family, causing severe annual yield loss worldwide. Expression of the VlsPLA2 gene, which encodes active PLA2 , is highly induced during the interaction of the fungus with the host plant Brassica napus. Heterologous expression of VlsPLA2 in Nicotiana benthamiana resulted in increased synthesis of certain phospholipids compared to plants in which enzymatically inactive PLA2 was expressed (VlsPLA2 ΔCD ). Moreover, VlsPLA2 suppresses the hypersensitive response triggered by the Cf4/Avr4 complex, thereby suppressing the chitin-induced reactive oxygen species burst. VlsPLA2 -overexpressing V. longisporum strains showed increased virulence in Arabidopsis plants, and transcriptomic analysis of this fungal strain revealed that the induction of the gene contributed to increased virulence. VlsPLA2 was initially localized to the host nucleus and then translocated to the chloroplasts at later time points. In addition, VlsPLA2 bound to the vesicle-associated membrane protein A (VAMPA) and was transported to the nuclear membrane. In the nucleus, VlsPLA2 caused major alterations in the expression levels of genes encoding transcription factors and subtilisin-like proteases, which play a role in plant immunity. In conclusion, our study showed that VlsPLA2 acts as a virulence factor, possibly by hydrolysing host nuclear envelope phospholipids, which, through a signal transduction cascade, may suppress basal plant immune responses.

Functional analysis of the apple fruit microbiome based on shotgun metagenomic sequencing of conventional and organic orchard samples

Fruits harbour abundant and diverse microbial communities that protect them from post-harvest pathogens. Identification of functional traits associated with a given microbiota can provide a better understanding of their potential influence. Here, we focused on the epiphytic microbiome of apple fruit. We suggest that shotgun metagenomic data can indicate specific functions carried out by different groups and provide information on their potential impact. Samples were collected from the surface of 'Golden Delicious' apples from four orchards that differ in their geographic location and management practice. Approximately 1 million metagenes were predicted based on a high-quality assembly. Functional profiling of the microbiome of fruits from orchards differing in their management practice revealed a functional shift in the microbiota. The organic orchard microbiome was enriched in pathways involved in plant defence activities; the conventional orchard microbiome was enriched in pathways related to the synthesis of antibiotics. The functional significance of the variations was explored using microbial network modelling algorithms to reveal the metabolic role of specific phylogenetic groups. The analysis identified several associations supported by other published studies. For example, the analysis revealed the nutritional dependencies of the Capnodiales group, including the Alternaria pathogen, on aromatic compounds.

A framework for the targeted recruitment of crop-beneficial soil taxa based on network analysis of metagenomics data

BACKGROUND

The design of ecologically sustainable and plant-beneficial soil systems is a key goal in actively manipulating root-associated microbiomes. Community engineering efforts commonly seek to harness the potential of the indigenous microbiome through substrate-mediated recruitment of beneficial members. In most sustainable practices, microbial recruitment mechanisms rely on the application of complex organic mixtures where the resources/metabolites that act as direct stimulants of beneficial groups are not characterized. Outcomes of such indirect amendments are unpredictable regarding engineering the microbiome and achieving a plant-beneficial environment.

RESULTS

This study applied network analysis of metagenomics data to explore amendment-derived transformations in the soil microbiome, which lead to the suppression of pathogens affecting apple root systems. Shotgun metagenomic analysis was conducted with data from 'sick' vs 'healthy/recovered' rhizosphere soil microbiomes. The data was then converted into community-level metabolic networks. Simulations examined the functional contribution of treatment-associated taxonomic groups and linked them with specific amendment-induced metabolites. This analysis enabled the selection of specific metabolites that were predicted to amplify or diminish the abundance of targeted microbes functional in the healthy soil system. Many of these predictions were corroborated by experimental evidence from the literature. The potential of two of these metabolites (dopamine and vitamin B₁₂) to either stimulate or suppress targeted microbial groups was evaluated in a follow-up set of soil microcosm experiments. The results corroborated the stimulant's potential (but not the suppressor) to act as a modulator of plant beneficial bacteria, paving the way for future development of knowledge-based (rather than trial and error) metabolic-defined amendments. Our pipeline for generating predictions for the selective targeting of microbial groups based on processing assembled and annotated metagenomics data is available at https://github.com/ot483/NetCom2 .

CONCLUSIONS

This research demonstrates how genomic-based algorithms can be used to formulate testable hypotheses for strategically engineering the rhizosphere microbiome by identifying specific compounds, which may act as selective modulators of microbial communities. Applying this framework to reduce unpredictable elements in amendment-based solutions promotes the development of ecologically-sound methods for re-establishing a functional microbiome in agro and other ecosystems. Video Abstract.

RNA silencing proteins and small RNAs in oomycete plant pathogens and biocontrol agents

Introduction

Oomycetes cause several damaging diseases of plants and animals, and some species also act as biocontrol agents on insects, fungi, and other oomycetes. RNA silencing is increasingly being shown to play a role in the pathogenicity of Phytophthora species, either through trans-boundary movement of small RNAs (sRNAs) or through expression regulation of infection promoting effectors.

Methods

To gain a wider understanding of RNA silencing in oomycete species with more diverse hosts, we mined genome assemblies for Dicer-like (DCL), Argonaute (AGO), and RNA dependent RNA polymerase (RDRP) proteins from Phytophthora plurivora, Ph. cactorum, Ph. colocasiae, Pythium oligandrum, Py. periplocum, and Lagenidium giganteum. Moreover, we sequenced small RNAs from the mycelium stage in each of these species.

Results and discussion

Each of the species possessed a single DCL protein, but they differed in the number and sequence of AGOs and RDRPs. SRNAs of 21nt, 25nt, and 26nt were prevalent in all oomycetes analyzed, but the relative abundance and 5' base preference of these classes differed markedly between genera. Most sRNAs mapped to transposons and other repeats, signifying that the major role for RNA silencing in oomycetes is to limit the expansion of these elements. We also found that sRNAs may act to regulate the expression of duplicated genes. Other sRNAs mapped to several gene families, and this number was higher in Pythium spp., suggesting a role of RNA silencing in regulating gene expression. Genes for most effector classes were the source of sRNAs of variable size, but some gene families showed a preference for specific classes of sRNAs, such as 25/26 nt sRNAs targeting RxLR effector genes in Phytophthora species. Novel miRNA-like RNAs (milRNAs) were discovered in all species, and two were predicted to target transcripts for RxLR effectors in Ph. plurivora and Ph. cactorum, indicating a putative role in regulating infection. Moreover, milRNAs from the biocontrol Pythium species had matches in the predicted transcriptome of Phytophthora infestans and Botrytis cinerea, and L. giganteum milRNAs matched candidate genes in the mosquito Aedes aegypti. This suggests that trans-boundary RNA silencing may have a role in the biocontrol action of these oomycetes.

Insights into the ecological generalist lifestyle of Clonostachys fungi through analysis of their predicted secretomes

Introduction

The fungal secretome comprise diverse proteins that are involved in various aspects of fungal lifestyles, including adaptation to ecological niches and environmental interactions. The aim of this study was to investigate the composition and activity of fungal secretomes in mycoparasitic and beneficial fungal-plant interactions.

Methods

We used six Clonostachys spp. that exhibit saprotrophic, mycotrophic and plant endophytic lifestyles. Genome-wide analyses was performed to investigate the composition, diversity, evolution and gene expression of Clonostachys secretomes in relation to their potential role in mycoparasitic and endophytic lifestyles.

Results and discussion

Our analyses showed that the predicted secretomes of the analyzed species comprised between 7 and 8% of the respective proteomes. Mining of transcriptome data collected during previous studies showed that 18% of the genes encoding predicted secreted proteins were upregulated during the interactions with the mycohosts Fusarium graminearum and Helminthosporium solani. Functional annotation of the predicted secretomes revealed that the most represented protease family was subclass S8A (11-14% of the total), which include members that are shown to be involved in the response to nematodes and mycohosts. Conversely, the most numerous lipases and carbohydrate-active enzyme (CAZyme) groups appeared to be potentially involved in eliciting defense responses in the plants. For example, analysis of gene family evolution identified nine CAZyme orthogroups evolving for gene gains (p ≤ 0.05), predicted to be involved in hemicellulose degradation, potentially producing plant defense-inducing oligomers. Moreover, 8-10% of the secretomes was composed of cysteine-enriched proteins, including hydrophobins, important for root colonization. Effectors were more numerous, comprising 35-37% of the secretomes, where certain members belonged to seven orthogroups evolving for gene gains and were induced during the C. rosea response to F. graminearum or H. solani. Furthermore, the considered Clonostachys spp. possessed high numbers of proteins containing Common in Fungal Extracellular Membranes (CFEM) modules, known for their role in fungal virulence. Overall, this study improves our understanding of Clonostachys spp. adaptation to diverse ecological niches and establishes a basis for future investigation aiming at sustainable biocontrol of plant diseases.

Computational Analysis of HTS Data and Its Application in Plant Pathology

High-throughput sequencing is a basic tool of biological research, and it is extensively used in plant pathology projects. Here, we describe how to handle data coming from a variety of sequencing experiments, focusing on the analysis of Illumina reads. We describe how to perform genome assembly and annotation with DNA reads, correctly analyze RNA-seq data to discover differentially expressed genes, handle amplicon sequencing data from microbial communities, and utilize small RNA sequencing data to predict miRNA sequences and their putative targets.

Role of Dicer-Dependent RNA Interference in Regulating Mycoparasitic Interactions

Dicer-like proteins (DCLs) play a vital role in RNA interference (RNAi), by cleaving RNA filament into small RNAs. Although DCL-mediated RNAi can regulate interspecific communication between pathogenic/mutualistic organisms and their hosts, its role in mycoparasitic interactions is yet to be investigated. In this study, we deleted dcl genes in the mycoparasitic fungus Clonostachys rosea and characterize the functions of DCL-dependent RNAi in mycoparasitism. Deletion of dcl2 resulted in a mutant with reduced secondary metabolite production, antagonism toward the plant-pathogenic fungus Botrytis cinerea, and reduced ability to control Fusarium foot rot disease on wheat, caused by Fusarium graminearum. Transcriptome sequencing of the in vitro interaction between the C. rosea Δdcl2 strain and B. cinerea or F. graminearum identified the downregulation of genes coding for transcription factors, membrane transporters, hydrolytic enzymes, and secondary metabolites biosynthesis enzymes putatively involved in antagonistic interactions, in comparison with the C. rosea wild-type interaction. A total of 61 putative novel microRNA-like RNAs (milRNAs) were identified in C. rosea, and 11 were downregulated in the Δdcl2 mutant. In addition to putative endogenous gene targets, these milRNAs were predicted to target B. cinerea and F. graminearum virulence factor genes, which showed an increased expression during interaction with the Δdcl2 mutant incapable of producing the targeting milRNAs. In summary, this study constitutes the first step in elucidating the role of RNAi in mycoparasitic interactions, with important implications for biological control of plant diseases, and poses the base for future studies focusing on the role of cross-species RNAi regulating mycoparasitic interactions. IMPORTANCE Small RNAs mediated RNA interference (RNAi) known to regulate several biological processes. Dicer-like endoribonucleases (DCLs) play a vital role in the RNAi pathway by generating sRNAs. In this study, we investigated a role of DCL-mediated RNAi in interference interactions between mycoparasitic fungus Clonostachys rosea and the two fungal pathogens Botrytis cinerea and Fusarium graminearum (here called mycohosts). We found that the dcl mutants were not able to produce 11 sRNAs predicted to finetune the regulatory network of genes known to be involved in production of hydrolytic enzymes, antifungal compounds, and membrane transporters needed for antagonistic action of C. rosea. We also found C. rosea sRNAs putatively targeting known virulence factors in the mycohosts, indicating RNAi-mediated cross-species communication. Our study expanded the understanding of underlying mechanisms of cross-species communication during interference interactions and poses a base for future works studying the role of DCL-based cross-species RNAi in fungal interactions.

CRISPR-Cas9-Based Discovery of the Verrucosidin Biosynthesis Gene Cluster in Penicillium polonicum

Penicillium polonicum, commonly found on food matrices, is a mycotoxigenic species able to produce a neurotoxin called verrucosidin. This methylated α-pyrone polyketide inhibits oxidative phosphorylation in mitochondria and thereby causes neurological diseases. Despite the importance of verrucosidin as a toxin, its biosynthetic genes have not been characterized yet. By similarity analysis with the polyketide synthase (PKS) genes for the α-pyrones aurovertin (AurA) and citreoviridin (CtvA), 16 PKS genes for putative α-pyrones were identified in the P. polonicum genome. A single PKS gene, verA, was found to be transcribed under verrucosidin-producing growth conditions. The annotated functions of the genes neighboring verA correspond to those required for verrucosidin biosynthesis. To prove the involvement of verA in verrucosidin biosynthesis, the clustered regularly interspaced short palindrome repeats (CRISPR) technology was applied to P. polonicum. In vitro reconstituted CRISPR-Cas9 was used to induce targeted gene deletions in P. polonicum. This approach allowed identifying and characterizing the verrucosidin biosynthetic gene cluster. VerA deletion mutants were no longer able to produce verrucosidin, whereas they were displaying morphological characteristics comparable with the wild-type strain. The available CRISPR-Cas9 technology allows characterizing the biosynthetic potential of P. polonicum as a valuable source of novel compounds.

Metagenomics Approaches for the Detection and Surveillance of Emerging and Recurrent Plant Pathogens

Globalization has a dramatic effect on the trade and movement of seeds, fruits and vegetables, with a corresponding increase in economic losses caused by the introduction of transboundary plant pathogens. Current diagnostic techniques provide a useful and precise tool to enact surveillance protocols regarding specific organisms, but this approach is strictly targeted, while metabarcoding and shotgun metagenomics could be used to simultaneously detect all known pathogens and potentially new ones. This review aims to present the current status of high-throughput sequencing (HTS) diagnostics of fungal and bacterial plant pathogens, discuss the challenges that need to be addressed, and provide direction for the development of methods for the detection of a restricted number of related taxa (specific surveillance) or all of the microorganisms present in a sample (general surveillance). HTS techniques, particularly metabarcoding, could be useful for the surveillance of soilborne, seedborne and airborne pathogens, as well as for identifying new pathogens and determining the origin of outbreaks. Metabarcoding and shotgun metagenomics still suffer from low precision, but this issue can be limited by carefully choosing primers and bioinformatic algorithms. Advances in bioinformatics will greatly accelerate the use of metagenomics to address critical aspects related to the detection and surveillance of plant pathogens in plant material and foodstuffs.

Elaborated regulation of griseofulvin biosynthesis in Penicillium griseofulvum and its role on conidiation and virulence

Penicilium griseofulvum, the causal agent of apple blue mold, is able to produce in vitro and on apple a broad spectrum of secondary metabolites (SM), including patulin, roquefortine C and griseofulvin. Among them, griseofulvin is known for its antifungal and antiproliferative activity, and has received interest in many sectors, from medicine to agriculture. The biosynthesis of SM is finely regulated by filamentous fungi and can involve global regulators and pathway specific regulators, which are usually encoded by genes present in the same gene cluster as the backbone gene and tailoring enzymes. In the griseofulvin gene cluster, two putative transcription factors were previously identified, encoded by genes gsfR1 and gsfR2, and their role has been investigated in the present work. Analysis of P. griseofulvum knockout mutants lacking either gene suggest that gsfR2 forms part of a different pathway and gsfR1 exhibits many spectra of action, acting as regulator of griseofulvin and patulin biosynthesis and influencing conidia production and virulence on apple. The analysis of gsfR1 promoter revealed that the regulation of griseofulvin biosynthesis is also controlled by global regulators in response to many environmental stimuli, such as carbon and nitrogen. The influence of carbon and nitrogen on griseofulvin production was further investigated and verified, revealing a complex network of response and confirming the central role of gsfR1 in many processes in P. griseofulvum.

Different Phenotypes, Similar Genomes: Three Newly Sequenced Fusarium fujikuroi Strains Induce Different Symptoms in Rice Depending on Temperature

Bakanae, caused by the hemibiotrophic fungus Fusarium fujikuroi, is one of the most important diseases of rice and is attributed to up to 75% of losses, depending on the strain and environmental conditions. Some strains cause elongation and thin leaves, whereas others induce stunting and chlorotic seedlings. Differences in symptoms are attributed to genetic differences in the strains. F. fujikuroi strains Augusto2, CSV1, and I1.3 were sequenced with Illumina MiSeq, and pathogenicity trials were conducted on rice cultivar Galileo, which is susceptible to bakanae. By performing gene prediction, single nucleotide polymorphism (SNP) calling, and structural variant analysis with a reference genome, we show how an extremely limited number of polymorphisms in genes not commonly associated with bakanae disease can cause strong differences in phenotype. CSV1 and Augusto2 were particularly close, with only 21,887 SNPs between them, but they differed in virulence, reaction to temperature, induced symptoms, colony morphology and color, growth speed, fumonisin, and gibberellin production. Genes potentially involved in the shift in phenotype were identified. Furthermore, we show how temperature variation may result in different symptoms even in rice plants inoculated with the same F. fujikuroi strain. Moreover, all of the F. fujikuroi strains became more virulent at higher temperatures. Significant differences were likewise observed in gibberellic acid production and in the expression of both fungal and plant gibberellin biosynthetic genes.

Genome Sequence, Assembly, and Characterization of the Antagonistic Yeast Candida oleophila Used as a Biocontrol Agent Against Post-harvest Diseases

Candida oleophila is an effective biocontrol agent used to control post-harvest diseases of fruits and vegetables. C. oleophila I-182 was the active agent used in the first-generation yeast-based commercial product, Aspire^®, for post-harvest disease management. Several action modes, like competition for nutrients and space, induction of pathogenesis-related genes in host tissues, and production of extracellular lytic enzymes, have been demonstrated for the biological control activity exhibited by C. oleophila through which it inhibits post-harvest pathogens. In the present study, the whole genome of C. oleophila I-182 was sequenced using PacBio and Illumina shotgun sequencing technologies, yielding an estimated genome size of 14.73 Mb. The genome size is similar in length to that of the model yeast strain Saccharomyces cerevisiae S288c. Based on the assembled genome, protein-coding sequences were identified and annotated. The predicted genes were further assigned with gene ontology terms and clustered in special functional groups. A comparative analysis of C. oleophila proteome with the proteomes of 11 representative yeasts revealed 2 unique and 124 expanded families of proteins in C. oleophila. Availability of the genome sequence will facilitate a better understanding the properties of biocontrol yeasts at the molecular level.

Genome Sequence, Assembly and Characterization of Two Metschnikowia fructicola Strains Used as Biocontrol Agents of Postharvest Diseases

The yeast Metschnikowia fructicola was reported as an efficient biological control agent of postharvest diseases of fruits and vegetables, and it is the bases of the commercial formulated product "Shemer." Several mechanisms of action by which M. fructicola inhibits postharvest pathogens were suggested including iron-binding compounds, induction of defense signaling genes, production of fungal cell wall degrading enzymes and relatively high amounts of superoxide anions. We assembled the whole genome sequence of two strains of M. fructicola using PacBio and Illumina shotgun sequencing technologies. Using the PacBio, a high-quality draft genome consisting of 93 contigs, with an estimated genome size of approximately 26 Mb, was obtained. Comparative analysis of M. fructicola proteins with the other three available closely related genomes revealed a shared core of homologous proteins coded by 5,776 genes. Comparing the genomes of the two M. fructicola strains using a SNP calling approach resulted in the identification of 564,302 homologous SNPs with 2,004 predicted high impact mutations. The size of the genome is exceptionally high when compared with those of available closely related organisms, and the high rate of homology among M. fructicola genes points toward a recent whole-genome duplication event as the cause of this large genome. Based on the assembled genome, sequences were annotated with a gene description and gene ontology (GO term) and clustered in functional groups. Analysis of CAZymes family genes revealed 1,145 putative genes, and transcriptomic analysis of CAZyme expression levels in M. fructicola during its interaction with either grapefruit peel tissue or Penicillium digitatum revealed a high level of CAZyme gene expression when the yeast was placed in wounded fruit tissue.