Fig. 7

Weighted gene co-expression network analysis (WGCNA) in lymphoid tissues for Control, M. bovis (MB), and Dual groups. (A) Module-treatment relationship graph where each row represents the module eigen value and each column represents infection status. The cells within the matrix show the correlation coefficient and p-value. Modules were found in mesenteric lymph node (MLN), tracheal-bronchial lymph node (TBLN), thymus, and spleen. (B) Correlation between module membership (MM) and gene significance (GS) of Module III genes, where MM represents the correlation between gene expression and the module eigen values and GS represents the correlation between gene expression and co-infection status. (C) Biological processes, pathways, and molecular functions enriched with co-expressed genes in Module III. The y-axis indicates the gene ontology (GO) term and the x-axis indicates the -log10 p-value. The size of the dot indicates the number of co-expressed genes enriched in the GO term. Enriched GO terms with a p-value < 0.05 (dashed line) were considered significant. IFI6 = interferon alpha inducible protein 6; BST2 = bone marrow stromal cell antigen 2; HERC6 = HECT and RLD Domain Containing E3 Ubiquitin Protein Ligase Family Member 6; ISG15 = interferon-stimulated gene 15; IFITM1/3 = interferon-induced transmembrane proteins 1/3