Skip to main content

Effects of sodium houttuyfonate on transcriptome of Pseudomonas aeruginosa



The purpose of this experiment is to analyze the changes of transcriptome in Pseudomonas aeruginosa under the action of sodium houttuyfonate (SH) to reveal the possible mechanism of SH inhibiting P. aeruginosa. We analyzed these data in order to compare the transcriptomic differences of P. aeruginosa in SH treatment and blank control groups.

Data description

In this project, RNA-seq of BGISEQ-500 platform was used to sequence the transcriptome of P. aeruginosa, and sequencing data of 8 samples of P. aeruginosa are generated as follows: SH treatment (SH1, SH2, SH3, SH4), negative control (Control 1, Control 2, Control 3, Control 4). Quality control is carried out on raw reads to determine whether the sequencing data is suitable for subsequent analysis. Totally 170.53 MB of transcriptome sequencing data is obtained. Then the filtered clean reads are aligned and compared to the reference genome to proceed second quality control. After completion, 5938 genes are assembled from sequencing data. Further quantitative analysis of genes and screening of differentially expressed genes based on gene expression level reveals that there are 2047 significantly differentially expressed genes under SH treatment, including 368 up-regulated genes and 1679 down-regulated genes.


Pseudomonas aeruginosa is a gram-negative bacterium, which can produce endotoxin, exotoxin, proteolytic enzyme and other substances and infect human and other organisms [1, 2]. At present, macrolide and aminoglycoside antibiotics are commonly used for curing the clinical infection of P. aeruginosa. However, with the emergence of drug resistance, P. aeruginosa are difficult to treat by common antibiotics. Thus, we are seeking for effective antimicrobial agents from traditional Chinese medicine to treat infection of P. aeruginosa. Previously, our research group have proved that sodium houttuyfonate (SH) can effectively inhibit the P. aeruginosa [3, 4]. Here, our aim is to investigate possible antimicrobial mechanism of SH by comparing the transcriptomic differences between SH drug and blank control groups.

The assemble transcriptome contains thousands of transcripts. Thus this study provides transcriptomic comparison between SH medication group and blank control group rather than comparisons of expression of several certain genes such as algD, algR, lasI, phzM, lasA and bdlA, in previous studies [4,5,6]. The difference of these transcriptome can be used as the basis for studying gene expression changes in SH treatment and control groups.

Data description

We cultured P. aeruginosa under two conditions, with 4 biological replications which were cultured independently under each condition:

  1. 1.

    ATCC 27853 was inoculated into LB liquid medium and cultured overnight at 37 °C. The culture were centrifuged for 1 min at 12,000 r/min, and pour out the supernatant, diluted with sterile water to 0.5 Maxwell colorimetric tube, and diluted to 107 times for later use. The SH was prepared according to our previous research [6]. The prepared SH samples in LB medium with 1 MIC (minimum inhibitory concentration) SH of 512 μg/ml were cultured for 24 h at 37 °C until OD600 was 0.6–0.8, and collected by centrifugation of 1 min at 12,000 r/min, and then rinsed with sterile water for 3 times. Then we placed the collected bacterial samples in the centrifuge tube, sealed it with sealing film, and sent the sample for RNA-seq with dry ice.

  2. 2.

    The samples of blank control were collected similarly as SH treatment samples except without drug treatment.

Totally 170.53 MB of transcriptome sequencing data is obtained after RNA-seq applying BGISEQ-500 platform. The original data of sequencing includes reads with low quality, linker contamination and high content of unknown base N are removed before data analysis to ensure the reliability of the results. This project used SOAPnuke [7], a filtering software independently developed by Huada Corporation, to make statistics and trimmomatic [8] to filter. Firstly, readers including connectors are removed. Then the reads with unknown base N content more than 5% are wiped off. Finally, the low-quality reads are removed (we define reads with a mass value of less than 10 and a proportion of more than 20% of the total number of bases in the reads as low-quality reads). The filtered “Clean Reads” are saved in FASTQ format. The file format corresponding to each sample is FASTQ format (Table 1).

Table 1 Overview of transcriptome data files

The original sequencing sequence data (fastq file), including reference genome information, can be obtained on NCBI. After the quality control of the original data, we used Bowtie2 [9] to compare clean reads to the reference gene sequence (Table S1), and then RSEM [10] was used to calculate the expression levels of genes and transcripts. After completion, totally 5938 genes are assembled from sequencing data.

In order to reflect the correlation of gene expression between samples, Pearson correlation coefficients of all gene expression amounts between every samples are calculated, and expression amount distribution analysis are performed. The obtained results are shown in Fig. S1A [13]. According to the gene expression level of each sample, a total of 2047 differentially expressed genes are detected by the threshold of fold changes > 2, Q value < 0.001, including 368 up-regulated genes and 1679 down-regulated genes. The results are shown in volcanic map of Fig. S1B [13].


The limitations of this data is that there is not a gradient comparison under multiple different concentrations of SH, and the transcriptome expression of P. aeruginosa under 24 h of culture is selected in this study, which may make the results inconsistent with other research results. In addition, in our previous studies [4, 5], the qRT-PCR results showed that genes of lasA, algD, algR are down regulated by SH treatment in P. aeruginosa. However, these genes are found to be below the detection threshold in this study. This may be due to the different technologies used in determination of the gene expression.

Availability of data materials

The data described in this data description can be freely and publicly accessed on the sequence reading file (SRA) with the item number of [11] and the gene expression omnibus (GEO) with the item number of [12]. We have used these data to conduct a large amount of additional analysis, and these results can be obtained from [13] and the corresponding authors according to reasonable requirements.



ribonucleic acid sequencing


sodium houttuyfonate


quantitative reverse transcription polymerase chain reaction


  1. 1.

    Denning GM, Iyer SS, Reszka KJ, et al. Phenazine-1-carboxylic acid, a secondary metabolite of Pseudomonas aeruginosa, alters expression of immunomodulatory proteins by human airway epithelial cells. Am J Physiol Lung Cell Mol Physiol. 2003;285(3):584–92.

    Article  Google Scholar 

  2. 2.

    Cui Q, Li F, Xing W, et al. Differential expression of two phenazine-producing loci mediated by deficiency of the global regulator rsmA in Psedomonas aeruginosa PAO1. Wei Sheng Wu Xue Bao. 2012;52:1326–34.

    CAS  PubMed  Google Scholar 

  3. 3.

    Cheng HJ, Wang CZ, Lu WB, et al. Effects of houttuyfonate sodium on eliminating adhesion of Psedomonas aeruginosa and forming biofilms. Zhongguo Zhong Yao Za Zhi. 2012;37:3297–301.

    CAS  PubMed  Google Scholar 

  4. 4.

    Wu DQ, Cheng H, Huang W, Duan Q. Sodium houttuyfonate in vitro inhibits the biofilm formation and alginate biosynthesis gene expression of Pseudomonas aeruginosa clinical strain. Exp Ther Med. 2015;10:753–8.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  5. 5.

    Wu D, Huang W, Duan Q, Li F, Cheng H. Sodium houttuyfonate affects production of N-acyl homoserine lactone and quorum sensing-regulated genes expression in Pseudomonas aeruginosa. Front microbiol. 2014;5:635.

    Article  PubMed  PubMed Central  Google Scholar 

  6. 6.

    Wang Tianming, Huang Weifeng, Duan Qiangjun, Wang Jian, Cheng Huijuan, Shao Jing, Li Fang, Daqiang Wu. Sodium houttuyfonate in vitro inhibits biofilm dispersion and expression of bdlA in Pseudomonas aeruginosa. Mol Biol Rep. 2019;46(1):471–7.

    CAS  Article  PubMed  Google Scholar 

  7. 7.

    Chen Y, Chen Y, Shi C, et al. SOAPnuke: a MapReduce acceleration-supported software for integrated quality control and preprocessing of high-throughput sequencing data. Gigascience. 2018;7(1):1–6.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  8. 8.

    Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for illumina sequence data. Bioinformatics. 2014;30(15):2114–20.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  9. 9.

    Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012;9(4):357–9.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  10. 10.

    Dewey CN, Bo L. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinform. 2011;12(1):323.

    CAS  Article  Google Scholar 

  11. 11.

    Wu D, Zhao Y. RNA-Seq of Pseudomonas aeruginosa. Sequence Read Archive. 2019. Accessed 11 Jun 2019.

  12. 12.

    Wu D, Zhao Y. Effect of sodium houttuyfonate on transcriptome of Pseudomonas aeruginosa. Accessed 28 Jun 2019.

  13. 13.

    Wu D, Zhao Y. SH transcriptome-S1.xlsx. Figshare. 2019.

Download references


We thank Professor Huang Xianqing of Shanghai Jiaotong University for providing P. aeruginosa strain and the sequencing facilities of Huada Gene Cloud Platform for the preparation and sequencing of the library.


The research is funded by National Natural Science Foundation of China (Grant No. 81503115) and the Discipline Construction Key Program of Anhui University of Chinese Medicine (DC18100042). The Funding is used to cover laboratory expenses, sample preparation and sequencing. These funding agencies have no role in research design, data collection, analysis of results, or manuscript writing.

Author information




DQW, JS and CZW designed and funded the experiment. YZ, YQS, LFM and JDW conducted experiments, and wrote the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Daqiang Wu.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Zhao, Y., Si, Y., Mei, L. et al. Effects of sodium houttuyfonate on transcriptome of Pseudomonas aeruginosa. BMC Res Notes 12, 685 (2019).

Download citation


  • RNA-seq
  • Pseudomonas aeruginosa
  • Sodium houttuyfonate
  • Transcriptome