Cecal microbiota composition differs under normal and high ambient temperatures in genetically distinct chicken lines

Modern broilers, selected for high growth rate, are more susceptible to heat stress (HS) as compared to their ancestral jungle fowl (JF). HS affects epithelia barrier integrity, which is associated with gut microbiota. The aim of this study was to determine the effect of HS on the cecal luminal (CeL) and cecal mucosal (CeM) microbiota in JF and three broiler populations: Athens Canadian Random Bred (ACRB), 1995 Random Bred (L1995), and Modern Random Bred (L2015). Broiler chicks were subjected to thermoneutral TN (24 °C) or chronic cyclic HS (8 h/day, 36 °C) condition from day 29 until day 56. HS affected richness in CeL microbiota in a line-dependent manner, decreasing richness in slow-growing JF and ACRB lines, while increasing richness in faster-growing L1995 and L2015. Microbiota were distinct between HS and TN conditions in CeL microbiota of all four lines and in CeM microbiota of L2015. Certain bacterial genera were also affected in a line-dependent manner, with HS tending to increase relative abundance in CeL microbiota of slow-growing lines, while decreases were common in fast-growing lines. Predictive functional analysis suggested a greater impact of HS on metabolic pathways in L2015 compared to other lines.


Animals and tissue sampling
The study was conducted in accordance with the recommendations in the guide for the care and use of laboratory animals of the National Institutes of Health and the protocols were approved by the University of Arkansas Animal Care and Use Committee under protocols 18,083 and 16,084.This study was performed and reported in accordance with ARRIVE guidelines (https:// arriv eguid elines.org/).Four lines of chicken were utilized for this study: Giant Jungle Fowl (JF), Athens Canadian Random Bred (ACRB), 1995 Arkansas Random Bred (L1995), and Modern Random Bred (L2015).Detailed characteristics of the lines were previously reported in Tabler et al. 3 .Protocols for incubation of embryonated eggs, hatching, housing conditions, and experimental design were as reported in Emami et al. 19 .In brief, to investigate the effect of heat stress (HS) on gut microbiota, chicks were separated by line and placed into twelve environmental chambers (24 pens total, 6 pens/line, 25 birds/pen, 0.09 m 2 /bird).Chamber temperatures were 32 °C at day 1 and gradually decreased to 20 °C at day 21.Water and 3-phase standard diets were provided ad libitum.At day 29, half of the pens for each line (3 pens/line) were raised under thermoneutral (TN, 24 °C) conditions, while 3 pens per line were subjected to chronic cyclic HS (8 h/day, 36 °C from 9 AM to 5 PM).On day 56, two birds per pen were selected based on the average pen weight and euthanized via cervical dislocation.The cecal contents (luminal, CeL) and cecal epithelial scrapings (mucosal, CeM) were collected for bacterial DNA sequencing.

DNA isolation, library preparation, and sequencing
DNA isolation, library preparation, and sequencing were as reported in Emami et al. 19 .A PowerSoil kit (Qiagen, Valencia, CA), PCR primers targeting the V3-V4 region of the 16S rRNA gene, and the Illumina MiSeq platform (Illumina, Inc) were utilized for these steps, respectively.The 16S rRNA gene sequences determined in this study were deposited in the NCBI Sequence Read Archive database (SRA accession no.PRJNA930873).

Bioinformatics and statistical analysis
Microbiota analyses on CeL and CeM microbiota were performed using the bioinformatics platform Quantitative Insights into Microbial Ecology 2 (QIIME 2) version 2022.8 22 .Quality control was performed through denoising with DADA2 23 via the q2-dada2 plugin, setting truncation parameters using a quality cutoff of 30.The SILVA rRNA database version 138 was utilized for taxonomic analyses by downloading reference sequences and taxonomy files pre-formatted for QIIME 2 (obtained at https:// docs.qiime2.org/ 2022.8/data-resou rces/) using RESCRIPt, a process that removes duplicate sequences assigned to different taxonomies to reduce inconsistencies and improve processing 24 .Reads were extracted from the reference sequences using the forward and reverse primers for the V3-V4 region, and the extracted reads were used to create a feature classifier via q2-featureclassifier 25 .Taxonomy was assigned to amplicon sequence variants (ASVs) from DADA2 via the q2-featureclassifier classify-sklearn naïve Bayes taxonomy classifier.Mitochondria, chloroplasts, and unassigned bacteria were filtered and excluded from the feature table.To construct a phylogeny, ASVs were aligned with MAFFT 26 via q2-alignment and passed to fasttree2 27 via q2-phylogeny.Rarefaction, or subsampling without replacement, was performed with sampling depths of 17,055 for CeL and 14,247 for CeM for alpha and beta diversity analyses via q2-diversity.Sampling depths were determined based on diversity captured at different depths (visualized by alpha rarefaction plots produced via q2-diversity) and the number of samples retained in the subset.
Alpha diversity metrics measure species richness and/or evenness within one sample.Shannon diversity index, observed features (ASVs), Faith's phylogenetic diversity (Faith PD) 28 , and evenness were the alpha diversity metrics measured, and differences in alpha diversity between the TN and HS groups of the four lines were analyzed using the non-parametric Kruskal-Wallis test.Beta diversity metrics, in particular unweighted and weighted UniFrac distance 29 , were used to analyze similarity or dissimilarity between microbiota of samples while considering phylogeny.Presence and absence of ASVs in samples is considered by unweighted UniFrac analysis, while the abundance of ASVs is considered in weighted UniFrac analysis 30 .For statistical analysis of UniFrac distances, the non-parametric permutational analysis of variance (PERMANOVA) test was used.Principal coordinates analysis (PCoA) was used to visualize distances between microbiota, where clustering of points may indicate similarities or differences between microbiota within treatment groups.PCoA results from QIIME 2 were imported to R 4.1.2 31using the package QIIME2R 0.99.35 32 .Within the tidyverse 1.3.0 33package, dplyr was used to select PCoA axes (PC1 and PC2) and join metadata, and ggplot2 was used to produce alpha diversity box plots and PCoA scatter plots.
The linear discriminant analysis effect size (LEfSe) algorithm 34 was used to analyze differential abundance of taxa up to the genus level between TN and HS birds in each line.The Huttenhower Lab Galaxy web server 35 was used to perform LEfSe analyses, using the default parameters.Phylogenetic Investigation of Communities by Reconstruction of Unobserved States 2 (PICRUSt2) version 2.4.2 software 36 was used to predict functional abundances based on marker gene sequences using the MetaCyc Metabolic Pathways Database 37 .PICRUSt2 output counts were transformed with center log-ratio transformation and inputted to STAMP 2.1.3 38to analyze and visualize predicted functional differences between TN and HS birds for each line, as well as compare all birds of the L2015 and JF lines.

Ethics approval and consent to participate
The study was conducted in accordance with the recommendations in the guide for the care and use of laboratory animals of the National Institutes of Health and the protocols were approved by the University of Arkansas Animal Care and Use Committee under protocols 18083 and 16084.All methods were carried out in accordance with relevant guidelines and regulations.This study was performed and reported in accordance with ARRIVE guidelines (https:// arriv eguid elines.org/).

Effects of heat stress on alpha diversity
Sequencing summaries of the CeL and CeM datasets are presented in Table 1.In CeL microbiota, Shannon diversity was significantly different based on groups (temperature condition and genetic line) overall (Kruskal-Wallis, H = 15.47,P = 0.03, Fig. 1A, and Additional file 1: Table S1), however, in pairwise comparisons between TN and HS groups for each genetic line, there were no significant differences (P > 0.05), though there was a trend of difference between TN and HS groups for the JF line (P = 0.055).Observed features were significantly different based on group (H = 39.81,P < 0.01, Fig. 1B, and Additional file 1: Table S2

Differences in alpha diversity in genetic lines
Pairwise comparisons were utilized where alpha diversity was significantly different based on groups to compare different genetic lines in the same temperature conditions.In CeL microbiota, Shannon diversity of the JF-TN group was lower compared to L1995-TN (H = 8.31, P < 0.01) and L2015-TN (H = 8.31, P < 0.01) and did not differ from ACRB-TN (P > 0.05, Fig. 1A, and Additional file 1: Table S1).Likewise, the ACRB-TN group was lower compared to L1995-TN (H = 4.03, P = 0.04) and L2015-TN (H = 5.63, P = 0.02).L1995-TN and L2015-TN did not differ in Shannon diversity (P > 0.05).Shannon diversity did not differ between genetic line under HS conditions (P > 0.05).Identical patterns were found in comparisons of TN groups for observed features and Faith PD (all P < 0.05, Fig. 1B,C, and Additional file 1: Tables S2 and S3), except for Faith PD only trending to be lower in JF-TN compared to L2015-TN (P = 0.055).JF-HS had greater observed features than L1995-HS and L2015-HS (both P < 0.05, Fig. 1B, and Additional file 1: Table S2), but it did not differ from ACRB-HS in observed features (P > 0.05).Observed features were different (P < 0.05) in all comparisons between ACRB-HS, L1995-HS, and L2015-HS, with observed features being lowest in L2015-HS and highest in ACRB-HS (Fig. 1B).JF-HS and ACRB-HS both had greater Faith PD than L1995-HS and L2015 HS (all P < 0.05), but Faith PD did www.nature.com/scientificreports/not differ between JF-HS and ACRB-HS (P > 0.05, Fig. 1C, and Additional file 1: Table S3).Faith PD did not differ between L1995-HS and L2015-HS (P > 0.05).Evenness was lower in JF-TN compared to L1995-TN and L2015-TN (both P < 0.05, Fig. 1D, and Additional file 1: Table S4) and did not differ from ACRB-TN (P > 0.05).ACRB-TN, L1995-TN, and L2015-TN did not differ in evenness from each other in all comparisons (P > 0.05).
Evenness was lower in JF-HS compared to L2015-HS (H = 6.53,P = 0.01), and there were no other significant comparisons between HS groups (all P > 0.05).

Effects of heat stress on beta diversity
CeL microbiota of different groups were distinct based on unweighted UniFrac distance (PERMANOVA, P < 0.01, visualized by PCoA in Fig. 2A, and statistics in Additional file 1: Table S6).This result included significant differences for comparisons between TN and HS birds for each genetic line: JF (P = 0.04), ACRB (P = 0.02), L1995 (P < 0.01), and L2015 (P < 0.01).When based on weighted UniFrac distance, microbiota of different groups were also distinct (P = 0.01, Fig. 2B, and Additional file 1: Table S7), however, only the JF line had significantly different microbiota (P = 0.02) between TN and HS birds.There were no significant differences (P > 0.05) between microbiota of TN and HS birds for the other three lines, though there was a trend of difference in the ACRB line (P = 0.06).
CeM microbiota of different groups were distinct (P < 0.01, Fig. 3A, and Additional file 1: Table S8) based on unweighted UniFrac distance.In comparisons between TN and HS birds for each genetic line, there were significant differences in microbiota for L2015 birds (P = 0.04), but no significant differences for the JF, ACRB, or L1995 lines (P > 0.05).According to weighted UniFrac distance, there were differences in CeM microbiota in groups overall (P = 0.03, Fig. 3B, and Additional file 1: Table S9), however there were no significant differences (P > 0.05) in comparisons of microbiota between TN and HS birds for each genetic line.

Beta diversity comparisons of genetic lines
Pairwise comparisons were utilized where beta diversity was significantly different based on groups to compare different genetic lines in the same temperature conditions.Based on unweighted UniFrac, the CeL microbiota of the JF-TN and ACRB-TN groups were distinct from L1995-TN and L2015-TN (all P < 0.05, Fig. 2A, and Additional file 1: Table S6), while JF-TN was not distinct from ACRB-TN and L1995-TN was not distinct from L2015-TN (both P > 0.05).ACRB-HS, L1995-HS, and L2015-HS were considered different to each other (all P < 0.05), while JF-HS was distinct from L1995-HS and L2015-HS, but not distinct from ACRB-HS.Based on weighted UniFrac, JF-TN were distinct from L1995-TN and L2015-TN (both P < 0.05, Fig. 2B, and Additional file 1: Table S7), but not from ACRB-TN (P > 0.05).ACRB-TN was distinct from L2015-TN (P = 0.03), but not from JF-TN or L1995-TN (P > 0.05), and L1995-TN and L2015-TN were not distinct (P > 0.05).Groups under HS conditions were not considered distinct from each other (all P > 0.05).
For CeM microbiota, JF-TN were not distinct from other TN groups based on unweighted UniFrac (all P > 0.05, Fig. 3A, and Additional file 1: Table S8).ACRB-TN was distinct from L1995-TN (P = 0.04), but not L2015-TN (P > 0.05), and L1995-TN and L2015-TN were not distinct (P > 0.05).JF-HS were not distinct from other HS groups (all P > 0.05), while ACRB-HS was distinct from L1995-HS (P = 0.01), but not distinct from www.nature.com/scientificreports/L2015-HS (P > 0.05).L1995-HS was distinct from L2015-HS (P = 0.01).Based on weighted UniFrac, all groups under the TN condition were not distinct from each other, and all groups under the HS condition were not distinct from each other (all P > 0.05, Fig. 3B, and Additional file 1: Table S9).

Differential bacterial abundance in heat-stressed and thermoneutral birds
The bacterial profiles of CeL and CeM microbiota by group (line and condition) are displayed in Fig. 4A,B, respectively.In CeL microbiota, the top five genera present in all samples were Bacteroides, Alistipes, Megamonas, Faecalibacterium, and Clostridia UCG-014 (Fig. 4A).In CeM microbiota, the top five genera present in all samples were Bacteroides, Helicobacter, Faecalibacterium, unclassified Lachnospiraceae, and Megamonas (Fig. 4B).
In CeM microbiota of JF, HS decreased the relative abundance of the of glycol metabolism and degradation, while it increased coenzyme M biosynthesis I (Fig. 8A).In ACRB, HS increased the relative abundance of the superpathway of sulfur oxidation (Acidianus ambivalens) (Fig. 8B).In L1995, HS decreased the relative abundance of the pathway l-lysine fermentation to acetate and butanoate (Fig. 8C).In L2015, HS decreased the relative abundance of 22 pathways and increased the relative abundance of 22 pathways (top 21 results by P-value displayed in Fig. 8D).Decreased pathways included l-lysine biosynthesis II, toluene degradation I and II, mevalonate pathway I, and superpathway of geranylgeranyldiphosphate biosynthesis I (via mevalonate).Increased pathways included mycolyl-arabinogalactan-peptidoglycan complex biosynthesis, superpathway of salicylate degradation, 4-methylcatechol degradation (ortho cleavage), aromatic compounds degradation via beta-ketoadipate, and catechol degradation III (ortho cleavage).
All L2015 birds were compared with JF birds to compare predicted functional abundance between a modern line and an ancestral line.In CeL microbiota, 18 pathways were in greater relative abundance in L2015 birds, and 13 pathways were in greater relative abundance in JF birds (Fig. 7E).Pathways in greater abundance in L2015 included ppGpp biosynthesis, protein N-glycosylation (bacterial), superpathway of demethylmenaquinol-6 biosynthesis II, superpathway of sulfolactate degradation, and l-methionine salvage cycle III.Pathways in greater abundance in JF included superpathway of 2,3-butanediol biosynthesis, isopropanol biosynthesis, superpathway of glycerol degradation to 1,3-propanediol, glycerol degradation to butanol, and mevalonate pathway I.In CeM microbiota, 17 pathways were in greater relative abundance in L2015, and 22 pathways were in greater relative abundance in JF birds (Fig. 8E).Pathways in greater abundance in L2015 included peptidoglycan biosynthesis V (beta; -lactam resistance), glycerol degradation to butanol, mevalonate pathway I, superpathway of geranylgeranyldiphosphate biosynthesis I (via mevalonate), and superpathway of glycerol degradation to 1,3-propanediol.

Discussion
The objectives of this study were to understand the effect of HS on factors such as alpha and beta diversity, bacterial abundance, and predicted metabolic function in the CeL and CeM microbiota, in addition to comparing these factors between four genetic lines of broilers and determining whether HS affected the microbiota differently in a line-dependent manner.Our findings demonstrate that HS had more prominent effects on the CeL microbiota compared to CeM and that the type of effects observed from HS were heavily dependent on genetic line.Lastly, predictive functional analysis showed that further studies in metabolic function of modern broilers may be of interest, as HS was predicted to have greater effects on metabolic pathways in modern L2015 compared to other lines.These findings are of importance in understanding both the roles of heat stress and host genetics on the microbiota, informing further research on development of effective microbiota-altering strategies (e.g.probiotics) to mitigate the negative effects of environmental stressors such as heat stress.
Our alpha diversity results showed that HS affected richness (observed features) in the CeL microbiota of all lines, however, the effects were in opposite directions in the slow-growing JF and ACRB lines compared to the moderate-growing L1995 and modern fast-growing L2015.Interestingly, richness increased in the two slowgrowing lines, while it decreased in the two faster-growing lines.Results followed the same pattern for another richness metric, Faith PD, except for HS having only a trend of increasing Faith PD in the JF line.These results suggest that potential microbiota-altering solutions would need to consider differences in genetic lines, including physiological differences and bacterial composition.It has been previously shown that HS increases core body temperature in the more modern L1995 and L2015 lines, but not in the JF or ACRB lines 3 , which could be one explanation for the opposing effects, as bacteria would be responding to different temperature conditions.Alternatively, genetic differences in lines have been found to influence intestinal microbiota composition 39,40 , therefore the differences in composition could lead to different responses to HS and contribute to the differential  changes in core temperature.Our data showed richness, as well as Shannon diversity, were significantly lower in the CeL microbiota of TN birds in the JF and ACRB lines compared to L1995 and L2015.These different starting points in diversity may influence competition dynamics of bacteria in the slower and faster-growing lines, affecting how richness changes during disruption by HS.Alpha diversity in CeM microbiota was not significantly affected by HS, and this result again demonstrates the differences between luminal and mucosal microbiota and the importance of understanding bacterial profiles of both regions, as they can be affected by stressors differently.Differing effects on luminal and mucosal microbiota were also observed in a previous experiment on the effect of HS on ileal microbiota, where we found HS had more prominent effects on Shannon diversity and taxonomic abundance in the ileal mucosa compared to the lumen 19 .
Beta diversity results supported alpha diversity findings showing that HS primarily altered CeL microbiota.Specifically, bacterial profiles tended to be affected by the presence and absence of bacterial taxa (unweighted UniFrac) as opposed to the abundances of bacterial taxa (weighted UniFrac).In all four lines, CeL microbiota were considered distinct between TN and HS birds in unweighted UniFrac analysis, while weighted UniFrac analysis determined CeL microbiota differed in HS birds only in the JF line.This was consistent with the alpha diversity results, as richness (presence and absence) was affected by HS.From visualization by PCoA, a puzzling phenomenon is observed where CeL bacterial profiles of faster-growing L1995 and L2015 birds under HS became more similar to TN bacterial profiles of slower-growing JF and ACRB birds.Meanwhile, CeL bacterial profiles of JF and ACRB birds under HS became more similar to TN bacterial profiles of L1995 and L2015 birds.CeM microbiota again displayed different patterns, where TN and HS microbiota differed only in L2015 birds, based on unweighted UniFrac.Future research to understand the underlying physiological and functional differences in different genetic lines may be of importance in developing probiotics or supplementation solutions that are suitable for the microbiota of specific breeds.
Although bacterial profiles primarily differed from presence and absence of taxa rather than abundance under HS, we could still analyze differential abundance in particular bacterial taxa at the genus level.Similar to how HS appeared to result in decreased richness in the CeL microbiota of more modern lines, HS appeared to decrease the relative abundance of certain genera in modern lines.Except for Parasutterella in L2015-HS birds, all other differentially abundant taxa were in greater relative abundance in the TN groups of the faster-growing L1995 and L2015 lines.In the slower-growing ACRB and JF lines, we observed the opposite pattern, where the majority of differentially abundant taxa were those that increased in relative abundance in the HS groups, while only five genera in JF birds were greater in the TN group.These results again demonstrate that genetic line plays a large role in how HS affects the CeL microbiota, as overall richness and abundance of certain taxa can be affected in opposite ways in different lines.In the L1995 and L2015 lines, genera within the orders Lachnospirales and Oscillospirales, which are typically obligate anaerobes, were the most likely to decrease in relative abundance from HS.The intestinal lumen is anoxic under healthy conditions, however, the oxidative stress that can occur during HS can affect the oxygen gradient, shifting the microbiota balance from obligate anaerobes to facultative anaerobes 41 .This could explain the decrease of these genera in HS birds of L1995 and L2015, though the dynamics may differ in JF and ACRB.Dynamics such as these may make certain groups of taxa more appropriate for inclusion in microbiota-altering solutions, or indicate the need to supplement metabolites to account for the reduced abundance of taxa which may be responsible in producing those metabolites.
While decreases in genera were the norm in L2015, Parasutterella was significantly higher under HS in both CeL and CeM microbiota of L2015, similar to results in another study on ileal mucosal microbiota of the same line 19 .Parasutterella is an obligate anaerobe of the phylum Proteobacteria, and potentially plays a role in bile acid maintenance and cholesterol metabolism due to its correlations with microbial-derived metabolites such as aromatic amino acid, bilirubin, purine, and bile acid derivatives 42 .Unlike our studies, Parasutterella has been positively correlated with performance in the CeL microbiota 13 , though this occurred due to ammonia exposure as opposed to HS.While HS tended to decrease relative abundance of other genera in CeL microbiota of L2015 birds, Parasutterella was the only genus in greater relative abundance, possibly indicating a unique interaction under HS.In contrast, Streptococcus and Tyzzerella decreased in relative abundance under HS in both the CeL and CeM microbiota of L2015 birds, and the same was true for Akkermansia in L1995 birds.Streptococcus has also been seen to decrease in the ileal microbiota of Arbor Acres broilers 6 and has been negatively correlated with broiler body weight in the ileum 12 , however, further investigation is required as some Streptococcus species may function as normal gut flora while others cause disease.In the same study, Akkermansia was found to negatively correlate with body weight in both the cecum and ileum 12 , however, the abundance of Akkermansia in our study was generally low and unlikely to have a significant functional impact.Unlike in the ileal mucosa of our previous study and in the CeL of Arbor Acres broilers in Liu et al. 43 , Tyzzerella instead decreased under HS.This genus has been linked to cardiovascular disease in humans 44 .As relative abundance of taxa may be affected by HS, future studies may seek to link changes in abundance of taxa to changes in metabolic function.To predict metabolic functional abundances, we utilized functional analysis.Although accuracy of functional analysis can be limited outside of human studies 45 , the tool provides a cost-free method of predicting functional abundances from 16S rRNA sequencing, which may assist in informing further research utilizing shotgun metagenome sequencing.Function of the gut microbiota is hypothesized to differ in modern broilers compared to ancestral broilers, as function is a potential link between the differences in bacterial composition in broiler lines observed in this study and the greater susceptibility of modern broilers to HS reported previously 3 .Our results from comparison of predicted function in CeL microbiota showed genes predicted to influence ppGpp biosynthesis were in greater relative abundance in the JF line.ppGpp biosynthesis is involved in the bacterial stringent response to environmental stresses, where a switch in transcription profile occurs to rapidly produce factors necessary for stress resistance, glycolysis, and amino acid synthesis 46 .Other pathways more frequent in JF included protein N-glycosylation and l-methionine salvage cycle III, which may indicate the CeL microbiota of the JF line has a strong response to environmental stressors.Multiple pathways related to production of alcohols, such as butanediol, isopropanol, and butanol, were in greater abundance in both the CeL and CeM microbiota of L2015 birds.
In addition to comparison of L2015 and JF birds, HS and TN birds were compared within each line.Overall, HS had more prominent effects on the predicted functional abundances in L2015 birds in both the CeL and CeM microbiota.Very few changes occurred in ACRB and L1995 birds, while there were some changes in the CeL microbiota of JF birds such as decreased relative abundance of pathways related to heme biosynthesis, TCA cycle IV, and usage of amino acids like tryptophan and methionine.Numerous predicted pathways were differentially abundant in L2015 birds, totaling 40 pathways in CeL and 44 pathways in CeM.In CeL microbiota, multiple ubiquinol biosynthesis pathways were in increased abundance under HS, which may be a response to deal with overproduction of free radicals under stressful conditions 47 .This prediction was also observed in the CeM microbiota of Eimeria tenella-infected broilers 11 .Unlike in JF, pathways related to heme biosynthesis were increased under HS.Lactose and galactose degradation I was decreased under HS in both CeL and CeM, which may have implications as prebiotics often contain lactose, galactose, and/or related sugars 9 .

Conclusions
HS affected the bacterial richness and resulted in distinct bacterial profiles in the cecal luminal microbiota in four genetic lines, decreasing richness in more modern lines and increasing richness in older lines.HS had limited effects in cecal mucosal microbiota, only resulting in distinct bacterial profiles in L2015.Differences in bacterial profiles in different genetic lines appeared to influence patterns of differential taxonomic abundance, potentially affecting likelihood of increases or decreases in certain taxa.Predicted functional abundances suggested that L2015 broilers may have functional differences to the ancestral JF line, and selection leading to modern L2015 broilers may have made modern broilers more susceptible to functional changes under HS.Continued research on GIT microbiota and the function of its members is of importance to alleviate the negative effects of HS on modern broiler chickens' performance and health.

Figure 2 .
Figure 2. Principal coordinate analysis (PCoA) of CeL microbiota based on: (A) unweighted UniFrac and (B) weighted UniFrac.Colors indicate difference in lines, while shape of points indicate temperature condition.Longer distances between points indicate microbiota profiles were different, while shorter distances between points indicate profiles were similar.

Figure 3 .
Figure 3. Principal coordinate analysis (PCoA) of CeM microbiota based on: (A) unweighted UniFrac and (B) weighted UniFrac.Colors indicate difference in lines, while shape of points indicate temperature condition.Longer distances between points indicate microbiota profiles were different, while shorter distances between points indicate profiles were similar.

Figure 4 .
Figure 4. Relative abundances of bacterial taxa at the genus level in groups (line and condition; n = 6 samples for each column) of (A) CeL microbiota and (B) CeM microbiota.The 20 most abundant taxa overall are listed in the legend.

Figure 5 .
Figure 5. Linear discriminant analysis effect size (LEfSe) in comparisons of CeL microbiota in HS and TN groups for each line: (A) Jungle Fowl (JF), (B) Athens Canadian Random Bred (ACRB), (C) 1995 Random Bred (L1995), (D) 2015 Modern Random Bred (L2015).Positive effect size indicates higher relative abundance in HS group, while negative effect size indicates higher relative abundance in TN group.

Figure 6 .
Figure 6.Linear discriminant analysis effect size (LEfSe) in comparisons of CeM microbiota in HS and TN groups for each line: (A) Jungle Fowl (JF), (B) Athens Canadian Random Bred (ACRB), (C) 1995 Random Bred (L1995), (D) 2015 Modern Random Bred (L2015).Positive effect size indicates higher relative abundance in HS group, while negative effect size indicates higher relative abundance in TN group.

Figure 7 .
Figure 7. Effect of HS on the mean proportion (%) of predicted MetaCyc pathways (up to the top 21 shown) in the CeL microbiota of each line: (A) Jungle Fowl (JF), (B) Athens Canadian Random Bred (ACRB), (C) 1995 Random Bred (L1995), (D) 2015 Modern Random Bred (L2015).In (E), L2015 and JF were compared to investigate the differences in a modern line and an ancestral line.

Figure 8 .
Figure 8.Effect of HS on the mean proportion (%) of predicted MetaCyc pathways (up to the top 21 shown) in the CeM microbiota of each line: (A) Jungle Fowl (JF), (B) Athens Canadian Random Bred (ACRB), (C) 1995 Random Bred (L1995), (D) 2015 Modern Random Bred (L2015).In (E), L2015 and JF were compared to investigate the differences in a modern line and an ancestral line.

), including significantly lower observedTable 1 .
Sequencing summary of CeL and CeM microbiota datasets processed in QIIME 2. QC quality control via DADA2, ASVs amplicon sequence variants.Reads after filtering indicates the number of reads after exclusion of mitochondria, chloroplasts, and unassigned bacteria.