key: cord-0794860-okahh9wg
authors: Park, Ga Young; Yu, Hyun Jin; Son, Jee Soo; Park, Sang Joon; Cha, Hee-Jae; Song, Kyoung Seob
title: Specific bacteriophage of Bordetella bronchiseptica regulates B. bronchiseptica-induced microRNA expression profiles to decrease inflammation in swine nasal turbinate cells
date: 2020-02-07
journal: Genes Genomics
DOI: 10.1007/s13258-019-00906-7
sha: 4462b7f646558c5a0eee6f02a45313e6585bdf52
doc_id: 794860
cord_uid: okahh9wg

BACKGROUND: Respiratory diseases in pigs are the main health concerns for swine producers. Similar to the diseases in human and other animals, respiratory diseases are primary related to morbidity and are the result of infection with bacteria, viruses, or both. B. bronchiseptica causes serious respiratory diseases in the swine airway track. However, the B. bronchiseptica-specific bacteriophage has diverse advantages such as decreasing antibiotic overuse and possible therapeutic potential against bacteria. OBJECTIVE: The objects of this study were to investigate the therapeutic effect of specific B. bronchiseptica bacteriophages and to identify genes related to bacteriophage signaling utilizing RNA microarrays in swine nasal turbinate cells. METHODS: Bor-BRP-1 phages were applied 24 h prior to B.bronchiseptica infection (1 × 10(7) cfu/ml) at several concentrations of bacterial infection. Cells were incubated to detect cytokines and 24 h to detect mucin production. And real-time quantitative PCR was performed to examine related genes expression. To determine the change of total gene expression based on B.bronchiseptica and Bor-BRP-1 treatment, we performed RNA sequencing experiments. RESULTS: The results showed that B. bronchiseptica induced increased expression of several inflammatory genes such as IL-1β, IL-6, and Muc1 in a dose-dependent manner. However, Bor-BRP-1 induced reduction of gene expression compared to the B. bronchiseptica induction group. In addition, microarrays detected Bor-BRP-1-altered inflammatory gene expression against B. bronchiseptica, reducing B. bronchiseptica-induced airway inflammation in swine epithelial cells. CONCLUSION: These results suggest that the specific bacteriophage has a therapeutic potential to defend against B. bronchiseptica infection by altering inflammatory gene expression profiles.

Bordetella bronchiseptica is a widespread gram-negative coccobacillus that causes respiratory diseases in wide range of mammalian species. Currently, it is a well-recognized etiological agent in many domestic and wild animal diseases including tracheobronchitis, conjunctivitis, rhinitis, and pneumonia (Garbal et al. 2016; Petrovic et al. 2017) . In particular, it is one of the primary causative agents of atrophic rhinitis and bronchopneumonia in swine (Moore and Kern 2001; Prüller et al. 2015; Chen et al. 2019) . The characteristic lesion of atrophic rhinitis is atrophy of the nasal septum and turbinate, leading to facial deformation (Petrovic et al. 2017; De Jong and Nielsen 1990; Brockmeier et al. 2002) . This impairs swine feeding and causes a reduction in growth performance resulting in great economic losses in the pig-rearing industry.

B. bronchiseptica is often isolated in mixtures with other bacterial or viral pathogens such as porcine reproductive and respiratory syndrome (PRRS) virus, porcine respiratory coronavirus, swine influenza virus (SIV), Haemophilus parasuis, Pasteurella multocida, and Streptococcus suis, from respiratory diseases in swine (Chen et al. 2019; Park et al. 2018; King et al. 2002) . B. bronchiseptica infection predisposes the upper respiratory tract to colonization with secondary pathogens including P. multocida and S. suis (Brockmeier 2004; Pósa et al. 2011) . P. multocida inefficiently colonizes by itself, while B. bronchiseptica acts as a precursor for infection with P. multocida. In addition, preinoculation with B. bronchiseptica prior to S. suis resulted in elevated isolation of S. suis as well as increased incidence of pneumonia and mortality in pigs (Mengelers et al. 1989; Vecht et al. 1992) . Presentation of B. bronchiseptica-infected mammalian animals can vary from asymptomatic cases to severe systemic infections. Generally, this infection leads to life-long chronic colonization in the host. Once infected, the infection is rarely eliminated from the host (Irie et al. 2005) .

Bacteriophages are tiny microorganisms that infect bacteria and proliferate within bacterial cells. After proliferation, the phages escape from the bacteria by destroying the bacterial cell wall (Forde and Fitzgerald 1999; Hermoso et al. 2007; Pepin et al. 2008) . Phages are very specific and can only selectively kill the target bacteria, not affecting other bacteria such as normal flora in animals (Leverentz et al. 2003; Rees and Dodd 2006) . Although B. bronchiseptica was first reported in 1911, several studies have been performed with its bacteriophages. Several reports have revealed that B. bronchiseptica is isolated through antibiotic resistance to β-lactams, macrolides, etc. (Petrovic et al. 2017; Prüller et al. 2015; Biswas et al. 2002; Dayao et al. 2014; Niemann et al. 2018) . Previously, we reported that Bor-BRP-1 down-regulated B. bronchiseptica-induced inflammation (Park et al. 2018) . However, the physiological mechanism by which Bor-BRP-1 suppresses B. bronchiseptica-induced inflammation remains unidentified.

In this report, we demonstrated that B. bronchisepticainduced inflammatory RNA expression can be effectively controlled by treatment with B. bronchiseptica-specific bacteriophage Bor-BRP-1 in vitro. Therefore, bacteriophages may function to suppress bacterial inflammation and may be a potential therapeutic candidate to eliminate B. bronchiseptica.

Swine nasal turbinate (PT-K75) cells were purchased from ATCC (CRL-2528; Manassas, VA) and cultured as previously described (Park et al. 2018) . Cells were prepared at 5 × 10 5 cells/well in 6-well plates and incubated for approximately one day to reach full confluence. Cells were serum starved for 24 h before treatment. bronchiseptica culture and Bor-BRP-1 bacteriophage production B. bronchiseptica strains (BB-01) were grown on Brain Heart Infusion (BHI) broth in a shaking incubator at 37 °C without antibiotics. The bacterial concentration was adjusted to 5 × 10 9 cfu/ml based on OD 600 nm. According to a previous study, the inoculum of B. bronchiseptica to PT-K75 cells was selected at 1 × 10 7 cfu/ml. Bor-BRP-1 bacteriophage was produced based on a previous study (Park et al. 2018 ). Bor-BRP-1 phages were applied 24 h prior to B. bronchiseptica-infection (1 × 10 7 cfu/ml) at 10 × concentration. One day after bacterial infection, the cells were lysed with TRIzol for RNA extraction.

Real-time PCR was carried out on an ABI 7300 Real-Time PCR system (Applied Biosystems, Foster City, CA) with TB Green™ Premix Ex Taq™ II (RR820A, Takara, Japan). Reactions were performed in a total volume of 20 µl including 10 µl of 2 × SYBR Green PCR Master Mix, 300 nM of each primer, and 1 µl of the previously reverse-transcribed cDNA template. The following primers were used: porcine IL-1β, forward 5′-GAG CAT CAG GCA GAT GGT GT -3′ and reverse 5′-CAA GGA TGA TGG GCT CTT CTTC -3′; porcine IL-6, forward 5′-GCT GCT TCT GGT GAT GGC TAC TGC C -3′ and reverse 5′-TGA AAC TCC ACA AGA CCG GTG GTG A -3′, porcine Muc1 forward 5′-GTG CCG ACG AAA GAA CTG -3′ and reverse 5′-TGC CAG GTT CGA GTA AGA G -3′, porcine GAPDH, forward 5′-AAG GAG TAA GAG CCC CTG GA -3′ and reverse 5′-TCT GGG ATG GAA ACT GGA AG -3′. All reactions were performed in triplicate. The relative quantity of mRNA was obtained using a comparative cycle threshold method and was normalized using GAPDH as an endogenous control.

Total RNA was quantified with an Agilent 2100 bioanalyzer using the RNA 6000 Nano Chip (Agilent Technologies, Amstelveen, The Netherlands). RNA quantification was determined using an ND-2000 Spectrophotometer (Thermo Inc., DE, USA). The library was constructed using a Quant-Seq 3′ mRNA-Seq Library Prep Kit (Lexogen, Inc., Austria) according to the manufacturer's instructions. In brief, 500 ng of total RNA was prepared, an oligo-dT primer containing an Illumina-compatible sequence at its 5′ end was hybridized to the RNA, and reverse transcription was applied. After RNA template degradation, second-strand synthesis was initiated by a random primer containing an Illumina-compatible linker sequence at its 5′ end. The double-stranded library was purified using magnetic beads for removal of all reaction components. The library was amplified to add the complete adapter sequences required for cluster generation. The finished library was purified from PCR components. High-throughput sequencing was performed as single-end 75 sequencing using NextSeq 500 (Illumina, Inc., USA). For data analysis, QuantSeq 3′ mRNA-Seq reads were aligned using Bowtie2 (Langmead and Salzberg 2012) . Bowtie2 indices were either generated from the genome assembly sequence or the representative transcript sequences for aligning to the genome and transcriptome. The alignment file was used for assembling transcripts, estimating their abundances, and detecting differential expression of genes. Differentially expressed genes were determined based on counts from unique and multiple alignments using coverage in Bedtools (Quinlan 2010) . The Read Count data were processed based on the quantile normalization method using EdgeR within R using Bioconductor (Gentleman et al. 2004 ). Gene classification was based on searches in Database for Annotation, Visualization, and Integrated Discovery (DAVID, http:// david .abcc.ncifc rf.gov/) and Medline databases (http:// www.ncbi.nlm.nih.gov/). A Porcine Quant-Seq microarray was performed using the customized service provided by eBiogen Inc, (Seoul, South Korea).

Data are presented as the mean ± SD of at least three independent experiments. Where appropriate, statistical differences among the groups were assessed by the ANOVA test. Differences with p value less than 0.05 was considered statistically significant.

First, we performed cell viability assays and confirmed that neither B. bronchiseptica nor Bor-BRP-1 exhibited cytotoxicity in a dose-or time-dependent manner (data not shown). B. bronchiseptica administration induced or 24 h (for mucin). RNA and cell lysates were isolated for specific quantitative real time-PCR (a) and ELISA (b). The figures are representative of three independent experiments. *p < 0.05 compared with B. bronchiseptica only inflammatory cytokine secretion and mucin gene expression to control the respiratory microenvironment in our previous study (Park et al. 2018 ). In addition, Bor-BRP-1 was challenged for 24 h in prior B. bronchiseptica infections. Bor-BRP-1 significantly decreased B. bronchiseptica-induced IL-1β, IL-6, and Muc1 expression (Fig. 1a,  b) . Muc1 is known as an anti-inflammatory mucin (Kato et al. 2017 ). Since the airway inflammatory environment was decreased by Bor-BRP-1, Muc1 gene expression decreased proportional to the status of anti-inflammation.

The results indicate that Bor-BRP-1 dramatically suppressed the B. bronchiseptica-induced physiological phenomena by eliminating the bacteria.

To examine the impact of gene expression according to Bor-BRP-1 treatment and B. bronchiseptic infection, RNA sequencing was performed on each group of PT-K 75 cells.

To investigate the effect of Bor-BRP-1 on B. bronchiseptica-infected cells, we classified the four sample groups as follows: (1) Control group with no treatment, (2) B. bronchiseptica infection group, (3) Bor-BRP-1 treatment group, and (4) Bor-BRP-1 treatment prior to B. bronchiseptica infection group. Cells were lysed using TRIzol to analyze the altered genes through a microarray.

Following normalization, a differential comparison of the expression profiles was performed. The normalized read counts were greater than 4, and the fold change (FC) was greater than 2, with a p value < 0.05. A total of 304 DEGs was identified. Firstly, hierarchical clustering of RNA-sequencing data was conducted to view the global features of the data. The samples from each treatment group were normalized with a z score, and the statistical significance for differential expression was displayed among the groups (Fig. 2) . Upregulation is indicated in red, and downregulation of a gene is indicated in blue. The results suggest that each comparison involved dramatic upor downregulation according to treatment. These results showed a dramatic increase in the group of B. bronchiseptica infection but also a significant decrease in the group of Bor-BRP-1 treatment prior to B. bronchiseptica infection.

Among a number of genes, we selected several inflammation-related genes that exhibited representative changes depending on the analysis (Fig. 3) . The results showed remarkable upregulation in the group of B. bronchiseptica infection but significant downregulation in the group with Bor-BRP-1 bacteriophage treatment prior to bacterial infection. These results agree with those in Fig. 1 . Taken together, these findings indicate that Bor-BRP-1 suppressed B. bronchiseptica-induced airway inflammation by altering the inflammation-related mRNA profile in swine epithelial cells.

To identify the changes in biological features due to bacterial infection and phage treatment, the identified DEGs were assembled and mapped into DAVID for functional enrichment analysis (Table 1 ). The DEGs were enriched mainly in biological processes associated with immune response. For further details regarding the identified DEGs, pathways were analyzed using KEGG. Pathway analysis showed that the metabolic pathways of PI3K-Akt signaling and TNF signaling were significantly enriched among the DEGs (p < 0.05, Table 2 ). This indicates that signaling of Bor-BRP-1 mainly affected the PI3K-Akt pathway, TNF signaling, and AMPK signaling to suppress specific B. bronchiseptica-induced inflammatory signaling in swine epithelial cells. Therefore, Fig. 2 Hierarchical clustering and heat map of RNA-seq data. Clustering of RNA-seq data generated in response to Control, B. bronchiseptica infection, Bor-BRP-1 treatment without B. bronchiseptica infection, and Bor-BRP-1 treatment 24 h prior to B. bronchiseptica induction groups. The heat map was generated using genes differentially expressed by at least twofold under two or more conditions. Downregulation of a gene is indicated in blue, and upregulation is indicated in red. Colors changing from blue to red indicate the differences in gene expression among the groups of control. Each row represents a single gene, and each column represents a group (color figure online)

Bor-BRP-1 may be a potential therapeutic candidate to eliminate B. bronchiseptica.

B. bronchiseptica-induced respiratory disease creates several serious problems for swine producers (Prüller et al. 2015; Nicholson et al. 2012 Nicholson et al. , 2014 . In this study, we evaluated the therapeutic potential of the specific bacteriophage Bor-BRP-1 on B. bronchiseptica-infected swine nasal turbinate cells. To the best of our knowledge, this is the first study that suggests alternation of micro RNA expression profiles in the Bor-BRP-1 bacteriophage effect on B. bronchisepticainfected cells. Treatment of Bor-BRP-1 for 24 h prior to bacterial infection showed significant inhibition of IL-1β and IL-6 secretion as well as Muc1 gene expression. Furthermore, we confirmed bacteriophage suppression of several genes of bacterial-induced increased expression by microarray analysis. The results of the study indicate that Bor-BRP-1 has therapeutic potential against B. bronchisepticaassociated respiratory disease.

As seen in Fig. 1 , IL-1β, IL-6, and Muc1 gene expression was upregulated in the B. bronchiseptica induction group, whereas the genes were decreased in a dose-dependent manner when the Bor-BRP-1 phage was applied for 24 h prior to B. bronchiseptica infection. This is interesting because IL-1β and IL-6 are important mediators of inflammation even though Muc1 is a membrane-bound mucin that serves an anti-inflammatory role (Kato et al. 2017) . It is implied that Bor-BRP-1 could suppress release of proinflammatory cytokines as well as anti-inflammatory responses. We still do not fully understand the mechanism by which the specific phage controls bacteria-induced airway inflammation in swine airway epithelial cells. Therefore, this study mainly focused on identifying the action mechanism of Bor-BRP-1.

We characterized micro RNA expression profiles known as epigenetic components that modulate gene expression post-transcriptionally. The hierarchical clustering map with 272 genes indicated differentially expressed gene cross-groups from downregulation to upregulation (Fig. 2) . Microarray revealed and qPCR confirmed increased gene expression in the B. bronchiseptica-infected group as well as decreased gene expression in the Bor-BRP-1-pretreated group, as shown in Figs. 1 and 3 . The result showed that the genes related to inflammation exhibited altered expression in the control, B. bronchiseptica induction, and Bor-BRP-1 treatment prior to B. bronchiseptica induction groups. In Fig. 3 Selected genes from RNA sequencing. The genes resulted in marked upregulation in the group of B. bronchiseptica induction, whereas dramatic downregulation was exhibited in Bor-BRP-1 bacteriophage treatment prior to bacterial infection. ACSM4, acyl-CoA synthetase medium-chain family member 4; RAB32, RAB32, member RAS oncogene family; GPN3, GPN-loop GTPase 3; BMPR1A, bone morphogenetic protein receptor type 1A; SERINC2, serine incorporator 2; LRP8, LDL receptor related protein 8; MRPS18B, mitochondrial ribosomal protein S18B; MPZL2, myelin protein zero like 2; MOSPD2, motile sperm domain containing 2. *p < 0.05 compared with the control; **p < 0.05 compared with B. bronchiseptica only addition, the selected genes exhibited similar traits assumed to be affected by the specific bacteriophage to suppress B. bronchiseptica-induced intracellular changes by killing the bacteria.

There are certain limitations in our study. First, we assessed the cellular viability of the host cells but could not confirm biosafety of the bacteriophage for clinical use. Second, we verified a protective effect of the bacteriophage as pretreatment to bacterial infection, but we did not evaluate the bacteriophage effect after bacterial infection. Further studies need to be performed to determine whether the bacteriophage could eliminate the bacteria after infection with or without antibiotic withdrawal.

In summary, we identified a protective effect of Bor-BRP-1 in B. bronchiseptica-infected swine nasal turbinate cells when applied 24 h prior to bacterial infection. Importantly, Bor-BRP-1 did not affect the viability of the host cells and could be a powerful alternative medicine to B. bronchiseptica-induced airway inflammation. Future studies will be performed to identify the intracellular mechanism of bacterial elimination according to phage challenge.

Bacteriophage therapy rescues mice bacteremic from a clinical isolate of vancomycin-resistant Enterococcus faecium

Prior infection with Bordetella bronchiseptica increases nasal colonization by Haemophilus parasuis in swine

Role of the dermonecrotic toxin of Bordetella bronchiseptica in the pathogenesis of respiratory disease in swine

Characterisation of a newly detected bacteriophage infecting Bordetella bronchiseptica in swine

Antimicrobial resistance in bacteria associated with porcine respiratory disease in Australia

Definition of progressive atrophic rhinitis

Bacteriophage defence systems in lactic acid bacteria

Occurence of Bordetella bronchiseptica in domestic cats with upper respiratory tract infections

Bioconductor: open software development for computational biology and bioinformatics

Taking aim on bacterial pathogens: from phage therapy to enzybiotics

Pseudomonas aeruginosa rhamnolipids disperse Bordetella bronchiseptica biofilms

MUC1: the first respiratory mucin with an anti-inflammatory function

Whatmore AM (2002) Development of a multilocus sequence typing scheme for the pig pathogen Streptococcus suis: identification of virulent clones and potential capsular serotype exchange

Fast gapped-read alignment with Bowtie 2

Biocontrol of Listeria monocytogenes on fresh-cut produce by treatment with lytic bacteriophages and a bacteriocin

In vitro antimicrobial activity of sulfonamides against some porcine pathogens

Atrophic rhinitis: a review of 242 cases

Phenotypic modulation of the virulent Bvg phase is not required for pathogenesis and transmission of Bordetella bronchiseptica in swine

The Bordetella bronchiseptica type III secretion system is required for persistence and disease severity but not transmission in swine

Antimicrobial susceptibility and genetic relatedness of respiratory tract pathogens in weaner pigs over a 12-month period

Bordetella bronchiseptica bacteriophage suppresses B. bronchiseptica-induced inflammation in swine nasal turbinate cells

Genomic evolution in a virus under specific selection for host recognition

The first Siphoviridae family bacteriophages infecting Bordetella bronchiseptica isolated from environment

Interaction of Bordetella bronchiseptica, Pasteurella multocida, and fumonisin B1 in the porcine respiratory tract as studied by computed tomography

Antimicrobial Susceptibility of Bordetella bronchiseptica isolates from Swine and companion animals and detection of resistance genes

BEDTools: a flexible suite of utilities for comparing genomic features

Phage for rapid detection and control of bacterial pathogens in food

Virulence of Streptococcus suis type 2 strains in newborn germfree pigs depends on phenotype

Conflict of interest Ga Young Park, Hyun Jin Yu, Jee Soo Son, Sang Joon Park, Hee-Jae Cha, and Kyoung Seob Song declare that we have no conflict of interest.Ethical approval This article does not contain any studies with human subjects or animals performed by any of the authors.