key: cord-286779-si3qml42
authors: Li, Hai-yan; Zhang, Hong-lei; zhao, Fu-jie; Wang, Shi-qiong; Wang, Zhi-xiang; Wei, Zhan-yong
title: Modulation of Gut Microbiota, Short-Chain Fatty Acid Production, and Inflammatory Cytokine Expression in the Cecum of Porcine Deltacoronavirus-Infected Chicks
date: 2020-06-04
journal: Front Microbiol
DOI: 10.3389/fmicb.2020.00897
sha: 
doc_id: 286779
cord_uid: si3qml42

Porcine deltacoronavirus (PDCoV) is a novel swine enteropathogenic coronavirus that causes watery diarrhea and induces proinflammatory cytokine responses in piglets. Our previous research showed that the specific-pathogen-free (SPF) chicks exhibited mild diarrhea and low fecal viral shedding, along with cecum lesions after PDCoV infection. Disturbances in the homeostasis of the gut microbiota have been associated with various diseases. We aimed to explore the effects of PDCoV infection on chick gut microbiota, short-chain fatty acid (SCFAs) production, and inflammatory cytokine expression in chicks, and also to investigate the relationship between gut microbiota and SCFAs or inflammatory cytokine expression of the PDCoV-infected chicks. Results obtained using 16S rRNA sequencing showed that infection with PDCoV strain HNZK-02 significantly altered the composition of chick gut microbiota, with the reduced abundance of Eisenbergiella and Anaerotruncus genera at 5 days post-inoculation (dpi) (P < 0.05), and an increased abundance of Alistipes genus at 17 dpi (P < 0.05). The production of SCFAs in the cecum of PDCoV HNZK-02–infected chicks, including acetic acid, propionic acid, and butyric acid, decreased in all cases. The expression of inflammatory cytokines (interferon-γ, tumor necrosis factor-α, and interleukin-10) was increased in the cecum tissue and serum of the PDCoV HNZK-02–infected chicks when detected by quantitative real-time polymerase chain reaction and enzyme-linked immunosorbent assay, respectively. Further analysis showed significant correlation between bacterial genera and SCFAs or inflammatory cytokines expression in cecum of the PDCoV infected chicks. These findings might provide new insight into the pathology and physiology of PDCoV in chicks.

Porcine deltacoronavirus (PDCoV) is an enveloped, positivesense single-stranded RNA virus, with a genome length of approximately 25 kb (Zhang, 2016) . This virus was originally found in Hong Kong in 2012 (Woo et al., 2012) , and the first outbreak of PDCoV-related diarrhea in swine herds was reported in the United States in 2014 (Wang et al., 2014) . As an important enteropathogen in pigs, PDCoV can cause acute diarrhea, vomiting, and dehydration in neonatal piglets (Chen et al., 2014; Hu et al., 2016) . Moreover, the infected pigs are characterized by thin and transparent intestinal walls and accumulation of large amounts of yellow fluid in the intestinal lumen (Jung et al., 2015; . Further, PDCoV-infected piglets show signs of proinflammatory cytokine responses during acute infection (Jung et al., 2018) . PDCoV has caused huge economic losses for the pig industries. Apart from infection of gnotobiotic (Gn) and conventional pigs, PDCoV was also reported to have limited ability to infect Gn calves (Jung et al., 2017) . Recent studies showed the susceptibility of specific-pathogen-free (SPF) chicks to PDCoV infection, with the clinical syndrome of mild diarrhea and slight lesions in cecum, which was weaker than that of piglets infected with PDCoV (Liang et al., 2019; Boley et al., 2020) . However, the reason why PDCoV infection can cause intestinal damage is still unknown.

The chicken gastrointestinal tract is an important site for immune cell development, which not only regulates gut microbiota, but also maintains extraintestinal immunity. Intestinal commensal microbes play a crucial role in gut homeostasis through mutually beneficial interactions with the host immune system (Deriu et al., 2016) . Once the intestinal mucosal barrier and microbiota are destroyed, intestinal inflammation will occur. In turn, the host's inflammatory response will elicit changes in the microbial community structure (Deriu et al., 2016) . The gut microbiota produces various metabolites, including short-chain fatty acids (SCFAs), which can regulate the host antiviral immune response (Budden et al., 2017) . Propionate and butyrate have also been reported to induce differentiation of T-regulatory cells, assisting in control of intestinal inflammation (Donohoe et al., 2014; Louis et al., 2014) .

In chicks, recent studies have emphasized the key roles of intestinal microbiota in shaping immunity against viral diseases, including avian influenza, Marek's disease, Newcastle disease (ND), and infectious bursal disease (IBD). These studies not only reflected connections between gut microbiota and distal organs in regulatory functions, but also emphasized the interaction between intestinal microbiota and viral infections and their impacts on immune regulations (Lillehoj and Trout, 1996; Balamurugan and Kataria, 2006; Li, 2007; Chen et al., 2017; Yitbarek et al., 2018) . It was reported that different commensal bacteria had their own unique role against viral infection by modulating diverse immune mechanisms (Abt et al., 2012) . In the H9N2 avian influenza virus infected chicks, increased levels of interferon gamma (IFN-γ). tumor necrosis factor alpha (TNFα), and interleukin-17A (IL-17A) led to intestinal microflora dysbiosis (Li H. et al., 2018) . However, there are no studies emphasizing how PDCoV affects the gut microbiota in chicks.

Thus, in the current study we aimed to observe the intestinal microbiota, SCFAs, and inflammation responses to PDCoV infection in chicks and elucidate the potential associations among them, which will provide the theoretical basis for investigation of the etiology and pathogenesis of PDCoV.

The virulent PDCoV strain-HNZK-02 (GenBank accession number MH708123) was isolated and identified by our laboratory and propagated in LLC-porcine kidney (LLC-PK) cells. Passage 5 of this strain was used in current study. Virus propagation and virus titers were performed as described previously (Liang et al., 2019) .

The SPF chicks were purchased from Jinan SIPAFAS Poultry Co. Ltd. in Jinan, China. Twenty-four-day-old SPF chicks with similar body weight (about 50 g/chick) were randomly divided into two groups, the control (mock) and PDCoV HNZK-02infected groups. All chicks were housed in the same facility but in different rooms under biosafety conditions and allowed free access to water and feed during the experiment. Feed was formulated to meet or exceed the National Research Council nutrient requirements for chicks (National Research Council [NRC] , 1994; Table 1 ). Each chick was intragastrically inoculated with PDCoV strain HNZK-02 or Eagle's Minimum Essential Medium (MEM) with 300 µL/chick (10.5 log10 GE/chick) (n = 10/per group).

Based on our previous study (Liang et al., 2019) , three chicks in each group were selected randomly for necropsy at 5 and 17 dpi, respectively. The cecal contents were collected for analysis of gut microbiota and SCFAs. The cecum tissue was collected for inflammatory cytokine mRNA measurements. The serum was collected for inflammatory cytokine detection by enzyme-linked immunosorbent assay (ELISA).

Viral RNA was extracted from cecal content suspensions using the TRizol TM method (Invitrogen, Carlsbad, CA, United States) according to the manufacturer's instructions. The viral RNA was further reverse transcribed into cDNA. Viral RNA titers were determined using quantitative real-time polymerase chain reaction (qRT-PCR) as reported previously (Liang et al., 2019) . The detection limit of the qRT-PCR was 4.6 log10 GE/mL. 

The SCFAs in cecum were determined using gas chromatography (GC) according to a previously described method (Stewart et al., 2010) . The samples were placed on a DB-WAX column (30 m long × 0.53 mm diameter and 1.00 µm film thickness) and were separated by using a TRACE TM 1300 GC with flame ionization detector (FID). The temperature program was 100 • C for 0.5 min, then raised to 170 • C at 8 • C/min and held for 1 min, then raised to 200 • C at 40 • C/min and held for 2 min. Samples were run with a 30:1 split ratio and a 5.0 mL/min column flow. High-purity hydrogen was used as the carrier gas. The temperatures of the injector and detector were 250 and 270 • C, respectively.

The levels of expression of the inflammatory cytokines (IFN-γ, TNF-α, and IL-10) in the cecum of SPF chicks were quantified by qRT-PCR. Cecum tissue samples were homogenized by the TISSUELYSER-24 (Jingxin, Shanghai, China). Total RNAs were extracted from the supernatant of the cecum tissue lysate by TRizol TM Reagent Kit (Invitrogen, Carlsbad, CA, United States) and were quantified by NanoDrop 2000 spectrophotometer (Thermo Fisher Scientific, Waltham, MA, United States), then further reverse transcribed into cDNA. The total reaction volume of 20 µL contained 2.0 µL of DNA, 10.0 µL of SYBR Green qPCR Mix (Takara Bio Inc., Tokyo, Japan), 7 µL of H 2 O, and 0.5 µL of each specific primer. Primers were synthesized according to previous reports (Emami et al., 2019) and are shown in Table 2 . Reactions were performed under the following conditions: one cycle of preincubated samples at 95 • C for 1 min and 39 cycles of amplification samples at 95 • C for 5 s, 60 • C for 30 s, 72 • C for 30 s. Target gene expression of each sample was normalized using glyceraldehyde 3-phosphate dehydrogenase (GAPDH). The relative expression of the target gene for each sample was calculated using the 2 − Ct method (Livak and Schmittgen, 2001) . The expression levels of the IFN-γ, TNF-α and IL-10 in serum were detected using ELISA Kit from Nanjing Jian cheng Biological Institute (Nanjing, China). Serum samples were obtained and measured as previously described (Huang et al., 2015) . Absorbance was read at 450 nm using a microplate spectrophotometer. The inflammatory cytokine concentrations were calculated according to a best-fit standard curve.

Raw fastq files were demultiplexed, quality-filtered by Trimmomatic, and merged by FLASH (Magoč and Salzberg, 2011) . Sequences with ambiguous bases, shorter than 200 bp, and sequences with an average mass less than 25 were removed. The high-quality reads were clustered into operational taxonomic units (OTUs) with 97% similarity cutoff using UPARSE (version 7.11) and used for further analysis of the Venn diagram, alpha diversity indices (Shannon and Simpson), and richness estimators (ACE and Chao). Beta diversity measures were based on weighted-uniFrac distance analysis and shown by principal coordinate analysis (PCoA), which was conducted to assess the relationships among the different groups. LEfSe was applied to recognize which bacterial taxa contribute to the differences between the two groups (Segata et al., 2011) , and the minimal threshold of linear discriminant analysis (LDA) was set at 3.5. The difference between specific taxa was analyzed using one-way analysis of variance (ANOVA) followed by post hoc t-tests. Data analysis of inflammatory cytokines and SCFAs was performed using GraphPad Prism (version 6.01, La Jolla, CA, United States). A P < 0.05 or 0.01 was considered statistically significant or highly significant, which is indicated as follows: * P < 0.05, * * P < 0.01.

SPF chicks, intragastrically inoculated with PDCoV HNZK-02 at 4 days old, showed mild diarrhea at 5 dpi. PDCoV RNA was detected in the intestinal contents by qRT-PCR at 5 and 17 dpi while non-infected control chicks appeared normal throughout the course of the experiments. PDCoV HNZK-02 infection had The operational taxonomic units (OTUs) were defined at a 97% similarity level. The coverage percentage (Good's) and richness estimators (ACE and Chao) and diversity indices (Shannon and Simpson) were calculated using Good's method and the mothur program, respectively. In the same column, values with the same superscript letter are not significantly different (P > 0.05); those with different superscript letters differ significantly (P < 0.05).

FIGURE 1 | Venn map. Venn map of shared OTUs based on the sequences with more than 97% similarity (n = 3). The mock group at 5 dpi (mock-5d), the PDCoV HNZK-02 group at 5 dpi (HNZK-02-5d), the mock at 17 dpi (mock-17d), and the PDCoV HNZK-02 group at 17 dpi (HNZK-02-17d).

no effect on feed consumption compared to controls during the whole experimental period.

A total of 12 cecal concent samples, the control chicks (n = 3), and PDCoV HNZK-02 groups (n = 3) at 5 and 17 dpi, respectively, were evaluated using 16S rRNA gene sequencing. A total of 663,736 sequences with a median read length of 429 bp (range from 274 to 494 bp) were collected. The total number of unique sequences from the four groups was 416. The sequence and species-level OTUs of each group are shown in Table 3 . The Venn diagram showed that 237 OTUs of the total gut microbial richness (416) were shared among all the sequenced samples. Two hundred and fifty-nine OTUs were shared between the samples of the mock and PDCoVHNZK-02 groups with a total of 343 at 5 dpi. Three hundred and seventeen OTUs of the total gut microbial richness (375) were shared between the samples of mock and PDCoV HNZK-02 groups at 17 dpi (Figure 1 ).

Alpha diversity estimators of community are shown in Table 3 . Gut microbial richness, according to the Ace and Chao indexes, was significantly increased in the PDCoV HNZK-02 group when compared with that of the mock group at 5 dpi (P < 0.05), while there was no significant difference between the two groups at 17 dpi (P > 0.05).

In Shannon and Simpson alpha diversity indexes, the differences between the PDCoV HNZK-02 and the mock groups at 5 and 17 dpi did not reach statistical significance (P > 0.05).

Beta diversity analysis showed that there were some similarities in the microbial composition between the mock and the PDCoV HNZK-02 groups at 5 dpi. However, the microbial composition from chicks with PDCoV inoculation and the controls could be divided into two different clusters at 17 dpi, indicating that the cecal microbial community structure and composition in the control group were significantly different from those in the PDCoV HNZK-02 group at 17 dpi (Figure 2) .

To further investigate whether PDCoV HNZK-02 infected chicks experienced any significant alteration in gut microbiota, we analyzed the relative abundance of microbiota at the phylum, family, and genus levels. At 5 dpi, genus level analysis showed the taxonomic abundance was altered in the PDCoV HNZK-02 group, which was characterized by lower abundance of Eisenbergiella and Anaerotruncus (P < 0.05). There were no significant differences at the levels of phylum and family between the mock and the PDCoV HNZk-02 groups (Figure 3A) . At 17 dpi, the ratio of Firmicutes to Bacteroidetes in the PDCoV HNZk-02 group was lower than that of the mock group. At the family level, the proportion of Rikenellaceac in the chicks inoculated with PDCoV HNZK-02 significantly increased (P < 0.05). At the genus level, Alistipes and the unclassified f-Ruminococcaceae were significantly increased (P < 0.05), whereas the proportion of Ruminococcaceae-UGG-014 was significantly decreased in the PDCoV HNZK-02 group when compared to that of the mock group at 17 dpi (P < 0.01) (Figure 3B) .

The characteristics of bacterial taxonomic abundant in the collected samples were further analyzed using LEfSe algorithm (Figure 4) . LEfSe analysis revealed that there existed main 8 and 33 discriminative features from phylum to genus (LDA score >3) between the mock and the PDCoV HNZK-02-inoculated chicks at 5 and 17 dpi (Figures 4A,B) , respectively. More specifically, Eisenbergiella and Anaerotruncus were the most differential microbiota in the mock at 5 dpi. Moraxellaceae, Acinetobacter, Pseudomonadates, Eubacterium, and norank-f Lachnospiraceae in the PDCoV HNZK-02 group increased when compared to that in the mock group. At 17 dpi, we identified 10 microbial biomarkers from the PDCoV HNZK-02 group, including Rikenellaceae, Bacteroidia, Bacteroidetes, Bacteroidales, Alistipes, unclassified-f-Ruminocococcaceae, Sphingomonadaceae, Novosphingobium, and Lachnoclostridium, FIGURE 2 | Beta diversity measures in cecal microbiota of chickens (n = 3). The scatterplot of the principal coordinate analysis (PCoA) scores shows four different clusters (weighted-uniFrac distance), the blue represents the PDCoV HNZK-02 group at 5 dpi (HNZK-02-5d), the red represents the mock group at 5 dpi (mock-5d), the purple represents the PDCoV HNZK-02 group at 17 dpi (HNZK-02-17d), and the green represents the mock group at 17 dpi (mock-17d).

Frontiers in Microbiology | www.frontiersin.org FIGURE 3 | Significant differences analysis in cecal microbiota of chickens. (A) Statistical comparison of the relative microbial abundance from phylum to genus between the mock and the PDCoV HNZK-02 groups at 5 dpi (n = 3). (B) Statistical comparison of the relative microbial abundance from phylum to genus between the mock and the PDCoV HNZK-02 groups at 17 dpi (n = 3). **P < 0.01; *P < 0.05.

LEfSe identifies the taxa with the greatest differences in abundance from phylum to genus between the mock and the PDCoV HNZK-02 groups at 5 dpi (n = 3). (B) LEfSe identifies the taxa with the greatest differences in abundance from phylum to genus between the mock and the PDCoV HNZK-02 groups at 17 dpi (n = 3). Linear discriminant analysis (LDA) effect size (LefSe) analysis showing those abundance from phylum to genus between groups; only taxa of an LDA significant threshold of 3.0 are shown.

which could be distinguished from the mock group. 13 clades were more abundant in the mock group than that in the PDCoV HNZK-02 group, including Firmicutes,

Ruminocococaceae-UCG-014, Clostridia, norank-f-Clostridiales-vadinBB60-group, Clostridiales-vadinBB60-group, Anaerostipe, Eisenbergiella, Coriobacteriaceae, 

Coriobacteriales, Actinobacteria, and Actinobacteria. There were bigger differences on microbiota compositions between the early and later stages of PDCoV infection (at 5 and 17 dpi, respectively).

Compared to the control, there was a decreasing trend in the amounts of acetic acid (AA), propionic acid (PA), and butyric acid (BA) in the PDCoV HNZK-02 group at 5 and 17 dpi. Furthermore, PA and AA levels in the PDCoV HNZK-02 group were greatly decreased at 5 and 17 dpi, respectively (P < 0.05) (Figure 5 ).

As shown in Figure 6 , the expression of IFN-γ (P < 0.01), TNFα (P < 0.05) and IL-10 (P < 0.01) was significantly upregulated at 5 dpi compared with that in the mock chicks (Figures 6A-C) . Furthermore, IFN-γ expression was significantly upregulated (Figure 6A ) at 17 dpi (P < 0.01).

To evaluate if PDCoV infection could affect IFN-γ, TFN-α, and IL-10 secretion in serum, the concentrations of IFN-γ, TFNα, and IL-10 in serum were tested with ELISA kits. Our results showed that the changes of the inflammatory cytokines in serum were similar with those in the cecum tissue. The levels of IFNγ and TNF-α were approximately two and three times as high as in the PDCoV HNZK-02 group when compared to that of the mock group at 5 dpi (P < 0.01). The level of IL-10 in the PDCoV infection group was also increased at 5 dpi (P < 0.05). The serum levels of IFN-γ in the PDCoV-infected chicks were higher than in the mock group at 17 dpi (P < 0.01), and no significant differences were observed for the secretion of TNFα and IL-10 between the mock and the PDCoV HNZK-02 groups at 17 dpi (Figures 7A-C) .

Spearman's correlation analysis was conducted between the top 15 bacterial genera and the environmental factors (SCFAs, proinflammatory cytokines), and was directly reflected by a heatmap (Figure 8) . The threshold | R| > 0.8 was considered as having a correlation. The results indicated that Eisenbergiella was negatively correlated with TNFα (P < 0.05) and IFN-γ (P < 0.01) expression levels, and was positively correlated with PA (P < 0.01) at 5 dpi ( Figure 8A) . The Ruminiclostridium_9 also showed a positively correlation with the IL-10 expression (P < 0.05), while Norank_f_ Ruminococcaceae was negatively correlated with BA (P < 0.001). At 17 dpi (Figure 8B) , Alistipes and unclassified_f_ Ruminococcaceae showed an extremely significant positive correlation with IFN-γ secretion (P < 0.05), but Eisenbergiella and Norank_f_Clostridiales_vadinbb60_group were negatively correlated with it (P < 0.05). The level of PA was associated with the increase in the abundance of Ruminiclostridium_5 (P < 0.05). Ruminococcus and Eisenbergiella were positively correlated with AA (P < 0.05), while Alistipes showed a negative relationship with it (P < 0.05). Eisenbergiella, unclassified_f_Lachnospiraceae, Ruminococcaceae_UCG-014, Norank_f_Clostridiales_vadinbb60_group, and Ruminiclostr idium_5 were positively correlated with BA.

Gut microflora play an important role in shaping immunity against viral diseases in host. Disruption of microbial homeostasis is associated with a variety of pathological states, which helps the establishment of acute viral infections in chickens (Abaidullah et al., 2019) . In this study, the infection model of PDCoV was developed on SPF chicks. Our results demonstrated similar clinical diseases, and the dynamics of the virus shedding in cecum were observed in all SPF chicks inoculated with the PDCoV HNZK-02 strains when compared with our previous experiment (Liang et al., 2019) . The influence of PDCoV infection on cecum microbiota of SPF chicken was evaluated, PDCoV HNZK-02 infection significantly altered the gut microbiota composition and decreased SCFAs products in chicks' cecum. Simultaneously, the expression levels of inflammatory cytokines (IFN-γ, IL-10, and TNF-α) in serum and cecum tissue were significantly increased, resulting in the inflammatory response in the chicken. The analysis showed a significant correlation between bacterial genera and SCFAs or inflammatory cytokines.

The diversity of indigenous intestinal microbiota is one of the key factors in resisting the colonization of invading pathogens (Keesing et al., 2010) . Several studies have reported that the diarrhea-relating coronaviruses could influence the intestinal microbiota diversity in pigs (Liu et al., 2015; Song et al., 2017 ; Huang et al., 2018; Tan et al., 2019) . Indeed, many viral agents have been shown to alter the intestinal microbiota in chickens, such as AIV (Yitbarek et al., 2018) , Marek's disease virus (MDV) (Lillehoj and Trout, 1996) , NDV (Li, 2007) , and IBDV (Li L. et al., 2018) . In PDCoV-infected chicks, the species richness indices were increased at 5 dpi, and there were no significant changes in the overall α-diversity of cecal microbiota at 17 dpi. It is possible that the reduced immune response may account for the change of microbial diversity. From the PCoA results, microbiota structure in the PDCoV infected group was clearly distinguished from that of the mock group at 17 dpi, which also reflected SPF chicks' resistance to PDCoV.

From microbial community profiling, the difference between the groups was analyzed at the phylum, family, and genus level. Eisenbergiella and Anaerotruncus genera were significantly decreased in the PDCoV HNZK-02 at 5 dpi, indicating that PDCoV HNZK-02 infection had a marked influence on chicks' cecal microbiota. According to previous studies, the decrease of the ratio of Firmicutes-to-Bacteroidetes was observed in mice with diabetes (Wen et al., 2008) and in some patients with Crohn's disease and ulcerative colitis (Frank et al., 2007) . In our results, a decreased ratio of Firmicutes-to-Bacteroidetes occurred in the PDCoV HNZK-02 group at 17 dpi. More remarkably, the Alistipes, considered an opportunistic pathogen (Pandit et al., 2018; Zhang et al., 2018) , was significantly increased in the PDCoV HNZK-02 group at 17 dpi. We used linear discriminant analysis (LDA) of effect size (LEfSe) to confirm the taxa that most likely explains the differences between the PDCoV HNZK-02 and mock groups. These findings provide compelling evidence that PDCoV HNZK-02 induced microbiota imbalance and showed that recovery was difficult in a short time period after PDCoV HNZK-02 infection.

Production of SCFAs, as the end products of protein and carbohydrate fermentation (Cummings and Macfarlane, 1991) , can easily be affected by the status of the gut microbiota (Lloyd-Price et al., 2019) . Studies have shown that Lactobacillus salivarius and L. agilis could increase propionate and butyrate contents in cecum of chicks (Meimandipour et al., 2010) . The changes in butyrate and propionate can influence intestinal physiology and immune function, while acetate can alter lipid metabolism (Macfarlane and Macfarlane, 2011) . Our data showed that the changes of gut microbiota composition contributed to the decreased SCFAs, which reduced this protective effect of SCFA in the cecum. However, the interaction of gut microbiota with SCFAs to respond to PDCoV infection needs further studies.

Inflammatory immune responses in the gut can alter the gut luminal environment in a way that favors dysbiosis (Winter et al., 2013) . Chicks with dysbiosis are more prone to acute viral infection (Abaidullah et al., 2019) , which may in turn further affect the pathogen infectivity. Eisenbergiella is strongly correlated with increased levels of TNF-α and IFN-γ, which help in modulating the functional activities of the cells of the immune system (Liu et al., 2010) , and to block virus replication and prevent clinical disease (Bradley, 2008; Praveena et al., 2010) . This may be the reason that Eisenbergiella was decreased in the PDCoV HNZK-02 group at 5 dpi. Moreover, the observation is similar to the increased proinflammatory (TNF-α) cytokine responses of 10-day-old Gn pigs to acute PDCoV infection (Jung et al., 2018) . Alistipes was strong positively correlated with IFN-γ levels, and the numbers of Ruminococcaceae_UCG-014 were strong negatively correlated with IFN-γ levels, indicating that these bacteria can be tolerated, with IFN-γ responses benefiting for those genera of bacteria. The anti−inflammation activity could be explained in the study regarding to the administration of Norank_f_Lachnospiraceae in gut, and was related to the IL−10 cytokine (Olszak et al., 2014) . The results showed the change of microbiota in the cecum was related to intestinal inflammation (directly or indirectly, positively or negatively).

In conclusion, the present study revealed that obvious changes were found in the microbial community of chicks infected with PDCoV HNZK-02. PDCoV-inoculated chicks showed a decrease in Eisenbergiella and Anaerotruncus populations at 5 dpi and an increase in pathogenetic Alistipes at 17 dpi with a reduction in SCFAs production. The changes in the microbial community may be related to changes in inflammatory cytokine expression. These results provide new insights into the pathophysiology of chicks infected with PDCoV.

The obtained raw sequencing reads have been added to the NCBI SRA page (Accession Number: PRJNA558384). Other data generated during the study are included in this article.

The experimental procedures were approved by the Laboratory Animals Ethics Committee at Hennan Qixiang Biological Technology Co. Ltd. (HNQX-2018-06) , and all husbandry practices and euthanasia were performed with full consideration of animal welfare.

HL and HZ performed the experiments. HL, FZ, and SW performed statistical analyses and wrote the manuscript. Z-xW and Z-yW designed the study and revised the manuscript. All authors read and approved the final manuscript. 

Current findings on gut microbiota mediated immune modulation against viral diseases in chicken

Commensal bacteria calibrate the activation threshold of innate antiviral immunity

Economically important non-oncogenic immunosuppressive viral diseases of chicken-current status

Porcine deltacoronavirus infection and transmission in poultry, United States

TNF-mediated inflammatory disease

Emerging pathogenic links between microbiota and the gut-lung axis

Role of the intestinal microbiota in the immunomodulation of influenza virus infection

Isolation and characterization of porcine epidemic diarrhea viruses associated with the 2013 disease outbreak among swine in the United States

The control and consequences of bacterial fermentation in the human colon

Influenza virus affects intestinal microbiota and secondary salmonella infection in the gut through type I interferons

A gnotobiotic mouse model demonstrates that dietary fiber protects against colorectal tumorigenesis in a microbiota-and butyrate-dependent manner

Necrotic enteritis in broiler chickens: the role of tight junctions and mucosal immune responses in alleviating the effect of the disease

Molecular-phylogenetic characterization of microbial community imbalances in human inflammatory bowel diseases

Experimental infection of gnotobiotic pigs with the cell-culture-adapted porcine deltacoronavirus strain OH-FD22

Eimeria maxima microneme protein 2 delivered as DNA vaccine and recombinant protein induces immunity against experimental homogenous challenge

Differences in the intestinal microbiota between uninfected piglets and piglets infected with porcine epidemic diarrhea virus

Pathogenicity of 2 porcine deltacoronavirus strains in gnotobiotic pigs

Calves are susceptible to infection with the newly emerged porcine deltacoronavirus, but not with the swine enteric alphacoronavirus, porcine epidemic diarrhea virus

Susceptibility of porcine IPEC-J2 intestinal epithelial cells to infection with porcine deltacoronavirus (PDCoV) and serum cytokine responses of gnotobiotic pigs to acute infection with IPEC-J2 cell culture-passaged PDCoV

Impacts of biodiversity on the emergence and transmission of infectious diseases

Porcine deltacoronavirus causes diarrhea in various ages of field-infected pigs in China

Avian influenza virus subtype H9N2 affects intestinal microbiota, barrier structure injury, and inflammatory intestinal disease in the chicken ileum

Infectious bursal disease virus infection leads to changes in the gut associated-lymphoid tissue and the microbiota composition

ERIC-PCR Based Fingerprint to Analyze the Structural Features of Intestinal Microbial Communities of Chicks Infected with Salmonella gallinarum, Newcastle Disease Virus and the Healthy Chicks

Susceptibility of chickens to porcine deltacoronavirus infection

Avian gut-associated lymphoid tissues and intestinal immune responses to Eimeria parasites

Comparison of the expression of cytokine genes in the bursal tissues of the chickens following challenge with infectious bursal disease viruses of varying virulence

Porcine epidemic diarrhea virus infection induced the unbalance of gut microbiota in piglets

Analysis of relative gene expression data using real-time quantitative PCR and the 2-CT method

Multi-omics of the gut microbial ecosystem in inflammatory bowel diseases

The gut microbiota, bacterial metabolites and colorectal cancer

Fermentation in the human large intestine: its physiologic consequences and the potential contribution of prebiotics

FLASH: fast length adjustment of short reads to improve genome assemblies

Gastrointestinal tract morphological alteration by unpleasant physical treatment and modulating role of Lactobacillus in broilers

Design and experimental application of a novel non-degenerate universal primer set that amplifies prokaryotic 16S rRNA genes with a low possibility to amplify eukaryotic rRNA genes

Nutritional Requirements of Poultry. 9th Rev

Protective mucosal immunity mediated by epithelial CD1d and IL-10

Microbial diversity and community composition of caecal microbiota in commercial and indigenous Indian chickens determined using 16s rDNA amplicon sequencing

Cytokine profiles, apoptosis and pathology of experimental Pasteurella multocida serotype A1 infection in mice

Metagenomic biomarker discovery and explanation

Altered gut microbiota profiles in sows and neonatal piglets associated with porcine epidemic diarrhea virus infection

Evaluation of the effect of four fibers on laxation, gastrointestinal tolerance and serum markers in healthy humans

Changes in cecal microbiota community of suckling piglets infected with porcine epidemic diarrhea virus

Detection and genetic characterization of deltacoronavirus in pigs

Innate immunity and intestinal microbiota in the development of Type 1 diabetes

The dynamics of gut-associated microbial communities during inflammation

Discovery of seven novel Mammalian and avian coronaviruses in the genus deltacoronavirus supports bat coronaviruses as the gene source of alphacoronavirus and betacoronavirus and avian coronaviruses as the gene source of gammacoronavirus and deltacoronavirus

Gut microbiota-mediated protection against influenza virus subtype H9N2 in chickens is associated with modulation of the innate responses

Porcine deltacoronavirus: overview of infection dynamics, diagnostic methods, prevalence and genetic evolution

The effect of exposure to high altitude and low oxygen on intestinal microbial communities in mice

The authors thank Dr. Hui Hu from Henan Agricultural University for designing the study and revising the manuscript.