Impact Factor 2.245 | CiteScore 2.6
More on impact ›

Original Research ARTICLE

Front. Vet. Sci., 11 August 2020 | https://doi.org/10.3389/fvets.2020.00452

Transcriptomic Analysis of Chicken Lungs Infected With Avian and Bovine Pasteurella multocida Serotype A

  • 1Chongqing Key Laboratory of Forage & Herbivore, College of Animal Science and Technology, Southwest University, Chongqing, China
  • 2The College of Life Sciences, Sichuan University, Chengdu, China
  • 3College of Veterinary Medicine, Kansas State University, Manhattan, KS, United States

Pasteurella multocida (P. multocida) is a common animal pathogen responsible for many animal diseases. Strains from different hosts exhibit disparate degrees of effect in other species. Here, we characterize an avian P. multocida serogroup A strain (PmQ) showing high lethality to chickens and a bovine P. multocida serogroup A strain (PmCQ2) with no lethality to chickens. We used RNA-seq to profile the transcriptomes of chicken lungs infected with PmQ and PmCQ2. A total of 1,649 differentially expressed genes (DEGs) due to PmQ infection (831 upregulated genes and 818 downregulated genes) and 1427 DEGs (633 upregulated genes and 794 downregulated genes) due to PmCQ2 infection were identified. Functional analysis of these DEGs demonstrated that the TNF signaling pathway, the toll-like receptor signaling pathway, complement and coagulation cascades, and cytokine–cytokine receptor interaction were both enriched in PmQ and PmCQ2 infection. STAT and apoptosis signaling pathways were uniquely enriched by PmQ infection, and the NOD-like receptor signaling pathway was enriched only by PmCQ2 infection. Cell-type enrichment analysis of the transcriptomes showed that immune cells, including macrophages and granulocytes, were enriched in both infection groups. Collectively, our study profiled the transcriptomic response of chicken lungs infected with P. multocida and provided valuable information to understand the chicken responses to P. multocida infection.

Introduction

Pasteurella multocida (P. multocida) is a Gram-negative bacteria, first characterized by Louis Pasteur as the causative agent of fowl cholera in 1881 (1). P. multocida causes many animal diseases, such as fowl cholera, swine atrophic rhinitis, rabbit septicemia, and bovine pneumonia (24). Fowl cholera, caused by avian P. multocida, is a highly contagious disease of various domestic and wild bird species, resulting in great economic losses worldwide (5). In the case of peracute or acute infection, the lung lesions dominated by hemorrhages can be characteristic of fowl cholera (6). Antibiotics are still the main treatment for pasteurellosis due to the lack of effective multi-serotype vaccines (7). In some cases, P. multocida can cause human infection via animal bites and scratches (8, 9). Thus, further study on the mechanism of its pathogenesis is still needed.

Pasteurella multocida can be classified into five serotypes, A, B, D, E, and F, according to the specificity of capsular antigens (10). Different types of P. multocida strains often associated with different diseases depend on the infected host; for example, avian cholera is usually caused by P. multocida of serotype A and F, and swine atrophic rhinitis usually results from serotype D (2). Moreover, P. multocida isolated from certain hosts shows differing pathogenicity to other animals (11, 12). A P. multocida B:2 strain can kill cattle and buffaloes at a low dose but has no affect on chickens even at very high doses (13).

In a previous study, we isolated two P. multocida strains of serotype A, PmQ from tissues of a dead duck and PmCQ2 from pneumonic lungs of a beef cow (14, 15). Both strains showed high virulence to mice, but only PmQ was virulent in chickens and ducks. Here, we have used RNA-seq to investigate the chicken lung transcriptional response to these two P. multocida strains. This study provides comprehensive insights into the chicken immune response to P. multocida and demonstrates that different host-isolated strains show different virulence and may induce diverse immune responses.

Materials and Methods

Bacteria Strains and Culture Conditions

PmQ (GenBank accession No. CP033597) used in this study was isolated from tissues of a dead duck, and PmCQ2 (GenBank accession No. CP033599) was isolated from pneumonic lungs of a beef cow (16). Colonies were picked from Martin's broth agar plate, cultured overnight in Martin's broth aerobically at 37°C with shaking at 220 rpm, and counted by plating on Martin agar plates.

Experimental Animals and Ethics Statement

A total of 150 commercial chickens (Ross 308, 7 days old), obtained from a commercial supplier (JinGang Village Farm, Chongqing), were used in our experiments. All animal experiments were carried out with the approval of the Animal Ethics and Research Committee of Southwest University (permit No. 11-1025), Chongqing, China.

Chicken Challenge

Chickens were kept in wire cages separately in groups (n = 10 per group) with free access to food and water. Chickens of experimental groups were intraperitoneally infected with PmQ (PmQ-P) at a dosage of 1 × 103 colony-forming units (CFU) in 200 μL saline and PmCQ2 (PmCQ2-P) at a dosage of 2 × 109 CFU in 200 μL saline and then monitored for 7 days. The control group (P-C) was treated in parallel with sterile saline of equal volume.

Chickens infected with PmQ were sacrificed at 16 h post-infection (hpi) and lungs were collected. Chickens infected with PmCQ2 were sacrificed at 4, 8, 16, 24, and 48 hpi. Bacterial burdens were determined by homogenizing tissues and diluting ten-fold serially in saline. The different dilutions were plated in triplicate on Martin's broth agar and were incubated at 37°C for 24 h.

Pathological Examinations

For histopathological examination, the lungs of chickens were collected and immediately fixed in 4% paraformaldehyde for 24 h, dehydrated in graded ethanol, and then embedded in paraffin wax. The tissues were sliced at 3 μm thickness and then stained with hematoxylin and eosin (H&E).

RNA Isolation, Library Construction, and Illumina Deep Sequencing

About 50–100 mg of each lung were homogenized in liquid nitrogen. Total RNA was extracted using Trizol reagent (Invitrogen Life Technologies, USA) following the manufacturer's protocol. RNA integrity was confirmed by agarose gel electrophoresis and quantified by NanoDrop (NanoDrop 2000, Thermo Scientific). Following that, the Ribo-Zero rRNA Removal Kit (Illumina, San Diego, CA, USA) was used to remove rRNA. Then, 1 μg of each RNA sample was used to construct a cDNA library according to the manufacturer's recommendations. The cDNA library size was analyzed by using an Agilent 2100 Bioanalyzer and determining the effective concentration using qPCR (StepOnePlus Real-Time PCR Systems, Thermo Scientific). Finally, samples were sequenced on the Illumina sequencing platform (HiSeq 4000), and 150 bp paired-end reads were generated.

Reads Filtering and Mapping

Read adaptors were removed by Cutadapt (Version 1.2.1), and reads with an average base quality lower than Q20 were discarded. Reads were mapped to the genome reference GCF_000002315.4, Gallus_gallus-5.0 from the Ensembl database (ftp://ftp.ebi.ac.uk/pub/databases/ena/wgs/public/aa/AADN05.fasta.gz). A genome index used Bowtie2, and then the clean reads were mapped to the genome by using Tophat2 (http://tophat.cbcb.umd.edu/). Mapping between reads and genomes that are within two mismatches were regarded as successful.

Qualification and Identification of Differently Expression Genes (DEGs)

HTSeq 0.6.1p2 (http://www-huber.embl.de/users/anders/HTSeq) were used to obtain the read count of genes. To eliminate the effects of gene length and sample differences, reads per kilo bases per million (RPKM) was used to qualify the gene expression. Differential expression analysis was performed using estimateSizeFactors in DESeq (version 1.18.0) R package, and p-value and fold-change value were identified using the nbinomTest in DESeq with the R package. P < 0.05 and |log2foldchange| >1 were set as the threshold for DEG identification. Hierarchical cluster analysis was used to identify DEGs with certain patterns of expression under three different challenge groups using Heatmap based on R.

Functional Annotation of DEGs

All DEGs were mapped to the gene ontology (GO) and KEGG ontology (KO) databases to obtain respective GO terms and KO classification and then enrichment analysis of the DEGs was performed using R based on hypergeometric distribution in the context of the entire reference genome. Gene clusters were subject to network analysis by using STRING (version 10.0) (17). Target genes were mapped to the STRING database to gain the relationship between each other and then drawn by Cytoscape (www.cytoscape.org).

To determine lung cell type enrichment, the transcript identifiers of the DEGs for each of the three major groups were converted into their associated gene names using BioMart software (http://uswest.ensembl.org/biomart/martview/). The three DEG lists were then input into the CTen database (18) to obtain the enriched cell types based on a highly expressed, cell-specific gene database. One-sided Fisher's exact test was used for enrichment analysis, and the CTen enrichment score was returned from the –log10 of the Benjamini-Hochberg (BH) adjusted P-values. The CTen database is originally based on human and mouse data and has previously been used for chickens (19, 20).

All sequence data discussed in this study have been deposited to NCBI's Sequence Read Archive (SRA) database, and the BioProject number is PRJNA597560.

Quantitative Real-Time PCR (qRT-PCR)

Total RNA of the lungs was extracted, and then cDNA was synthesized using FastKing gDNA Dispelling RT SuperMix (Tiangen, Beijing, China). Thirteen immune response–related genes, including Il6, Il1β, Il22, Il4i1, S100A9, Socs1, Socs3, AVD, AvBD1, AvBD2, AvBD4, TLR2A, and TLR4, were selected to confirm RNAseq results. Criteria for gene selection were based on immune response function and significance in the RNAseq study. Beta-actin was used as a reference gene to normalize transcript levels of the target gene. Specific primers were designed according to the reference sequences in NCBI with Primer-BLAST under the criteria: (a) PCR product size 70–200 bp, (b) Melting temperatures 60 ± 2°C, and (c) primer must span an exon–exon junction (21). Primer sequences are listed in Table 1. The relative expression levels of genes are shown as a ratio of the target gene to the reference gene using the formula 2−(ΔΔCt) (22).

TABLE 1
www.frontiersin.org

Table 1. Primers used for validation of RNA-seq data.

Statistical Analysis

All data are displayed as means ± standard error of mean (SEM). GraphPad Prism 6.0 software was applied for all statistical analysis. Survival rates of chicks were evaluated using Kaplan-Meier analysis. All the other data were evaluated using unpaired, two-tailed Student's t-test. Significant differences were considered at P < 0.05. *P < 0.05, **P < 0.01, and ***P < 0.001.

Results

Mortality and Histopathology of Infected Chickens

Chickens infected with PmCQ2 and saline all survived, and PmQ-infected chickens all died rapidly within 24 h (Figure 1A). The bacterial burden of PmCQ2 in chicken lungs was approximately 105 CFU/mg at 4 and 8 hpi, maintained a high density at 16 hpi, then decreased to 101 CFU/mg at 24 hpi (Figure 1B). The bacterial load of PmQ in chicken lungs was about 104 CFU/mg, on average, at 16 hpi (Figure 1C). Although chickens infected with PmQ had many inflammatory cell infiltrates, PmCQ2 infected chickens had few inflammatory cells detected (Figure 1D). PmCQ2-infected chickens all recovered with 7 days with no clinical signs of infection.

FIGURE 1
www.frontiersin.org

Figure 1. Chicken model infected with PmCQ2 and PmQ. (A) Survival curve of chickens challenged intraperitoneally with PmQ (1 × 103 CFU), PmCQ2 (2 × 109 CFU), and control (saline); (B) Bacterial burden of chicken lungs infected with PmCQ2 at 4, 8, 16, and 24 h; (C) Bacterial burden comparison of chicken lungs infected PmQ and PmCQ2 at 16 h; (D) The pulmonary histopathological examination of chicken at 16 h, stained with H&E (magnification 400X with a bar of 100 μm labeled). **p < 0.001, ****p < 0.00001, and ns, no significance.

Evaluation of Transcriptomic Sequencing Quality

A total of nine samples were generated from three biological replicates under three treatments; 44,374,891, 95,465,238, and 99,426,423 average clean reads were obtained from control and two infection group samples, respectively. The percentages mapped to the reference genome were more than 80% in all samples. A total of 19,119 genes were annotated in the Ensembl database for chicken genome. The average number of detected genes for individuals are 14,363, accounting for more than 75% of the reference genes. The percentage of uniquely mapped reads in all these samples was above 96%. Sequence data statistics are shown in Table 2.

TABLE 2
www.frontiersin.org

Table 2. Statistic of RNA-seq output.

Identification of DEGs in the Lungs After P. multocida Challenge

Based on the criteria: |log2foldchange| >1 and P < 0.05, there were 1,649 DEGs in P-C vs. PmQ-P, of which 831 were upregulated and 818 were downregulated; 1,427 DEGs (633 upregulated genes and 794 downregulated genes) in P-C vs. PmCQ2-P; and 318 DEGs with 232 upregulated genes and 86 downregulated genes in PmCQ2-P vs. PmQ-P (Figure 2A). The description of total DEGs is displayed in Supplementary Table 1. A total of 899 DEGs were common in P-C vs. PmQ-P and P-C vs. PmCQ2-P (Figure 2B), of which 31 DEGs were also found in PmQ infection vs. PmCQ2. The full gene list is displayed in Supplementary Table 2. Unsupervised hierarchical clustering analysis of global gene expression profiles revealed that each of the three biological triplicates clustered together (Figure 2C).

FIGURE 2
www.frontiersin.org

Figure 2. Differentially expressed genes (DEGs) results. (A) Numbers of up-/downregulated DEGs in three contrasts (P-C vs. PmQ-P, P-C vs. PmCQ2-P, and PmCQ2-P vs. PmQ-P). (B) Venn diagram of three contrasts. (C) The clustering of DEGs in a heat map.

Experimental Validation of DEGs

To validate the results of RNA-seq, 13 immune response–related genes, including Il6, Il1β, Il22, Il4i1, S100A9, Socs1, Socs3, AVD, AvBD1, AvBD2, AvBD4, TLR2A, and TLR4, were chosen for qPCR detection. The gene expression patterns detected by qPCR were similar to those obtained from the RNA-seq results (Table 3).

TABLE 3
www.frontiersin.org

Table 3. qRT-PCR validation.

Functional Analysis of DEGs

GO enrichment revealed that there were 34 GO terms enriched (of 3,864 annotated) in P-C vs. PmQ-P; 21 GO terms enriched (of 3,448 annotated) in P-C vs. PmCQ2-P (Table 4). The detailed information is listed in Supplementary Table 3. In both P-C vs. PmQ-P and P-C vs. PmCQ2-P, “defense response” and “immune system process” categories were predominant. Although categories including “inflammatory response,” “response to cytokine,” and “response to chemical” are only enriched in P-C vs. PmQ-P, the enriched GO terms in P-C vs. PmCQ2-P are more often associated with leukocytes and granulocytes, such as “leukocyte migration,” “granulocyte chemotaxis,” “leukocyte chemotaxis,” and “granulocyte migration.”

TABLE 4
www.frontiersin.org

Table 4. Enriched GO terms analysis.

KEGG pathway enrichment showed that 33 pathways were enriched in P-C vs. PmQ-P, of which the top 20 enriched pathways are displayed (Figure 3). Similarly, 30 pathways were enriched in P-C vs. PmCQ2-P (Figure 4). Many of these pathways are related to the immune system; signal transduction; the digestive system; and metabolism of lipids, amino acids, and xenobiotics. Both of the two contrasts clustered several conserved pathways, including “TNF signaling pathway,” “Toll-like receptor signaling pathway,” “Complement and coagulation cascades,” “Cytokine–cytokine receptor interaction,” “NF-kappa B signaling pathway,” “Fc gamma R-mediated phagocytosis,” and “Chemokine signaling pathway,” indicating that these pathways are related to P. multocida infection (Supplementary Table 4).

FIGURE 3
www.frontiersin.org

Figure 3. Top 20 enriched KEGG pathway in contrasts, P-C vs. PmQ-P. The numbers of DEGs in each pathway are counted, and rich factors and P-values are displayed.

FIGURE 4
www.frontiersin.org

Figure 4. Top 20 enriched KEGG pathway in contrasts, P-C vs. PmCQ2-P. The numbers of DEGs in each pathway are counted, and rich factors and P-values are displayed.

Furthermore, cell type enrichment analysis based on upregulated DEGs showed that different types of macrophages are the most enriched immune cell type in both PmQ and PmCQ2 infected chicken vs. uninfected chicken (Figures 5A,B). Other immune related cell types, such as microglia, osteoclasts, granulocytes, lymph nodes, CD8a+ Dend. cells myeloid, bone marrow, and mast cells, were all significantly enriched (Supplementary Table 3). Specifically, comparing the PmCQ2 infected group with PmQ infection, the most enriched cell type was also focused on macrophages (Figure 5C), which suggests that macrophage may be the key immune cell upon P. multocida infection in chicken. Moreover, the expression level of macrophage related genes was displayed with a heat map (Figure 6). The clustered gene group of each gene showed that TLR4, TLR2A, Macro, Ccr2, Ccr5, MRC1L-B, and MRC1L-D were upregulated in both PmCQ2 and PmQ infected groups. Moreover, genes such as Ifr9, Ifr8, Ifr7, CSF2RA, STAT1, and STAT2 were expressed to a greater extent in the PmQ infection group than PmCQ2, and Il22, Il6, Nos2, Ccl4, Ccl5, Ifr1, Socs1, and Socs3 were increased only in PmQ infected chickens rather than PmCQ2.

FIGURE 5
www.frontiersin.org

Figure 5. The 20 most enriched cell types. The enrichment scores (labeled inside the circle) is displayed in the –log10 of the Benjamini-Hochberg (BH) adjusted P-values. Cell/tissue types (labeled outside the circle) with an enrichment score (red line) greater than two were considered significant. (A) P-C vs. PmQ-P. (B) P-C vs. PmCQ2-P. (C) PmCQ2-P vs. PmQ-P.

FIGURE 6
www.frontiersin.org

Figure 6. The expression level of macrophage-related genes displayed with heat map. The expression level (RPKM) of each gene were picked from Supplementary Table 7 and clustered by z-score.

Discussion

Pasteurella multocida is a very common animal pathogen and is classified into five serotypes (A, B, D, E, and F) depending on the capsular antigen (10). Among them, serotype A tends to have a wide host range (2, 3). However strains isolated from non-avian hosts generally do not infect poultry, and the underling mechanisms are poorly understood (23). An avian P. multocida of serotype A, PmQ, isolated from tissues of a dead duck shows high virulence to poultry as well as to mammals. A bovine strain of serotype A, PmCQ2, from pneumonic lungs of a beef cow can cause severe pneumonia in mice (14, 16) but is non-lethal to chickens. It has been suggested that some of the differences in host susceptibility to P. multocida infection may be due to differences in the host response expressed in the lung during the early phase of infection (24). In the present study, we infected chickens with PmCQ2 at an extremely high dosage (2 × 109 CFU), via the abdominal cavity to artificially create a chicken “infection model.” However, the PmCQ2 is quickly eliminated by the hosts, and a dramatic drop between 16 and 24 hpi of the bacterial load occurs. Whereas, the PmQ infected chickens died rapidly within 24 hpi. Therefore, transcripts of the chicken lungs at 16 hpi were used for RNA-seq.

The innate immune system is the first line in host defense and protects the host from bacterial infection (25). Previously, we found TLRs and inflammatory cytokines such as IL17, IL6, TNF-α, IFN-γ were significantly upregulated in mouse lungs infected with PmCQ2 (14). Another study also reports that IL6, IL17, and IL22 expression levels were significantly increased in the lungs of chicken challenged with avian P. multocida strain X73, which succumbed to acute infection (26). Consistent with these studies, the transcriptome of chicken lungs infected with PmQ also demonstrate that cytokines, such as IL6 and IL22, were significantly upregulated in PmQ-P vs. P-C, which implies PmQ infection induces a strong inflammatory reaction in the lungs. IL6 is a pleiotropic cytokine and can be both pro- and anti-inflammatory depending on the context (27). Impaired IL6 function causes enhanced susceptibility of mice to infection of Listeria monocytogenes (28), but excessive and sustained production of IL6 can be life-threatening, including systemic inflammatory response syndrome (SIRS) and chronic inflammatory diseases (29). IL22 is a member of the IL10 family and targets epithelial cells to induce various innate immune responses, which is essential for host defense against pathogen invasion (30). It can not only protect the tissue integrity and maintain mucosal homeostasis, but also act as a pro-inflammatory cytokine capable of amplifying the inflammatory reaction, and overproduction of IL22 may result in tissue damage (31). IL22 can mediate heme-associated iron scavenging away from pathogens during systematic infection (32, 33).

In PmCQ2 infected chickens, the expression level of IL22 and IL6 were unaffected. Similarly, chickens challenged with capsular-deficient strain 775 with low virulence also fail to increase IL6 and IL22 expression (26). Further verification using q-PCR indicated that IL6 and IL22 expression levels were unchanged during PmCQ2 infection (Supplementary Figure 1), which implies that PmCQ2 causes a lower inflammatory response. CCL4, CCL5, and CCL28 were upregulated in PmQ and PmCQ2 infected chickens, but the expression level was increased to a high level in PmQ-P. CCL4, also known as macrophage inflammatory protein-1β (MIP-1β), is a chemoattractant for natural killer cells, monocytes, and a variety of other immune cells (34); CCL5 is chemotactic for T cells, eosinophils, and basophils and plays an active role in recruiting leukocytes into inflammatory sites (35); and CCL28 shows both antimicrobial activity and functions as a chemoattractant for T and B lymphocytes and migration of eosinophils. Together, these chemokines are ascribed to more immune cell infiltration, such as lymphocytes and eosinophils, in the PmQ infected chicken lungs. Whereas, the expression level of CCL19 and CCL26 were only increased in PmQ-P vs. P-C. CCL19 binding to its receptor CCR7 attracts dendritic cells to lymph nodes (36), and CCL26, also called macrophage inflammatory protein 4-alpha (MIP-4-alpha), can drive eosinophils and basophil migration by binding to CCR3 (37). In sum, these data suggest a different immune response of the host infected with the two P. multocida strains.

In this study, a number of enrichment GO function categories were identified based on the DEGs of PmQ-P vs. P-C and PmCQ2-P vs. P-C. The protein–protein interactions based on DEGs annotated in immune system response of two infection groups (Figures 7A,B) indicated a different immune response interaction network. In PmQ infected chicken, IL6 and PTPRC (CD45) are in the center of the net. Upregulation of IL6 combined with STAT1 and STAT2/3, attributed to the key points of the Jak-STAT pathway, demonstrated that the STAT pathway was activated during P. multocida infection, which is consistent with our previous study (14). The activation of the STAT pathway can then trigger a signaling cascade required for immune response and inflammatory reaction (38), resulting in production of IFN-γ, contributing to bacterial clearance (39, 40). CD45, expressed on leukocytes, is a crucial player in the activation of T-cell receptor signaling by controlling the activation of the Src family protein-tyrosine kinases Lck and Fyn. CD45 deficiency can result in T- and B-lymphocyte dysfunction in the form of severe combined immune deficiency (41). Some immunomodulatory proteins, such as UL11 of human cytomegalovirus, can interact with CD45, resulting in function disruption of T-cells (42).

FIGURE 7
www.frontiersin.org

Figure 7. STRING network analysis on DEGs annotated in GO: immune system response. (A) P-C vs. PmQ-P. (B) P-C vs. PmCQ2-P. The nod represents proteins, bigger means more interaction with others. The red color means upregulated DEGs, and the green means downregulated DEGs.

Different from P. multocida toxin (PMT) (43, 44), lipopolysaccharides isolated from P. multocida serotype A and B stains can induce apoptosis of bovine leukocytes, including lymphocytes, neutrophils, monocytes, and macrophages, resulting in tissue damage (45, 46). In the current study, genes (FAS, APAF1, TNFRSF1A, CASP7, and CASP8) related to apoptosis were significantly upregulated in PmQ infected chickens but not in chickens infected with PmCQ2. In contrast to our previous study, the apoptosis pathway is enriched in PmCQ2 infected mice with a deadly outcome (14). Together, these results hint that the PmQ infected chicken attempt to cause cell death of the host, leading to tissue lesions of the lungs.

The NOD-like receptors (NLRs) are a specialized group of intracellular innate immune receptors that mediate IL1β production against bacterial infection (47). In the present study, the NOD-like receptor signaling pathway was enriched alone in the PmCQ2 infected chicken, which is consistent with the previous study of PmCQ2 infected mice (14). NLRC5, which can act as a negative modulator of inflammatory pathways, is critical for LPS-induced IL-10 production in macrophages (48), which upregulated in PmCQ2 infection, which suggests it might contribute to the lower inflammatory reaction in PmCQ2 infection. Our previous study revealed the NLRP3 inflammasome of mouse peritoneal macrophages infected with PmCQ2 played a role in mediating IL1β against P. multocida infection (49). In line with the study, IL1β was also significantly increased in PmCQ2 infected chicken, which might responsible for PmCQ2 clearance.

Macrophages play a crucial role in sensing P. multocida and secretion of inflammatory factors (49). We found that L-serine supplementation can decrease the macrophage-mediated inflammatory reaction in PmCQ2 infected mice and reduce the lung damage (50). In the current study, cell type enrichment analysis based on DEGs showed that macrophages were the most significantly enriched cell types, suggesting that macrophages were also activated in chickens infected with two different host-isolated P. multocida strains (Figure 5). Moreover, when analyzing the DEGs between PmCQ2 and PmQ infected chicken, macrophage was still the most enriched cell type. However, the different outcomes of chicken challenged with PmCQ2 and PmQ hint that macrophage may play different roles. For example, NOS2, a key marker, and PTX3 and IRF5, two important regulators of M1 macrophage activation (5153), were significantly overexpressed in PmQ infected chicken, tending to activate to more pro-inflammatory activation. Whereas, in the PmCQ2 infected chicken lung, M2 markers (MRC1LD, MRC1LB, and MACRO) and M1 markers (CD86 and CD80) were upregulated, implying a mixed activation of both M1 and M2 macrophage occurred. In Marek's disease virus (MDV) infection, macrophages derived from differentially sensitive chicken lines undergo different immune strategies, and the macrophages from resistant chicken lines activate different biological pathways and suppression of oncogenic potential (54). Macrophages can be activated into different states depending upon the microenvironment and stimulus encountered (55, 56), and the activation styles are often related to its morphology (57).

In addition to the immune response, metabolism was also affected upon P. multocida infection. For example, GAD1, a key enzyme that converts L-glutamate to GABA, used for succinic semialdehyde synthesis and eventually succinate, which is the main route of M1 macrophage and fundamental for stabilization of HIFα (58, 59), was significantly upregulated in PmQ infected chicken while GLUL (also known as glutamine synthetase, GS), highly expressed in M2 macrophages and important to M2 phenotype (60), is increased in PmCQ2 infected chicken. However, PmCQ2 infection also caused reduction of all NADH dehydrogenase subunits (ND1~ND6), implying the oxidative phosphorylation of chicken was impaired, which was unfavorable to M2 macrophage activation (61). A recent study reported that arachidonic acid can bind to myeloid differentiation factor-2 (MD2) and prevent MD2/TLR4 signaling activation, thus, it inhibits inflammatory response (62). Ten genes associated with arachidonic acid metabolism are upregulated in PmQ infected chickens, resulting in overproduction of a number of metabolites, for example, PGD2 (Prostaglandin D2), which can bind to its receptor DP2 and amplify inflammatory response in asthma (63). Overall, these results suggest macrophages might undergo different polarization and contribute to opposite outcomes of chickens.

In conclusion, different host-isolated serotype A P. multocida strains, bovine PmCQ2 and avian PmQ, exhibited opposite pathogenicity to birds. Additionally, this study also revealed several crucial pathways related to anti-infection and metabolisms were enriched and that macrophages might play a key role in P. multocida infection.

Data Availability Statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found at: https://www.ncbi.nlm.nih.gov/, PRJNA597560.

Ethics Statement

This animal study was reviewed and approved by Animal Ethics and Research Committee of Southwest University.

Author Contributions

YP and NL conceived the experiments. PL conducted the experiments and wrote the paper. FH, CW, and GZ contributed to the experiments and data analysis. PH, YP, and NL supervised the project and revised the paper. All authors discussed the results and approved the final manuscript.

Funding

This work was supported by the earmarked fund of the China Agriculture Research System (Beef/Yak Cattle, CARS-37), Chongqing Science & Technology Commission (cstc2017shms-zdyfx0036 and cstc2017jcyjAX0288). The funding organizations had no role in study design, data collection or analysis, decision to publish, or preparation.

Conflict of Interest

The sequencing service was provided by Personal Biotechnology Co., Ltd., Shanghai, China.

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We thank the sequencing service provided by Personal Biotechnology Co., Ltd., Shanghai, China.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fvets.2020.00452/full#supplementary-material

References

1. Kubatzky KF. Pasteurella multocida and immune cells. Curr Top Microbiol Immunol. (2012) 361:53–72. doi: 10.1007/82_2012_204

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Wilkie IW, Harper M, Boyce JD, Adler B. Pasteurella multocida: diseases and pathogenesis. Curr Top Microbiol Immunol. (2012) 361:1–22. doi: 10.1007/82_2012_216

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Wilson BA, Ho M. Pasteurella multocida: from zoonosis to cellular microbiology. Clin Microbiol Rev. (2013) 26:631–55. doi: 10.1128/CMR.00024-13

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Peng Z, Wang XR, Zhou R, Chen HC, Wilson BA, Wu B. Pasteurella multocida: Genotypes and Genomics. Microbiol Mol Biol Rev. (2019) 83:e00014–00019. doi: 10.1128/MMBR.00014-19

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Sellyei B, Thuma A, Volokhov D, Varga Z. Comparative analysis of Pasteurella multocida isolates from acute and chronic fowl cholera cases in Hungary during the period 2005 through 2010. Avian Dis. (2017) 61:457–65. doi: 10.1637/11674-051817-Reg.1

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Christensen JP, Bisgaard M. Fowl cholera. Rev Sci Tech. (2000) 19:626–37. doi: 10.20506/rst.19.2.1236

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Liu Q, Hu Y, Li P, Kong Q. Identification of fur in Pasteurella multocida and the potential of its mutant as an attenuated live vaccine. Front Vet Sci. (2019) 6:5. doi: 10.3389/fvets.2019.00005

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Honnorat E, Seng P, Savini H, Pinelli PO, Simon F, Stein A. Prosthetic joint infection caused by Pasteurella multocida: a case series and review of literature. BMC Infect Dis. (2016) 16:435. doi: 10.1186/s12879-016-1763-0

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Ryan JM, Feder HM Jr. Dog licks baby. Baby gets Pasteurella multocida meningitis. Lancet. (2019) 393:e41. doi: 10.1016/S0140-6736(19)30953-5

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Carter GR. Studies on Pasteurella multocida. I. A hemagglutination test for the identification of serological types. Am J Vet Res. (1955) 16:481–484.

PubMed Abstract | Google Scholar

11. Okay S, Kurt Kizildogan A. Comparative genome analysis of five Pasteurella multocida strains to decipher the diversification in pathogenicity and host specialization. Gene. (2015) 567:58–72. doi: 10.1016/j.gene.2015.04.063

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Peng Z, Liang W, Wang F, Xu Z, Xie Z, Lian Z, et al. Genetic and phylogenetic characteristics of Pasteurella multocida isolates from different host species. Front Microbiol. (2018) 9:1408. doi: 10.3389/fmicb.2018.01408

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Aktories K, Orth J, Adler B. Pasteurella multocida. Berlin; Heidelberg: Verlag; Springer (2012). p. 5

Google Scholar

14. Wu C, Qin X, Li P, Pan T, Ren W, Li N, et al. Transcriptomic analysis on responses of murine lungs to Pasteurella multocida infection. Front Cell Infect Microbiol. (2017) 7:251. doi: 10.3389/fcimb.2017.00251

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Du H, Fang R, Pan T, Li T, Li N, He Q, et al. Comparative genomics analysis of two different virulent bovine Pasteurella multocida isolates. Int J Genomics. (2016) 2016:4512493. doi: 10.1155/2016/4512493

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Li N, Long Q, Du H, Zhang J, Pan T, Wu C, et al. High and low-virulent bovine Pasteurella multocida capsular type A isolates exhibit different virulence gene expression patterns in vitro and in vivo. Vet Microbiol. (2016) 196:44–49. doi: 10.1016/j.vetmic.2016.10.017

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Szklarczyk D, Franceschini A, Wyder S, Forslund K, Heller D, Huerta-Cepas J, et al. STRING v10: protein-protein interaction networks, integrated over the tree of life. Nucleic Acids Res. (2015) 43:D447–52. doi: 10.1093/nar/gku1003

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Shoemaker JE, Lopes TJS, Ghosh S, Matsuoka Y, Kawaoka Y, Kitano H. CTen: a web-based platform for identifying enriched cell types from heterogeneous microarray data. BMC Genomics. (2012) 13:460. doi: 10.1186/1471-2164-13-460

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Deist MS, Gallardo RA, Bunn DA, Dekkers JCM, Zhou H, Lamont SJ. Resistant and susceptible chicken lines show distinctive responses to Newcastle disease virus infection in the lung transcriptome. BMC Genomics. (2017) 18:989. doi: 10.1186/s12864-017-4380-4

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Zhang J, Kaiser MG, Deist MS, Gallardo RA, Bunn DA, Kelly TR, et al. Transcriptome analysis in spleen reveals differential regulation of response to Newcastle disease virus in two chicken lines. Sci Rep. (2018) 8:1278. doi: 10.1038/s41598-018-19754-8

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Ye J, Coulouris G, Zaretskaya I, Cutcutache I, Rozen S, Madden TL. Primer-BLAST: a tool to design target-specific primers for polymerase chain reaction. BMC Bioinform. (2012 13:134. doi: 10.1186/1471-2105-13-134

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Lee SI, Jang HJ, Jeon MH, Lee MO, Kim JS, Jeon IS, et al. Transcriptional regulation of cathelicidin genes in chicken bone marrow cells. Poult Sci. (2016) 95:912–19. doi: 10.3382/ps/pev361

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Ahir VB, Roy A, Jhala MK, Bhanderi BB, Mathakiya RA, Bhatt VD, et al. Genome sequence of Pasteurella multocida subsp gallicida Anand1_poultry. J Bacteriol. (2011) 193:5604–4. doi: 10.1128/Jb.05706-11

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Jens P, Christensen, Miki Bojesen A, Bisgaard M. Fowl Cholera. Philadelphia, PA: Elsevier. (2008). p. 149–54.

Google Scholar

25. Mogensen TH. Pathogen recognition and inflammatory signaling in innate immune defenses. Clin Microbiol Rev. (2009) 22:240–73. doi: 10.1128/CMR.00046-08

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Petruzzi B, Dalloul RA, LeRoith T, Evans NP, Pierson FW, Inzana TJ. Biofilm formation and avian immune response following experimental acute and chronic avian cholera due to Pasteurella multocida. Vet Microbiol. (2018) 222:114–23. doi: 10.1016/j.vetmic.2018.07.005

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Hunter CA, Jones SA. IL-6 as a keystone cytokine in health and disease. Nat Immunol. (2015) 16:448–57. doi: 10.1038/ni.3153

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Hoge J, Yan I, Janner N, Schumacher V, Chalaris A, Steinmetz OM, et al. IL-6 controls the innate immune response against listeria monocytogenes via classical IL-6 signaling. J Immunol. (2013) 190:703–11. doi: 10.4049/jimmunol.1201044

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Kang S, Tanaka T, Kishimoto T. Therapeutic uses of anti-interleukin-6 receptor antibody. Int Immunol. (2015) 27:21–9. doi: 10.1093/intimm/dxu081

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Kim S, Faris L, Cox CM, Sumners LH, Jenkins MC, Fetterer RH, et al. Molecular characterization and immunological roles of avian IL-22 and its soluble receptor IL-22 binding protein. Cytokine. (2012) 60:815–27. doi: 10.1016/j.cyto.2012.08.005

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Eidenschenk C, Rutz S, Liesenfeld O, Ouyang W. Role of IL-22 in microbial host defense. Curr Top Microbiol Immunol. (2014) 380:213–36. doi: 10.1007/978-3-662-43492-5_10

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Kugelberg E. Infection: IL-22 controls iron scavenging. Nat Rev Immunol. (2017) 17:146–7. doi: 10.1038/nri.2017.15

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Sakamoto K, Kim YG, Hara H, Kamada N, Caballero-Flores G, Tolosano E, et al. IL-22 controls iron-dependent nutritional immunity against systemic bacterial infections. Sci Immunol. (2017) 2:eaai8371. doi: 10.1126/sciimmunol.aai8371

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Bystry RS, Aluvihare V, Welch KA, Kallikourdis M, Betz AG. B cells and professional APCs recruit regulatory T cells via CCL4. Nat Immunol. (2001) 2:1126–32. doi: 10.1038/ni735

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Maghazachi AA, AlAoukaty A, Schall TJ. CC chemokines induce the generation of killer cells from CD56(+) cells. Eur J Immunol. (1996) 26:315–9. doi: 10.1002/eji.1830260207

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Robbiani DF, Finch RA, Jager D, Muller WA, Sartorelli AC, Randolph GJ. The leukotriene C(4) transporter MRP1 regulates CCL19 (MIP-3beta, ELC)-dependent mobilization of dendritic cells to lymph nodes. Cell. (2000) 103:757–68. doi: 10.1016/s0092-8674(00)00179-3

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Kitaura M, Suzuki N, Imai T, Takagi S, Suzuki R, Nakajima T, et al. Molecular cloning of a novel human CC chemokine (eotaxin-3) that is a functional ligand of CC chemokine receptor 3. J Biol Chem. (1999) 274:27975–80. doi: 10.1074/jbc.274.39.27975

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Kiu H, Nicholson SE. Biology and significance of the JAK/STAT signalling pathways. Growth Factors. (2012) 30:88–106. doi: 10.3109/08977194.2012.660936

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Zeng X, Gu H, Peng LS, Yang Y, Wang N, Shi Y, et al. Transcriptome profiling of lung innate immune responses potentially associated with the pathogenesis of Acinetobacter baumannii acute lethal Pneumonia. Front Immunol. (2020) 11:708. doi: 10.3389/fimmu.2020.00708

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Hu X, Ivashkiv LB. Cross-regulation of signaling pathways by interferon-gamma: implications for immune responses and autoimmune diseases. Immunity. (2009) 31:539–50. doi: 10.1016/j.immuni.2009.09.002

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Rheinlander A, Schraven B, Bommhardt U. CD45 in human physiology and clinical medicine. Immunol Lett. (2018) 196:22–32. doi: 10.1016/j.imlet.2018.01.009

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Gabaev I, Steinbruck L, Pokoyski C, Pich A, Stanton RJ, Schwinzer R, et al. The human cytomegalovirus UL11 protein interacts with the receptor tyrosine phosphatase CD45, resulting in functional paralysis of T cells. PLoS Pathog. (2011) 7:e1002432. doi: 10.1371/journal.ppat.1002432

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Hildebrand D, Heeg K, Kubatzky KF. Pasteurella multocida toxin manipulates T cell differentiation. Front Microbiol. (2015) 6:1273. doi: 10.3389/fmicb.2015.01273

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Preuss I, Hildebrand D, Orth JH, Aktories K, Kubatzky KF. Pasteurella multocida toxin is a potent activator of anti-apoptotic signalling pathways. Cell Microbiol. (2010) 12:1174–185. doi: 10.1111/j.1462-5822.2010.01462.x

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Periasamy S, Praveena PE, Singh N. Effects of Pasteurella multocida lipopolysaccharides on bovine leukocytes. Microb Pathog. (2018) 119:225–32. doi: 10.1016/j.micpath.2018.04.030

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Praveena PE, Periasamy S, Kumar AA, Singh N. Cytokine profiles, apoptosis and pathology of experimental Pasteurella multocida serotype A1 infection in mice. Res Vet Sci. (2010) 89:332–9. doi: 10.1016/j.rvsc.2010.04.012

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Chen G, Shaw MH, Kim YG, Nunez G. NOD-like receptors: role in innate immunity and inflammatory disease. Ann Rev Pathol Mech Dis. (2009) 4:365–98. doi: 10.1146/annurev.pathol.4.110807.092239

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Benko S, Magalhaes JG, Philpott DJ, Girardin SE. NLRC5 limits the activation of inflammatory pathways. J Immunol. (2010) 185:1681–91. doi: 10.4049/jimmunol.0903900

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Fang R, Du H, Lei G, Liu Y, Feng S, Ye C, et al. NLRP3 inflammasome plays an important role in caspase-1 activation and IL-1beta secretion in macrophages infected with Pasteurella multocida. Vet Microbiol. (2019) 231:207–13. doi: 10.1016/j.vetmic.2019.03.019

PubMed Abstract | CrossRef Full Text | Google Scholar

50. He F, Yin Z, Wu C, Xia Y, Wu M, Li P, et al. l-serine lowers the inflammatory responses during Pasteurella multocida infection. Infect Immun. (2019) 87:e00677–19. doi: 10.1128/IAI.00677-19

PubMed Abstract | CrossRef Full Text | Google Scholar

51. Murray PJ. Macrophage polarization. Annu Rev Physiol. (2017) 79:541–66. doi: 10.1146/annurev-physiol-022516-034339

PubMed Abstract | CrossRef Full Text | Google Scholar

52. Ferrer MF, Thomas P, Lopez Ortiz AO, Errasti AE, Charo N, Romanowski V, et al. Junin virus triggers macrophage activation and modulates polarization according to viral strain pathogenicity. Front Immunol. (2019) 10:2499. doi: 10.3389/fimmu.2019.02499

PubMed Abstract | CrossRef Full Text | Google Scholar

53. Wei J, Tang D, Lu C, Yang J, Lu Y, Wang Y, et al. Irf5 deficiency in myeloid cells prevents necrotizing enterocolitis by inhibiting M1 macrophage polarization. Mucosal Immunol. (2019) 12:888–96. doi: 10.1038/s41385-019-0169-x

PubMed Abstract | CrossRef Full Text | Google Scholar

54. Chakraborty P, Kuo R, Vervelde L, Dutia BM, Kaiser P, Smith J. Macrophages from susceptible and resistant chicken lines have different transcriptomes following Marek's disease virus infection. Genes. (2019) 10:74. doi: 10.3390/genes10020074

PubMed Abstract | CrossRef Full Text | Google Scholar

55. Arora S, Dev K, Agarwal B, Das P, Syed MA. Macrophages: their role, activation and polarization in pulmonary diseases. Immunobiology. (2018) 223:383–96. doi: 10.1016/j.imbio.2017.11.001

PubMed Abstract | CrossRef Full Text | Google Scholar

56. Murray PJ, Wynn TA. Protective and pathogenic functions of macrophage subsets. Nat Rev Immunol. (2011) 11:723–37. doi: 10.1038/nri3073

PubMed Abstract | CrossRef Full Text | Google Scholar

57. McWhorter FY, Wang TT, Nguyen P, Chung T, Liu WF. Modulation of macrophage phenotype by cell shape. Proc Natl Acad Sci USA. (2013) 110:17253–58. doi: 10.1073/pnas.1308887110

PubMed Abstract | CrossRef Full Text | Google Scholar

58. Tannahill GM, Curtis AM, Adamik J, Palsson-McDermott EM, McGettrick AF, Goel G, et al. Succinate is an inflammatory signal that induces IL-1 beta through HIF-1 alpha. Nature. (2013) 496:238. doi: 10.1038/nature11986

PubMed Abstract | CrossRef Full Text | Google Scholar

59. Meiser J, Kramer L, Sapcariu SC, Battello N, Ghelfi J, D'Herouel AF, et al. Pro-inflammatory macrophages sustain pyruvate oxidation through pyruvate dehydrogenase for the synthesis of itaconate and to enable cytokine expression. J Biol Chem. (2016) 291:3932–46. doi: 10.1074/jbc.M115.676817

PubMed Abstract | CrossRef Full Text | Google Scholar

60. Palmieri EM, Menga A, Martin-Perez R, Quinto A, Riera-Domingo C, De Tullio G, et al. Pharmacologic or genetic targeting of glutamine synthetase skews macrophages toward an M1-like phenotype and inhibits tumor metastasis. Cell Rep. (2017) 20:1654–66. doi: 10.1016/j.celrep.2017.07.054

PubMed Abstract | CrossRef Full Text | Google Scholar

61. Viola A, Munari F, Sanchez-Rodriguez R, Scolaro T, Castegna A. The metabolic signature of macrophage responses. Front Immunol. (2019) 10:1462. doi: 10.3389/fimmu.2019.01462

PubMed Abstract | CrossRef Full Text | Google Scholar

62. Zhang YL, Chen HJ, Zhang WX, Cai Y, Shan PR, Wu D, et al. Arachidonic acid inhibits inflammatory responses by binding to myeloid differentiation factor-2 (MD2) and preventing MD2/toll-like receptor 4 signaling activation. Biochim Biophys Acta Mol Basis Dis. (2020) 1866:165683. doi: 10.1016/j.bbadis.2020.165683

PubMed Abstract | CrossRef Full Text | Google Scholar

63. Brightling CE, Brusselle G, Altman P. The impact of the prostaglandin D2 receptor 2 and its downstream effects on the pathophysiology of asthma. Allergy. (2020) 75:761–8. doi: 10.1111/all.14001

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: Pasteurella multocida, chicken lung, infection, transcriptome, immune response

Citation: Li P, He F, Wu C, Zhao G, Hardwidge PR, Li N and Peng Y (2020) Transcriptomic Analysis of Chicken Lungs Infected With Avian and Bovine Pasteurella multocida Serotype A. Front. Vet. Sci. 7:452. doi: 10.3389/fvets.2020.00452

Received: 04 April 2020; Accepted: 22 June 2020;
Published: 11 August 2020.

Edited by:

Subhash Verma, Chaudhary Sarwan Kumar Himachal Pradesh Krishi Vishvavidyalaya, India

Reviewed by:

Katharina F. Kubatzky, Heidelberg University, Germany
Md. Aminul Islam, Tohoku University, Japan

Copyright © 2020 Li, He, Wu, Zhao, Hardwidge, Li and Peng. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Nengzhang Li, lich2001020@163.com; Yuanyi Peng, pyy2002@sina.com