Explore

Communities in English

Advertise on Engormix

Host transcriptome response to heat stress and Eimeria maxima infection in meat-type chickens

Published: May 2, 2024
By: Ahmed F. A. Ghareeb 1, James C. Foutz 1, Gustavo H. Schneiders 1, Jennifer N. Richter 1, Marie C. Milfort 1, Alberta L. Fuller 1, Romdhane Rekaya 2, Samuel E. Aggrey 1 / 1 NutriGenomics Laboratory, Department of Poultry Science, University of Georgia, Athens, Georgia, USA; 2 Department of Animal and Dairy Science, University of Georgia, Athens, Georgia, USA.
Introduction
Heat stress (HS) is one of the most challenging environmental stressors despite the modern climate control equipment in broiler chickens’ houses. Broiler flocks may be seasonally exposed to HS that varies in intensity according to the relative humidity of the region [1]. Multiple studies have identified the various adverse effects of acute or chronic HS on chickens, such as a severe reduction in feed intake and growth [2–4], systemic alkalosis [5], immune suppression [6, 7], and increased intestinal permeability [8]. During HS the cellular production of reactive oxygen species (ROS) is dramatically increased causing oxidative stress that impacts cellular functions and even viability [5, 9, 10].
Eimeria (E.) maxima is one of the most common apicomplexan species infecting chickens causing coccidiosis. Eimeria maxima preferably invades the upper part of the jejunum and whole ileum and produces a huge number of merogonic and gametogenic stages, damaging the intestinal mucosal lining and resulting in severe inflammation, hemorrhage, and diarrhea. Consequently, areas of erosion and ulceration would be developed exposing the intestine to secondary pathogenic invasion, such as Clostridium perfringens causing necrotic enteritis. E. maxima spp infection is associated with significant production losses resulting from weight loss, poor feed conversion, high morbidity, and mortality [11–13]. Exposing E. maxima-infected chickens to HS-suppressed E. maxima gametogony. Broiler chickens exposed to a combination of HS and E. maxima infection exhibited higher values for the apparent ileal digestibility and less enterocytic damage than the E. maxima-infected chickens raised under thermoneutral condition [14–16].
Identifying the differentially expressed genes (DEGs) and integrating them into known molecular pathways and biological functions could lead to the mechanisms by which cells respond to different conditions [17]. This study provides an unprecedented analysis of the ileal transcriptome under HS, E. maxima infection, and their combination to elucidate the molecular mechanisms controlling the ileal tissue response to each stressor and key molecular functions by which HS limits the E. maxima life cycle.
Materials and methods
This study was implemented under the Animal Use Proposal (AUP) A2015 04–005 approved by the University of Georgia Animal Care and Use Committee (IACUC).
Experimental design
Two hundred and forty 14-day-old Ross 708 male broilers were randomly allocated into four treatment groups: Thermoneutral control (TNc), thermoneutral infected (TNi), heat stress control (HSc), and heat stress infected (HSi). The thermoneutral and HS groups were raised at 20°C and 35°C, respectively. The infected groups received via oral gavage 2x105 E. maxima sporulated oocysts/bird, while the control groups were mock-infected with water. E. maxima infective doses were acquired from a single oocyst cloning performed as described by Schneider et al. [14]. The infected groups showed obvious coccidiosis clinical signs started at 4 day-post-infection (dpi), and the beak of growth depression was at 6 dpi [16]. Only HSi group had two recorded mortalities, one was found dead at 3dpi and the other was humanly euthanized at 6 dpi. To confirm the infection, oocyst shedding was detected in the infected groups (TNi and HSi) groups, but not in control groups (TNc and HSc). All the chickens were housed in wired-floor cages and ad libitum supplied with water and a non-medicated standard grower diet. The human end point was considered to euthanize the bird upon showing inability to reach feed and water, sever respiratory distress, or severely emaciation. Any of these factors led to humane euthanization by cervical dislocation.
Five chickens were randomly selected from each treatment group at 6 (dpi) and euthanized by cervical dislocation to sample ~1 cm of ileum tissue from Meckel’s diverticulum. Ileum samples were snap-frozen in liquid nitrogen and then, stored at -80°C for later RNA extraction and sequencing using the NGS Illumina sequencing platform.
Nucleic acid extraction
The trizol-chloroform method was used to extract the mRNA. Briefly, 100 mg of frozen ileum tissue was homogenized in 1 mL of trizol using Benchmark BeadBlaster Microtube Homogenizer® (OVERSTOCK Lab Equipment, NH, USA), and then phase separated by mixing 0.2 mL of chloroform, 3 minutes incubation at room temperature (RT), and centrifugation (13500 xg, 4°C, for 15 min). About 550 μL of the aqueous phase was collected into a new tube containing 0.5 mL isopropanol to precipitate the RNA, mixing by vortex and 10 min incubation at RT. The mixture was centrifugated at 13500 xg, 4°C for 10 min, followed by supernatant disposal. The pelleted RNA was washed by vortex mixing with 1 mL of 75% ethanol, then reprecipitated by centrifugation at 13500 xg, 4°C for 5 min, followed by discarding the supernatants. The tubes containing the pelleted mRNA were left for 10 min at RT for drying, followed by resuspending the mRNA pellet in 100μL of nuclease-free distilled water and then incubating at 55–60˚C for 10 min before being stored at -80°C.
cDNA library construction and sequencing
The extracted RNA was purified using the RNeasy Mini Kit (Qiagen, Hilden, Germany) following the manufacturer’s protocol. The concentration of the purified RNA was determined using a NanoDrop 2000 Spectrophotometer (Thermo Fischer Scientific, DE, USA), showing OD260/280 ratios for all samples > 1.9. The RNA integrity number (RIN) of all samples was ≥ 9. cDNA libraries were prepared with 4 μg total RNA using the TruSeq RNA Sample Preparation Kit to acquire cDNA fragments with an average size of 229 bp, or 355 bp inclusive of the adapter sequences. The prepared cDNA libraries were sequenced by 150 bp paired-end read chemistry using the Illumina HiSeq 2000 system.
Sequence quality control, mapping, and annotation
The quality of raw reads was assessed by FastQC, and the trimming of the low-quality reads was conducted using Trimmomatic v.0.36 [18]. The reads were aligned to the reference genome of chicken (Gallus gallus 5.0.90, Ensembl) using STAR aligner v.2.5.2b [19]. Hit counts were accounted for using the feature Counts v.1.5.2 package [20], considering only the unique reads for downstream analysis. Before the mapping step, the quality of reads was inspected -and trimmed when required- for the adapter sequence using Flexbar version 2.4 [21].
Sequence reads mapped to the transcriptome were reported according to their genome-equivalent coordinates. The genome coverage data were acquired using the Bamtools version 2.5.1 ’stats’ command [22]. Cufflinks version 2.2.1 [23] was utilized to detect the expressed genes and transcripts by contrasting the Tophat2 read-mapping to the Ensembl gene models for the Gallus gallus 5.0.90 assembly of the chicken genome, including the known- and predicted genes. Afterward, the fragments per kilobase of exon per million mapped fragments (FPKM) expression values for each sample were collected using Cuffdiff version 2.1.1 [23] both at the Gene (G) and gene Isoform (I) levels.
Gene differential expression, gene ontology, and Kyoto encyclopedia of genes and genomes pathway analysis
DESeq2 [24] was used to ascertain DEGs. The Wald test within the DESeq2 package was used to acquire the p-value and log2 fold-change of the DEGs. The p-value then was adjusted for multiple testing with the Benjamini-Hochberg method [25]. The DEGs with false discovery rate adjusted (adj. p-value) ≤0.05 and absolute log2 fold change > 1.2 were considered as significantly differentially expressed. The gene ontogeny (GO) analysis was conducted and functionally visualized using the Cytoscape v3.7.2 software platform provided with the ClueGO Plugin v2.5.4 application [26, 27]. The pathways analysis was performed using the Kyoto Encyclopedia of Genes and Genomes (KEGG) terms [28, 29]. The GO pathways/networks were used to classify DEGs into cellular, molecular, and biological functions.
Hi-Seq data confirmation
The Hi-Seq expression data was confirmed using RT-qPCR. Genes were randomly selected from the common DEGs shared between different group comparisons (n = 10 genes for each comparison). The relative mRNA expression values were calculated according to 2-ΔΔCt method [30], and then compared with the fold change values yielded from the Hi-Seq data (S1A–S4B Figs).
Results
The sequencing generated an average of 38,692,865 reads per sample, with an average of 32,541,154 reads uniquely mapped. The average total mapped read and unique mapped reads were 88 and 84%, respectively. The number of DEGs resulting from each pairwise comparison is presented in Fig 1. The top 100 DEGs (50 downregulated and 50 upregulated) of each pairwise comparison are provided in the supplementary Excel sheets (S1–S4 Tables).
Host transcriptome response to heat stress and Eimeria maxima infection in meat-type chickens - Image 1
Fig 1. Number of DEGs of each pairwise comparison at 6 dpi of chickens infected with Eimeria maxima and their uninfected controls that are raised in a thermoneutral or heat stress environment: TNc = thermoneutral control, TNi = thermoneutral infected, HSc = heat stress control, HSi = heat stress infected.
The effect of heat stress on the ileal transcriptome
A total of 413 DEGs were identified when the HSc group was compared with the TNc group (Fig 1). The genes C-X-C motif chemokine ligand 13 and ligand 13-like 2 (CXCL13, CXCL13L2); interleukins 22, 8-like 1, and 1 receptor type 2 (IL-22, IL-8L1, and IL-1R2); TNF receptor superfamily member 13B (TNFRSF13B); heat shock proteins family A (HSP70) member 8 and member 4 like (HSPA8 and HSPA4L), family H (Hsp110) member 1 (HSPH1), and 90 alpha family class A member 1 (HSP90AA1) were among the most downregulated (Fig 2A) in HSc compared to the TNc group. The highest upregulated genes in the HSc group compared with the TNc group included amidohydrolase domain containing 1 (AMDHD1), histidine ammonia-lyase (HAL), cholecystokinin A receptor (CCKAR), cholinergic receptor nicotinic alpha 7 subunit (CHRNA7), polypeptide N-acetylgalactosaminyltransferase 13 (GALNT13), glutathione S-transferase alpha 2 (GSTA2), and lipase G, endothelial type (LIPG) (Fig 2A).
Host transcriptome response to heat stress and Eimeria maxima infection in meat-type chickens - Image 2
Host transcriptome response to heat stress and Eimeria maxima infection in meat-type chickens - Image 3
Host transcriptome response to heat stress and Eimeria maxima infection in meat-type chickens - Image 4
Fig 2.
A. Heatmap of the top ~60 DEGs of TNc vs. HSc, calculated as Log 2 relative hit counts of the comparison between the chickens exposed to HS (HSc) and their thermoneutral control (TNc) at 6 day-post-treatment. B. The top KEGG pathways of TNc vs. HSc, based on the pathway term significance (α≤0.1) of the comparison between the chickens exposed to HS (HSc) and their thermoneutral control (TNc) at 6 day-post-treatment. The green squares and arrowheads depict the downregulated pathways and genes, respectively. The red circles and arrowheads depict the upregulated pathways and genes, respectively. The bigger size squares and circles denote the lower p-value. An increase in the color transparency (lighter color) of the squares and circles indicates an increase in the percentage of the included genes with expression that counter the regulation direction of the pathway. C. The top significant KEGG pathways of TNc vs. HSc, the comparison between the chickens exposed to HS (HSc) and their thermoneutral control (TNc) at 6 day-post-treatment. D. The top significant Gene Ontology terms of TNc vs. HSc, the comparison between the chickens exposed to HS (HSc) and their thermoneutral control (TNc) at 6 day-post-treatment.
The KEGG pathway network analysis depicted in Fig 2B and 2C shows downregulation of the “cytokine-cytokine receptor interaction”, “intestinal immune network for IgA production”, “Toll- and nod-like receptor signaling”, and “protein processing in the endoplasmic reticulum (ER)” pathways, whereas the pathways involved in metabolism and degradation of many amino acids such as histidine, glycine, tyrosine, valine, and leucine were upregulated. The GO terms are exhibited in Fig 2D. The Cellular compartments showed downregulation of the terms related to “endoplasmic reticulum compartments”, “basolateral plasma membrane”, and “immunological synapse”, and upregulation of the “apical plasma membrane” (Fig 2D). Consistently, the molecular function and biological process terms related to immune response and inflammation, and protein synthesis in the endoplasmic reticulum were downregulated in the HSc group compared to the TNc group, while the terms involved in the amino acids and organic acids metabolic processing were upregulated (Fig 2D).
The effect of E. maxima infection on the ileal transcriptome
A total of 3,377 DEGs were identified when the TNi group was compared with the TNc group (Fig 1). The most downregulated genes in the TNi group compared to the TNc group included alcohol dehydrogenase 1C (class I); gamma polypeptide (ADH1C); aldehyde dehydrogenase 1 family member A1 (ALDH1A1), cytochrome P450 family 1 subfamily A members 1 and 2, family 2 subfamily AC polypeptide 1, family 2 subfamily C polypeptide 23, family 2 subfamily D member 6 (CYP1A1, CYP1A2, CYP2AC1, CYP2C23b, and CYP2D6); fatty acids binding proteins 1 and 6 (FABP1 and FABP6); and phosphoenolpyruvate carboxykinase 1 (PCK1) (Fig 3A). Among the top upregulated genes were identified interleukins (IL-1β and IL12B); interleukin 1 receptor type 2 (IL1R2); C-C motif chemokine ligand 4 and 26 (CCL4 and CCAH221); nitric oxide synthase 2 (NOS2); colony-stimulating factor 3 (CSF3), integrin subunit alpha 2 (ITGA2); and secreted phosphoprotein 1 (SPP1) (Fig 3A).
Host transcriptome response to heat stress and Eimeria maxima infection in meat-type chickens - Image 5
Host transcriptome response to heat stress and Eimeria maxima infection in meat-type chickens - Image 6
Host transcriptome response to heat stress and Eimeria maxima infection in meat-type chickens - Image 7
Fig 3.
A. Heatmap of the top ~60 DEGs of TNc vs. TNi, as Log 2 relative hit counts of the comparison between the chickens infected with Eimeria maxima (TNi) and their uninfected thermoneutral control (TNc) at 6 day-post-treatment. B. The top KEGG pathways of TNc vs. TNi, based on the pathway term significance (α≤0.1) of the comparison between the chickens infected with Eimeria maxima (TNi) and their uninfected thermoneutral control (TNc) at 6 day-post-treatment. The green squares and arrowheads depict the downregulated pathways and genes, respectively. The red circles and arrowheads depict the upregulated pathways and genes, respectively. The bigger size squares and circles denote the lower p-value. An increase in the color transparency (lighter color) of the squares and circles indicates an increase in the percentage of the included genes with expression that counter the regulation direction of the pathway. C. The top significant KEGG pathways of TNc vs. TNi, the comparison between the chickens infected with Eimeria maxima (TNi) and their uninfected thermoneutral control (TNc) at 6 day-post-treatment. D. The top significant Gene Ontology terms of TNc vs. TNi, the comparison between the chickens infected with Eimeria maxima (TNi) and their uninfected thermoneutral control (TNc) at 6 day-post-treatment.
T
he KEGG pathway analysis shown in Fig 3B and 3C shows the downregulation of the metabolic pathways involved in cellular respiration, carbohydrate, lipids, and amino acids catabolism, while the “cytokine-cytokine receptor interaction”, “toll-like receptor signaling, cell cycle”, and “necroptosis pathways” pathways were enriched. The GO terms are depicted in Fig 3D. The cellular compartments and molecular function showed downregulation of the terms related to absorption and transportation, especially for fatty acids, and enrichment of the terms involved in producing, binding, and transporting the inflammatory signal molecules. The biological process terms related to fatty acid metabolism and lipid homeostasis were downregulated, while the terms of immunity and inflammation, including immune regulatory and signaling processes, were upregulated.
The effect of E. maxima infection on the ileal transcriptome of heat-stressed chickens
A total of 1,908 DEGs were identified when the HSi group was compared with the HSc group. The top ~60 DEGs of the HSi group compared with the HSc group are exhibited in Fig 4A. The most downregulated genes in the HSi group compared with the HSc group included alcohol dehydrogenases 6 (class V) and 1C (class I), gamma polypeptide (ADH6 and ADH1C), glucose-6-phosphatase catalytic subunit 1 (G6PC), monoamine oxidase B (MAOB), D-aspartate oxidase (DDO), glycogen synthase 2 (GYS2) (Fig 4A). While the genes for IL-8, IL-22, IL-4I1, CCL4, CXCL13L2, TNFRSF6B, CTLA4, secreted phosphoprotein 1 (SPP1), Programmed cell death 1 ligand 1 (CD274), and claudin 1 (CLDN1) were of the top upregulated (Fig 4A).
Host transcriptome response to heat stress and Eimeria maxima infection in meat-type chickens - Image 8
Host transcriptome response to heat stress and Eimeria maxima infection in meat-type chickens - Image 9
Host transcriptome response to heat stress and Eimeria maxima infection in meat-type chickens - Image 10
Fig 4.
A. Heatmap of the top ~60 DEGs of HSc vs. HSi, as Log 2 relative hit counts of the comparison between the chickens infected with Eimeria maxima raised under HS (HSi) their uninfected HS control (HSc) at 6 day-post-treatment. B. The top KEGG pathways of HSc vs. HSi, based on the pathway term significance (α≤0.1) of the comparison between the chickens infected with Eimeria maxima raised under HS (HSi) and their uninfected HS control (HSc) at 6 day-post-treatment. The green squares and arrowheads depict the downregulated pathways and genes, respectively. The red circles and arrowheads depict the upregulated pathways and genes, respectively. The bigger size squares and circles denote the lower p-value. An increase in the color transparency (lighter color) of the squares and circles indicates an increase in the percentage of the included genes with expression that counter the regulation direction of the pathway. C. The top significant KEGG pathways of HSc vs. HSi, the comparison between the chickens infected with Eimeria maxima raised under HS (HSi) and their uninfected HS control (HSc) at 6 day-post-treatment. D. The top significant Gene Ontology terms of HSc vs. HSi, the comparison between the chickens infected with Eimeria maxima raised under HS (HSi) and their uninfected HS control (HSc) at 6 day-post-treatment.
The KEGG analysis depicted in Fig 4B and 4C shows downregulation of the pathways “fatty acids degradation”, “tryptophan metabolism”, “ABC transporters”, and “calcium signaling”. At the same time, the immune-related pathways “toll-like receptor” and “cytokine-cytokine receptor interaction”, “cell cycle” and its related pathways, “proteasome”, and “oxidative phosphorylation pathway” were among the significantly enriched pathways in the HSi group compared with the HSc group. The GO terms are presented in Fig 4D. The molecular function and biological process of the HSi showed downregulation of “lipid absorption and homeostasis”, “calcium transmembrane transporter activity”, and “anion transport, and ion channels”. The top enriched terms were the “defense”, ‘innate immune”, and “inflammatory responses” and related terms such as “positive regulation of cytokine production” and “activity and leukocyte activation”. Also, the “DNA polymerase”, “cell cycle”, “fructose 1,6-biphosphate”, and “fatty-acyl-CoA metabolic processes” were of the top enriched terms (Fig 4D). Consistently, the GO-cellular component showed downregulation of cell junction terms, such as “cell-cell adherence” and “anchoring junctions”, and “intercalated disc”. While, the “mitochondrial compartments”, “endoplasmic reticulum”, and “nucleus” terms were enriched.
The effect of HS on the ileal transcriptome of E. maxima-infected chickens
A total of 2,304 DEGs were identified when the HSi group was compared with the TNi groups. The genes interleukins 8, 1B, and 12 B (IL8, IL1B, and IL12B); the interleukin receptors 1R2 and 20RA (IL1R2 and IL20RA), secreted phosphoprotein 1 (SPP1), hepatocyte growth factor, and nitric oxide synthase 2 (NOS2) were among the most downregulated in the HSi group compared with the TNi group (Fig 5A). The top upregulated genes included apolipoproteins A1, A5, and C3 (APOA1, APOA5, and APOC3); fatty acid-binding proteins 1, 2, and 6 (FABP1, FABP2, and FABP6); cytochrome P450 family 1 A polypeptides 1 and 2 (CYP1A1 and CYP1A2); and glutamate-ammonia ligase (GLUL).
Host transcriptome response to heat stress and Eimeria maxima infection in meat-type chickens - Image 11
Host transcriptome response to heat stress and Eimeria maxima infection in meat-type chickens - Image 12
Host transcriptome response to heat stress and Eimeria maxima infection in meat-type chickens - Image 13
Fig 5.
A. Heatmap of the top ~60 DEGs of TNi vs. HSi, as Log 2 relative hit counts of the comparison between the chickens infected with Eimeria maxima raised either under HS (HSi) or thermoneutral condition (TNi) at 6-day-post-treatmen. B. The top KEGG pathways TNi vs. HSi, based on the pathway term significance (α≤0.05) of the comparison between the chickens infected with Eimeria maxima raised either under HS (HSi) or thermoneutral condition (TNi) at 6 day-post-treatment. The green squares and arrowheads depict the downregulated pathways and genes, respectively. The red circles and arrowheads depict the upregulated pathways and genes, respectively. The bigger size squares and circles denote the lower p-value. An increase in the color transparency (lighter color) of the squares and circles indicates an increase in the percentage of the included genes with expression that counter the regulation direction of the pathway. C. The top significant KEGG pathways of TNi vs. HSi, the comparison between the chickens infected with Eimeria maxima raised either under HS (HSi) or thermoneutral condition (TNi) at 6 day-post-treatment. D. The top significant Gene Ontology terms of TNi vs. HSi, the comparison between the chickens infected with Eimeria maxima raised either under HS (HSi) or thermoneutral condition (TNi) at 6 day-post-treatment.
The KEGG pathway analysis of the HSi group compared with the TNi group is depicted in Fig 5B and 5C. The top downregulated terms were “toll-like receptor signaling”, “cytokine-cytokine receptor interaction”, “focal adhesion”, and “calcium signaling” pathways. In contrast, the highest enriched pathways were the “citrate cycle”, “glycolysis/gluconeogenesis”, “oxidative phosphorylation”, “fatty acid degradation”, “glutathione metabolism”, “tryptophan metabolism”, and “steroid biosynthesis”. The GO terms are exhibited in Fig 5D. The top downregulated molecular function and biological process terms were the “calcium ion transmembrane regulation”, “inflammatory response”, “cytokine secretion and receptors”, “collagen and actin-binding”, and “leukocytes activation and differentiation”. Whereas, “antioxidant activity”, “glutathione metabolic process”, “cellular respiration, electron transport chain”, “lipid absorption, and fatty acid metabolism” were among the top enriched terms (Fig 5D). Concordantly, the GO-cellular component analysis showed downregulation of the terms “actin cytoskeleton”, “focal adhesion”, and “membrane raft”. While the terms “mitochondrial respiratory chain complex 1”, “endoplasmic reticulum”, “lysosome, and peroxisome” were the top enriched compartments (Fig 5D).
The RT-qPCR results of the target genes selected to confirm the Hi-Seq data showed a very high similarity in the expression levels to the fold change derived from the Hi-Seq analysis (S1–S4 Figs). The correlation coefficients (r) of the RT-qPCR values versus Hi-Seq yielded fold change were 0.758, 0.857, 0.891, and 0.999 for the pair comparisons TNc vs. HSc, TNc vs. TNi, HSc vs. HSi, and TNi vs. HSi; respectively. The values reflect a very strong correlation confirming the Hi-Seq results.
Discussion
Ileum tissue transcriptome change in response to HS
The ileal transcriptome was reshaped to maintain resources for cellular molecular functions to cope with HS-induced oxidative stress. Ileum cells upregulated the metabolism and degradation of multiple amino acids, such as lysine, histidine, valine, leucine, and isoleucine. This may plausibly be a compensatory mechanism to maintain cellular functions and energy production under low nutrient availability created by HS [15, 16, 31–33]. The results show the upregulation of multiple aldehyde dehydrogenase (ALDH) genes. KEGG pathway analysis shows that ALDH is connecting multiple amino acid degradations and metabolism pathways (Fig 2B). The ALDH7A1 gene encodes antiquitin enzyme (α-aminoadipic semialdehyde (α-AASA) dehydrogenase) that plays a fundamental role in lysine catabolism, the process by which lysine is broken down for energy production [34]. The enzyme α-AASA oxidizes multiple aldehyde containing molecules produced from metabolizing various amino acids by converting NAD+ to NADH [34]. The energy yield from each NADH is identified as 2.5 ATP produced in the electron transport chain [35]. The upregulated genes HAL and AMDHD1 encode enzymes involved in histidine metabolism producing L-glutamate [36]. The BCKDHB gene encodes the 2-oxoisovalerate dehydrogenase E1 component beta subunit involved in valine, leucine, and isoleucine degradation by producing acetoacetate and acetyl-CoA utilized in the tricarboxylic acid cycle for energy production [37]. The upregulated gene LIPG possesses both phospholipase and triglyceride lipase roles that mainly hydrolyzes the phospholipids part of lipoproteins and are generally involved in glycolipid metabolism [38]. Chickens under HS conditions endure visceral ischemia along with feed intake reduction resulting in severe nutrient scarcity at the cellular level [39]. The degradation and catabolism of lipoproteins and amino acids may be a potential source for producing energy and nutrient molecules in the ileum cells of chickens reared under chronic HS.
Stress hormones (mainly glucocorticoids and catecholamine) are released under the control of the neuroendocrine system in response to HS and were found to negatively impact both innate and adaptive immune functions [7]. Furthermore, nutritional limitations induced by chronic HS results putatively in a trade-off between maintaining the cellular basic functions and the immune response [40, 41]. Heat-stressed chickens expressed downregulation in the “cytokine-cytokine receptor interaction” KEGG pathway (Fig 2B) including many chemokines, cytokines, and cytokine receptor genes such as (CXCL13, CXCL13L2, IL22, IL8L1, IL1R2, and TNFRSF13B) that play a fundamental role in the innate immune response by attracting the immune cells to the affected site and initiating the inflammatory response (Fig 2B) [42]. The “cytokine-cytokine receptor interaction” pathway includes other downregulated genes that affect the cell-mediated and adaptive immune response such as FAS, IL18, and IFNAR1 [43, 44]. Also, the HSc group showed downregulation in gene CCR6, which is the key connection between the immature dendritic cells and the B cell lineage maturation for developing an adaptive immune response [45]. Consequently, the HS group show downregulation of KEGG pathways “influenza A” and “intestinal immune network for IgA production” indicating the negative regulation of both adaptive immune response compartments, cell-mediated and humeral immunity (Fig 2C & 2D).
Additionally, ER responds to cellular stress by multiple mechanisms including sensing the unfolded and misfolded proteins and then either correcting the protein folding or degrading it [46]. Chronic HS results in an increase the ROS production causing severe cellular oxidative stress that increases the rate of unfolded and misfolded protein production [47]. Eventually, ER reduces the protein synthesis to decrease the unfolded and misfolded protein production rate and also in response to reduced resources [48]. The HS group showed downregulation of the “protein processing in the ER” KEGG pathway including many downregulated genes that are responsible for recognizing the misfolded and unfolded proteins such as DNAJC3, HYOU1, and HSP90B1 [29, 49]. The downregulated pathway also includes many genes related to ER-associated protein degradation such as HSPH1, HSP90AA1, HSPA8, and HSPA4L (Fig 2A and 2B) [29, 49]. The downregulation of the ER lumen and chaperone machinery was observed in the GO analysis terms (cellular compartments, biological function, and molecular process) (Fig 2D). Exposing chickens to chronic HS resulted in amino acid degradation for energy, and a reduction in protein synthesis and overall ER functions.
Ileum tissue transcriptome change in response to E. maxima infection
The E. maxima-infected chickens expectedly expressed upregulation of the immune-related pathways and biological functions (Fig 3B–3D). Immune-related KEGG pathways such as “cytokine-cytokine receptor interaction”, “proteasome”, and “homologs recombination” were upregulated to establish the required innate and adaptive immune responses against E. maxima parasite (Fig 3B) [50]. The immune response is associated with increasing the transcription rate, especially for the genes encoded for inflammatory mediators, cytokines, and chemokines [51]. Therefore, there was upregulation of KEGG pathways “mRNA surveillance”, “spliceosome”, and “RNA transport” reflecting the cellular attempts for eukaryotic mRNA maturation as an initial step for protein synthesis [52]. The adaptive immune response against coccidiosis involves lymphocyte activation and homologous recombination that showed upregulation in the KEGG pathway and GO terms analysis [53].
Additionally, the “glutathione metabolism” KEGG pathway was downregulated, and some pathway-related genes were upregulated including GPX1 and 2 which encode major enzymes that catalyze the reduction of hydrogen peroxide into water. The GPX1 and 2 enzymes contribute to protecting the cell from oxidative damage induced by E. maxima infection by catalyzing the scavenger of the oxygen radical by glutathione [54–56]. The downregulation of multiple alcohol dehydrogenase genes such as ADH1C, ADH6, ADH4, ACADSB, and ACAA2, acyl-CoA synthetase genes (ACSL1 and 3), acyl-CoA dehydrogenase (ACADS and ACOX3) that encode enzymes contributing to long and short chain fatty acid metabolism suggests the downregulation of fatty acid metabolic pathways and disturbances of cellular lipid homeostasis [57, 58]. Alteration of lipid metabolism was also reported in broiler chickens infected with E. acervulina [59] and sub-clinical necrotic enteritis [60].
The effect of combining HS and E. maxima infection on the cellular immune functions
Comparing the HSi group with the HSc group showed a shift in the immune response related pathways triggered by E. maxima infection, since both groups were heat stressed. There was a significant upregulation in KEGG pathways “Toll-like receptor signaling”, “cytokine-cytokine receptor interaction”, and “Salmonella infection” suggesting that the infection might alter the previously explained immune suppression induced by HS (Fig 4B) [7]. However, the immune response seemed to be of lower expression capacity in the HSi chickens compared with the TNi chickens due to the downregulation of the “homologous combination” pathway required for lymphocyte differentiation to develop an adaptive immune response. Meanwhile, the chickens showed upregulation of the GO molecular function term for regulation of wound healing when compared with HSc (Fig 4D) but not in the TNi when compared with TNc (Fig 3D). Taken together, HSi showed an altered immune response and tissue repair functions than TNi when each was compared to its uninfected control at the same temperature level.
Comparing the HSi group against the TNi group reveals how HS altered the cellular response to the E. maxima infection and showed downregulation of all the immune response and signaling-related pathways and GO terms (Fig 5B–5D). The downregulation of the “Toll-like receptor signaling” pathway may have resulted in the reduction in the ability of the immune system to identify the E. maxima antigen and trigger an inflammatory and immune response. Consequently, all the cellular and molecular functions such as cytokine secretion and leukocytic activation contributing to the downstream cascade of immune response were downregulated. Similarly, HS diminished the immune response to Salmonella Typhimurium infection in broiler chickens [61]. The role of macrophages and cytotoxic (CD 8+) and peripheral lymphocytes in Eimeria spp sporozoite transportation through intestinal tissue laminae has been documented [62–64]. Schneider et al. [65] reported that HS curtailed the E. maxima life cycle and reduced gametogony in broiler chickens. Thus, it can be hypothesized that suppressing the lymphocytes and monocyte migration induced by HS might contribute to reducing the sporozoite transportation which disrupted the parasite’s life cycle leading to reducing its gametogony capacity.
The effect of combining HS and E. maxima infection on the cellular signaling and metabolic functions
The upregulation of the “glutamate metabolism”, “pyridine metabolism”, “steroid biosynthesis”, and cell respiration-related KEGG pathways in the HSi group compared with either HSc or TNi groups may suggest improving the enterocytes functions and proliferation when the two stressors (HS and E. maxima infection) are combined rather than individually applied [66, 67]. The metabolism of both the sulfur-containing amino acids (cysteine and methionine) and glutathione was upregulated in the HSi compared to TNi (Fig 5B). The genes GSS and GPX1 encoding enzymes glutathione synthetase and glutathione peroxidase1, respectively, were upregulated. Glutathione synthetase catalyzes the glutathione production form γ-glutamyl cysteine and glycine in an ATP-dependent reaction [68]. Glutathione peroxidase 1 catalyzes the glutathione antioxidant activity by scavenger ROS from hydrogen peroxide [69]. This indicates that the cellular antioxidation machinery was possibly enhanced in the HSi group compared with the TNi group. Interestingly, the molecular function GO terms showed responses to oxidative stress and low oxygen level (Fig 5D) regularly associated with HS were downregulated in the HSi group compared with the HSc group suggesting that the local hyperemia induced by E. maxima infection reversed the HS-induced intestinal ischemia.
Calcium signaling is required for the development of the protozoan Plasmodium spp which is closely classified with Eimeria spp in the phylum Apicomplexa [70]. Depleting the cell culture calcium suppressed the Plasmodium asexual stage which is restored by excess calcium supplementation [71, 72]. The results showed downregulation of molecular function GO terms “positive regulation of calcium ion transport” and “regulation of calcium-mediated signaling”, and KEGG term “calcium signaling pathway” in HSi compared to both HSc and TNi (Figs 4B, 5B and 5D). That includes the downregulation of genes RYR2, CACNA1C, and ORAI3 encoded for ryanodine receptor 2, voltage-dependent L-type calcium channel subunit alpha-1C, and calcium release-activated calcium modulator 3, respectively. RYR is responsible for a calcium-induced-calcium release from the endoplasmic reticulum stimulated by increasing cytoplasmic calcium concentration [73]. CACNA1C is responsible for gating calcium ions current into the cell, which may activate RYR receptors [37, 74]. ORAI3 mediates calcium ions influx into the cell in response to the depletion of stored cellular calcium [75]. The catabolism of the essential amino acid tryptophan via the kynurenine pathway is necessary for the optimum development of E. falciformis in mice [65]. Tryptophan metabolism was also revealed as a metabolomic signature in E. acervulina infection [59]. The current study showed that the metabolism of the essential amino acid tryptophan was upregulated in the HSi group compared with the TNi. However, the upregulated genes mainly contribute to tryptophan metabolism through pathways other than kynurenine. For instance, the genes DDC, CYP1A1, and CYP1A2 participate in the serotonin pathway, and the genes MAOA, ALDH3A2, AOX2, and IL4I1 contribute to the indole pathway [76]. Taken together, reducing the cellular pool of nutrients and ions (tryptophan and calcium) that are scavenged by E. maxima for development and replication may explain the suppression of the parasite gametogony under HS [14, 77, 78].
The effect of combining HS and E. maxima infection on cellular proliferation and protein synthesis
When HSi is compared with HSc, cell cycle arrest was more prominent due to the upregulation of the “P53 signaling” KEGG pathway which is one of the major cell arrest inducers (Fig 4B). P53 protein is phosphorylated under stress conditions resulting in activating the transcription of the cyclins suppressors such as CDKN1A and GADD45A that arrest the cell cycle at the S, G2 and M phases [79, 80]. That negative feedback regulation seems to be accompanied with upregulation in the genes involved in cycle regulation (Fig 4B) to enhance the immune response and tissue repair in the HSi chickens compared with HSc group [16]. The upregulation of genes upstream of the cellular apoptosis such as FAS and SHISA5 enhanced the programmed death of the infected cell resulting in reducing inflammation, infection elimination, and prompt tissue repair [81]. The cell cycle pathway was downregulated in the HSi group compared with the TNi group. The upregulation of GO molecular function “DNA repair” accompanied with the upregulation of “histone kinase” which grants a negative charge to histone resulting in a higher level of chromatin opening conformation [82] suggesting that the DNA repair is enhanced in the HSi over the other groups with a single stressor.
However, KEGG pathways “ribosome biogenesis in eukaryotes” and “RNA degradation” were downregulated, the GO analysis showed upregulation of the cellular components related to the ER and the molecular function term “positive regulation of protein exit from the ER” in HSi when compared with TNi. Comparing HSi with HSc showed upregulation of KEGG pathways “ribosome” and “protein export”, and GO cellular components terms “endoplasmic reticulum chaperon complex, “endoplasmic reticulum lumen”, “endoplasmic reticulum inner membrane”, and GO molecular functions “positive regulation of protein exit from ER” suggesting that the TNi group might express a higher level of cellular and molecular functions related to protein synthesis when compared either to the TNi or HSc group.
Conclusion
The ileal transcriptome profile was modified in response to either HS or E. maxima infection. Heat stress significantly affects the ileal transcriptome by shifting the cellular metabolic pathways from oxidative phosphorylation toward pathways exploiting the amino acids while reducing other pathways of crucial functions such as immune response and protein synthesis. Under E. maxima infection, cellular and molecular pathways involved in immune functions were enhanced to counter the parasitic invasion. However, the analysis was conducted at 6 dpi, the upregulation of pathways related to the primary immune response, such as PRR, were still prominent. Surprisingly, chickens exposed to a combination of HS and E. maxima infection showed a unique response. Both HS and E. maxima interact to modulate the cellular and molecular responses to each stressor. The HSi group showed lower expression levels of the oxidative stress-related pathways compared to the HSc group. Furthermore, the HSi group showed modifications in pathways involved in nutrient metabolism and calcium signaling compared with other treatment groups. We hypothesized that limiting the nutrients available to be scavenged by the parasite and/or hindering the sporozoite transportation by macrophages and lymphocytes may be two possible mechanisms by which HS interferes with the E. maxima life cycle.
     
This article was originally published in . PLoS ONE 19(2): e0296350. https://doi.org/10.1371/journal. pone.0296350. This is an Open Access article distributed under the terms of the Creative Commons Attribution License.

1. Ashraf S, Zaneb H, Yousaf MS, Ijaz A, Sohail MU, Muti S, et al. Effect of dietary supplementation of prebiotics and probiotics on intestinal microarchitecture in broilers reared under cyclic heat stress. J Anim Physiol Anim Nutr (Berl). 2013;97 Suppl 1:68–73. Epub 2013/05/10. https://doi.org/10.1111/jpn.12041 PMID: 23639019.

2. Nardone A, Ronchi B, Lacetera N, Ranieri MS, Bernabucci U. Effects of climate changes on animal production and sustainability of livestock systems. Livestock Science. 2010; 130(1–3):57–69.

3. Ghazi S, Habibian M, Moeini M, Abdolmohammadi A. Effects of different levels of organic and inorganic chromium on growth performance and immunocompetence of broilers under heat stress. Biological trace element research. 2012; 146:309–17. https://doi.org/10.1007/s12011-011-9260-1 PMID: 22127829

4. Imik H, Ozlu H, Gumus R, Atasever MA, Urcar S, Atasever M. Effects of ascorbic acid and α-lipoic acid on performance and meat quality of broilers subjected to heat stress. British Poultry Science. 2012; 53 (6):800–8.

5. Marder J, Arad Z. Panting and acid-base regulation in heat stressed birds. Comp Biochem Physiol A Comp Physiol. 1989; 94(3):395–400. Epub 1989/01/01. https://doi.org/10.1016/0300-9629(89)90112-6 PMID: 2574090.

6. Butts CL, Sternberg EM. Neuroendocrine factors alter host defense by modulating immune function. Cell Immunol. 2008; 252(1–2):7–15. Epub 2008/03/11. https://doi.org/10.1016/j.cellimm.2007.09.009 PMID: 18329009; PubMed Central PMCID: PMC2590632.

7. Marketon JIW, Glaser RJCi. Stress hormones and immune function. 2008; 252(1–2):16–26. https://doi. org/10.1016/j.cellimm.2007.09.006 PMID: 18279846

8. Ruff J, Barros TL, Tellez G Jr, Blankenship J, Lester H, Graham BD, et al. Research Note: Evaluation of a heat stress model to induce gastrointestinal leakage in broiler chickens. Poultry Science. 2020; 99 (3):1687–92. https://doi.org/10.1016/j.psj.2019.10.075 PMID: 32115037

9. Renaudeau D, Collin A, Yahav S, De Basilio V, Gourdine J-L, Collier R. Adaptation to hot climate and strategies to alleviate heat stress in livestock production. Animal. 2012; 6(5):707–28. https://doi.org/10. 1017/S1751731111002448 PMID: 22558920

10. Habashy WS, Milfort MC, Rekaya R, Aggrey SE. Cellular antioxidant enzyme activity and biomarkers for oxidative stress are affected by heat stress. International journal of biometeorology. 2019; 63:1569– 84. https://doi.org/10.1007/s00484-019-01769-z PMID: 31352522

11. Dubey JP, Jenkins MC. Re-evaluation of the life cycle of Eimeria maxima Tyzzer, 1929 in chickens (Gallus domesticus). Parasitology. 2018; 145(8):1051–8. Epub 2017/12/15. https://doi.org/10.1017/ S0031182017002153 PMID: 29239290.

12. Schnitzler BE, Shirley MW. Immunological aspects of infections with Eimeria maxima: a short review. Avian Pathology. 1999; 28(6):537–43.

13. Schneiders G, Foutz J, Milfort M, Ghareeb A, Sorhue U, Richter J, et al. Ontogeny of intestinal permeability in chickens infected with Eimeria maxima: implications for intestinal health. J Adv Parasitol. 2019; 6(3):41–50.

14. Schneiders GH, Foutz JC, Milfort MC, Ghareeb AF, Fuller AL, Rekaya R, et al. Heat stress reduces sexual development and affects pathogenesis of Eimeria maxima in meat-type chickens. Scientific Reports. 2020; 10(1):10736.

15. Ghareeb AF, Schneiders GH, Foutz JC, Milfort MC, Fuller AL, Yuan J, et al. Heat Stress Alters the Effect of Eimeria maxima Infection on Ileal Amino Acids Digestibility and Transporters Expression in Meat-Type Chickens. Animals. 2022; 12(12):1554.

16. Ghareeb AFA, Schneiders GH, Richter JN, Foutz JC, Milfort MC, Fuller AL, et al. Heat stress modulates the disruptive effects of Eimeria maxima infection on the ileum nutrient digestibility, molecular transporters, and tissue morphology in meat-type chickens. PLoS One. 2022; 17(6):e0269131. Epub 2022/0 04. https://doi.org/10.1371/journal.pone.0269131 PMID: 35657942; PubMed Central PMCID: PMC9165794.

17. Adams JJNE. Transcriptome: connecting the genome to gene function. Nature Education. 2008; 1 (1):195.

18. Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014; 30(15):2114–20. https://doi.org/10.1093/bioinformatics/btu170 P

19. Dobin A, Davis CA, Schlesinger F, Drenkow J, Zaleski C, Jha S, et al. STAR: ultrafast universal RNA seq aligner. Bioinformatics. 2013; 29(1):15–21. https://doi.org/10.1093/bioinformatics/bts635 PMID: 23104886

20. Liao Y, Smyth GK, Shi W. featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics. 2014; 30(7):923–30. https://doi.org/10.1093/bioinformatics/ btt656 PMID: 24227677

21. Dodt M, Roehr JT, Ahmed R, Dieterich C. FLEXBAR—flexible barcode and adapter processing for next-generation sequencing platforms. Biology. 2012; 1(3):895–905. https://doi.org/10.3390/ biology1030895 PMID: 24832523

22. Barnett DW, Garrison EK, Quinlan AR, Stro¨mberg MP, Marth GT. BamTools: a C++ API and toolkit for analyzing and managing BAM files. Bioinformatics. 2011; 27(12):1691–2. https://doi.org/10.1093/ bioinformatics/btr174 PMID: 21493652

23. Trapnell C, Roberts A, Goff L, Pertea G, Kim D, Kelley DR, et al. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nature protocols. 2012; 7(3):562–78. https://doi.org/10.1038/nprot.2012.016 PMID: 22383036

24. Anders S, Huber W. Differential expression analysis for sequence count data. Nature Precedings. 2010:1-. https://doi.org/10.1186/gb-2010-11-10-r106 PMID: 20979621

25. Benjamini Y, Hochberg Y. Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. Journal of the Royal Statistical Society: Series B (Methodological). 1995; 57(1):289– 300. https://doi.org/10.1111/j.2517-6161.1995.tb02031.x

26. Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome research. 2003; 13 (11):2498–504. https://doi.org/10.1101/gr.1239303 PMID: 14597658

27. Bindea G, Mlecnik B, Hackl H, Charoentong P, Tosolini M, Kirilovsky A, et al. ClueGO: a Cytoscape plug-in to decipher functionally grouped gene ontology and pathway annotation networks. Bioinformatics. 2009; 25(8):1091–3. https://doi.org/10.1093/bioinformatics/btp101 PMID: 19237447

28. Ogata H, Goto S, Sato K, Fujibuchi W, Bono H, Kanehisa M. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic acids research. 1999; 27(1):29–34. https://doi.org/10.1093/nar/27.1.29 PMID: 9847135

29. Kanehisa M, Furumichi M, Sato Y, Ishiguro-Watanabe M, Tanabe M. KEGG: integrating viruses and cellular organisms. Nucleic acids research. 2021; 49(D1):D545–D51. https://doi.org/10.1093/nar/ gkaa970 PMID: 33125081

30. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2− ΔΔCT method. methods. 2001; 25(4):402–8.

31. Syafwan S, Kwakkel R, Verstegen M. Heat stress and feeding strategies in meat-type chickens. World’s Poultry Science Journal. 2011; 67(4):653–74.

32. Habashy W, Milfort M, Adomako K, Attia Y, Rekaya R, Aggrey S. Effect of heat stress on amino acid digestibility and transporters in meat-type chickens. Poultry science. 2017; 96(7):2312–9. https://doi. org/10.3382/ps/pex027 PMID: 28339933

33. Habashy WS, Milfort MC, Fuller AL, Attia YA, Rekaya R, Aggrey SE. Effect of heat stress on protein utilization and nutrient transporters in meat-type chickens. International Journal of Biometeorology. 2017; 61:2111–8. https://doi.org/10.1007/s00484-017-1414-1 PMID: 28799035

34. Brocker C, Lassen N, Estey T, Pappa A, Cantore M, Orlova VV, et al. Aldehyde dehydrogenase 7A1 (ALDH7A1) is a novel enzyme involved in cellular defense against hyperosmotic stress. Journal of Biological Chemistry. 2010; 285(24):18452–63. https://doi.org/10.1074/jbc.M109.077925 PMID: 20207735

35. Feher JJ. Quantitative human physiology: an introduction: Academic press; 2017.

36. Petrova B, Kanarek N. Potential Benefits and Pitfalls of Histidine Supplementation for Cancer Therapy Enhancement. J Nutr. 2020; 150(Suppl 1):2580S–7S. Epub 2020/10/02. https://doi.org/10.1093/jn/ nxaa132 PMID: 33000153.

37. Rougier JS, Abriel H. Cardiac voltage-gated calcium channel macromolecular complexes. Biochim Biophys Acta. 2016; 1863(7 Pt B):1806–12. Epub 2015/12/29. https://doi.org/10.1016/j.bbamcr.2015.12. 014 PMID: 26707467.

38. Jin W, Broedl UC, Monajemi H, Glick JM, Rader DJ. Lipase H, a new member of the triglyceride lipase family synthesized by the intestine. Genomics. 2002; 80(3):268–73. https://doi.org/10.1006/geno.2002. 6837 PMID: 12213196

39. Cronje PJRAiANiA. Heat stress in livestock—The role of the gut in its aetiology and a potential role for betaine in its alleviation. 2005; 15:107–22.

40. Norris K, Evans MRJBE. Ecological immunology: life history trade-offs and immune defense in birds. 2000; 11(1):19–26.

41. Lucas LD, French SSJPo. Stress-induced tradeoffs in a free-living lizard across a variable landscape: consequences for individuals and populations. 2012; 7(11):e49895.

42. Borish LC, Steinke JW. 2. Cytokines and chemokines. Journal of Allergy and Clinical Immunology. 2003; 111(2):S460–S75. https://doi.org/10.1067/mai.2003.108 PMID: 12592293

43. Van Parijs L, Abbas AK. Role of Fas-mediated cell death in the regulation of immune responses. Current opinion in immunology. 1996; 8(3):355–61. https://doi.org/10.1016/s0952-7915(96)80125-7 PMID: 8793992

44. Lebratti T, Lim YS, Cofie A, Andhey P, Jiang X, Scott J, et al. A sustained type I IFN-neutrophil-IL-18 axis drives pathology during mucosal viral infection. Elife. 2021; 10:e65762. https://doi.org/10.7554/ eLife.65762 PMID: 34047696

45. Ito T, Carson IV WF, Cavassani KA, Connett JM, Kunkel SL. CCR6 as a mediator of immunity in the lung and gut. Experimental cell research. 2011; 317(5):613–9. https://doi.org/10.1016/j.yexcr.2010.12. 018 PMID: 21376174

46. Schwarz DS, Blower MD. The endoplasmic reticulum: structure, function and response to cellular signaling. Cellular and molecular life sciences. 2016; 73:79–94. https://doi.org/10.1007/s00018-015-2052- 6 PMID: 26433683

47. Azad Kalam M, Motoi K, Hoque Azharul M, Masaaki T. Effect of Chronic Heat Stress on Performance and Oxidative Damage in Different Strains of Chickens. 日本家禽学会誌. 2010; 47(4):333–7.

48. Lee M, Park H, Heo JM, Choi HJ, Seo S. Multi-tissue transcriptomic analysis reveals that L-methionine supplementation maintains the physiological homeostasis of broiler chickens than D-methionine under acute heat stress. PLoS One. 2021; 16(1):e0246063. Epub 2021/01/28. https://doi.org/10.1371/journal. pone.0246063 PMID: 33503037; PubMed Central PMCID: PMC7840013.

49. Kanehisa M. Toward understanding the origin and evolution of cellular organisms. Protein Science. 2019; 28(11):1947–51. https://doi.org/10.1002/pro.3715 PMID: 31441146

50. Memon FU, Yang Y, Leghari IH, Lv F, Soliman AM, Zhang W, et al. Transcriptome analysis revealed ameliorative effects of Bacillus based probiotic on immunity, gut barrier system, and metabolism of chicken under an experimentally induced Eimeria tenella infection. Genes. 2021; 12(4):536. https://doi. org/10.3390/genes12040536 PMID: 33917156

51. Powell N, Canavan J, MacDonald T, Lord G. Transcriptional regulation of the mucosal immune system mediated by T-bet. Mucosal immunology. 2010; 3(6):567–77. https://doi.org/10.1038/mi.2010.53 PMID: 20844482

52. Bremner A, Kim S, Morris KM, Nolan MJ, Borowska D, Wu Z, et al. Kinetics of the cellular and transcriptomic response to Eimeria maxima in relatively resistant and susceptible chicken lines. Frontiers in Immunology. 2021; 12:653085.

53. Lillehoj H, Min W, Dalloul R. Recent progress on the cytokine regulation of intestinal immune responses to Eimeria. Poultry Science. 2004; 83(4):611–23. https://doi.org/10.1093/ps/83.4.611 PMID: 15109059

54. Galli GM, Baldissera MD, Griss LG, Souza CF, Fortuoso BF, Boiago MM, et al. Intestinal injury caused by Eimeria spp. impairs the phosphotransfer network and gain weight in experimentally infected chicken chicks. Parasitology Research. 2019; 118:1573–9. https://doi.org/10.1007/s00436-019-06221-0 PMID: 30815727

55. Khatlab AdS Del Vesco AP, de Oliveira Neto AR, Fernandes RPM, Gasparino E. Dietary supplementation with free methionine or methionine dipeptide mitigates intestinal oxidative stress induced by Eimeria spp. challenge in broiler chickens. Journal of Animal Science and Biotechnology. 2019; 10(1):1–17.

56. Alehagen U, Opstad TB, Alexander J, Larsson A, Aaseth J. Impact of selenium on biomarkers and clinical aspects related to ageing. A review. Biomolecules. 2021; 11(10):1478. https://doi.org/10.3390/ biom11101478 PMID: 34680111

57. Hunt MC, Alexson SE. The role Acyl-CoA thioesterases play in mediating intracellular lipid metabolism. Progress in lipid research. 2002; 41(2):99–130. https://doi.org/10.1016/s0163-7827(01)00017-0 PMID: 11755680

58. Ho¨o¨g J-O, Stro¨mberg P, Hedberg JJ, Griffiths WJ. The mammalian alcohol dehydrogenases interact in several metabolic pathways. Chemico-biological interactions. 2003; 143:175–81. https://doi.org/10. 1016/s0009-2797(02)00225-9 PMID: 12604202

59. Aggrey SE, Milfort MC, Fuller AL, Yuan J, Rekaya R. Effect of host genotype and Eimeria acervulina infection on the metabolome of meat-type chickens. PLoS One. 2019; 14(10):e0223417. Epub 2019/ 10/17. https://doi.org/10.1371/journal.pone.0223417 PMID: 31618222; PubMed Central PMCID: PMC6795442.

60. Qing X, Zeng D, Wang H, Ni X, Lai J, Liu L, et al. Analysis of hepatic transcriptome demonstrates altered lipid metabolism following Lactobacillus johnsonii BS15 prevention in chickens with subclinical necrotic enteritis. Lipids in Health and Disease. 2018; 17:1–10.

61. Tang L-P, Li W-H, Liu Y-L, Lun J-C, He Y-M. Heat stress inhibits expression of the cytokines, and NF κB-NLRP3 signaling pathway in broiler chickens infected with salmonella typhimurium. Journal of Thermal Biology. 2021; 98:102945.

62. DORAN DJ. The migration of Eimeria acervulina sporozoites to the duodenal glands of Lieberku¨hn. The Journal of Protozoology. 1966; 13(1):27–33.

63. Trout J, Lillehoi H. Evidence of a role for intestinal CD8+ lymphocytes and macrophages in transport of Eimeria acervulina sporozoites. The Journal of parasitology. 1993; 79(5):790–2. PMID: 8105047

64. Beattie SE, Barta JR, Fernando M. Involvement of CD 8+ and CD 3+ lymphocytes in the transport of Eimeria necatrix sporozoites within the intestinal mucosa of chickens. Parasitology research. 2001; 87:405–8. https://doi.org/10.1007/s004360000371 PMID: 11403384

65. Schmid M, Lehmann MJ, Lucius R, Gupta N. Apicomplexan parasite, Eimeria falciformis, co-opts host tryptophan catabolism for life cycle progression in mouse. Journal of Biological Chemistry. 2012; 287 (24):20197–207. https://doi.org/10.1074/jbc.M112.351999 PMID: 22535959

66. McCauley R, Kong SE, Hall J. Glutamine and nucleotide metabolism within enterocytes. Journal of parenteral and enteral nutrition. 1998; 22(2):105–11. https://doi.org/10.1177/0148607198022002105 PMID: 9527969

67. Sakai R, Ooba Y, Watanabe A, Nakamura H, Kawamata Y, Shimada T, et al. Glutamate metabolism in a human intestinal epithelial cell layer model. Amino Acids. 2020; 52:1505–19. https://doi.org/10.1007/ s00726-020-02908-2 PMID: 33180203

68. Njålsson R. Glutathione synthetase deficiency. Cellular and Molecular Life Sciences CMLS. 2005; 62:1938–45. https://doi.org/10.1007/s00018-005-5163-7 PMID: 15990954

69. Lei XG, Cheng W-H, McClung JP. Metabolic regulation and function of glutathione peroxidase-1. Annu Rev Nutr. 2007; 27:41–61. https://doi.org/10.1146/annurev.nutr.27.061406.093716 PMID: 17465855

70. Zipprer EM, Neggers M, Kushwaha A, Rayavara K, Desai SA. A kinetic fluorescence assay reveals unusual features of Ca++ uptake in Plasmodium falciparum-infected erythrocytes. Malaria journal. 2014; 13(1):1–11.

71. Wasserman M, Alarco´n C, Mendoza PM. Effects of Ca++ depletion on the asexual cell cycle of Plasmodium falciparum. The American journal of tropical medicine and hygiene. 1982; 31(4):711–7. https://doi. org/10.4269/ajtmh.1982.31.711 PMID: 6808848

72. Gupta Y, Goicoechea S, Pearce CM, Mathur R, Romero JG, Kwofie SK, et al. The emerging paradigm of calcium homeostasis as a new therapeutic target for protozoan parasites. Medicinal Research Reviews. 2022; 42(1):56–82. https://doi.org/10.1002/med.21804 PMID: 33851452

73.  Van Petegem FJJoBC. Ryanodine receptors: structure and function. 2012; 287(38):31624–32

74. Uurasmaa T-M, Streng T, Alkio M, Heinonen I, Anttila K. Short-term exercise affects cardiac function ex vivo partially via changes in calcium channel levels, without influencing hypoxia sensitivity. Journal of physiology and biochemistry. 2021; 77(4):639–51. https://doi.org/10.1007/s13105-021-00830-z PMID: 34449060

75. Hoth M, Niemeyer BAJCtim. The neglected CRAC proteins: Orai2, Orai3, and STIM2. 2013; 71:237– 71. https://doi.org/10.1016/B978-0-12-407870-3.00010-X PMID: 23890118

76. Roth W, Zadeh K, Vekariya R, Ge Y, Mohamadzadeh M. Tryptophan metabolism and gut-brain homeostasis. International journal of molecular sciences. 2021; 22(6):2973. https://doi.org/10.3390/ ijms22062973 PMID: 33804088

77. Saliba KJ, Kirk K. Nutrient acquisition by intracellular apicomplexan parasites: staying in for dinner. International journal for parasitology. 2001; 31(12):1321–30. https://doi.org/10.1016/s0020-7519(01) 00258-2 PMID: 11566300

78. Landfear SM. Nutrient transport and pathogenesis in selected parasitic protozoa. Eukaryotic cell. 2011; 10(4):483–93. https://doi.org/10.1128/EC.00287-10 PMID: 21216940

79. Gentile M, Latonen L, Laiho M. Cell cycle arrest and apoptosis provoked by UV radiation-induced DNA damage are transcriptionally highly divergent responses. Nucleic acids research. 2003; 31(16):4779– 90. https://doi.org/10.1093/nar/gkg675 PMID: 12907719

80. Kleinsimon S, Longmuss E, Rolff J, Ja¨ger S, Eggert A, Delebinski C, et al. GADD45A and CDKN1A are involved in apoptosis and cell cycle modulatory effects of viscumTT with further inactivation of the STAT3 pathway. Scientific reports. 2018; 8(1):5750. https://doi.org/10.1038/s41598-018-24075-x PMID: 29636527

81. Engeland KJCD, Differentiation. Cell cycle regulation: p53-p21-RB signaling. 2022; 29(5):946–60.

82. Singh D, Nishi K, Khambata K, Balasinor N. Introduction to epigenetics: basic concepts and advancements in the field. Elsevier; 2020. p. xxv–xliv

Related topics:
Related Questions

Stress hormones (mainly glucocorticoids and catecholamine) are released under the control of the neuroendocrine system in response to HS and were found to negatively impact both innate and adaptive immune functions [7]. Furthermore, nutritional limitations induced by chronic HS results putatively in a trade-off between maintaining the cellular basic functions and the immune response [40, 41]. Heat-stressed chickens expressed downregulation in the “cytokine-cytokine receptor interaction” KEGG pathway (Fig 2B) including many chemokines, cytokines, and cytokine receptor genes such as (CXCL13, CXCL13L2, IL22, IL8L1, IL1R2, and TNFRSF13B) that play a fundamental role in the innate immune response by attracting the immune cells to the affected site and initiating the inflammatory response (Fig 2B) [42]. The “cytokine-cytokine receptor interaction” pathway includes other downregulated genes that affect the cell-mediated and adaptive immune response such as FAS, IL18, and IFNAR1 [43, 44]. Also

Additionally, ER responds to cellular stress by multiple mechanisms including sensing the unfolded and misfolded proteins and then either correcting the protein folding or degrading it [46]. Chronic HS results in an increase the ROS production causing severe cellular oxidative stress that increases the rate of unfolded and misfolded protein production [47]. Eventually, ER reduces the protein synthesis to decrease the unfolded and misfolded protein production rate and also in response to reduced resources [48]. The HS group showed downregulation of the “protein processing in the ER” KEGG pathway including many downregulated genes that are responsible for recognizing the misfolded and unfolded proteins such as DNAJC3, HYOU1, and HSP90B1 [29, 49]. The downregulated pathway also includes many genes related to ER-associated protein degradation such as HSPH1, HSP90AA1, HSPA8, and HSPA4L (Fig 2A and 2B) [29, 49]. The downregulation of the ER lumen and chaperone machinery was observed in the G
Authors:
James Foutz
Gustavo Schneiders
Jennifer Richter
University of Georgia
University of Georgia
Marie Milfort
University of Georgia
University of Georgia
Lorraine Fuller
University of Georgia
University of Georgia
Romdhane Rekaya
University of Georgia
University of Georgia
Samuel Aggrey
University of Georgia
University of Georgia
Show more
Recommend
Comment
Share
Profile picture
Would you like to discuss another topic? Create a new post to engage with experts in the community.
Featured users in Poultry Industry
Padma Pillai
Padma Pillai
Cargill
United States
Kendra Waldbusser
Kendra Waldbusser
Pilgrim´s
United States
Karen Christensen
Karen Christensen
Tyson
Tyson
PhD, senior director of animal welfare at Tyson Foods
United States
Join Engormix and be part of the largest agribusiness social network in the world.