Construction and analysis of gene co-expression network in the pathogenic fungus Ustilago maydis

Biological systems respond to environmental disturbances and a wide range of compounds through complex gene interaction networks. The enormous growth of experimental information obtained using large-scale genomic techniques such as microarrays and RNA sequencing led to the construction of a wide variety of gene co-expression networks in recent years. These networks allow the discovery of clusters of co-expressed genes that potentially work in the same process linking them to biological processes often of interest to industrial, medicinal, and academic research.

In this study, we built the gene co-expression network of Ustilago maydis from the gene expression data of 168 samples belonging to 19 series, which correspond to the GPL3681 platform deposited in the NCBI using WGCNA software. This network was analyzed to identify clusters of co-expressed genes, gene hubs and Gene Ontology terms. Additionally, we identified relevant modules through a hypergeometric approach based on a predicted set of transcription factors and virulence genes.

We identified 13 modules in the gene co-expression network of U. maydis. The TFs enriched in the modules of interest belong to the superfamilies of Nucleic acid-binding proteins, Winged helix DNA-binding, and Zn2/Cys6 DNA-binding. On the other hand, the modules enriched with virulence genes were classified into diseases related to corn smut, Invasive candidiasis, among others. Finally, a large number of hypothetical, a large number of hypothetical genes were identified as highly co-expressed with virulence genes, making them possible experimental targets.