Computational analysis of microarray data of Arabidopsis thaliana challenged with Alternaria brassicicola for identification of key genes in Brassica

J Genet Eng Biotechnol. 2020 Jul 1;18(1):17. doi: 10.1186/s43141-020-00032-y.

Abstract

Background: Alternaria blight, a recalcitrant disease caused by Alternaria brassicae and Alternaria brassicicola, has been recognized for significant losses of oilseed crops especially rapeseed-mustard throughout the world. Till date, no resistance source is available against the disease; hence, plant breeding methods cannot be used to develop disease-resistant varieties. Therefore, in the present study, efforts have been made to identify resistance and defense-related genes as well as key components of JA-SA-ET-mediated pathway involved in resistance against Alternaria brasscicola through computational analysis of microarray data and network biology approach. Microarray profiling data from wild type and mutant Arabidopsis plants challenged with Alternaria brassicicola along with control plant were obtained from the Gene Expression Omnibus (GEO) database. The data analysis, including DEGs extraction, functional enrichment, annotation, and network analysis, was used to identify genes associated with disease resistance and defense response.

Results: A total of 2854 genes were differentially expressed in WT9C9; among them, 1327 genes were upregulated and 1527 genes were downregulated. A total of 1159 genes were differentially expressed in JAM9C9; among them, 809 were upregulated and 350 were downregulated. A total of 2516 genes were differentially expressed in SAM9C9; among them, 1355 were upregulated and 1161 were downregulated. A total of 1567 genes were differentially expressed in ETM9C9; among them, 917 were upregulated and 650 were downregulated. Besides, a total of 2965 genes were differentially expressed in contrast WT24C24; among them, 1510 genes were upregulated and 1455 genes were downregulated. A total of 4598 genes were differentially expressed in JAM24C24; among them, 2201 were upregulated and 2397 were downregulated. A total of 3803 genes were differentially expressed in SAM24C24; among them, 1819 were upregulated and 1984 were downregulated. A total of 4164 genes were differentially expressed in ETM24C24; among them, 1895 were upregulated and 2269 were downregulated. The upregulated genes of Arabidopsis thaliana were mapped and annotated with CDS sequences of Brassica rapa obtained from PlantGDB database. Additionally, PPI network of these genes were constructed to investigate the key components of hormone-mediated pathway involved in resistance during pathogenesis.

Conclusion: The obtained information from present study can be used to engineer resistance to Alternaria blight caused by Alternaria brasscicola through molecular breeding or genetic manipulation-based approaches for improving Brassica oilseed productivity.

Keywords: Brassica; GO analysis; Hubs; Microarray data analysis; Network analysis; Oilseeds.