- Research Article
- Open Access
Cecal MicroRNAome response to Salmonella enterica serovar Enteritidis infection in White Leghorn Layer
BMC Genomics volume 18, Article number: 77 (2017)
Salmonella enterica serovar Enteritidis (SE) is a food-borne pathogen and of great threat to human health through consuming the contaminated poultry products. MicroRNAs (miRNAs) play an important role in different biological activities and have been shown to regulate the innate immunity in the bacterial infection. The objective of this study is to identify miRNAs associated with SE infection in laying chicken cecum.
Average number of reads of three libraries constructed from infected and non-infected chickens was 12,476,156 and 10,866,976, respectively. There were 598 miRNAs including 194 potential novel miRNAs identified in which 37 miRNAs were significantly differentially expressed between infected and non-infected chickens. In total, 2897 unique target genes regulated by differentially expressed miRNAs were predicted, in which, 841 genes were uniquely regulated by up-regulated miRNAs (G1), 636 genes were uniquely regulated by down-regulated miRNAs (G2), and 1420 were co-regulated by both up and down- regulated miRNAs (G3). There were 118, 73 and 178 GO (Gene ontology) BP (Biological process) terms significantly enriched in G1, G2 and G3 groups, respectively. More immune-related GO BP terms than metabolism-related terms were found in G1. Expression of 12 immune-related genes of four differentially expressed miRNAs was detected through qRT-PCR. The regulatory direction of gga-miR-1416-5p, gga-miR-1662, and gga-miR-34a-5p were opposite with the target genes of TLR21, BCL10, TLR1LA, NOTCH2 and THBS1, respectively.
The miRNAs contribute to the response to SE infection at the onset of egg laying through regulating the homeostasis between metabolism and immunity. The gga-miR-125b-5p, gga-miR-34a-5p, gga-miR-1416-5p and gga-miR-1662 could play an important role in SE infection through regulating their target genes. The finding herein will pave the foundation for the studies of microRNA regulation in SE infection in laying hens.
Salmonellosis is one of the most important food-borne diseases, with an estimated one million human cases and more than 350 deaths occurring each year in the United States . Non-typhoidal Salmonella is the main cause of foodborne disease in the United States . Salmonella enterica serotype Enteritidis (SE) is one of the most common serotypes of Salmonella bacteria reported worldwide which is the major source of human intestinal infections . The animal or animal product including poultry, eggs, milk and milk products could be the carriers of SE which threaten the human health [4–6]. Salmonella organisms could penetrate the intestines mucosal epithelium to outcompete the resident microbial flora and cause the pathological reaction . Egg-related salmonellosis is costing $44 million per year in Australia .
MicroRNAs (miRNAs) are small (19 ~ 24 nt long), non-coding, single-stranded and conserved RNAs. It performs multiple functions through regulating gene expression mainly at the post-transcriptional levels . It plays an important role in different biological activities such as the development, cell differentiation and disease [10, 11]. MiRNA could regulate innate immunity caused by viruses, bacteria, fungi, and protozoa infection [12, 13]. MiR-21, miR-146a/b and miR-155 were obviously up-regulated in rat’s mononuclear cells after Salmonella infection [14, 15]. Let-7 was down-regulated to induce the release of cytokine IL6 (interleukin 6) and IL10 to participate in the regulation of immune response to Salmonella infection in macrophages . MiRNAs also could buffer and alter the variance of relatively lowly expressed genes in the response to Salmonella infection in pig . However, the responsive miRNAs in laying hen Salmonella infection is still unclear.
The objective of the present study is to discover the miRNAs in the response to Salmonella infection in chicken cecum at the onset of egg laying. Next generation sequencing (NGS) has been widely used to analyze the miRNA expression profile in many studies [17–19]. In the current study, the Solexa Sequencing method was used to detect differentially expressed miRNAs in the response to SE infection in White Leghorn at the onset of egg laying. Our results will expand the list of miRNAs related to the host responses to SE infection in chickens.
Preliminary analysis of the raw data
Six libraries were constructed from infected (I) and non-infected (N) chickens, three in each group. The average number of total reads of three libraries obtained from I and N chickens were 12,476,156 and 10,866,976, respectively, and the filtered clean reads were 5,078,218 and 2,411,757, respectively (Table 1). In the I group, 3,456,099 clean reads were exactly matched to the chicken genome, 525,400 1-mismatched and 108,131 2-mismatched. In the N group, 1,722,678 clean reads were exactly matched to chicken reference genome, 149,372 1-mismatched and 47,761 2-mismatched. Altogether, 80.5% (4,089,630) and 79.6% (1,919,811) of clean reads were mapped to the genome in the I and N group, respectively. The miRNAs with 21 nt in length were most abundant followed by 22 nt in both I and N groups.
Number of reads in each chromosome (Chr) was counted (Additional file 1 and Fig. 1). In general, number of read counts in I group was more than that in N group in each chromosome. For the I group, there were 769,514, 929,207 and 810,295 clean reads matched on Chr1, Chr3, and Chr13, which account for 19.81, 18.82 and 22.72% of total clean reads, respectively. For the N group, there were 367,058 and 693,717 clean reads matched on Chr1 and Chr3, which account for 19.12 and 36.13% of total clean reads, respectively. The count of matched clean reads on Chr13 in I group was eight folds more than that in N group.
MiRNA identification and genome distribution
To identify the miRNAs, small RNA sequences obtained by deep sequencing were aligned with known miRNA sequences in the miRBase and Rfam (http://pfam.xfam.org). The novel miRNAs were predicted through miRDeep. There were 598 miRNAs including 194 novel miRNAs identified. The number of identified miRNAs in each chromosome was counted (Fig. 2). The chromosomes could be divided into three clusters based on the number of mapped miRNAs: (1) Chr18, Chr22, Chr24, Chr25, Chr27, less than ten miRNAs in each chromosome, (2) Chr5, Chr6, Chr8-12, Chr19, Chr21, Chr23, Chr26 and Chr28, 10–20 miRNAs in each chromosome, (3) Chr1-4, Chr7, Chr13-17, Chr20 and ChrZ, more than 20 miRNAs in each chromosome. There were 77 miRNAs mapped on Chr1 and only two miRNAs on Chr22, respectively. The density distribution of miRNAs across chromosomes showed that number of miRNAs in 1Mbp DNA in each chromosome ranged from 0.19 (ChrW) to 2.77 (Chr21) (Additional file 2). The density of miRNAs in Chr1-12 was lower than one miRNA per 1 Mbp DNA. The density of miRNAs in Chr21, Chr23, Chr26 and Chr28 were more than two miRNAs per 1 Mbp DNA. There were 20–60% of known miRNAs in the miRBase observed in each chromosome except for the Chr21 (91.67%), Chr23 (73.33%), Chr18 (15.79%) and ChrUn (12.5%).
Differentially expressed miRNAs responding to SE infection
The differentially expressed miRNAs between I and N groups were identified through edgeR package. There were 37 miRNAs significantly differentially expressed between I and N groups including 19 known miRNAs and 18 potentially novel miRNAs (P < 0.05 with false discovery rate (FDR) of 0.49 and fold change > 2) (Table 2). There were 22 miRNAs up-regulated post SE infection, which included 15 known and seven novel miRNAs. The highest fold change (71.24) was observed for gga-miR-490- 5p and the lowest fold change (2.06) for gga-miR-193b-3p. Fifteen miRNAs were down-regulated in the comparison of I/N, which included four known and 11 novel miRNAs. The highest fold change (15.13) was observed for gga-miR-chr17_13654 and the lowest (2.01) for gga-miR-chr15_12378. Gga-miR-490-5p, gga-miR-chr13_10137 and gga-miR-chrUn_AADN03024004_45551 were unique to the infected chickens. More than half of significantly differentially expressed miRNAs (17/37) were located on Chr1-10.
The heat map and hierarchical clustering demonstrated that the miRNA profiles from I and N groups were distinct (Fig. 3). The infected and non-infected groups were separated distinctly. The differentially expressed miRNAs were clustered into two groups. MiRNAs in group one showed higher expression in the N group but lower expression in I group, which including gga-miR-34a-5p, gga-miR-215-5p, and gga-miR-1662. MiRNAs in group two showed higher expression in I group but lower expression in N group including gga-miR-125b-3p, gga-miR-1416-5p and gga-miR-125b-5p.
MiRNA target gene identification
To investigate the biological role of the differentially expressed miRNAs, the potential target genes were predicted using Miranda algorithm with the Vienna package. In total, 2897 unique target genes regulated by those differentially expressed miRNAs were predicted, in which, 841 genes were uniquely regulated by up-regulated miRNAs (G1), 636 genes were uniquely regulated by down-regulated miRNAs (G2), and 1420 were co-regulated by both up and down-regulated miRNAs (G3). One hundred and seventy-six immune-related genes were retrieved from the Ensembl BioMart database . In total, 587 pairs of interaction between miRNAs and immune-related genes were obtained (Additional file 3). Gga-miR-34a-5p interacted with the greatest number of immune-related genes (46), while the gga-miR-215-5p only interacted with one immune-related gene.
Gene ontology (GO) analysis for target genes
To demonstrate the function of target genes of differentially expressed miRNAs, functional annotation was performed through DAVID (The Database for Annotation, Visualization and Integrated Discovery) for the target genes in the G1, G2 and G3 groups, respectively (Additional file 4). There were 118, 73 and 178 GO BP (biological process) terms significantly enriched in G1, G2 and G3, respectively (P < 0.05). The GO BP terms were categorized by CateGOrizer using “Immune System Gene Classes” GO classification (Fig. 4). For target genes in G1, the enriched GO BP terms were categorized into ten ancestral classes, which could be divided into three groups: 1) immune-related function including death, apoptosis, cell adhesion, stress response, response to abiotic stimulus, regulation of apoptosis, and response to external stimulus occupied 61.54% of all ancestral classes, 2) metabolism-related including carbohydrate metabolism, protein metabolism, and catabolism, occupied 30.76% of all ancestral classes, 3) cell adhesion occupied 7.69% of all ancestral classes. For target genes in G2, the enriched GO BP terms were categorized into 12 ancestral classes. The immune-related classes including stress response, death, response to abiotic stimulus, lymphocyte activation, apoptosis, T cell activation, lymphocyte differentiation, regulation of apoptosis, occupied 36.83% of all ancestral classes. The metabolism-related classes including catabolism, protein metabolism, carbohydrate metabolism, occupied 60.53% of all ancestral classes. For target genes in G3, the enriched BP terms were categorized into 20 ancestral classes, the immune-related classes occupied 80.58% of all ancestral classes, and metabolism-related classes occupied 15.28%.
KEGG Pathway analysis
For the target genes in the G1, only one pathway of endocytosis was significantly enriched (P < 0.05) with the fold enrichment of 1.65, which included 18 target genes. For target genes in the G2, five pathways were significantly enriched, which were ECM-receptor interaction, glycolysis/gluconeogenesis, focal adhesion, melanogenesis and proteasome. Proteasome had the highest fold enrichment of 3.53 and focal adhesion included the largest number of target genes (16). For target genes in the G3, Cell cycle and SNARE interactions in vesicular transport were significantly enriched and associated with 27 and 9 genes, respectively (Table 3).
Protein-protein Interacting Network of the immune-related target genes
The Reactome FIViz (http://f1000research.com/articles/3-146/v2) app in Cytoscape was used to examine the potential protein-protein interactions of the 176 immune-related genes encoded proteins . Through Cytoscape, there were 98 proteins accepted by Reactome FIViz and clustered into six network modules (Fig. 5). The top ten proteins interacted with more than 15 other proteins were illustrated in the center of the network. JUN protein interacted with 32 proteins, and LYN and JAK with 21 proteins, respectively.
Validation of differentially expressed miRNAs
The relative expression level of 12 differentially expressed miRNAs regulating immune-related genes was validated through quantitative real-time PCR (qRT-PCR). The specific primers for those 12 miRNAs (seven up-regulated miRNAs (including two novel miRNAs) and five down-regulated miRNAs (including two novel miRNAs)) were listed in Additional file 5. The results showed that all the miRNAs selected for qRT-PCR validation were significantly differentially expressed. All the results in qRT-PCR were consistent with those obtained from the sequencing except for gga-miR-193b-3p which had reverse regulatory direction (Table 4). The gga-miR-1416-5p and gga-miR-125b-5p were up-regulated. The gga-miR-1662 and gga-miR-34a-5p were down-regulated after SE infection.
Target genes expression
Twelve immune-related target genes of 4 differentially expressed miRNAs (gga-miR-1416-5p, gga-miR-1662, gga-miR-125b-5p and gga-miR-34a-5p) were selected to detect the relative expression using qRT-PCR. The potential regulation was listed in Additional file 6. The specific primers were listed in Additional file 7. The results showed that eight target genes were significantly differentially expressed between I and N groups (Fig. 6). NOTCH2, THBS1, RIPK2, IGJ and TLR1LA were significantly up-regulated (P < 0.05) and CCL4, TLR21 and BCL10 were significantly down-regulated following SE infection. Five of eight target genes had reversely regulatory direction with their regulating miRNAs following SE infection.
MiRNAs play a vital role in regulating immunity induced by bacteria [22–24]. Regulation of miRNA expression responding to bacterial infection is emerging as a crucial part of the host response to infection . The onset of egg laying is a complex physiological process when physiological and hormonal stress occurs. Previously, we have studied the splenic gene expression in laying hens through microarray . In the current study, cecal miRNA profiling in the response to SE infection in laying hens has been identified through Solexa sequencing. In total, 404 known miRNAs and 194 potential novel miRNAs were expressed in six libraries.
MiRNA distribution varies widely in different chromosome in different species, and chromosomes with higher number of miRNAs most likely play regulatory roles in several cellular processes including disease. MiRNAs associated with cancer and cardiovascular disease are mostly located in Chr1, 14, 19, and X in human . There is a significant association between the chromosomal location of miRNAs and those of cancer susceptibility loci in mouse . There are 39.66% miRNA genes located on chicken chromosome 1, 2, 3, and 4 . The skew distribution of miRNAs has been observed following SE infection in the current study. Majority of miRNAs responding to SE infection located on Chr1, 2, 3, 10, 13, and 20. There was no miRNAs observed in chicken Chr16 which is consistent with the current results from miRBase database (http://www.mirbase.org).
SE infection stimulates the miRNAs expression. Thirty-seven miRNAs were differentially expressed following SE infection. Read counts in the infected group were higher than those in the non-infected group (Fig. 1). In the current study, more up-regulated miRNAs were observed at 7 days post infection (dpi) (22 up-regulated vs. 15 down-regulated). More miRNAs are up-regulated at 2 days post infection with S. Typhimurium in piglet which is consistent with current miRNA profiling . The similar response to Gram-negative bacterial infection could be regulated by miRNAs. MiRNA regulation in immune response is influenced by genetic background. It has been reported that opposite miRNA profiling was found between SPF layer and broiler chickens infected with AIV (Avian Influenza Virus) [19, 29]. Functionally, the miRNAs target mRNAs and trigger either translation repression or RNA degradation in animal . It has been reported that the complicated interaction between immune system and metabolism exists during SE infection in laying hens . Maintaining metabolic homeostasis requires a balanced immune response. The perturbation of this equilibrium could result in pathological situations . It has been reported that miRNAs are fine-tuners of metabolic processes  and could generate a bidirectional functional link between metabolism and pathogenesis . In the current study, target genes of up-regulated miRNAs were associated with stronger immune function than metabolism function. Target genes of down-regulated miRNAs were associated with stronger metabolism compared with immune function. The target genes regulated by both up- and down-regulated miRNAs were also associated with stronger immune function than metabolism (Fig. 4). The target genes of up- and down regulated miRNAs may buffer the homeostasis between immune response and metabolism. Those miRNAs may contribute to buffer the gene expression of target genes [16, 32]. The immune function was repressed at the onset of egg laying [25, 34, 35]. The results herein suggested that miRNAs could regulate the interaction between immune response and metabolism on day 7 post SE infection in laying hens.
It is noteworthy that the Proteasome pathway had the highest fold enrichment in the enriched pathways. The proteasome is a protein-destroying apparatus involved in many essential cellular functions including antigen processing for appropriate immune responses and inflammatory responses [36–38]. The nonclassical class Ib molecule Qa-1b is a dominant restricting element to the recognition of CD8+ T Cells which is proteasome-dependent during Salmonella infection . Type III protein secretion systems (T3SS) effector protein activities of Salmonella are temporally regulated by proteasome-dependent protein degradation . The LMP2 protein, a subunit of the cytosolic proteasome complex, is up-regulated after invasion of HLA-B27-transfected HeLa cells by Salmonella typhimurium . Both Salmonella enteritidis and Salmonella typhimurium are gram-negative bacteria. This could suggest that proteasome is an important defense system to protect against bacterial colonization regulated by miRNAs in chicken.
MiR-34a-5p was down-regulated and potentially targeted 46 immune-related genes in chicken after SE infection. MiR-34a was associated with the cellular senescence, it was strongly up-regulated and targeted the important proto-oncogene MYC during B-RAF-induced senescence . MiR-34a may play an important role in the loss of oxidative defense in rat liver during aging through suppressing the expression of Sirt1, Mgst1, Sp1 and Nrf2 . Bcl-2 and the Foxp1 transcription factor, which were required for early B cell development , are the direct target of miR-34a, so it was hypothesized that miR-34a plays an important role in B cell development [45–47]. MiR-34a causes a partial block in B cell development, whereas its knockdown results in increased B cells development in rat bone marrow . CCL4, IL8L1 and CDC42 were contained in the Salmonella infection pathway and were all targeted by miR-34a. CCL4 gene is significantly up-regulated at different time points (1, 2, 4, and 8 h) after SE endotoxin treatment . CDC42 could activate the Arp2/3 by N-WASP (Neural Wiskott-Aldrich syndrome protein) to mediate actin polymerization in the invasion of SE [50–53]. It has been reported that mir-34a exclusively expressed in AIV infected chicken lung with targeting 14 immune-related genes and four AIV genes . Although there were no SE genes regulated by mir-34a predicted, the down-regulation of miR-34a may conduce to B cell development in the response to SE infection in chickens.
Gga-miR-1416 is located in the intron of MAMDC 2 (MAM domain containing 2) on ChrZ. Gga-miR-1416 is up-regulated in both layer and broiler following AIV infection suggests that it responds to AIV infection across diverse genetic lines . Six immune-related genes were potentially targeted by gga-miR-1416-5p including BCL10, NFKBIA and TLR21, which were important in the response to bacterial infection. BCL10 mediates the LPS-induced activation of NF-κB and IL8 in normal human intestinal epithelial cells . Chicken TLR21 is involved in the recognition of bacterial components and Salmonella in the innate immune response . The expression of TLR 21 was decreased in the duodenum, jejunum, ileum, ceca and large intestine of broilers following SE infection . In the current study, the expression of TLR 21 was significantly down-regulated and had reversely regulatory direction with the regulating gga-miR-1416-5p post SE infection. Gga-miR-1416-5p could play an important role in response to SE infection through regulating the target genes.
The interaction of miRNA with target genes plays an important role in many biological functions. Toll-like receptors (TLR) recognizes specific patterns of microbial components and participate in the innate immunity and antigen-specific adaptive immunity . As one of the main pattern recognition receptors, TLR can identify the pathogen to activate the immune cells response and induce to produce type I interferon and a series of proinflammatory cytokines . TLR1LA could combine with TLR2 to form a dimer and efficiently identify bacterial peptidoglycan and lipoprotein . TLR1LA is significantly up-regulated at the ileum, cecum and colon at 24 h post SE infection in day-old chickens . Gga-miR-1662 was significantly up-regulated in chicken lung following AIV infection . Gga-miR-1662 was down-regulated while its target gene TLR1LA was up-regulated in the response to SE infection in the current study. The interaction of gga-miR-1662 with TLR1LA could be important in the respond to SE infection in chicken. Further study is warranted to verify the function of interaction between gga-miR-1662 and TLR1LA.
MiRNAs mediate the homeostasis between metabolism and immunity in the response to SE infection at the onset of egg laying. The gga-miR-34a-5p, gga-miR-1416-5p and gga-miR-1662 could play important roles in SE infection through regulating their target genes. The finding herein will pave the foundation for the studies of microRNA regulation in SE infection in laying hens.
Animals and SE challenges
One hundred 20-week old SE negative White Leghorn layers were used in the current study. Chickens were randomly divided into two groups with equal number, infected group (I) and non-infected (N) group. At the first day, chickens in the infected group were fed with 30 g feed containing 5.8 × 108 cfu SE at a time and the non-infected birds were fed with regular feed as described previously (Wu et al., ). Chickens in each group were sacrificed by cervical dislocation at 7 days post infection (dpi). The cecum samples were aseptically harvested from each chicken and put into RNAlater solution (Life technologies, Grand Island, USA) and stored at –20 °C. All animal procedures were approved by Shandong Agricultural University Animal Care and Use Committee.
Small RNAs library construction and deep sequencing
In total, nine infected and six non-infected birds were randomly selected for further RNA isolation. Total RNA was isolated from each individual sample using Trizol following the manufacture’s protocol (Life technologies, Grand Island, USA). The RNA integrity and concentration were checked and measured using gel electrophoresis and Nanodrop (Thermo Fisher Scientific, Wilmington, USA).
Six small RNA libraries were pooled from I and N groups, three in each group. Each small RNA library from I (tp1, tp2, tp3) and N (cp1, cp2, cp3) group consisted of three and two individual samples with equal amount of RNA, respectively. The small RNA pools were purified and enriched using denaturing polyacrylamide gel electrophoresis. A pair of Illumina proprietary adaptors was ligated to their 5’ and 3’ ends, followed by reverse transcription and cluster generation using TruSeq Small RNA Library Preparation Kit (Illumina Inc., San Diego, USA). Subsequently, the libraries were sequenced by Illumina Hiseq 2500 according to the manufacturer’s instructions at Genergy Inc (Shanghai, China).
Basic data processing
The small RNA sequence reads were pre-processed using FASTX-Toolkit to filter low-quality reads and trim adaptor. After filtering adaptor sequences and removing contaminated reads, the clean reads were matched to chicken reference genome using the Bowtie . The first 15 bp of the read was exactly matched. Two mismatched bases were accepted on the rest nucleotides of the read. The filtered sequences were matched with miRBase (http://mirbase.org/) to search for known miRNAs with exact matches. The unmatched data sets were aligned with chicken genomic sequence (Gallus_gallus-4.0) and predict the novel miRNAs using miRDeep .
Differential expression analysis of miRNAs
To compare miRNAs expression level between I and N groups, read count of each identified miRNA was normalized to the total number of reads in each given sample. The differentially expressed miRNAs were identified through edgeR package with False discovery rate multiple testing correction . The P < 0.05 and fold change > 2 was considered as significant difference.
Target prediction and Gene Ontology (GO) enrichment analysis
The target genes of those significantly differentially expressed miRNAs were predicted using Miranda algorithm through the Vienna package [63, 64]. Functional annotation of GO and pathway analysis for those target genes of miRNAs were performed through DAVID 6.7 [65–67]. Significant over-representation is based on a Fisher Exact statistical methodology similar to that described by Al-Shahrour et al . CateGOrizer was used to categorize the significantly enriched GO BP (biological process) terms .
Protein-protein Interaction Network Analysis of the immune-related target genes
We converted the gene IDs of 176 immune-related target genes to the gene symbols of the encoded proteins using bioDBnet software (www.biodbnet.abcc.ncifcrf.gov/db/db2 db.php#biodb). Then the gene symbols were mapped onto the chicken functional interaction network found in the Reactome database using the Reactome FI network plug-in in the Cytoscape software .
Quantitative real-time PCR of miRNAs and target genes
The same individual RNA samples used for sequencing were used for quantitative real-time PCR (qRT-PCR). In total, nine infected and six non-infected samples were used. Twelve differentially expressed miRNAs were validated and characterized using qRT-PCR. In brief, 1 μg of total RNA was reverse transcribed using Step PrimeScript® miRNA cDNA Synthesis Kit (Perfect Real Time) and amplified using Stratagene MX3000 real-time PCR System with miRNA specific primers (Additional file 5). The qRT-PCR was performed with SYBR green PCR master mix. Small nucleolar RNA U6 was used as endogenous control to normalize RNA input. The 20 μL PCR reactions included 10 μL SYBR® Premix Ex TaqTM II (2×), 0.4 μL ROX Reference Dye (50×), 0.8 μL Uni-miR qPCR Primer (10 μM), 0.8 μL forward qPCR Primer, 2 μL cDNA and 6 μL ddH2O. The qPCR amplification conditions were: 1 cycle of 95 °C for 30s, 40 cycles of 95 °C for 5 s and 60 °C for 30s. All qRT-PCR reactions were performed in triplicate. The relative expression were calculated using 2-△△CT method. The student’s T-test was performed to examine the significance of miRNA expression between I and N groups.
Relative expression level of 12 immune-related target genes were quantified using qRT-PCR. The RNA was reverse transcribed to cDNA using TaKaRa Primer Script™ RT reagent kit (Perfect Real Time) (TaKaRa, Dalian, China) according to the manufacturer’s manual. The qRT-PCR was performed by Stratagene MX3000 real-time PCR System with SYBR green method. The specific primers were designed by primer premier 5.0 according to gene sequence (Additional file 7). The 20 μL PCR reactions included 10 μl SYBR Primer Ex TaqTM (2×), 0.4 μl forward primer (10 μM), 0.4 μl reverse primer (10 μM), 0.4 μL ROX Reference Dye II (50×), 2 μL cDNA, 6.8 μL ddH2O. The conditions of qRT-PCR amplification and data analysis were the same as those used in miRNA qRT-PCR. The chicken β-actin was used as the internal control.
Salmonella enterica serovar Enteritidis
Genes were uniquely regulated by up-regulated miRNAs
Genes were uniquely regulated by down-regulated miRNAs
Genes co-regulated by both up and down- regulated miRNAs
Next generation sequencing
Chromosome; ChrUn contains clone contigs that can’t be confidently placed on a specific chromosome
The Database for Annotation, Visualization and Integrated Discovery
Kyoto Encyclopedia of Genes and Genomes
Quantitative real-time PCR
Avian influenza virus
Type III protein secretion systems
Toll like receptor
False discovery rate
Scallan E, Hoekstra RM, Angulo FJ, Tauxe RV, Widdowson MA, Roy SL, Jones JL, Griffin PM. Foodborne illness acquired in the United States--major pathogens. Emerg Infect Dis. 2011;17(1):7–15.
Scallan E, Hoekstra RM, Mahon BE, Jones TF, Griffin PM. An assessment of the human health impact of seven leading foodborne pathogens in the United States using disability adjusted life years. Epidemiol Infect. 2015;143(13):2795–804.
Karaffova V, Bobikova K, Husakova E, Levkut M, Herich R, Revajova V, Levkutova M, Levkut M. Interaction of TGF-beta4 and IL-17 with IgA secretion in the intestine of chickens fed with E. faecium AL41 and challenged with S. Enteritidis. Res Vet Sci. 2015;(100):75–9.
Gomez TM, Motarjemi Y, Miyagawa S, Kaferstein FK, Stohr K. Foodborne salmonellosis. World Health Stat Q. 1997;50(1-2):81–9.
Foley SL, Lynne AM. Food animal-associated Salmonella challenges: pathogenicity and antimicrobial resistance. J Anim Sci. 2008;86(14 Suppl):E173–87.
Bosilevac JM, Guerini MN, Kalchayanand N, Koohmaraie M. Prevalence and characterization of salmonellae in commercial ground beef in the United States. Appl Environ Microbiol. 2009;75(7):1892–900.
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. Infect Immun. 2007;75(12):5993–6007.
Samiullah, Chousalkar KK, Roberts JR, Sexton M, May D, Kiermeier A. Effects of egg shell quality and washing on Salmonella Infantis penetration. Int J Food Microbiol. 2013;165(2):77–83.
Guarnieri DJ, DiLeone RJ. MicroRNAs: a new class of gene regulators. Ann Med. 2008;40(3):197–208.
Ebert MS, Sharp PA. Roles for microRNAs in conferring robustness to biological processes. Cell. 2012;149(3):515–24.
Yang C, Wei W. The miRNA expression profile of the uveal melanoma. Sci China Life Sci. 2011;54(4):351–8.
Zhou R, O’Hara SP, Chen XM. MicroRNA regulation of innate immune responses in epithelial cells. Cell Mol Immunol. 2011;8(5):371–9.
Takeuchi O, Akira S. Pattern recognition receptors and inflammation. CELL. 2010;140(6):805–20.
Schulte LN, Eulalio A, Mollenkopf HJ, Reinhardt R, Vogel J. Analysis of the host microRNA response to Salmonella uncovers the control of major cytokines by the let-7 family. EMBO J. 2011;30(10):1977–89.
Sharbati S, Sharbati J, Hoeke L, Bohmer M, Einspanier R. Quantification and accurate normalisation of small RNAs through new custom RT-qPCR arrays demonstrates Salmonella-induced microRNAs in human monocytes. BMC Genomics. 2012;13:23.
Bao H, Kommadath A, Plastow GS, Tuggle CK, Guan LL, Stothard P. MicroRNA buffering and altered variance of gene expression in response to Salmonella infection. PLoS ONE. 2014;9(4):e94352.
Lawless N, Foroushani AB, McCabe MS, O’Farrelly C, Lynn DJ. Next generation sequencing reveals the expression of a unique miRNA profile in response to a gram-positive bacterial infection. PLoS ONE. 2013;8(3):e57543.
Rathjen T, Pais H, Sweetman D, Moulton V, Munsterberg A, Dalmay T. High throughput sequencing of microRNAs in chicken somites. Febs Lett. 2009;583(9):1422–6.
Wang Y, Brahmakshatriya V, Zhu H, Lupiani B, Reddy SM, Yoon BJ, Gunaratne PH, Kim JH, Chen R, Wang J, et al. Identification of differentially expressed miRNAs in chicken lung and trachea with avian influenza virus infection by a deep sequencing approach. BMC Genomics. 2009;10:512.
Kinsella RJ, Kahari A, Haider S, Zamora J, Proctor G, Spudich G, Almeida-King J, Staines D, Derwent P, Kerhornou A, et al. Ensembl BioMarts: a hub for data retrieval across taxonomic space. Database (Oxford). 2011;2011:r30.
Wu G, Dawson E, Duong A, Haw R, Stein L. ReactomeFIViz: a Cytoscape app for pathway and network-based data analysis. F1000Res. 2014;3:146.
Maudet C, Mano M, Eulalio A. MicroRNAs in the interaction between host and bacterial pathogens. FEBS Lett. 2014;588(22):4140–7.
Maudet C, Mano M, Sunkavalli U, Sharan M, Giacca M, Forstner KU, Eulalio A. Functional high-throughput screening identifies the miR-15 microRNA family as cellular restriction factors for Salmonella infection. Nat Commun. 2014;5:4718.
Staedel C, Darfeuille F. MicroRNAs and bacterial infection. Cell Microbiol. 2013;15(9):1496–507.
Wu G, Liu L, Qi Y, Sun Y, Yang N, Xu G, Zhou H, Li X. Splenic gene expression profiling in White Leghorn layer inoculated with the Salmonella enterica serovar Enteritidis. Anim Genet. 2015;46(6):617–26.
Ghorai A, Ghosh U. miRNA gene counts in chromosomes vary widely in a species and biogenesis of miRNA largely depends on transcription or post-transcriptional processing of coding genes. Front Genet. 2014;5:100.
Sevignani C, Calin GA, Nnadi SC, Shimizu M, Davuluri RV, Hyslop T, Demant P, Croce CM, Siracusa LD. MicroRNA genes are frequently located near mouse cancer susceptibility loci. Proc Natl Acad Sci U S A. 2007;104(19):8017–22.
Uribe JH, Collado-Romero M, Zaldivar-Lopez S, Arce C, Bautista R, Carvajal A, Cirera S, Claros MG, Garrido JJ. Transcriptional analysis of porcine intestinal mucosa infected with Salmonella Typhimurium revealed a massive inflammatory response and disruption of bile acid absorption in ileum. Vet Res. 2016;47:11.
Wang Y, Brahmakshatriya V, Lupiani B, Reddy SM, Soibam B, Benham AL, Gunaratne P, Liu HC, Trakooljul N, Ing N, et al. Integrated analysis of microRNA expression and mRNA transcriptome in lungs of avian influenza virus infected broilers. BMC Genomics. 2012;13:278.
Lim LP, Lau NC, Garrett-Engele P, Grimson A, Schelter JM, Castle J, Bartel DP, Linsley PS, Johnson JM. Microarray analysis shows that some microRNAs downregulate large numbers of target mRNAs. Nat. 2005;433(7027):769–73.
Matarese G, La Cava A. The intricate interface between immune system and metabolism. Trends Immunol. 2004;25(4):193–200.
Dumortier O, Hinault C, Van Obberghen E. MicroRNAs and metabolism crosstalk in energy homeostasis. Cell Metab. 2013;18(3):312–24.
Rottiers V, Naar AM. MicroRNAs in metabolism and metabolic disorders. Nat Rev Mol Cell Biol. 2012;13(4):239–50.
Johnston CE, Hartley C, Salisbury AM, Wigley P. Immunological changes at point-of-lay increase susceptibility to Salmonella enterica Serovar enteritidis infection in vaccinated chickens. PLoS ONE. 2012;7(10):e48195.
Wigley P, Hulme SD, Powers C, Beal RK, Berchieri AJ, Smith A, Barrow P. Infection of the reproductive tract and eggs with Salmonella enterica serovar pullorum in the chicken is associated with suppression of cellular immunity at sexual maturity. Infect Immun. 2005;73(5):2986–90.
Adams J. The proteasome: structure, function, and role in the cell. Cancer Treat Rev. 2003;29 (Suppl 1):3–9.
Song R, Harding CV. Roles of proteasomes, transporter for antigen presentation (TAP), and beta 2-microglobulin in the processing of bacterial or particulate antigens via an alternate class I MHC processing pathway. J Immunol. 1996;156(11):4182–90.
Kruger E, Kuckelkorn U, Sijts A, Kloetzel PM. The components of the proteasome system and their role in MHC class I antigen processing. Rev Physiol Biochem Pharmacol. 2003;148:81–104.
Lo WF, Ong H, Metcalf ES, Soloski MJ. T cell responses to Gram-negative intracellular bacterial pathogens: a role for CD8+ T cells in immunity to Salmonella infection and the involvement of MHC class Ib molecules. J Immunol. 1999;162(9):5398–406.
Kubori T, Galan JE. Temporal regulation of salmonella virulence effector function by proteasome-dependent protein degradation. Cell. 2003;115(3):333–42.
Maksymowych WP, Ikawa T, Yamaguchi A, Ikeda M, McDonald D, Laouar L, Lahesmaa R, Tamura N, Khuong A, Yu DT, et al. Invasion by Salmonella typhimurium induces increased expression of the LMP, MECL, and PA28 proteasome genes and changes in the peptide repertoire of HLA-B27. Infect Immun. 1998;66(10):4624–32.
Christoffersen NR, Shalgi R, Frankel LB, Leucci E, Lees M, Klausen M, Pilpel Y, Nielsen FC, Oren M, Lund AH. p53-independent upregulation of miR-34a during oncogene-induced senescence represses MYC. Cell Death Differ. 2010;17(2):236–45.
Li N, Muthusamy S, Liang R, Sarojini H, Wang E. Increased expression of miR-34a and miR-93 in rat liver during aging, and their impact on the expression of Mgst1 and Sirt1. Mech Ageing Dev. 2011;132(3):75–85.
Hu H, Wang B, Borde M, Nardone J, Maika S, Allred L, Tucker PW, Rao A. Foxp1 is an essential transcriptional regulator of B cell development. Nat Immunol. 2006;7(8):819–26.
Shi C, Zhang X, Chen Z, Sulaiman K, Feinberg MW, Ballantyne CM, Jain MK, Simon DI. Integrin engagement regulates monocyte differentiation through the forkhead transcription factor Foxp1. J Clin Invest. 2004;114(3):408–18.
Shi C, Sakuma M, Mooroka T, Liscoe A, Gao H, Croce KJ, Sharma A, Kaplan D, Greaves DR, Wang Y, et al. Down-regulation of the forkhead transcription factor Foxp1 is required for monocyte differentiation and macrophage function. Blood. 2008;112(12):4699–711.
Bommer GT, Gerin I, Feng Y, Kaczorowski AJ, Kuick R, Love RE, Zhai Y, Giordano TJ, Qin ZS, Moore BB, et al. p53-mediated activation of miRNA34 candidate tumor-suppressor genes. Curr Biol. 2007;17(15):1298–307.
Rao DS, O’Connell RM, Chaudhuri AA, Garcia-Flores Y, Geiger TL, Baltimore D. MicroRNA-34a perturbs B lymphocyte development by repressing the forkhead box transcription factor Foxp1. Immun. 2010;33(1):48–59.
Ciraci C, Tuggle CK, Wannemuehler MJ, Nettleton D, Lamont SJ. Unique genome-wide transcriptome profiles of chicken macrophages exposed to Salmonella-derived endotoxin. BMC Genomics. 2010;11:545.
Ma L, Cantley LC, Janmey PA, Kirschner MW. Corequirement of specific phosphoinositides and small GTP-binding protein Cdc42 in inducing actin assembly in Xenopus egg extracts. J Cell Biol. 1998;140(5):1125–36.
Ma L, Rohatgi R, Kirschner MW. The Arp2/3 complex mediates actin polymerization induced by the small GTP-binding protein Cdc42. Proc Natl Acad Sci U S A. 1998;95(26):15362–7.
Rohatgi R, Ma L, Miki H, Lopez M, Kirchhausen T, Takenawa T, Kirschner MW. The interaction between N-WASP and the Arp2/3 complex links Cdc42-dependent signals to actin assembly. Cell. 1999;97(2):221–31.
Stender S, Friebel A, Linder S, Rohde M, Mirold S, Hardt WD. Identification of SopE2 from Salmonella typhimurium, a conserved guanine nucleotide exchange factor for Cdc42 of the host cell. Mol Microbiol. 2000;36(6):1206–21.
Bhattacharyya S, Borthakur A, Pant N, Dudeja PK, Tobacman JK. Bcl10 mediates LPS-induced activation of NF-kappaB and IL-8 in human intestinal epithelial cells. Am J Physiol Gastrointest Liver Physiol. 2007;293(2):G429–37.
Ruan WK, Wu YH, An J, Zheng SJ. Polymorphisms of chicken Toll-like receptors 4, 15, and 21 in different breeds. Poult Sci. 2012;91(10):2512–6.
MacKinnon KM, He H, Nerren JR, Swaggerty CL, Genovese KJ, Kogut MH. Expression profile of toll-like receptors within the gastrointestinal tract of 2-day-old Salmonella enteriditis-infected broiler chickens. Vet Microbiol. 2009;137(3-4):313–9.
Takeda K, Akira S. Toll-like receptors. Curr Protoc Immunol. 2015;109:12–4.
Manicassamy S, Pulendran B. Modulation of adaptive immunity with Toll-like receptors. Semin Immunol. 2009;21(4):185–93.
Higuchi M, Matsuo A, Shingai M, Shida K, Ishii A, Funami K, Suzuki Y, Oshiumi H, Matsumoto M, Seya T. Combinational recognition of bacterial lipoproteins and peptidoglycan by chicken Toll-like receptor 2 subfamily. Dev Comp Immunol. 2008;32(2):147–55.
Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012;9(4):357–9.
An J, Lai J, Lehman ML, Nelson CC. miRDeep*: an integrated application tool for miRNA identification from RNA sequencing data. Nucleic Acids Res. 2013;41(2):727–37.
Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26(1):139–40.
John B, Enright AJ, Aravin A, Tuschl T, Sander C, Marks DS. Human MicroRNA targets. PLoS Biol. 2004;2(11):e363.
Enright AJ, John B, Gaul U, Tuschl T, Sander C, Marks DS. MicroRNA targets in Drosophila. Genome Biol. 2003;5(1):R1.
Dennis GJ, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, Lempicki RA. DAVID: Database for Annotation, Visualization, and Integrated Discovery. Genome Biol. 2003;4(5):3.
Huang DW, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4(1):44–57.
Huang DW, Sherman BT, Lempicki RA. Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009;37(1):1–13.
Al-Shahrour F, Diaz-Uriarte R, Dopazo J. FatiGO: a web tool for finding significant associations of Gene Ontology terms with groups of genes. Bioinformatics. 2004;20(4):578–80.
Joslyn CA, Mniszewski SM, Fulmer A, Heaton G. The gene ontology categorizer. Bioinformatics. 2004;20 (Suppl 1):i169–77.
Smoot ME, Ono K, Ruscheinski J, Wang PL, Ideker T. Cytoscape 2.8: new features for data integration and network visualization. Bioinformatics. 2011;27(3):431–2.
This project was supported by the National High Technology Research and Development Program of China (863 Program) (2011AA100305), the National Natural Science Foundation of China (31172204), Shandong Modern Agricultural Industry & Technology System (SDAIT-11-02), Shandong Province Agricultural Seed Project, and China Scholarship Council (201508370050).
Availability of data and materials
The expression data has been deposited into the Gene Expression Omnibus with the accession number GSE59748.
GW analyzed data, qPCR validation and drafted the manuscript. YQ carried out animal inoculation and sample collection. XL assisted in data analysis and PCR validation. NY participated in the experiment design. GX assisted in animal experiment design and sample collection. LL participated in the development of the experimental design and assisted in data analysis. XL provided the concepts of the study, and revised the manuscript. All authors read, edited and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
All animal procedures were approved by Shandong Agricultural University Animal Care and Use Committee.
The distribution of clean reads in the genome. (DOCX 13 kb)
Density distribution of miRNAs across chromosomes. Note: Densities are shown as number of miRNAs per megabase of DNA. (TIFF 256 kb)
Immune-related target genes and interaction with miRNAs. (XLS 128 kb)
Enriched GO BP terms in each group. (XLS 100 kb)
Primer information for miRNAs validated by qRT-PCR. (DOCX 17 kb)
The potential regulation between miRNAs and target genes. (DOCX 15 kb)
The designed primers for genes validated by qRT-PCR. (DOC 43 kb)
About this article
Cite this article
Wu, G., Qi, Y., Liu, X. et al. Cecal MicroRNAome response to Salmonella enterica serovar Enteritidis infection in White Leghorn Layer. BMC Genomics 18, 77 (2017). https://doi.org/10.1186/s12864-016-3413-8
- Laying chicken
- Salmonella enterica serovar Enteritidis
- Next generation sequencing