Identification of key genes in bovine mammary epithelial cells challenged with Escherichia coli and Staphylococcus aureus by integrated bioinformatics analysis

ABSTRACT This study integrated four microarray datasets by Robust Rank Aggregation (RRA) method to identify the differentially expressed genes (DEG) in bovine mammary epithelial (BME) cells in response to Escherichia coli and Staphylococcus aureus infection. Furthermore, the GO function and KEGG pathway enrichment analysis of the integrated DEG were performed. Finally, the protein-protein interaction (PPI) network was constructed. A total of 72 integrated DEG were identified from the four datasets. The most significantly enriched terms within the integrated DEG were mainly involved in the immune response. The PPI network of DEG was constructed with 53 nodes. Seventeen genes, which constitute a significant module, were identified as hub genes. Among them, CD40, CXCL6, and NFKBIZ were further screened as the key genes and have the potential to become biomarkers of E. coli and S. aureus mastitis, considering the specificity of biomarkers for diseases. The identified key genes and pathways in this study can assist in the search for biomarkers for mastitis diagnosis and disease resistance breeding.