Genome-Wide Transcriptome Analysis Reveals the Comprehensive Response of Two Susceptible Poplar Sections to Marssonina brunnea Infection

Transcriptome Analysis and Reactive Oxygen Species Detection Suggest Contrasting Molecular Mechanisms in Populus canadensis' Response to Different Formae Speciales of Marssonina brunnea

Revealing plant-pathogen interactions is important for resistance breeding, but it remains a complex process that presents many challenges. Marssonina leaf spot of poplars (MLSP) is the main disease in poplars; in China, its pathogens consist of two formae speciales, namely, Marssonina brunnea f. sp. Monogermtubi (MO) and M. brunnea f. sp. Multigermtubi (MU). However, the mechanism of the molecular interaction between poplars and the two formae speciales, especially for an incompatible system, remains unclear. In this study, we conducted transcriptome sequencing and reactive oxygen species (ROS) staining based on the interactions between Populus canadensis and the two formae speciales. The results show that the gene expression patterns of P. canadensis induced by MO and MU were significantly different, especially for the genes associated with biotic stress. Furthermore, MO and MU also triggered distinct ROS reactions of P. canadensis, and ROS (mainly H2O2) burst was only observed around the cells penetrated by MU. In conclusion, this study suggested that P. canadensis experienced different resistance reactions in response to the two formae speciales of M. brunnea, providing valuable insights for further understanding the host-pathogen interactions of MLSP.

The responses of poplars to fungal pathogens: A review of the defensive pathway

Long-lived tree species need to cope with changing environments and pathogens during their lifetime. Fungal diseases cause damage to trees growth and forest nurseries. As model system for woody plants, poplars are also hosts of a large variety of fungus. The defense strategies to fungus are generally associated with the type of fungus, therefore, the defense strategies of poplar against necrotrophic and biotrophic fungus are different. Poplars initiate constitutive defenses and induced defenses based on recognition of the fungus, hormone signaling network cascades, activation of defense-related genes and transcription factors and production of phytochemicals. The means of sensing fungus invasion in poplars are similar with herbs, both of which are mediated by receptor proteins and resistance (R) proteins, leading to pattern-triggered immunity (PTI) and effector-triggered immunity (ETI), but poplars have evolved some unique defense mechanisms compared with Arabidopsis due to their longevity. In this paper, current researches on poplar defensive responses to necrotrophic and biotrophic fungus, which mainly include the physiological and genetic aspects, and the role of noncoding RNA (ncRNA) in fungal resistance are reviewed. This review also provides strategies to enhance poplar disease resistance and some new insights into future research directions.

Effects of a Microbial Restoration Substrate on Plant Growth and Rhizosphere Microbial Community in a Continuous Cropping Poplar

In poplar cultivation, continuous cropping obstacles affect wood yield and soil-borne diseases, primarily due to structural changes in microbes and fungus infection. The bacterium Bacillus cereus BJS-1-3 has strong antagonistic properties against pathogens that were isolated from the rhizosphere soil of poplars. Poplar rhizospheres were investigated for the effects of Bacillus cereus BJS-1-3 on microbial communities. Three successive generations of soil were used to replant poplar seedlings. BJS-1-3 inoculated poplars were larger, had higher plant height and breast height diameter, and had a greater number of total and culturable bacteria than non-inoculated controls. B. cereus BJS-1-3 inoculated poplar rhizospheres were sequenced, utilizing the Illumina MiSeq platform to analyze changes in diversity and structure. The fungi abundance and diversity in the BJS-1-3 rhizosphere were significantly lower than in the control rhizosphere. In comparison to the control group, Bacillus sp. constituted 2.87% and 2.38% of the total bacterial community, while Rhizoctonia sp. constituted 2.06% and 6.00% of the total fungal community. Among the potential benefits of B. cereus BJS-1-3 in poplar cultivation is that it enhances rhizosphere microbial community structure and facilitates the growth of trees.

Systematic identification and functional characterization of the CFEM proteins in poplar fungus Marssonina brunnea

Proteins containing Common in Fungal Extracellular Membrane (CFEM) domains uniquely exist in fungi and play significant roles in their whole life history. In this study, a total of 11 MbCFEM proteins were identified from Marssonina brunnea f. sp. multigermtubi (MULT), a hemibiotrophic pathogenic fungus on poplars that causes severe leaf diseases. Phylogenic analysis showed that the 11 proteins (MbCFEM1-11) were divided into three clades based on the trans-membrane domain and the CFEM domain. Sequence alignment and WebLogo analysis of CFEM domains verified the amino acids conservatism therein. All of them possess eight cysteines except MbCFEM4 and MbCFEM11, which lack two cysteines each. Six MbCFEM proteins with a signal peptide and without trans-membrane domain were considered as candidate effectors for further functional analysis. Three-dimensional (3D) models of their CFEM domains presented a helical-basket structure homologous to the crucial virulence factor Csa2 of Candida albicans. Afterward, four (MbCFEM1, 6, 8, and 9) out of six candidate effectors were successfully cloned and a yeast signal sequence trap (YSST) assay confirmed their secretion activity. Pathogen challenge assays demonstrated that the transient expression of four candidate MbCFEM effectors in Nicotiana benthamiana promoted Fusarium proliferatum infection, respectively. In an N. benthamiana heterogeneous expression system, MbCFEM1, MbCFEM6, and MbCFEM9 appeared to suppress both BAX/INF1-triggered PCD, whereas MbCFEM8 could only defeat BAX-triggered PCD. Additionally, subcellular localization analysis indicated that the four candidate MbCFEM effectors accumulate in the cell membrane, nucleus, chloroplast, and cytosolic bodies. These results demonstrate that MbCFEM1, MbCFEM6, MbCFEM8, and MbCFEM9 are effectors of M. brunnea and provide valuable targets for further dissection of the molecular mechanisms underlying the poplar-M. brunnea interaction.

Transcriptome analysis of sugarcane reveals differential switching of major defense signaling pathways in response to Sporisorium scitamineum isolates with varying virulent attributes

Sugarcane smut caused by the basidiomycetous fungus Sporisorium scitamineum is one of the most devastating diseases that affect sugarcane production, globally. At present, the most practical and effective management strategy for the disease is the cultivation of resistant cultivars. In this connection, a detailed understanding of the host’s defense mechanism in response to smut isolates with varying degrees of virulence at the molecular level would facilitate the development of reliable and durable smut-resistant sugarcane varieties. Hence, in this study, a comparative whole transcriptome analysis was performed employing Illumina RNA-seq in the smut susceptible cultivar Co 97009 inoculated with two distinct S. scitamineum isolates, Ss97009 (high-virulent) and SsV89101 (low-virulent) during the early phases of infection (2 dpi and 5 dpi) and at the phase of sporogenesis (whip emergence) (60 dpi). Though the differential gene expression profiling identified significant transcriptional changes during the early phase of infection in response to both the isolates, the number of differentially expressed genes (DEGs) were more abundant at 60 dpi during interaction with the high virulent isolate Ss97009, as compared to the low virulent isolate SsV89101. Functional analysis of these DEGs revealed that a majority of them were associated with hormone signaling and the synthesis of defense-related metabolites, suggesting a complex network of defense mechanisms is being operated in response to specific isolates of the smut pathogen. For instance, up-regulation of hormone-related genes, transcription factors, and flavonoid biosynthesis pathway genes was observed in response to both the isolates in the early phase of interaction. In comparison to early phases of infection, only a few pathogenesis-related proteins were up-regulated at 60 dpi in response to Ss97009, which might have rendered the host susceptible to infection. Strikingly, few other carbohydrate metabolism-associated genes like invertases were up-regulated in Ss97009 inoculated plants during the whip emergence stage, representing a shift from sucrose storage to smut symptoms. Altogether, this study established the major switching of defense signaling pathways in response to S. scitamineum isolates with different virulence attributes and provided novel insights into the molecular mechanisms of sugarcane-smut interaction.

Advanced Breeding for Biotic Stress Resistance in Poplar

Poplar is one of the most important forest trees because of its high economic value. Thanks to the fast-growing rate, easy vegetative propagation and transformation, and availability of genomic resources, poplar has been considered the model species for forest genetics, genomics, and breeding. Being a field-growing tree, poplar is exposed to environmental threats, including biotic stresses that are becoming more intense and diffused because of global warming. Current poplar farming is mainly based on monocultures of a few elite clones and the expensive and long-term conventional breeding programmes of perennial tree species cannot face current climate-change challenges. Consequently, new tools and methods are necessary to reduce the limits of traditional breeding related to the long generation time and to discover new sources of resistance. Recent advances in genomics, marker-assisted selection, genomic prediction, and genome editing offer powerful tools to efficiently exploit the Populus genetic diversity and allow enabling molecular breeding to support accurate early selection, increasing the efficiency, and reducing the time and costs of poplar breeding, that, in turn, will improve our capacity to face or prevent the emergence of new diseases or pests.

Deciphering the Genome-Wide Transcriptomic Changes during Interactions of Resistant and Susceptible Genotypes of American Elm with Ophiostoma novo-ulmi

Dutch elm disease (DED), caused by Ophiostoma novo-ulmi (Onu), is a destructive disease of American elm (Ulmus americana L.). The molecular mechanisms of resistance and susceptibility against DED in American elm are still largely uncharacterized. In the present study, we performed a de novo transcriptome (RNA-sequencing; RNA-Seq) assembly of U. americana and compared the gene expression in a resistant genotype, 'Valley Forge', and a susceptible (S) elm genotype at 0 and 96 h post-inoculation of Onu. A total of 85,863 non-redundant unigenes were identified. Compared to the previously characterized U. minor transcriptome, U. americana has 35,290 similar and 55,499 unique genes. The transcriptomic variations between 'Valley Forge' and 'S' were found primarily in the photosynthesis and primary metabolism, which were highly upregulated in the susceptible genotype irrespective of the Onu inoculation. The resistance to DED was associated with the activation of RPM1-mediated effector-triggered immunity that was demonstrated by the upregulation of genes involved in the phenylpropanoids biosynthesis and PR genes. The most significantly enriched gene ontology (GO) terms in response to Onu were response to stimulus (GO:0006950), response to stress (GO:0050896), and secondary metabolic process (GO:0008152) in both genotypes. However, only in the resistant genotype, the defense response (GO:0006952) was among the topmost significantly enriched GO terms. Our findings revealed the molecular regulations of DED resistance and susceptibility and provide a platform for marker-assisted breeding of resistant American elm genotypes.

Contrasting transcriptional responses to Fusarium virguliforme colonization in symptomatic and asymptomatic hosts

The broad host range of Fusarium virguliforme represents a unique comparative system to identify and define differentially induced responses between an asymptomatic monocot host, maize (Zea mays), and a symptomatic eudicot host, soybean (Glycine max). Using a temporal, comparative transcriptome-based approach, we observed that early gene expression profiles of root tissue from infected maize suggest that pathogen tolerance coincides with the rapid induction of senescence dampening transcriptional regulators, including ANACs (Arabidopsis thaliana NAM/ATAF/CUC protein) and Ethylene-Responsive Factors. In contrast, the expression of senescence-associated processes in soybean was coincident with the appearance of disease symptom development, suggesting pathogen-induced senescence as a key pathway driving pathogen susceptibility in soybean. Based on the analyses described herein, we posit that root senescence is a primary contributing factor underlying colonization and disease progression in symptomatic versus asymptomatic host-fungal interactions. This process also supports the lifestyle and virulence of F. virguliforme during biotrophy to necrotrophy transitions. Further support for this hypothesis lies in comprehensive co-expression and comparative transcriptome analyses, and in total, supports the emerging concept of necrotrophy-activated senescence. We propose that F. virguliforme conditions an environment within symptomatic hosts, which favors susceptibility through transcriptomic reprogramming, and as described herein, the induction of pathways associated with senescence during the necrotrophic stage of fungal development.

Jasmonic Acid- and Ethylene-Induced Mitochondrial Alternative Oxidase Stimulates Marssonina brunnea Defense in Poplar

Mitochondrial processes are implicated in plant response to biotic stress caused by viruses, actinomyces, bacteria and pests, but their function in defense against fungal invasion remains unclear. Here, we investigated the role and regulation of mitochondrial alternative oxidase (AOX) in response to black spot disease caused by the hemibiotrophic fungus Marssonina brunnea in poplar. M. brunnea inoculation induced the transcription of the AOX1a gene in the mitochondrial electron transport chain and of jasmonic acid (JA) and ethylene (ET) biosynthetic genes, with the accumulation of these phytohormones in poplar leaf, while inhibiting the transcript amount of the mitochondrial cytochrome c oxidase gene (COX6b) and genes related to salicylic acid (SA). Enhanced AOX reduced poplar susceptibility to M. brunnea with a higher ATP/ADP ratio while the repressed AOX caused the reverse effect. Exogenous JA and 1-aminocyclopropane-1-carboxylic acid (ACC, a biosynthetic precursor of ET) inhibited the transcript amount of COX6b and consequently increased the ratio of AOX pathway to total respiration. Furthermore, the transcription of CYS C1 and CYS D1 genes catalyzing cyanide metabolism was induced, while the cysteine (CYS) substrate levels reduced upon M. brunnea inoculation; exogenous JA and ACC mimicked the effect of M. brunnea infection on cysteine. Exogenous SA enhanced, while JA and ACC reduced, poplar susceptibility to M. brunnea. Moreover, inhibiting AOX completely prohibited JA- and ET-increased tolerance to M. brunnea in poplar. These observations indicate that the JA- and ET-induced mitochondrial AOX pathway triggers defense against M. brunnea in poplar. This effect probably involves cyanide. These findings deepen our understanding of plant-pathogenic fungi interactions.

High-Throughput Sequencing Reveals the Regulatory Networks of Transcriptome and Small RNAs During the Defense Against Marssonina brunnea in Poplar

MicroRNAs are implicated in the adjustment of gene expression in plant response to biotic stresses. However, the regulatory networks of transcriptome and miRNAs are still poorly understood. In the present study, we ascertained the induction of genes for small RNA biosynthesis in poplar defense to a hemibiotrophic fungus Marssonina brunnea and afterward investigated the molecular regulatory networks by performing comprehensive sequencing analysis of mRNAs and small RNAs in M. brunnea-inoculated leaves. Differentially expressed genes in M. brunnea-infected poplar are mainly involved in secondary metabolisms, phytohormone pathways, the recognition of pathogens, and MAPK pathway in the plant, with real-time quantitative PCR (qPCR) validating the mRNA-seq results. Furthermore, differentially expressed miRNAs, such as MIR167_1-6, MIR167_1-12, MIR171_2-3, MIR395-13, MIR396-3, MIR396-16, MIR398-8, and MIR477-6, were identified. Through psRobot and TargetFinder programs, MIR167-1-6, MIR395-13, MIR396-3, MIR396-16, and MIR398-8 were annotated to modulate the expression of genes implicated in transportation, signaling, and biological responses of phytohormones and activation of antioxidants for plant immunity. Besides, validated differentially expressed genes involved in lignin generation, which were phenylalanine ammonia-lyase, ferulate-5-hydroxylase, cinnamyl alcohol dehydrogenase, and peroxidase 11, were selected as targets for the identification of novel miRNAs. Correspondingly, novel miRNAs, such as Novel MIR8567, Novel MIR3228, Novel MIR5913, and Novel MIR6493, were identified using the Mireap online program, which functions in the transcriptional regulation of lignin biosynthesis for poplar anti-fungal response. The present study underlines the roles of miRNAs in the regulation of transcriptome in the anti-fungal response of poplar and provides a new idea for molecular breeding of woody plants.

Key Genes and Genetic Interactions of Plant-Pathogen Functional Modules in Poplar Infected by Marssonina brunnea

Marssonina brunnea, the causative pathogen of Marssonina leaf spot of poplars (MLSP), devastates poplar plantations by forming black spots on leaves and defoliating trees. Although MLSP has been studied for over 30 years, the key genes that function during M. brunnea infection and their effects on plant growth are poorly understood. Here, we used multigene association studies to investigate the effects of key genes in the plant-pathogen interaction pathway, as revealed by transcriptome analysis, on photosynthesis and growth in a natural population of 435 Populus tomentosa individuals. By analyzing transcriptomic changes during three stages of infection, we detected 628 transcription factor genes among the 7,611 differentially expressed genes that might be associated with basal defense responses. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analyses revealed that transcriptomic changes across different stages of infection lead to the reprogramming of metabolic processes possibly related to defense activation. We identified 29,399 common single-nucleotide polymorphisms (SNPs) within 221 full-length genes in plant-pathogen interaction pathways that were significantly associated with photosynthetic and growth traits. We also detected 4,460 significant epistatic pairs associated with stomatal conductance, tree diameter, and tree height. Epistasis analysis uncovered significant interactions between 2,561 SNP-SNP pairs from different functional modules in the plant-pathogen interaction pathway, revealing possible genetic interactions. This analysis revealed many key genes that function during M. brunnea infection and their potential roles in mediating photosynthesis and plant growth, shedding light on genetic interactions between functional modules in the plant-pathogen interaction pathway.

Time-series expression profiling of sugarcane leaves infected with Puccinia kuehnii reveals an ineffective defense system leading to susceptibility

Successful orange rust development on sugarcane can potentially be explained as suppression of the plant immune system by the pathogen or delayed plant signaling to trigger defense responses. Puccinia kuehnii is an obligate biotrophic fungus that infects sugarcane leaves causing a disease called orange rust. It spread out to other countries resulting in reduction of crop yield since its first outbreak. One of the knowledge gaps of that pathosystem is to understand the molecular mechanisms altered in susceptible plants by this biotic stress. Here, we investigated the changes in temporal expression of transcripts in pathways associated with the immune system. To achieve this purpose, we used RNA-Seq to analyze infected leaf samples collected at five time points after inoculation. Differential expression analyses of adjacent time points revealed substantial changes at 12, 48 h after inoculation and 12 days after inoculation, coinciding with the events of spore germination, haustoria post-penetration and post-sporulation, respectively. During the first 24 h, a lack of transcripts involved with resistance mechanisms was revealed by underrepresentation of hypersensitive and defense response related genes. However, two days after inoculation, upregulation of genes involved with immune response regulation provided evidence of some potential defense response. Events related to biotic stress responses were predominantly downregulated in the initial time points, but expression was later restored to basal levels. Genes involved in carbohydrate metabolism showed evidence of repression followed by upregulation, possibly to ensure the pathogen nutritional requirements were met. Our results support the hypothesis that P. kuehnii initially suppressed sugarcane genes involved in plant defense systems. Late overexpression of specific regulatory pathways also suggests the possibility of an inefficient recognition system by a susceptible sugarcane genotype.

Distinctive Gene Expression Profiles and Effectors Consistent With Host Specificity in Two Formae Speciales of Marssonina brunnea

The knowledge on the host specificity of a pathogen underlying an interaction is becoming an urgent necessity for global warming. In this study, the gene expression profiles and the roles of effectors in host specificity were integrally characterized with two formae speciales, multigermtubi and monogermtubi, of a hemibiotrophic pathogen Marssonina brunnea when they were infecting respective susceptible poplar hosts. With a functional genome comparison referring to a de novo transcriptome of M. brunnea and Pathogen-Host Interaction database functional annotations, the multigermtubi strain showed abundant and significant differentially expressed unigenes (DEGs) (more than 40%) in colonizing the initial invasion stage and in the necrotrophic stage. The monogermtubi strain induced less than 10% of DEGs in the initial invasion stage but which abruptly increased to more than 80% DEGs in the necrotrophic stage. Both strains induced the least DEGs in the biotrophic stage compared to the initial invasion and necrotrophic stages. The orthologs of the effector genes Ecp6, PemG1, XEG1, ACE1, and Mg3LysM were exclusively induced by one of the two formae speciales depending on the infection stages. Some unigenes homologous to carbohydrate lytic enzyme genes, especially pectate lyases, were notably induced with multigermtubi forma specialis infection but not expressed in the monogermtubi forma specialis at an earlier infection stage. The extraordinary differences in the functional genome level between the two formae speciales of M. brunnea could be fundamental to exploring their host specificity determinant and evolution. This study also firstly provided the fungal transcriptome of the monogermtubi forma specialis for M. brunnea.

Fusarium virguliform e Transcriptional Plasticity Is Revealed by Host Colonization of Maize versus Soybean

We exploited the broad host range of Fusarium virguliforme to identify differential fungal responses leading to either an endophytic or a pathogenic lifestyle during colonization of maize (Zea mays) and soybean (Glycine max), respectively. To provide a foundation to survey the transcriptomic landscape, we produced an improved de novo genome assembly and annotation of F. virguliforme using PacBio sequencing. Next, we conducted a high-resolution time course of F. virguliforme colonization and infection of both soybean, a symptomatic host, and maize, an asymptomatic host. Comparative transcriptomic analyses uncovered a nearly complete network rewiring, with less than 8% average gene coexpression module overlap upon colonizing the different plant hosts. Divergence of transcriptomes originating from host specific temporal induction genes is central to infection and colonization, including carbohydrate-active enzymes (CAZymes) and necrosis inducing effectors. Upregulation of Zn(II)-Cys6 transcription factors were uniquely induced in soybean at 2 d postinoculation, suggestive of enhanced pathogen virulence on soybean. In total, the data described herein suggest that F. virguliforme modulates divergent infection profiles through transcriptional plasticity.

Fungal canker pathogens trigger carbon starvation by inhibiting carbon metabolism in poplar stems

Carbon starvation is the current leading hypothesis of plant mortality mechanisms under drought stress; recently, it is also used to explain tree die-off in plant diseases. However, the molecular biology of the carbon starvation pathway is unclear. Here, using a punch inoculation system, we conducted transcriptome and physiological assays to investigate pathogen response in poplar stems at the early stages of Botryosphaeria and Valsa canker diseases. Transcriptome assays showed that the majority of differentially expressed genes (DEGs) in stem phloem and xylem, such as genes involved in carbon metabolism and transportation, aquaporin genes (in xylem) and genes related to the biosynthesis of secondary metabolites and the phenylpropanoid pathway (related to lignin synthesis), were downregulated at 7 days after inoculation (DAI). Results also showed that the expression of the majority of disease-resistance genes upregulated in poplar stems, which may be connected with the downregulation expression of the majority of WRKY family genes. Physiological assays showed that transpiration rate decreased but WUE (water use efficiency) increased the 3 and 7 DAI, while the net photosynthetic rate decreased at 11 DAI in Botryosphaeria infected poplars (ANOVA, P < 0.05). The NSC (non-structural carbohydrates) content assays showed that the soluble sugar content of stem phloem samples increased at 3, 7, and 11 DAI that might due to the impede of pathogen infection. However, soluble sugar content of stem xylem and root samples decreased at 11 DAI; in contrast, the starch content unchanged. Therefore, results revealed a chronological order of carbon related molecular and physiological performance: declination of genes involved in carbon and starch metabolism first (at least at 7 DAI), declination of assimilation and carbon reserve (at 11 DAI) second. Results implied a potential mechanism that affects the host carbon reserve, by directly inhibiting the expression of genes involved in carbon metabolism and transport.