Tanscriptomic Study of the Soybean-Fusarium virguliforme Interaction Revealed a Novel Ankyrin-Repeat Containing Defense Gene, Expression of Whose during Infection Led to Enhanced Resistance to the Fungal Pathogen in Transgenic Soybean Plants

The soybean plasma membrane GmDR1 protein conferring broad-spectrum disease and pest resistance regulates several receptor kinases and NLR proteins

Overexpression of Glycine max disease resistant 1 (GmDR1) exhibits broad-spectrum resistance against Fusarium virguliforme, Heterodera glycines (soybean cyst nematode), Tetranychus urticae (Koch) (spider mites), and Aphis glycines Matsumura (soybean aphids) in soybean. To understand the mechanisms of broad-spectrum immunity mediated by GmDR1, the transcriptomes of a strong and a weak GmDR1-overexpressor following treatment with chitin, a pathogen- and pest-associated molecular pattern (PAMP) common to these organisms, were investigated. The strong and weak GmDR1-overexpressors exhibited altered expression of 6098 and 992 genes, respectively, as compared to the nontransgenic control following chitin treatment. However, only 192 chitin- and 115 buffer-responsive genes exhibited over two-fold changes in expression levels in both strong and weak GmDR1-overexpressors as compared to the control. MapMan analysis of the 192 chitin-responsive genes revealed 64 biotic stress-related genes, of which 53 were induced and 11 repressed as compared to the control. The 53 chitin-induced genes include nine genes that encode receptor kinases, 13 encode nucleotide-binding leucine-rich repeat (NLR) receptor proteins, seven encode WRKY transcription factors, four ethylene response factors, and three MYB-like transcription factors. Investigation of a subset of these genes revealed three receptor protein kinases, seven NLR proteins, and one WRKY transcription factor genes that are induced following F. virguliforme and H. glycines infection. The integral plasma membrane GmDR1 protein most likely recognizes PAMPs including chitin and activates transcription of genes encoding receptor kinases, NLR proteins and defense-related genes. GmDR1 could be a pattern recognition receptor that regulates the expression of several NLRs for expression of PAMP-triggered immunity and/or priming the effector triggered immunity.

Influence of Integrated Management Strategies on Soybean Sudden Death Syndrome (SDS) Root Infection, Foliar Symptoms, Yield and Net Returns

Three soybean field trials were conducted in Indiana to evaluate the integration of seed treatment, cultivar selection, and seeding rate on sudden death syndrome (SDS) root rot, pathogen load in the root, foliar symptoms, yield, and net return. Two soybean cultivars, one moderately resistant and one susceptible to SDS, were planted at three seeding rates (272,277 seeds/ha, 346,535 seeds/ha, and 420,792 seeds/ha). Fluopyram and pydiflumetofen seed treatments were applied to both cultivars, and the cultivars were then compared with a control. Low foliar SDS disease pressure was observed in our study. Seed treatment with either fluopyram or pydiflumetofen and the use of a moderately resistant cultivar decreased Fusarium virguliforme DNA concentration in the root relative to the control and the use of a susceptible cultivar. Fluopyram significantly reduced visual root rot severity by 8.8% and increased yield by 105 kg/ha relative to the control but was not different from pydiflumetofen. However, pydiflumetofen performed the same as the control with respect to root rot severity and yield. Findings from this study support the use of a seed treatment to protect roots from infection and the use of a moderately resistant cultivar planted at a seeding rate of 346,535 seeds/ha to protect yield and maximize net returns when a field has low foliar SDS pressure.

Genome-Wide Association Study (GWAS) of White Mold Resistance in Snap Bean

White mold can result in snap bean yield losses of 90 to 100% when field conditions favor the pathogen. A genome-wide association study (GWAS) was conducted to detect loci significantly associated with white mold resistance in a panel of snap bean (Phaseolus vulgaris L.) cultivars. Two populations of snap bean were used in this study. The first population was the BeanCAP (Coordinated Agriculture Project) Snap Bean Diversity Panel (SBDP) (n = 136), and the second population was the Snap Bean Association Panel (SnAP) (n = 378). SBDP was evaluated for white mold reaction in the field in 2012 and 2013, and SnAP was screened in a greenhouse only using the seedling straw test in 2016. Two reference genomes representing the Andean and Middle American centers of domestication were utilized to align the genotyping-by-sequencing (GBS) data. A GWAS was performed using FarmCPU with one principal component after comparing five models. Thirty-four single-nucleotide polymorphisms (SNPs) significantly associated with white mold resistance were detected. Eleven significant SNPs were identified by the seedling straw test, and 23 significant SNPs were identified by field data. Fifteen SNPs were identified within a 100 kb window containing pentatricopeptide repeat (PPR)-encoding genes, and eleven were close to leucine-rich repeat (LRR)-encoding genes, suggesting that these two classes are of outsized importance for snap bean resistance to white mold.

Comparative Transcriptome Analysis of Onion in Response to Infection by Alternaria porri (Ellis) Cifferi

Purple blotch (PB) is one of the most destructive foliar diseases of onion and other alliums, caused by a necrotrophic fungal pathogen Alternaria porri. There are no reports on the molecular response of onion to PB infection. To elucidate the response of onion to A. porri infection, we consequently carried out an RNAseq analysis of the resistant (Arka Kalyan; AK) and susceptible (Agrifound rose; AFR) genotype after an artificial infection. Through differential expression analyses between control and pathogen-treated plants, we identified 8,064 upregulated and 248 downregulated genes in AFR, while 832 upregulated and 564 downregulated genes were identified in AK. A further significant reprogramming in the gene expression profile was also demonstrated by a functional annotation analysis. Gene ontology (GO) terms, which are particularly involved in defense responses and signaling, are overrepresented in current analyses such as "oxidoreductase activity," "chitin catabolic processes," and "defense response." Several key plant defense genes were differentially expressed on A. porri infection, which includes pathogenesis-related (PR) proteins, receptor-like kinases, phytohormone signaling, cell-wall integrity, cytochrome P450 monooxygenases, and transcription factors. Some of the genes were exclusively overexpressed in resistant genotype, namely, GABA transporter1, ankyrin repeat domain-containing protein, xyloglucan endotransglucosylase/hydrolase, and PR-5 (thaumatin-like). Antioxidant enzyme activities were observed to be increased after infection in both genotypes but higher activity was found in the resistant genotype, AK. This is the first report of transcriptome profiling in onion in response to PB infection and will serve as a resource for future studies to elucidate the molecular mechanism of onion-A. porri interaction and to improve PB resistance in onions.

Current recommendations and novel strategies for sustainable management of soybean sudden death syndrome

The increase in food production requires reduction of the damage caused by plant pathogens, minimizing the environmental impact of management practices. Soil-borne pathogens are among the most relevant pathogens that affect soybean crop yield. Soybean sudden death syndrome (SDS), caused by several distinct species of Fusarium, produces significant yield losses in the leading soybean-producing countries in North and South America. Current management strategies for SDS are scarce since there are no highly resistant cultivars and only a few fungicide seed treatments are available. Because of this, innovative approaches for SDS management need to be developed. Here, we summarize recently explored strategies based on plant nutrition, biological control, priming of plant defenses, host-induced gene silencing, and the development of new SDS-resistance cultivars using precision breeding techniques. Finally, sustainable management of SDS should also consider cultural control practices with minimal environmental impact. © 2021 Society of Chemical Industry.

Arabidopsis non-host resistance PSS30 gene enhances broad-spectrum disease resistance in the soybean cultivar Williams 82

Non-host resistance (NHR), which protects all members of a plant species from non-adapted or non-host plant pathogens, is the most common form of plant immunity. NHR provides the most durable and robust form of broad-spectrum immunity against non-adaptive pathogens pathogenic to other crop species. In a mutant screen for loss of Arabidopsis (Arabidopsis thaliana) NHR against the soybean (Glycine max (L.) Merr.) pathogen Phytophthora sojae, the Phytophthora sojae-susceptible 30 (pss30) mutant was identified. The pss30 mutant is also susceptible to the soybean pathogen Fusarium virguliforme. PSS30 encodes a folate transporter, AtFOLT1, which was previously localized to chloroplasts and implicated in the transport of folate from the cytosol to plastids. We show that two Arabidopsis folate biosynthesis mutants with reduced folate levels exhibit a loss of non-host immunity against P. sojae. As compared to the wild-type Col-0 ecotype, the steady-state folate levels are reduced in the pss1, atfolt1 and two folate biosynthesis mutants, suggesting that folate is required for non-host immunity. Overexpression of AtFOLT1 enhances immunity of transgenic soybean lines against two serious soybean pathogens, the fungal pathogen F. virguliforme and the soybean cyst nematode (SCN) Heterodera glycines. Transgenic lines showing enhanced SCN resistance also showed increased levels of folate accumulation. This study thus suggests that folate contributes to non-host plant immunity and that overexpression of a non-host resistance gene could be a suitable strategy for generating broad-spectrum disease resistance in crop plants.

Diverse Trajectories Drive the Expression of a Giant Virus in the Oomycete Plant Pathogen Phytophthora parasitica

Giant viruses of amoebas, recently classified in the class Megaviricetes, are a group of viruses that can infect major eukaryotic lineages. We previously identified a set of giant virus sequences in the genome of Phytophthora parasitica, an oomycete and a devastating major plant pathogen. How viral insertions shape the structure and evolution of the invaded genomes is unclear, but it is known that the unprecedented functional potential of giant viruses is the result of an intense genetic interplay with their hosts. We previously identified a set of giant virus sequences in the genome of P. parasitica, an oomycete and a devastating major plant pathogen. Here, we show that viral pieces are found in a 550-kb locus and are organized in three main clusters. Viral sequences, namely RNA polymerases I and II and a major capsid protein, were identified, along with orphan sequences, as a hallmark of giant viruses insertions. Mining of public databases and phylogenetic reconstructions suggest an ancient association of oomycetes and giant viruses of amoeba, including faustoviruses, African swine fever virus (ASFV) and pandoraviruses, and that a single viral insertion occurred early in the evolutionary history of oomycetes prior to the Phytophthora–Pythium radiation, estimated at ∼80 million years ago. Functional annotation reveals that the viral insertions are located in a gene sparse region of the Phytophthora genome, characterized by a plethora of transposable elements (TEs), effectors and other genes potentially involved in virulence. Transcription of viral genes was investigated through analysis of RNA-Seq data and qPCR experiments. We show that most viral genes are not expressed, and that a variety of mechanisms, including deletions, TEs insertions and RNA interference may contribute to transcriptional repression. However, a gene coding a truncated copy of RNA polymerase II along a set of neighboring sequences have been shown to be expressed in a wide range of physiological conditions, including responses to stress. These results, which describe for the first time the endogenization of a giant virus in an oomycete, contribute to challenge our view of Phytophthora evolution.

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.

Overexpression of a plasma membrane protein generated broad-spectrum immunity in soybean

Plants fight-off pathogens and pests by manifesting an array of defence responses using their innate immunity mechanisms. Here we report the identification of a novel soybean gene encoding a plasma membrane protein, transcription of which is suppressed following infection with the fungal pathogen, Fusarium virguliforme. Overexpression of the protein led to enhanced resistance against not only against F. virguliforme, but also against spider mites (Tetranychus urticae, Koch), soybean aphids (Aphis glycines, Matsumura) and soybean cyst nematode (Heterodera glycines). We, therefore, name this protein as Glycine max disease resistance 1 (GmDR1; Glyma.10g094800). The homologues of GmDR1 have been detected only in legumes, cocoa, jute and cotton. The deduced GmDR1 protein contains 73 amino acids. GmDR1 is predicted to contain an ecto- and two transmembrane domains. Transient expression of the green fluorescent protein fused GmDR1 protein in soybean leaves showed that it is a plasma membrane protein. We investigated if chitin, a pathogen-associated molecular pattern (PAMP), common to all pathogen and pests considered in this study, can significantly enhance defence pathways among the GmDR1-overexpressed transgenic soybean lines. Chitin induces marker genes of the salicylic- and jasmonic acid-mediated defence pathways, but suppresses the defence pathway regulated by ethylene. Chitin induced SA- and JA-regulated defence pathways may be one of the mechanisms involved in generating broad-spectrum resistance among the GmDR1-overexpressed transgenic soybean lines against two serious pathogens and two pests including spider mites, against which no known resistance genes have been identified in soybean and among the most other crop species.

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.

Comparative genomic analyses of two segregating mutants reveal seven genes likely involved in resistance to Fusarium equiseti in soybean via whole genome re-sequencing

KEY MESSAGE

The candidate genes involved in resistance to Fusarium equiseti in soybean PI 437654 were identified through comparative genomic analyses of mutants via whole genome re-sequencing. The fungus Fusarium infects each stage of the growth and development of soybean and causes soybean (Glycine max (L.)) seed and root rot and seedling damping-off and wilt with a large quantity of annual yield loss worldwide. It is very important to identify the resistant genes in soybean to prevent and control this pathogen. One Fusarium equiseti isolate was previously identified to be incompatible with 'PI 437654' but compatible with a Chinese soybean cultivar 'Zhonghuang 13'. In this study, with the infection of this isolate on the seedling roots of developed PI 437654 mutants, 6 mutants were identified from 500 mutants to significantly alter their phenotypes to F. equiseti compared to wild-type PI 437654. Then, two identified segregating mutants were selected to directly perform whole genome re-sequencing. Finally, through comparative genomic analyses 7 genes including one cluster of 4 nucleotide binding site-leucine-rich repeat genes on one genomic region of chromosome 7, a 60S ribosomal protein L12 gene and 2 uncharacterized genes were identified to be likely involved in the resistance to F. equiseti. These genes will facilitate the breeding of resistant germplasm resources and the identification of resistance of soybean to Fusarium spp.

The Soybean Laccase Gene Family: Evolution and Possible Roles in Plant Defense and Stem Strength Selection

Laccase is a widely used industrial oxidase for food processing, dye synthesis, paper making, and pollution remediation. At present, laccases used by industries come mainly from fungi. Plants contain numerous genes encoding laccase enzymes that show properties which are distinct from that of the fungal laccases. These plant-specific laccases may have better potential for industrial purposes. The aim of this work was to conduct a genome-wide search for the soybean laccase genes and analyze their characteristics and specific functions. A total of 93 putative laccase genes (GmLac) were identified from the soybean genome. All 93 GmLac enzymes contain three typical Cu-oxidase domains, and they were classified into five groups based on phylogenetic analysis. Although adjacent members on the tree showed highly similar exon/intron organization and motif composition, there were differences among the members within a class for both conserved and differentiated functions. Based on the expression patterns, some members of laccase were expressed in specific tissues/organs, while some exhibited a constitutive expression pattern. Analysis of the transcriptome revealed that some laccase genes might be involved in providing resistance to oomycetes. Analysis of the selective pressures acting on the laccase gene family in the process of soybean domestication revealed that 10 genes could have been under artificial selection during the domestication process. Four of these genes may have contributed to the transition of the soft and thin stem of wild soybean species into strong, thick, and erect stems of the cultivated soybean species. Our study provides a foundation for future functional studies of the soybean laccase gene family.

Effective identification of soybean candidate genes involved in resistance to soybean cyst nematode via direct whole genome re-sequencing of two segregating mutants

KEY MESSAGE

Three soybean candidate genes involved in resistance to soybean cyst nematode race 4 were identified via direct whole genome re-sequencing of two segregating mutants. The genes conferring resistance to soybean cyst nematode (SCN) race 4 (Hg type 1.2.3.5.7) in soybean (Glycine max L. Merr.) remains unknown. Next generation sequencing-based methods identify a wide range of targets, it is difficult to identify genes underlying traits. Use of the MutMap and QTL-seq methods to identify trait candidate genes needs backcrossing and is very time-consuming. Here we report a simple method to effectively identify candidate genes involved in resistance to SCN race 4. Two ethane methylsulfonate mutagenized mutants of soybean 'PI 437654', whose SCN race 4-infection phenotype altered, were selected. Six relevant whole genomes were re-sequenced, and then calling of genomic variants (SNPs and InDels) was conducted and compared to 'Williams 82'. The comparison eliminated many genomic variants from the mutant lines that overlapped two non-phenotypic but mutant progeny plants, wild-type PI 437654 and 'Zhonghuang 13'. Finally, only 27 mutations were found among 10 genes. Of these 10 genes, 3 genes, Glyma.09g054000, Glyma.16g065700 and Glyma.18g192200 were overlapped between two phenotypic mutant progeny plants. Therefore, the three genes may be the candidate genes involved in resistance of PI 437654 to soybean cyst nematode race 4. This method simplifies the effective identification of candidate genes.

Mapping Gene Markers for Apple Fruit Ring Rot Disease Resistance Using a Multi-omics Approach

Apple fruit ring rot (FRR), caused by Botryosphaeria dothidea, is a worldwide disease that impacts Asian apple production regions. However, no substantial progress has thus far been made toward the mapping of candidate genes or the development of effective genetic makers. In this five-year study, the resistance of 1,733 F1 hybrids from the cross ‘Jonathan’ × ‘Golden Delicious’ was phenotyped by non-wounding inoculation with four B. dothidea isolates. We first conducted systematic comparison of different analytic strategies for bulk segregant analysis by re-sequencing (BSA-Seq) and obtained suitable one for outbreeding species such as Malus. Forty-six quantitative trait loci (QTL) for resistance/susceptibility to the four isolates, including one QTL ‘hotspot’ on chromosome 14, were identified via BSA-Seq. Using integrated multi-omics strategies including RNA-sequencing, parental re-sequencing, BSA-Seq and meta-analysis of RNA-sequencing, fifty-seven candidate genes and corresponding functional mutations from the QTL were predicted. Functional mutations located on the candidate genes were validated using kompetitive allele-specific PCR in hybrids and Malus germplasm accessions with extremely resistant/susceptible phenotypes. Ten effective markers for apple ring rot were developed. The results provide an example of rapid candidate gene mapping for complex traits in outbreeding species.

Arabidopsis Novel Glycine-Rich Plasma Membrane PSS1 Protein Enhances Disease Resistance in Transgenic Soybean Plants

Nonhost resistance is defined as the immunity of a plant species to all nonadapted pathogen species. Arabidopsis (Arabidopsis thaliana) ecotype Columbia-0 is nonhost to the oomycete plant pathogen Phytophthora sojae and the fungal plant pathogen Fusarium virguliforme that are pathogenic to soybean (Glycine max). Previously, we reported generating the pss1 mutation in the pen1-1 genetic background as well as genetic mapping and characterization of the Arabidopsis nonhost resistance Phytophthora sojae-susceptible gene locus, PSS1 In this study, we identified six candidate PSS1 genes by comparing single-nucleotide polymorphisms of (1) the bulked DNA sample of seven F2:3 families homozygous for the pss1 allele and (2) the pen1-1 mutant with Columbia-0. Analyses of T-DNA insertion mutants for each of these candidate PSS1 genes identified the At3g59640 gene encoding a glycine-rich protein as the putative PSS1 gene. Later, complementation analysis confirmed the identity of At3g59640 as the PSS1 gene. PSS1 is induced following P. sojae infection as well as expressed in an organ-specific manner. Coexpression analysis of the available transcriptomic data followed by reverse transcriptase-polymerase chain reaction suggested that PSS1 is coregulated with ATG8a (At4g21980), a core gene in autophagy. PSS1 contains a predicted single membrane-spanning domain. Subcellular localization study indicated that it is an integral plasma membrane protein. Sequence analysis suggested that soybean is unlikely to contain a PSS1-like defense function. Following the introduction of PSS1 into the soybean cultivar Williams 82, the transgenic plants exhibited enhanced resistance to F. virguliforme, the pathogen that causes sudden death syndrome.