Genome-wide gene expression profiling allows for identification of genes involved in the defense response of the host against pathogens. As presented here, transcriptomic analysis and bioinformatics tools were applied in order to identify genes expressed in the mammary gland parenchyma of cows naturally infected with coagulase-positive and coagulase-negative Staphylococci.In cows infected with coagulase-positive Staphylococci, being in 1st or 2nd lactation, 1700 differentially expressed genes (DEGs) were identified. However, examination of the 3rd or 4th lactations revealed 2200 DEGs. Gene ontology functional classification showed the molecular functions of the DEGs overrepresented the activity of cytokines, chemokines, and their receptors. In cows infected with coagulase-negative Staphylococci, in the 1st or 2nd lactations 418 DEGs, while in the 3rd or 4th lactations, 1200 DEGs were identified that involved in molecular functions such as protein, calcium ion and lipid binding, chemokine activity, and protein homodimerization. Gene network analysis showed DEGs associated with inflammation, cell migration, and immune response to infection, development of cells and tissues, and humoral responses to infections caused by both types of Staphylococci.A coagulase-positive Staphylococci infection caused a markedly stronger host response than that of coagulase-negative, resulting in vastly increased DEGs. A significant increase in the expression of the FOS, TNF, and genes encoding the major histocompatibility complex proteins (MHC) was observed. It suggests these genes play a key role in the synchronization of the immune response of the cow's parenchyma against mastitis-causing bacteria. Moreover, the following genes that belong to several physiological pathways (KEGG pathways) were selected for further studies as candidate genes of mammary gland immune response for use in Marker Assisted Selection (MAS): chemokine signaling pathway (CCL2, CXCL5, HCK, CCR1), cell adhesion molecules (CAMs) pathway (BOLA-DQA2, BOLA-DQA1, F11R, ITGAL, CD86), antigen processing and presentation pathway (CD8A, PDIA3, LGMN, IFI30, HSPA1A), and NOD-like receptor signaling pathway (TNF, IL8, IL18, NFKBIA).