Skip to main content
Download PDF

Comparative transcriptomic and proteomic analyses reveal upregulated expression of virulence and iron transport factors of Aeromonas hydrophila under iron limitation

BMC Microbiology volume 18, Article number: 52 (2018) Cite this article

Abstract

Background

Iron plays important roles in the growth, reproduction and pathogenicity of Aeromonas hydrophila. In this study, we detected and compared the mRNA and protein expression profiles of A. hydrophila under normal and iron restricted medium with 200 μM 2,2-Dipyridyl using RNA Sequencing (RNA-seq) and isobaric tags for relative and absolute quantification (iTRAQ) analyses.

Results

There were 1204 genes (601 up- and 603 down-regulated) and 236 proteins (90 up- and 146 down-regulated) shown to be differentially expressed, and 167 genes and proteins that showed consistent expression. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses revealed that the differentially expressed genes and proteins were mainly involved in iron ion transport, protein activity, energy metabolism and virulence processes. Further validation of the RNA-seq and iTRAQ results by quantitative real-time PCR (qPCR) revealed that 18 of the 20 selected genes were consistently expressed. The iron-ion absorption and concentration of A. hydrophila under iron-limited conditions were enhanced, and most virulence factors (protease activity, hemolytic activity, lipase activity, and swimming ability) were also increased. Artificial A. hydrophila infection caused higher mortality in cyprinid Megalobrama amblycephala under iron-limited conditions.

Conclusion

Understanding the responses of pathogenic Aeromonas hydrophila within the hostile environment of the fish host, devoid of free iron, is important to reveal bacterial infection and pathogenesis. This study further confirmed the previous finding that iron-limitation efficiently enhanced the virulence of A. hydrophila using multi-omics analyses. We identified differentially expressed genes and proteins, related to enterobactin synthesis and virulence establishment, that play important roles in addressing iron scarcity.

Background

Aeromonas hydrophila is an opportunistic pathogenic bacterium that is ubiquitous in aquatic environments and causes serious infections worldwide in cultured fishes, amphibians, reptiles, and even mammals [1,2,3,4]. The pathogenesis of A. hydrophila is multifactorial, causing disease with virulence factors, such as adhesins, cytotoxins, hemolysins, and proteases, and it has the capacity to form biofilms and alter metabolic pathways and gene expression under various host environments [5, 6]. Its virulence expression is closely related to the environment in which the bacteria live (in vivo and in vitro), nutrients, and so on [7]. For example, the nutrient iron deficiency in the host environment has been thoroughly documented as having a pronounced effect on the virulence of pathogens [8].

Iron is an indispensable element of most living cells that is involved in many cellular functions, including electron transportation and oxygen transportation. The quantity of iron has a great impact on biological processes, for instance, iron overload will result in iron toxicity to cellular components [9], especially for DNA damage, owing to the reactions between hydroxyl radicals and other biomolecules [10, 11]. However, iron deficiency can also cause malnutrition cell death in severe cases [12]. In vivo, iron is usually oxidized to an insoluble form due to its special physico-chemical properties, bonding with heme, ferritin, hemoglobin, and transferrin within the cells, and thus is not readily accessible to bacteria [13]. In response to this iron deficiency predicament, microorganisms have evolved a series of sophisticated mechanisms to compete against the host, such as the secretion of siderophores [14], to grab iron from transferrin, hemoglobin, and ferritin and maintain iron dynamic balance for bacterial growth, proliferation, and toxin secretion [15,16,17]. During the past decades, the bacterial iron acquisition system and virulence have attracted much attention. For example, CaFTR1-mediated iron-uptake was proven to be an important virulence factor of Candida albicans [18], iron-responsive transcriptional repressor PerR was required for full virulence in Staphylococcus aureus [19], and FeoB was determined to play an important role in Fe acquisition expression of virulence of Helicobacter pylori [20].

Pathogenic bacteria virulence factors under iron-restricted growth conditions have previously been published [21,22,23,24]. Proteomes and transcriptomes reflect gene expressions from two different levels, and their joint analysis provides more complete expression information about bacteria. Therefore, in this study, an iron stress model was established to maximize the simulation of iron deficiency environment in vivo, and the effects of iron-restricted stress on the growth and virulence of A. hydrophila were evaluated comprehensively by combining transcriptome and proteomics data.

Methods

Selection of iron chelator concentration and growth of A. hydrophila

A. hydrophila (NJ-35) was isolated from dead cultured cyprinid in Jiangsu Province, China [25], and kindly provided by Professor Yongjie Liu from the College of Veterinary Medicine, Nanjing Agricultural University, P.R. China. We selected 2,2’-Bipyridyl (Bip) (Sinopharm Chemical Reagent Co., Ltd., Shanghai, China) as the ferrous iron chelating agent because of its high cell membrane permeation and intracellular iron sequestering ability [26,27,28]. The accuracy and virulence of A. hydrophila NJ-35 were confirmed by 16S rRNA gene sequencing (Biological Engineering Technology Co., Shanghai, China) and lab infection assays, respectively. Six concentrations (0, 100, 200, 300, 400, and 500 μM Bip in normal tryptic soy broth medium (TSB; BD; final pH = 7.3)) were set to detect the optimal concentration according to the growth curve of A. hydrophila NJ-35. A. hydrophila NJ-35 was inoculated in 5 ml of normal TSB and incubated (28 °C, 24 h); bacteria cells were collected via centrifugation, washed three times with PBS, and then diluted to an optical density at 600 nm (OD 600) of 0.01 in 100 mL of normal TSB to culture (180 rpm, 28 °C).

Sample collection

A. hydrophila NJ-35 cells (OD 600 ≅ 0.8) in normal and iron-limited groups were collected by centrifugation (5000 rpm, 4 °C, 10 mins). The pellet was rinsed twice with saline and stored immediately at − 80 °C until further transcriptomic and proteomic analyses. The supernatant was retained, filtered (MILLEX®GP filter unit, 0.22 μm), and frozen at − 20 °C, and it was used for the following proteolytic and hemolytic activity analyses.

Determination of iron concentration

The atomic absorption spectrophotometry (GB/T 5009.90–2003) method [29] was used to the measure variations in the intracellular iron of A. hydrophila NJ-35 in normal and iron-limited groups, as well as the iron concentration in the broth. Samples were analyzed by the Jiangsu Provincial Food Safety Testing Co., Ltd.

Quantitative transcriptomics (RNA-seq)

(i) RNA isolation and mRNA purification

Total RNA was purified using an RNAqueous kit (Thermo Fisher Scientific, San Jose, CA, USA) according to the manufacturer’s instructions. The RNA concentration and integrity (RIN) were measured following the previous description of Wang et al. [30]. The mRNA was enriched using a MICROBExpress Kit (Ambion, USA) [31], and determined on Agilent 2100 Bioanalyzer.

(ii) cDNA Synthesis, Illumina sequencing and library construction

Bacterial mRNA was fragmented using an RNA fragmentation kit (Illumina, San Diego, CA, USA). Double-stranded cDNA was synthesized using SuperScript II Reverse Transcriptase (Invitrogen, Carlsbad, CA) according to the manufacturer’s recommendations. Libraries were prepared with the standard protocol of the TruSeq RNA Sample Prep v2 Low Throughput (LT) kit. Paired-end sequencing was processed by the Hiseq™2000 (Illumina, San Diego, CA, USA) sequencer.

(iii) Bioinformatics Analyses

The assembled reads were mapped to the complete genome of the A. hydrophila NJ-35 strain (http://www.ncbi.nlm.nih.gov/nuccore/CP006870.1). The QC of alignment was produced based on the standard generated by Qin et al. [31]. The gene expression level was calculated using the RPKM method (fragments per kb per million reads) [32]. Differentially expressed genes (DEGs) were identified with EdgeR software [33], and used to generate statistical information such as expression level, fold change, p-value and FDR (false discovery rate). The specific filter conditions of DEGs were: log2(fold change) ≥ 2, p < 0.05 and bcv (biological coefficient of variation) = 0.01.

GO enrichment analyses of DEGs were performed on website (http://www.geneontology.org/). The calculation method, p-value formula and enrichment score were analyzed according to the method reported by Yan et al. [34].

Additionally, the DEGs were subjected to KEGG enrichment analyses [35] to identify their main metabolic pathways. The formula used for calculation was the same as that in the GO analyses.

Quantitative proteomics (iTRAQ)

(i) Protein extraction, quantization, and SDS-PAGE electrophoresis

The extract of whole cellular protein was conducted according to Isaacson et al. [36] with some modification. The bacterial cells pellets were suspended in cooled acetone (1 h, − 20 °C), centrifuged (15,000×g, 15 mins, 4 °C), and dried with a vacuum freeze dryer. The samples were resuspended in cold saturated-phenol (pH 7.5) and shaken (30 mins, 4 °C). The upper phenolic phase was collected by centrifugation (5000×g, 30 mins, 4 °C), 5 volumes of cold 0.1 M ammonium acetate in methanol was added, and then it was stored (1 h, − 20 °C). After centrifugation (5000×g, 30 mins, 4 °C), the pellets were washed and mixed with 2 volumes of ice-cold methanol. The pellets were centrifuged, dried and dissolved in lysis solution (1 h, 30 °C). The supernatants were isolated by centrifugation (15,000×g, 15 mins). The protein concentrations were measured with the BCA method [37], after which they were stored at − 80 °C for iTRAQ analyses. Additionally, 10 μg samples were subjected to 12% SDS-PAGE, visualized and then scanned according to Candiano’s protocol [38].

(ii) protein samples preparation and labeling

The filter-aided sample preparation (FASP) method [39] was adopted for enzymatic hydrolysis of the proteins (100 μg). After 50 μL trypsin (50 ng/μL) digestion, peptides were labeled according to the manufacture’s protocol for 8-plex iTRAQ reagent (AB SCIEX, USA).

(iii) 2D-LC-MSMS analyses

RPLC analyses

The dried samples were resuspended with 100 μL buffer A, after which reversed-phase liquid chromatography (RPLC) was employed on an Agilent 1200 HPLC System (Agilent). Separation was conducted according to the method of You et al. [40]. The first segment was collected from 0 to 5 mins, after which each additional segment was collected at a 4.5 min interval for 6–45 min, while the last segment was collected from 46 to 50 mins for a total of 10 segments. Each segment was dried and used for subsequent RPLC-MSMS analyses.

RPLC-MSMS analyses

In brief, samples were resuspended with Nano-RPLC buffer, filtered through a C18 nanoLC trap column, and a Chromxp C18 column (75 μm × 15 cm, C18, 3 μm 120 Å). The Eksigent nanoLC-Ultra™ 2D System (AB SCIEX) was used to perform the online Nano-RPLC. Triple TOF 5600 system (AB SCIEX, USA) was used to analyze MS data combined with Nanospray III source (AB SCIEX, USA).

(iiii) protein identification and quantification

Data were processed with the Protein Pilot Software v. 5.0 (AB SCIEX, USA) against the NCBI database using the Paragon algorithm [41]. The results of protein quantification were obtained by the matching of tandem mass spectrometry (MS) data and theoretical data, and was performed with the search option: emphasis on biological modifications.

An Orbitrap Elite high-resolution mass spectrometer (Thermo Fisher Scientific, USA) was used for ITRAQ quantitative proteomic analyses. Normalized high-energy collision dissociation (HCD) was performed, with the collision energy set at 30%. A protein database search and quantification were performed using Maxquant 1.5.1.0 (Thermo Fisher Scientific, USA). The protein database contained 4119 proteins (https://www.ncbi.nlm.nih.gov/genome/?term=Aeromonas+hydrophila, GCF_000014805.1_ASM1480v1_protein.faa). Oxidation (M) and acetyl (protein N-term) were used as the variable modifications and carbamidomethyl (C) was the fixed modification. The MS/MS tol. (FTMS) was 20 ppm. The protein quantitation, peptides matching and the functional annotations of DEPs were performed according to the method reported by Yao et al. [24].

Primer design, quantitative real-time PCR (qRT-PCR) validation

All of the sequence-specific primers of the target genes for qRT-PCR analyses were designed using Primer 5.0 based on the obtained fragment (Table 3). The mRNA level of rpoB was used as an internal reference because of its stable expression according to Zhang et al. [42].

Total RNA from A. hydrophila was extracted using RNAiso Plus (TaKaRa, Japan), and measured using a Nanodrop 2000 (Thermo Fisher Scientific, USA), the RNA concentration of each sample were diluted to 40 ng/μL, and then 2 μg of the total RNA was subjected to the following quantitative analysis with a One Step SYBR® PrimeScript® Plus RT-PCR Kit (TaKaRa, Dalian). Triplicate quantitative assays were performed on each type of cDNA using the ABI 7500 Real-time PCR System (Applied Biosystems, Foster City, CA, USA) and analyzed with the two-standard curve method.

Proteolytic activity

Proteolytic activity was measured by an azocasein assay method of Swift et al. [43] and Chu et al. [44], with some modifications. Briefly, 150 μL of normal group and iron-limitation group NJ-35 culture supernatants were added to 1 ml of 0.3% azocasein (Sigma, St. Louis, USA) in 0.05 M Tris-HC1 and 0.5 mM CaCl2 (pH 7.5), then they were incubated (37 °C, 30 mins) respectively. Precooling trichloroacetic acid (l0%, 0.5 ml) was then added to stop the reaction, after which the samples were allowed to stand for 15 mins at room temperature, then they were centrifuged (12,000 rpm, 10 mins, 4 °C) to remove the precipitate. Next, 500 μL of the supernatants were added to an equal volume of NaOH (1 mol/L). The supernatants (200 μL) were subsequently transferred to a 96-well tissue culture plate, after which the absorbance (OD400) of the supernatant was measured. The proteolytic activity was calculated using the following equation: proteolytic activity = OD400nm sample – OD400nm blank control (normal TSB/iron limitation TSB).

Hemolytic activity

Hemolytic activity was determined as previously described [45, 46], and sheep blood (Ping Rui Biotechnology, China) was prepared by washing thrice with PBS. Washed sheep blood (10 μL) was added to 490 μL of the experiment supernatants (sample), normal TSB/iron limitation TSB (blank control), 1% (v/v) Trinton X-100 (positive control), or PBS (phosphate buffer solution, negative control). After 30 mins of incubation at 37 °C, all of the samples were centrifuged (5000 rpm, 10 mins) at room temperature. The supernatants (200 μL) were then transferred to a 96-well tissue culture plate, after which the absorbance of hemoglobin released for each solution at 540 nm was measured. The percentage of hemolysis was calculated using the following equation: hemolysis (%) = (OD540nm sample - OD540nm blank control)/ (OD540nm positive control Trinton X-100 - OD540nm negative control PBS).

Lipase activity

Bacterial cells were centrifuged and washed with PBS, after which 5 μL of bacterial fluid was used to inoculate the LB medium containing a 1% mass fraction of Tween 80. Samples were then incubated at 28 °C for 24 h, after which they were observed for lipase production, which was indicated by a white precipitate zone around the colony.

Motility

The target bacteria were centrifuged and washed with sterilized PBS. Next, 5 μL of bacterial fluid was dropped onto LB semisolid agar plates containing 0.3% agar (to determine swimming ability) and 0.5% agar (to determine swarming motility). The LB plates were subsequently sealed with parafilm and incubated at 28 °C for 24 h (three parallel groups were set up for each group). At the end of the culture period, the migration distance from the colony edge to the colony center was determined. The experiment was repeated three times.

Infection assays in vivo

A health check was conducted and healthy M. amblycephala (50 ± 5 g) were obtained from the Nanquan Experimental Station of the Freshwater Fisheries Research Center (Chinese Academy of Fishery Sciences, China) and acclimatized in circulating water system with thermo-control for 2 weeks before use. Fish were given commercial feed. The water temperature fluctuated between 27.5–28.5 °C, with a pH between 7.2–7.8, and the DO was about 5.5 mg/L.

Strain NJ-35 was inoculated aseptically into normal TSB medium and iron-limitation medium and then incubated for 18 h at 28 °C while shaking at 180 rpm. The artificial challenge experiment was performed as the previous report [47]. To determine the 50% lethal dose (LD50) [48], five groups of 20 M. amblycephala each were injected intraperitoneally with 150 μL of serial ten-fold diluted bacterial suspensions (1 × 109, 108, 107, 106, and 105 CFU·mL-1 measured by turbidimeter (Yue Fung Instrument Co., Ltd., Shanghai, China)), which were diluted with 0.9% saline. Next, an experimental group and a control group were injected intraperitoneally with 150 μL A. hydrophila (LD50) iron-limited and A. hydrophila (LD50) basal, respectively, and the virulence was compared. Three replicate tanks per challenge isolate (containing 20 fish each) were used to calculate survival (from a total of 60 fish per isolate). The mortality of the fish of experimental groups and control groups were monitored (7 days), and the activity and behavior were recorded daily; pathogenic bacteria were isolated and identified from the lesion tissues of dead fish as the judging standard.

Results

Growth of A. hydrophila under different iron-limitation medium

The effects of different concentrations of Bip on the growth of A. hydrophila are shown in Fig. 1. When compared with the control group, inhibitory effects were observed in the Bip addition groups, and higher Bip concentrations delayed the time of entering the logarithmic phase and reduced the maximum. When the Bip concentration was 500 μM, the growth of A. hydrophila was totally inhibited for at least 24 h. Due to the significant inhibition and higher cells concentration, 200 μM Bip was chosen as the proper iron-limitation concentration for subsequent analyses.

Fig. 1
figure 1

Effect of Bip supplementation on A. hydrophila growth. Growth curve (OD600) of A. hydrophila NJ-35 grown in TSB medium in the presence of 0, 100, 200, 300, 400, and 500 μM Bip

Full size image

Expression profile of iron-limited A. hydrophila

Based on the transcripts of A. hydrophila, 4327 genes were identified and quantified (Table 1). After filtering with FDR, 1204 genes were found to be differentially expressed between the control and iron-limitation groups. Detailed information for most of the DEGs is shown in Table 2. In comparison, the quantity of down-regulated DEGs detected (603) was greater than that of the up-regulated genes (601). A total of 2244 proteins were identified; 2012 were quantified and 1946 were correlated with the transcripts. Additionally, while compared with the control group, a total of 236 DEPs (90 up-regulated and 146 down-regulated) were identified in the iron-limitation groups with an at least 2-fold difference, and 167 of the DEPs were correlated to the corresponding DEGs, which have the same trends. Fewer DEPs are probably due to the removal of some proteins that were secreted by A. hydrophila NJ-35 in the supernatant of the experimental design.

Table 1 Overall features of the iron-limitation responsive expression profile
Full size table
Table 2 List of differentially expressed genes under iron restriction
Full size table

Integration analyses of transcriptome and proteome

To identify robust pathways that were corroborated by both datasets, we integrated the differentially expressed transcripts and proteins to find the corresponding genes and proteins, and the results are listed in Additional file 1: Excel S1.

The distribution of the corresponding mRNA: protein ratios is shown in a scatterplot of the log2-transformed ratios. As shown in Fig. 2, almost all of the log2 mRNA: log2 protein ratios are concentrated at the center of the plot, where mRNA and protein levels did not vary above 2-fold. Integration analyses of transcriptome and proteome data revealed that 67 genes and their corresponding proteins were up-regulated, while 94 were down-regulated, reflecting significant changes and showing a strong correlation between the transcripts and proteins. Overall, 680 transcriptomes showed DEGs with no difference in proteins, while 35 transcriptomes showed different proteins but no difference in genes. Conversely, the expression of the following six genes and proteins was opposite (e.g., when the gene was upregulated, the protein was downregulated and vice versa): (U876_04575, YP_857861.1), (U876_17130, YP_855747.1), (U876_17135, YP_855746.1), (U876_19295, YP_855421.1), (U876_20135, YP_855265.1), and (U876_21295, YP_855025.1). This exception can be caused by regulation at several levels, such as post transcriptional processing, degradation of the transcript, translation, post-translational processing and modification. In summary, most of the trends in DEP abundance were consistent with the DEG data.

Fig. 2
figure 2

Relationship patterns of all of the quantitative mRNA and protein. In the nine-quadrant diagram, the abscissa is the protein expression and the ordinate is the gene expression. Each color denotes a log2 mRNA ratio and a log2 protein ratio. Gray (filtered) represents genes and proteins with no significant difference, red (Cor_up) indicates up-regulated genes and proteins, green (Cor_down) indicates down-regulated genes and proteins, purple (Opposite_Sig) indicates that DEGs and DEPs show opposite up- and down- regulation and blue (Single_Sig) indicates that one of the genes and proteins differ

Full size image

Functional classification of enriched DEGs and DEPs by GO and KEGG

GO enrichment analyses were used to classify the enriched DEGs and DEPs between the control and iron-limitation groups using bioinformatics methods, and the results are listed in Additional file 2: Excel S2 and Additional file 3: Excel S3, respectively. As shown in Fig. 3, the following three ontologies (molecular function, cellular component and biological process) were observed.

Fig. 3
figure 3

GO enrichment analyses of DEGs and DEPs Control group vs Iron-Limitation group. GO term analyses of transcriptomics (a) and proteomics (b) that were catalogued as Biological Process, Cellular Component, and Molecular Function

Full size image

DEGs were distributed in up to 1460 GO terms, while DEPs were classified into 402 GO terms. In this case, GO terms related to bacteria energy metabolism, iron ion transport, and virulence. Based on the ‘−log10Pvalue’, most of the GO terms in the biological process category were associated with energy metabolism (Fig. 3a and b). Additionally, six genes were categorized as ‘glycerol catabolic process’ (GO: 0019563), three as ‘propionate catabolic process, 2-methylcitrate cycle’ (GO: 0019629), five as ‘oxidative phosphorylation’ (GO: 0006119), and five as ‘respiratory electron transport chain’ (GO: 0022904). Regarding proteomics, DEPs were mainly involved in the synthesis and transport of iron ions and proteins, particularly the following GO terms: ‘iron assimilation’ (GO: 0033212), ‘ion transport’ (GO: 0006811), ‘enterobactin biosynthetic process’ (GO: 0009239), ‘protein secretion’ (GO: 0009306), ‘protein transport’ (GO: 0015031), and ‘electron transport chain’ (GO: 0022900).

In the cellular component category (Fig. 3a and b), three genes were categorized as ‘glycerol-3-phosphate dehydrogenase complex’ (GO: 0009331), five as ‘proton-transporting ATP synthase complex, catalytic core F(1)’ (GO: 0045261), four as ‘proton-transporting ATP synthase complex, coupling factor F(o)’ (GO: 0045263), and seven as ‘bacterial-type flagellum hook’ (GO: 0009424). Regarding proteomics, DEPs were mainly classified in the cell membrane and cytoplasm of GO terms, including ‘integral component of membrane’ (GO: 0016021), ‘plasma membrane’ (GO: 0005886), ‘cell outer membrane’ (GO: 0009279), ‘cytosol’ (GO: 0005829), and ‘cytoplasm’ (GO: 0005737).

In the molecular function category (Fig. 3a and b), 11 genes were categorized as ‘receptor activity’ (GO: 0004872), three as ‘energy transducer activity’ (GO: 0031992), three as ‘cytochrome o ubiquinol oxidase activity’ (GO: 0008827), four as ‘siderophore uptake transmembrane transporter activity’ (GO: 0015344), and three as ‘siderophore transmembrane transporter activity’ (GO: 0015343). Regarding proteomics, DEPs were mainly related to protein activity and binding capacity, including ‘siderophore transmembrane transporter activity’ (GO: 0015343), ‘receptor activity’ (GO: 0004872), ‘iron ion binding’ (GO: 0005506), ‘heme binding’ (GO: 0020037), ‘metal ion binding’ (GO: 0046872), and ‘porin activity’ (GO: 0015288). In summary, GO term enrichment analyses further explained that metabolism, biosynthesis, transmembrane transport and redox homeostasis should be tightly regulated.

Enriched KEGG terms are listed under Additional file 4: Excel S4 and Additional file 5: Excel S5, as transcriptomics and proteomics, respectively. When compared with the whole genome, a total of 624 genes were present in the 139 KEGG pathways as DEGs, and we selected the 20 most critical KEGG pathways according to the enrichment scores (Fig. 4a). The up-regulated KEGG pathways included 78 genes under the category of ‘ABC transporters’ (ko02010), 20 genes under ‘TCA cycle’ (ko00020), and 38 genes under ‘quorum sensing’ (ko02024). We inferred that transport, energy production and bacteria interact with each other and may play important roles via stress responses that are regulated through several pathways. The down-regulated KEGG pathways included 47 genes categorized as ‘Ribosome’ (ko03010), 71 as ‘Carbon metabolism’ (ko01200), 31 as ‘Pyruvate metabolism’ (ko00620), and 35 genes as ‘Oxidative phosphorylation’ (ko00190), which confirmed that bacteria slowed down material synthesis and life activities. With respect to proteomics, a total of 41 proteins were detected in the 34 KEGG pathways by DEP, while only eight pathways were found to be significantly enriched by filtration (Fig. 4b). The up-regulated KEGG pathways included three that were labeled under ‘biosynthesis of siderophore group nonribosomal peptides’ (aha01053) and 10 that were labeled under ‘ABC transporters’ (aha02010), indicating clear changes in synthesis and transportation of siderophores. The down-regulated KEGG pathways included 11 proteins that were classified as ‘oxidative phosphorytation’ (aha00190), six as ‘butanoate metabolism’ (aha00650), five proteins as ‘TCA cycle’ (aha00020), five as ‘pyruvate metabolism’ (aha00620), seven as ‘carbon metabolism’ (aha01200), and six as ‘two-component system’ (aha02020), indicating the bacteria repress energy metabolize to adaptive constraint environment. Conversely, the total number of DEPs among them was far smaller than that of the DEGs, and most DEGs and DEPs were down-regulated.

Fig. 4
figure 4

KEGG enrichment analyses of DEGs and DEPs Control group vs Iron-Limitation group. KEGG enrichment analyses of transcriptomics (a) and proteomics (b)

Full size image

Clustering of virulence genes and proteins in A. hydrophila in iron-limited medium

According to the bioinformatics analyses, we found that there were 60 virulence factors in the differential genes, which mainly fell under the category of synthesis of iron carriers (U876_01620, U876_18555, U876_21285, U876_21455, U876_23515, and U876_24445), motility of flagella (U876_20435, U876_07265, U876_07270, and U876_07305), and generation of hemolysin (U876_04005, U876_15265, U876_16300, and U876_16315). Heat map analyses (Fig. 5) were used to visualize genes and proteins, and the results indicated a comprehensive impact and clear changes in the regulation of virulence factors.

Fig. 5
figure 5

Clustering of 60 mainly related virulence genes and proteins. Numbers are listed as the log2 value of difference multiples. Expression differences are shown in different colors; red indicates up-regulation, while green indicates down-regulation. A heatmap was used to visualize the genes and proteins that were related to virulence factor (hemolysis, secretion system, lipase, phospholipid, serine-type peptidase, metallopeptidase, flagellum, polysaccharides, siderophore transporter, quorum sensing, and outer membrane)

Full size image

Validation of selected DEGs/DEPs by qRT-PCR analyses

To further evaluate the expression of genes in an iron-limited environment, 20 virulence genes (13 up-regulated and seven down-regulated genes) together with reference genes (rpoB) were selected for investigation based on their expressions, which were measured by real-time quantitative PCR (RT-qPCR) (Table 3) according to the results of the GO analyses. These selected genes were involved in virulence factors, hemolysis, secretion systems, lipases, phospholipids, serine-type peptidases, metallopeptidases, flagella, polysaccharides, siderophore transporters, quorum sensing, and outer membrane production.

Table 3 Primers and sequences used in this study for q-PCR
Full size table

The results of qPCR showed that the majority of the selected virulence factors (90%, 18/20) were consistent with the transcriptome data. Notably, five virulence-related factors, U876_15265 (hemolysin, log2FC = 3.80), U876_15575 (secretin, log2FC = 5.00), U876_18585 (hemin ABC transporter substrate-binding protein, log2FC = 4.71), U876_20975 (transcriptional activator protein AhyR/AsaR, log2FC = 3.33), and U876_09860 (2,3-dihydroxybenzoate-AMP ligase, log2FC = 6.20) were shown to be significantly up-regulated (log2FC > 3.00) under iron-limited conditions. Moreover, two selected genes, U876_07270 (flagellar hook protein FlgE) and U876_12225 (murein transglycosylase A), showed appositive results to the RNA-seq data, which might have been due to differences in the analyses methods.

Determination of iron concentration

Atomic absorption spectrophotometry revealed that the medium iron concentration of 0.44 mg/100 g in the normal TSB group was higher than 0.28 mg/100 g in the iron-limitation group, indicating that iron scavenger 2,2- bipyridine has a higher efficiency. After bacterial growth, the medium iron content of the normal TSB group was higher than that of the iron-limitation group. Surprisingly, the concentration of 0.664 mg/100 g in the normal TSB group strain cell was lower than 0.998 mg/100 g in the iron-limitation group strain cell. All of the results are shown in Table 4.

Table 4 Determination of iron concentration under two culture conditions
Full size table

Effect of iron-limitation on virulence factors production in A. hydrophila

As shown in Table 5, the total protease activity in supernatants from A. hydrophila NJ-35 growing without Bip was 0.105 (OD400 nm), whereas the presence of Bip resulted in a significant increase in protease activity to 0.36 (OD400 nm) (Fig. 6a). When compared with the control group, the hemolytic activity of A. hydrophila NJ-35 was significantly enhanced under iron limitation, indicating that NJ-35 produced 83.8% more hemolysin (Fig. 6b). To observe the hemolysis ability, sheep blood agar plates were used for rough detection. A. hydrophila NJ-35 under iron limitation generated a large hemolytic zone on the blood agar plates compared to the control group, but the lipase activity and swarming motility did not differ significantly (Table 5). Interestingly, the swimming ability of the bacteria was strong under iron limiting conditions, which could reflect attempts to move to areas with more suitable conditions (Table 5).

Table 5 Effect of iron limitation on A. hydrophila extracellular enzyme activity and motility
Full size table
Fig. 6
figure 6

Effect of control and iron-limitation conditions on A. hydrophila NJ-35. a Total protease, and (b) hemolytic activity. The data represent the mean values of three independent experiments and are presented as the means ± SD

Full size image

Infection assays

The isolated pathogenic bacteria were A. hydrophila after morphological, physiological and biochemical, molecular identification. Megalobrama amblycephala injected with A. hydrophila NJ-35 showed distinct mortality rates under iron and non-iron limited conditions (Fig. 7). Although the difference was not significant, the survival rate in the group injected with A. hydrophila was substantially higher (by 19.77%) than that of the iron-limitation group at four days post-challenge.

Fig. 7
figure 7

Kaplan-Meier survival analyses of Megalobrama amblycephala challenged with A. hydrophila NJ-35 from normal and iron-limited media. Data represent accumulative fish mortality in three replicates

Full size image

Discussion

Comparative transcriptomic and proteomic analyses

The survival and proliferation of bacteria was sensitive to environment factors. Many environmental stress factors, e.g., pH, temperature, oxygen, acidity and salinity [49, 50] significantly affected the expression of virulence. Iron limitation is an important external stimulus [51] that has profound impacts on almost all bacteria. The culturability and growth rate of A. hydrophila were reduced under iron-limited conditions [52]; however, bacterial multiplication was enhanced after injecting exogenous iron into experimentally infected animals, and the virulence of pathogens including Vibrio cholerae, Pseudomonas aeruginosa, Klebsiella pneumoniae, and Mycobacterium tuberculosis was established with sufficient iron [8, 53]. A. hydrophila establishes virulence through many mechanisms [54], including iron-binding systems, secretion systems, biofilm formation, flagella and pili adhesion, structural proteins, phospholipids, polysaccharides, hemolysis, collagenase, serine protease, metalloprotease, enolase, lipase, and nucleases [5, 6]. The pathogenesis of diseases involves most virulence factors [1], beginning with molecular changes at the micro level and progressing to phenotypic changes at the macro level [55]. Under iron-limited conditions, virulence genes and proteins were up-regulated more than down-regulated (Fig. 5), suggesting that virulence expression was enhanced in A. hydrophila to compensate for iron insufficiency, which was confirmed in F. tularensis [56]. These virulence factors exerted synergistic effects [57] and contributed to the production of toxins. The results of the infection assays further confirmed this conclusion (Fig. 7). The ferric uptake regulator (Fur) is a negative regulator in iron acquisition systems [58] that controls the expression of 90 virulence and metabolic genes [7, 15, 59]. For example, the biosynthesis of rhizoferring, an iron siderophore in F. tularensis, is regulated by operon fslABCDEF [60]. In this study, the expression of the fur gene (U876_15170) was up-regulated (log2FC = 0.3187). This phenomenon could be explained by the higher iron concentration in bacterial cells of the iron-limited group. At the sampling time-point, more iron was stored in the iron-limited group, after which fur was up-regulated to reduce the iron absorption [61].

Iron homeostasis was coordinated by the absorption, transport, utilization, and storage of iron ions [62]. A. hydrophila utilized multiple iron sequestration systems to hijack host iron ions [63]. Under an iron deficient environment, A. hydrophila secreted large amounts of iron transporters and iron-specific scavenger-siderophores. The same results were confirmed by transcriptome analyses of Bacillus cereus ATCC 10987, which showed the upregulation of predicted iron transporters in the presence of 2,2-Bipyridine [64]. As an important virulence characteristic of pathogens to both animals and plants [65], siderophores were formed and played a major role in microbial iron acquisition. Siderophore-assisted iron uptake and reductive iron assimilation are both induced upon iron starvation [58]. In previous studies, A. hydrophila was found to secrete siderophores to compete with transferrin in vivo to meet the iron demand required for growth and virulence [66]. Measurement of the iron concentration confirmed that the iron chelating ability of bacterial siderophores was notable (Table 4), because E. coli [67] and A. hydrophila synthesize and secrete enterobactin siderophores [68] in response to iron starvation. Enterobactin synthase subunit E (entE), which is encoded by entA, entB, and entC genes, is a key enzyme involved in the synthesis of isochorismate synthase. In both E. coli and A. hydrophila, a 22 kB gene cluster including entD-fepA-fes-entD-fepE-fepC-fepG-fepD-fepB-entC-entE-entB-entA-ybdA genes encodes proteins responsible for the synthesis and transport of enterobactin [69]. During this process, the entE polypeptide is responsible for activating the DHBA carboxylate group with ATP by forming the enzyme-bound 2,3-dihydroxybenzoyadenylate as an intermediary in the biosynthetic pathway [70]. Genes with similar enterobactin transport functions (iroN, fepC, cirA, fepC, and iroC) were also found in Salmonella enterica [71]. After differential analyses of the genes and proteins, we found that the entE expression level of gene U876_09860 (log2FC = 9.39) and protein YP_856992.1 (log2FC = 15.46) had increased significantly during the biosynthesis of the siderophore subunits (ko01053). Upon GO term analyses of the DEGs, the entE gene and protein expression levels were not increased significantly, which may have been because of differences in the analyses methods and software. Ferritin is the major iron storage protein in A. hydrophila [72]. The data demonstrated that ferritin (U876_00270, log2FC = − 0.4088) and bacterioferritin (U876_02285, log2FC = 13.4043) participated in iron ion transport and storage, which may benefit the survival of bacteria. The up-regulation of this protein may be responsible for the increased intracellular iron concentration in A. hydrophila. The expression levels of bacterioferritin in different isolates, including F. tularensis, also varied [73, 74]. The TonB mechanism is essential to the virulence of avian pathogenic E. coli [75], indicating that a specific TonB-dependent outer membrane receptor might be involved in the transport of iron from transferrin [76]. TonB-dependent outer membrane receptors TonB-2 (U876_00270, log2FC = 5.9844), AHA_4249 (YP_858666.1, log2FC = 6.1718), AHA_4250 (YP_858667.1, log2FC = 7.4891), and AHA_4251 (YP_858668.1, log2FC = 10.7778) were found to be required for the transfer of iron chelators and heme to the periplasm, followed by transport to the cytoplasm by ATP-binding cassette (ABC)-type transporters. Inorganic iron in the periplasm is transported to the cytoplasm by membrane transporters, such as Sfu ABC [77].

Iron influences a number of catalytic reactions involving cell energy metabolism in vivo, including respiration and nucleic acid replication [78]. Overall, when iron demand is not met, some enzymes related to metabolism, the regulation of protein synthesis, and the ability of A. hydrophila to utilize nutrients, such as carbohydrates, decreased. It has been hypothesized that decreased virulence might be caused by the loss of metabolic activity and the lack of toxin production [79, 80]. According to bioinformatic analyses conducted in this study, the energy generation system and electron respiration chain appeared to be depressed under iron starvation, which is consistent with previous quantitative proteomic analyses of A. hydrophila [24]. Iron scarcity reduces iron utilization in iron nonessential pathways, and limited iron is used for the synthesis of iron-containing enzymes involved in the citric acid cycle and the electron transport chain [81]. For example, the expression of NADP-dependent glyceraldehyde-3-phosphate significantly altered the antioxidant activity of bacteria, and NADPH is involved in the transformation of Fe3+ into Fe2+ in some of the identified bacteria [16]. Similar to S. pneumoniae in manganese limited environments [82], the metabolic activity of bacteria will become inert, so bacteria can survive in these environments for a long time [83]. Based on the high-throughput data analyses, it is apparent that 969 genes decreased, 905 genes increased, 146 proteins decreased, 90 proteins increased, the gene and protein ratio was down-regulated, and the regulation of bacteria itself was also used to interpret iron starvation.

Virulence evaluation of A. hydrophila under iron-limited environment

Many studies have shown that the virulence of A. hydrophila increased in response to iron deficiency [52]. Two aspects may contribute to the establishment of bacterial pathogenicity: invasiveness and toxin production [84]. The invasive ability of A. hydrophila is closely related to their motility, as well as the secretion of toxins, including aerolysin, hemolysin, and enterotoxin, and extracellular protease. To evaluate the virulence of A. hydrophila more comprehensively, we conducted an encompassing study of A. hydrophila hemolytic and enzymatic activity in vitro and lethality rate in vivo.

A. hydrophila pilus is an important coagulation factor and a major colonization factor that enables bacteria to adhere to host digestive epithelial cells during the invasion process. In terms of virulence establishment, pili-assisted adhesion bacteria were 10–200 times more effective than bacteria that do not express pili [85]. Flagella-mediated motility also promotes the initial stages of adhesion [86]. In this study, although the swimming ability of the control group was significantly stronger than that of the iron restriction group, swimming ability was enhanced under iron-limited conditions, indicating that A. hydrophila can overcome unfavorable conditions by accelerating their swimming and adhesion abilities, thereby enhancing their resilience to environmental restraints. Alternatively, these findings demonstrate the complexity of Aeromonas sp. virulence.

Lethal pathogenic extracellular products (ECPs) of A. hydrophila are produced to compete with rivals for limited iron resources [87]. After removal of ECPs by repeated washing with normal saline, the invasion and pathogenicity of the pathogenic bacteria to the host cells was reduced or even completely lost. As a typical ECP, hemolysin that is synthesized and secreted into the organism’s environment can dissolve various sources of iron by destroying intracellular red blood cells or hydrolyzing hemoglobin. Hemolytic activity was detected both in vivo and in vitro in septic animals, and beta hemolysins isolated from protease deficient strains of A. hydrophila were found to lead to the death of catfish [88]. Blood-plate hemolysis results showed that the hemolytic ability of A. hydrophila under iron deficiency was stronger than that under normal conditions, and it caused greater toxicity and damage to the host. The results also showed that iron exerted an inhibitory effect on extracellular hemolysin and protease activity. Notability, the presence of hemolysin alone does not cause disease [89].

The invasion of pathogenic bacteria was found to be significantly correlated with the level of corresponding enzyme production, and protease activity [90], which is consistent with the results of this trial. Not only can proteases degrade a variety of proteins to provide amino acids for bacterial survival and growth, but they can also directly cause tissue injury, resulting in the spread through the defense mechanism and evasion of the immune system of the host [91]. In addition, the A. hydrophila family of extracellular proteases can cooperate with other virulence factors [92] to activate other pathogenic factors. In this study, A. hydrophila NJ-35 under low-iron growth conditions were detected with higher protease activity than the control, demonstrating that iron scarcity can promote NJ-35 virulence factor expression.

Conclusion

In this paper, we simulated the iron restriction environment in the fish host, coalition analyzed the transcriptome and proteomics data of A. hydrophila, and identified the changes of enzyme activity, comprehensively revealed the pathogenicity of A. hydrophila increased. This study also provide a profound theoretical basis for the effect of exogenous iron preparation on the toxicity of bacteria.

Abbreviations

ABC:

ATPbinding cassette

Bip:

2,2’-Bipyridyl

BLAST:

Basic Local Alignment Search Tool

CA:

Citric acid

CID:

Collision-induced dissociation

DEGs:

Differentially expressed genes

DEPs:

Differentially expressed proteins

ECP:

Extracellular products

FDR:

False discovery rate

GO:

Gene Ontology

HCD:

High-energy collision dissociation

IDA:

Information dependent acquisition

iTRAQ:

Isobaric tags for relative and absolute quantification

KEGG:

Kyoto Encyclopedia of Genes and Genomes

LC MS/MS:

Liquid chromatography tandem mass spectrometry

LD50 :

50% lethal dose

MS:

Mass spectrometry

NCBI:

National Centre for Biotechnology Information

PBS:

Phosphate buffer solution

PRIDE:

PRoteomics IDEntifications

QC:

Quality control

qPCR:

Real-time Quantitative polymerase chain reaction

RIN:

RNA integrity value

RNA-seq:

RNA Sequencing

SRA:

Sequence Read Archive

TCA:

Tricarboxylic acid

TSB:

Tryptic soy broth

References

  1. Janda JM, Abbott SL. The genus Aeromonas: taxonomy, pathogenicity, and infection. Clin Microbiol Rev. 2010;23:35–73.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  2. Feelders RA, Vreugdenhil G, Eggermont AM, Kuiper-Kramer PA, van Eijk HG, Swaak AJ. Regulation of iron metabolism in the acute-phase response: interferon gamma and tumour necrosis factor alpha induce hypoferraemia, ferritin production and a decrease in circulating transferrin receptors in cancer patients. Eur J Clin Investig. 1998;28:520–7.

    Article  CAS  Google Scholar 

  3. Reines HD, Cook FV. Pneumonia and bacteremia due to Aeromonas hydrophila. Chest. 1981;80:264–7.

    Article  PubMed  CAS  Google Scholar 

  4. Brenden RA, Huizinga HW. Pathophysiology of experimental Aeromonas hydrophila infection in mice. J Med Microbiol. 1986;21:311–7.

    Article  PubMed  CAS  Google Scholar 

  5. Rasmussen-Ivey CR, Figueras MJ, Mcgarey D, Liles MR. Virulence factors of Aeromonas hydrophila: in the wake of reclassification. Front Microbiol. 2016;7:1337.

    PubMed  PubMed Central  Google Scholar 

  6. Toma'S JM. The main Aeromonas pathogenic factors. ISRN Microbiol. 2012;2012:256261.

  7. Mekalanos JJ. Environmental signals controlling expression of virulence determinants in bacteria. J Bacteriol. 1992;174:1–7.

  8. Sritharan M. Iron as a candidate in virulence and pathogenesis in mycobacteria and other microorganisms. World J Microbiol Biotechnol. 2000;16:769–80.

    Article  CAS  Google Scholar 

  9. Halliwell B, Gutteridge JM. Oxygen toxicity, oxygen radicals, transition metals and disease. Biochem J. 1984;219:1–14.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  10. Braun V. Avoidance of iron toxicity through regulation of bacterial iron transport. Biol Chem. 1997;378:779–86.

    PubMed  CAS  Google Scholar 

  11. Miller RA, Britigan BE. Role of oxidants in microbial pathophysiology. Clin Microbiol Rev. 1997;10:1–18.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  12. Zhang YC, Shen YY, Yan XH, Wang FD. Molecular mechanisms of mammalian iron homeostasis. Chin. J Cell Biol. 2011;33:1179–90.

    CAS  Google Scholar 

  13. Teng T, Xi BW, Xie J, Chen K, Pao X, Pan LK. Molecular cloning and expression analysis of Megalobrama amblycephala transferrin gene and effects of exposure to iron and infection with Aeromonas hydrophila. Fish Physiol Biochem. 2017;43:987–97.

  14. Telford JR, Raymond KN. Amonabactin: a family of novel siderophores from a pathogenic bacterium. J Biol Inorg Chem. 1997;2:750–61.

    Article  CAS  Google Scholar 

  15. Litwin CM, Calderwood SB. Role of iron in regulation of virulence genes. Clin Microbiol Rev. 1993;6:137–49.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  16. Ratledge C, Dover LG. Iron metabolism in pathogenic bacteria. Annu Rev Microbiol. 2003;54:881–941.

    Article  Google Scholar 

  17. Miethke M, Marahiel M. Siderophore-based iron acquisition and pathogen control. Microbiol Mol Biol Rev. 2007;71:413–51.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  18. Ramanan N, Wang Y. A high-affinity iron permease essential for candida albicans virulence. Science. 2000;288:1062–4.

    Article  PubMed  CAS  Google Scholar 

  19. Horsburgh MJ, Clements MO, Crossley H, Ingham E, Foster SJ. Perr controls oxidative stress resistance and iron storage proteins and is required for virulence in staphylococcus aureus. Infect Immun. 2001;69:3744.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  20. Velayudhan J, Hughes NJ, Mccolm AA, Bagshaw J, Clayton CL, Andrews SC, Kelly DJ. Iron acquisition and virulence in helicobacter pylori : a major role for feob, a high-affinity ferrous iron transporter. Mol Microbiol. 2000;37:274–86.

    Article  PubMed  CAS  Google Scholar 

  21. Deng K, Blick RJ, Liu W, Hansen EJ. Identification of Francisella tularensis genes affected by iron limitation. Infect Immun. 2006;74:4224–36.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  22. Lenco J, Hubálek M, Larsson P, Fucíková A, Brychta M, Macela A, Stulík J. Proteomics analysis of the Francisella tularensis LVS response to iron restriction: induction of the F. tularensis pathogenicity island proteins IglABC. FEMS Microbiol Lett. 2007;269:11–21.

    Article  PubMed  CAS  Google Scholar 

  23. Folsom JP, Parker AE, Carlson RP. Physiological and proteomic analysis of Escherichia coli iron-limited chemostat growth. J Bacteriol. 2014;196:2748–61.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  24. Yao Z, Wang Z, Sun L, Li W, Shi Y, Lin L, Lin WX, Lin XM. Quantitative proteomic analysis of cell envelope preparations under iron starvation stress in Aeromonas hydrophila. BMC Microbiol. 2016;16:161.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  25. Pang MD, Jiang JW, Xie X, Wu YF, Dong YH, Kwok AH, Zhang W, Yao HC, Lu CP, Leung FC, Liu YJ. Novel insights into the pathogenicity of epidemic Aeromonas hydrophila ST251 clones from comparative genomics. Sci Rep. 2015;5:9833.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Caliaperumal J, Wowk S, Jones S, Ma YL, Colbourne F. Bipyridine, an iron chelator, does not lessen intracerebral iron-induced damage or improve outcome after intracerebral hemorrhagic stroke in rats. Transl Stroke Res. 2013;4:719–28.

    Article  PubMed  CAS  Google Scholar 

  27. Alencar TD, Wilmart-Gonçalves TC, Vidal LS, Fortunato RS, Leitão AC, Lage C. Bipyridine (2,2′-dipyridyl) potentiates Escherichia coli lethality induced by nitrogen mustard mechlorethamine. Mutat Res. 2014;765:40.

    Article  PubMed  CAS  Google Scholar 

  28. Lee P, Tan KS. Effects of epigallocatechin gallate against Enterococcus faecalis biofilm and virulence. Arch Oral Biol. 2015;60:393.

    Article  PubMed  CAS  Google Scholar 

  29. GB/T 5009.90–2003. Determination of iron, magnesium and manganese in foods. Beijing: Standardization Administration of the People’s republic of China; 2003.

    Google Scholar 

  30. Wang XK, Yang RQ, Zhou YL, Gu ZX. A comparative transcriptome and proteomics analysis reveals the positive effect of supplementary Ca2+ on soybean sprout yield and nutritional qualities. J Proteome. 2016;143:161.

    Article  CAS  Google Scholar 

  31. Qin N, Tan X, Jiao Y, Liu L, Zhao W, Yang S, Jia AQ. RNA-Seq-based transcriptome analysis of methicillin-resistant Staphylococcus aureus biofilm inhibition by ursolic acid and resveratrol. Sci Rep. 2014;4:5467.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  32. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215:403–10.

    Article  PubMed  CAS  Google Scholar 

  33. Robinson MD, Mccarthy DJ, Smyth GK. Edger: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26:139–40.

    Article  PubMed  CAS  Google Scholar 

  34. Yan MX, Dai WJ, Cai EP, Deng YZ, Chang CQ, Jiang ZD, Zhang LH. Transcriptome analysis of Sporisorium scitamineum reveals critical environmental signals for fungal sexual mating and filamentous growth. BMC Genomics. 2016;17:354.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  35. Kanehisa M, Araki M, Goto S, Hattori M, Hirakawa M, Itoh M, Katayama T, Kawashima S, Okuda S, Tokimatsu T, Yamanishi Y. KEGG for linking genomes to life and the environment. Nucleic Acids Res. 2008;36:480–4.

    Article  CAS  Google Scholar 

  36. Isaacson T, Damasceno CM, Saravanan RS, He Y, Catalá C, Saladié M, Rose JKC. Sample extraction techniques for enhanced proteomic analysis of plant tissues. Nat Protoc. 2006;1:769–74.

    Article  PubMed  CAS  Google Scholar 

  37. Smith PK, Krohn RIG, Hermanson G, Mallia AKFD, Gartner FJH, Provenzano MD, Fujimoto EK, Goeke NM, Olson BJ, Klenk DC. Measurement of protein using Bicinchoninic acid. Anal Biochem. 1985;150:76–85.

    Article  PubMed  CAS  Google Scholar 

  38. Candiano G, Bruschi M, Musante L, Santucci L, Ghiggeri GM, Carnemolla B, Orecchia P, Zardi L, Righetti PG. Blue silver: a very sensitive colloidal Coomassie G-250 staining for proteome analysis. Electrophoresis. 2004;25:1327–33.

    Article  PubMed  CAS  Google Scholar 

  39. Wisniewski J, Zougman A, Nagaraj N, Mann M. Universal sample preparation method for proteome analysis. Nat Methods. 2009;6:359–62.

    Article  PubMed  CAS  Google Scholar 

  40. You C, Lin C. He H, et al. iTRAQ-based proteome profile analysis of superior and inferior Spikelets at early grain filling stage in japonica Rice. BMC Plant Biol. 2017;17(1):100.

    Article  PubMed  PubMed Central  Google Scholar 

  41. Shilov IV, Seymour SL, Patel AA, Loboda A, Tang WH, Keating SP, Hunter CL, Nuwaysir LM, Schaeffer DA. The paragon algorithm, a next generation search engine that uses sequence temperature values and feature probabilities to identify peptides from tandem mass spectra. Mol Cell Proteomics. 2007;6:1638–55.

    Article  PubMed  CAS  Google Scholar 

  42. Zhang MC, Cao YN, Yao B, Bai DQ, Zhou ZG. Characteristics of quenching enzyme AiiO-AIO6 and its effect on Aeromonas hydrophila virulence factors expression. J Fish China. 2011;35:1720–8.

    CAS  Google Scholar 

  43. Swift S, Lynch MJ, Fish L, Kirke DF, Tomás JM, Stewart GSAB, Williams P. Quorum sensing-dependent regulation and blockade of exoprotease production in Aeromonas hydrophila. Infect Immun. 1999;67:5192–9.

    PubMed  PubMed Central  CAS  Google Scholar 

  44. Chu W, Zhou S, Zhu W, Zhuang X. Quorum quenching bacteria Bacillus sp. QSI-1 protect zebrafish (Danio rerio) from Aeromonas hydrophila infection. Sci Rep. 2014;4:5446.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  45. Gang L, Huang L, Su Y, Qin Y, Xu X, Zhao L, Yan Q. Flra, flrb and flrc regulate adhesion by controlling the expression of critical virulence genes in Vibrio alginolyticus. Emerging Microbes Infect. 2016;5:e85.

    Article  CAS  Google Scholar 

  46. Tsou AM, Zhu J. Quorum sensing negatively regulates hemolysin transcriptionally and posttranslationally in Vibrio cholerae. Infect Immun. 2010;78:461–7.

    Article  PubMed  CAS  Google Scholar 

  47. Teng T, Liang LG, Chen K, Xi BW, Xie J, Xu P. Isolation, identification and phenotypic and molecular characterization of pathogenic Vibrio vulnificus isolated from Litopenaeus vannamei. PLoS One. 2017;12:e0186135.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  48. Saganuwan SA. A modified arithmetical method of reed and Muench for determination of a relatively ideal median lethal dose (LD 50). Afr J Pharm Pharmacol. 2011;5:1543–6.

    Article  Google Scholar 

  49. Liu W, Dong H, Li J, Ou Q, Lv Y, Wang X, Xiang Z, He Y, Wu Q. RNA-seq reveals the critical role of OtpR in regulating Brucella melitensis metabolism and virulence under acidic stress. Sci Rep. 2015;5:10864.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  50. Nagar V, Bandekar JR, Shashidhar R. Expression of virulence and stress response genes in Aeromonas hydrophila under various stress conditions. J Basic Microbiol. 2016;56:1132–7.

    Article  PubMed  CAS  Google Scholar 

  51. Teixeiragomes AP, Cloeckaert A, Zygmunt MS. Characterization of heat, oxidative, and acid stress responses in Brucella melitensis. Infect Immun. 2000;68:2954–61.

    Article  CAS  Google Scholar 

  52. Casabianca A, Orlandi C, Barbieri F, Sabatini L, Cesare AD, Sisti D, Pasquaroli S, Magnani M, Citterio B. Effect of starvation on survival and virulence expression of Aeromonas hydrophila from different sources. Arch Microbiol. 2015;197:431–8.

    Article  PubMed  CAS  Google Scholar 

  53. Payne SM, Lawlor KM. Chapter 11 – molecular studies on iron acquisition by non- Escherichia coli, species. Bacteria. 1990:225–48.

  54. Allan BJ, Stevenson RM. Extracellular virulence factors of Aeromonas hydrophila in fish infections. Can J Microbiol. 1981;27:1114–22.

    Article  PubMed  CAS  Google Scholar 

  55. Lv J, Yu GC, Sun ZH, Wang NJ, Zhu Y, Wang HC, Sun XS. Proteomic analysis of effects of iron depletion on Streptococcus pyogenes MGAS5005. Microbiology China. 2012;39:515–25.

    Google Scholar 

  56. Bhatnagar, Elkins, Fortier. Heat stress alters the virulence of a rifampin-resistant mutant of Francisella tularensis LVS. I nfect Immun 1995; 63:154–159.

  57. Ellis AE, Burrows AS, Stapleton KJ. Lack of relationship between virulence of Aeromonas salmonicida and the putative virulence factors: A-layer, extracellular proteases and extracellular haemolysins. J Fish Dis. 1988;11:309–23.

    Article  CAS  Google Scholar 

  58. Brickman TJ, Cummings CA, Liew SY, Relman DA, Armstrong SK. Transcriptional profiling of the iron starvation response in Bordetella pertussis provides new insights into siderophore utilization and virulence gene expression. J Bacteriol. 2011;193:4798–812.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  59. Sheikh MA, Taylor GL. Crystal structure of the Vibrio cholerae, ferric uptake regulator (Fur) reveals insights into metal co-ordination. Mol Microbiol. 2009;72:1208–20.

    Article  PubMed  CAS  Google Scholar 

  60. Girija R. Iron and virulence in Francisella tularensis. Front Cell Infect Microbiol. 2017;7:107.

    Google Scholar 

  61. Skaar EP. The battle for Iron between bacterial pathogens and their vertebrate hosts. PLoS Pathog. 2010;6:e1000949.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  62. Wrighting DM, Andrews NC. Iron homeostasis and erythropoiesis. Curr Top Dev Biol. 2008;82:141–67.

    Article  PubMed  CAS  Google Scholar 

  63. Maltz M, Levarge BL, Graf J. Identification of iron and heme utilization genes in Aeromonas and their role in the colonization of the leech digestive tract. Front Microbiol. 2015;6:763.

    Article  PubMed  PubMed Central  Google Scholar 

  64. Hayrapetyan H, Siezen R, Abee T, Groot MN. Comparative genomics of Iron-transporting systems in Bacillus cereus strains and impact of Iron sources on growth and biofilm formation. Front Microbiol. 2016;7:842.

    Article  PubMed  PubMed Central  Google Scholar 

  65. Neilands JB. Siderophores: structure and function of microbial iron transport compounds. J Biol Chem. 1995;270:26723–6.

    Article  PubMed  CAS  Google Scholar 

  66. Long H, Zeng Y. Studies on resistance property of fish serum Transferrins against Aeromonas hydrophila. Journal of Hubei Agricultural College. 2004;24:119–23.

    Google Scholar 

  67. Gehring AM, Mori I, Walsh CT. Reconstitution and characterization of the Escherichia coli Enterobactin Synthetase from EntB, EntE, and EntF. Biochemistry. 1998;37:2648–59.

    Article  PubMed  CAS  Google Scholar 

  68. Neilands JB. Molecular aspects of regulation of high affinity iron absorption in microorganisms. Adv Inorg Biochem. 1990;8:63–90.

    PubMed  CAS  Google Scholar 

  69. Crosa JH, Walsh CT. Genetics and assembly line enzymology of Siderophore biosynthesis in Bacteria. Microbiol Mol Biol Rev. 2002;66:223–49.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  70. Franza T, Enard C, Van GF, Expert D. Genetic analysis of the Erwinia chrysanthemi 3937 chrysobactin iron-transport system: characterization of a gene cluster involved in uptake and biosynthetic pathways. Mol Microbiol. 1991;5:1319–29.

    Article  PubMed  CAS  Google Scholar 

  71. Bearson BL, Bearson SM, Uthe JJ, Dowd SE, Houghton JO, Lee I, Toscano MJ, Lay Jr DC. Iron regulated genes of Salmonella enterica serovar typhimurium in response to norepinephrine and the requirement of fepDGC for norepinephrine-enhanced growth. Microbes Infect. 2008;10:807–16.

    Article  PubMed  CAS  Google Scholar 

  72. Bou-Abdallah F. The iron redox and hydrolysis chemistry of the ferritins. Biochim Biophys Acta. 2010;1800:719–31.

    Article  PubMed  CAS  Google Scholar 

  73. Hubálek M, Hernychová L, Havlasová J, Kasalová I, Neubauerová V, Stulík J, Macela A, Lundqvist M, Larsson P. Towards proteome database of Francisella tularensis. J Chromatogr B. 2003;787:149–77.

    Article  CAS  Google Scholar 

  74. Hubálek M, Hernychová L, Brychta M, Lenco J, Zechovská J, Stulík J. Comparative proteome analysis of cellular proteins extracted from highly virulent Francisella tularensis ssp. tularensis and less virulent F. tularensis ssp. holarctica and F. tularensis ssp. mediaasiatica. Proteomics. 2004;4:3048–60.

    Article  PubMed  CAS  Google Scholar 

  75. Holden KM, Browning GF, Noormohammadi AH, Markham PF, Marenda MS. TonB is essential for virulence in avian pathogenic Escherichia coli. Comparative Immunology Microbiology and Infectious Diseases. 2012;35:129–38.

    Article  Google Scholar 

  76. Dong Y, Liu J, Pang M, Du H, Wang N, Awan F, Lu C, Liu Y. Catecholamine-stimulated growth of Aeromonas hydrophila requires the TonB2 energy transduction system but is independent of the Amonabactin Siderophore. Front Cell Infect Microbiol. 2016;6:183.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  77. Angerer A, Gaisser S, Braun V. Nucleotide sequences of the sfuA, sfuB, and sfuC genes of Serratia marcescens suggest a periplasmic-binding-protein-dependent iron transport mechanism. J Bacteriol. 1990;172:572–8.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  78. Andrews NC. Iron homeostasis: insights from genetics and animal models. Nat Rev Genet. 2000;1:208–17.

    Article  PubMed  CAS  Google Scholar 

  79. Rahman MH, Suzuki S, Kawai K. Formation of viable but non-culturable state (VBNC) of Aeromonas hydrophila, and its virulence in goldfish, Carassius auratus. Microbiol Res. 2001;156:103–6.

    Article  PubMed  CAS  Google Scholar 

  80. Maalej S, Gdoura R, Dukan S, Hammami A, Bouain A. Maintenance of pathogenicity during entry into and resuscitation from viable but nonculturable state in Aeromonas hydrophila, exposed to natural seawater at low temperature. J Appl Microbiol. 2004;97:557–65.

    Article  PubMed  CAS  Google Scholar 

  81. Mchugh JP, Rodríguezquinoñes F, Abdultehrani H, Svistunenko DA, Poole RK, Cooper CE, Andrews SC. Global iron-dependent gene regulation in Escherichia coli. A new mechanism for iron homeostasis J Biol Chem. 2003;278:29478–86.

    PubMed  CAS  Google Scholar 

  82. He X. Proteomic analysis of biological affects on Streptococcus pneumoniae induced by manganese depression. Jinan University 2011.

  83. Cunningham AB, Sharp RR, Jr FC, Gerlach R. Effects of starvation on bacterial transport through porous media. Adv Water Resour. 2007;30:1583–92.

    Article  Google Scholar 

  84. Elgaml A, Miyoshi SI. Regulation systems of protease and hemolysin production in Vibrio vulnificus. Microbiol Immunol. 2017;61:1–11.

    Article  PubMed  CAS  Google Scholar 

  85. Tang TS, Lu CP. An Acinetobacter baumannii strain isolated from mandarin fish possesses type 4 pili. J Nanjing Agric Univ. 1997;20:114–6.

    Google Scholar 

  86. Chen XQ, Cai HY, Zhang W, Yan MY, Gao H, Duan GX, Yan ZY. The Role of Fur Involved in Bioflim Formation of Vibrio Cholerae. Prog Mod Biomed. 2013;13:2841–5.

  87. Vasil ML, Ochsner UA. The response of Pseudomonas aeruginosa, to iron: genetics, biochemistry and virulence. Mol Microbiol. 1999;34:399–413.

    Article  PubMed  CAS  Google Scholar 

  88. Thune RL, Johnson MC, Graham TE, Amborski RL. Aeromonas hydrophila B-haemolysin: purification and examination of its role in virulence in 0-group channel catfish, Ictalurus punctatus (Rafinesque). J Fish Dis. 1986;9:55–61.

    Article  CAS  Google Scholar 

  89. Zhu DL, Li AH, Wang JG, Li M, Cai TZ, Hu G. The correlation between the distribution pattern of virulence genes and the virulence of Aeromonas hydrophila strains. Acta Sci Nat Univ Sunyatseni. 2006;45:82–5.

    CAS  Google Scholar 

  90. Chu WH. Invasion mechanism and extracellular protease of Aeromonas hydrophila. Nanjing Agric Univ. 2002;

  91. Cao Q, Zhang XY, Pang MD, Wang NN. Identification and molecular typing of the epidemic Aeromonas hydrophila strains in one farm of Nanjing. J Fish China. 2017;41:134–41.

    CAS  Google Scholar 

  92. Cascón A, Yugueros J, Temprano A, Sánchez M, Hernanz C, Luengo JM, Naharro GN. A major secreted elastase is essential for pathogenicity of Aeromonas hydrophila. Infect Immun. 2000;68:3233–41.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Funding

This study was supported by The earmarked fund for China Agriculture Research System (CARS-45), Natural Science Foundation of China (31572662), Jiangsu Natural Science Foundation (BK20171152), and Postgraduate Research & Practice Innovation Program of Jiangsu Province (KYLX16_1082). The funding agencies have not been involved in the design of research and collection, analysis and interpretation of data and in writing of manuscripts. We also thank OEbiotech.co.ltd for technical support in our transcriptome and proteomics analyses.

Availability of data and materials

The RNA-seq data and analyses discussed in this publication were deposited in the NCBI Sequence Read Archive (SRA) database under accession number SRR5894319. The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE partner repository with the dataset identifier PXD007641.

Author information

Author notes

  1. Tao Teng and Bingwen Xi contributed equally to this work.

Authors and Affiliations

  1. Wuxi Fisheries College, Nanjing Agricultural University, Wuxi, 214081, China

    Tao Teng, Bingwen Xi, Jun Xie & Pao Xu

  2. Key Laboratory of Freshwater Fisheries and Germplasm Resources Utilization, Ministry of Agriculture, Freshwater Fisheries Research Center, Chinese Academy of Fishery Sciences, Wuxi, 214081, China

    Tao Teng, Bingwen Xi, Kai Chen, Liangkun Pan, Jun Xie & Pao Xu

Authors
  1. Tao Teng
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Bingwen Xi
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Kai Chen
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Liangkun Pan
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Jun Xie
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Pao Xu
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

BX, JX and PX conceived and designed the experiments; BX guided the experiments, TT performed the experiments, analyzed the data; TT and BX wrote the paper, revised the paper, they contributed equally to this work; LP, KC participated in the collection of samples, planning and coordination of the study, provided general supervision. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Jun Xie or Pao Xu.

Ethics declarations

Ethics approval and consent to participate

The study protocol was granted by the Research Ethics Committee, Wuxi Fisheries College of Nanjing Agriculture University (Permit No. NJYY20160929–1), and all methods were performed in accordance with the approved guidelines and regulations.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

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

Additional files

Additional file 1:

Excel S1. The results of differentially expressed transcripts and proteins to find the corresponding genes and proteins. (XLSX 87 kb)

Additional file 2:

Excel S2. The enriched DEGs GO terms between control and iron-limitation groups using bioinformatics methods. (XLSX 42 kb)

Additional file 3:

Excel S3. The enriched DEPs GO terms between control and iron-limitation groups using bioinformatics methods. (XLSX 42 kb)

Additional file 4:

Excel S4. The enriched DEGs KEGG terms between control and iron-limitation groups using bioinformatics methods. (XLSX 25 kb)

Additional file 5:

Excel S5. The enriched DEPs KEGG terms between control and iron-limitation groups using bioinformatics methods. (XLSX 13 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), 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 (http://creativecommons.org/publicdomain/zero/1.0/) 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

Teng, T., Xi, B., Chen, K. et al. Comparative transcriptomic and proteomic analyses reveal upregulated expression of virulence and iron transport factors of Aeromonas hydrophila under iron limitation. BMC Microbiol 18, 52 (2018). https://doi.org/10.1186/s12866-018-1178-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12866-018-1178-8

Keywords

  • Transcriptomic
  • Proteomic
  • Virulence
  • Iron
  • Aeromonas hydrophila

BMC Microbiology

ISSN: 1471-2180

Contact us