Gene coexpression network analysis combined with metabonomics reveals the resistance responses to powdery mildew in Tibetan hulless barley

Powdery mildew is a fungal disease that represents a ubiquitous threat to crop plants. Transcriptomic and metabolomic analyses were used to identify molecular and physiological changes in Tibetan hulless barley in response to powdery mildew. There were 3418 genes and 405 metabolites differentially expressed between the complete resistance cultivar G7 and the sensitive cultivar Z13. Weighted gene coexpression network analysis was carried out, and the differentially expressed genes were enriched in five and four major network modules in G7 and Z13, respectively. Further analyses showed that phytohormones, photosynthesis, phenylpropanoid biosynthesis, and flavonoid biosynthesis pathways were altered during Qingke-Blumeria graminis (DC.) f.sp. hordei (Bgh) interaction. Comparative analyses showed a correspondence between gene expression and metabolite profiles, and the activated defenses resulted in changes of metabolites involved in plant defense response, such as phytohormones, lipids, flavone and flavonoids, phenolamides, and phenylpropanoids. This study enabled the identification of Bgh responsive genes and provided new insights into the dynamic physiological changes that occur in Qingke during response to powdery mildew. These findings greatly improve our understanding of the mechanisms of induced defense response in Qingke and will provide new clues for the development of resistant Tibetan hulless barley varieties.