Gene regulatory network inference in soybean upon infection by Phytophthora sojae

Phytophthora sojae is a soil-borne oomycete and the causal agent of Phytophthora root and stem rot (PRR) in soybean (Glycine max [L.] Merrill). Yield losses attributed to P. sojae are devastating in disease-conducive environments, with global estimates surpassing 1.1 million tonnes annually. Historically, management of PRR has entailed host genetic resistance (both vertical and horizontal) complemented by disease-suppressive cultural practices (e.g., oomicide application). However, the vast expansion of complex and/or diverse P. sojae pathotypes necessitates developing novel technologies to attenuate PRR in field environments. Therefore, the objective of the present study was to couple high-throughput sequencing data and deep learning to elucidate molecular features in soybean following infection by P. sojae. In doing so, we generated transcriptomes to identify differentially expressed genes (DEGs) during compatible and incompatible interactions with P. sojae and a mock inoculation. The expression data were then used to select two defense-related transcription factors (TFs) belonging to WRKY and RAV families. DNA Affinity Purification and sequencing (DAP-seq) data were obtained for each TF, providing putative DNA binding sites in the soybean genome. These bound sites were used to train Deep Neural Networks with convolutional and recurrent layers to predict new target sites of WRKY and RAV family members in the DEG set. Moreover, we leveraged publicly available Arabidopsis (Arabidopsis thaliana) DAP-seq data for five TF families enriched in our transcriptome analysis to train similar models. These Arabidopsis data-based models were used for cross-species TF binding site prediction on soybean. Finally, we created a gene regulatory network depicting TF-target gene interactions that orchestrate an immune response against P. sojae. Information herein provides novel insight into molecular plant-pathogen interaction and may prove useful in developing soybean cultivars with more durable resistance to P. sojae. Author SummaryGlobal food security is threatened continually by plant pathogens. One approach to circumvent these disease-causing agents entails understanding how hosts balance primary growth and defense upon pathogen perception. Molecular signatures of perception-rendered defense may be leveraged subsequently to develop resistant/tolerant crop plants. Additionally, evidence suggests that the plant immune system is characterized by tuning primary and secondary metabolic activity via transcription factor-mediated transcriptional reprogramming. Therefore, we investigated transcription factor-target gene interactions in soybean upon infection by compatible and incompatible races of Phytophthora sojae. Through transcriptome analysis, we found that the interactions elicited vast, overlapping transcriptional responses and identified overrepresented, defense-related transcription factor families. We then generated/acquired DNA-protein interactome data for the most represented transcription factor families in the transcriptome analysis and trained deep learning-based models to predict novel transcription factor targets. Transcription factor/target gene metrics were used to construct a gene regulatory network with prioritized components. We identified hub transcription factors belonging to WRKY and ERF families, the majority of which function in response to various biotic and abiotic stressors. These findings propose novel regulators in the soybean defense response to Phytophthora sojae and provide an avenue for the investigation of transcription factor-target gene interactions in plants.