Next Article in Journal
Characterization and Functional Divergence of a Novel DUF668 Gene Family in Rice Based on Comprehensive Expression Patterns
Next Article in Special Issue
Complete Range of the Universal mtDNA Gene Pool and High Genetic Diversity in the Thai Dog Population
Previous Article in Journal
Artificial Intelligence (AI) in Rare Diseases: Is the Future Brighter?
Previous Article in Special Issue
Evidence for Early European Neolithic Dog Dispersal: New Data on Southeastern European Subfossil Dogs from the Prehistoric and Antiquity Ages
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 10
Issue 12
10.3390/genes10120979
genes-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Transcriptome Analysis of the Cecal Tonsil of Jingxing Yellow Chickens Revealed the Mechanism of Differential Resistance to Salmonella

by
Fei Wang
1,†,
Jin Zhang
1,†,
Bo Zhu
1,
Jie Wang
1,
Qiao Wang
1,
Maiqing Zheng
1,
Jie Wen
1,2,
Qinghe Li
1 and
Guiping Zhao
1,3,*
1
Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing 100193, China
2
State Key Laboratory of Animal Nutrition, Beijing 100193, China
3
School of Life Science and Engineering, Foshan University, Foshan 528000, China
*
Author to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Genes 2019, 10(12), 979; https://doi.org/10.3390/genes10120979
Submission received: 9 October 2019 / Revised: 21 November 2019 / Accepted: 26 November 2019 / Published: 28 November 2019
(This article belongs to the Special Issue Domestication Genetics)
Browse Figures
Versions Notes

Abstract

:
Salmonella is one of the most common food-borne pathogens. It can be transmitted between chickens, as well as to people by contaminated poultry products. In our study, we distinguished chickens with different resistances mainly based on bacterial loads. We compared the cecal tonsil transcriptomes between the susceptible and resistant chickens after Salmonella infection, aiming to identify the crucial genes participating in the antibacterial activity in the cecal tonsil. A total of 3214 differentially expressed genes (DEGs), including 2092 upregulated and 1122 downregulated genes, were identified between the two groups (fold change ≥ 2.0, padj < 0.05). Many DEGs were mainly involved in the regulation of two biological processes: crosstalk between the cecal tonsil epithelium and pathogenic bacteria, such as focal adhesion, extracellular-matrix–receptor interaction, and regulation of the actin cytoskeleton and host immune response including the cytokine–receptor interaction. In particular, the challenged resistant birds exhibited strong activation of the intestinal immune network for IgA production, which perhaps contributed to the resistance to Salmonella infection. These findings give insight into the mRNA profile of the cecal tonsil between the two groups after initial Salmonella stimulation, which may extend the known complexity of molecular mechanisms in chicken immune response to Salmonella.

1. Introduction

Salmonella is an important cause of food-borne and zoonotic disease, which can colonize in chickens. It also poses a serious threat to people’s health via the consumption of contaminated meat and eggs [1]. Many measures have been taken to reduce salmonellosis during poultry production, including the improvement of the breeding environment, sterilization and vaccination [2]. However, cases of salmonellosis happen occasionally [3]. Therefore, new efficient and permanent strategies are necessary to control this disease. These strategies require an understanding of the interaction between chickens and bacteria.
Chickens are normally infected by Salmonella via the fecal–oral route. In young chickens, the innate immune mechanism is particularly crucial during infection with Salmonella, because of the immature adaptive immunity [4]. The intestinal epithelium is the first physical and immunological barrier against enteric pathogens, because the epithelial cell layer can reduce invasion and localization of Salmonella to the liver and spleen tissues [5]. Besides the physical barrier, Toll-like receptors (TLRs) on the cell membrane, especially TLR4, play an essential role in the recognition and defense against Salmonella infection [6]. The activation of the TLR signal pathway induces the secretion of pro-inflammatory cytokines to control bacterial growth [7]. Previous research has shown that the expression of genes related to the TLR4-MyD88-dependent pathway slightly increased on the first and third day after Salmonella infection in the caecum [8]. In addition to cytokines, IgA and IgG, secreted by the lymphocytes of the mucosa, participate in the host’s intestinal immune response, and the concentrations increased after the Salmonella invasion [9,10,11,12]. During bacterial–epithelial crosstalk, a series of changes occur in the gut, including intestinal flora and metabolism. Salmonella relies on the virulence factors, T3SS-1 and T3SS-2, to trigger intestinal inflammation in adult mice [13]. During invasion, Salmonella reduces the intestinal concentration of butyrate by the depletion of Clostridium [14]. Inflammation and butyrate decreases transform the metabolism of epithelial cells from mitochondrial β-oxidation of fatty acids to anaerobic glycolysis [15]. Gut morphology, such as villous height and crypt depth, can be affected by Salmonella stimulation, because Salmonella could lead to apoptosis of cells [9].
The cecal tonsil is a crucial gut-associated lymphoid tissue and plays a major role in controlling the entry of bacteria and other pathogens into the cecum [16]. Many studies have focused on the splenic transcriptome after infection with Salmonella [17,18] to reveal the host’s immune-related genes and pathways, because the spleen primarily participates in the recognition and clearance of bacteria. In their research, cytokine–cytokine receptor interaction pathways were significantly enriched after the Salmonella challenge. The genes involved in Forkhead box O (FoxO) and mitogen-activated protein kinase (MAPK) signaling pathways were identified as the potential markers related to host resistance against Salmonella. So far, little is known about the cecal tonsil transcriptome change after infection with Salmonella, especially between resistant and susceptible chickens. In this study, we compared the mRNA expression of the cecal tonsils between resistant and susceptible chickens to elucidate effective mechanisms of host resistance against Salmonella.

2. Materials and Methods

2.1. Ethics Statement

Ethical approval on animal care and experimental procedures were performed in accordance with the Animal Ethics Committee of the Institute of Animal Sciences, Chinese Academy of Agricultural Sciences (IAS-CAAS, Beijing, China).

2.2. Animals and Experimental Design

The Salmonella Typhimurium (ST, 21484 standard strain) was purchased from China Industrial Microbial Culture Preservation Center (Beijing, China). The bacteria were resuscitated overnight in Luria–Bertani (LB) broth (Amresco, Washington, DC, USA) at 37 °C in an orbital shaking incubator at 150 rpm. After recovery, ST was cultured for 12 h, and concentrated in a centrifuge. The final number of colony forming units (CFU) was determined by plating serial dilutions. Jingxing Yellow Chickens were obtained from the Changping Experimental Base of Institute of Animal Sciences (Beijing, China). All the chickens were checked for the presence of Salmonella by culturing faecal samples in buffered peptone water overnight with shaking at 150 rpm and spreading the samples on brilliant green agar (37 °C, 18–24 h) [8]. According to the results, the positive chickens were eliminated. A total of 146 1-day-old chicks were raised in separate cages at the experimental center of China Agricultural University (Beijing, China) with free access to feed and water. At 7 days of age, the chicks were orally inoculated with 1 mL culture containing 2.5 × 1010 CFU Salmonella Typhimurium. The birds’ blood samples, livers and cecal tonsils were collected at 3 days post infection, respectively. The cecal tonsils were collected and placed in an −80 °C freezer for short-term storage, the blood was allowed to clot and the serum was stored, and the livers were used for the later measurement of bacterial loads.
The lysozyme concentration in the serum was measured by Chicken lysozyme (LZM) ELISA Kit (Cusabio, Wuhan, China) [19]. In brief, standards and serum diluent with biotin-conjugated lysozyme were pipetted into the wells, where the lysozyme antibody had been pre-coated. After incubation, the wells were washed and avidin conjugated Horseradish Peroxidase (HRP) was added to the wells. Following incubation and wash, a substrate solution was added to the wells, and color developed. The intensity of the color was measured by a microplate reader. According to the standards, the lysozyme in the serum was calculated. For the determination of bacterial loads, we first took tissues of the same weight, then ground and diluted these tissues to different concentrations, and finally plated the same amount of serial dilutions onto a MacConkey agar medium (Thermo Fisher Scientific, Waltham, MA, USA) to determine the presence of ST and calculate the bacterial loads. Finally, according to the liver bacterial loads and clinical signs including diarrhea, drooping wings and dying, 14 chicks were identified as susceptible (severe clinical signs and liver loads >104 CFU) or resistant (slight clinical sighs and liver loads <104 CFU) birds, from which the cecal tonsils were selected for RNA sequencing (RNA-seq).

2.3. Total RNA Isolation, cDNA Library Construction, and Sequencing

The cecal tonsil RNA was extracted from the 14 chickens using a QIAGEN kit (Qiagen, Hilden, Germany). The quality and quantity of the total RNA were assessed by a 2100 Bioanalyzer and RNA 6000 Nano kit (Agilent, Santa Clara, CA, USA). For the mRNA library construction and deep sequencing, 3 µg total of RNA was prepared using the TruSeq RNA Sample Preparation Kit (Illumina, San Diego, CA, USA) to capture the coding transcriptome. After purification, the RNA was fragmented using divalent cations at 95 °C. The cleaved RNA fragments were reversely transcribed into first-strand cDNA using TruSeq RNA Library Preparation Kit, followed by second-strand cDNA synthesis. After cDNA fragments purification and adaptor ligation, sequencing was performed on the HiSeq X Ten platform (Illumina, San Diego, CA, USA).
The sequence data reported in this paper was deposited in the Genome Sequence Archive of BIG Data Center, Beijing Institute of Genomics (BIG), Chinese Academy of Sciences and are publicly accessible at http://bigd.big.ac.cn/gsa/s/ZzF2353M (accession number: CRA002119).

2.4. Analysis of Differentially Expressed Genes (DEGs) and Pathway Enrichment

The quality control of reads was analyzed by FastQC software (v0.10.1) [20]. Sequence adapters and low-quality reads were eliminated. The clean reads were mapped to the chicken reference genome (Gallus gallus 5.0) by Hisat2 [21]. The expression quantities of the mapped transcripts were calculated using Htseq [22]. Analysis of the DEGs were conducted by DESeq2 [23]. Genes with fold-change ≥2 and padj <0.5 were considered to be DEGs. To assess the variation between samples, principal component analysis (PCA) was conducted by gmodels (based on the gene transcripts) (https://CRAN.R-project.org/package=gmodels) and hierarchical clustering (based on DEGs) was conducted by pheatmap (https://CRAN.R-project.org/package=pheatmap), respectively. Volcano plots were performed by ggplot2 to provide an overview of the DEGs [24]. Based on the DEGs, the gene enrichment analysis was conducted by Gene Ontology (GO) functional enrichment analysis [25] and KOBAS [26].

2.5. Quantitative Real-time PCR

After reverse transcription, qPCR was carried out using a SYBR Fast qPCR Master Mix (KAPA, Wilmington, MA, USA). The qPCR amplification system was as follows: 3 µL of cDNA (10-fold dilution), 5 µL of 2× SYBR Master Mix, 0.2 µL of ROX, 0.25 µL of each primer, and adding water to 10 µL. The samples were amplified using the real-time PCR Detection System ABI 7500 (Applied Biosystems, Foster City, CA, USA). The qPCR cycle parameters were as follows: 95 °C for 3 min, 40 cycles of 95 °C for 3 s, and 60 °C for 34 s. Three independent replicates were used for each assay. The 2−ΔΔCT method was used to calculate the relative abundance of transcripts [27]. The correlation analysis was conducted between the relative expression calculated by the q-PCR and the fold-change by RNA-seq.

3. Results

3.1. Overview of the Immune Response between Susceptible and Resistant Chickens

To choose that candidate individuals of the chickens for RNA-seq, the liver bacterial load was measured by counting the number of colonies on the MacConkey agar. In addition, the lysozyme content in the serum was evaluated by ELISA. Seven individuals were selected in combination with the clinical symptoms in each group. The bacterial loads of the resistant chickens were significantly lower than those in the susceptible chickens, as seen in Figure 1A. However, as seen in Figure 1B, the lysozyme contents of the resistant chickens were higher than those of the susceptible chickens, which indicated that the resistant chickens cleared ST effectively.

3.2. Sequencing of Cecal Tonsil Transcriptomes

RNA-seq of the cecal tonsils yielded >40 Mb in the 14 samples. Around 93% of the clean reads had quality scores exceeding the Q30 value. The data demonstrated the reliability of the RNA-seq and could be used for data analysis. After eliminating the interference reads, the clean reads accounted for >93%, as shown in Supplementary Table S1. A total of 17,296 genes were detected, and 4864 genes were novel.

3.3. Differentially Expressed mRNAs Responding to ST Infection in Cecal Tonsil

The differentially expressed mRNAs between susceptible and resistant chickens were identified using DESeq2. As seen in Figure 2A, there were 3214 genes significantly differentially expressed after ST stimulation. As shown in Supplementary Table S2, 2092 genes were upregulated and 1122 were downregulated. PCA indicated that two groups of birds were distinctly clustered, as shown in Figure 2B. As seen in Figure 2C, hierarchical clustering (based on all DEGs) was consistent with PCA, and showed that there were more upregulated than downregulated genes.

3.4. Quantitative Real-Time PCR Validation

To verify the accuracy of the sequencing, qPCR was carried out. A total of 13 DEGs were randomly selected from RNA-seq. The specific primers shown in Supplementary Table S3 were designed using Oligo 6.0 software, and were subsequently synthesized by BGI (Beijing, China). As seen in Figure 3, the correlation coefficient between RNA-seq and qPCR was 0.8824, which indicated that the RNA-seq data was reliable for subsequent analysis.

3.5. Functional Enrichment Analysis of the DEGs

Potential function analysis of all the DEGs was performed using GO and KEGG enrichment. Some significantly enriched GO terms were mainly involved in cell communication and adhesion, integral components of membranes, response to stimulus and the immune system, as shown in Table 1. After KEGG pathway analysis, 20 pathways were significantly enriched, as seen in Figure 4 and Supplementary Table S4. Several of these pathways were related to metabolism, including sphingolipid, nicotinate and nicotinamide metabolism. Many of the pathways were involved in the immune response, including the intestinal immune network for IgA production and the cytokine–receptor interaction. In addition, four enriched pathways were related to the interaction and adhesion between intestinal microorganisms and epithelial cells, including focal adhesion, regulation of actin cytoskeleton, cell adhesion molecules (CAMs) and the extracellular matrix (ECM)–receptor interaction. The DEGs relating to the cytokine–receptor interaction and focal adhesion were examined further. The relative DEGs of each pathway are shown in the Supplementary Tables S5 and S6.

4. Discussion

An increasing number of studies focus on the transcriptome changes in chickens after Salmonella stimulation, especially the spleen transcriptome. However, there are few studies on the changes in the cecal tonsil transcriptome. In order to understand the function and transcriptome changes in the cecal tonsil after Salmonella infection, we compared the cecal tonsil expression profiles between susceptible and resistant chickens. Apart from clinical symptoms, the bacterial burden in the livers and the lysozyme content in the serum were taken into consideration to evaluate the susceptibility and resistance of the chickens. The cecal tonsil expression profile was first detected by RNA-seq. A total of 3217 DEGs were identified, some of which were related to IgA production including AICDA, TNFRSF13C, CD86, TNFSF13B, CD80, ICOS, CCR9 and CD28. The significantly changed pathways related to the interaction between the bacteria and the epithelia and the immune response were identified in the cecal tonsil.
Many significantly enriched pathways, including focal adhesion, sphingolipid metabolism, regulation of actin cytoskeleton, CAMs and the ECM–receptor interaction were related to the bacterial–epithelial crosstalk, which occurs during the bacterial invasion of the intestine. In this process, a necessary step is the adhesion of bacterial pathogens to the intestinal epithelium, which is the basis for successful colonization and ultimate production of disease [28]. Recent studies suggest that adhesion relies on glycoconjugates induced by indigenous intestinal microflora on the host epithelial cells [29]. During this time, type III secretion, which Salmonella possesses [30], mediates the transfection of the bacterial proteins into the host cell membrane [31]. These secreted soluble molecules activate epithelial cells by a series of signaling pathways including Ca2+ transportation, resulting in actin cytoskeleton rearrangement and the entry of the bacteria into the cells. Sphingolipids, produced by both the host and specific bacteria, are closely related to metabolic and inflammatory pathways in host cells. The cellular levels and distribution of sphingolipids differ between inflamed and noninflamed intestinal tissue. Sphingolipids from a symbiotic microbe regulate homeostasis of host intestinal natural killer T cells, resulting in inflammatory bowel disease [32,33].
Correspondingly, the host has specialized strategies to resist such invasion. In our research, the intestinal immune network for IgA production was enriched. IgA is the predominant antibody isotype produced at mucosal surfaces, and plays a critical role in the intestinal immune response and the prevention of tissue damage in inflammation [34,35,36]. When the invasion of pathogenic microorganisms is recognized by the intestinal immune system, IgA is induced and transported into the lumen [37]. It has been confirmed that pathogen-binding IgA regulates bacterial motility and protects the host from the pathogens [38]. Furthermore, IgA can also influence the bacterial invasion of epithelial cells by inhibiting the type III secretion system [39]. In our research, the genes crucial for IgA production are more highly expressed in resistant chickens, as shown in Figure 5, which means that the IgA signaling pathway was more active in this group.
Li et al. characterized the splenic transcriptomes of susceptible and resistant chickens [18]. Compared with his results, we found that 58 DEGs were the same, some of which were mainly involved in metabolism, as shown in Supplementary Table S7. Metabolism is essential for the host to produce an effective immune response against pathogenic microorganisms. Previous research has revealed that during the recognition of microbial ligands, macrophages initiate upregulation of glycolysis. The energy of glycolysis supports antimicrobial inflammation and secretion of cytokines [40,41]. It has been shown that a glucose-rich diet improves host survival rate in systemic fungal infection [42]. The spleen and the cecal tonsil are two critical immune organs that may have metabolic differences between susceptible and resistant chickens after Salmonella infection.

5. Conclusions

Our research was aimed at analyzing the transcriptome of cecal tonsils of susceptible and resistant chickens after Salmonella infection. A total of 3214 DEGs were identified between them. DEGs were mainly involved in two biological processes, the interaction and crosstalk between the cecal tonsil epithelium and the host immune response, and the immune response of the host. Importantly, the adhesion between Salmonella and the intestine is a crucial signaling pathway, which plays an important role in the gut pathological changes after Salmonella infection. In addition, there was stronger activation of the intestinal immune network for the IgA production signaling pathway in resistant chickens, which may contribute to the protection of chicken from Salmonella invasion.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4425/10/12/979/s1, Table S1: Sample grouping and sequencing data quality control (XLSX), Table S2: The screened DEGs (FC ≥ 2; padi < 0.05) in this study (XLSX), Table S3: The specific primers for qRT-PCR in this study (XLSX), Table S4: The enriched pathways (XLSX), Table S5: The DEGs related to focal adhesion in this study (XLSX), Table S6: The DEGs related to cytokine-cytokine receptor interaction in this study (XLSX), Table S7: The same DEGs and some of which mainly involved in metabolism.

Author Contributions

Conceptualization and methodology, F.W., J.Z., B.Z., J.W. (Jie Wang), and G.Z.; validation, Q.L. and G.Z.; formal analysis, F.W., J.Z. and B.Z.; investigation, Q.L. and G.Z.; resources, J.W. (Jie Wen) and G.Z.; data curation, F.W., J.Z. and B.Z.; writing—original draft preparation, F.W. and J.Z.; writing—review & editing, Q.L. and G.Z.; visualization, F.W. and J.Z.; supervision, Q.L. and G.Z.; project administration, J.W. (Jie Wen) and G.Z.; funding acquisition, Q.W. and M.Z.

Funding

This research was supported by grants from National Natural Science Foundation of China (No. 31572393), the Natural Science Foundation of Beijing (6182032), National Nonprofit Institute Research Grant (2017ywf-zd-5) and the Agricultural Science and Technology Innovation Program (ASTIPIAS04) of Chinese Academy of Agricultural Sciences.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

References

  1. Barrow, P.A.; Jones, M.A.; Smith, A.L.; Wigley, P. The long view: Salmonella—The last forty years. Avian Pathol. 2012, 41, 413–420. [Google Scholar] [CrossRef] [PubMed]
  2. Zhang-Barber, L.; Turner, A.K.; Barrow, P.A. Vaccination for control of Salmonella in poultry. Vaccine 1999, 17, 2538–2545. [Google Scholar] [CrossRef]
  3. Snow, L.C.; Davies, R.H.; Christiansen, K.H.; Carrique-Mas, J.J.; Wales, A.D.; O’Connor, J.L.; Cook, A.J.C.; Evans, S.J. Survey of the prevalence of Salmonella species on commercial laying farms in the United Kingdom. Vet. Rec. 2007, 161, 471. [Google Scholar] [CrossRef] [PubMed]
  4. Bar-Shira, E.; Sklan, D.; Friedman, A. Establishment of immune competence in the avian GALT during the immediate post-hatch period. Dev. Comp. Immunol. 2003, 27, 147–157. [Google Scholar] [CrossRef]
  5. Bryan, T.; Nicholas, P.; Caitlyn, D.; Rafaela, P.; Mary, M.; Hassan, H.M.; Koci, M.D. Poultry Body Temperature Contributes to Invasion Control through Reduced Expression of Salmonella Pathogenicity Island 1 Genes in Salmonella enterica Serovars Typhimurium and Enteritidis. Appl. Environ. Microbiol. 2015, 81, 8192. [Google Scholar]
  6. Monack, D.M.; Anne, M.; Stanley, F. Persistent bacterial infections: The interface of the pathogen and the host immune system. Nat. Rev. Microbiol. 2004, 2, 747–765. [Google Scholar] [CrossRef]
  7. Weiss, D.S.; Raupach, B.; Takeda, K.; Akira, S.; Zychlinsky, A. Toll-like receptors are temporally involved in host defense. J. Immunol. 2004, 172, 4463–4469. [Google Scholar] [CrossRef]
  8. Li, P.; Xia, P.; Wen, J.; Zheng, M.; Chen, J.; Zhao, J.; Jiang, R.; Liu, R.; Zhao, G. Up-regulation of the MyD88-dependent pathway of TLR signaling in spleen and caecum of young chickens infected with Salmonella serovar Pullorum. Vet. Microbiol. 2010, 143, 346–351. [Google Scholar] [CrossRef]
  9. Qiu, J.W.; Xiao, C.Z.; Tian, W.; Zhang, T.Y. Effect of dietary oridonin supplementation on growth performance, gut health, and immune response of broilers infected with Salmonella pullorum. Irish Vet. J. 2018, 71, 16. [Google Scholar]
  10. Wieland, W.H.; Diego, O.; Aart, L.; Parmentier, H.K.; Verstegen, M.W.A.; Arjen, S. A functional polymeric immunoglobulin receptor in chicken (Gallus gallus) indicates ancient role of secretory IgA in mucosal immunity. Biochem. J. 2004, 380, 669–676. [Google Scholar] [CrossRef]
  11. Webb, T.; Goodman, H.C. The structure and function of immunoglobulins. Mod. Trends Immunol. 1967, 2, 151–187. [Google Scholar] [PubMed]
  12. Wieland, W.H.; Kruijt, L.; Jansma, A. Successive immunoglobulin and cytokine expression in the small intestine of juvenile chicken. Dev. Comp. Immunol. 2010, 34, 1254–1262. [Google Scholar]
  13. Bryan, C.; Yuling, L.; David, O.; Vallance, B.A.; B Brett, F. Salmonella enterica serovar Typhimurium pathogenicity island 2 is necessary for complete virulence in a mouse model of infectious enterocolitis. Infect. Immun. 2005, 73, 3219–3227. [Google Scholar]
  14. Rivera-Chávez, F.; Zhang, L.F.; Faber, F.; Lopez, C.A.; Byndloss, M.X.; Olsan, E.E.; Xu, G.; Velazquez, E.M.; Lebrilla, C.B.; Winter, S.E. Depletion of Butyrate-Producing Clostridia from the Gut Microbiota Drives an Aerobic Luminal Expansion of Salmonella. Cell Host Microbe 2016, 19, 443–454. [Google Scholar] [CrossRef] [PubMed]
  15. Byndloss, M.X.; Olsan, E.E.; Rivera-Chávez, F.; Tiffany, C.R.; Cevallos, S.A.; Lokken, K.L.; Torres, T.P.; Byndloss, A.J.; Faber, F.; Gao, Y. Microbiota-activated PPAR-γ signaling inhibits dysbiotic Enterobacteriaceae expansion. Science 2017, 357, 570–575. [Google Scholar] [CrossRef] [PubMed]
  16. Setta, A.M.; Barrow, P.A.; Kaiser, P.; Jones, M.A. Early immune dynamics following infection with Salmonella enterica serovars Enteritidis, Infantis, Pullorum and Gallinarum: Cytokine and chemokine gene expression profile and cellular changes of chicken cecal tonsils. Comp. Immunol. Microbiol. Infect. Dis. 2012, 35, 397–410. [Google Scholar] [CrossRef] [PubMed]
  17. Zhou, H.; Lamont, S.J. Global gene expression profile after Salmonella enterica Serovar enteritidis challenge in two F8 advanced intercross chicken lines. Cytogenet. Genome Res. 2007, 117, 131–138. [Google Scholar] [CrossRef]
  18. Peng, L.; Fan, W.; Everaert, N.; Liu, R.; Li, Q.; Zheng, M.; Cui, H.; Zhao, G.; Jie, W. Messenger RNA Sequencing and Pathway Analysis Provide Novel Insights into the Susceptibility toSalmonella enteritidisInfection in Chickens. Front. Genet. 2018, 9, 256. [Google Scholar]
  19. Lollike, K.; Kjeldsen, L.; Sengelov, H.; Borregaard, N. Purification of lysozyme from human neutrophils, and development of an ELISA for quantification in cells and plasma. Leukemia 1995, 9, 206–209. [Google Scholar]
  20. Chen, S.; Zhou, Y.; Chen, Y.; Gu, J. fastp: An ultra-fast all-in-one FASTQ preprocessor. Bioinformatics 2018, 34, i884–i890. [Google Scholar] [CrossRef]
  21. Kim, D.; Langmead, B.; Salzberg, S.L. HISAT: A fast spliced aligner with low memory requirements. Nat. Methods 2015, 12, 357–360. [Google Scholar] [CrossRef] [PubMed]
  22. Anders, S.; Pyl, P.T.; Huber, W. HTSeq--A Python framework to work with high-throughput sequencing data. Bioinformatics 2015, 31, 166–169. [Google Scholar] [CrossRef] [PubMed]
  23. Love, M.I.; Huber, W.; Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014, 15, 550. [Google Scholar] [CrossRef] [PubMed]
  24. Ginestet, C. ggplot2: Elegant Graphics for Data Analysis. J. R. Stat. Soc. A Stat. 2011, 174, 245–246. [Google Scholar] [CrossRef]
  25. Young, M.D.; Wakefield, M.J.; Smyth, G.K.; Oshlack, A. Gene ontology analysis for RNA-seq: Accounting for selection bias. Genome Biol. 2010, 11, R14. [Google Scholar] [CrossRef]
  26. Xie, C.; Mao, X.; Huang, J.; Ding, Y.; Wu, J.; Dong, S.; Kong, L.; Gao, G.; Li, C.Y.; Wei, L. KOBAS 2.0: A web server for annotation and identification of enriched pathways and diseases. Nucleic Acids Res. 2011, 39, W316–W322. [Google Scholar] [CrossRef]
  27. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef]
  28. Lu, L.; Walker, W. Pathologic and physiologic interactions of bacteria with the gastrointestinal epithelium. Am. J. Clin. Nutr. 2001, 73, 1124S. [Google Scholar] [CrossRef]
  29. Schauer, D.B. Indigenous microflora: Paving the way for pathogens? Curr. Biol. 1997, 7, R75–R77. [Google Scholar] [CrossRef]
  30. Goosney, D.L.; Knoechel, D.G.; Finlay, B.B. Enteropathogenic, E. coli, Salmonella, and Shigella: Masters of host cell cytoskeletal exploitation. Emerg. Infect. Dis. 1999, 5, 216–223. [Google Scholar] [CrossRef]
  31. Gauthier, A.; Finlay, B.B. Protein translocation: Delivering virulence into the host cell. Curr. Biol. 1998, 8, R768–R770. [Google Scholar] [CrossRef]
  32. An, D.; Oh, S.; Olszak, T.; Neves, J.; Avci, F.; Erturk-Hasdemir, D.; Lu, X.; Zeissig, S.; Blumberg, R.; Kasper, D. Sphingolipids from a Symbiotic Microbe Regulate Homeostasis of Host Intestinal Natural Killer T Cells. Cell 2014, 156, 123–133. [Google Scholar] [CrossRef] [PubMed]
  33. Abdel, H.L.; Di, V.C.; Riboni, L. Fostering Inflammatory Bowel Disease: Sphingolipid Strategies to Join Forces. Mediat. Inflamm. 2016, 2016, 3827684. [Google Scholar] [CrossRef] [PubMed]
  34. Oliver, P. New concepts in the generation and functions of IgA. Nat. Rev. Immunol. 2012, 12, 821–832. [Google Scholar]
  35. Slack, E.; Balmer, M.L.; Fritz, J.H.; Hapfelmeier, S. Functional flexibility of intestinal IgA-broadening the fine line. Front. Immunol. 2012, 3, 100. [Google Scholar] [CrossRef] [PubMed]
  36. Mkaddem, S.B.; Christou, I.; Rossato, E.; Berthelot, L.; Lehuen, A.; Monteiro, R.C. IgA, IgA Receptors, and Their Anti-inflammatory Properties. Curr. Top. Microbiol. Immunol. 2014, 382, 221. [Google Scholar]
  37. Palm, N.; Dezoete, M.; Cullen, T.; Barry, N.; Stefanowski, J.; Hao, L.; Degnan, P.; Hu, J.; Peter, I.; Zhang, W. Immunoglobulin A Coating Identifies Colitogenic Bacteria in Inflammatory Bowel Disease. Cell 2014, 158, 1000–1010. [Google Scholar] [CrossRef]
  38. Forbes, S.J.; Marisa, E.; Mantis, N.J. Inhibition of Salmonella enterica serovar typhimurium motility and entry into epithelial cells by a protective antilipopolysaccharide monoclonal immunoglobulin A antibody. Infect. Immun. 2008, 76, 4137–4144. [Google Scholar] [CrossRef]
  39. Forbes, S.J.; Tia, B.; Mccarthy, E.A.; Blaise, C.; Mantis, N.J. Transient suppression of Shigella flexneri type 3 secretion by a protective O-antigen-specific monoclonal IgA. Mbio 2011, 2, 00042-11. [Google Scholar] [CrossRef] [Green Version]
  40. Freemerman, A.J.; Johnson, A.R.; Sacks, G.N.; Milner, J.J.; Kirk, E.L.; Troester, M.A.; Macintyre, A.N.; Pankuri, G.H.; Rathmell, J.C.; Liza, M. Metabolic reprogramming of macrophages: Glucose transporter 1 (GLUT1)-mediated glucose metabolism drives a proinflammatory phenotype. J. Biol. Chem. 2014, 289, 7884. [Google Scholar] [CrossRef] [Green Version]
  41. Van den Bossche, J.; O’Neill, L.A.; Menon, D. Macrophage Immunometabolism: Where Are We (Going)? Trends Immunol. 2017, 38, 395. [Google Scholar] [CrossRef] [PubMed]
  42. Tucey, T.M.; Verma, J.; Harrison, P.F.; Snelgrove, S.L.; Lo, T.L.; Scherer, A.K.; Barugahare, A.A.; Powell, D.R.; Wheeler, R.T.; Hickey, M.J. Glucose Homeostasis Is Important for Immune Cell Viability during Candida Challenge and Host Survival of Systemic Fungal Infection. Cell Metab. 2018, 27, 988. [Google Scholar] [CrossRef] [PubMed]
Genes 10 00979 g001 550
Figure 1. Bacterial burden and lysozyme content in chickens. (A) The liver was plated onto the MacConkey agar medium and the bacterial burden was determined by counting. (B) Lysozyme concentration in the serum was measured by ELISA. Data are shown as means ± SD. Data with asterisks were statistically significant (p ≤ 0.05).
Figure 1. Bacterial burden and lysozyme content in chickens. (A) The liver was plated onto the MacConkey agar medium and the bacterial burden was determined by counting. (B) Lysozyme concentration in the serum was measured by ELISA. Data are shown as means ± SD. Data with asterisks were statistically significant (p ≤ 0.05).
Genes 10 00979 g001
Genes 10 00979 g002 550
Figure 2. RNA-seq analysis and identification of genes differentially expressed between susceptible and resistant chickens. (A) Volcano plot showing DEGs. The red dots represent significantly upregulated genes (log2 FC ≥1 and padj <0.05); green dots represent significantly downregulated genes; and black dots represent genes with no significant change. (B) PCA of the RNA-seq data based on all the identified genes. The susceptible and resistant chickens were distinctly clustered. (C) The heat map of the RNA-seq data based on DEGs. S, susceptible chickens; R, resistant chickens.
Figure 2. RNA-seq analysis and identification of genes differentially expressed between susceptible and resistant chickens. (A) Volcano plot showing DEGs. The red dots represent significantly upregulated genes (log2 FC ≥1 and padj <0.05); green dots represent significantly downregulated genes; and black dots represent genes with no significant change. (B) PCA of the RNA-seq data based on all the identified genes. The susceptible and resistant chickens were distinctly clustered. (C) The heat map of the RNA-seq data based on DEGs. S, susceptible chickens; R, resistant chickens.
Genes 10 00979 g002
Genes 10 00979 g003 550
Figure 3. Correlation analysis of the relative expression levels of 13 DEGs between the RNA-seq and qPCR.
Figure 3. Correlation analysis of the relative expression levels of 13 DEGs between the RNA-seq and qPCR.
Genes 10 00979 g003
Genes 10 00979 g004 550
Figure 4. KEGG pathway analysis of DEGs. The y axis represents the significantly enriched KEGG terms based on the DEGs. The x axis represents rich factors (rich factor = number of DEGs enriched in each term/number of all genes in each term). Color represents significance, and size of the bubble represents the number of DEG enriched in the pathway (p-Value: before correction).
Figure 4. KEGG pathway analysis of DEGs. The y axis represents the significantly enriched KEGG terms based on the DEGs. The x axis represents rich factors (rich factor = number of DEGs enriched in each term/number of all genes in each term). Color represents significance, and size of the bubble represents the number of DEG enriched in the pathway (p-Value: before correction).
Genes 10 00979 g004
Genes 10 00979 g005 550
Figure 5. IgA production signaling pathway. The DEGs identified in this research are highlighted in purple.
Figure 5. IgA production signaling pathway. The DEGs identified in this research are highlighted in purple.
Genes 10 00979 g005
Table
Table 1. Significant enriched biological processes by gene ontology analysis based on the DEGs.
Table 1. Significant enriched biological processes by gene ontology analysis based on the DEGs.
TermDescriptionCountCorrected p-Value
GO:0023052signaling2370.002
GO:0007154cell communication2370.002
GO:0050896response to stimulus3040.002
GO:0007165signal transduction2190.003
GO:0002376immune system process980.004
GO:0007155cell adhesion720.007
GO:0022610biological adhesion720.007
GO:0016020membrane3550.007
GO:0016021Integral component of membrane2360.009
GO:0044425membrane part2760.009
GO:0098602single organism cell adhesion460.100
GO:0031224intrinsic component of membrane2380.012
GO:0016337single organismal cell–cell adhesion430.023
GO:0051716cellular response to stimulus2480.028
GO:0006955immune response520.032
GO:0005886plasma membrane1650.032
GO:0071944cell periphery1680.043
GO:0032501multicellular organismal process2360.049
All DEGs between susceptible and resistant chickens were used to identify enriched biological functions (p < 0.05).

Share and Cite

MDPI and ACS Style

Wang, F.; Zhang, J.; Zhu, B.; Wang, J.; Wang, Q.; Zheng, M.; Wen, J.; Li, Q.; Zhao, G. Transcriptome Analysis of the Cecal Tonsil of Jingxing Yellow Chickens Revealed the Mechanism of Differential Resistance to Salmonella. Genes 2019, 10, 979. https://doi.org/10.3390/genes10120979

AMA Style

Wang F, Zhang J, Zhu B, Wang J, Wang Q, Zheng M, Wen J, Li Q, Zhao G. Transcriptome Analysis of the Cecal Tonsil of Jingxing Yellow Chickens Revealed the Mechanism of Differential Resistance to Salmonella. Genes. 2019; 10(12):979. https://doi.org/10.3390/genes10120979

Chicago/Turabian Style

Wang, Fei, Jin Zhang, Bo Zhu, Jie Wang, Qiao Wang, Maiqing Zheng, Jie Wen, Qinghe Li, and Guiping Zhao. 2019. "Transcriptome Analysis of the Cecal Tonsil of Jingxing Yellow Chickens Revealed the Mechanism of Differential Resistance to Salmonella" Genes 10, no. 12: 979. https://doi.org/10.3390/genes10120979

APA Style

Wang, F., Zhang, J., Zhu, B., Wang, J., Wang, Q., Zheng, M., Wen, J., Li, Q., & Zhao, G. (2019). Transcriptome Analysis of the Cecal Tonsil of Jingxing Yellow Chickens Revealed the Mechanism of Differential Resistance to Salmonella. Genes, 10(12), 979. https://doi.org/10.3390/genes10120979

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop