ORIGINAL RESEARCH article
Transcriptome Profiling of Pacu (Piaractus mesopotamicus) Challenged With Pathogenic Aeromonas hydrophila: Inference on Immune Gene Response
- 1Aquaculture Center, São Paulo State University (Unesp), Jaboticabal, Brazil
- 2Faculty of Agricultural and Veterinary Sciences, São Paulo State University (Unesp), Jaboticabal, Brazil
- 3The Royal (Dick) School of Veterinary Studies, The University of Edinburgh, Edinburgh, United Kingdom
Pacu (Piaractus mesopotamicus) is a Neotropical fish of major importance for South American aquaculture. Septicemia caused by Aeromonas hydrophila bacteria is currently considered a substantial threat for pacu aquaculture that have provoked infectious disease outbreaks with high economic losses. The understanding of molecular aspects on progress of A. hydrophila infection and pacu immune response is scarce, which have limited the development of genomic selection for resistance to this infection. The present study aimed to generate information on transcriptome of pacu in face of A. hydrophila infection, and compare the transcriptomic responses between two groups of time-series belonging to a disease resistance challenge, peak mortality (HM) and mortality plateau (PM) groups of individuals. Nine RNA sequencing (RNA-Seq) libraries were prepared from liver tissue of challenged individuals, generating ∼160 million 150 bp pair-end reads. After quality trimming/cleanup, these reads were assembled de novo generating 211,259 contigs. When the expression of genes from individuals of HM group were compared to individuals from control group, a total of 4,413 differentially expressed transcripts were found (2,000 upregulated and 2,413 downregulated candidate genes). Additionally, 433 transcripts were differentially expressed when individuals from MP group were compared with those in the control group (155 upregulated and 278 downregulated candidate genes). The resulting differentially expressed transcripts were clustered into the following functional categories: cytokines and signaling, epithelial protection, antigen processing and presentation, apoptosis, phagocytosis, complement system cascades and pattern recognition receptors. The proposed results revealing relevant differential gene expression on HM and PM groups which will contribute to a better understanding of the molecular defense mechanisms during A. hydrophila infection.
Disease outbreaks are considered a primary constraint and a continuous challenge to the growth of aquaculture sector (FAO, 2018). Detrimental effects due improper management practices, high-density culture systems and fluctuating water parameters have facilitated the spreading of bacteria or other pathogens, resulting in frequent and often severe disease outbreaks. This scenario has led to the billions of dollars in losses to aquaculture every year, that may exceed 40% of global capacity of production (Lafferty et al., 2015; Stentiford et al., 2017; Owens, 2019). While economic impact of disease outbreaks are under-reported for fish production in South America, a quarter of Brazilian large producers have reported moderate to high production losses due to disease outbreaks in the major aquaculture regions (Benavides and Maciel, 2015).
The incidence of sanitary problems linked to intensification of fish production from South America have resulted in large mortalities of Pacu (Piaractus mesopotamicus), one of the most important native fish species produced in Brazil (IBGE, 2016; Woynárovich and Van Anrooy, 2019). Pacu is considered a promising species for Brazilian aquaculture due to its large size, high quality meat and low-cost diet (Woynárovich and Van Anrooy, 2019). However, one of the main causes of disease associated with large mortalities for this species is the infection by Aeromonas hydrophila (Carraschi et al., 2012; Farias et al., 2016; Mastrochirico-Filho et al., 2019).
Aeromonas hydrophila is considered the etiological agent of Motile Aeromonas Septicemia (MAS), and is generally associated with symptoms such as reddened or rotten fins, small ulcerative lesions on the skin and hemorrhage, often followed by mortality of individuals within hours after manifestation of disease (Silva et al., 2012; Stratev and Odeyemi, 2017). Despite the lack of official data on the economic losses in fish production attributed to A. hydrophila infection, Brazilian fish farmers have reported that MAS prevalence is responsible for 20–30% of annual mortalities on pacu farming. However, insufficient knowledge on preventive measures and immune mechanisms of susceptibility and resistance are available to MAS (Belém-Costa and Cyrino, 2006; Farias et al., 2016; Zanuzzo et al., 2017; Mastrochirico-Filho et al., 2019).
RNA-Seq has been considered an effective technique for understanding the transcriptomic changes that occur during infection, exploring the activation of defense mechanisms involved in fast responses of infected hosts, even in organisms without a reference genome (Sudhagar et al., 2018). Over the last years, an increasing effort has been accomplished to characterize host-pathogen interaction by differentially expressed genes (DEGs) of fish species submitted to A. hydrophila resistance challenges, such as blunt snout bream Megalobrama amblycephala (Tran et al., 2015); grass carp Ctenopharyngodon idella (Dang et al., 2016); darkbarbel catfish Pelteobagrus vachellii (Qin et al., 2017); hybrid sturgeon (Huso dauricus × Acipenser schrenckii) (Jiang et al., 2018); and other species (Maekawa et al., 2019). These expression profiling have been considered useful to obtain information about the genomic makeup of immune responses at a molecular level, contributing to a better understanding of genetic resistance to disease, and guiding selective breeding strategies to produce more resistant strains. Additionally, in order to develop knowledge to suppress the disease, a better understanding of host-pathogen interaction between fish and bacteria, during critical periods of a disease outbreak is needed.
Therefore, the aim of the study was to characterize the transcriptome of pacu after A. hydrophila infection. A de novo assembly was necessary, since pacu is considered a non-model species and does not have a reference genome until the present moment. Further, the identification of differentially expressed genes involved in the response to this bacteria was performed aiming immune system genes, comparing differences between the actions of genomic mechanisms under two progressive time trends of a disease resistance challenge, peak mortality (HM) and mortality plateau (PM).
Materials and Methods
This study was conducted in strict accordance with the recommendations of the National Council for Control of Animal Experimentation (CONCEA) (Brazilian Ministry for Science, Technology and Innovation) and was approved by the Ethics Committee on Animal Use (CEUA number 19.005/17) of Faculdade de Ciências Agrárias e Veterinárias, UNESP, Campus Jaboticabal, SP, Brazil. All biosafety procedures related to the management of bacteria and infected animals were carried out. In addition, all measures related to the safety of the participants of the present study in handling the strains of bacteria were ensured.
Generation of Pacu Families
Individuals used for RNA sequencing were obtained from an experimental challenge carried out in 14 full-sibling families of pacu, generated by a hierarchical mating scheme using 4 dams and 14 sires (3 or 4 sires per dam). Induced spawning was performed using carp pituitary extract dissolved in saline solution (0.9% NaCl) and applied in two dosages, with a 12 h interval (first and second dosage of 0.6 and 5.4 mg/kg, respectively). For sires, a single dosage was used, at the same time of the second dosage for dams, equivalent to 1.5 mg/kg of carp pituitary extract.
After hatching in 20 l conical fiberglass incubators, the larvae were fed with artemia nauplii for 20 days. Gradually, the feed was replaced by 50% of crude protein. In the fingerling stage, 1.2 mm pelleted feeds were used (40% of crude protein), being gradually replaced by 2 to 3 mm pelleted feeds (36% of crude protein) provided twice daily in 60 l tanks. Animals used in the experiment were pit-tagged to maintain the pedigree information during the challenge experiments. Laterally, fish were kept in 800 l fiberglass tanks at the Laboratory of Genetics in Aquaculture and Conservation (LaGeAC), at the São Paulo State University (UNESP), Jaboticabal city (São Paulo, Brazil).
Experimental Challenge by Aeromonas hydrophila
Bacterial challenge was performed using a strain of A. hydrophila isolated from an outbreak of pacu aeromoniasis in a commercial facility from the São Paulo State, by the Laboratory of Microbiology and Parasitology of Aquatic Organisms, at UNESP, Jaboticabal city (São Paulo, Brazil). The strain of A. hydrophila was the same of that used by Mastrochirico-Filho et al. (2019). The strain was cultured in Agar Trypticase Soy (ATS), Vegitone (Sigma-Aldrich) for 24 h (28°C). The colony was then transferred to a nutrient tryptic soy broth (TSB) (Sigma-Aldrich) and cultured for 24 h (28°C). After bacteria growth, the culture was centrifuged at 5,000 × g for 10 min (4°C, in the Eppendorf Centrifuge 5810), forming a bacterial pellet suspended in a saline solution (PBS) and washed twice.
The challenge test performed in this study was carried out by intraperitoneal inoculation, according to protocols of Mastrochirico-Filho et al. (2019). The LD50 (lethal dose in 50% of individuals) was previously tested in 60 randomly chosen individuals from the same pacu families using concentrations adjusted by optical density of the solution at 0.400, 0.600, and 0.800 at 625 nm in spectrophotometer (2100 Unico, Japan). A sample of 100 μl of the LD50 was removed from the inoculum to perform serial dilutions and plate counts in duplicate on Trypticase Soy Agar (TSA).
Prior to challenge, the presence of A. hydrophila and other pathogens, such as Flavobacterium columnare and Streptococcus agalactiae in the populations was discarded checking subsamples of the populations by routine microbiological tests. In the experimental design of the challenge test, the individuals were distributed into three communal tanks of 2 m3 (length = 2 m, width = 1 m, depth = 1 m), where approximately ten individuals from each family were randomly distributed into each treatment tank. Prior to the inoculation of bacteria, fish were anesthetized with benzocaine (0.1 mg/l) and weighed. The mean weight of animals prior to the bacterial inoculation time was 23.0 ± 9.06 g. Individual fish were injected by intraperitoneal inoculation of the predefined LD50 dose of live cells of A. hydrophila (8 × 105 CFU/g body weight), according to protocols carried out by Mastrochirico-Filho et al. (2019). Moreover, approximately ten fish from each family were also used as control and kept in a separated tank (called as control tank) of 2 m3 (length = 2 m, width = 1 m, depth = 1 m). Individuals of the control tank were injected by intraperitoneal inoculation of saline solution (PBS).
Each treatment and control tank were maintained with an independent water recirculation system, fitted with mechanical and biological filters, external aeration system, and controlled temperature using thermal controller connected to heaters (2 × 500 w). Water quality parameters were determined daily during 14 days of the challenge. Temperature, dissolved oxygen and pH were measured with a Multiparameter Water Quality Checker U-50 (Horiba, Kyoto, Japan). Water temperature was maintained at 30°C (SD = 0.5), with dissolved oxygen at 5.8 mg/L (SD = 1.0) and pH at 7.0 (SD = 1.1) during the challenge period. No water was exchanged during the challenge, but the tanks were topped off to compensate for evaporation. The control tank had the same condition when compared with treatment tanks, except for an antibacterial treatment through a UV filter.
Fish mortality was observed throughout the entire day (24 h) in the initial 3 days of challenge, and in intervals of 8 h in the remaining days of the challenge experiment. Fish that showed mortality by clinical signs of A. hydrophila infection (e.g., disequilibrium, hemorrhage, isolation from the group) were recorded and removed immediately from the tanks. Necropsy examination and microbiological tests were performed in a sub-sample of dead fish in order to confirm mortality by A. hydrophila and discard other pathogens. The isolation of bacteria was performed in a specific growth medium (phenol red agar and ampicillin) incubated at 28°C for 48 h. All surviving fish were examined externally for clinical signs of disease.
Tissue Collection, RNA Isolation and Sequencing
According to experimental challenges by A. hydrophila in different fish species, the response of the immune system is very fast to this bacteria infection, occurring when mortality rates are generally detected in the first 24 h after exposure, whereas the mortality plateau generally occurs 72 h after infection (Mahapatra et al., 2008; Ødegård et al., 2010; Xiong et al., 2017; Mastrochirico-Filho et al., 2019). Therefore, fish presenting classical symptoms of A. hydrophila infection at 24 h after bacteria inoculation were collected and considered as belonging to HM (period of high mortality) group. When the mortality plateau occurred at 72 h after bacteria inoculation, fish were also collected and considered as belonging to PM (period of mortality plateau) group.
Laterally, three randomly selected families were sampled for RNA sequencing. Therefore, liver tissue from three individuals belonging to the HM group and three individuals belonging to the PM group were sampled from each family. Three control fish at each time trends of mortality were collected from the same families. Liver was collected from all individuals and stored in RNAlater®. This organ was chosen because it plays a critical role in the coordination of innate and adaptive immune responses (Kumar et al., 2017; Rebl and Goldammer, 2018), and it has been used in several studies involving transcriptome analysis in fish for a better understanding of the behavior of the immune system in face of bacterial infections (Tran et al., 2015; Maekawa et al., 2019; Jiang et al., 2020).
RNA was extracted using RNeasy Plus Universal Tissue Kit (Qiagen, Valencia, CA, United States) following the manufacturer’s specifications. After extraction, total RNA was treated with DNAse I to remove any DNA contamination. Total RNA was quantified by Qubit RNA Assay Kit in a Qubit 2.0 Fluorometer (Life Technologies, Carlsbad, CA, United States) and RNA integrity verified using a 2100 Bioanalyzer (Agilent Technologies). After mRNA extraction by poly-A tail selection, cDNA libraries were constructed using the TruSeq RNA Library Prep Kit protocol (Illumina). Six cDNA libraries were created from the combination of three pooled samples per family (family A, V, and R) diluted to the same concentration, for each time-series of mortality (HM and PM). At the same time, three cDNA libraries (one per family) were created for control samples, pooling three samples belonging to time-series (HM + PM), totalizing 9 cDNA libraries (Table 1). Sequencing was performed on the Illumina HiseqTM 2500 platform.
Table 1. Preparation of cDNA pooled libraries of P. mesopotamicus families (A, V, R) challenged on A. hydrophila infection.
Data Processing and Analysis
All RNA-Seq data generated for this study were deposited in the National Center for Biotechnology Information Sequence Read Archive under BioProject ID PRJNA632934. Bioinformatics pipelines are available in Supplementary Table S1. The quality of raw data was initially checked using FASTQC v. 0.11.9 (Andrews, 2010). Quality filtering and removal of residual Illumina adaptors were conducted on read pairs using TRIMMOMATIC v.0.36 (Bolger et al., 2014). Leading and trailing bases with a Phred score less than 20 were removed. Additionally, the read was trimmed if a sliding window average Phred score over four bases was less than 20. Following other trimming strategies adopted on studies approaching transcriptomic changes in fish submitted to disease infections, only pair reads greater than 36 bp were maintained (Robledo et al., 2014, 2018).
After filtering out low quality reads, the remained reads were used to assemble a de novo transcriptome using TRINITY v.2.9.1, adopting default parameters (Grabherr et al., 2011) as kmer-length value of 25 and contiguous sequences (contigs) longer than 200 bp. Read representation was evaluated by mapping the cleaned reads back to generated assembly with BOWTIE2 v.220.127.116.11 (Langmead and Salzberg, 2012). A. hydrophila genome (Yang et al., 2016), zebrafish ribosomal RNA (rRNA) sequences and pacu mitochondrial genome (da Pimentel et al., 2014) were then downloaded from the NCBI (National Center for Biotechnology Information) database and mapped against the assembly using BOWTIE2 v.18.104.22.168 (Langmead and Salzberg, 2012) to remove any contamination. Redundancy was reduced using CD-HIT-EST v.4.6.8 software by clustering sequences at 95% identity (Li and Godzik, 2006). Transcriptome statistics were generated using TRANSRATE v1.0.2 (Smith-Unna et al., 2016).
Annotation and homology searches of assembled transcripts were conducted using the TRINOTATE pipeline (Bryant et al., 2017). The assembled transcripts were translated into long open reading frames by TRANSDECODER (Haas et al., 2013). Homology searches were performed against the UniprotKB/Swiss-Prot database (BLASTX - E-value 1e–3) to identify known protein sequences (Altschul et al., 1990). Gene Ontology (GO) (Ashburner et al., 2000), eggNOG/COG (Non-supervised orthologous groups) (Powell et al., 2012), and Kyoto Encyclopedia of Genes and Genomes (KEGG) (Kanehisa et al., 2012) databases were used to assign functions and biology pathways to transcripts.
KALLISTO v0.43.1 software was used to compute expression levels for each pooled samples (HM, PM and control) by pseudoalignments (Bray et al., 2016). Initially, an index was built with default k-mer size based on the contigs assembled by TRINITY. After, in order to quantify abundances of transcripts, TPM (transcripts per million reads) expression levels of genes for each of the nine pooled samples were computed. DESeq2 (Love et al., 2014) used negative binomial distribution, data normalization and FDR (false recovery rate) correction through Benjamini-Hochberg (BH) method, to identify genes differentially expressed comparing the gene expression from pooled samples belonging to HM and PM groups against pooled samples from control groups. To control false discovery rates of multiple-hypothesis testing, the p-values calculated by DESeq2 were processes to q-values. Then, the list of differentially expressed transcripts were generated with cutoff q-values less than 0.05 and log2FC (fold change (condition/control) for a gene) > 1 or log2FC < −1) as described in similar transcriptome studies (Marancik et al., 2015; Polinski et al., 2016; Jiang et al., 2020).
Disease Challenge Experiment
During the A. hydrophila resistance challenge, moribund individuals demonstrated lethargic behavior, erratic swimming and red spotted skin lesions on the operculum, head, fins and gills. No symptoms related to infection were reported in individuals in the unchallenged control group. High mortality rates (93.1%) were observed mainly at the first day after inoculation (HM period). Resistance rates of sampled families for RNA sequencing ranged from 16 to 32%. A. hydrophila was successfully isolated from infected individuals, confirming the bacterial infection. The last mortality was recorded at approximately 65 h post inoculation (PM period).
Sequencing and de novo Assembly
RNA sequencing, generated from liver samples, yielded between 15 and 20 million paired-end reads per pooled library. Approximately 12% of low-quality reads (quality score < 20) were removed. This resulted in approximately 143 million (87.8%) high quality paired-end reads, with a total length of 39.81 × 109 bp to be explored in the analysis of the pacu transcriptome (Table 2). De novo assembly produced 238,054 contigs longer than 200 bp, totalizing 169,652,373 assembled bases. Additionally, the average contig length was 712.6 bp with an N50 value of 1,296 bp. Approximately 92% of the filtered reads mapped back to the assembly, of which 76% mapped adequately (both pair reads aligned in the correct orientation and in the same contig) (Table 3).
Table 2. RNA-seq statistics in P. mesopotamicus families (A, V, R) challenged on A. hydrophila infection.
Table 3. Summary of de novo assembly statistics of P. mesopotamicus liver transcriptome in A. hydrophila disease challenge.
Posteriorly, 115 sequences mapped to the A. hydrophila genome, 20 sequences mapped to pacu mitochondrial genome and 90 transcripts identified as part of zebrafish ribosomal RNA (rRNA) were identified as contaminants, and therefore were excluded of the analysis. In addition, CD-HIT-EST clustered around 11.2% (26,570) of the transcripts resulting in a non-redundant assembly containing 211,259 sequences.
A total of 25,807 (12.2%) transcripts were annotated based on currently known proteins in the UniprotKB/Swiss-Prot database (Supplementary Table S2). Most of the annotate transcripts (15,547; 60%) showed significant homology with Astyanax mexicanus proteins, followed by Danio rerio, Cyprinus carpio, Clupea harengus, and Salmo salar proteins with 2,340 (9%); 1,368 (5.3%); 574 (2.2%); and 422 (1,6%) matched transcripts, respectively. Laterally, 23,420 transcripts were assigned into one or more GO terms classifying 20,851 transcripts into biological processes, 20,604 transcripts into molecular function and 19,884 transcripts into cellular component category (Figure 1). The most frequent GO term subcategories were cellular process, metabolic process and biological regulation for biological process; binding and catalytic activity for molecular function; and cell, cell part and organelle for cellular component. Furthermore, GO annotation revealed pertinent information on 2,726 transcripts assigned to immune system processes.
Figure 1. Functional classification of the liver transcriptome of pacu (P. mesopotamicus) exposed to A. hydrophila infection. The transcripts were classified in the x-axis, according to Gene Ontology (GO) categories at level 2: cellular component, molecular function and biological process (E-value < 1e-3). The primary y-axis represents the percentage of genes classified in each subcategory, whereas the secondary y-axis represents the number of classified genes.
Cluster of Orthologous Groups (COG) using the EggNOG database assigned 4,931 transcripts to 25 COG functional categories (Supplementary Table S3 and Figure 2). The largest clusters were composed by transcripts of unknown function (1,075 transcripts, 21.8%), followed by clusters related to signal transduction mechanisms (884 transcripts, 17.9%), post-translational modification and protein turnover (chaperones systems) (510 transcripts, 10.3%). Other clusters related to transcripts involved in intracellular trafficking secretion and vesicular transport (296 transcripts, 6.0%), RNA processing and modification (187 transcripts, 3.8%) and defense mechanisms (61 transcripts, 1.2%) were well represented. Further, 8,422 transcripts were mapped to 405 KEGG pathways (Supplementary Table S4). The most represented KEGG pathways are shown in Figure 3. “Signal Transduction” pathway had the largest number of matched transcripts (1,540, 18.3%), followed by “Global and Overview Maps” (1,374, 16.3%), “Infectious Disease (viral)” (814, 9.6%), “Endocrine System” (732, 8.7%), and “Immune System” (703, 8.3%) pathways.
Figure 2. Functional classification of the pacu liver transcriptome exposed to A. hydrophila infection based on eggnog/COG orthologous groups (E-value < 1e-3). 4,931 pacu transcripts were classified in 25 functional categories arranged in the x-axis. The y-axis presents the frequency of each represented category.
Figure 3. Number of pacu transcripts (x-axis) assigned to each KEGG pathway arranged in the y-axis. The categories GIP and EIP stand for “Genetic Information Processing” and “Environmental Information Processing.”
Differential Gene Expression Analysis
Differential expression analyses between pooled samples belonging to control and HM group showed 4,413 differentially expressed transcripts, where 2,000 transcripts were up-regulated (45.3%) and 2,413 transcripts were down-regulated (54.7%) in HM group (q-value < 0.05, log2FC > | 1|) (Figure 4A, Table 4, and Supplementary Table S5). Additionally, 433 transcripts were differentially expressed between PM group and control fish, with 155 up-regulated (35.8%) and 278 down-regulated (64.2%) transcripts considering the PM group (q-value < 0.05, log2FC > |1|) (Figure 4B, Table 4, and Supplementary Table S6).
Table 4. Number of differentially expressed transcripts of pacu (P. mesopotamicus) after infection with A. hydrophila.
Figure 4. Volcano plots between controls and HM (A) and PM (B) pacu (P. mesopotamicus) challenged with A. hydrophila. Each dot represents the expression (x-axis, Log2 fold change), and significance (y-axis, −Log10 Q-value) of each transcript. Differentially expressed transcripts (q-value < 0.05) are red.
A significant fraction of the differentially expressed transcripts showed log2FC > | 5|, specifically 400 genes (9%) differentially expressed for HM group, and 5 genes (1.2%) differentially expressed in PM group. Further, a total of 4,271 genes (96.8%) were differentially expressed only in HM group, whereas 291 genes (67.2%) were differentially expressed exclusively in PM group (Figures 4A,B and Table 4).
In this study, only genes specifically related to the immune responses against A. hydrophila and differentially expressed in HM (Table 5) and PM (Table 6) groups were considered. Therefore, 57 transcripts identified as immune transcripts were considered differentially expressed in the HM group (31 up-regulated and 26 down-regulated transcripts) (q-value < 0.09) (Table 5), whereas 23 immune transcripts were considered differentially expressed in the PM group (15 up-regulated and 8 down-regulated transcripts) (q-value < 0.07) (Table 6). For didactic purposes, these transcripts were grouped and classified according to the important immune mechanisms that may be related to fighting infections such as those caused by A. hydrophila, such as: cytokines and signaling, epithelial protection, antigen processing and presentation, apoptosis, phagocytosis, complement system cascades and pattern recognition receptors (Tables 5, 6).
Table 5. List of the most significant differentially expressed immune system transcripts in HM group of P. mesopotamicus exposed to A. hydrophila infection vs. the controls.
Table 6. List of most significant differentially expressed immune system transcripts in PM group of P. mesopotamicus exposed to A. hydrophila infection vs. the controls.
De novo assembly from pacu liver RNA sequencing generated a high-quality transcriptome containing 211,259 transcripts, presenting an N50 value of 1,296 bp and average length of 713 bp. Moreover, default settings (kmer-length fixed to 25 and minimum contig length of 200 bp) used by the Trinity software provided an optimal assembly, as they subsequently allowed 92% of the reads to be mapped back to the transcriptome, having few transcriptomic information missing. The evident quality of the transcriptome was similar to other RNA-seq studies involving infections in fish (Wang et al., 2013; Kumar et al., 2017; Song et al., 2017). Therefore, this immune-enriched de novo pacu liver transcriptome can be used as framework to investigate transcriptomic changes upon infection by A. hydrophila and other pathogens.
This is the first study on gene expression of pacu infected by A. hydrophila. The choice to focus only on genes related to the immune system was due to the generation of subsidies for future studies on application of immunostimulating treatments that could increase the pacu resistance when exposed to infectious outbreaks of A. hydrophila. Additionally, this results will contribute to develop effective strategies to protect pacu production from this threatening bacterial disease, and could be exploited through future development of breeding programs to increase the genetic resistance of pacu to A. hydrophila infection. Our results were similar to previous transcriptomic analyses in different species of fish infected by A. hydrophila (Tran et al., 2015; Dang et al., 2016; Qin et al., 2017; Jiang et al., 2018; Maekawa et al., 2019) showing transcripts related to innate immune system, which consist in a complex set of pathways that mediate the activation of phagocytosis, inflammation and destruction of the pathogen by antimicrobial substances, stimulating adaptive immune responses (Tort et al., 2003). In order to facilitate the discussion we have grouped these differentially expressed genes in categories based on their putative molecular mechanisms.
Cytokines and Signaling
Cytokines are small proteins released by immune cells acting as signaling molecules that regulate immunity, inflammation and hematopoiesis. We found the Interleukin-1 receptor type 1 (IL1R1) and Interleukin-1 receptor type 2 (IL1R2) transcripts up-regulated in HM and PM groups. These transcripts are part of the Interleukin-1 receptor complex, involved in many cytokine-induced inflammatory responses.
We also found four chemokines upregulated and four downregulated in HM group. Chemokines are cytokines that regulate the leukocyte migration and the differentiation of recruited cells under inflammatory conditions, in order to orchestrate the first steps of both innate and adaptive immune responses (Bird and Tafalla, 2015). Further, we found downregulated genes involved in chemokine activation in HM group, and one transcript related to the negative regulation of cytokines (Suppressor of cytokine signaling, SOCS1), downregulated in PM group. Excess production of anti-inflammatory cytokines can compromise the ability to clear microorganisms through suppression of immune cell function. If a balance is not maintained, the result is either an excessive proinflammatory response or an immunosuppression (Chaudhry et al., 2015).
Phagocytosis is an essential cellular process that captures and destroys foreign particles, such as small unicellular organisms. One of the major effects of A. hydrophila infection in fish is characterized by the production of toxins that can cause severe tissue destruction, affecting the epithelial cells due to the immune phagocytic activity in response to the bacteria (Cipriano et al., 2001; Abdelhamed et al., 2017). According to the results, Natural resistance-associated macrophage protein 2 (NRAMP2) transcript was upregulated only in PM group. This protein has been previously considered an important regulator of macrophages in the response to A. hydrophila (Jiang et al., 2018). Further, a previous study on cytochrome B-245 beta chain (CYBB) transcription, upregulated in HM group, revealed that deficiency in the production of the protein impairs phagocytosis-mediated defense against A. hydrophila. This means that even if the phagocytes are able to phagocytize the bacteria, A. hydrophila can survive into the phagocytic vacuoles (Jiang et al., 2018).
Pattern Recognition Receptors
Interleukin-1 receptor-associated kinase 1 (IRAK1) transcript is directly involved in the recognition of pathogen-derived products and the consequent initiation of immune and inflammatory responses. This protein is considered part of Toll-like receptor complex, and it is usually expressed on sentinel cells, such as macrophages and dendritic cells, which recognize structurally conserved molecules derived from pathogens. Interleukin 1 receptor-associated kinase 3 (IRAK3), which inhibits the dissociation of IRAK1 from the Toll-like receptor signaling complex (Shan et al., 2015), is upregulated in HM group. Therefore, these results can suggest that these genes are involved in the induction of the inflammatory response in response to A. hydrophila infection in pacu.
The complement system is a key component of the innate immune system of fish. The activation of the complement pathway leads to opsonization of pathogens, recruitment of inflammatory and immune competent cells, and direct killing of pathogens through lysis (Merle et al., 2015). In the present study, practically all complement transcripts identified were downregulated in HM group (especially C1, C3, and C4), suggesting that complement pathways is not activated following Aeromonas hydrophila infection, which is different than reported in previous studies considering other fish species (Tran et al., 2015; Dang et al., 2016; Jiang et al., 2018). Complement system inactivation in HM group could be an important factor leading to the failure of the defense system, and consequent host death during the PM time-series that we found in the present study. However, further studies need to be performed in order to validate these genes and corroborate this hypothesis. In the same way, several coagulation factors were also downregulated in HM group, such as Coagulation factor VII (F7), Coagulation Factor IX (F9), and Coagulation Factor XI (F11), indicating the absence of blood coagulation, and favoring the occurrence of hemorrhages, as observed in HM group with individuals presenting clinical signs of A. hydrophila infection.
Transcriptome analysis by RNA-seq allowed the identification and comparative quantification of gene expression in pacu infected by A. hydrophila. De novo liver transcriptome assembled in this study enables capturing expression levels of specific transcripts associated with important immune pathways, such as the coagulation factor and complement system cascades. Using this transcriptome, we have described for the first time the response of pacu infected by A. hydrophila in two critical periods of infection (PM and HM) comparing individuals belonging two distinct groups based on different phases of bacterial infection (acute mortality and plateau of mortality). To our knowledge this is the first descriptive study on transcriptome analysis of infected pacu that may serve as a reference to understand the genetic mechanisms underlying the host-pathogen interaction. The differences in the transcriptomic responses to a standardized A. hydrophila challenge was able to provide a initial data set to begin to understand the impact of selective breeding on genetic basis of disease resistance. Therefore, further studies trying to connect the genomic and genetic basis of resistance to A. hydrophila are fundamental to enable the implementation of A. hydrophila resistance into pacu breeding programs. Despite these limitations, these results may suggest possible suppression of genes belonging to the immune system that may contribute to susceptibility to infection, as for example genes related with complement system and coagulation factor cascades. Such supression of important immune system genes will serve as subsidies for future studies aiming at the use of prophylactic and immunostimulating treatments that could increase the pacu resistance when exposed to infectious outbreaks of A. hydrophila.
Data Availability Statement
All RNA-Seq data generated for this study were deposited in the National Center for Biotechnology Information Sequence Read Archive under BioProject ID PRJNA632934. Bioinformatics pipelines are available in Supplementary Table S1.
This study was conducted in strict accordance with the recommendations of the National Council for Control of Animal Experimentation (CONCEA) (Brazilian Ministry for Science, Technology and Innovation) and was approved by the Ethics Committee on Animal Use (CEUA number 19.005/17) of Faculdade de Ciências Agrárias e Veterinárias, UNESP, Campus Jaboticabal, SP, Brazil.
VM-F and DH: conceptualization, funding acquisition, and roles/writing – original draft. VM-F, DR, RH, and DH: formal analysis. VM-F, MH, RK, MF, RA, and DP: investigation and methodology. DH: project administration. DH, DR, and RH: supervision. All authors: writing – review and editing.
This work was supported by the São Paulo Research Foundation (FAPESP grants 2016/21011-9, 2017/26900-9, and 2016/18294-9), the Comisión Nacional de Investigación Científica y Tecnológica (International Call FAPESP-CONICYT 2018/08416-5), the National Council for Scientific and Technological Development (CNPq grant 311559/2018-2), the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior – Brazil (CAPES – Finance Code 001), and the BBSRC Institute Strategic Program Grants (BB/P013759/1 and BB/P013740/1).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2020.00604/full#supplementary-material
Abdelhamed, H., Ibrahim, I., Baumgartner, W., Lawrence, M. L., and Karsi, A. (2017). Characterization of histopathological and ultrastructural changes in channel catfish experimentally infected with virulent Aeromonas hydrophila. Front. Microbiol. 8:1519. doi: 10.3389/fmicb.2017.01519
Andrews, S. (2010). FastQC: A Quality Control Tool for High Throughput Sequence Data. Available online at: http://www.bioinformatics.babraham.ac.uk/projects/fastqc (accessed April 15, 2018).
Belém-Costa, A., and Cyrino, J. E. P. (2006). Antibiotic resistence of Aeromonas hydrophila isolated from Piaractus mesopotamicus (Holmberg, 1887) and Oreochromis niloticus (Linnaeus, 1758). Sci. Agric. 63, 281–284. doi: 10.1590/s0103-90162006000300011
Bryant, D. M., Johnson, K., DiTommaso, T., Tickle, T., Couger, M. B., Payzin-Dogru, D., et al. (2017). A tissue-mapped axolotl de novo transcriptome enables identification of limb regeneration factors. Cell Rep. 18, 762–776. doi: 10.1016/j.celrep.2016.12.063
Carraschi, S. P., Cruz, C., Neto, J. G. M., Moraes, F., Rossi, O. D., Neto, A. N., et al. (2012). Evaluation of experimental infection with Aeromonas hydrophila in pacu (Piaractus mesopotamicus) (Holmberg, 1887). Int. J. Fish. Aquac. 4, 81–84. doi: 10.5897/IJFA11.008
da Pimentel, J. S. M., do Carmo, A. O., de Maciel, D. C. L., de Faria Siqueira, F., and Kalapothakis, E. (2014). Complete mitochondrial genome sequence of Piaractus mesopotamicus (Holmberg, 1887). Mitochondrial DNA 1736, 1–2. doi: 10.3109/19401736.2014.971297
Dang, Y., Xu, X., Shen, Y., Hu, M., Zhang, M., Li, L., et al. (2016). Transcriptome analysis of the innate immunity-related complement system in spleen tissue of Ctenopharyngodon idella infected with Aeromonas hydrophila. PLoS One 11:e0157413. doi: 10.1371/journal.pone.0157413
Farias, T. H. V., Levy-Pereira, N., de Oliveira Alves, L., de Carla Dias, D., Tachibana, L., Pilarski, F., et al. (2016). Probiotic feeding improves the immunity of pacus, Piaractus mesopotamicus, during Aeromonas hydrophila infection. Anim. Feed Sci. Tech. 211, 137–144. doi: 10.1016/j.anifeedsci.2015.11.004
Grabherr, M. G., Haas, B. J., Yassour, M., Levin, J. Z., Thompson, D. A., Amit, I., et al. (2011). Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat. Biotechnol. 29, 644–652. doi: 10.1038/nbt.1883
Haas, B. J., Papanicolaou, A., Yassour, M., Grabherr, M., Blood, P. D., Bowden, J., et al. (2013). De novo transcript sequence reconstruction from RNA-seq using the Trinity platform for reference generation and analysis. Nat. Protoc. 8, 1494–1512. doi: 10.1038/nprot.2013.084
Jiang, H., Wang, M., Fu, L., Zhong, L., Liu, G., Zheng, Y., et al. (2020). Liver transcriptome analysis and cortisol immune-response modulation in lipopolysaccharide-stimulated in channel catfish (Ictalurus punctatus). Fish Shellfish Immun. 101, 19–50. doi: 10.1016/j.fsi.2020.03.024
Jiang, N., Fan, Y., Zhou, Y., Wang, W., Ma, J., and Zeng, L. (2018). Transcriptome analysis of Aeromonas hydrophila infected hybrid sturgeon (Huso dauricus × Acipenser schrenckii). Sci. Rep. 8, 1–18. doi: 10.1038/s41598-018-36376-2
Kanehisa, M., Goto, S., Sato, Y., Furumichi, M., and Tanabe, M. (2012). KEGG for integration and interpretation of large-scale molecular data sets. Nucleic Acids Res. 40, 109–114. doi: 10.1093/nar/gkr988
Kumar, R., Sahoo, P. K., and Barat, A. (2017). Transcriptome profiling and expression analysis of immune responsive genes in the liver of Golden mahseer (Tor putitora) challenged with Aeromonas hydrophila. Fish Shellfish Immun. 67, 655–666. doi: 10.1016/j.fsi.2017.06.053
Lafferty, K. D., Harvell, C. D., Conrad, J. M., Friedman, C. S., Kent, M. L., Kuris, A. M., et al. (2015). Infectious diseases affect marine fisheries and aquaculture economics. Ann. Rev. Mar. Sci. 7, 471–496. doi: 10.1146/annurev-marine-010814-015646
Mahapatra, K., Das Gjerde, B., Sahoo, P. K., Saha, J. N., Barat, A., Sahoo, M., et al. (2008). Genetic variations in survival of rohu carp (Labeo rohita, Hamilton) after Aeromonas hydrophila infection in challenge tests. Aquaculture 279, 29–34. doi: 10.1016/j.aquaculture.2008.03.054
Marancik, D., Gao, G., Paneru, B., Ma, H., Hernandez, A. G., Salem, M., et al. (2015). Whole-body transcriptome of selectively bred, resistant-, control-, and susceptible-line rainbow trout following experimental challenge with Flavobacterium psychrophilum. Front. Genet. 5:453. doi: 10.3389/fgene.2014.00453
Mastrochirico-Filho, V. A., Ariede, R. B., Freitas, M. V., Lira, L. V. G., Agudelo, J. F. G., Pilarski, F., et al. (2019). Genetic parameters for resistance to Aeromonas hydrophila in the Neotropical fish pacu (Piaractus mesopotamicus). Aquaculture 513:734442. doi: 10.1016/j.aquaculture.2019.734442
Ødegård, J., Olesen, I., Dixon, P., Jeney, Z., Nielsen, H. M., Way, K., et al. (2010). Genetic analysis of common carp (Cyprinus carpio) strains. II: resistance to koi herpesvirus and Aeromonas hydrophila and their relationship with pond survival. Aquaculture 304, 7–13. doi: 10.1016/j.aquaculture.2010.03.017
Polinski, M. P., Bradshaw, J. C., Inkpen, S. M., Richard, J., Fritsvold, C., Poppe, T. T., et al. (2016). De novo assembly of Sockeye salmon kidney transcriptomes reveal a limited early response to piscine reovirus with or without infectious hematopoietic necrosis virus superinfection. BMC Genomics 17:848. doi: 10.1186/s12864-016-3196-y
Powell, S., Szklarczyk, D., Trachana, K., Roth, A., Kuhn, M., Muller, J., et al. (2012). eggNOG v3.0: orthologous groups covering 1133 organisms at 41 different taxonomic ranges. Nucleic Acids Res. 40, 284–289. doi: 10.1093/nar/gkr1060
Qin, C., Gong, Q., Wen, Z., Yuan, D., Shao, T., Wang, J., et al. (2017). Transcriptome analysis of the spleen of the darkbarbel catfish Pelteobagrus vachellii in response to Aeromonas hydrophila infection. Fish Shellfish Immun. 70, 498–506. doi: 10.1016/j.fsi.2017.09.042
Robledo, D., Gutiérrez, A. P., Barría, A., Yáñez, J. M., and Houston, R. D. (2018). Gene expression response to sea lice in Atlantic salmon skin: RNA sequencing comparison between resistant and susceptible animals. Front. Genet. 9:287. doi: 10.3389/fgene.2018.00287
Robledo, D., Ronza, P., Harrison, P. W., Losada, P. P., Bermúdez, R., Pardo, B. G., et al. (2014). RNA-seq analysis reveals significant transcriptome changes in turbot (Scophthalmus maximus) suffering severe enteromyxosis. BMC Genomics 15:1149. doi: 10.1186/1471-2164-15-1149
Shan, S. J., Liu, D. Z., Wang, L., Zhu, Y. Y., Zhang, F. M., Li, T., et al. (2015). Identification and expression analysis of irak1 gene in common carp Cyprinus carpio L.: indications for a role of antibacterial and antiviral immunity. J. Fish Biol. 87, 241–255. doi: 10.1111/jfb.12714
Silva, B. C., Mouriño, J. L. P., Vieira, F. N., Jatobá, A., Seiffert, W. Q., and Martins, M. L. (2012). Haemorrhagic septicaemia in the hybrid surubim (Pseudoplatystoma corruscans×Pseudoplatystoma fasciatum) caused by Aeromonas hydrophila. Aquac. Res. 43, 908–916. doi: 10.1111/j.1365-2109.2011.02905.x
Smith-Unna, R., Boursnell, C., Patro, R., Hibberd, J. M., and Kelly, S. (2016). TransRate: reference free quality assessment of de novo transcriptome assemblies. Genome Res. 26:021626. doi: 10.1101/gr.1964 69.115
Song, X., Hu, X., Sun, B., Bo, Y., Wu, K., Xiao, L., et al. (2017). A transcriptome analysis focusing on inflammation-related genes of grass carp intestines following infection with Aeromonas hydrophila. Sci. Rep. 7, 1–12. doi: 10.1038/srep40777
Stentiford, G. D., Sritunyalucksana, K., Flegel, T. W., Williams, B. A. P., Withyachumnarnkul, B., Itsathitphaisarn, O., et al. (2017). New paradigms to help solve the global aquaculture disease crisis. PLoS Pathog. 13:e1006160. doi: 10.1371/journal.ppat.1006160
Sudhagar, A., Kumar, G., and El-Matbouli, M. (2018). Transcriptome analysis based on RNA-Seq in understanding pathogenic mechanisms of diseases and the immune system of fish: a comprehensive review. Int. J. Mol. Sci. 19, 245. doi: 10.3390/ijms19010245
Tran, N. T., Gao, Z., Zhao, H., Yi, S., Chen, B., Zhao, Y., et al. (2015). Transcriptome analysis and microsatellite discovery in the blunt snout bream (Megalobrama amblycephala) after challenge with Aeromonas hydrophila. Fish Shellfish Immun. 45, 72–82. doi: 10.1016/j.fsi.2015.01.034
Wang, R., Sun, L., Bao, L., Zhang, J., Jiang, Y., Yao, J., et al. (2013). Bulk segregant RNA-seq reveals expression and positional candidate genes and allele-specific expression for disease resistance against enteric septicemia of catfish. BMC Genomics 14:929. doi: 10.1186/1471-2164-14-929
Xiong, X., Chen, Y., Liu, L., Wang, W., Robinson, N. A., and Gao, Z. (2017). Estimation of genetic parameters for resistance to Aeromonas hydrophila in blunt snout bream (Megalobrama amblycephala). Aquaculture 479, 768–773. doi: 10.1016/j.aquaculture.2017.07.011
Zanuzzo, F. S., Sabioni, R. E., Montoya, L. N. F., Favero, G., and Urbinati, E. C. (2017). Aloe vera enhances the innate immune response of pacu (Piaractus mesopotamicus) after transport stress and combined heat killed Aeromonas hydrophila infection. Fish Shellfish Immun. 65, 198–205. doi: 10.1016/j.fsi.2017.04.013
Keywords: disease resistance, complement system, de novo assembly, RNA-Seq, aquaculture
Citation: Mastrochirico-Filho VA, Hata ME, Kuradomi RY, de Freitas MV, Ariede RB, Pinheiro DG, Robledo D, Houston R and Hashimoto DT (2020) Transcriptome Profiling of Pacu (Piaractus mesopotamicus) Challenged With Pathogenic Aeromonas hydrophila: Inference on Immune Gene Response. Front. Genet. 11:604. doi: 10.3389/fgene.2020.00604
Received: 28 February 2020; Accepted: 18 May 2020;
Published: 09 June 2020.
Edited by:Argyris Papantonis, University Medical Center Göttingen, Germany
Reviewed by:Bhanwar Lal Puniya, University of Nebraska-Lincoln, United States
Leda Torres, National Institute of Pediatrics of Mexico, Mexico
Copyright © 2020 Mastrochirico-Filho, Hata, Kuradomi, de Freitas, Ariede, Pinheiro, Robledo, Houston and Hashimoto. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Diogo Teruo Hashimoto, email@example.com