The lemon industry in the Mediterranean basin is strongly threatened by "mal secco" disease (MSD) caused by the fungus Plenodomus tracheiphlilus. Leaf pretreatments with Pseudomonas mediterranea 3C have been proposed as innovative tools for eco-sustainable interventions aimed at controlling the disease. In this study, by exploiting the results of previously performed RNAseq analysis, WCGNA was conducted among gene expression patterns in both inoculated (Pt) and pretreated and fungus-inoculated lemon plants (Citrus limon L.) (3CPt), and two indicators of fungal infection, i.e., the amount of fungus DNA measured in planta and the disease index (DI). The aims of this work were (a) to identify gene modules significantly associated with those traits, (b) to construct co-expression networks related to mal secco disease; (c) to define the effect and action mechanisms of P. mediterranea by comparing the networks. The results led to the identification of nine hub genes in the networks, with three of them belonging to receptor-like kinases (RLK), such as HERK1, CLAVATA1 and LRR, which play crucial roles in plant-pathogen interaction. Moreover, the comparison between networks indicated that the expression of those receptors is not induced in the presence of P. mediterranea, suggesting how powerful WCGNA is in discovering crucial genes that must undergo further investigation and be eventually knocked out.
Keywords: Plenodomus tracheiphilus; Pseudomonas mediterranea; RNAseq; WGCNA; biocontrol agent; lemon; mal secco; transcriptome.