Discovery of leaf region and time point related modules and genes in maize (Zea mays L.) leaves by Weighted Gene Co-expression Network analysis (WGCNA) of gene expression profiles of carbon metabolism

Abstract

Maize (Zea mays L.) yield depends not only on the conversion and accumulation of carbohydrates in kernels, but also on the supply of carbohydrates by leaves. However, the carbon metabolism process in leaves can vary across different leaf regions and during the day and night. Hence, we used Weighted Gene Co-expression Network analysis (WGCNA) with the gene expression profiles of carbon metabolism to identify the modules and genes that may associate with particular regions in a leaf and time of day. There were a total of 45 samples of maize leaves that were taken from three different regions of a growing maize leaf at five time points. Robust Multi-array Average analysis was used to pre-process the raw data of GSE85963 (accession number), and quality control of data was based on Pearson correlation coefficients. We obtained eight co-expression network modules. The modules with the highest significance of association with LeafRegion and TimePoint were selected. Functional enrichment and gene-gene interaction analyses were conducted to acquire the hub genes and pathways in these significant modules. These results can support the findings of similar studies by providing evidence of potential module genes and enriched pathways associated with leaf development in maize.