Discriminant haplotypes of avirulence genes of Phytophthora sojae lead to a molecular assay to predict phenotypes

The soybean–Phytophthora sojae interaction operates on a gene-for-gene relationship, where the product of a resistance gene (Rps) in the host recognizes that of an avirulence gene (Avr) in the pathogen to generate an incompatible reaction. To exploit this form of resistance, one must match with precision the appropriate Rps gene with the corresponding Avr gene. Currently, this association is evaluated by phenotyping assays that are labour-intensive and often imprecise. To circumvent this limitation, we sought to develop a molecular assay that would reveal the avirulence allele of the seven main Avr genes (Avr1a, Avr1b, Avr1c, Avr1d, Avr1k, Avr3a, and Avr6) in order to diagnose with precision the pathotypes of P. sojae isolates. For this purpose, we analysed the genomic regions of these Avr genes in 31 recently sequenced isolates with different virulence profiles and identified discriminant mutations between avirulence and virulence alleles. Specific primers were designed to generate amplicons of a distinct size, and polymerase chain reaction conditions were optimized in a final assay of two parallel runs. When tested on the 31 isolates of known virulence, the assay accurately revealed all avirulence alleles. The test was further assessed and compared to a phenotyping assay on 25 isolates of unknown virulence. The two assays matched in 97% (170/175) of the interactions studied. Interestingly, the sole cases of discrepancy were obtained with Avr3a, which suggests a possible imperfect interaction with Rps3a. This molecular assay offers a powerful and reliable tool to exploit and study with greater precision soybean resistance against P. sojae.     

A CRISPR/Cas9-mediated in situ complementation method for Phytophthora sojae mutants

Phytophthora sojae is an important model species for oomycete functional genomics research. Recently, a CRISPR/Cas9-mediated genome-editing technology has been successfully established in P. sojae, which has been rapidly and widely applied in oomycete research. However, there is an emerging consensus in the biological community that a complete functional gene research system is needed such as developed in the investigations in functional complementation carried out in this study. We report the development of an in situ complementation method for accurate restoration of the mutated gene. We targeted a regulatory B-subunit of protein phosphatase 2A (PsPP2Ab1) to verify this knockout and subsequent complementation system. We found that the deletion of PsPP2Ab1 in P. sojae leads to severe defects in vegetative hyphal growth, soybean infection, and loss of the ability to produce sporangia. Subsequently, the reintroduction of PsPP2Ab1 into the knockout mutant remedied all of the deficiencies. This study demonstrates the successful implementation of an in situ complementation system by CRISPR/Cas9, which will greatly accelerate functional genomics research of oomycetes in the post-genomic era.     

The Pectin Methylesterase Gene Complement of Phytophthora sojae: Structural and Functional Analyses,and the Evolutionary Relationships with Its Oomycete Homologs

Phytophthora sojae is an oomycete pathogen that causes the disease known as root and stem rot in soybean plants, frequently leading to massive economic damage. Additionally, P. sojae is increasingly being utilized as a model for phytopathogenic oomycete research. Despite the economic and scientific importance of P. sojae, the mechanism by which it penetrates the host roots is not yet fully understood. It has been found that oomycetes are not capable of penetrating the cell wall solely through mechanical force, suggesting that alternative factors facilitate breakdown of the host cell wall. Pectin methylesterases have been suggested to be important for Phytophthora pathogenicity, but no data exist on their role in the P. sojae infection process. We have scanned the newly revised version of the annotated P. sojae genome for the presence of putative pectin methylesterases genes and conducted a sequence analysis of all gene models found. We also searched for potential regulatory motifs in the promoter region of the proposed P. sojae models, and investigated the gene expression levels throughout the early course of infection on soybean plants. We found that P. sojae contains a large repertoire of pectin methylesterase-coding genes and that most of these genes display similar motifs in the promoter region, indicating the possibility of a shared regulatory mechanism. Phylogenetic analyses confirmed the evolutionary relatedness of the pectin methylesterase-coding genes within and across Phytophthora spp. In addition, the gene duplication events that led to the emergence of this gene family appear to have occurred prior to many speciation events in the genus Phytophthora. Our results also indicate that the highest levels of expression occurred in the first 24 hours post inoculation, with expression falling after this time. Our study provides evidence that pectin methylesterases may be important for the early action of the P. sojae infection process.     

Zinc finger nuclease‐mediated targeting of multiple transgenes to an endogenous soybean genomic locus via non‐homologous end joining

Emerging genome editing technologies hold great promise for the improvement of agricultural crops. Several related genome editing methods currently in development utilize engineered, sequence‐specific endonucleases to generate DNA double strand breaks (DSBs) at user‐specified genomic loci. These DSBs subsequently result in small insertions/deletions (indels), base substitutions or incorporation of exogenous donor sequences at the target site, depending on the application. Targeted mutagenesis in soybean (Glycine max) via non‐homologous end joining (NHEJ)‐mediated repair of such DSBs has been previously demonstrated with multiple nucleases, as has homology‐directed repair (HDR)‐mediated integration of a single transgene into target endogenous soybean loci using CRISPR/Cas9. Here we report targeted integration of multiple transgenes into a single soybean locus using a zinc finger nuclease (ZFN). First, we demonstrate targeted integration of biolistically delivered DNA via either HDR or NHEJ to the FATTY ACID DESATURASE 2‐1a (FAD2‐1a) locus of embryogenic cells in tissue culture. We then describe ZFN‐ and NHEJ‐mediated, targeted integration of two different multigene donors to the FAD2‐1a locus of immature embryos. The largest donor delivered was 16.2 kb, carried four transgenes, and was successfully transmitted to T₁ progeny of mature targeted plants obtained via somatic embryogenesis. The insertions in most plants with a targeted, 7.1 kb, NHEJ‐integrated donor were perfect or near‐perfect, demonstrating that NHEJ is a viable alternative to HDR for gene targeting in soybean. Taken together, these results show that ZFNs can be used to generate fertile transgenic soybean plants with NHEJ‐mediated targeted insertions of multigene donors at an endogenous genomic locus.     

Roles of small RNAs in soybean defense against Phytophthora sojae infection

The genus Phytophthora consists of many notorious pathogens of crops and forestry trees. At present, battling Phytophthora diseases is challenging due to a lack of understanding of their pathogenesis. We investigated the role of small RNAs in regulating soybean defense in response to infection by Phytophthora sojae, the second most destructive pathogen of soybean. Small RNAs, including microRNAs (miRNAs) and small interfering RNAs (siRNAs), are universal regulators that repress target gene expression in eukaryotes. We identified known and novel small RNAs that differentially accumulated during P. sojae infection in soybean roots. Among them, miR393 and miR166 were induced by heat‐inactivated P. sojae hyphae, indicating that they may be involved in soybean basal defense. Indeed, knocking down the level of mature miR393 led to enhanced susceptibility of soybean to P. sojae; furthermore, the expression of isoflavonoid biosynthetic genes was drastically reduced in miR393 knockdown roots. These data suggest that miR393 promotes soybean defense against P. sojae. In addition to miRNAs, P. sojae infection also resulted in increased accumulation of phased siRNAs (phasiRNAs) that are predominantly generated from canonical resistance genes encoding nucleotide binding‐leucine rich repeat proteins and genes encoding pentatricopeptide repeat‐containing proteins. This work identifies specific miRNAs and phasiRNAs that regulate defense‐associated genes in soybean during Phytophthora infection.     

Resistance to Phytophthora pathogens is dependent on gene silencing pathways in plants

RNA silencing is one of the main defence mechanisms employed by plants to fight pathogens. p19 protein encoded by the tomato bushy stunt virus (TBSVp19) is known as a suppressor of RNA silencing via siRNA sequestration to prevent the assembly of RISC. To better understand the impact of TBSVp19 on silencing and its roles in Phytophthora pathogens, we used the transient expression assay in Nicotiana benthamiana and found that the leaves expressing TBSVp19 were more susceptible to Phytophthora parasitica. Furthermore, we demonstrated that TBSVp19‐mediated plant susceptibility in N. benthamiana is dependent on RNA‐dependent RNA polymerase 6 (RDR6). We also tested the role of RNA silencing in resistance of soybean hairy roots to Phytophthora. The lesion size induced by P. sojae on TBSVp19‐expressing soybean hairy roots was slightly, but significantly larger than GFP‐expressing soybean hairy roots. Finally, the Arabidopsis gene silencing mutants ago1‐27, zip‐1, sgs3‐11 and rdr6‐11 were also examined for their resistance to P. parasitica. The results clearly showed that resistance levels of the mutants were visibly reduced compared with the wild type. Taken together, these results suggest that the gene silencing system in plants is essential for resistance to Phytophthora pathogens.     

The soybean-Phytophthora resistance locus <Emphasis Type="Italic">Rps1</Emphasis>-k encompasses coiled coil-nucleotide binding-leucine rich repeat-like genes and repetitive sequences

Background  

A series of Rps (resistance to P ytophthora s ojae) genes have been protecting soybean from the root and stem rot disease caused by the Oomycete pathogen, Phytophthora sojae. Five Rps genes were mapped to the Rps1 locus located near the 28 cM map position on molecular linkage group N of the composite genetic soybean map. Among these five genes, Rps1-k was introgressed from the cultivar, Kingwa. Rps1-k has been providing stable and broad-spectrum Phytophthora resistance in the major soybean-producing regions of the United States. Rps1-k has been mapped and isolated. More than one functional Rps1-k gene was identified from the Rps1-k locus. The clustering feature at the Rps1-k locus might have facilitated the expansion of Rps1-k gene numbers and the generation of new recognition specificities. The Rps1-k region was sequenced to understand the possible evolutionary steps that shaped the generation of Phytophthora resistance genes in soybean.     

Molecular mapping of two genes conferring resistance to Phytophthora sojae in a soybean landrace PI 567139B

Phytophthora root and stem rot (PRR), caused by the soil-borne oomycete pathogen Phytophthora sojae, is one of the most destructive diseases of soybean. PRR can be effectively controlled by race-specific genes conferring resistance to P. sojae (Rps). However, the Rps genes are usually non-durable, as populations of P. sojae are highly diverse and quick to adapt, and can be overcome 8–15 years after deployment. Thus, it is important to identify novel Rps genes for development of resistant soybean cultivars. PI 567139B is a soybean landrace carrying excellent resistance to nearly all predominant P. sojae races in Indiana. A mapping population consisting of 245 F₂ individuals and 403 F_2:3 families was developed from a cross between PI 567139B and the susceptible cultivar ‘Williams’, and used to dissect the resistance carried by PI 567139B. We found that the resistance in PI 567139B was conferred by two independent Rps genes, designated RpsUN1 and RpsUN2. The former was mapped to a 6.5 cM region between SSR markers Satt159 and BARCSOYSSR_03_0250 that spans the Rps1 locus on chromosome 3, while the latter was mapped to a 3.0 cM region between BARCSOYSSR_16_1275 and Sat_144, approximately 3.0–3.4 cM upstream of Rps2 on chromosome 16. According to the ‘Williams 82’ reference genome sequence, both regions are highly enriched with NBS-LRR genes. Marker assisted resistance spectrum analyses of these genes with 16 isolates of P. sojae, in combination with the mapping results, suggested that RpsUN1 was likely to be a novel allele at the Rps1 locus, while RpsUN2 was more likely to be a novel Rps gene.     

Peroxidase Activity in Soybeans Following Inoculation with Phytophthora sojae

The effects of race-specific resistance as conditioned by Rps genes (rps, Rps1-k, Rps2, Rps3, Rps6) in two genetic backgrounds (Williams & Harosoy) on accumulation of soluble peroxidases were determined by a soybean peroxidase capture assay (SPCA) after inoculation with P. sojae races 2, 7, or 25. Peroxidase activity increased in all isolines during the 72 h after inoculation, but reactions varied depending on time after inoculation, genetic background, Rps gene and P. sojae race. Peroxidase activity was higher in race-specific resistant than in susceptible reactions at 72 h. after inoculation, except for plants with the Rps2 gene which confers a unique form of root resistance in addition to the whole plant race-specific resistance. Williams isolines had larger increases in peroxidase activity than Harosoy isolines when data were averaged across Rps genes, and was most evident when plants were inoculated with race 2. When soybeans were inoculated with race 7 Rps1-k resistant plants had the highest increase in peroxidase activity, but Rps2 susceptible plants had a significantly higher peroxidase activity than plants with rps, Rps3, and Rps6 that were also susceptible. Results from inoculations with race 25 were somewhat different, Rps2 resistant plants had the highest increase in peroxidase activity; however, plants with the Rps3 or Rps6 gene that were also resistant did not have a significantly higher peroxidase activity than susceptible plants with the rps or Rps1-k gene.     

Editing of an effector gene promoter sequence impacts plant‐Phytophthora interaction

Pathogen avirulence (Avr) effectors interplay with corresponding plant resistance (R) proteins and activate robust plant immune responses. Although the expression pattern of Avr genes has been tied to their functions for a long time, it is still not clear how Avr gene expression patterns impact plant‐microbe interactions. Here, we selected PsAvr3b, which shows a typical effector gene expression pattern from a soybean root pathogen Phytophthora sojae. To modulate gene expression, we engineered PsAvr3b promoter sequences by in situ substitution with promoter sequences from Actin (constitutive expression), PsXEG1 (early expression), and PsNLP1 (later expression) using the CRISPR/Cas9. PsAvr3b driven by different promoters resulted in distinct expression levels across all the tested infection time points. Importantly, those mutants with low PsAvr3b expression successfully colonized soybean plants carrying the cognate R gene Rps3b. To dissect the difference in plant responses to the PsAvr3b expression level, we conducted RNA‐sequencing of different infection samples at 24 h postinfection and found soybean immune genes, including a few previously unknown genes that are associated with resistance. Our study highlights that fine‐tuning in Avr gene expression impacts the compatibility of plant disease and provides clues to improve crop resistance in disease control management.     

The Activation of Phytophthora Effector Avr3b by Plant Cyclophilin is Required for the Nudix Hydrolase Activity of Avr3b

Plant pathogens secrete an arsenal of effector proteins to impair host immunity. Some effectors possess enzymatic activities that can modify their host targets. Previously, we demonstrated that a Phytophthora sojae RXLR effector Avr3b acts as a Nudix hydrolase when expressed in planta; and this enzymatic activity is required for full virulence of P. sojae strain P6497 in soybean (Glycine max). Interestingly, recombinant Avr3b produced by E. coli does not have the hydrolase activity unless it was incubated with plant protein extracts. Here, we report the activation of Avr3b by a prolyl-peptidyl isomerase (PPIase), cyclophilin, in plant cells. Avr3b directly interacts with soybean cyclophilin GmCYP1, which activates the hydrolase activity of Avr3b in a PPIase activity-dependent manner. Avr3b contains a putative Glycine-Proline (GP) motif; which is known to confer cyclophilin-binding in other protein substrates. Substitution of the Proline (P132) in the putative GP motif impaired the interaction of Avr3b with GmCYP1; as a result, the mutant Avr3b^P132A can no longer be activated by GmCYP1, and is also unable to promote Phytophthora infection. Avr3b elicits hypersensitive response (HR) in soybean cultivars producing the resistance protein Rps3b, but Avr3b^P132A lost its ability to trigger HR. Furthermore, silencing of GmCYP1 rendered reduced cell death triggered by Avr3b, suggesting that GmCYP1-mediated Avr3b maturation is also required for Rps3b recognition. Finally, cyclophilins of Nicotiana benthamiana can also interact with Avr3b and activate its enzymatic activity. Overall, our results demonstrate that cyclophilin is a “helper” that activates the enzymatic activity of Avr3b after it is delivered into plant cells; as such, cyclophilin is required for the avirulence and virulence functions of Avr3b.     

Novel quantitative trait loci for partial resistance to Phytophthora sojae in soybean PI 398841

Phytophthora root and stem rot caused by Phytophthora sojae Kaufmann and Gerdemann is one of the most severe soybean [Glycine max (L.) Merr] diseases in the USA. Partial resistance is as effective in managing this disease as single-gene (Rps gene)-mediated resistance and is more durable. The objective of this study was to identify quantitative trait loci (QTL) associated with partial resistance to P. sojae in PI 398841, which originated from South Korea. A population of 305 F_7:8 recombinant inbred lines derived from a cross of OX20-8 × PI 398841 was used to evaluate partial resistance against P. sojae isolate C2S1 using a tray test. Composite interval mapping using a genome-wide logarithm of odd (LOD) threshold detected three QTL on chromosomes 1, 13, and 18, which individually explained 4–16 % of the phenotypic variance. Seven additional QTL, accounting for 2–3 % of phenotypic variance each, were identified using chromosome-wide LOD thresholds. Seven of the ten QTL for resistance to P. sojae were contributed by PI 398841. Seven QTL co-localized with known Rps genes and previously reported QTL for soil-borne root pathogens, isoflavone, and seed oil. Three QTL on chromosomes 3, 13, and 18 co-localized with known Rps genes, but PI 398841 did not exhibit an Rps gene-mediated resistance response following inoculation with 48 different isolates of P. sojae. PI 398841 is potentially a source of novel genes for improving soybean cultivars for partial resistance to P. sojae.     

Virulence of Phytophthora sojae in the Pampeana Subregion of Argentina from 1998 to 2004

Phytophthora root rot is one of the most serious diseases of soybeans in Argentina. Surveys of commercial fields and trial plots of soybean were conducted throughout the northern Pampeana subregion (Argentina) between 1998 and 2004. A total of 193 isolates of Phytophthora sojae were collected and classified into races or virulence formulae. Among the 173 isolates tested on 8 differentials, 42 different pathotypes were detected, including 18 described races. Races 1, 4, 5, 7, 9, 13, 23 and 24 were found in both plants and soils, whereas races 2, 3, 6, 8, 11, 14, 15, 17, 43 and 44 were only isolated from plants. An additional 19 pathotypes were described from 20 isolates tested in Canada on the expanded set of 14 differential cultivars. Currently, all Rps genes/alleles associated with resistance have been defeated, indicating an increased complexity of virulence within the P. sojae populations in the region. The great increase in virulence complexity found in this study is most likely a result of a long period of continuous production of soybean cultivars with Rps genes and the extensive adoption of the no‐tillage system.     

CRISPR‐Cas9‐mediated efficient directed mutagenesis and RAD51‐dependent and RAD51‐independent gene targeting in the moss Physcomitrella patens

The ability to address the CRISPR‐Cas9 nuclease complex to any target DNA using customizable single‐guide RNAs has now permitted genome engineering in many species. Here, we report its first successful use in a nonvascular plant, the moss Physcomitrella patens. Single‐guide RNAs (sgRNAs) were designed to target an endogenous reporter gene, PpAPT, whose inactivation confers resistance to 2‐fluoroadenine. Transformation of moss protoplasts with these sgRNAs and the Cas9 coding sequence from Streptococcus pyogenes triggered mutagenesis at the PpAPT target in about 2% of the regenerated plants. Mainly, deletions were observed, most of them resulting from alternative end‐joining (alt‐EJ)‐driven repair. We further demonstrate that, in the presence of a donor DNA sharing sequence homology with the PpAPT gene, most transgene integration events occur by homology‐driven repair (HDR) at the target locus but also that Cas9‐induced double‐strand breaks are repaired with almost equal frequencies by mutagenic illegitimate recombination. Finally, we establish that a significant fraction of HDR‐mediated gene targeting events (30%) is still possible in the absence of PpRAD51 protein, indicating that CRISPR‐induced HDR is only partially mediated by the classical homologous recombination pathway.     

PsHint1, associated with the G‐protein α subunit PsGPA1, is required for the chemotaxis and pathogenicity of Phytophthora sojae

Zoospore chemotaxis to soybean isoflavones is essential in the early stages of infection by the oomycete pathogen Phytophthora sojae. Previously, we have identified a G‐protein α subunit encoded by PsGPA1 which regulates the chemotaxis and pathogenicity of P. sojae. In the present study, we used affinity purification to identify PsGPA1‐interacting proteins, including PsHint1, a histidine triad (HIT) domain‐containing protein orthologous to human HIT nucleotide‐binding protein 1 (HINT1). PsHint1 interacted with both the guanosine triphosphate (GTP)‐ and guanosine diphosphate (GDP)‐bound forms of PsGPA1. An analysis of the gene‐silenced transformants revealed that PsHint1 was involved in the chemotropic response of zoospores to the isoflavone daidzein. During interaction with a susceptible soybean cultivar, PsHint1‐silenced transformants displayed significantly reduced infectious hyphal extension and caused a strong cell death in plants. In addition, the transformants displayed defective cyst germination, forming abnormal germ tubes that were highly branched and exhibited apical swelling. These results suggest that PsHint1 not only regulates chemotaxis by interacting with PsGPA1, but also participates in a Gα‐independent pathway involved in the pathogenicity of P. sojae.     

Highly efficient homology‐directed repair using CRISPR/Cpf1‐geminiviral replicon in tomato

Genome editing via the homology‐directed repair (HDR) pathway in somatic plant cells is very inefficient compared with error‐prone repair by nonhomologous end joining (NHEJ). Here, we increased HDR‐based genome editing efficiency approximately threefold compared with a Cas9‐based single‐replicon system via the use of de novo multi‐replicon systems equipped with CRISPR/LbCpf1 in tomato and obtained replicon‐free but stable HDR alleles. The efficiency of CRISPR/LbCpf1‐based HDR was significantly modulated by physical culture conditions such as temperature and light. Ten days of incubation at 31 °C under a light/dark cycle after Agrobacterium‐mediated transformation resulted in the best performance among the tested conditions. Furthermore, we developed our single‐replicon system into a multi‐replicon system that effectively increased HDR efficiency. Although this approach is still challenging, we showed the feasibility of HDR‐based genome editing of a salt‐tolerant SlHKT1;2 allele without genomic integration of antibiotic markers or any phenotypic selection. Self‐pollinated offspring plants carrying the HKT1;2 HDR allele showed stable inheritance and germination tolerance in the presence of 100 mm NaCl. Our work may pave the way for transgene‐free editing of alleles of interest in asexually and sexually reproducing plants.