Identification of key transcription factors and their functional role involved in Salmonella typhimurium infection in chicken using integrated transcriptome analysis and bioinformatics approach
BMC Genomics volume 24, Article number: 214 (2023)
Salmonella enterica serovar typhimurium is the cause of significant morbidity and mortality worldwide that causes economic losses to poultry and is able to cause infection in humans. Indigenous chicken breeds are a potential source of animal protein and have the added advantage of being disease resistant. An indigenous chicken, Kashmir favorella and commercial broiler were selected for understanding the mechanism of disease resistance. Following infection in Kashmir favorella, three differentially expressed genes Nuclear Factor Kappa B (NF-κB1), Forkhead Box Protein O3 (FOXO3) and Paired box 5 (Pax5) were identified. FOXO3, a transcriptional activator, is the potential marker of host resistance in Salmonella infection. NF-κB1 is an inducible transcription factor which lays the foundation for studying gene network of the innate immune response of Salmonella infection in chicken. Pax5 is essential for differentiation of pre-B cells into mature B cell. The real time PCR analysis showed that in response to Salmonella Typhimurium infection a remarkable increase of NF-κB1 (P˂0.01), FOXO3 (P˂0.01) gene expression in liver and Pax5 (P˂0.01) gene expression in spleen of Kashmir favorella was observed. The protein–protein interaction (PPI) and protein-TF interaction network by STRINGDB analysis suggests that FOXO3 is a hub gene in the network and is closely related to Salmonella infection along with NF-κB1. All the three differentially expressed genes (NF-κB1, FOXO3 and PaX5) showed their influence on 12 interacting proteins and 16 TFs, where cyclic adenosine monophosphate Response Element Binding protein (CREBBP), erythroblast transformation-specific (ETSI), Tumour-protein 53(TP53I), IKKBK, lymphoid enhancer-binding factor-1 (LEF1), and interferon regulatory factor-4 (IRF4) play role in immune responses. This study shall pave the way for newer strategies for treatment and prevention of Salmonella infection and may help in increasing the innate disease resistance.
Salmonella is a gram-negative, flagellated, rod shaped, facultative anaerobe, and a food borne pathogen which is commonly associated with contaminated food and water. Salmonella enterica accounts for approximately 99% of salmonella infections . It has been detected in cattle, pigs, poultry and food products . S. enterica serovar typhimurium is zoonotic and the most common cause of human food poisoning (11%) which iscaused by poultry products infected by S. enterica serovar Typhimurium also written as S. enterica serovar typhimurium or simply Salmonella typhimurium. . Salmonella typhimurium can be transferred to humans via undercooked or raw infected food, particularly meat and eggs. Salmonella cause an illness called salmonellosis, and some salmonella types can cause typhoid fever, paratyphoid fever as well as gastroenteritis. The lack of an effective immunoprophylactic strategy warrants the need for strict biosecurity in addition to proper cooking and cleaning of food. A So, analysis on Salmonella typhimurium infection is highly looked-for, from the perspective of both human and animal health.
Transcription factors play an important role in response to bacterial, viral, and fungal infections, and have been reported as differentially expressed genes involved in disease resistance . Over the last few years, genetic selection of birds is considered to be a systematic and continuous way to control Salmonella infections . Selecting chicken that are genetically resistant to Salmonella infection offers an alternative environment-friendly control measure. Jammu and Kashmir state in northern India is home to the native breed of chicken known as the Kashmir favorella. It is raised primarily for the production of meat and eggs, but is also regarded as a significant source of animal protein . This native breed is well adapted to the local climatic conditions, feed, and stress management. In addition to this, this breed also has a high disease resistance . In this study, the differentially expressed transcriptional factors NF-κB1, FOXO3 and Pax5 in liver and spleen were examined. After infecting chicken with Salmonella Typhimurium, pathogen was first found at the intestinal tract followed by liver and spleen . Salmonella enterica can quickly enter the bloodstream, colonize the internal organs including liver, spleen and heart. In this paper identification of immune related genes (NF-κB1, FOXO3 and Pax5) was done. They confer resistance to infection and may have an impact on potential attempts to control salmonellosis in chicken and subsequently, help in preventing human salmonellosis. These findings will provide newer approaches to develop strategies for Salmonella typhimurium prevention and treatment.
Materials and methods
The tissue was collected after approval from Institutional Animal Ethics Committee of Faculty of Veterinary Sciences & Animal Husbandry (SKUAST-Kashmir) ethical standards in animal experimentation (AU/FVSc/PS-57/16021). All the guidelines of the Institutional Animal Ethics Committee of SKUAST Kashmir were strictly followed during all the animal experiments.
Experimental birds and Bacterial strain
Susceptible and resistant one day old chicks were initially screened for the presence of any infection. After the initial screening, one day old chicks from Kashmir favorellaand broiler (Cobb) were procured from Division of Livestock Production and Management, SKUAST-Kashmir,India. The experimental birds were nurtured under standard hygienic conditions at the animal house facility at Faculty of Veterinary Sciences and Animal Husbandry, SKUAST-Kashmir, India. The birds were kept under surveillance daily and had ad-libitum access to feed and water. Salmonella Typhimurium strain (ATCC 14,028) from the American Type Culture Collection was procured.
A pilot study in Kashmir favorella and broiler breeds was carried out where both the breeds were infected at day 4 of age and at the outset were evaluated for disease resistance up to 10 days post infection. The breeds showed optimal differences in clinical symptoms and bacterial loads at day 5 post infection. The birds were divided into two groups i.e., resistant and susceptible groups based on the slight or severity of differences in clinical symptoms, pathological manifestations and bacterial loads. Kashmir favorella was found to be resistant while the broilers were found to be susceptible to Salmonella infection based on visible clinical signs and bacterial load. Both the breeds were divided into control and infected groups and were scrutinized for the existence of Salmonella before challenge. Only Salmonella negative birds were taken for the study. After 3 days of habituation, the infected group was challenged orally with 2 × 108 CFU/ml of S. typhimurium (14,028) strain. Chicks in the control group were artificially challenged by 1 ml of nutrient broth orally.
The experimental birds were maintained under standard hygienic, temperature and pressure conditions at the animal house facility centre at Faculty of Veterinary Sciences and Animal Husbandry (FVSc & AH), SKUAST-Kashmir, India.. The birds were monitored daily and had ad libitum access to feed and water. Salmonella typhimurium strain (ATCC 14,028) was obtained from the American Type Culture Collection. The culture was revived in tetrathionate broth (TTB) after incubating at 37 °C for 18 h. The overnight growth from TTB was streaked on MacConkey and Brilliant green agar (BGA) plates, incubated at 37 °C for 24 h and finally examined for typical Salmonella colonies.
Two experiments were carried out in the current study. A pilot study was conducted to assess the Salmonella resistance and susceptibility between the two chicken breeds (Kashmir favorella and commercial broiler). On the basis of the pilot study, the actual experiment was carried out.
Assessment of Salmonella resistance and susceptibility
For the initial screening of the poultry birds for disease resistance and susceptibility, a pilot study was carried in which one-day old chicks were procured from two poultry breeds i.e. Kashmir favorella (n = 100) and broiler (Cobb) (n = 100). The birds were orally challenged with 2 × 108 CFU/ml of Salmonella typhimurium (14,028) strain at day 4 of age and were initially assessed for disease resistance up to 10 days post infection. For confirmation of infection, faecal swabs were taken after 12 h post infection from all the birds and were incubated at 42 °C for 18 h in tetrathionate broth (TTB). Clinical symptoms and gross pathology were also observed. The infection was further confirmed by amplification of Salmonella specific 16S rRNA and Salmonella typhimurium specific primers.
The birds were divided into two groups i.e. resistant and susceptible groups based on the differences in clinical symptoms, pathological manifestations and bacterial loads. The chicks that exhibited severe clinical symptoms (progressive weakness, anorexia, diarrhoea and lowering of the head), severe liver pathology and higher bacterial loads in faecal swabs compared with others were identified as challenged susceptible group. Chicks with slight clinical and pathological manifestations and lower bacterial loads were identified as challengedresistant birds. Based on the above clinical symptoms and bacterial loads, Kashmir favorella was found to be resistant and broilers were found to be susceptible to Salmonella infection. Day 5 post infection the infected chicken showed optimal differences in clinical symptoms and bacterial loads.
After the initial screening of chicks (resistant and susceptible), one-day-old chicks from both the breeds i.e. Kashmir favorella (n = 50) and broiler (n = 50) procured from the Division of Livestock Production and management, SKUAST-Kashmir-India, were maintained under standard hygienic, temperature and pressure conditions at the animal house facility centre. The chicks were monitored daily and had ad libitum access to feed and water. The chicks from both breeds were divided into control (n = 25) and infected groups (n = 25) and were screened for the presence of Salmonella before challenge. All the chicks were found to be Salmonella negative. After 3 days of acclimatization, the infected group was challenged orally with 2 × 108 CFU/ml of Salmonella Typhimurium (14,028) strain. Chicks in the control group were mock challenged and were given 1 ml of nutrient broth orally.
The birds were euthanized with CO2 and tissue samples (liver and spleen) were taken from twelve infected chicks (resistant 6 and susceptible 6) at day 5 post infection and processed for pathological studies, bacterial enumeration and RNA isolation. Six respective controls were taken from each group. The selection of the chicks for RNA sequencing was based on the highest clinical score, bacterial counts, gross pathology, and lesion scores in their respective groups. The clinical scores were recorded twice daily following points-based scoring system .
For global identification of genes and transcription factorss interacting with NF-κB1, FOXO3 and PAX5, we utilized RNA-Sequencing data from our previously published dataset NCBI GEO ID: GSE168060 . The standard protocols/pipeline was utilized for RNA extraction, cDNA preparation and RNA-Sequencing . Total RNA was extracted using Trizol method (Ambion, USA). RNA quality was assessed by spectrophotometer (ThermoFisher, USA) and integrity was assessed by Bioanalyzer (Agilent, USA). Only those RNA samples whose RIN values ≥ 8 were used to construct libraries. Illumina TruSeq stranded mRNA sample preparation kit was used to generate cDNA libraries following manufacturer’s protocol. 4 µg/sample of total RNA was used. Poly-T attached magnetic beads were used for the purification of poly-A containing mRNA molecules. Once the purification was done, the RNA was fragmented into smaller pieces under high temperature using divalent cations. The RNA fragments thus generated, were used for the preparation of first strand cDNA using enzyme reverse transcriptase and random primers (Illumina, USA). Using DNA polymerase 1 and RNase H, the second strand cDNA was synthesized. DNA fragments (at 3ʹ ends) were adenylated and hybridized by ligating Illumina paired-end adapter and index. Using the Illumina PCR primer cocktail, the cDNA fragments (150 bp) were generated and selectively enriched to prepare the final sequencing paired end cDNA library . The High Throughput Model flow cell on an Illumina HiSeq 2500 platform were used to pool the libraries in equimolar amounts and paired end sequenced by SciGenom, Cochin, Kerala-India.
Total RNA Isolation and cDNA synthesis
Total RNA was extracted from the tissue samples of liver and spleen by Trizol™ method (Invitrogen, USA) as per manufacturer's protocol. The quantity and quality of isolated RNA was checked with UV–Visible spectrophotometer. RNA samples were run on 1% agarose gel to check the quality of RNA. DNase treatment using DNase1 kit (Sigma, USA) was given to rule out genomic DNA contamination. cDNA synthesis was done with equal concentration of RNA (1.5 μg/μl) in all the samples using Thermo Scientific RevertAid First Strand cDNA Synthesis Kit™ (Lithuania) using oligodT primers following manufacturer protocol.
Quantitative real-time PCR
Real-time PCR was performed to verify the expression of three key transcription factors (NF-κB 1, FOXO3 and PAX5) which are involved in different immunological pathways. Following the manufacturer’s protocol, total RNA was used to synthesize cDNA using revert Aid cDNA synthesis kit (Thermo Scientific, USA). The primer details are enlisted in (Table 1). All the qPCR reactions were run in replicates using Light Cycler II (Roche). β-actin and GAPDH were used as internal controls. The reaction mixture (20 µl) consisted of 8.4 µl nuclease free water, 0.3 µl of 10 µM each forward and reverse primers, SYBR green master mix (Roche) 10 µl and 70 µg/µl of cDNA. The reactions were amplified by 40 cycles of denaturation at 95 °C for 5 min, annealing at variable temperatures (specific for each gene) for 15 s and extension at 72 °C for 15 s. For the confirmation of product specificity, melting curve analysis was done. The relative mRNA expression was determined by using 2−∆∆CT method .
Network and enrichment analysis
STRINGDB was utilized to identify the possible protein–protein interaction (PPI) interaction network. Two topographical features “Betweenness between nodes” and “Degree” were used to identify key interacting partners in the PPI network. Higher the quantitative value of a node, the more important is the node in the network. The TF were subjected to functional using g: Profiler [a] web server using Gallus gallus (chicken) as reference organism.
Once the chicks were inoculated with Salmonella typhimurium, typical clinical signs of progressive weakness, anorexia, diarrhoea, exhaustion, dullness, ruffled feathers, lowering of head, dropping of wings, inappetence, and reluctance to move with both eyes closed were observed (Fig. 1). For further confirmation of infection, amplification of Salmonella specific 16S rRNA and Salmonella typhimurium specific primers was done using PCR.
The resistant group i.e. Kashmiri favorella displayed mild clinical symptoms with low bacterial load while as the birds in the susceptible group i.e. broiler displayed severe clinical signs with high bacterial load. The greatest differences in clinical symptoms, pathology, and bacterial load were observed on day 5 after inoculation. The resistant birds had minor liver lesions compared to the susceptible birds, which had major necrotic lesions. (Fig. 2 and Fig. 3).
Gene expression study
Transcriptome data was cleaned using pipeline described . In this study we performed a multi-group comparison between liver samples (Boiler Liver Infected (A) Vs Boiler Liver Control (B) Vs Kashmiri Liver Control (E)) and spleen samples (Boiler Spleen Infected (C) Vs Boiler Spleen Control (D) Vs Kashmiri Spleen Control (G) Vs Kashmiri Spleen Infected (F) to identify more genes and TFs associated with infection resistance. In accordance with qPCR results, RNA-Sequencing data shows similar trend in the expression of all three TFs in spleen and liver Samples. FOXO3 and NF-κB1 were expressed more in liver and PAX5 in spleen. Fold increase in liver of Kashmir favorella liver and spleen is given in Table 2. The PCA plot shows clear separation between samples (Supplementary Fig. 1).
2-Dimensional hierarchical clustering of top DEGs in spleen samples (multi-group comparison) shows the key role of PAX5 in regulating infection resistance in Kashmiri favorella (Fig. 4). The protein–protein and protein-TFs interaction analysis (Fig. 5) suggests that FOXO3 (Degree = 23 after thresholds) was hub gene in the network and closely related in Salmonella infection with NF-κB1, which is in accordance with the RNA-Seq and q-PCR results.
To gain insight into the biological processes and pathways mediated by PAX5, FOXO3 and NF-κB1Gene ontology (GO) and KEGG enrichment analysis was performed. Pathway analyses suggest that NF-κB1 and FOXO3 were involved in Gallus Gallus Cellular senescence signalling pathway (FDR adjusted P-value = 9.107 × 10–3). The GO analysis suggest that all three TFs were involved in molecular function; DNA-binding transcription activator activity (padj = 2.378 × 10–3) and biological process negative regulation of transcription (padj = 3.064 × 10–2). The transcription factors of interest NF-κB1, FOXO3 and PAX5 were significantly upregulated in Kashmir favorella. All the three showed their influence on 12 interacting proteins (CREBBP, AR, ETS1, TP53, IKBKB, UBC, CUL4B, SKP2, ESR1, CHUK, BTRC, TSC22D3) and 16 TFs (LEF1, NFRKB, IRF4, REST, GFI1B, GAGAD1, BCL11B, TEADE, EGR1, HMG20, BCTNNB1, MXD3, ZNF263, ZNF8, ZNF83, ZNF24), where CREBBP, ETSI, TP53I, IKKBK, LEF1, IRF4 play role in immune responses as shown in Fig. 5. Furthermore, the differentially expressed genes that were mainly related to immune function, such as TNFRSF11B, SOCS1, IL-18, IL6, IL1β, IL10RB, IKBKB, FOS, and EGR1 were also upregulated in Kashmir favorella.
Real-time quantitative PCR validation of the RNA-seq expression results
The differentially expressed transcription factors NF-κB1, FOXO3 and PAX5 in Kashmir favorella were significantly upregulated when compared to commercial broiler where the expression of NF-κB1 was 12.04-fold with P˂ 0.01, FOXO3 was 74.54-fold with P˂ 0.01 in liver and PAX5 was 61.39-fold with P˂0.01 in spleen as given in Fig. 1 and Table 2. Similar trend in the expression of NF-κB1 and PaX5 were identified in RNA seq data. The comparative RNA sequencing expression and qPCR expression is shown (Fig. 6).
The genetic mechanism of disease resistance in poultry against Salmonella infection have been described by several studies. Our study elucidates the genetic basis of disease resistance in Kashmir favorella and commercial Broiler against Salmonella typhimurium based on the severity of clinical symptoms, bacterial load post infection and the pathological manifestations. Following Salmonella typhimurium challenge in chicken liver and spleen, we identified certain transcription factors that played an important role in resistance and susceptibility of two different chicken breeds (Commercial broiler and Kashmir favorella). Nuclear Factor kappa B subunit (NF-κB1), Fork head box transcription factor 3 (FOXO3) and Paired box 5 (PaX5), superoxide dismutases (SOD1), Immunoglobulin Superfamily Containing Leucine Rich Repeat (ISLR), Pyruvate Dehydrogenase E1 Subunit Beta (PDHB), Zinc finger protein (GLI2) transcription factors were identified. FOXO3 is involved in a variety of cellular functions including cell-cycle control, DNA damage repair, cell growth and differentiation, response to oxidative stress, and apoptosis . Chicken FOXO3 reportedly shares more than 85.4% identity with humans . The FOXO3 gene controls lymphoid homeostasis in the host immune system [20, 21]. FOXO3 serves as a marker of host resistance to Salmonella infection. In previous studies, JAK-STAT signalling pathway, MAPK-related pathways and cytokine-cytokine receptor interaction were activated significantly in resistant chickens . In our study, mRNA levels of FOXO3 were quantified in liver and spleen of Kashmir favorella. FOXO3 is highly upregulated in liver of Kashmir favorella. It was found that MAPK pathway, cytokine-cytokine and JAK-STAT pathways are upregulated in Kashmir favorella. TNFRSF11B is one of the members of tumor necrosis factor receptor (TNFR) molecular sub-family was also upregulated in Kashmir favorella and has role in RANK activation . The ligation of RANK by RANKL plays an important role in the immune response regulating the interactions between T cells and dendritic cells. The JAK-STAT pathway is also activated involving genes IL10RA, IL10RB, SOCS1 positively regulated. The JAK-STAT pathway ensures T and B-cell development . In the present study, various genes involved in the JAK-STAT signalling pathway had increased expression including IL10RB that activate monocytes and neutrophils and SOCS1 that regulates M1 macrophage activation. FOS and IRG1 have role in cell proliferation and differentiation and promote endotoxin tolerance respectively. Also, IL18 has role in IFN-γ production, promote Th1 immunity host defense inflammation (innate and acquired immunity) is highly upregulated in Kashmir favorella.
NF-κB1 was more expressed in liver of Kashmir favorella than commercial broiler. It is a primary regulator of the innate immune response to bacterial invasion. It is activated by MyD88-dependent TLR pathway signalling . It has five related proteins (subunits) that can bind to DNA: p50 (a 50 kDa protein also known as NF-KB1), p52 (also known as NF-B2), p65 (also known as RelA), c-Rel, and RelB. . According to our findings, Salmonella infection causes inflammatory expression of NF-κB genes. NF-κB is the master regulator of the inflammatory response in the complex inflammatory network, serving as both a host defence against pathogens and a cell protector against apoptosis. This occurs as a result of the activation of proinflammatory genes . It is crucial in the activation of chemokines/cytokines (IL18) . The NF-κB1 pathway was successfully activated in Kashmir favorella. Our results showed that NF-κB1 and FOXO3 were upregulated in Kashmir favorella play their role in expression of interacting proteins and transcription factors as given in Fig. 6. FOXO3 and NF-κB1 were expressed more in liver, as liver sinusoids monocytes are present in bulk whereas in spleen only red pulp contains monocytes [19, 28]. FOXO3 binds NF-KBRelA, which provides cytoplasmic stability while also inhibiting proinflammatory signalling pathways. NF-κB1 was able to shuttle to the nucleus and turn on proinflammatory target genes (IL-6, IL-12, and IL-15) after FOXO3 expression was silenced . The signalling adapter molecule IKBKB was significantly upregulated in Kashmir favorella after Salmonella infection that showed its influence on FOXO3 expression as well. Furthermore, EGR1 and IRF4 were also upregulated. EGR1 plays role in inflammatory responses whereas IRF4 in development of NK cells, B and T cells. Bcl11b is essential for development of CD4+ and CD8+ T cells [30,31,32].
The expression of PAX5 (BSAP) was remarkably more in spleen of Kashmir favorella. PAX5 expressed more in spleen as it was continuously required for B-cell function . It has 94% similarity with mice PAX5 and is solely expressed from the pro-B to mature B cell stages and is downregulated in the course of terminal differentiation into plasma cells . It is a vital transcription factor in the process of B cell development and intricated in the activation of the chromatin of essential genes involved in B-cell signalling, adhesion, migration, inflammatory and immune responses [35,36,37]. It plays a principal role in the upregulation of B lymphocyte-induced maturation protein 1 (Blimp-1) and X-box binding protein 1 (XBP-1) [38, 39] obligatory for plasma cell differentiation [3, 40]. Blimp-1 is enough to shut down the B cell gene expression and initiate the development of plasma cells. Apart from this it also plays a significant role in activation of CD19, CD79A and repression of FMS, Notch1, CSF1 receptor [16, 41]. Notch 2 is expressed in liver of Kashmir favorella. PAX5 either act as activator or repressor.
Both NF-κB1 and PAX5 play role in the expression of CREBBP, AR, LEF1 and TP53. Lymphoid Enhancer Binding Factor-1 (LEF1) is a Wnt/-catenin signalling pathway mediator that interacts with -catenin to regulate Wnt target gene expression [42, 43]. When the Wnt/-catenin pathway is activated, catenin enters the nucleus and interacts with transcription factors T-cell and lymphoid enhancer factors (TCF/LEF1) to activate and/or repress transcription of specific target genes . LEF1 in chickens is expressed in pre-B and T-cells, T-helper 1 cell differentiation and B-cell proliferation. Cyclic adenosine monophosphate Response Element Binding Protein Binding Protein (CREBBP or CREB) is present in liver of Kashmir favorella. The CREB promotes activation and proliferation of T and B cells and differentially regulates Th1, Th2, and Th17 responses . In Kashmir favorella IL10, TNF are highly upregulated. Tp53 also regulates innate immune and adaptive immune responses. FOXO3 and PaX5 expression also show its effect on ETS1. ETS1 deficiency has a significant impact on B cells, T cells, and NK cells, all of which express ETS1 at high levels under normal physiological conditions . PAX5 and FOXO3 have been shown to interact with ETS1. It is manifested in spleen. The ETS1 transcription factor belongs to the ETS gene family and has been found to be highly conserved throughout evolution (600–700 million years ago). The ETS gene family includes 28 genes in humans (441 amino acids) and 27 genes in mice (440 amino acids). Mice ETS1 has 95% similarity with that of chicken ETS1. Adult humans and chicken show a similar pattern where high expression of ETS1 is found mainly in lymphoid tissues [47, 48]. ETS1 is expressed in B cells (highly expressed in naïve mature B-cells and memory B-cells), T cells, NK cells, and NK T cells [49, 50] and is inducible in other, non-lymphoid cell types in response to certain stimuli .
SOD1 in presence of copper and zinc catalyse the transfer of superoxide anions into less dangerous free radical (hydrogen peroxide). It decreases ROS formation . It has antibacterial and anti-fungal activity. SOD1 also show its effect on NF-κB and p38 pathways involved in influenza A virus pathogenesis and pro inflammatory responses. SOD1 overexpression significantly attenuates transcription of several pro inflammatory cytokines like CCL2, IL6, TNFα, IL-8 . ISLR is an immunoglobulin superfamily containing leucine rich repeats that stabilizes Wnt signalling and regulates muscle atrophy (skeletal muscles) via IGF1-PI3K/Akt-FOXO signalling pathway  PDHB plays role in metabolic processes TCA cycle, gluconeogenesis acetyl CoA synthesis and canonical pathways . GLI-2 is a zinc finger protein which acts as transcriptional activator in hedgehog (hh) signalling. Expression of GLI-2 has an effect on activation of NF-κB pathway as well as cell apoptosis. Upregulation of GLI-2 decreases the proinflammatory cytokines, interleukin-6, IFN-γ, FASL whereas increases expression IL-10 .
Based on our RNA Seq data, we discovered that Kashmir favorella was more resistant to Salmonella typhimurium infection than commercial broiler. A consistency in gross pathology, histopathology and lesion score with bacterial load in liver and spleen was found.
Transcriptome profiling of the liver and spleen in two chicken breeds/strains i.e., Kashmir favorella and commercial broiler manifest that Kashmir favorella shows a significant improvement in the majority of the genes and signalling pathways involved in protecting the host from bacterial infections. The differentially expressed genes NF-κB1, FOXO3 and PAX5 are involved in defending the host against bacterial infections were significantly enriched in the indigenous breed Kashmir favorella. Also, the genes that are found to be involved in immune response and disease tolerance were differentially expressed in Kashmir favorella. The protein–protein interaction (PPI) and protein -TF interaction networks by STRINGDB analysis suggests that FOXO3 was hub gene in the network and closely related in Salmonella infection with NF-κB1. By using model animals like chicken, these findings will help us better understand the phenomenon of disease resistance and susceptibility to other bacterial and viral infections through pathogen-induced genes. In addition, this research will be useful for developing fresh plans and methods for preventing and treating Salmonella infections. It may also increase innate disease resistance through genetic selection in poultry, which would lower the risk of human Salmonellosis.
Availability of data and materials
The sequencing data is available in NCBI under accession number GSE168060.
Brenner FW, Villar RG, Angulo FJ, Tauxe R, Swaminathan B. Salmonella nomenclature. J Clin Microbiol. 2000;38:2465–7.
Andino A, Hanning I. Salmonella enterica: survival, colonization, and virulence differences among serovars. Scientific World J. 2015;2015:520179.
Rabsch W, Tschape H, Baumler AJ. Non-typhoidal salmonellosis: emerging problems emerging problems. Microbes Infect. 2001;3:237–324.
Dar MA, Urwat U, Ahmad SM, Ahmad R, Kashoo ZA, Dar TA, Heidari. Gene expression and antibody response in chicken against Salmonella Typhimurium challenge. Poultry Sci. 2019;98:2008–13.
Li P, Fan W, Everaert N, Liu R, Li Q, Zheng M, Wen J. Messenger RNA sequencing and pathway analysis provide novel insights into the susceptibility to Salmonella enteritidis infection in chickens. Front Gen. 2018;9:256.
Wani H, Darzi MM, Kamil SA, Wani SA, Munshi ZH, Shakoor A, Shah A. Histological and histochemical studies on the reproductive tract of Kashmir favorella chicken. J Etnomol Zool Stud. 2017;5:2256–62.
Iqbal S, Pampori ZA. Production potential and qualitative traits of indigenous chicken of Kashmir. Livest Res Rural Dev. 2008;20:14.
Lee GM, Jackson GD, Cooper GN. Infection and immune responses in chickens exposed to Salmonella typhimurium. Avian Dis. 1983;27:577–83.
Dulwich KL, Asfor A, Gray A, Giotis ES, Skinner MA, Broadbent AJ. The stronger downregulation of in vitro and in vivo innate antiviral responses by a very virulent strain of infectious bursal disease virus (IBDV), compared to a classical strain, is mediated, in part, by the VP4 protein. Front Cell Infect Microbiol. 2020;10:315.
Dar MA, Ahmad SM, Bhat BA, Dar TA, UlHaq Z, Wani BA, Shabir N, Kashoo ZA, Shah RA, Ganai NA, Heidari M. Comparative RNA-Seq analysis reveals insights in Salmonella disease resistance of chicken; and database development as resource for gene expression in poultry. Genomics. 2022;114(5):110475.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) method. Methods. 2001;25:402–8.
He XM, Fang MX, Zhang ZT, Hu YS, Jia XZ, He DL, Liang SD, Nie QH, Zhang XQ. Characterization of chicken natural resistance-associated macrophage protein encoding genes (Nramp1 and Nramp2) and association with salmonellosis resistance. Genet Mol Res. 2013;12:618-30.
Berndt A, Wilhelm A, Jugert C, Pieper J, Sachse K, Methner U. Chicken cecum immune response to Salmonella enterica serovars of different levels of invasiveness. Infection and Immunity. 2007;75(12):5993-6007.
Chiang HI, Berghman LR, Zhou H. Inhibition of NF-kB 1 (NF-kBp50) by RNA interference in chicken macrophage HD11 cell line challenged with Salmonellaenteritidis. Genet Mol Biol. 2009;32(3):507-15.
Chen B, Guo L, Chen X, El-Senousey HK, Ma M, Jebessa E, Nie Q. Cellular function of chicken FOXO3 and its associations with chicken growth. Poultry Science. 2019;98(10):5109-17.
Nera KP, Lassila O. Pax5–a critical inhibitor of plasma cell fate. Scand J Immunol. 2006;64:190–9.
Bhat B, Yaseen M, Singh A, Ahmad SM, Ganai NA. Identification of potential key genes and pathways associated with the pashmina fiber initiation using RNA-seq and integrated bioinformatics analysis. Sci Rep. 2021;11:1–9.
Farhan M, Wang H, Gaur U, Little PJ, Xu J, Zheng W. FOXO signalling pathways as therapeutic targets in cancer. Int J Biol Sci. 2017;13:815–27.
Chen B, Xu J, He X, Xu H, Li G, Du H, Zhang. A genome-wide mRNA screen and functional analysis reveal FOXO3 as a candidate gene for chicken growth. PloS one. 2015;10:0137087.
Lin L, Hron JD, Peng SL. Regulation of NF-κB, The activation, and autoinflammation by the forkhead transcription factor Foxo3a. Immunity. 2004;21:203–13.
Lu L, Barbi J, Pan F. The regulation of immune tolerance by FOXP3. Nat Rev Immunol. 2017;17:703–17.
Chistiakov DA, Sobenin IA, Orekhov AN. Regulatory T cells in atherosclerosis and strategies to induce the endogenous atheroprotective immune response. Immunol Lett. 2013;151:10–22.
Rawlings JS, Rosler KM, Harrison DA. The JAK/STAT signaling pathway. J Cell Sci. 2004;117:1281–3.
Chiang H I. Genomic approaches to study innate immune response to Salmonella Enteritidis infection in chickens. Texas A&M University.
Sehnert B, Burkhardt H, Dübel S, Voll RE. Cell-type targeted NF-kappaB inhibition for the treatment of inflammatory diseases. Cells. 2020;9:1627.
Rius-Pérez S, Pérez S, Martí-Andrés P, Monsalve M & Sastre J . NF-κBSignaling complexes in acute inflammation. Antioxid. Redox Signal 2019
Kaiser P & Stäheli P. Avian cytokines and chemokines,[w:] Davison F, Kaspers B, Schat KA (red.): Avian Immunology.
Balan KV, Babu US. Comparative responses of chicken macrophages to infection with Salmonella enterica serovars. Poult Sci. 2017;96:1849–54.
Thompson MG, Larson M, Vidrine A, Barrios K, Navarro F, Meyers K, Watkins SK. FOXO3–NF-κB RelA protein complexes reduce proinflammatory cell signaling and function. J Immunol. 2015;195:5637–47.
Albu DI, Feng D, Bhattacharya D, Jenkins NA, Copeland N, Liu P, Avram D. BCL11B is required for positive selection and survival of double-positive thymocytes. J Exp Med. 2007;204:3003–15.
Kastner P, Chan S, Vogel WK, Zhang LJ, Topark Ngarm A. Golonzhka, & Leid M Bcl11b represses a mature T-cell gene expression program in immature CD4+ CD8+ thymocytes. Eur J Immunol. 2010;40:2143–54.
Kojo S, Tanaka H, Endo TA, Muroi S, Liu Y, Seo W, Taniuchi I. Priming of lineage specifying genes by Bcl11b is required for lineage choice in post-selection thymocytes. Nat Commun. 2017;8:1–14.
Schaniel C, Bruno L, Melchers F, Rolink AG. Multiple hematopoietic cell lineages develop in vivo from transplanted Pax5- deficient pre-B I-cell clones. Blood. 2002;99:472–8.
Nutt SL, Morrison AM, Dörfler P, Rolink A, Busslinger M. Identification of BSAP (Pax-5) target genes in early B-cell development by loss-and gain-of-function experiments. EMBO J. 1998;17:2319–33.
Medvedovic J, Ebert A, Tagoh H, Busslinger M. Pax5: a master regulator of B cell development and leukemogenesis. Adv Immunol. 2011;111:179–82.
Schebesta A, McManus S, Salvagiotto G, Delogu A, Busslinger GA, Busslinger M. Transcription factor Pax5 activates the chromatin of key genes involved in B cell signalling, adhesion, migration, and immune function. Immunity. 2007;27:49–63.
Wang S, Xu X, Xu F, Meng Y, Sun C, Shi L, Zhao E. Combined expression of c-jun, c-fos, and p53 improves estimation of prognosis in oral squamous cell carcinoma. Cancer Invest. 2016;34:393–400.
Lin KI, Tunyaplin C, Calame K. Transcriptional regulatory cascades controlling plasma cell differentiation. Immunol Rev. 2003;194:19–28.
Shapiro-Shelef M, Calame K. Regulation of plasma-cell development. Nat Rev Immunol. 2005;5:230–42.
Abd El-Aziz DM. Detection of Salmonella typhimurium in retail chicken meat and chicken giblets. Asian Pac J Trop Biomed. 2013;3:678–81.
Yasuda T, Hayakawa F, Kurahashi S, Sugimoto K, Minami Y, Tomita A, Naoe T. B cell receptor-ERK1/2 signal cancels PAX5-dependent repression of BLIMP1 through PAX5 phosphorylation: a mechanism of antigen-triggering plasma cell differentiation. J Immunol. 2012;188:6127–34.
Teo JL, Kahn M. The Wnt signalling pathway in cellular proliferation and differentiation: a tale of two coactivators. Adv Drug Deliv Rev. 2010;62:1149–55.
Mu F, Huang J, Xing T, Jing Y, Cui T, Guo Y, Wang N. The Wnt/β-catenin/Lef1 pathway promotes cell proliferation at least in part through direct upregulation of miR-17-92 cluster. Front Genet. 2019;10:525.
Clevers H. Wnt/β-catenin signalling in development and disease. Cell. 2006;127:469–80.
Wen AY, Sakamoto KM, Miller LS. The role of the transcription factor CREB in immune function. J Immunol. 2010;185:6413–9.
Garrett-Sinha LA. Review of Ets1 structure, function, and roles in immunity. Cell Mol Life Sci. 2013;70:3375–90.
Maroulakou IG, Bowe DB. Expression and function of Ets transcription factors in mammalian development: a regulatory network. Oncogene. 2000;19:6432–42.
Ghysdael J, Gegonne A, Pognonec P, Dernis D, Leprince D, Stehelin D. Identification and preferential expression in thymic and bursal lymphocytes of a c-ets oncogene-encoded Mr 54,000 cytoplasmic protein. Proc Natl Acad Sci. 1986;83:1714–8.
Fisher CL, Ghysdael J, Cambier JC. Ligation of membrane Ig leads to calcium-mediated phosphorylation of the proto-oncogene product, Ets-1. J Immunol. 1991;146:1743–9.
John SA, Clements JL, Russell LM, Garrett-Sinha LA. Ets-1 regulates plasma cell differentiation by interfering with the activity of the transcription factor Blimp-1. J Biol Chem. 2008;283:951–62.
Inaba H, Kuboniwa M, Sugita H, Lamont RJ, Amano A. Identification of signalling pathways mediating cell cycle arrest and apoptosis induced by Porphyromonas gingivalis in human trophoblasts. Infect Immun. 2012;80:2847–57.
Abdel-Latif MA, Elbestawy AR, El-Far AH, Noreldin AE, Emam M, Baty RS, Albadrani GM, Abdel-Daim MM, El-Hamid A, Hatem S. Quercetin dietary supplementation advances growth performance, gut microbiota, and intestinal mrna expression genes in broiler chicken. Animals. 2021;11:2302.
Lin X, Wang R, Zou W, Sun X, Liu X, Zhao L, Wang S, Jin M. The influenza virus H5N1 infection can induce ROS production for viral replication and host cell death in A549 cells modulated by human Cu/Zn Superoxide Dismutase (SOD1) Overexpression. Viruses. 2016;8:13.
Cui C, Han S, Shen X, He H, Chen Y, Zhao J, Wei Y, Wang Y, Zhu Q, Li D, Yin H. ISLR regulates skeletal muscle atrophy via IGF1-PI3K/Akt-Foxo signaling pathway. Cell Tissue Res. 2020;381:479–92.
Peng M, Han J, Li L, Ma H. Suppression of fat deposition in broiler chickens by (-)-hydroxycitric acid supplementation: a proteomics perspective. Sci Rep. 2016;6:1–11.
Liu Z, Lai K, Xie Y, He X, Zhou X. Gli2 mediated activation of hedgehog signaling attenuates acute pancreatitis via balancing inflammatory cytokines in mice. Cell Physiol Biochem. 2018;48:120–30.
The Science and Engineering Research Board, Department of Science and Technology, Government of India is duly acknowledged for supporting this study under grant EMR/2017/000580.
Ethics approval and consent to participate
The tissue collection was accomplished only after accredited from Institutional Animal Ethics Committee on ethical standards in animal experimentation (AU/FVSc/PS-57/16021). All the guidelines of the Institutional Animal Ethics Committee of SKUAST Kashmir were strictly followed during all the experimental procedures. All the experiments were performed in accordance with ARRIVE guidelines.
Consent for publication
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary figure 1. The PCA plot shows between different samples. This plot shows clear separation between samples. Supplementary table 1. Downloaded RNA-Seq data from Mashooq et al 2022 UR_Genomics_submission number (NCBI Accession number GSE 168060).
About this article
Cite this article
Ahmad, S.M., Bhat, S.S., Shafi, S. et al. Identification of key transcription factors and their functional role involved in Salmonella typhimurium infection in chicken using integrated transcriptome analysis and bioinformatics approach. BMC Genomics 24, 214 (2023). https://doi.org/10.1186/s12864-023-09315-3