A Comprehensive Analysis of the Transcriptomes of Marssonina brunnea and Infected Poplar Leaves to Capture Vital Events in Host-Pathogen Interactions

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.

Asparagine-rich protein (NRP) mediates stress response by regulating biosynthesis of plant secondary metabolites in Arabidopsis

The plant-specific stress response protein NRP (asparagine-rich protein) is characterized by an asparagine-rich domain at its N-terminus and a conserved development and cell death (DCD) domain at its C-terminus. Previous transcriptional studies and phenotypic analyses have demonstrated the involvement of NRP in response to severe stress conditions, such as high salt and ER Endoplasmic reticulum-stress. We have recently identified distinct roles for NRP in biotic- and abiotic-stress signaling pathways, in which NRP interacts with different signaling proteins to change their subcellular localizations and stability. Here, to further explore the function of NRP, a transcriptome analysis was carried out on nrp1nrp2 knock-out lines at different life stages or under different growing conditions. The most significant changes in the transcriptome at both stages and conditions turned out to be the induction of the synthesis of secondary metabolites (SMs). Such an observation implicates that NRP is a general stress-responsive protein involved in various challenges faced by plants during their life cycle, which might involve a broad alteration in the distribution of SMs.

Characterizations of an Emerging Disease: Apple Blotch Caused by Diplocarpon coronariae (syn. Marssonina coronaria) in the Mid-Atlantic United States

Apple orchards with minimal or reduced fungicide inputs in the Mid-Atlantic region of the United States have experienced outbreaks of severe premature defoliation with symptoms that matched those of apple blotch disease (ABD) caused by Diplocarpon coronariae. Fungal isolates obtained from symptomatic apple leaves and fruit produced uniform slow-growing, dark-gray colonies on peptone potato dextrose agar and had conidia. Internal transcribed spacer DNA sequences matched with D. coronariae and Koch's postulates were fulfilled when typical ABD symptoms occurred when reinoculated onto apple leaves and fruit. Spore dispersal in nonfungicide-treated orchards detected with quantitative PCR was low in early spring and dropped to undetectable levels in late May and early June before rising exponentially to highs in July and August, which coincided with symptom development. Only low spore numbers were detected in fungicide-treated orchards and nearby forests. In preliminary fungicide tests, fluxapyroxad, thiophanate methyl, and difenoconazole effectively inhibited mycelial growth of isolates in vitro. When apple cultivars Fuji and Honeycrisp were inoculated with D. coronariae, Honeycrisp showed delayed onset of symptoms and lower disease severity, and the transcription profile of seven host defense-related genes showed that PR-2, PR-8, LYK4, and CERK1 were highly induced in Honeycrisp at 2 and 5 days postinoculation. This is the first report of ABD in the Mid-Atlantic United States, which includes studies of seasonal D. coronariae spore dispersal patterns, preliminary fungicide efficacy, and host defense-related gene expression to assist development of best ABD management practices.

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.

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.

Unified feature association networks through integration of transcriptomic and proteomic data

High-throughput multi-omics studies and corresponding network analyses of multi-omic data have rapidly expanded their impact over the last 10 years. As biological features of different types (e.g. transcripts, proteins, metabolites) interact within cellular systems, the greatest amount of knowledge can be gained from networks that incorporate multiple types of -omic data. However, biological and technical sources of variation diminish the ability to detect cross-type associations, yielding networks dominated by communities comprised of nodes of the same type. We describe here network building methods that can maximize edges between nodes of different data types leading to integrated networks, networks that have a large number of edges that link nodes of different-omic types (transcripts, proteins, lipids etc). We systematically rank several network inference methods and demonstrate that, in many cases, using a random forest method, GENIE3, produces the most integrated networks. This increase in integration does not come at the cost of accuracy as GENIE3 produces networks of approximately the same quality as the other network inference methods tested here. Using GENIE3, we also infer networks representing antibody-mediated Dengue virus cell invasion and receptor-mediated Dengue virus invasion. A number of functional pathways showed centrality differences between the two networks including genes responding to both GM-CSF and IL-4, which had a higher centrality value in an antibody-mediated vs. receptor-mediated Dengue network. Because a biological system involves the interplay of many different types of molecules, incorporating multiple data types into networks will improve their use as models of biological systems. The methods explored here are some of the first to specifically highlight and address the challenges associated with how such multi-omic networks can be assembled and how the greatest number of interactions can be inferred from different data types. The resulting networks can lead to the discovery of new host response patterns and interactions during viral infection, generate new hypotheses of pathogenic mechanisms and confirm mechanisms of disease.

Integrated Transcriptome Analysis Reveals Plant Hormones Jasmonic Acid and Salicylic Acid Coordinate Growth and Defense Responses upon Fungal Infection in Poplar

Plants have evolved a sophisticated system to respond to various stresses. Fungal attack or infection is one of the most important biotic stresses for most plants. During the defense response to fungal infection, the plant hormones jasmonic acid (JA) and salicylic acid (SA) play critical roles. Here, gene expression data on JA/SA treatments and Melampsora larici-populina (MLP) infection were generated. Integrated transcriptome analyses of these data were performed, and 943 genes in total were identified as common responsive genes (CRG). Gene ontology (GO) term analysis revealed that the genes from CRG are generally involved in the processes of stress responses, metabolism, and growth and development. The further cluster analysis of the CRG identified a set of core genes that are involved in the JA/SA-mediated response to fungal defense with distinct gene expression profiles upon JA/SA treatment, which highlighted the different effects of these two hormones on plant fungal defenses. The modifications of several pathways relative to metabolism, biotic stress, and plant hormone signal pathways suggest the possible roles of JA/SA on the regulation of growth and defense responses. Co-expression modules (CMs) were also constructed using the poplar expression data on JA, SA, M. larici-populina, Septoria musiva, and Marssonina brunnea treatment or infection. A total of 23 CMs were constructed, and different CMs clearly exhibited distinct biological functions, which conformably regulated the concerted processes in response to fungal defense. Furthermore, the GO term analysis of different CMs confirmed the roles of JA and SA in regulating growth and defense responses, and their expression profiles suggested that the growth ability was reduced when poplar deployed defense responses. Several transcription factors (TFs) among the CRG in the co-expression network were proposed as hub genes in regulating these processes. According to this study, our data finely uncovered the possible roles of JA/SA in regulating the balance between growth and defense responses by integrating multiple hormone signaling pathways. We were also able to provide more knowledge on how the plant hormones JA/SA are involved in the regulation of the balance between growth and plant defense.

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

Marssonina leaf spot disease of poplar (MLDP), caused by the hemibiotrophic pathogen Marssonina brunnea, frequently results in damage to many poplar species. In nature, two formae speciales of M. brunnea exist that are susceptible to different poplar subgenera. Marssonina brunnea f. sp. monogermtubi infects poplar hosts from Populus sect. Aigeiros (Aig), while M. brunnea f. sp. multigermtubi always infects poplar hosts from Populus sect. Leuce Duby (Leu). Based on the fungal penetration structures, a comprehensive transcriptomic approach was used to investigate the gene expression patterns of these two poplar subgenera at three crucial infection stages. MLDP significantly altered the expression patterns of many genes involved in mitogen activated protein kinase (MAPKs) and calcium signaling, transcription factors, primary and secondary metabolism, and other processes in both poplar subgenera. However, major differences in gene expression were also observed between the two poplar subgenera. Aig was most responsive at the initial infection stage, while Leu largely interacted with M. brunnea at the necrotrophic phase. Furthermore, the differentially expressed genes (DEGs) involved in pathways related to biotic stress also differed substantially between the two poplar subgenera. Further analysis indicated that the genes involved in cell wall metabolism and phenylpropanoid metabolism were differentially expressed in the progression of the disease. By examining the expression patterns of genes related to the defense against disease, we found that several genes annotated with causing hypersensitive cell death were upregulated at the necrotrophic phase of MLDP, inferring that plant immune response potentially happened at this infection stage. The present research elucidated the potential molecular differences between the two susceptible interaction systems in MLDP and provided novel insight into the temporal regulation of genes during the susceptible response. To the best of our knowledge, this study also constitutes the first to reveal the molecular mechanisms of poplar in response to the transition of hemibiotrophic fungal pathogens from the biotrophic phase to the necrotrophic phase.