Soybean leaves transcriptomic data dissects the phenylpropanoid pathway genes as a defence response against Phakopsora pachyrhizi

Soybean-Phakopsora pachyrhizi interactions: towards the development of next-generation disease-resistant plants

Soybean rust (SBR), caused by the obligate biotrophic fungus Phakopsora pachyrhizi, is a devastating foliar disease threatening soybean production. To date, no commercial cultivars conferring durable resistance to SBR are available. The development of long-lasting SBR resistance has been hindered by the lack of understanding of this complex pathosystem, encompassing challenges posed by intricate genetic structures in both the host and pathogen, leading to a gap in the knowledge of gene-for-gene interactions between soybean and P. pachyrhizi. In this review, we focus on recent advancements and emerging technologies that can be used to improve our understanding of the P. pachyrhizi-soybean molecular interactions. We further explore approaches used to combat SBR, including conventional breeding, transgenic approaches and RNA interference, and how advances in our understanding of plant immune networks, the availability of new molecular tools, and the recent sequencing of the P. pachyrhizi genome could be used to aid in the development of better genetic resistance against SBR. Lastly, we discuss the research gaps of this pathosystem and how new technologies can be used to shed light on these questions and to develop durable next-generation SBR-resistant soybean plants.

RNA-Seq and Comparative Transcriptomic Analyses of Asian Soybean Rust Resistant and Susceptible Soybean Genotypes Provide Insights into Identifying Disease Resistance Genes

Asian soybean rust (ASR), caused by Phakopsora pachyrhizi, is one of the most destructive foliar diseases that affect soybeans. Developing resistant cultivars is the most cost-effective, environmentally friendly, and easy strategy for controlling the disease. However, the current understanding of the mechanisms underlying soybean resistance to P. pachyrhizi remains limited, which poses a significant challenge in devising effective control strategies. In this study, comparative transcriptomic profiling using one resistant genotype and one susceptible genotype was performed under infected and control conditions to understand the regulatory network operating between soybean and P. pachyrhizi. RNA-Seq analysis identified a total of 6540 differentially expressed genes (DEGs), which were shared by all four genotypes. The DEGs are involved in defense responses, stress responses, stimulus responses, flavonoid metabolism, and biosynthesis after infection with P. pachyrhizi. A total of 25,377 genes were divided into 33 modules using weighted gene co-expression network analysis (WGCNA). Two modules were significantly associated with pathogen defense. The DEGs were mainly enriched in RNA processing, plant-type hypersensitive response, negative regulation of cell growth, and a programmed cell death process. In conclusion, these results will provide an important resource for mining resistant genes to P. pachyrhizi infection and valuable resources to potentially pyramid quantitative resistance loci for improving soybean germplasm.

Arogenate dehydratases: unique roles in light-directed development during the seed-to-seedling transition in Arabidopsis thaliana

The seed-to-seedling transition is impacted by changes in nutrient availability and light profiles, but is still poorly understood. Phenylalanine affects early seedling development; thus, the roles of arogenate dehydratases (ADTs), which catalyze phenylalanine formation, were studied in germination and during the seed-to-seedling transition by exploring the impact of light conditions and specific hormone responses in adt mutants of Arabidopsis thaliana. ADT gene expression was assessed in distinct tissues and for light-quality dependence in seedlings for each of the six-member ADT gene family. Mutant adt seedlings were evaluated relative to wild type for germination, photomorphogenesis (blue, red, far red, white light, and dark conditions), anthocyanin accumulation, and plastid development-related phenotypes. ADT proteins are expressed in a light- and tissue-specific manner in transgenic seedlings. Among the analyzed adt mutants, adt3, adt5, and adt6 exhibit significant defects in germination, hypocotyl elongation, and root development responses during the seed-to-seedling transition. Interestingly, adt5 exhibits a light-dependent disruption in plastid development, similar to a phyA mutant. These data indicate interactions between photoreceptors, hormones, and regulation of phenylalanine pools in the process of seedling establishment. ADT5 and ADT6 may play important roles in coordinating hormone and light signals for normal early seedling development.

Identifying Soybean Pod Borer (Leguminivora glycinivorella) Resistance QTLs and the Mechanism of Induced Defense Using Linkage Mapping and RNA-Seq Analysis

The soybean pod borer (Leguminivora glycinivorella) (SPB) is a major cause of soybean (Glycine max L.) yield losses in northeast Asia, thus it is desirable to elucidate the resistance mechanisms involved in soybean response to the SPB. However, few studies have mapped SPB-resistant quantitative trait loci (QTLs) and deciphered the response mechanism in soybean. Here, we selected two soybean varieties, JY93 (SPB-resistant) and K6 (SPB-sensitive), to construct F₂ and F_2:3 populations for QTL mapping and collected pod shells before and after SPB larvae chewed on the two parents to perform RNA-Seq, which can identify stable QTLs and explore the response mechanism of soybean to the SPB. The results show that four QTLs underlying SPB damage to seeds were detected on chromosomes 4, 9, 13, and 15. Among them, qESP-9-1 was scanned in all environments, hence it can be considered a stable QTL. All QTLs explained 0.79 to 6.09% of the phenotypic variation. Meanwhile, 2298 and 3509 DEGs were identified for JY93 and K6, respectively, after the SPB attack, and most of these genes were upregulated. Gene Ontology enrichment results indicated that the SPB-induced and differently expressed genes in both parents are involved in biological processes such as wound response, signal transduction, immune response, and phytohormone pathways. Interestingly, secondary metabolic processes such as flavonoid synthesis were only significantly enriched in the upregulated genes of JY93 after SPB chewing compared with K6. Finally, we identified 18 candidate genes related to soybean pod borer resistance through the integration of QTL mapping and RNA-Seq analysis. Seven of these genes had similar expression patterns to the mapping parents in four additional soybean germplasm after feeding by the SPB. These results provide additional knowledge of the early response and induced defense mechanisms against the SPB in soybean, which could help in breeding SPB-resistant soybean accessions.

A Phakopsora pachyrhizi Effector Suppresses PAMP-Triggered Immunity and Interacts with a Soybean Glucan Endo-1,3-β-Glucosidase to Promote Virulence

Asian soybean rust, caused by the fungus Phakopsora pachyrhizi, is one of the most important diseases affecting soybean production in tropical areas. During infection, P. pachyrhizi secretes proteins from haustoria that are transferred into plant cells to promote virulence. To date, only one candidate P. pachyrhizi effector protein has been characterized in detail to understand the mechanism by which it suppresses plant defenses to enhance infection. Here, we aimed to extend understanding of the pathogenic mechanisms of P. pachyrhizi based on the discovery of host proteins that interact with the effector candidate Phapa-7431740. We demonstrated that Phapa-7431740 suppresses pathogen-associated molecular pattern-triggered immunity (PTI) and that it interacts with a soybean glucan endo-1,3-β-glucosidase (GmβGLU), a pathogenesis-related (PR) protein belonging to the PR-2 family. Structural and phylogenetic characterization of the PR-2 protein family predicted in the soybean genome and comparison to PR-2 family members in Arabidopsis thaliana and cotton, demonstrated that GmβGLU is a type IV β-1,3-glucanase. Transcriptional profiling during an infection time course showed that the GmβGLU mRNA is highly induced during the initial hours after infection, coinciding with peak of expression of Phapa-7431740. The effector was able to interfere with the activity of GmβGLU in vitro, with a dose-dependent inhibition. Our results suggest that Phapa-7431740 may suppress PTI by interfering with glucan endo-1,3-β-glucosidase activity. [Formula: see text] The author(s) have dedicated the work to the public domain under the Creative Commons CC0 "No Rights Reserved" license by waiving all of his or her rights to the work worldwide under copyright law, including all related and neighboring rights, to the extent allowed by law, 2022.

Transcriptomic analysis reveals potential candidate pathways and genes involved in toxin biosynthesis in true toads

Synthesized chemical defenses have broadly evolved across countless taxa and are important in shaping evolutionary and ecological interactions within ecosystems. However, the underlying genomic mechanisms by which these organisms synthesize and utilize their toxins are relatively unknown. Herein, we use comparative transcriptomics to uncover potential toxin synthesizing genes and pathways, as well as interspecific patterns of toxin synthesizing genes across ten species of North American true toads (Bufonidae). Upon assembly and annotation of the ten transcriptomes, we explored patterns of relative gene expression and possible protein-protein interactions across the species to determine what genes and/or pathways may be responsible for toxin synthesis. We also tested our transcriptome dataset for signatures of positive selection to reveal how selection may be acting upon potential toxin producing genes. We assembled high quality transcriptomes of the bufonid parotoid gland, a tissue not often investigated in other bufonid related RNAseq studies. We found several genes involved in metabolic and biosynthetic pathways (e.g. steroid biosynthesis, terpenoid backbone biosynthesis, isoquinoline biosynthesis, glucosinolate biosynthesis) that were functionally enriched and/or relatively expressed across the ten focal species that may be involved in the synthesis of alkaloid and steroid toxins, as well as other small metabolic compounds that cause distastefulness in bufonids. We hope that our study lays a foundation for future studies to explore the genomic underpinnings and specific pathways of toxin synthesis in toads, as well as at the macroevolutionary scale across numerous taxa that produce their own defensive toxins.

Genome-wide identification of ascorbate-glutathione cycle gene families in soybean (Glycine max) reveals gene duplication events and specificity of gene members linked to development and stress conditions

Ascorbate-glutathione (AsA-GSH) cycle plays an important role in tuning beneficial ROS accumulation for intracellular signals and imparts plant tolerance to oxidative stress by detoxifying excess of ROS. Here, we present genome-wide identification of AsA-GSH cycle genes (APX, MDHAR, DHAR, and GR) in several leguminous species and expression analyses in G. max during stress, germination and tissue development. Our data revealed 24 genes in Glycine genus against the maximum of 15 in other leguminous species, which was due to 9 pars of duplicated genes mostly originated from sub/neofunctionalization. Cytosolic APX and MDHAR genes were highly expressed in different tissues and physiological conditions. Germination induced genes encoding AsA-GSH proteins from different cell compartments, whereas vegetative phase (leaves) stimulated predominantly genes related to chloroplast/mitochondria proteins. Moreover, cytosolic APX-1, 2, MDHAR-1a, 1b and GR genes were the primary genes linked to senescence and biotic stresses, while stAPX-a, b and GR (from organelles) were the most abiotic stress related genes. Biotic and abiotic stress tolerant genotypes generally showed increased MDHAR, DHAR and/or GR mRNA levels compared to susceptible genotypes. Overall, these data clarified evolutionary events in leguminous plants and point to the functional specificity of duplicated genes of the AsA-GSH cycle in G. max.

Secondary Metabolites of Pseudomonas aeruginosa LV Strain Decrease Asian Soybean Rust Severity in Experimentally Infected Plants

Asian Soybean Rust (ASR), a disease caused by Phakopsora pachyrhizi, causing yield losses up to 90%. The control is based on the fungicides which may generate resistant fungi. The activation of the plant defense system, should help on ASR control. In this study, secondary metabolites of Pseudomonas aeruginosa LV strain were applied on spore germination and the expression of defense genes in infected soybean plants. The F4A fraction and the pure metabolites were used. In vitro, 10 µg mL⁻¹ of F4A reduced spore germination by 54%, while 100 µg mL⁻¹ completely inhibited. Overexpression of phenylalanine ammonia lyase (PAL), O-methyltransferase (OMT) and pathogenesis related protein-2 (PR-2; glucanases) defense-related genes were detected 24 and 72 h after soybean sprouts were sprayed with an organocopper antimicrobial compound (OAC). Under greenhouse conditions, the best control was observed in plants treated with 60 µg mL⁻¹ of PCA, which reduced ASR severity and lesion frequency by 75% and 43%, respectively. Plants sprayed with 2 and 20 µg mL⁻¹ of F4A also decreased severity (41%) and lesion frequency (32%). The significant reduction in spore germination ASR in plant suggested that the strain of these metabolites are effective against P. pachyrhizi, and they can be used for ASR control.

Untargeted Metabolomics Analysis by UHPLC-MS/MS of Soybean Plant in a Compatible Response to Phakopsora pachyrhizi Infection

Phakopsora pachyrhizi is a biotrophic fungus, causer of the disease Asian Soybean Rust, a severe crop disease of soybean and one that demands greater investment from producers. Thus, research efforts to control this disease are still needed. We investigated the expression of metabolites in soybean plants presenting a resistant genotype inoculated with P. pachyrhizi through the untargeted metabolomics approach. The analysis was performed in control and inoculated plants with P. pachyrhizi using UHPLC-MS/MS. Principal component analysis (PCA) and the partial least squares discriminant analysis (PLS-DA), was applied to the data analysis. PCA and PLS-DA resulted in a clear separation and classification of groups between control and inoculated plants. The metabolites were putative classified and identified using the Global Natural Products Social Molecular Networking platform in flavonoids, isoflavonoids, lipids, fatty acyls, terpenes, and carboxylic acids. Flavonoids and isoflavonoids were up-regulation, while terpenes were down-regulated in response to the soybean–P. pachyrhizi interaction. Our data provide insights into the potential role of some metabolites as flavonoids and isoflavonoids in the plant resistance to ASR. This information could result in the development of resistant genotypes of soybean to P. pachyrhizi, and effective and specific products against the pathogen.

Transcriptome Analysis Reveals the Symbiotic Mechanism of Ustilago esculenta-Induced Gall Formation of Zizania latifolia

Zizania latifolia is a perennial aquatic vegetable, whose symbiosis with the fungus Ustilago esculenta (member of Basidiomycota, class Ustilaginaceae) results in the establishment of swollen gall formations. Here, we analyzed symbiotic relations of Z. latifolia and U. esculenta, using a triadimefon (TDF) treatment and transcriptome sequencing (RNA-seq). Specifically, accurately identify the whole growth cycle of Z. latifolia. Microstructure observations showed that the presence of U. esculenta could be clearly observed after gall formation but was absent after the TDF treatment. A total of 17,541 differentially expressed genes (DEGs) were identified, based on the transcriptome. According to gene ontology term and Kyoto Encyclopedia of Genes and Genomes pathway results, plant hormone signal transduction, and cell wall-loosening factors were all significantly enriched due to U. esculenta infecting Z. latifolia; relative expression levels of hormone-related genes were identified, of which downregulation of indole 3-acetic acid (IAA)-related DEGs was most pronounced in JB_D versus JB_B. The ultra-high performance liquid chromatography analysis revealed that IAA, zeatin+trans zeatin riboside, and gibberellin 3 were increased under U. esculenta infection. Based on our results, we proposed a hormone-cell wall loosening model to study the symbiotic mechanism of gall formation after U. esculenta infects Z. latifolia. Our study thus provides a new perspective for studying the physiological and molecular mechanisms of U. esculenta infection of Z. latifolia causing swollen gall formations as well as a theoretical basis for enhancing future yields of cultivated Z. latifolia.[Formula: see text] The author(s) have dedicated the work to the public domain under the Creative Commons CC0 "No Rights Reserved" license by waiving all of his or her rights to the work worldwide under copyright law, including all related and neighboring rights, to the extent allowed by law. 2021.

Covering Soybean Leaves With Cellulose Nanofiber Changes Leaf Surface Hydrophobicity and Confers Resistance Against Phakopsora pachyrhizi

Asian soybean rust (ASR) caused by Phakopsora pachyrhizi, an obligate biotrophic fungal pathogen, is the most devastating soybean production disease worldwide. Currently, timely fungicide application is the only means to control ASR in the field. We investigated cellulose nanofiber (CNF) application on ASR disease management. CNF-treated leaves showed reduced lesion number after P. pachyrhizi inoculation compared to control leaves, indicating that covering soybean leaves with CNF confers P. pachyrhizi resistance. We also demonstrated that formation of P. pachyrhizi appressoria, and also gene expression related to these formations, such as chitin synthases (CHSs), were significantly suppressed in CNF-treated soybean leaves compared to control leaves. Moreover, contact angle measurement revealed that CNF converts soybean leaf surface properties from hydrophobic to hydrophilic. These results suggest that CNF can change soybean leaf surface hydrophobicity, conferring resistance against P. pachyrhizi, based on the reduced expression of CHSs, as well as reduced formation of pre-infection structures. This is the first study to investigate CNF application to control field disease.

Comparative RNA-Seq profiling of a resistant and susceptible peanut (Arachis hypogaea) genotypes in response to leaf rust infection caused by Puccinia arachidis

The goal of this study was to identify differentially expressed genes (DEGs) responsible for peanut plant (Arachis hypogaea) defence against Puccinia arachidis (causative agent of rust disease). Genes were identified using a high-throughput RNA-sequencing strategy. In total, 86,380,930 reads were generated from RNA-Seq data of two peanut genotypes, JL-24 (susceptible), and GPBD-4 (resistant). Gene Ontology (GO) and KEGG analysis of DEGs revealed essential genes and their pathways responsible for defence response to P. arachidis. DEGs uniquely upregulated in resistant genotype included pathogenesis-related (PR) proteins, MLO such as protein, ethylene-responsive factor, thaumatin, and F-box, whereas, other genes down-regulated in susceptible genotype were Caffeate O-methyltransferase, beta-glucosidase, and transcription factors (WRKY, bZIP, MYB). Moreover, various genes, such as Chitinase, Cytochrome P450, Glutathione S-transferase, and R genes such as NBS-LRR were highly up-regulated in the resistant genotype, indicating their involvement in the plant defence mechanism. RNA-Seq analysis data were validated by RT-qPCR using 15 primer sets derived from DEGs producing high correlation value (R ² = 0.82). A total of 4511 EST-SSRs were identified from the unigenes, which can be useful in evaluating genetic diversity among genotypes, QTL mapping, and plant variety improvement through marker-assisted breeding. These findings will help to understand the molecular defence mechanisms of the peanut plant in response to P. arachidis infection.

The Spread and Transmission of Sweet Potato Virus Disease (SPVD) and Its Effect on the Gene Expression Profile in Sweet Potato

Sweet potato virus disease (SPVD) is the most devastating viral disease in sweet potato (Ipomoea batatas (L.) Lam.), causing substantial yield losses worldwide. We conducted a systemic investigation on the spread, transmission, and pathogenesis of SPVD. Field experiments conducted over two years on ten sweet potato varieties showed that SPVD symptoms first occurred in newly developed top leaves, and spread from adjacent to distant plants in the field. The SPVD incidence was mainly (but not only) determined by the resistance of the varieties planted, and each variety exhibited a characteristic subset of SPVD symptoms. SPVD was not robustly transmitted through friction inoculation, but friction of the main stem might contribute to a higher SPVD incidence rate compared to friction of the leaf and branch tissues. Furthermore, our results suggested that SPVD might be latent in the storage root. Therefore, using virus-free storage roots and cuttings, purposeful monitoring for SPVD according to variety-specific symptoms, and swiftly removing infected plants (especially during the later growth stages) would help control and prevent SPVD during sweet potato production. Comparative transcriptome analysis revealed that numerous genes involved in photosynthesis, starch and sucrose metabolism, flavonoid biosynthesis, and carotenoid biosynthesis were downregulated following SPVD, whereas those involved in monolignol biosynthesis, zeatin biosynthesis, trehalose metabolism, and linoleic acid metabolism were upregulated. Notably, critical genes involved in pathogenesis and plant defense were significantly induced or suppressed following SPVD. These data provide insights into the molecular changes of sweet potato in response to SPVD and elucidate potential SPVD pathogenesis and defense mechanisms in sweet potato. Our study provides important information that can be used to tailor sustainable SPVD control strategies and guide the molecular breeding of SPVD-resistant sweet potato varieties.

Unraveling Asian Soybean Rust metabolomics using mass spectrometry and Molecular Networking approach

Asian Soybean Rust (ASR), caused by the biotrophic fungus Phakopsora pachyrhizi, is a devastating disease with an estimated crop yield loss of up to 90%. Yet, there is a nerf of information on the metabolic response of soybean plants to the pathogen Untargeted metabolomics and Global Natural Products Social Molecular Networking platform approach was used to explore soybean metabolome modulation to P. pachyrhizi infection. Soybean plants susceptible to ASR was inoculated with P. pachyrhizi spore suspension and non-inoculated plants were used as controls. Leaves from both groups were collected 14 days post-inoculation and extracted using different extractor solvent mixtures. The extracts were analyzed on an ultra-high performance liquid chromatography system coupled to high-definition electrospray ionization-mass spectrometry. There was a significant production of defense secondary metabolites (phenylpropanoids, terpenoids and flavonoids) when P. pachyrhizi infected soybean plants, such as putatively identified liquiritigenin, coumestrol, formononetin, pisatin, medicarpin, biochanin A, glyoceollidin I, glyoceollidin II, glyoceollin I, glyoceolidin II, glyoceolidin III, glyoceolidin IV, glyoceolidin VI. Primary metabolites (amino acids, peptides and lipids) also were putatively identified. This is the first report using untargeted metabolomics and GNPS-Molecular Networking approach to explore ASR in soybean plants. Our data provide insights into the potential role of some metabolites in the plant resistance to ASR, which could result in the development of resistant genotypes of soybean to P. pachyrhizi, and effective and specific products against the pathogen.

Network Analysis Combining Proteomics and Metabolomics Reveals New Insights Into Early Responses of Eucalyptus grandis During Rust Infection

Eucalyptus rust is caused by the biotrophic fungus, Austropuccinia psidii, which affects commercial plantations of Eucalyptus, a major raw material for the pulp and paper industry in Brazil. In this manuscript we aimed to uncover the molecular mechanisms involved in rust resistance and susceptibility in Eucalyptus grandis. Epifluorescence microscopy was used to follow the fungus development inside the leaves of two contrasting half-sibling genotypes (rust-resistance and rust-susceptible), and also determine the comparative time-course of changes in metabolites and proteins in plants inoculated with rust. Within 24 h of complete fungal invasion, the analysis of 709 metabolomic features showed the suppression of many metabolites 6 h after inoculation (hai) in the rust-resistant genotype, with responses being induced after 12 hai. In contrast, the rust-susceptible genotype displayed more induced metabolites from 0 to 18 hai time-points, but a strong suppression occurred at 24 hai. Multivariate analyses of genotypes and time points were used to select 16 differential metabolites mostly classified as phenylpropanoid-related compounds. Applying the Weighted Gene Co-Expression Network Analysis (WGCNA), rust-resistant and rust-susceptible genotypes had, respectively, 871 and 852 proteins grouped into 5 and 6 modules, of which 5 and 4 of them were significantly correlated to the selected metabolites. Functional analyses revealed roles for photosynthesis and oxidative-dependent responses leading to temporal activity of metabolites and related enzymes after 12 hai in rust-resistance; while the initial over-accumulation of those molecules and suppression of supporting mechanisms at 12 hai caused a lack of progressive metabolite-enzyme responses after 12 hai in rust-susceptible genotype. This study provides some insights on how E. grandis plants are functionally modulated to integrate secondary metabolites and related enzymes from phenylpropanoid pathway and lead to temporal divergences of resistance and susceptibility responses to rust.