
==== Front
Commun Biol
Commun Biol
Communications Biology
2399-3642
Nature Publishing Group UK London

39300162
6815
10.1038/s42003-024-06815-0
Article
The gut microbiota and its metabolite butyrate shape metabolism and antiviral immunity along the gut-lung axis in the chicken
Saint-Martin Vincent 1
Guillory Vanaique 1
Chollot Mélanie 1
Fleurot Isabelle 1
Kut Emmanuel 1
Roesch Ferdinand 1
Caballero Ignacio 1
http://orcid.org/0000-0001-5886-4205
Helloin Emmanuelle 1
Chambellon Emilie 1
http://orcid.org/0000-0002-6873-1032
Ferguson Brian 2
Velge Philippe 1
Kempf Florent 1
http://orcid.org/0000-0002-2735-0217
Trapp Sascha 1
http://orcid.org/0000-0003-4005-1753
Guabiraba Rodrigo rodrigo.guabiraba-brito@inrae.fr

1
1 https://ror.org/02wwzvj46 grid.12366.30 0000 0001 2182 6141 INRAE, ISP, Université de Tours, 37380 Nouzilly, France
2 https://ror.org/013meh722 grid.5335.0 0000 0001 2188 5934 Department of Pathology, University of Cambridge, Cambridge, United Kingdom
20 9 2024
20 9 2024
2024
7 118530 11 2023
30 8 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by/4.0/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
The gut microbiota exerts profound influence on poultry immunity and metabolism through mechanisms that yet need to be elucidated. Here we used conventional and germ-free chickens to explore the influence of the gut microbiota on transcriptomic and metabolic signatures along the gut-lung axis in poultry. Our results demonstrated a differential regulation of certain metabolites and genes associated with innate immunity and metabolism in peripheral tissues of germ-free birds. Furthermore, we evidenced the gut microbiota’s capacity to regulate mucosal immunity in the chicken lung during avian influenza virus infection. Finally, by fine-analysing the antiviral pathways triggered by the short-chain fatty acid (SCFA) butyrate in chicken respiratory epithelial cells, we found that it regulates interferon-stimulated genes (ISGs), notably OASL, via the transcription factor Sp1. These findings emphasize the pivotal role of the gut microbiota and its metabolites in shaping homeostasis and immunity in poultry, offering crucial insights into the mechanisms governing the communication between the gut and lungs in birds.

Comparison of conventional and germ-free chicken show gut microbiota's significant impact on immunity and metabolism. Key findings include its role in lung mucosal immunity during avian influenza and butyrate's influence on antiviral pathways

Subject terms

Innate immunity
Microbiome
Agriculture
This work was funded by the INRAE, France (specifically Santé Animale – SA and Microbiologie et Chaine Alimentaire – MICA divisions), by "EuroFéRi" (FEDER-FSE Centre Val de Loire 2014-2020, N° EX 010233) and by VetBioNet, EU H2020 (Grant Agreement N° 731014).issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

The gut microbiota (GM) is a complex ecosystem that coexists in a symbiotic relationship with the host. These microorganisms constituting the GM actively participate in a broad range of physiological functions, including digestion, nutrient absorption, and regulation of the immune system1,2. Disruptions in the GM composition have been shown to impact the host’s health in multiple ways, ranging from metabolic disorders to immune system dysfunctions3. The GM also modulates the host immune activity, enhancing its function against pathogens while limiting its reactivity to self-antigens1. One of the key ways in which the GM modulates the host’s physiology is via secreted molecules such as bacterial metabolites, including short-chain fatty acids (SCFAs). SCFAs are the end-products of bacterial fermentation of dietary fibre, and can play a crucial role in regulating immunity and metabolism in the gut and peripheral tissues4. SCFAs can also improve intestinal barrier function, limiting the translocation of pathogenic bacteria and their toxins into the bloodstream5.

Recent research suggests that SCFAs have broad, extra-intestinal, and long-lasting biological activity beyond the gut4,6. They have specific receptors in other tissues and organs, indicating their relevance to the regulation of whole-body physiology. The periphery, including the liver, adipose tissue, lungs, and brain, is particularly sensitive to SCFAs and other microbial metabolites, leading to the coining of terms such as “gut-lung axis” and “gut-brain axis”7,8. Therefore, the regulation of microbiota composition has well-recognised consequences for immune development and overall health in humans and mammalian model species.

Poultry farming is a central pillar of global animal protein production, and poultry health is paramount to the success of this industry. The GM plays a crucial role in chicken health and well-being9,10, and its acquisition and development are particularly important in early life11–13. After hatching, the microbiota in the chicken’s gut begins to develop, with colonisation by bacteria such as Lactobacillus and Bifidobacterium14. The GM is also seen as essential for the development of the immune system in poultry, and studies have shown that an unregulated or poorly enriched microbiota in early life can retard immune system development and functioning15. Moreover, we have previously shown that GM can influence the immune response in the chicken respiratory tract, suggesting the functional existence of a gut-lung axis in the chicken16. The chicken respiratory tract is particularly sensitive to viral infections, especially those caused by avian influenza viruses (AIV), a significant cause of morbidity and mortality in commercial and backyard poultry17. Understanding the mechanisms through which the GM and SCFAs influence metabolism and immunity in the respiratory tract and other peripheral tissues is an area of active research and would allow the development of probiotics/postbiotics with optimised capacity to enhance immune system functions beyond the gut. Nevertheless, the gut-lung axis in birds is largely unexplored, despite most pathogens in commercial poultry targeting the gut and/or airways.

To define how the chicken GM regulates metabolism and immunity in the gastrointestinal tract and in the periphery, we analysed transcriptomic and metabolic signatures in conventional (CV) and germ-free (GF) chickens. GF birds exhibited significant physiological changes in the caecal compartment, lacking central GM-derived metabolites like SCFAs. SCFA detection in peripheral compartments of CV birds and differential expression of innate immunity genes in GF birds’ peripheral tissues further corroborated the functional existence of a gut-lung axis in birds. In an AIV H7N1 challenge, GF chickens showed different immunoregulatory responses to infection, though clinical outcomes were similar. Finally, butyrate, an SCFA, displayed antiviral effects in chicken respiratory epithelial cells, likely via the regulation of interferon-stimulated genes, particularly OASL. The observed antiviral signalling mechanism would involve HDAC inhibition and Sp1-dependent regulation of the OASL promoter, shedding light on new mechanisms through which the GM regulates poultry mucosal immunity along the gut-lung axis in the chicken.

Methods

Animals

GF chickens were produced following an established protocol18. Eggs were obtained from a histocompatible B21/B21 White Leghorn layer chicken line (PA12), sourced from SPF hens at INRAE, Plateforme d’Infectiologie Expérimentale (PFIE, Nouzilly, France). Following surface disinfection, the eggs were incubated in a sterile hatching incubator. Following hatch, GF chicks were transferred to sterile isolators in BSL2 containment. For CV birds, hatched under SPF conditions without specific disinfection, a similar isolator environment was provided from day 1 post-hatching. Both groups received ad libitum access to a commercial diet (Safe 115, Safe-Diets, France). GF bird feed was sterilised by gamma irradiation, and tap water (autoclaved for GF birds) was provided through water dispensers. These conditions were maintained for 21 days. Sterility tests, a combination of molecular methods (broad-range 16S ribosomal RNA gene polymerase chain reaction (PCR)19) and microbiological methods (thioglycolate and Brain Heart Infusion (BHI) broths to assess the absence of bacteria and fungi up to 14 days after sampling)18, were routinely performed on pooled fecal contents from birds within a given isolator. As for the microbiological analyses, 1 mL of stool was mixed with 9 mL of thioglycolate broth with resazurin, while the remaining sample was mixed with 9 mL of BHI. The tubes were then incubated at 37 °C without shaking for 18–48 hours, favouring the growth of aerobic, facultative aerobic, and non-fastidious anaerobic species, including bacteria, fungi, and yeast. Thioglycolate broth with resazurin, designed for detecting non-fastidious anaerobic bacteria, also facilitated the detection of aerobic bacteria, yeast, and fungi, in compliance with the sterility testing standards of the European, American, and Japanese pharmacopoeias20–22. Perturbations in the growth media were visually examined after 18 hours of incubation. After 48 hours, a drop from the BHI fecal-broth media was placed on a glass slide and examined under a microscope (×40 magnification) for the presence or absence of bacteria, yeast, or fungi. If suspicion arose regarding the presence of these microorganisms, a sample from the BHI culture was inoculated onto BHI agar plates and incubated at 37 °C for 18 to 48 hours for further analysis. Swabs from isolators were periodically tested during hatching and post-hatching periods until the end of the experiment, employing the same detection methods for bacteria and fungi described above. All experiments were conducted in compliance with the relevant European and national regulations and authorised by the local ethics committee (Comité d’Ethique en Expérimentation Animale Val de Loire) under the reference number APAFIS#30655-2021042114338110 v1.

Tissue sampling

At day 21 post-hatching, CV and GF birds (n = 8 per group) were sacrificed by cervical dislocation followed by blood withdrawal from the occipital venous plexus. Blood samples were stored in heparin-coated tubes for subsequent plasma preparation. Caecal contents, caeca, spleens, and lungs were immediately collected from each chicken, snap-frozen in liquid nitrogen, and stored at −80 °C until further use for metabolomics analysis and/or RNA and DNA extraction and sequencing.

Avian influenza virus infection model

The low-pathogenic avian influenza virus (LPAIV) strain A/Turkey/Italy/977/1999 (H7N1)23, which causes moderate to severe disease in experimentally infected chickens, including interstitial pneumonia and intestinal infection24,25, was used in the AIV challenge infection experiment. Virus stocks were prepared by one round of amplification in 10-day-old embryonated chicken eggs at a multiplicity of infection (M.O.I.) of 2 plaque-forming units (PFU) per egg. Three-week-old immunocompetent/immunoresponsive26 GF or CV chickens were housed in BSL2 cabinets under negative pressure with HEPA-filtered air. Briefly, a group of 24 GF and 24 CV birds were equally inoculated intra-tracheally (i.t.) with 0.2 ml of 106 EID550 of H7N1, and 15 GF and 15 CV control (mock) birds were inoculated i.t. with 0.2 ml of PBS. At days 1, 2, and 3 p.i., caecal contents, tissues, and plasma were recovered from a group of animals (n = 5–8) for different downstream analyses. Birds were carefully monitored (twice daily) during the course of infection, and diseased birds presenting at least 4 of the following symptoms for 24/48 hours were killed by intravenous pentobarbital injection (humane endpoint): ruffled feathers, tremors, coughing/sneezing, diarrhoea, lethargy or anorexia. Clinical signs were evaluated according to the following score: 0 (no clinical signs), 1 (mild clinical signs), 2 (temporary severe clinical signs), or 3 (dead or euthanized due to persistent severe clinical signs). A scoring system was used to evaluate macroscopic lung lesions as follows: 1 (mild, localised oedema and fibrinous exudate), 2 (moderate oedema and with haemorrhage and fibrinous exudate over ~1/4 of the lung), or 3 (severe haemorrhage and extensive oedema over ~1/2 of the lung).

Cell culture and in vitro treatment and infection

The spontaneously immortalised chicken respiratory epithelial cell line CLEC213 has previously been obtained from the lungs of histocompatible B21/B21 white leghorn layer chickens (PA12)27. Cells were cultured in Dulbecco’s Modified Eagles’s Medium (DMEM) supplemented with 10% fetal calf serum (FCS) (Gibco, UK), 2 mM L-Glutamine (Gibco, UK), and 100 UI/ml penicillin and 100 μg/ml streptomycin (PS, Gibco, UK) at 41 °C in a 5% CO2 incubator. For SCFA biological effects, cells were plated at 5 × 106 cells/ml in 12-well plates, incubated with acetate, propionate, or butyrate (sodium salt form, Sigma-Aldrich-UK) at varying concentrations and times. Some experiments included chicken interferon alpha (Yeast-derived, Kingfisher Biotech, USA) and Mithramycin A (Merck Millipore, USA). Cell lysates and RNA were analysed for gene expression and metabolomics. Cell viability and counts were assessed using Spark® Cyto (Tecan, Switzerland) after Trypan Blue (Sigma-Aldrich, UK) addition.

AIV infection experiments using CLEC213 were performed to assess the effects of butyrate pre-treatment for 16 h on viral replication and innate immune gene regulation. To this end, subconfluent CLEC213 monolayers (3 × 105 cells/well in six-wells plates) were pre-treated with butyrate (3 mM) for 16 h and then infected with the AIV strains A/Mallard/Marquenterre/Z237/83 (H1N1)28 or A/Turkey/Italy/977/1999 (H7N1)23 at a MOI of 0.1 or 0.01 by incubating the cells with a virus dilution prepared in serum-free DMEM for 1 h at 37 °C. Following this 1 h adsorption period, cells were washed with PBS and then overlayed with serum-free DMEM containing TPCK-treated trypsin (0.4 μg/ml final). Viral replication was assessed in supernatants or cell lysates (RNA) from 12 to 48 h post-inoculation via titration of supernatants on MDCK cells or by qRT-PCR targeting the M1 coding sequence of gene segment 729. Primers M52C and M253R span conserved sequences in gene segment 7 of influenza A virus and have no homology to nucleotide sequences from other species available from GenBank. Viral load quantification was expressed as PFU/mL or M1 genomic copy numbers. Cellular RNA extractions from infection experiments with the H1N1 and H7N1 strains were also used for gene expression analysis.

HDAC activity assay

The histone deacetylase (HDAC) inhibitory properties of butyrate are extensively documented in mammals30,31. The HDAC-Glo™ I/II assay kit (Promega, UK) was used to measure class I and II HDAC activity in CLEC213 cells according to the manufacture’s protocol. Briefly, the acetylated peptide substrate was added to the cell culture medium in the presence or absence of butyrate, and peptide cleavage requiring deacetylated protease activity was quantified by measuring the release of luminescent aminoluciferin using a GloMax plate reader (Promega, UK). Reduction in relative light units (RLU) correlated with decreased HDAC activity due to the inhibitory properties of the tested molecule. Trichostatin A served as a positive control. Data are presented as RLU or as percentage of HDAC activity, with untreated cells as the reference (100%).

Metabolomics analysis

Caecal contents, tissues, plasma and cell lysates were specifically prepared for optimised SCFA quantification through proton nuclear magnetic resonance (H1-NMR) according to Hauser et al.32. Caecal contents samples (20–50 mg) were homogenised in 1.2 ml of phosphate buffer (0.2 M, pH 7) and a CDCl3/CD3OD (2:1, v/v) solution. After centrifugation (5500 × g, 10 min, 4 °C) twice, the aqueous phases were collected, dried using a SpeedVac Vacuum Concentrator, reconstituted in phosphate buffer (0.2 M, pH 7), vortexed, and centrifuged (5500 × g, 15 min, 4 °C). Subsequently, 500 µl of supernatant were transferred into 5 mm NMR tubes, and 100 µl of TSP (1 mM) were added. Tissues (100–200 mg), plasma, and cell lysates were homogenised in 4.85 ml of methanol/water (80/20, v/v) per gram of tissue for 40 s. Next, 2 ml of dichloromethane per gram of tissue were added, and the sample was vortexed (5 s). Finally, 2 ml of dichloromethane and 2 ml of water per gram of tissue were added, the sample was vortexed (10 s), and kept standing at 4 °C for 15 min. After centrifugation (5500 × g, 15 min, 4 °C), the aqueous phases were collected. Samples were dried using a SpeedVac Vacuum Concentrator, reconstituted in 200 µl of phosphate buffer (0.2 M, pH 7), vortexed, and centrifuged (5500 × g, 15 min, 4 °C). Then, 150 µL of supernatant were transferred into 3 mm NMR tubes, and 50 µL of TSP solution (1 mM) were added to the NMR tube. 1H NMR spectra were obtained at 300 K on a Bruker Avance III HD 600 MHz NMR spectrometer. The “noesypr1d” pulse sequence was used for water signal suppression with a mixing time of 100 ms. A total of 1024 and 256 transients were collected for tissue and caecal samples, respectively, into 64k data points using a spectral width of 12 ppm, a relaxation delay of 15 s, and an acquisition time of 4.5 s. Prior to Fourier transform, an exponential line broadening function of 0.3 Hz was applied to the FID. All NMR spectra were phase- and baseline-corrected and referenced to the chemical shift of TSP (0 ppm) using Topspin (V3.2, Bruker Biospin, Germany). Metabolite concentrations were calculated using the TSP signal at 0 ppm integrating for 9 protons and with known concentrations of different metabolites.

Gene expression analysis

Total RNA was extracted from cells or tissues using the NucleoSpin RNA II kit (Macherey-Nagel, Germany) with rDNAse treatment to eliminate contaminant DNA. RNA quality and concentration were assessed with a NanoDrop (Thermo Scientific, USA). Reverse transcription of 0.5–1 μg total RNA per reaction was performed using the iScript cDNA synthesis kit (Bio-Rad, USA). Quantitative real-time PCR (qPCR) on a CFX96 machine (Bio-Rad, USA) used iQ SYBR Green Supermix, cDNA, primers (250 nM, Eurogentec, Belgium), and nuclease-free water (Sigma-Aldrich, UK) in a 20 μL reaction volume. CFX Manager software 3.1 (Bio-Rad, USA) analysed the qPCR data. Amplicon size was confirmed by 2% agarose gel. Target gene expression was normalised to chicken glyceraldehyde 3-phosphate dehydrogenase and β-2-microglobulin. Relative normalised expression represented as log2 fold change compared to the CV chicken group was calculated using the 2ΔΔCt method. Primer pairs are listed in Supplementary Table 1.

High-throughput qPCR using the 48.48 Dynamic Array

Quantitative high-throughput qPCR was conducted on the BioMark HD instrument using a 48.48 Dynamic Array (Fluidigm). Total lung RNA (20 ng) was reverse transcribed using iSCRIPTtm Reverse Transcription Supermix (Bio-Rad). Pre-amplification used a 1:10 dilution of cDNA in water with Preamp Master Mix (Fluidigm, #100-5581) in a T100 Thermal Cycler (BioRad), following thermal cycling conditions. Exonuclease I-treated pre-amplified cDNAs were diluted 1:5 in TE. The qPCR employed primers, 2x SsoFast EvaGreen with Low Rox (BioRad 172-5211), and 20× DNA binding Dye (Fluidigm PN 100-7609) on a 48.48 Dynamic Array Integrated Fluidic Circuit (Fluidigm). Thermal cycling conditions were 1 min at 95 °C, followed by 30 cycles. The analysis software 3.1.3 (Fluidigm) was used with a quality threshold, linear baseline correction, and auto (global) quantitation cycle (Cq) threshold method. Data were analysed using the 2ΔΔCt method, expressing results as relative fold change (Fc) compared to a given control. Primer pairs are listed in Supplementary Table 1.

Knockdown of OASL and SP1 expression in CLEC213 pulmonary epithelial cells using siRNA

To investigate the roles played by OASL and SP1 in the cellular and molecular mechanisms underlying the biological effects of butyrate CLEC213 cells, we opted to knockdown (KD) their expression in these cells using RNA interference methods33,34. Optimal siRNA sequences for knocking-down OASL (GenBank: NM_001397447.1) and SP1 (Genbank: NM_204604.2) were determined using a tool available on the website of the siRNA manufacturer (Horizon Discovery, USA). The transfection involved adding 10 pmol of siRNA to the Lipofectamine 2000 (Thermo, USA) and OPTI-MEM 1× (Gibco, USA) mixture, vortexing, and incubating for 15 minutes at room temperature. Then, 50 µL of the mixture was added per well in a 24-wells plate or 200 µL per well for a six-wells plate, one hour before butyrate addition. siRNA efficiency was verified through qPCR by measuring SP1 and OASL expression. siRNA sequences used are listed in Supplementary Table 2.

RNA sequencing analyses

The total RNA from tissue and CLEC213 cells was extracted using the NucleoSpin RNA Plus, Mini kit for RNA purification with a DNA removal column (Macherey-Nagel, Germany). RNA quality was checked by the Bioanalyzer 2100 system using the RNA Nano 6000 Assay kit (Agilent Technologies, CA, USA). The NEBNext® Ultra™ II RNA Library Prep Kit for Illumina® was used to prepare the library. Raw reads were trimmed with TrimGalore (v0.4.1, parameters: –clip_R2 2)35. Trimmed reads were mapped and quantified using STAR (v2.6.1c) and RSEM (v1.3.1) using the function rsem-calculate-expression (parameters: –star –sort-bam-by-coordinate) and the reference file Ensembl annotation release GRCg6a, Ensemble annotation release 98, genome-build-accession NCBI:GCA_000002315.536. Differential gene expression was calculated using DESeq2 (v1.30.0)37. The Wald test was used to generate p values and log2 fold changes. Genes with an adjusted p value < 0.05 and absolute log2 fold change >1 were called as differentially expressed genes (DEGs) for each comparison. GO and KEGG pathways enrichment analysis were done using g:Profiler38. An FDR threshold of 0.05 was used for determining GO category overrepresentation. For significant enrichment, Benjamini-corrected value P ≤ 0.05 was considered. Gene ontology enrichment analysis was conducted using Ingenuity Pathway Analysis software (IPA, Qiagen, Germany). For such, a confidence-boosting maximum adjusted p value threshold was applied: 0.45 for caeca, which included the analyses of the top 244 genes from the caeca RNAseq analyses; 0.4 for lungs, which included the analyses of the top 280 genes from the lungs RNAseq analyses; and 0.1 for spleen, which included the analyses of the top 84 genes from the spleen RNAseq analyses. From these datasets, histograms were generated with GraphPad Prism software (version 6.0, GraphPad software, USA), illustrating the relative expression of a chosen set of genes. The complete list of DEG from RNAseq in CV and GF chickens at 21 days of age are shown in Supplementary Data 1. The complete list of DEG from RNAseq in CLEC213 chicken lung epithelial cells are shown in Supplementary Data 2. Finally, selected DEG from tissue RNAseq used for IPA analysis are listed in Supplementary Data 3.

DNA extraction and 16S sequencing

Total DNA was extracted from caecal contents using the Dneasy PowerSoil kit (Qiagen, UK) following the manufacturer’s protocol, and 16S (V3–V4 region) amplification was performed as described by Zheng et al.39. In the limited cycle 2nd round PCR, sample-specific barcodes were added for multiplexing. Pooled libraries underwent QC before Illumina 2 × 250bp paired-end sequencing for comprehensive V3 and V4 coverage. The FROGS analysis pipeline40 was used for clustering paired-end reads with a 0.1 mismatch rate, selecting reads based on an expected size of 300 bp and a total amplicon size of 400–600 bp (mean 460 bp). Swarm41 was employed for sequence clustering with aggregation distance parameters of 1 and 1 for denoising and final clustering steps. VSearch removed OTUs, including chimeric sequences (69), and quality control measures were applied, excluding rare OTUs (relative abundance < 0.0005% of total read numbers) and those matching phiX sequences in a specific data bank. OTUs were classified using NCBI BLAST+ search within the Silva SSU 123 database42. Diversity assessment included Chao1, Shannon α-diversity, and Bray–Curtis β-diversity indexes. Differential abundances were assessed using the R-package DESeq237, with significant logarithmic fold change ratios determined by Wald tests and Benjamini-Hochberg adjustment for multiple testing (P < 0.01).

Flow cytometry analysis

Spleen samples were obtained from both CV and GF chickens, then homogenised using 40-μm sterile strainers and resuspended in PBS supplemented with 2% FCS and 2 mm EDTA. Next, nucleated red blood cells were partly eliminated using the ammonium chloride-based RBC lysis buffer Hybri-Max™ (Sigma-Aldrich, UK). The resulting leucocytes were carefully layered onto Histopaque 1077 (Sigma-Aldrich, UK) and centrifuged for 30 minutes at 750 g without braking, at room temperature, for enhanced purification of leucocytes and elimination of persistent erythrocytes. The leucocytes were subsequently collected from the cellular interface ring, washed in PBS containing 2% FCS and 2 mm EDTA, and then quantified. Viability of purified leucocytes in each sample was determined using trypan blue dye exclusion with the Spark® Cyto (Tecan, Switzerland) system. Cells underwent staining using a series of chicken-specific monoclonal antibodies alongside their corresponding isotype controls, following established protocols43. All antibodies were sourced from Biorad, France. To distinguish monocytes and macrophages in the spleen, we utilised an APC-conjugated antibody (clone Kul01, IgG1) targeting the chicken mannose receptor C-type 1 like B (MRCL1-B) of avian monocytes and macrophages. Additionally, we incorporated staining for MHC-II (clone 21-1A6, IgG1, FITC-conjugated) in this process. T lymphocytes were identified based on the expression of their T-cell receptor (TCR) and coreceptors. For this purpose, we employed monoclonal antibodies against CD4 (clone CT-4, FITC-conjugated), CD8α (clone CT-8, PE-conjugated), or the γδ TCR (clone TCR1, IgG1, APC-conjugated). Chicken thrombocytes were characterised using antibodies recognising the fibrinogen receptor CD41/61 (clone 11C3, IgG1, PE-conjugated). To discern chicken B cells, we utilised an APC-conjugated antibody recognising the antigen Bu-1 (clone AV20, IgG1). Quadrant markers were established based on negative populations and isotype controls. Finally, cell viability was assessed using the fluorescent DNA intercalator 7-aminoactinomycin D (7-AAD, BD Biosciences). Cell acquisition (5 × 105 events) was performed using a BD FACSCanto II cytometer (BD Biosciences), with subsequent analysis conducted using FlowJo 7.5.3 software (TreeStar Inc., Ashland, OR). The percentage of the analysed population relative to total acquired events was employed in graph construction.

Statistical analysis

Unless otherwise specified above or in the figure legends, data are presented as the median or mean ± SEM. The unpaired Student’s t test was employed to compare the means of two independent groups. For three or more independent groups, one-way analysis of variance (ANOVA) was conducted to ascertain statistically significant differences among the means, followed by a Tukey post-hoc test utilising a Studentized range statistic to perform pairwise comparisons between groups. Graph Pad Prism 8.0 software (GraphPad, San Diego, CA, USA) was employed for statistical analysis. Significance levels were set at p < 0.05 for most analyses, unless otherwise indicated.

Results

The GM regulates a wide array of transcriptional signatures governing innate immunity along the gut-lung axis in the chicken

We started by assessing the caecal microbiota profile of our CV chickens during their first weeks of life using 16S DNA sequencing. Diversity patterns across caecal samples collected at days 7, 14 and 21 post-hatch were summarised using four α‐diversity indices. Notably, a significant increase in α‐diversity from day 7 to day 21 was observed (Supplementary Fig. 1a). On day 21, caecal samples exhibited the highest average α‐diversity values (Supplementary Fig. 1a). Further examination of bacterial family relative abundances between 7 and 14 days highlights the dynamics of taxa appearance and disappearance (Supplementary Fig. 1b). Among these families, Lachnospiraceae, whose role in the fermentation of polysaccharides into SCFA is well-documented in mammals44, is present at both 7 and 14 days. Lactobacillaceae, abundant in very young animals, decreases by 14 days. Few operational taxonomic units (OTUs) show notable enrichment between 14 and 21 days, indicating a slower acquisition of new bacterial taxa. Enriched rare OTUs after 14 days belong to Ruminococcaceae, Lachnospiraceae, and Lactobacillaceae families (Supplementary Fig. 1c), suggesting similarities with bacterial taxa found in other studies in the chicken13,14,45,46. β‐diversity measures revealed similar taxonomic profiles within samples collected at each time-point post-hatch, and significant associations were noted between sample type and all four β‐diversity indices (P < 0.001; refer to Supplementary Fig. 1d for the distribution of Bray–Curtis distances). Overall, these data highlight that the isolators employed in this study did not impede the development of a complex caecal microbiota in our chicken line.

Next, we assessed the impact of GM on chicken immune system development using 3-week-old immunocompetent/immunoresponsive26 GF chickens and CV chickens of the same line. RNAseq analysis (padj <0.05, log2-fold change ≥1) revealed significant transcriptional profile changes in GF compared to CV chickens, particularly in the caeca, with ~3000 DEGs showing a downregulation profile, compared to 50 in the spleen and 20 in the lungs. Gene ontology enrichment analysis illustrated disruptions in biological processes, including the cell cycle, biosynthesis, and response pathways to stimuli, emphasising the substantial influence of the GM in the caecal physiology (Fig. 1a). Notable genes that are negatively regulated in GF chickens include JCHAIN, a critical factor in mucosal intestinal protection47, NMRAL1, which regulates innate immunity through NF-κb inhibition and RLR pathways48, and AREG, which encodes amphiregulin, crucial for maintaining the integrity of the intestinal epithelial barrier49. Only a few pathways related to cellular stress or transcriptional regulation showed positive regulation in GF chickens, notably through the sirtuin signalling pathway. Nevertheless, the exact function or expression pattern of these genes or immunological pathways in the chicken are not yet completely resolved, with the exception of JCHAIN or RIGI (which is absent from the chicken genome)50.Fig. 1 The influence of the gut microbiota on gene expression profiles varies according to the tissue studied.

Differently expressed genes (DEG, padj < 0.05, log2 fold change ≥1) from RNAseq data from germ-free (GF) chickens, compared to conventionally raised (CV) chickens at 21 days post-hatch were analysed using Ingenuity Analysis Pathway software, visualised as a “bubble chart.” Circles in each row represent groups of genes associated with a biological function, with size indicating gene count and colour reflecting an intensity gradient ranging from overexpression (orange pattern) to downregulation (blue pattern) in GF chickens using CV chickens as the normalisation controls. The software enables tissue or cell selection and provides inferences on unrepresented genes. A confidence-boosting maximum adjusted p value threshold was applied: 0.45 for caeca (a), including the analyses of the top 244 genes from the caeca RNAseq analyses; 0.4 for lungs (b), including the analyses of the top 280 genes from the lungs RNAseq analyses; and 0.1 for spleen (c), including the analyses of the top 84 genes from the spleen RNAseq analyses. Each group comprised caeca, lungs or spleen samples from 5 GF chickens or 5 CV chickens (n = 5 biological replicates).

In the lung tissue, the variations involved a general downregulation of pathways related to immune responses, development and cell survival in GF chickens (Fig. 1b). Genes affected included those coding for myosin subunits (MYH1C, MYL3, MYL10), and those associated with B lymphocyte development or recruitment, such as CXCL13, the immunoglobulin Lambda gene, and FAB4, implicated in neutrophil recruitment by alveolar macrophages in mouse bacterial pneumonia51. However, chickens have rigid parabronchial lungs, harbouring no alveoli, and birds have no neutrophils, possessing heterophils instead, thus making any parallel to mammalian immunological literature a delicate issue50,52,53. In the spleen, the main secondary lymphoid organ in chickens, we observed a varying degree of downregulation in immune response processes, including cellular and humoral responses, pathogen-dependent signalling pathways, and cytokine pathways in GF compared to CV chickens (Fig. 1c). Particularly, we observed a significant downregulation of the MZB1 gene, encoding a protein highly expressed by marginal zone B lymphocytes in the mouse spleen54. Additionally, the KK34 gene, which encodes a cytokine expressed by T lymphocytes in chickens55, exhibited downregulation. Nevertheless, the particularities of avian lymphoid organs compared to mammalian species50 and how these genes operate within compartments harbouring lymphocytes in the chicken are yet to be described in detail. Furthermore, the antiviral innate response was compromised, with notable downregulation of key genes such as OASL, IFIT5, and IFI6, all involved in the interferon (IFN) stimulation pathway56.

By closely examining a selection of innate immune-related genes regulated by the GM (padj <0.05, log2 fold change ≥1), we observed a certain overrepresentation of IFN-related genes. In the GF caecal tissue, there was significant downregulation of OASL, IRF1, IFI6, IFIT5, and MX1, along with genes related to mucosal immune responses (IL22RA2, SOCS1/3), tissue repair (MMP1), and apoptosis (CASP8) compared to CV chickens (Fig. 2a). The lung tissue shows a similar trend with significant downregulation of IFN-stimulated genes (ISGs) IRF1, IFIT5, IFI6, and OASL and NF-κB binding activity regulators (NFKBIA), and the canonical MHC class I gene BF1 compared to CV chickens (Fig. 2b). There was a modest impact on genes within the splenic compartment, with ISGs like IRF7, IFI6, IFIT5, and OASL significantly downregulated in GF compared to CV chickens (Fig. 2c). Finally, when aggregating DEGs identified through RNAseq analysis, we identified a select group of genes that exhibit consistent downregulation in GF chickens (padj < 0.05, log2 fold change ≥1) (Fig. 2d). Indeed, OASL emerges as a primary standout DEG in all three analysed organs.Fig. 2 The absence of a gut microbiota results in the regulation of several immune-related genes along the gut-lung axis in the chicken.

Differently expressed genes (DEG, padj <0.05, log2 fold change ≥1) associated with immunity and inflammation were identified in the RNAseq analyses of caeca (a), lungs (b), and spleen (c) samples of 21 days-old chickens. This selection included standout genes coding for transcription factors, cytokines, enzymes and other molecules with broad functions in innate and adaptive immunity of vertebrates. Histograms present the relative expression of each gene in log2 format for germ-free (GF) animals compared to conventionally raised (CV) animals. A downregulation trend is particularly pronounced in the caeca and moderately in the spleen. In the lungs, contrasting regulatory phenomena with similar intensities coexist. d Certain DEGs identified using RNAseq analyses followed by an assessment through the Ingenuity Analysis Pathway software (in which a confidence-boosting maximum adjusted p value threshold was applied) are shared by the three organs studied. A Venn diagram illustrates DEG counts per organ, with caeca in blue, lungs in grey, and spleen in orange. Genes shared between organs, where available in the chicken genome, are annotated in the intersections. Immune-related genes, indicated in bold, reveal OASL as the sole immune-related gene significantly regulated in all three examined organs. Each group comprised caeca, lungs or spleen samples from 5 GF chickens or 5 CV chickens (n = 5 biological replicates). Data are represented as the mean ± SEM.

Attempts were made to complement these data with information on the dynamics of immune-cells populations likely to be regulated by the GM in our model. Specifically, we have tried to perform deconvolution analyses using the R package “granulator”57 to access the heterogeneity in the cellular composition of our bulk RNAseq from caeca and lungs, a technique that has been widely used to track compositional alterations of cell types in gene expression data. However, the lack of adequate calibrators for the chicken available in public databases led to our analyses yielding no significant results. We therefore opted to perform a flow cytometry analysis of spleen samples from CV and GF chickens at 7, 14 and 21 days post-hatching. As the main secondary lymphoid organ in chickens, the spleen’s leucocyte provides insightful information regarding the GM’s impact on systemic immunity. Total leucocyte number showed no significant changes between 7 and 21 days-post hatching between CV and GF chickens (Supplementary Fig. 2a). We then applied a staining strategy combining antibodies targeting lymphocytes (CD4, CD8α, γδTCR) and antigen-presenting avian phagocytes (MRCL1-B and MHC-II) (Supplementary Fig. 2b). T CD4 lymphocytes or gamma-delta T cells (CD4+γδTCR+) showed no differences between CV and GF animals independent of the time-point analysed (Supplementary Fig. 2c). The same was observed for CD8α+ cells, which may not only include cytotoxic T cells but also NK cells and dendritic cells43. However, MHC-II negative CD8α+ cells, possibly including CTLs and NK cells, presented a significant higher percentage in GF animals at 21 days post-hatching (Supplementary Fig. 2c). Finally, no changes were observed in antigen-presenting avian phagocytes (MRCL1-B+MHC-II+), regardless if they presented high or low levels of expression for MHC-II (MHC-IIhi or MHC-IIlo) (Supplementary Fig. 2d).

Overall, these findings highlight the GM’s significant impact on the transcriptional regulation of IFN-related genes along the gut-lung axis in the chicken, while cellular heterogeneity, at least in the spleen, the main systemic lymphoid organ, remains largely unaffected by the lack of a GM.

SCFA distribution in the chicken is significantly reduced in peripheral tissues in the absence of a GM

We have previously established that SCFAs are detectable beyond the gastrointestinal tract in CV chickens and that their concentrations are noticeably reduced or absent in GF animals16. Here, our metabolomic analysis corroborates this observation in both caecal contents and lungs (Fig. 3a–d) and extends it to the spleen and serum (Supplementary Fig. 3a–d). More specifically, the concentrations of SCFAs in the caecal contents ranged from 1000–2000 nM/mg for acetate, 100–300 nM/mg for butyrate, and 100–400 nM for propionate. In the lungs, the concentrations found ranged from 0.2–04 nM/mg for acetate, 0.02–0.06 nM/mg for butyrate, and 0.025–0.075 nM/mg for propionate. As expected, the absence of SCFAs-producing microbes in the gut leads to a systemic loss of butyrate, propionate, and acetate across the gut-lung axis, with more pronounced differences in caecal contents (Fig. 3a), although acetate is still present to some extent, probably from residual food sources.Fig. 3 Germ-free chickens show an altered metabolic landscape along the gut-lung axis.

We used H1-NMR to analyse metabolite concentrations (nM/mg) in caecal contents (a) and lungs (c) of 21-day-old conventional (CV) and germ-free (GF) chickens. Metabolic pathway enrichments for caecal contents (b) and lungs (d) were determined using MetaboAnalyst103. This tool identifies essential metabolic pathways (e.g., amino-acid metabolism, biosynthesis, catabolism) across organs. Ratios for each pathway were calculated as “number of metabolites in pathway X”/“Total number of metabolites in pathway X in the database,” and enrichment testing used the “globaltest” method104, with resulting p values generating colour scales for the histograms. No prior transformations (normalisation, exclusion threshold) were applied to the presented data. Each group comprised caecal contents or lung samples from 6 GF chickens or 6 CV chickens (n = 6 biological replicates). Unpaired student’s t test was employed for statistical analyses (a, c). Data are represented as the mean ± SEM.

These differences also extend to other metabolites, including amino acids and Krebs cycle molecules (e.g., alanine, aspartate, leucine), abundant in CV animals’ caecal contents, spleen, and lungs but diminished in GF animals (Fig. 3a–d and Supplementary Fig. 3a–d). In caecal contents of GF chickens, pathways related to SCFA metabolism and amino acid biosynthesis are prominently downregulated (Fig. 3b). Serum shows similar profiles, with identified pathways such as vitamin K metabolism and ethanol degradation being altered (Supplementary Fig. 3d). Finally, lung pathways primarily involving amino acid synthesis and essential compounds, such as folate and glutathione, are significantly altered in GF chickens (Fig. 3d), while spleen pathways exclusively encompass alterations in the biosynthesis or metabolism of sugars, nitrogenous bases, and membrane components in the absence of a GM from birth (Supplementary Fig. 3b). The complete dataset of all quantified metabolites is provided in Supplementary Figs. 4 and  5. These observed disparities underscore the GM’s profound influence on chicken metabolism and immune response, providing solid evidence to further explore these regulatory mechanisms in the context of infectious challenges.

The transcriptional immune response to avian influenza infection in the lung is profoundly regulated by the GM

To define the putative role of the GM in mucosal immunity, we used a H7N1 LPAIV infection model in immunocompetent/immunoresponsive26 CV and GF chickens24,25. Intratracheal administration of 5 × 105 EID50 resulted in nuanced disease with interstitial pneumonia and an innate antiviral response in the lungs. Over a 3-day post-infection period, both CV and GF chickens were monitored closely. Only 3 out of 24 infected animals per group died during the 3-day infection period. There were thus no significant differences between CV or GF animals in terms of mortality. Macroscopic lung lesions were readily observed at 2–3 days post-infection, showing necrosis, haemorrhage, and inflammatory tissue, with low to no mortality observed as mentioned above (Fig. 4a). Despite individual variations, viral replication occurred in the lungs and, to a minor extent, in the caeca (Fig. 4b, c). Surprisingly, we did not observe significant differences in lung pathology or viral burden between CV and GF chickens. Nevertheless, at day 2 post-infection, four out of seven GF birds had undetectable viral RNA loads in the lungs, whereas in the CV group, only one bird showed an absence of detectable viral replication. This suggests that the majority of GF birds may better deal with viral clearance at this specific time-point, although statistical analyses did not reveal significant differences between the two groups (Fig. 4b).Fig. 4 The transcriptional immune response to avian influenza infection in the chicken lung is profoundly regulated by the gut microbiota.

Germ-free (GF) animals infected via the tracheal route with a 5 × 105 EID50 dose of LPAIV H7N1 exhibited a disease profile indistinguishable from their conventional (CV) counterparts. Macroscopic lung lesion scores (haemorrhage, necrosis, oedema) (a) and viral genome copies in the lungs (b) or in the caeca (c) at Days 1–3 post-infection (p.i.) show no significant differences between the two groups. d Examination of infection kinetics in GF chickens’ lungs using Medium-Throughput qPCR (Fluidigm) indicates sustained expression of genes associated with inflammation (IL8L1, IL1B) and innate antiviral immunity (IFNA, IFNB, IRF7, and OASL – upper-right box) at Day 3 p.i. The heatmap illustrates relative expression levels in GF-infected animals compared to CV-infected animals on the same day. Each group comprised caeca or lungs samples from 5 mock and 8 infected GF chickens or 5 mock and 8 infected CV chickens (n = 5–8 biological replicates), depending on the time-point. a, c, d, each biological replicate is the mean of three technical replicates. Statistical analysis employed unpaired student’s t test (a–c) or one-way ANOVA followed by Tukey multiple comparison test (d), where *p < 0.05 and ***p < 0.001. Data are represented as the median (a–c) or the mean ± SEM (d, OASL expression).

By analysing pulmonary gene expression dynamics during H7N1 infection, we tracked 48 immune-related genes58 over a three-day post-infection (p.i.) period in CV and GF animals (Fig. 4d). Disparities emerged in the lungs of GF animals, especially in key antiviral innate immunity genes such as IRF7, IFNA, IFNB, IFNL, and STAT1, with OASL particularly showing a very distinctive regulatory pattern in the absence of a GM (Fig. 4d, upper-right box). These genes showed lower expression at 24 and 48 hours compared to CV animals but a significant elevation at 72 hours, indicating a sustained antiviral response beyond 48 hours in GF animals. These changes also influenced mucosal immunity, as seen in the distinct regulation pattern for IL22 and its receptors (IL22RA1, IL22RA2) in GF animals. Furthermore, reduced expression of anti-inflammatory genes like TGFB and SOCS3 at 72 hours p.i. in GF animals possibly indicates an unresolved inflammatory state. In summary, while the absence of GM does not hinder antiviral responses against AIV in chickens, it influences the quality and amplitude of the immune response at the transcriptional level.

Butyrate regulates cellular and molecular processes in chicken respiratory epithelial cells

Shifting our focus to investigations in our unique chicken respiratory epithelial cell line CLEC21327, we aimed to model events in the communication between SCFAs and the respiratory tissue. Cytotoxicity tests on this cell line revealed a wide tolerance to SCFAs (Fig. 5a). Despite concentrations exceeding physiological lung levels in chickens (3–10,000 µM), this concentration range facilitated the examination of SCFA-induced effects on metabolism and transcription. Acetate and propionate alone induced moderate transcriptional responses in CLEC213 (Fig. 5b, d). For instance, the expression of IFNB, OASL, STAT1 (with propionate), and IL1B (with acetate) increased by approximately three-fold compared to untreated cells. In contrast, butyrate showed a mean 5-fold increase for the analysed genes (Fig. 5c), with a strong concentration-dependent induction of OASL (Fig. 5e).Fig. 5 Short-chain fatty acids elicit cellular and molecular responses in chicken lung epithelial cells.

The CLEC213 chicken lung epithelial cell line was treated with different concentrations of SCFA for 16 h. Cytotoxicity of butyrate, propionate, and acetate was assessed using a Spark® Cyto and expressed as the percentage of cell viability compared to untreated controls (a). qPCR analysis demonstrates that acetate (b), butyrate (c), and propionate (d) elicit different responses in terms of selected innate immune gene expression in CLEC213 cells. Notably, the induction of OASL by butyrate exhibits a concentration-dependent relationship (e). Statistical analysis was performed using One-way ANOVA followed by Tukey multiple comparison test, with significance levels indicated as follows: *p < 0.05, **p < 0.01, ***p < 0.005, ****p < 0.001. a, b, an n = 3 biological replicates is shown, in which each biological replicate is the mean of three technical replicates. c–e an n = 5–7 biological replicates is shown, in which each biological replicate is the mean of three technical replicates. Data are represented as the mean ± SEM.

Next, we further defined the effects of butyrate by analysing the metabolic and transcriptional profiles of CLEC213 cells treated with 3 mM of butyrate for 16 hours. While concentrations of detectable metabolites increased in untreated cells, only butyrate concentrations showed statistically significant differences (Fig. 6a). Enriched metabolic pathways revealed that butyrate functions in CLEC213 cells by enhancing the biosynthesis of essential compounds (e.g., vitamin K, pyruvate, folate, etc.) and amino acids, particularly tryptophan. This suggests a broad impact of butyrate on chicken lung epithelial cells, enhancing their overall metabolic activity (Fig. 6b), consistent with the well-documented effects of butyrate on cellular metabolism59.Fig. 6 Butyrate regulates metabolic and transcriptional signatures in chicken lung epithelial cells.

CLEC213 chicken lung epithelial cells were treated with 3 mM butyrate for 16 h, resulting in an increase in metabolites associated with various biosynthetic pathways (tryptophan, glycogenesis, pyruvate), as analysed by H1-NMR (a). Statistical significance was assessed using unpaired student’s t test. b Metabolic pathway enrichment analysis, employing the “globaltest” method104, revealed pathway ratios, with colour scales in histograms corresponding to p values. c Butyrate induced positive regulation across diverse biological functions in chicken lung epithelial cells, as depicted in a ‘bubble chart’ generated from RNAseq data analysis using Ingenuity Analysis Pathway software. Circles in each row represent groups of genes associated with a biological function, with size indicating gene count and colour reflecting an intensity gradient ranging from overexpression (orange pattern) to downregulation (blue pattern) in butyrate (3 mM) treated cells using untreated cells as the normalisation controls. A confidence-boosting maximum adjusted p value threshold of 045 was applied. d The log2 relative expression of a group of selected immune-related genes extracted from all DEG from RNAseq data (padj < 0.05, log2 fold change ≥1) in cells receiving butyrate (3 mM) compared to untreated cells at 16 h. a, b an n = 3 biological replicates is shown, while in c, d an n = 6 biological replicates. Data in a are represented as the median. Data in d are represented as the mean ± SEM.

Ontology enrichment and gene expression analyses revealed that butyrate significantly influences cellular functions in CLEC213 cells, encompassing cytokine signalling, pathogen-influenced pathways, and cellular immunity (Fig. 6c). In addition, a cluster of negatively regulated genes associated with cell growth aligns with butyrate’s role as a histone deacetylase inhibitor (HDACi)31,60, inducing cell cycle slowdown or arrest (Fig. 6c). Immune-related gene regulation by butyrate shows a bidirectional pattern, balancing anti-inflammatory and pro-inflammatory signatures. This includes the downregulation of NFKB1 and CD40, alongside overexpressed genes like SOCS3 and SOCS4, known for negative regulation of cytokine signalling (Fig. 6d). Butyrate also exhibits pro-inflammatory actions, as seen in the upregulation of CCL4, IL6R, and ISGs (OASL, IFIT5, IRF1, EIF2AK2), mirroring observations in GF animals (which lack systemic butyrate) and emphasising the role of SCFA in regulating ISGs in the respiratory mucosa.

To verify butyrate’s HDACi functions, CLEC213 cells were treated with 0.3 or 3 mM for 16–48 hours. Initial cell proliferation remained unaffected, but between 16 and 24 hours, treatment with 3 mM butyrate inhibited proliferation and led to 15% reduction in cell viability. A similar growth inhibitory profile was also found for 300 µM of butyrate at later timepoints, whilst cell viability remained unchanged at this concentration (Fig. 7a, b). Overall, in the presence of butyrate, the cells reached a steady state where neither the total cell count nor the number of dead cells increased, reflecting a net cessation of cell proliferation (Fig. 7b). Butyrate was shown to inhibit mammalian HDACs from classes I and II31. Using the HDAC-Glo I/II assay kit, we measured enzymatic activity in CLEC213 cells. Luminescence counts (Fig. 7c) and percentage of inhibition (Fig. 7d) indicate a significant reduction in HDAC activity following treatment with 3 and 10 mM of butyrate, emphasising its impact on histone deacetylase activity in chickens. Trichostatin A, a known broad-range HDAC inhibitor, exhibits a more pronounced decrease, highlighting butyrate’s HDACi function in CLEC213 cells.Fig. 7 Butyrate possess HDACi functions and strongly regulates Sp1 in chicken respiratory epithelial cells.

Incubation of CLEC213 cells with 0.3 and 3 mM of butyrate for 16 hours resulted in reduced cell counts (per ml) (a) with no major impact on cell viability (b), as assessed using a Spark® Cyto and expressed as the percentage of cell viability compared to untreated controls. c The exposure of CLEC213 cells to 3 mM of butyrate for 16 hours led to diminished histone deacetylase (HDAC) activity as assessed using the HDAC-Glo I/II assay. Data are presented as luminescence ratio (RLU) (c) or as the percentage relative to the control group (d), with decreased luminescence indicating HDAC inhibition. Trichostatin A (50–200 nM) was used as a positive control. e Gene ontology enrichment analysis in butyrate-treated cells (3 mM for 16 h) reveals a significant enrichment (green box) of two Sp1 transcription factor motifs among the top 1000 differentially expressed genes obtained through RNAseq (padj < 0.05, log2 fold change ≥1). Enrichment analysis originally used g:profiler105 for molecular function, biological processes, cellular components, various databases, and transcription factors. f Identification of potential binding sites for Sp1 in the OASL promoter of selected bird and mammalian species. One-way ANOVA followed by Tukey multiple comparison tests (a–d) was employed for statistical analyses, where *p < 0.05, **p < 0.01. a, b an n = 3 biological replicates are shown. c, d an n = 3 biological replicates is shown, in which each biological replicate is the mean of three technical replicates. Data in e are representative of an n = 6 biological replicates. Data are represented as the mean ± SEM.

The transcription factor Sp1 is involved in the regulation of OASL by butyrate in chicken respiratory epithelial cells

Gene ontology analysis from RNAseq data of cLEC213 cells treated with 3 mM butyrate for 16 hours significantly revealed two Sp1 transcription factor motifs (Fig. 7e). In mammals, inhibition of Sp1 with mithramycin A halts TGF-β induction by butyrate in intestinal epithelial cells61. We then asked whether this observation could apply to our system. Exploring Sp1-binding sites in the chicken OASL promoter, we observed a clear pattern where domesticated bird species have 1–3 Sp1-binding sites, while humans and pigs show none to 1 (Fig. 7f and Supplementary Fig. 6). To test if inhibiting Sp1 mitigates butyrate-induced OASL overexpression, we used siRNAs targeting distinct chicken Sp1 transcript regions. siRNA n°2 achieves significant SP1 expression knockdown (Supplementary Fig. 7a), even when cells are treated with butyrate before transfection. In addition, RNAseq analysis of siRNA-transfected CLEC213 cells versus controls highlights SP1 as the prominently downregulated transcription factor among Sp family members (Supplementary Fig. 7b). SP1 knockdown significantly downregulated canonical Sp1-regulated genes LMCD1 and AQP162,63, without affecting vital cellular functions in the CLEC213 cell line. To our knowledge, the function of these genes in the chicken are still not resolved, although their conservation among several vertebrate species suggests that they may exhibit overall conserved biological functions.

Transfecting CLEC213 cells with a Sp1-targeting siRNA and treating them with butyrate for 16 hours resulted in a 50% reduction in OASL expression, indicating a partial dependence on Sp1 binding to the OASL promoter (Supplementary Fig. 7d). This effect did not impact the type I IFN stimulation pathway, as demonstrated by the similar response to IFN-α in transfected and stimulated cells compared to non-transfected cells (Supplementary Fig. 7e). Additionally, mithramycin A, a competitive inhibitor of Sp1 and Sp3, showed cytotoxicity only beyond 0.1 mM, and while it exhibited a non-statistically significant reduction in OASL induction by butyrate, it also affected OASL induction by IFN-α, suggesting an unexpected interaction between these supposedly independent pathways (Supplementary Fig. 7c).

Antiviral effects of butyrate in chicken respiratory epithelial cells are partially dependent on the induction of OASL expression

To explore butyrate’s potential antiviral effects in CLEC213 cells, we applied our well-established LPAIV H1N1 (A/Mallard/Marquenterre/Z237/83) infection model64, which yields productive infection in CLEC213 cells with mild to moderate cytopathic effects and a low-grade but consistent IFN-I response. Cells pre-treated with 3 mM butyrate for 16 hours, and infected at a MOI of 0.01 or 0.1, exhibited significantly reduced infectious particles, especially at an MOI of 0.1 for 6 h, despite lower infection titre (104−105 PFU/ml) (Fig. 8a). H1N1 genome copies were also significantly decreased (Fig. 8b), while transcription of OASL and, to a lesser extent, IFNB were significantly enhanced by infection, particularly at an MOI of 0.1 (Fig. 8c). Overall, butyrate pre-treatment led to a reduction in viral load by 50% and amplified OASL expression induced by H1N1 infection (MOI 0.1) in chicken respiratory epithelial cells. Interestingly, butyrate was not effective in enhancing IFNB expression induced by H1N1 infection (Fig. 8c), which was already low without treatment. We further confirmed the reproducibility of this pattern by repeating the antiviral activity assay with the LPAIV H7N1 (A/Turkey/Italy/977/1999) strain used in the in vivo challenge infection experiment. Infections at an MOI of 0.01 or 0.1 yielded markedly higher viral RNA loads when compared to LPAIV H1N1 infection (Supplementary Fig. 8a). Likewise, butyrate pre-treatment (3 mM) led to a significant reduction in H7N1 genome copies at an MOI of 0.1. Interestingly, H7N1 infection led to enhanced OASL and IFNB expression levels at both MOIs similar to mock control groups pre-treated with butyrate (Supplementary Fig. 8b, c). Nevertheless, butyrate pre-treatment significantly enhanced both OASL and IFNB expression in the infected cells to levels higher than those seen in the H1N1 infection experiment at the same conditions (Supplementary Fig. 8b, c), thus confirming the overall capacity of butyrate to limit AIV replication and enhance IFN-I/ISG responses in chicken respiratory epithelial cells.Fig. 8 Butyrate enhances viral clearance in chicken lung epithelial cells upon infection with a low-pathogenicity H1N1 avian influenza strain via increased type I IFN/ISG responses.

Butyrate treatment (3 mM) improves viral clearance in CLEC213 cells infected with an LPAIV H1N1 avian influenza strain. Incubation with butyrate prior to infection (MOI 0.1 for 6 hours or MOI 0.01 for 16 hours) significantly reduces viral titre (PFU/ml) (a) and genomic copies (b) compared to untreated cells. Additionally, qPCR analysis showed increased expression of IFNB and OASL in butyrate-treated infected cells (c). Silencing OASL expression abolishes butyrate’s antiviral effects, resulting in infectious titre identical to untreated cells, concomitantly with reduced OASL expression (d). Unpaired student’s t test (a, b) and One-way ANOVA followed by Tukey multiple comparison test (c, d) were employed for statistical analyses, where p < 0.05, **p < 0.01, ***p < 0.001. a an n = 3 biological replicates are shown. b–d an n = 3 biological replicates is shown, in which each biological replicate is the mean of three technical replicates. a, b are represented as the median. Data in c, d are represented as the mean ± SEM.

RNAseq analysis from LPAIV H1N1-infected CLEC213 cells further revealed butyrate’s modulation of infection-related pathways, notably by suppressing HIFs, TGF-β, p53, and NF-κB signalling pathways, emphasising its role as a negative regulator of inflammation during H1N1 infection in chicken lung epithelial cells (Supplementary Fig. 9a). Finally, we developed three siRNAs targeting chicken OASL transcripts and used them in CLEC213 cells to target OASL expression induced by butyrate (Supplementary Fig. 9b). Subsequent H1N1 infection highlighted OASL’s pivotal role in butyrate-induced antiviral effects. qPCR was used to confirm diminished OASL expression following siRNA transduction in butyrate-treated cells. In these cells, infectious titres of H1N1 virus were comparable to untreated cells (Supplementary Fig. 9c). Thus, the overall antiviral effects of butyrate would involve enhanced metabolic activity, regulatory influence on inflammatory signals, and pronounced overexpression of ISGs, notably OASL.

Discussion

The establishment of the GM is a critical early event in vertebrate life. Specific bacterial genera and species influence various aspects of poultry physiology2,12,15,65. Using a specific pathogen-free (SPF) laying hen model, we aimed to ensure the development of a biologically significant GM before exploring the gut-lung axis mechanisms. Between days 7 and 21, the caecal microbiota saw an enrichment of numerous new species, predominantly within the Firmicutes phylum, including bacteria from the Lachnospiraceae, Lactobacilleae, and Ruminococcacae families. These findings align with Glendinning’s study on broiler chickens (Ross 308) raised without maternal influence from birth and followed until their 35th day of life66. Notably, the Lachnospiraceae genus Roseburia, known for butyrate production, is more prevalent in our animals even at a young age, ensuring significant SCFA production, particularly butyrate. This also emphasises the crucial role of Roseburia in poultry mucosal homoeostasis, mainly attributed to the anti-inflammatory properties of butyrate67–69.

Comparing metabolites and transcriptional profiles in CV and GF animals highlighted the GM’s profound influence on mucosal homoeostasis in birds. Most of the data available in the literature to discuss these aspects was performed in mammalian model species, especially mice70–72. Virtually no mechanistic work has been performed in detail in poultry birds, especially without resorting to antibiotic-induced microbiota suppression73. In GF animals, absence of GM leads to altered gene expression and metabolite concentrations in the caecal compartment, signalling changes in intestinal physiology, probably represented by reduced villi, thinner walls, and a more acidic pH as already discussed elsewhere74,75. Notably, the near total absence of butyrate and propionate in the caecal contents and in peripheral organs of GF animals confirms their exclusive bacterial origin in our model, since they are unlikely to be produced in significant amounts by the host metabolic machinery59,76. Furthermore, the acetate concentrations we observed in GF birds likely originates from dietary sources77. Finally, while we consistently observed SCFAs at varying gradients and compartments along the gut-lung axis in chickens as it was observed for mammals7,78, we could not elucidate the cellular and molecular mechanisms underlying these systemic distribution. These observations underscore the importance of examining the roles of SCFA in various cell or tissue types along this biological gradient in birds, notably in regard to their epithelial utilisation, integrity, splanchnic utilisation, and peripheral availability, mirroring work performed in mammalian species78.

Additionally, the absence of GM results in decreased amino acid concentrations in the periphery, which could directly influence systemic immunity. For instance, the 85% drop in alanine concentration in the caecal contents of GF animals would be crucial for T-lymphocyte activation as suggested elsewhere79. Amino acids play pivotal roles in various biological processes essential for immune cell development and function79,80. Despite these metabolic disparities, GF animals exhibit no zootechnical differences up to 4 weeks of age16,18, suggesting potential compensatory mechanisms for lower nutrient absorption.

The level of knowledge for the molecular mechanisms underlying the communication between the GM and innate immunity in birds does not yet match the level achieved for humans and mice7,81–83. Our investigation revealed significant alterations in gene expression signatures due to the absence of a GM in the chicken. In caeca, lungs, and spleen of GF chickens, we observed a substantial reduction in the expression of ISGs, such as IRF1, IFIT, OASL, and MX1, aligning with findings in GF mice84. Particularly, OASL, a pivotal player in chicken antiviral immunity, for example during Newcastle disease virus (NDV) infection34, exhibits significant downregulation along the gut-lung axis in our model. In-depth analysis of the chicken OASL previously revealed a unique motif, LRLRGG, absent in mammals, which might represent a compensatory mechanism for the lack of ISG15 in birds85,86. Analysis of OASL sequences from six diverse bird species indicate that these features are likely conserved among avian OASLs. Furthermore, the covalent bonding of avian OASLs to themselves and other immunological proteins, along with structural similarities between avian OASLs’ tandem ubiquitin-like (Ubl) domain and ISG15, suggests a distinctive approach to immune signalling pathways in avian species85,86. ISG15 is an important regulator of the mammalian antiviral innate immune response by regulating RIG-I, NF-κB, cytokine and chemokine production, and immune cell activation87. However, studies linking the GM and SCFA to the regulation of ISG15 or OASL are currently lacking. Overall, OASL’s regulation in GF animals appears markedly altered both at homoeostasis and infection, emphasising its dependence upon positive microbial signals coming from the GM in the chicken. Our comparative study using GF and CV chickens thus provides robust evidence of OASL regulation by the microbiota, shedding new light on evolutionary biology and comparative immunology.

Over the past decade, the link between the microbiota and type I IFNs has been extensively documented in humans and mice88–90. Recent research has revealed that commensal microbes play a crucial role in maintaining tonic type I IFN levels, which are essential for mounting an effective antiviral immune response. Although specific bacterial strains associated with IFN-I signalling vary among studies, it is noteworthy that the effect of the microbiota is not reliant on a single species. Rather, a diverse community of microbes can reproduce type I IFN responses84. These findings collectively support the hypothesis that microbiota-induced type I IFN signalling contributes to the maintenance of mucosal homoeostasis and immune tolerance. Likewise, our findings suggest that at homoeostasis, the innate immune system, including a consistently large panel of IFN-I related genes (e.g., OASL, IFI6, IFIT5), in GF chickens is systematically tuned down due to the physical absence of bacteria in the intestinal tract and the potential lack of microbiota-derived metabolites in the bloodstream. This leads to compromised communication with distant organs such as the lungs, and possibly affects the development of tertiary lymphoid structures in mucosal tissues15 and the adequate coordination of antiviral responses.

GF animals infected with an H7N1 AIV strain exhibited clinical manifestations similar to CV animals. Moreover, the virus replicated in the caeca, as described for AIV in domesticated birds91,92. Surprisingly, no differences were observed between GF and CV animals with respect to viral replication and clinical outcomes, particularly in regard to lung pathology. This contradicts mouse literature suggesting the GM’s protective role against respiratory viruses88,93. Disparities may result from a less effective innate antiviral immune response in the context of GM imbalance, especially in IFN-related genes (IFNB, IRF7, OASL, and MX1). At 24 and 48 hours p.i., GF animals showed downregulation of most tested genes in the lungs, possibly due to a lack of GM-mediated signals that otherwise prime innate immunity development and triggering threshold. At 72 hours p.i., GF chickens exhibited significant overexpression of ISGs (OASL or IRF7), potentially aiding in viral clearance. Interestingly, at day 2 post-infection, the lungs of 4 out of 7 GF birds showed nearly undetectable viral loads, while this observation was only made for a single bird in the CV group. This points to a superior viral clearance in the majority of GF birds at this time-point, although statistical analyses did not reveal significant differences between the groups (Fig. 4b). Subsequent analysis of Fluidigm data aimed to explore within-group differences in the relationship between innate immune regulation and viral load, considering the possibility that any lack of significance might stem from a division between responders and non-responders. However, correlations between the two analyses were not identified, possibly due to the limited number of animals used and significant intra-group variability. Despite the observed altered dynamics, both in viral load and gene expression, viral load in GF animals at day 3 post-infection is nearly identical to CV animals. This contrasts with previous findings in antibiotic-treated chickens challenged with a H9N2 LPAIV, where decreased IFNA expression and increased viral particle excretion were observed94. Regarding OASL, its expression decreased by 2 log2 at homoeostasis and surged to a level 60 times higher 24 hours post-infection. Among tested genes, only IRF7 followed a similar kinetics. Therefore, it appears that a delay in triggering the antiviral innate response in the lungs of GF chickens does not impair the proper control of viral load compared to CV chickens. This may be related to the fact that, at homoeostasis, the selected innate immune genes in GF chickens are already at a lower basal level than those in CV counterparts. This could potentially attenuate the exacerbated inflammatory response that leads to severe pneumonia and increased mortality. However, this is not what we observed, as overall lung pathology and clinical scores are similar between CV and GF chickens. Further studies with a broader analysis beyond day 3 are necessary to determine if GF animals can sustain viral load control and if their innate immune response will peak and resolve similarly to CV animals. This remains undefined at present and could lead to the discovery of new molecular mechanisms through which the chicken GM directs innate immunity toward effective and pro-resolutive viral clearance.

A previous report indicated that butyrate influences the response to type I IFN in A549 adenocarcinoma-derived human alveolar basal epithelial cells by modulating the induction of specific ISGs, both positively and negatively, depending on the gene95. Unlike the previously documented mechanism of butyrate’s regulation of the type II IFN response96, it does not seem to inhibit type I IFN-mediated STAT activation or nuclear translocation. Furthermore, baseline expression of over 30 ISGs was significantly upregulated by at least 4-fold in the presence of butyrate alone, including antiviral restriction factors OASL and BST2. Consistent with our findings, other researchers also observed that butyrate alters the induction level of several ISGs95, revealing a novel mechanism through which butyrate may influence mucosal immunity in the lungs and respiratory virus infections in both mammals and birds. Finally, these results suggest that different subsets of ISGs would have varying requirements for HDAC activity, unveiling a previously unappreciated layer of complexity in the regulation of these genes. The authors, therefore, do not dismiss the possibility that butyrate’s effect on ISGs stems from increased acetylation of non-histone proteins, which are also substrates of HDACs97.

Indeed, butyrate’s HDACi properties appears to contribute to OASL expression in our CLEC213 cells, albeit showing modest effects in our chicken respiratory epithelial cell model compared to a bona fide HDAC synthetic inhibitor. Unique chicken-specific butyrate functions and the absence of well-characterised GPCR receptors for SCFAs in chickens may contribute to this phenomenon. It was previously reported that more than 20 genes encoding FFAR2/GPR43 paralogs exist in the chicken genome. Although all paralogs seem to encode full-length proteins, the functionality of those genes remains unknown. Moreover, this expansion among FFAR2/GPR43 genes seems to be chicken-specific, as it was not found in other galliform birds such as turkeys and quails98. Inspired by studies on butyrate and TGF-β production in mouse intestinal cells61, we found Sp1 transcription factor involvement in its mode of action in chicken respiratory epithelial cells. Comparing OASL promoter sequences in birds and mammals revealed abundant Sp1 binding sites in birds, a contrast to the scarcity in mammals, emphasising avian immune system peculiarities. siRNA-mediated Sp1 knockdown in CLEC213 cells partially reduced OASL expression by 50% upon butyrate treatment, confirming Sp1’s involvement and suggesting a hitherto unidentified mechanism. Sp1 tightly regulates numerous mammalian genes, including both coding and noncoding RNAs, with its expression/activity controlled by various post-translational modifications, such as phosphorylation, acetylation, glycosylation, ubiquitination, and sumoylation99. Genetic disruption of Sp1 in mice results in embryonic lethality, highlighting its crucial role in developmental processes100.

Finally, butyrate exhibits minor yet noteworthy antiviral effects on chicken respiratory epithelial cells. Butyrate has been shown to directly influence virus replication in certain cases, either independently or in addition to its impact on the type I IFN response, as observed in the case of HIV-1101. Moreover, GM-derived butyrate has been recognised for its protective role against influenza infection, enhancing the activity of macrophages and virus-specific CD8+ lymphocytes81. Our study corroborates the antiviral activity of butyrate by demonstrating its modest, yet reproducible effect on H1N1 and H7N1 infection in CLEC213 cells, notably through a boost in OASL expression. Furthermore, transcriptome sequencing shows butyrate induces wide cellular changes, which possibly contribute to protection against influenza by modifying active infection pathways. Notably, butyrate inhibits the HIF pathway, previously linked to inflammation during H5N1 infection in macaques102, potentially mitigating inflammation by preventing HIF activation in pulmonary epithelium. Beyond Sp1, butyrate may regulate OASL expression through additional metabolic pathways yet to be defined, highlighting its dual role in shaping immunity and metabolism in our respiratory epithelial system. This unveils a novel mechanism through which a GM-derived SCFA modulates antiviral immunity in peripheral organs of birds. Further investigation into the multifaceted regulation of the type I IFN response by butyrate and its varying effects on virus infection in vertebrates is warranted.

Conclusions

Our research revealed the intricate interaction between caecal microbiota-derived SCFAs, particularly butyrate, and their role in regulating protective immunity and metabolism across the gut-lung axis in chickens. These insights have broad implications for understanding the gut-lung axis and developing prophylactic strategies in commercial poultry through microbiota modulation. Further exploration of butyrate’s multifaceted effects and its interactions with specific cellular pathways is crucial to harnessing the full potential of SCFAs as next-generation postbiotics for poultry birds.

Supplementary information

Supplementary Material

Description of additional supplementary files

Supplementary Data 1

Supplementary Data 2

Supplementary Data 3

Supplementary information

The online version contains supplementary material available at 10.1038/s42003-024-06815-0.

Acknowledgements

We express our gratitude to the staff of the Plate-forme d’Infectiologie Expérimentale (UE PFIE, INRAE) for providing the necessary infrastructure and technical support essential for conducting the animal experiments in this study. This work was funded by the INRAE, France (specifically Santé Animale—SA and Microbiologie et Chaine Alimentaire—MICA divisions), by “EuroFéRi” (FEDER-FSE Centre Val de Loire 2014-2020, no. EX 010233) and by VetBioNet, EU H2020 (Grant Agreement no. 731014).

Author contributions

V.S.M. acquired, analysed, and interpreted data. V.G., M.C., I.F., and E.K. contributed to data acquisition and analysis. F.R. and I.C. provided molecular and bioinformatics tools and participated in data analysis. E.H. and E.C. conducted microbiological analyses. B.F., P.V., and F.K. analysed and interpreted data using bioinformatics. S.T. provided tools, infrastructure, and funding. R.G. conceived the experiments, acquired, analysed, and interpreted data, and provided funding. All authors contributed to writing the article and approved the submitted version.

Peer review

Peer review information

Communications Biology thanks Mohamed Faizal Abdul-Careem, Dirkjan Schokker, and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Primary Handling Editors: Sabina Leanti La Rosa and Tobias Goris.

Data availability

The authors declare that all data supporting the findings in this study are available within the article, within Supplementary Data files, or from the corresponding author on reasonable request. All RNAseq data and 16S sequencing data were deposited in the ELIXIR Deposition Database ArrayExpress under the accession numbers E-MTAB-14262, E-MTAB-14264, E-MTAB-14264, E-MTAB-14270, and E-MTAB-14346.

Competing interests

The authors declare no competing interests.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Thaiss CA Zmora N Levy M Elinav E The microbiome and innate immunity Nature 2016 535 65 74 10.1038/nature18847 27383981
Thaiss, C. A., Zmora, N., Levy, M. & Elinav, E. The microbiome and innate immunity. Nature 535, 65–74 (2016).27383981
2. Rowland I Gut microbiota functions: metabolism of nutrients and other food components Eur. J. Nutr. 2018 57 1 24 10.1007/s00394-017-1445-8 28393285
Rowland, I. et al. Gut microbiota functions: metabolism of nutrients and other food components. Eur. J. Nutr. 57, 1–24 (2018).28393285
3. Willing BP Russell SL Finlay BB Shifting the balance: antibiotic effects on host–microbiota mutualism Nat. Rev. Microbiol. 2011 9 233 243 10.1038/nrmicro2536 21358670
Willing, B. P., Russell, S. L. & Finlay, B. B. Shifting the balance: antibiotic effects on host–microbiota mutualism. Nat. Rev. Microbiol. 9, 233–243 (2011).21358670
4. Koh A De Vadder F Kovatcheva-Datchary P Bäckhed F From dietary fiber to host physiology: short-chain fatty acids as key bacterial metabolites Cell 2016 165 1332 1345 10.1016/j.cell.2016.05.041 27259147
Koh, A., De Vadder, F., Kovatcheva-Datchary, P. & Bäckhed, F. From dietary fiber to host physiology: short-chain fatty acids as key bacterial metabolites. Cell 165, 1332–1345 (2016).27259147
5. Fachi JL Butyrate protects mice from clostridium difficile-induced colitis through an HIF-1-dependent mechanism Cell Rep. 2019 27 750 761.e7 10.1016/j.celrep.2019.03.054 30995474
Fachi, J. L. et al. Butyrate protects mice from clostridium difficile-induced colitis through an HIF-1-dependent mechanism. Cell Rep. 27, 750–761.e7 (2019).30995474
6. van der Hee B Wells JM Microbial regulation of host physiology by short-chain fatty acids Trends Microbiol 2021 29 700 712 10.1016/j.tim.2021.02.001 33674141
van der Hee, B. & Wells, J. M. Microbial regulation of host physiology by short-chain fatty acids. Trends Microbiol. 29, 700–712 (2021).33674141
7. Dang AT Marsland BJ Microbes, metabolites, and the gut-lung axis Mucosal Immunol. 2019 12 843 850 10.1038/s41385-019-0160-6 30976087
Dang, A. T. & Marsland, B. J. Microbes, metabolites, and the gut-lung axis. Mucosal Immunol. 12, 843–850 (2019).30976087
8. Rhee SH Pothoulakis C Mayer EA Principles and clinical implications of the brain–gut–enteric microbiota axis Nat. Rev. Gastroenterol. Hepatol. 2009 6 306 314 10.1038/nrgastro.2009.35 19404271
Rhee, S. H., Pothoulakis, C. & Mayer, E. A. Principles and clinical implications of the brain–gut–enteric microbiota axis. Nat. Rev. Gastroenterol. Hepatol. 6, 306–314 (2009).19404271
9. Bindari YR Gerber PF Centennial review: factors affecting the chicken gastrointestinal microbial composition and their association with gut health and productive performance Poult. Sci. 2022 101 101612 10.1016/j.psj.2021.101612 34872745
Bindari, Y. R. & Gerber, P. F. Centennial review: factors affecting the chicken gastrointestinal microbial composition and their association with gut health and productive performance. Poult. Sci. 101, 101612 (2022).34872745
10. Oviedo-Rondón EO Holistic view of intestinal health in poultry Anim. Feed Sci. Technol. 2019 250 1 8 10.1016/j.anifeedsci.2019.01.009
Oviedo-Rondón, E. O. Holistic view of intestinal health in poultry. Anim. Feed Sci. Technol. 250, 1–8 (2019).
11. Varmuzova, K. et al. Composition of gut microbiota influences resistance of newly hatched chickens to salmonella enteritidis infection. Front. Microbiol. 7, 957 (2016).
12. Diaz Carrasco JM Casanova NA Fernández Miyakawa ME Microbiota, gut health and chicken productivity: what is the connection? Microorganisms 2019 7 374 10.3390/microorganisms7100374 31547108
Diaz Carrasco, J. M., Casanova, N. A. & Fernández Miyakawa, M. E. Microbiota, gut health and chicken productivity: what is the connection? Microorganisms 7, 374 (2019).31547108
13. Rubio LA Possibilities of early life programming in broiler chickens via intestinal microbiota modulation Poult. Sci. 2019 98 695 706 10.3382/ps/pey416 30247675
Rubio, L. A. Possibilities of early life programming in broiler chickens via intestinal microbiota modulation. Poult. Sci. 98, 695–706 (2019).30247675
14. Sergeant MJ Extensive microbial and functional diversity within the chicken cecal microbiome PLOS One 2014 9 e91941 10.1371/journal.pone.0091941 24657972
Sergeant, M. J. et al. Extensive microbial and functional diversity within the chicken cecal microbiome. PLOS One 9, e91941 (2014).24657972
15. Broom LJ Kogut MH The role of the gut microbiome in shaping the immune system of chickens Vet. Immunol. Immunopathol. 2018 204 44 51 10.1016/j.vetimm.2018.10.002 30596380
Broom, L. J. & Kogut, M. H. The role of the gut microbiome in shaping the immune system of chickens. Vet. Immunol. Immunopathol. 204, 44–51 (2018).30596380
16. Saint-Martin V Quere P Trapp S Guabiraba R Uncovering the core principles of the gut-lung axis to enhance innate immunity in the chicken Front. Immunol. 2022 13 956670 10.3389/fimmu.2022.956670 36268022
Saint-Martin, V., Quere, P., Trapp, S. & Guabiraba, R. Uncovering the core principles of the gut-lung axis to enhance innate immunity in the chicken. Front. Immunol. 13, 956670 (2022).36268022
17. Verhagen JH Fouchier RAM Lewis N Highly pathogenic avian influenza viruses at the wild–domestic bird interface in europe: future directions for research and surveillance Viruses 2021 13 212 10.3390/v13020212 33573231
Verhagen, J. H., Fouchier, R. A. M. & Lewis, N. Highly pathogenic avian influenza viruses at the wild–domestic bird interface in europe: future directions for research and surveillance. Viruses 13, 212 (2021).33573231
18. Guitton, E. et al. Production of germ-free fast-growing broilers from a commercial line for microbiota studies. J. Vis. Exp. e61148 10.3791/61148 (2020)
19. Rampini SK Broad-range 16S rRNA gene polymerase chain reaction for diagnosis of culture-negative bacterial infections Clin. Infect. Dis. 2011 53 1245 1251 10.1093/cid/cir692 21976460
Rampini, S. K. et al. Broad-range 16S rRNA gene polymerase chain reaction for diagnosis of culture-negative bacterial infections. Clin. Infect. Dis. 53, 1245–1251 (2011).21976460
20. European Pharmacopoeia. European Pharmacopoeia (Ph. Eur.): Supplement 6.3, Sterility., Reference 01/2009:20601 (2009).
21. Japanese Pharmacopoeia, Japanese Pharmacopoeia. Japanese Pharmacopoeia (JP): the 4.06 Sterility Test as It Appeared in the Partial Revision of the JP 15th Edition, Made Official March 31, 2009, by the Ministry of Health, Labour and Welfare Ministerial Notification No. 190 (2009).
22. United States Pharmacopeia (USP). United States Pharmacopeia (USP): sterility tests as presented in pharmacopeial forum, interim revision announcement No. 6. 34 (6), December 1, 2008, Official on May 1, 2009 (2008).
23. Banks J Changes in the haemagglutinin and the neuraminidase genes prior to the emergence of highly pathogenic H7N1 avian influenza viruses in Italy Arch. Virol. 2001 146 963 973 10.1007/s007050170128 11448033
Banks, J. et al. Changes in the haemagglutinin and the neuraminidase genes prior to the emergence of highly pathogenic H7N1 avian influenza viruses in Italy. Arch. Virol. 146, 963–973 (2001).11448033
24. Rebel JM Highly pathogenic or low pathogenic avian influenza virus subtype H7N1 infection in chicken lungs: small differences in general acute responses Vet. Res. 2011 42 10 10.1186/1297-9716-42-10 21314972
Rebel, J. M. et al. Highly pathogenic or low pathogenic avian influenza virus subtype H7N1 infection in chicken lungs: small differences in general acute responses. Vet. Res. 42, 10 (2011).21314972
25. Trapp S Major contribution of the RNA-binding domain of NS1 in the pathogenicity and replication potential of an avian H7N1 influenza virus in chickens Virol. J. 2018 15 55 10.1186/s12985-018-0960-4 29587792
Trapp, S. et al. Major contribution of the RNA-binding domain of NS1 in the pathogenicity and replication potential of an avian H7N1 influenza virus in chickens. Virol. J. 15, 55 (2018).29587792
26. Ask B van der Waaij EH Glass EJ Bishop SC Modeling immunocompetence development and immunoresponsiveness to challenge in chicks Poult. Sci. 2007 86 1336 1350 10.1093/ps/86.7.1336 17575180
Ask, B., van der Waaij, E. H., Glass, E. J. & Bishop, S. C. Modeling immunocompetence development and immunoresponsiveness to challenge in chicks. Poult. Sci. 86, 1336–1350 (2007).17575180
27. Esnault E A novel chicken lung epithelial cell line: characterization and response to low pathogenicity avian influenza virus Virus Res. 2011 159 32 42 10.1016/j.virusres.2011.04.022 21557972
Esnault, E. et al. A novel chicken lung epithelial cell line: characterization and response to low pathogenicity avian influenza virus. Virus Res. 159, 32–42 (2011).21557972
28. Trapp S Shortening the unstructured, interdomain region of the non-structural protein NS1 of an avian H1N1 influenza virus increases its replication and pathogenicity in chickens J. Gen. Virol. 2014 95 1233 1243 10.1099/vir.0.063776-0 24694396
Trapp, S. et al. Shortening the unstructured, interdomain region of the non-structural protein NS1 of an avian H1N1 influenza virus increases its replication and pathogenicity in chickens. J. Gen. Virol. 95, 1233–1243 (2014).24694396
29. Fouchier RAM Detection of influenza A viruses from different species by PCR amplification of conserved sequences in the matrix gene J. Clin. Microbiol. 2000 38 4096 4101 10.1128/JCM.38.11.4096-4101.2000 11060074
Fouchier, R. A. M. et al. Detection of influenza A viruses from different species by PCR amplification of conserved sequences in the matrix gene. J. Clin. Microbiol. 38, 4096–4101 (2000).11060074
30. Chang PV Hao L Offermanns S Medzhitov R The microbial metabolite butyrate regulates intestinal macrophage function via histone deacetylase inhibition PNAS 2014 111 2247 2252 10.1073/pnas.1322269111 24390544
Chang, P. V., Hao, L., Offermanns, S. & Medzhitov, R. The microbial metabolite butyrate regulates intestinal macrophage function via histone deacetylase inhibition. PNAS 111, 2247–2252 (2014).24390544
31. Davie JR Inhibition of histone deacetylase activity by butyrate J. Nutr. 2003 133 2485S 2493S 10.1093/jn/133.7.2485S 12840228
Davie, J. R. Inhibition of histone deacetylase activity by butyrate. J. Nutr. 133, 2485S–2493S (2003).12840228
32. Hauser A Efficient extraction from mice feces for NMR metabolomics measurements with special emphasis on SCFAs Metabolites 2019 9 E55 10.3390/metabo9030055
Hauser, A. et al. Efficient extraction from mice feces for NMR metabolomics measurements with special emphasis on SCFAs. Metabolites 9, E55 (2019).
33. Paddison PJ Caudy AA Hannon GJ Stable suppression of gene expression by RNAi in mammalian cells Proc. Natl Acad. Sci. 2002 99 1443 1448 10.1073/pnas.032652399 11818553
Paddison, P. J., Caudy, A. A. & Hannon, G. J. Stable suppression of gene expression by RNAi in mammalian cells. Proc. Natl Acad. Sci. 99, 1443–1448 (2002).11818553
34. Del Vesco AP Jang HJ Monson MS Lamont SJ Role of the chicken oligoadenylate synthase-like gene during in vitro Newcastle disease virus infection Poult. Sci. 2021 100 101067 10.1016/j.psj.2021.101067 33752069
Del Vesco, A. P., Jang, H. J., Monson, M. S. & Lamont, S. J. Role of the chicken oligoadenylate synthase-like gene during in vitro Newcastle disease virus infection. Poult. Sci. 100, 101067 (2021).33752069
35. Martin M Cutadapt removes adapter sequences from high-throughput sequencing reads EMBnet. J. 2011 17 10 12 10.14806/ej.17.1.200
Martin, M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet. J. 17, 10–12 (2011).
36. Li B Dewey CN RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome BMC Bioinform. 2011 12 323 10.1186/1471-2105-12-323
Li, B. & Dewey, C. N. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinform. 12, 323 (2011).
37. Love MI Huber W Anders S Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2 Genome Biol. 2014 15 550 10.1186/s13059-014-0550-8 25516281
Love, M. I., Huber, W. & Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15, 550 (2014).25516281
38. Raudvere U g:Profiler: a web server for functional enrichment analysis and conversions of gene lists (2019 update) Nucleic Acids Res. 2019 47 W191 W198 10.1093/nar/gkz369 31066453
Raudvere, U. et al. g:Profiler: a web server for functional enrichment analysis and conversions of gene lists (2019 update). Nucleic Acids Res. 47, W191–W198 (2019).31066453
39. Zheng W An accurate and efficient experimental approach for characterization of the complex oral microbiota Microbiome 2015 3 48 10.1186/s40168-015-0110-9 26437933
Zheng, W. et al. An accurate and efficient experimental approach for characterization of the complex oral microbiota. Microbiome 3, 48 (2015).26437933
40. Escudié F FROGS: find, rapidly, OTUs with galaxy solution Bioinformatics 2018 34 1287 1294 10.1093/bioinformatics/btx791 29228191
Escudié, F. et al. FROGS: find, rapidly, OTUs with galaxy solution. Bioinformatics 34, 1287–1294 (2018).29228191
41. Mahé F Rognes T Quince C de Vargas C Dunthorn M Swarm: robust and fast clustering method for amplicon-based studies PeerJ 2014 2 e593 10.7717/peerj.593 25276506
Mahé, F., Rognes, T., Quince, C., de Vargas, C. & Dunthorn, M. Swarm: robust and fast clustering method for amplicon-based studies. PeerJ 2, e593 (2014).25276506
42. Quast C The SILVA ribosomal RNA gene database project: improved data processing and web-based tools Nucleic Acids Res. 2013 41 D590 D596 10.1093/nar/gks1219 23193283
Quast, C. et al. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 41, D590–D596 (2013).23193283
43. Hofmann T Schmucker S Characterization of chicken leukocyte subsets from lymphatic tissue by flow cytometry Cytom. Pt A 2021 99 289 300 10.1002/cyto.a.24214
Hofmann, T. & Schmucker, S. Characterization of chicken leukocyte subsets from lymphatic tissue by flow cytometry. Cytom. Pt A 99, 289–300 (2021).
44. Biddle A Stewart L Blanchard J Leschine S Untangling the genetic basis of fibrolytic specialization by Lachnospiraceae and Ruminococcaceae in diverse gut communities Diversity 2013 5 627 640 10.3390/d5030627
Biddle, A., Stewart, L., Blanchard, J. & Leschine, S. Untangling the genetic basis of fibrolytic specialization by Lachnospiraceae and Ruminococcaceae in diverse gut communities. Diversity 5, 627–640 (2013).
45. Glendinning L Chintoan-Uta C Stevens MP Watson M Effect of cecal microbiota transplantation between different broiler breeds on the chick flora in the first week of life Poult. Sci. 2022 101 101624 10.1016/j.psj.2021.101624 34936955
Glendinning, L., Chintoan-Uta, C., Stevens, M. P. & Watson, M. Effect of cecal microbiota transplantation between different broiler breeds on the chick flora in the first week of life. Poult. Sci. 101, 101624 (2022).34936955
46. Kubasova T Contact with adult hen affects development of caecal microbiota in newly hatched chicks PLoS One 2019 14 e0212446 10.1371/journal.pone.0212446 30840648
Kubasova, T. et al. Contact with adult hen affects development of caecal microbiota in newly hatched chicks. PLoS One 14, e0212446 (2019).30840648
47. Li Y Jin L Chen T The effects of secretory IgA in the mucosal immune system Biomed. Res. Int 2020 2020 2032057 31998782
Li, Y., Jin, L. & Chen, T. The effects of secretory IgA in the mucosal immune system. Biomed. Res. Int 2020, 2032057 (2020).31998782
48. Zang W Zheng X Structure and functions of cellular redox sensor HSCARG/NMRAL1, a linkage among redox status, innate immunity, DNA damage response, and cancer Free Radic. Biol. Med. 2020 160 768 774 10.1016/j.freeradbiomed.2020.09.016 32950687
Zang, W. & Zheng, X. Structure and functions of cellular redox sensor HSCARG/NMRAL1, a linkage among redox status, innate immunity, DNA damage response, and cancer. Free Radic. Biol. Med. 160, 768–774 (2020).32950687
49. Shao J Sheng H Amphiregulin promotes intestinal epithelial regeneration: roles of intestinal subepithelial myofibroblasts Endocrinology 2010 151 3728 3737 10.1210/en.2010-0319 20534719
Shao, J. & Sheng, H. Amphiregulin promotes intestinal epithelial regeneration: roles of intestinal subepithelial myofibroblasts. Endocrinology 151, 3728–3737 (2010).20534719
50. Kaspers, B., Schat, K. A., Göbel, T. & Vervelde, L. Avian Immunology. 3rd edn. (Elsevier Academic Press, London, 2022).
51. Liang X Macrophage FABP4 is required for neutrophil recruitment and bacterial clearance in Pseudomonas aeruginosa pneumonia FASEB J. 2019 33 3562 3574 10.1096/fj.201802002R 30462529
Liang, X. et al. Macrophage FABP4 is required for neutrophil recruitment and bacterial clearance in Pseudomonas aeruginosa pneumonia. FASEB J. 33, 3562–3574 (2019).30462529
52. Sutton K Visualisation and characterisation of mononuclear phagocytes in the chicken respiratory tract using CSF1R-transgenic chickens Vet. Res. 2018 49 104 10.1186/s13567-018-0598-7 30305141
Sutton, K. et al. Visualisation and characterisation of mononuclear phagocytes in the chicken respiratory tract using CSF1R-transgenic chickens. Vet. Res. 49, 104 (2018).30305141
53. Reese S Dalamani G Kaspers B The avian lung-associated immune system: a review Vet. Res. 2006 37 311 324 10.1051/vetres:2006003 16611550
Reese, S., Dalamani, G. & Kaspers, B. The avian lung-associated immune system: a review. Vet. Res. 37, 311–324 (2006).16611550
54. Flach H Mzb1 protein regulates calcium homeostasis, antibody secretion, and integrin activation in innate-like B cells Immunity 2010 33 723 735 10.1016/j.immuni.2010.11.013 21093319
Flach, H. et al. Mzb1 protein regulates calcium homeostasis, antibody secretion, and integrin activation in innate-like B cells. Immunity 33, 723–735 (2010).21093319
55. Koskela K Identification of a novel cytokine-like transcript differentially expressed in avian γδ T cells Immunogenetics 2004 55 845 854 10.1007/s00251-004-0643-8 14760511
Koskela, K. et al. Identification of a novel cytokine-like transcript differentially expressed in avian γδ T cells. Immunogenetics 55, 845–854 (2004).14760511
56. Mazewski, C., Perez, R. E., Fish, E. N. & Platanias, L. C. Type I interferon (IFN)-regulated activation of canonical and non-canonical signaling pathways. Front. Immunol. 11, 606456 (2020).
57. Pfister S., Kuettel V., Ferrero E. granulator: Rapid benchmarking of methods for *in silico* deconvolution of bulk RNA-seq data. R package version 1.12.0. Granulator. ([object Object]). 10.18129/B9.BIOC.GRANULATOR (2024).
58. Borowska D Highly multiplexed quantitative PCR-based platform for evaluation of chicken immune responses PLoS One 2019 14 e0225658 10.1371/journal.pone.0225658 31794562
Borowska, D. et al. Highly multiplexed quantitative PCR-based platform for evaluation of chicken immune responses. PLoS One 14, e0225658 (2019).31794562
59. Zhang L Liu C Jiang Q Yin Y Butyrate in energy metabolism: there is still more to learn Trends Endocrinol. Metab. 2021 32 159 169 10.1016/j.tem.2020.12.003 33461886
Zhang, L., Liu, C., Jiang, Q. & Yin, Y. Butyrate in energy metabolism: there is still more to learn. Trends Endocrinol. Metab. 32, 159–169 (2021).33461886
60. Salvi PS Cowles RA Butyrate and the intestinal epithelium: modulation of proliferation and inflammation in homeostasis and disease Cells 2021 10 1775 10.3390/cells10071775 34359944
Salvi, P. S. & Cowles, R. A. Butyrate and the intestinal epithelium: modulation of proliferation and inflammation in homeostasis and disease. Cells 10, 1775 (2021).34359944
61. Martin-Gallausiaux, C. et al. Butyrate produced by gut commensal bacteria activates TGF-beta1 expression through the transcription factor SP1 in human intestinal epithelial cells. Sci. Rep. 9, 2838 (2018).
62. Jiang Y Endothelial aquaporin-1 (AQP1) expression is regulated by transcription factor Mef2c Mol. Cells 2016 39 292 298 10.14348/molcells.2016.2223 26923194
Jiang, Y. et al. Endothelial aquaporin-1 (AQP1) expression is regulated by transcription factor Mef2c. Mol. Cells 39, 292–298 (2016).26923194
63. He J Li D Zhou J Zhu Y Yu B SP1-mediated upregulation of lncRNA LMCD1-AS1 functions a ceRNA for miR-106b-5p to facilitate osteosarcoma progression Biochem. Biophys. Res. Commun. 2020 526 670 677 10.1016/j.bbrc.2020.03.151 32248969
He, J., Li, D., Zhou, J., Zhu, Y. & Yu, B. SP1-mediated upregulation of lncRNA LMCD1-AS1 functions a ceRNA for miR-106b-5p to facilitate osteosarcoma progression. Biochem. Biophys. Res. Commun. 526, 670–677 (2020).32248969
64. Munier S A genetically engineered waterfowl influenza virus with a deletion in the stalk of the neuraminidase has increased virulence for chickens J. Virol. 2010 84 940 952 10.1128/JVI.01581-09 19889765
Munier, S. et al. A genetically engineered waterfowl influenza virus with a deletion in the stalk of the neuraminidase has increased virulence for chickens. J. Virol. 84, 940–952 (2010).19889765
65. Rooks MG Garrett WS Gut microbiota, metabolites and host immunity Nat. Rev. Immunol. 2016 16 341 352 10.1038/nri.2016.42 27231050
Rooks, M. G. & Garrett, W. S. Gut microbiota, metabolites and host immunity. Nat. Rev. Immunol. 16, 341–352 (2016).27231050
66. Glendinning L Watson KA Watson M Development of the duodenal, ileal, jejunal and caecal microbiota in chickens Anim. Microbiome 2019 1 17 10.1186/s42523-019-0017-z 33499941
Glendinning, L., Watson, K. A. & Watson, M. Development of the duodenal, ileal, jejunal and caecal microbiota in chickens. Anim. Microbiome 1, 17 (2019).33499941
67. Nie, K. et al. Roseburia intestinalis: a beneficial gut organism from the discoveries in genus and species. Front. Cell. Infection Microbiology 11, (2021).
68. Ruan G Roseburia intestinalis and its metabolite butyrate inhibit colitis and upregulate TLR5 through the SP3 signaling Pathway Nutrients 2022 14 3041 10.3390/nu14153041 35893896
Ruan, G. et al. Roseburia intestinalis and its metabolite butyrate inhibit colitis and upregulate TLR5 through the SP3 signaling Pathway. Nutrients 14, 3041 (2022).35893896
69. Kang X Roseburia intestinalis generated butyrate boosts anti-PD-1 efficacy in colorectal cancer by activating cytotoxic CD8 + T cells Gut 2023 72 2112 2122 10.1136/gutjnl-2023-330291 37491158
Kang, X. et al. Roseburia intestinalis generated butyrate boosts anti-PD-1 efficacy in colorectal cancer by activating cytotoxic CD8 + T cells. Gut. 72, 2112–2122 (2023).37491158
70. Trompette A Gut microbiota metabolism of dietary fiber influences allergic airway disease and hematopoiesis Nat. Med. 2014 20 159 166 10.1038/nm.3444 24390308
Trompette, A. et al. Gut microbiota metabolism of dietary fiber influences allergic airway disease and hematopoiesis. Nat. Med. 20, 159–166 (2014).24390308
71. Schroeder BO Bäckhed F Signals from the gut microbiota to distant organs in physiology and disease Nat. Med. 2016 22 1079 1089 10.1038/nm.4185 27711063
Schroeder, B. O. & Bäckhed, F. Signals from the gut microbiota to distant organs in physiology and disease. Nat. Med. 22, 1079–1089 (2016).27711063
72. Kogut MH The gut microbiota and host innate immunity: regulators of host metabolism and metabolic diseases in poultry?1 1Presented as a part of the Informal Nutrition Symposium “Metabolic Responses to Nutrition and Modifiers” at the Poultry Science Association’s annual meeting in Athens, Georgia, July 9, 2012 J. Appl. Poult. Res. 2013 22 637 646 10.3382/japr.2013-00741
Kogut, M. H. The gut microbiota and host innate immunity: regulators of host metabolism and metabolic diseases in poultry?1 1Presented as a part of the Informal Nutrition Symposium “Metabolic Responses to Nutrition and Modifiers” at the Poultry Science Association’s annual meeting in Athens, Georgia, July 9, 2012. J. Appl. Poult. Res. 22, 637–646 (2013).
73. Yitbarek A Commensal gut microbiota can modulate adaptive immune responses in chickens vaccinated with whole inactivated avian influenza virus subtype H9N2 Vaccine 2019 37 6640 6647 10.1016/j.vaccine.2019.09.046 31542262
Yitbarek, A. et al. Commensal gut microbiota can modulate adaptive immune responses in chickens vaccinated with whole inactivated avian influenza virus subtype H9N2. Vaccine 37, 6640–6647 (2019).31542262
74. Cowieson AJ Comparative biology of germ-free and conventional poultry Poult. Sci. 2022 101 102105 10.1016/j.psj.2022.102105 36057194
Cowieson, A. J. Comparative biology of germ-free and conventional poultry. Poult. Sci. 101, 102105 (2022).36057194
75. Mitsuhiro F Jun-ichi O Nutritional and physiological characteristics in germ-free chickens Comp. Biochem. Physiol. Part A: Physiol. 1994 109 547 556 10.1016/0300-9629(94)90193-7
Mitsuhiro, F. & Jun-ichi, O. Nutritional and physiological characteristics in germ-free chickens. Comp. Biochem. Physiol. Part A: Physiol. 109, 547–556 (1994).
76. den Besten G The role of short-chain fatty acids in the interplay between diet, gut microbiota, and host energy metabolism J. Lipid Res. 2013 54 2325 2340 10.1194/jlr.R036012 23821742
den Besten, G. et al. The role of short-chain fatty acids in the interplay between diet, gut microbiota, and host energy metabolism. J. Lipid Res. 54, 2325–2340 (2013).23821742
77. González Hernández MA Canfora EE Jocken JWE Blaak EE The short-chain fatty acid acetate in body weight control and insulin sensitivity Nutrients 2019 11 1943 10.3390/nu11081943 31426593
González Hernández, M. A., Canfora, E. E., Jocken, J. W. E. & Blaak, E. E. The short-chain fatty acid acetate in body weight control and insulin sensitivity. Nutrients 11, 1943 (2019).31426593
78. Morrison DJ Preston T Formation of short chain fatty acids by the gut microbiota and their impact on human metabolism Gut. Microbes 2016 7 189 200 10.1080/19490976.2015.1134082 26963409
Morrison, D. J. & Preston, T. Formation of short chain fatty acids by the gut microbiota and their impact on human metabolism. Gut. Microbes 7, 189–200 (2016).26963409
79. Ron-Harel N T cell activation depends on extracellular alanine Cell Rep. 2019 28 3011 3021.e4 10.1016/j.celrep.2019.08.034 31533027
Ron-Harel, N. et al. T cell activation depends on extracellular alanine. Cell Rep. 28, 3011–3021.e4 (2019).31533027
80. Li P Yin Y-L Li D Kim SW Wu G Amino acids and immune function Br. J. Nutr. 2007 98 237 252 10.1017/S000711450769936X 17403271
Li, P., Yin, Y.-L., Li, D., Kim, S. W. & Wu, G. Amino acids and immune function. Br. J. Nutr. 98, 237–252 (2007).17403271
81. Trompette A Dietary fiber confers protection against Flu by shaping Ly6c-patrolling monocyte hematopoiesis and CD8+T cell metabolism Immunity 2018 48 992 1005.e8 10.1016/j.immuni.2018.04.022 29768180
Trompette, A. et al. Dietary fiber confers protection against Flu by shaping Ly6c-patrolling monocyte hematopoiesis and CD8+T cell metabolism. Immunity 48, 992–1005.e8 (2018).29768180
82. Hanson C The relationship between dietary fiber intake and lung function in the National Health and Nutrition Examination surveys Ann. ATS 2016 13 643 650 10.1513/AnnalsATS.201509-609OC
Hanson, C. et al. The relationship between dietary fiber intake and lung function in the National Health and Nutrition Examination surveys. Ann. ATS 13, 643–650 (2016).
83. Marsland BJ Trompette A Gollwitzer ES The gut–lung axis in respiratory disease Ann. ATS 2015 12 S150 S156 10.1513/AnnalsATS.201503-133AW
Marsland, B. J., Trompette, A. & Gollwitzer, E. S. The gut–lung axis in respiratory disease. Ann. ATS 12, S150–S156 (2015).
84. Ayala, A. V. et al. Commensal bacteria promote type I interferon signaling to maintain immune tolerance. https://www.biorxiv.org/content/10.1101/2021.10.21.464743v2 (2022).
85. Shepard, J. D. et al. The Structure and immune regulatory implications of the ubiquitin-like tandem domain within an Avian 2’−5’ oligoadenylate synthetase-like protein. Front. Immunol. 12, 1664–3224 (2022).
86. Magor KE Defense genes missing from the flight division Dev. Comp. Immunol. 2013 41 377 388 10.1016/j.dci.2013.04.010 23624185
Magor, K. E. et al. Defense genes missing from the flight division. Dev. Comp. Immunol. 41, 377–388 (2013).23624185
87. dos Santos PF Mansur DS Beyond ISGlylation: functions of free intracellular and extracellular ISG15 J. Interferon Cytokine Res. 2017 37 246 253 10.1089/jir.2016.0103 28467275
dos Santos, P. F. & Mansur, D. S. Beyond ISGlylation: functions of free intracellular and extracellular ISG15. J. Interferon Cytokine Res. 37, 246–253 (2017).28467275
88. Abt MC Commensal bacteria calibrate the activation threshold of innate antiviral Immun. Immun. 2012 37 158 170
Abt, M. C. et al. Commensal bacteria calibrate the activation threshold of innate antiviral. Immun. Immun. 37, 158–170 (2012).
89. Bradley KC Microbiota-driven tonic interferon signals in lung stromal cells protect from influenza virus infection Cell Rep. 2019 28 245 256.e4 10.1016/j.celrep.2019.05.105 31269444
Bradley, K. C. et al. Microbiota-driven tonic interferon signals in lung stromal cells protect from influenza virus infection. Cell Rep. 28, 245–256.e4 (2019).31269444
90. Wirusanti NI Baldridge MT Harris VC Microbiota regulation of viral infections through interferon signaling Trends Microbiol. 2022 30 778 792 10.1016/j.tim.2022.01.007 35135717
Wirusanti, N. I., Baldridge, M. T. & Harris, V. C. Microbiota regulation of viral infections through interferon signaling. Trends Microbiol. 30, 778–792 (2022).35135717
91. Evseev D Magor KE Innate immune responses to avian influenza viruses in ducks and chickens Vet. Sci. 2019 6 5 30634569
Evseev, D. & Magor, K. E. Innate immune responses to avian influenza viruses in ducks and chickens. Vet. Sci. 6, 5 (2019).30634569
92. Han X Bertzbach LD Veit M Mimicking the passage of avian influenza viruses through the gastrointestinal tract of chickens Vet. Microbiol. 2019 239 108462 10.1016/j.vetmic.2019.108462 31767100
Han, X., Bertzbach, L. D. & Veit, M. Mimicking the passage of avian influenza viruses through the gastrointestinal tract of chickens. Vet. Microbiol. 239, 108462 (2019).31767100
93. Ichinohe T Microbiota regulates immune defense against respiratory tract influenza A virus infection Proc. Natl Acad. Sci. USA 2011 108 5354 5359 10.1073/pnas.1019378108 21402903
Ichinohe, T. et al. Microbiota regulates immune defense against respiratory tract influenza A virus infection. Proc. Natl Acad. Sci. USA 108, 5354–5359 (2011).21402903
94. Yitbarek A Gut microbiota modulates type I interferon and antibody-mediated immune responses in chickens infected with influenza virus subtype H9N2 Benef. Microbes 2018 9 417 427 10.3920/BM2017.0088 29380643
Yitbarek, A. et al. Gut microbiota modulates type I interferon and antibody-mediated immune responses in chickens infected with influenza virus subtype H9N2. Benef. Microbes 9, 417–427 (2018).29380643
95. Chemudupati M Butyrate reprograms expression of specific interferon-stimulated genes J. Virol. 2020 94 e00326 20 10.1128/JVI.00326-20 32461320
Chemudupati, M. et al. Butyrate reprograms expression of specific interferon-stimulated genes. J. Virol. 94, e00326–20 (2020).32461320
96. Klampfer L Huang J Sasazuki T Shirasawa S Augenlicht L Inhibition of interferon gamma signaling by the short chain fatty acid butyrate Mol. Cancer Res. 2003 1 855 862 14517348
Klampfer, L., Huang, J., Sasazuki, T., Shirasawa, S. & Augenlicht, L. Inhibition of interferon gamma signaling by the short chain fatty acid butyrate. Mol. Cancer Res. 1, 855–862 (2003).14517348
97. Narita T Weinert BT Choudhary C Functions and mechanisms of non-histone protein acetylation Nat. Rev. Mol. Cell Biol. 2019 20 156 174 10.1038/s41580-018-0081-3 30467427
Narita, T., Weinert, B. T. & Choudhary, C. Functions and mechanisms of non-histone protein acetylation. Nat. Rev. Mol. Cell Biol. 20, 156–174 (2019).30467427
98. Meslin C Expanding duplication of free fatty acid receptor-2 (GPR43) genes in the chicken genome Genome Biol. Evol. 2015 7 1332 1348 10.1093/gbe/evv072 25912043
Meslin, C. et al. Expanding duplication of free fatty acid receptor-2 (GPR43) genes in the chicken genome. Genome Biol. Evol. 7, 1332–1348 (2015).25912043
99. Tan NY Khachigian LM Sp1 phosphorylation and its regulation of gene transcription Mol. Cell Biol. 2009 29 2483 2488 10.1128/MCB.01828-08 19273606
Tan, N. Y. & Khachigian, L. M. Sp1 phosphorylation and its regulation of gene transcription. Mol. Cell Biol. 29, 2483–2488 (2009).19273606
100. Marin M Karis A Visser P Grosveld F Philipsen S Transcription factor Sp1 is essential for early embryonic development but dispensable for cell growth and differentiation Cell 1997 89 619 628 10.1016/S0092-8674(00)80243-3 9160753
Marin, M., Karis, A., Visser, P., Grosveld, F. & Philipsen, S. Transcription factor Sp1 is essential for early embryonic development but dispensable for cell growth and differentiation. Cell 89, 619–628 (1997).9160753
101. Kantor B Ma H Webster-Cyriaque J Monahan PE Kafri T Epigenetic activation of unintegrated HIV-1 genomes by gut-associated short chain fatty acids and its implications for HIV infection Proc. Natl Acad. Sci. USA 2009 106 18786 18791 10.1073/pnas.0905859106 19843699
Kantor, B., Ma, H., Webster-Cyriaque, J., Monahan, P. E. & Kafri, T. Epigenetic activation of unintegrated HIV-1 genomes by gut-associated short chain fatty acids and its implications for HIV infection. Proc. Natl Acad. Sci. USA 106, 18786–18791 (2009).19843699
102. Tolnay, A.-E. et al. Extrapulmonary tissue responses in cynomolgus macaques (Macaca fascicularis) infected with highly pathogenic avian influenza A (H5N1) virus. Arch. Virol. 155, 905–914 (2010).
103. Pang Z MetaboAnalyst 5.0: narrowing the gap between raw spectra and functional insights Nucleic Acids Res. 2021 49 W388 W396 10.1093/nar/gkab382 34019663
Pang, Z. et al. MetaboAnalyst 5.0: narrowing the gap between raw spectra and functional insights. Nucleic Acids Res. 49, W388–W396 (2021).34019663
104. Hulsegge I Kommadath A Smits MA Globaltest and GOEAST: two different approaches for Gene Ontology analysis BMC Proc. 2009 3 S10 10.1186/1753-6561-3-S4-S10 19615110
Hulsegge, I., Kommadath, A. & Smits, M. A. Globaltest and GOEAST: two different approaches for Gene Ontology analysis. BMC Proc. 3, S10 (2009).19615110
105. Kolberg L g:Profiler—interoperable web service for functional enrichment analysis and gene identifier mapping (2023 update) Nucleic Acids Res. 2023 51 W207 W212 10.1093/nar/gkad347 37144459
Kolberg, L. et al. g:Profiler—interoperable web service for functional enrichment analysis and gene identifier mapping (2023 update). Nucleic Acids Res. 51, W207–W212 (2023).37144459
