RNA-Seq based transcriptome of whole blood from immunocompetent pigs (Sus scrofa) experimentally infected with Mycoplasma suis strain Illinois
Veterinary Research volume 49, Article number: 49 (2018)
Pigs are popular animal models in biomedical research. RNA-Seq is becoming the predominant tool to investigate transcriptional changes of the pig’s response to infection. The high sensitivity of this tool requires a strict control of the study design beginning with the selection of healthy animals to provide accurate interpretation of research data. Pigs chronically infected with Mycoplasma suis often show no obvious clinical signs, however the infection may affect the validity of animal research. The goal of this study was to investigate whether or not this silent infection is also silent at the host transcriptional level. Therefore, immunocompetent pigs were experimentally infected with M. suis and transcriptional profiles of whole blood, generated by RNA-Seq, were analyzed and compared to non-infected animals. RNA-Seq showed 55 differentially expressed (DE) genes in the M. suis infected pigs. Down-regulation of genes related to innate immunity (tlr8, chemokines, chemokines receptors) and genes containing IFN gamma-activated sequence (gbp1, gbp2, il15, cxcl10, casp1, cd274) suggests a general suppression of the immune response in the infected animals. Sixteen (29.09%) of the DE genes were involved in two protein interaction networks: one involving chemokines, chemokine receptors and interleukin-15 and another involving the complement cascade. Genes related to vascular permeability, blood coagulation, and endothelium integrity were also DE in infected pigs. These findings suggest that M. suis subclinical infection causes significant alterations in blood mRNA levels, which could impact data interpretation of research using pigs. Screening of pigs for M. suis infection before initiating animal studies is strongly recommended.
The pig (Sus scrofa) is an important meat-production animal and the most commonly used large animal in research, particularly as a mammalian model for the study of human diseases . Their size, anatomy, and physiology as well as the progression of many diseases, such as diabetes, atherosclerosis, and cardiovascular disease, are similar to humans . Further, the improvement of techniques for gene modification has resulted in an increase in the number of genetically engineered swine models [2, 3]. Because of these factors, the National Institutes of Health (NIH) has established the National Swine Research and Resource Center to help meet the needs of the biomedical community .
The sequencing of the complete pig genome  confirmed that its gene content, sequence and chromosomal structure are highly conserved relative to the human genome [5, 6], making the pig a particularly useful model for genome-based studies. Accordingly, RNA sequencing (RNA-Seq) using next-generation sequencing (NGS) technology is becoming the predominant tool for large-scale gene expression analyses. Numerous studies using RNA-Seq methodology have been performed to investigate global gene expression in various tissues or cells from pigs of different breeds or phenotypes [7,8,9,10,11,12,13,14]. Since transcriptome analyses are extremely sensitive and specific, they may result in misleading and inaccurate functional genomic data if environmental and host factors are not tightly controlled in the study design. The selection of healthy pigs is the first step toward the goal of generating valid research data that will provide results that are more predictive of human diseases.
One of the most important host factors that may affect the validity of animal research is the presence of pre-existing subclinical infectious diseases in test and/or control animals. Mycoplasma suis is notorious for causing a chronic infection without obvious clinical signs. This organism typically attaches to porcine erythrocytes, and is the only hemotrophic mycoplasma shown to invade its host cell . When chronically infected animals become stressed or undergo splenectomy, overt infection develops and they may succumb to an acute, life-threatening disease characterized by hemolytic anemia, hypoglycemia, and bleeding accompanied by intravascular coagulation and endothelial damage [16, 17]. Some reports also suggest that pigs chronically infected with M. suis may have impaired immune defenses, being more susceptible to respiratory and enteric diseases [18, 19]. Research pigs, however, are not routinely tested for the presence of M. suis in their bloodstream. We have shown that levels of bacteria may fluctuate from 103 to 109/mL of blood for months in naturally infected pigs . Considering the high within-herd prevalence (up to 86%) and wide geographic distribution [21,22,23], M. suis constitutes a major threat to the use of pigs in research.
Hemotrophic mycoplasmas, such as M. suis, have been described in several mammalian species and the recrudescence of acute disease in chronically infected animals, especially when stressed, has been reported to have a deleterious effect on the course of experimental research in splenectomized pigs and dogs, non-human primates (Saimiri sciureus), rodents, and sheep [24,25,26,27,28]. These findings strengthen our hypothesis that chronic infections, while clinically silent, are unlikely to be silent at the host transcriptional level and should be further explored. Therefore, experimental M. suis infection of immunocompetent (non-splenectomized) pigs, a model for mimicking natural subclinical infection, was performed to study its effects on the blood transcriptional profile of the porcine host.
Materials and methods
The Purdue Animal Care and Use Committee, Protocol Number 1111000223, approved this study. Six domestic piglets (Sus scrofa domesticus), all females from the same litter, were purchased from the Animal Sciences Research and Education Center-ASREC, Purdue University, West Lafayette, IN, USA. The animals were received when they were 20 days old and kept until turning 99 days old at Purdue animal housing facilities. Piglets were fed antibiotic-free feed, as well as water ad libitum for the entire study.
The M. suis strain Illinois, obtained in 1985 from a naturally infected pig having acute infection and high loads of bacteremia, was used in an experimental infection study . After 14 years of storage (−80 °C), M. suis strain Illinois, as demonstrated by experimental infection of pigs, retained its pathogenicity . Aliquots collected from a specific pathogen free (SPF) pig inoculated with M. suis and at the peak of bacteremia , representing the second passage of this strain, were frozen at −80 °C and used in the experimental infection described herein.
Study design and experimental infection
Piglets were considered M. suis—free by testing blood samples using qPCR  13 days after birth at ASREC and twice (at 33 and 53 days old) prior to the infection day while at Purdue animal housing facilities. The mother also tested negative for M. suis by qPCR while pregnant, 50 days prior to farrowing. They were randomly assigned to either a control group (n = 3; pigs #1, #2, and #3) or an experimental M. suis infection group (n = 3; pigs #5, #6, and #7). The control group was housed in a separate room from the infection group to avoid M. suis transmission. Thirty-three days after arrival (day 0), experimental infection was performed. Briefly, animals were sedated with TX (telazol reconstituted with 500 mg of xylazine, 1.0 mL/100 kg of body weight); three piglets in the infection group were intravenously inoculated through the ear vein with 1.0 mL of cryopreserved blood containing M. suis strain Illinois at a concentration of 108 to 109 organisms/mL of blood. After inoculation, 1.0 mL of heparin sodium was flushed through the butterfly catheter. Pigs in the control group (sham-infected) were also intravenously inoculated, using 1.0 mL of blood collected from another control pig, collected at the same time. The animals were monitored daily (minimum twice a day) for direct observation of behavior (BAR status) and body temperature. Blood was collected through the anterior venous cava vein and added to EDTA tubes every 3–4 days for monitoring infection by qPCR , hematocrit, and detection of specific antibodies for M. suis recombinant GrpE (rGrpE) antigen using microbead immunoassay (MIA) . Pigs were considered seropositive when their median fluorescence intensity (MFI) crossed the cut-off defined as the mean plus 3 standard deviation of the pre-inoculation MFIs (day 0) for all animals.
RNA preparation and RNA-sequencing (RNA-Seq) procedure
Six aliquots of 400 L EDTA-whole blood from each animal were immediately transferred after collection to 2 mL tubes containing TRIzol® reagent (Life Technologies Corp., Carlsbad, CA, USA) and stored at −80 °C until RNA extraction. Total RNA was extracted using TRIzol® (Life Technologies Corp.) followed by RNeasy Mini Kit, including on-column DNase treatment step (QIAGEN® Inc., Valencia, CA, USA) . Concentration and quality of the RNA samples were assessed using Agilent 2100 Bioanalyzer (Agilent Technologies, Inc. Headquarters, CA, USA), and samples with highest RIN (RNA integrity number) and concentration were chosen for mRNA purification. Enriched mRNA (purified with Globin-Zero™ Gold Kit, Epicentre®, by Illumina® company, Madison, WI, USA) was sent to the Genomics Core Facility at Purdue University for library construction and RNA-sequencing. One RNA sample from each pig was selected for transcriptome analysis based on the first, highest load of M. suis detected (peak of bacteremia; shown in results) in the peripheral blood of the infected pigs according to qPCR results. The days selected for the infected group were: pig #5, 30 days post-inoculation (dpi); pig #6, 23 dpi, and pig #7, 23 dpi; the highest peaks of bacteremia for each animal within 7–8 days from each other. For control pigs (#1, #2, and #3), 27 dpi was chosen for transcriptome analysis, as an intermediate day among days 23 and 30. Samples were sequenced from paired-end libraries (TruSeq DNA sample preparation kit, Illumina Inc., San Diego, CA, USA) barcoded with different tags using one lane of Illumina HiSeq 2000.
Differential gene expression analysis of pigs infected with M. suis and control pigs
Illumina reads were pre-processed for mapping by assessing quality using FastQC (v 0.10.1) , and trimming was done using FASTX toolkit (v 0.0.13). The quality trimmed reads were mapped against the bowtie2-indexed Sus scrofa genome (version 75 from ENSEMBL) using Tophat (v 2.0.9)  with default parameters. HTSeq (v 0.5.3p7) was used to generate raw read counts from each sample for each gene feature using Tophat output and the known Sus scrofa annotations. Genes with zero counts were removed from the count matrix and genes where the count was zero in some samples and not in others were changed to one in order to prevent infinite values during the statistical analysis.
Three different methods were used to perform pairwise analyses of differential gene expression between control and M. suis-infected groups. In the first two methods, edgeR (v 3.0.3) and DESeq 2 (v 1.0.19), “R” package (Version 2.15.2)  was used to statistically analyze the results based on the count matrix mentioned above. And finally, in the third method (Cufflinks), bam files from mapping to the genome were analyzed using the Cufflinks (v 2.0.2) suite of programs (Cufflinks, Cuffmerge and Cuffdiff) . The lists of significant differentially expressed (DE) genes were generated by using the adjusted p-value or FDR values. DE genes identified by at least two methods (edgeR, DESeq 2 and Cufflinks) (cut-off of p-value ≤ 0.1) were defined as truly DE genes.
Quantitative RT-PCR validation of selected genes
Selected up and down DE genes were validated using two step quantitative RT-PCR (qRT-PCR). Briefly, independent cDNA synthesis was performed for all samples (three infected and three control pigs) starting from 100 ng of total RNA, same total RNA used for RNA-Seq procedure, using SuperScript® VILO™ cDNA Synthesis Kit (Invitrogen, Carlsbad, CA, USA). qPCR was carried out in a 7300 Real-Time PCR System (Applied Biosystems®, Foster City, CA, USA) using Power SYBR® Green PCR Master Mix (Applied Biosystems®) as chemistry according to manufacturer’s instructions. Primers for each target gene selected for validation by qRT-PCR, in addition to GAPDH, used as reference gene, were designed using Primer3web version 4.0.0  (Table 1). In addition to the biological replicates (cDNA from three infected and three control pigs) three technical replicates for each target gene (DE gene selected from RNA-Seq results) were included in the qPCR assay. GAPDH replicates were added and analyzed in every assay for correct comparison of target genes between infected and control pigs. Relative quantification was evaluated using the Comparative Ct method with GAPDH gene as reference gene . Fold change was calculated as the difference between infected and control group. Graphics of relative expression display the mean ± standard deviation of biological replicates.
The relative expression of the target genes tested by qRT-PCR between infected and control groups was compared by t-test (p-value ≤ 0.1) using JMP® Software (JMP®, SAS Institute Inc., Cary, NC, 1989–2007).
Network associations and gene ontology annotation of DE genes
DE genes identified in M. suis infected pigs were analyzed using STRING database with default settings for prediction of network associations between proteins . In addition, Blast2Go software was used to assign functional categories to each protein encoded by DE genes based on gene ontology (GO) terms .
Hematocrit, rectal temperature, M. suis loads, and humoral response to rGrpE
Results of hematocrit and rectal temperature of M. suis-infected and sham-inoculated groups are shown in Figure 1. The hematocrit of M. suis infected pigs ranged from 23 to 35%, while in sham-inoculated animals this number varied from 26 to 29% (Figure 1A). Interestingly, the lowest hematocrit values (pig #5: 23%; pig #6: 24%; pig #7: 23%) were observed 10 dpi in M. suis infected-pigs. The rectal temperature of M. suis-infected and sham-infected pigs ranged from 36.8 to 40.4 °C and from 37.7 to 39.8 °F, respectively (pig reference interval 38.7–39.8 °C) (Figure 1B). High rectal temperatures (> 39.8 °C, maximum of 40.4 °C) were observed intermittently during the first 7 dpi only in pigs #6 and #7 from the M. suis-infected group (Figure 1B). Despite these minor increases in rectal temperatures of M. suis-infected pigs, none of the six animals displayed obvious clinical signs of infection.
Results of M. suis copy number/mL of blood and humoral response (MIA) of M. suis-infected and sham-inoculated pigs are shown in Figure 2. All three pigs inoculated with M. suis strain Illinois became infected (as shown by qPCR results) and seroconverted (as shown by MIA results). Infected pigs showed the first qPCR positive result for M. suis on different days following experimental infection, varying from day 10 (pig #6) to day 27 (pig #5). Bacterial loads varied from approximately 103 to 109 organisms/mL of blood. Interestingly, pig #7 was qPCR positive for M. suis on only two separate occasions, having a relatively low bacterial load (103 organisms/mL of blood). According to the established MFI cut-off (74.07), all three pigs (infected group) seroconverted by days 14 (pig #6), 34 (pig #5) or 37 (pig #7) post inoculation and remained seropositive (pigs #5 and #6) until the last day evaluated (37 dpi) (Figure 2; MFI values varied from 20.5 to 939; individual MFI values are given in Additional file 1). It is important to note that, although pigs #5 and #7 seroconverted at a later day post-inoculation, both animals showed a constant increase on their individual MFI values at earlier days post-inoculation (Additional file 1), which likely represents earlier seroconversion. Sham-infected animals remained M. suis qPCR and MIA negative throughout the experiment (MFI values varied from 25.5 to 58.2; Additional file 1).
RNA-Seq results and validation by qRT-PCR
Total RNA concentrations and quality of the samples were as follows: pig #1, 156.5 ng/µL, RIN 8.8; pig #2, 77.7 ng/µL, RIN 8.5; pig #3, 107.9 ng/µL, RIN 8.9; pig #5, 113.1 ng/µL, RIN 8; pig #6 30.9 ng/µL, RIN 8; pig #7 115.7 ng/µL, RIN 8.5. A Venn diagram showing the number of DE genes identified (p-value ≤ 0.1) in the M. suis-infected group when compared to the sham-inoculated group by each of the three analysis methods (DESeq 2, edgeR and Cufflinks) is shown in Figure 3. A list of all DE genes and their respective fold change and p-values are available in Additional file 2. Twenty DE genes were detected by all three methods, while 55 DE genes, 3 up and 52 down-regulated, were identified by at least two of the three methods (p-value ≤ 0.1) (Figure 3). Forty-seven out of the 55 DE genes (85.45%) were predicted to encode proteins with assigned functions, while 8 (14.55%) encode proteins with uncharacterized functions (Additional file 2). The raw RNA-Seq data was submitted to the Sequence Read Archive (SRA) database at the National Center for Biotechnology Information (NCBI) under Bioproject Accession Number PRJNA321932.
Approximately 20% of the 55 DE genes (n = 10) were chosen for qRT-PCR validation based on their role in immune response and/or fold change of 2× or more (Table 2). The relative gene expression results (mean and standard deviation) of DE genes validated by qRT-PCR are shown on Additional file 3. All genes were confirmed to be up or down regulated based on qRT-PCR results (p-value ≤ 0.1) (Table 2).
Network analysis of the relationships between proteins encoded by DE genes
STRING database predictions of functional association networks for all 55 proteins encoded by DE genes identified by 2 or more methods are shown in Figure 4. A main network involving chemokines, chemokine receptors and interleukin-15 was observed (encoded by genes CXCL10, CCL8, PPBP, CCR5 and IL15). This network is part of pathways ID GO:0006954 and GO:0006955 (biological process, gene ontology), and 04060 and 04062 (KEGG pathways), corresponding to inflammatory response, immune response, cytokine–cytokine receptor interaction, and chemokine signaling pathways, respectively. Interactions among genes involved in programmed cell death (CASP4, BCL2L14 and CASP1) were also noted. Interestingly, genes containing IFN gamma-activated sequence (GAS) (e.g. GBP1, GBP2, IL15, CXCL10, CASP1, and CD274) are scattered throughout different networks, and a strong interaction between C2, the gene encoding for Complement 2 protein, and the C1 inhibitor (C1-INH, encoded by SERPING1 gene), both involved in the complement cascade, was detected. Further, many proteins (28/55, 50.9%) that showed no direct relationships to other proteins were also identified in the analysis.
Gene ontology analysis of differentially expressed (DE) genes of M. suis infected pigs
DE genes of M. suis infected pigs were analyzed using Blast2GO bioinformatics tool, which provides functional annotation of sequences by mapping the Blast hits of each sequence with gene ontology terms assigning molecular functions, biological processes and/or cellular components to each sequence. A total of 314 different functional categories, identified as biological processes, were assigned to the proteins encoded by DE genes (Additional file 4). Many proteins were assigned more than one “biological processes” roles. The top five identified functions under the “biological processes” category are listed on Table 3, which includes 18 of the 55 (32.7%) identified genes that were DE by 2 methods or more. These 5 functions of the “biological processes” category had at least 4 genes each. The following GO categories were also important (representing 3 or 4 genes each): chemokine-mediated signaling pathway; chemotaxis defense response to virus; G-protein coupled receptor signaling pathway; ion transport; negative regulation of T cell proliferation; oxidation–reduction process; positive regulation of T cell proliferation; protein phosphorylation; regulation of apoptotic process; response to lipopolysaccharide; transport; signal transduction (n = 3 each); ion transport (n = 4). These functional roles indicate that DE genes were mostly associated with immune and inflammatory responses.
In this study, immunocompetent pigs experimentally infected with M. suis did not show clinical signs of overt acute disease during the whole experiment, except for mild, transient temperature elevations and low hematocrit values (23–24%) at 10 dpi. Thus, a subclinical infection, resembling the chronic type of the disease seen in the field, developed, with bacteremia first detected by qPCR between days 10 and 27 post M. suis inoculation. It is unclear whether one pig (pig #7) eliminated the infection after day 23 or whether its bacterial loads fell below the detection limit of our assay. This variability in the course of M. suis infection may have important implications related to the host–pathogen interactions identified herein. Nevertheless, the transcriptome profile obtained in this study was performed on samples collected at the first highest level (peak) of bacteremia for each individual animal within the same week, which may mitigate possible biases. The study reveals that pigs subclinically infected with M. suis show altered patterns of gene expression when compared to non-infected pigs, confirming our main hypothesis.
Immune mediated hemolytic anemia (IMHA) has been described in M. suis-infected pigs. Reactive antibodies can be detected as early as 10 days following M. suis inoculation in splenectomized pigs . The low hematocrit values (23–24%) detected 10 dpi in the infected pigs may be a result of a M. suis-induced IMHA. It is unlikely a direct consequence of the pathogen, which was still undetectable by qPCR. This finding warrants further studies in the pathophysiology of M. suis-induced IMHA in non-splenectomized pigs and its impact on pig health.
Pigs #5 and #6 showed consistently high M. suis blood loads (between 108 and 109 organisms/mL of blood) and showed no clinical signs or significant alterations in hematocrit values during these “peaks” of bacteremia. A similar phenomenon (up to 8.23 × 108 organisms/mL of blood, no anemia and negative blood smears) has been observed in a naturally infected, non-splenectomized pig that was followed-up for 80 days by our research group . In that same study, a splenectomized pig showed profound anemia when M. suis loads reached between 1011 and 1012 organisms/mL of blood, which were easily detected on blood smears. On the other hand, anemia and a significant inverse correlation between M. suis loads in the blood and hematocrit values have been described in non-splenectomized, feeder pigs from Germany that have showed up to approximately 109 organisms/mL of blood . The discrepancy in blood loads in anemic and non-anemic animals observed in different studies may be due to differences in M. suis PCR detection protocols and has been previously discussed . Nevertheless, the extent of individual hematological parameter variations of pigs with high M. suis loads (> 108 organisms/mL of blood), and the influence of stress, splenectomy, genetics, co-infections and age on M. suis infection have never been fully explored. Whether or not some of these factors have affected the clinical course of M. suis infection in the studied pigs or in pigs from different studies is unknown.
Due to recent advances in NGS platforms, RNA-Seq studies have become a cost-effective and efficient approach to study global gene expression based on RNA levels within a biological sample. Unfortunately, there is no consensus on which analysis method should be used to guarantee maximum sensitivity and specificity in the detection of DE genes. A recent study that compared the performance of the three most frequently used software tools for RNA-Seq analysis (Cufflinks, DESeq and edgeR) identified that edgeR performs slightly better than the other two methods to detect true positives (i.e. truly DE genes) . Nevertheless, in order to decrease the possibility of false positives, it is still recommended that DE genes should be detected based on their identification in two or more analysis tools . Considering that the objective of our study was to determine if M. suis infection truly causes alterations in pigs’ gene expression profiles, we considered only genes that were identified as DE by two or more methods, thus decreasing the chance of detecting false positive DE genes.
As whole blood (i.e. multiple cell types) was the RNA source of this study, comparisons regarding the detected DE genes with studies using the same system in pigs are limited. However, the significance of each detected DE gene for the pathogenesis of the disease can be evaluated using information available in the literature about the expected function of these genes in each white blood cell (e.g. lymphocytes, monocytes, neutrophils). Interestingly, most DE genes in this study were down-regulated. While several mycoplasma infections are associated with up-regulation of pro-inflammatory genes related with overt disease [42,43,44], our study focused on subclinical infection; thus, more subtle changes were expected. A general suppression of genes related to the innate immune response was observed in the M. suis-infected animals. Accordingly, the expression of TLR8 (Toll like receptor 8) and CD180 genes, the physiological inhibitor of TLR4, was down-regulated in infected pigs. The activation of TLRs by pathogen associated molecular patterns (PAMPs) is critical to the induction of an innate immune response. TLR8 recognizes viral or bacterial single-stranded RNA and is primarily expressed in monocytes/macrophages on intracellular vesicular membranes [45,46,47,48,49]. On the other hand, CD180 expressed in monocytes/macrophages, inhibits TLR4 signaling and thus, low levels can increase TLR4 mediated inflammatory response upon stimulation . TLR4 mRNA expression, however, was not elevated in M. suis infected pigs despite the down-regulation of CD180. Another way by which M. suis may be modulating the innate response is through the down-regulation of cytochrome b-245 beta polypeptide (CYBB) and indoleamine 2,3-dioxygenase (IDO1) genes. CYBB is part of the NADPH oxidase 2 enzyme (NOX2) responsible for the production of superoxide in phagocytes. Low levels of NOX2 may dampen the microbicidal activity of these cells . IDO1 is a tryptophan-degrading enzyme mainly present in antigen presenting cells (APC). It can have opposing effects against intracellular pathogens by either suppressing or increasing their replication . Thus, the down-regulation of TLR8, CD180, CYBB and IDO1 suggests mechanisms by which M. suis may modulate innate immune responses. Considering that many biomedical research projects using pig models focus on the evaluation of innate immune responses, alterations in such genes are likely to interfere with subsequent research results.
The biology of inflammasomes and programmed cell death has been increasingly evaluated in biomedical research projects. The observed down-regulation of the genes that encode caspase-1 (CASP1), caspase-4 (CASP4), apoptosis facilitator Bcl-2-like protein 14 (BCL2L14), and programmed death-ligand 1 (PD-L1 or CD274) suggest the possibility that M. suis is able to interfere with the host’s programmed cell death. Pyroptosis, a type of programmed cell death, is a caspase-1-mediated process that drives the host response to an inflammatory state through the recruitment of the inflammasome’s multiprotein oligomers . Caspase-4 has also been linked to inflammasome signaling involving caspase-1-dependent and caspase-1-independent outcomes [54,55,56]. Intracellular pathogens and LPS have been described to directly interfere with these caspases [56,57,58]. BCL2L14 protein, on the other hand, are linked to the p53-dependent intrinsic apoptosis pathway . Whether or not M. suis may be acting to inhibit the programmed cell death of peripheral blood leukocytes needs to be further explored.
Cytokine production has been studied in porcine models of infectious and metabolic diseases for many years. Numerous commercial assays exist to measure cytokines and chemokines in pig samples. Herein, the down-regulation of several IFN-γ-regulated genes (a type II-IFN) could indicate a mechanism of impaired IFN-γ control of M. suis infection, which could directly interfere with cytokine analyses of pigs in research studies if animals are infected. It is important to note that IFN genes were not differentially expressed in our study, however only one time point was evaluated and the expression of these genes could have varied over the course of infection. The only interleukin gene with altered mRNA levels in M. suis infected pigs was IL15, which was down-regulated. Deficiency of IL15 may enhance erythropoiesis , but at the same time hamper the maintenance of memory CD8+ T cells and NK proliferation . The failure in maintaining proper levels of memory CD8+ T cells may contribute to impaired host protection to re-infection. While cats can be resistant to re-infection with M. haemofelis , this feature is still unclear for M. suis.
The bradykinin, complement, coagulation and fibrinolytic cascades play important roles in inflammation [63, 64]. Of particular interest, C1 inhibitor (C1-INH, encoded by SERPING1 gene), which was down-regulated in M. suis infected pigs, is a major biological inhibitor regulating several proteins of these pathways. The biological activities of C1-INH can be broadly divided into regulation of vascular permeability, coagulation and anti-inflammatory functions . The clinical manifestations of the C1-INH deficiency are dramatically demonstrated in hereditary angioedema (HAE) in humans, where the principle mediator of subcutaneous and submucosal angioedema is excessive bradykinin production . Vascular edema has been suggested to occur in pigs infected with M. suis  and Arthus-like vascular lesions have been reported in cattle infected with M. wenyoni, another hemotrophic mycoplasma . In addition, intravascular thrombosis, superficial hemorrhage (petechia and ecchymosis), prolonged partial thromboplastin and prothrombin times, and decreased platelet count have been observed in splenectomized M. suis-infected pigs [16, 17, 67]. Although other molecules are likely involved in the disturbance of hemostasis in acute infection, the absence of C1-INH could establish a pro-coagulation and pro-thrombotic state in subclinical infections. This could more easily trigger the coagulation cascade leading to clot formation and contributing to angioedema formation by increasing vascular permeability . Although M. suis chronically infected animals show a less severe disease when compared to acutely infected pigs , transcriptional changes observed herein suggest mechanisms that help explain clinical signs and laboratory findings observed during acute attacks, as well as the host’s impaired ability to stop the development of intravascular coagulation. While it is possible that these lesions are related to deficiencies of C1-INH, further studies are needed to clearly define the mechanism involved.
Surprisingly, several DE genes related to endothelial cell physiology and angiogenesis were identified in our analyses. The expression of the gene encoding plasma prolylcarboxypeptidase (PRCP or angiotensinase C), a membrane protein of endothelial and other cells that cleaves a variety of bioactive peptides, including bradykinin, angiotensin II/III, and prekallikrein, was down-regulated at the time point evaluated during subclinical M. suis infection. PRCP activity also has been identified in blood leukocytes in humans, principally monocytes/macrophages . It has been reported that this serine protease acts like an endothelial growth factor. A series of experiments has shown that PRCP deficiency in mice is associated with decreased endothelial cell growth, reduced angiogenesis and impaired wound healing following vascular injury . Given the recently recognized cellular tropism of M. suis for vascular endothelium and believed propagation at this site , the possibility that decreased PRCP may be associated with disruption of endothelial homeostasis is intriguing. The down-regulation of calgranulin C (S100A12) and transglutaminase 1 (TGM1) mRNAs also may have effects related to vascular functions. Acute infection with M. suis in pigs has been shown to result in disturbed endothelial homeostasis and endothelial denudation . While direct effects of M. suis infection on the endothelial cells of pigs were not evaluated in this study, our results suggest that vascular homeostasis may be indirectly affected. It is important to note that the sources of these mRNAs in this study are nucleated cells in the peripheral blood, and its effects on endothelial cells are speculative. Nevertheless, possible changes in vascular homeostasis could have a significant effect on porcine models, such as the use of pigs to study diabetes mellitus, a metabolic disorder known to cause vascular alterations.
In summary, subclinical, possibly chronic, M. suis infection was capable of causing significant transcriptional alterations in host pathways commonly studied in biomedical research (e.g. innate immune response, programmed cell death, vascular homeostasis), and may severely affect scientific results using porcine models. We strongly recommend that pigs be screened for M. suis infection by PCR before being used in research studies.
Lunney JK (2007) Advances in swine biomedical model genomics. Int J Biol Sci 3:179–184
Fan N, Lai L (2013) Genetically modified pig models for human diseases. J Genet Genomics 40:67–73
Gün G, Kues WA (2014) Current progress of genetically engineered pig models for biomedical research. Bioresour Open Access 3:255–264
National Swine Research and Resource Center. http://grants.nih.gov/grants/guide/rfa-files/RFA-RR-03-003.html. Accessed 03 Mar 2015
Groenen MA, Archibald AL, Uenishi H, Tuggle CK, Takeuchi Y, Rothschild MF, Rogel-Gaillard C, Park C, Milan D, Megens HJ, Li S, Larkin DM, Kim H, Frantz LA, Caccamo M, Ahn H, Aken BL, Anselmo A, Anthon C, Auvil L, Badaoui B, Beattie CW, Bendixen C, Berman D, Blecha F, Blomberg J, Bolund L, Bosse M, Botti S, Bujie Z, et al. (2012) Analyses of pig genomes provide insight into porcine demography and evolution. Nature 491:393–398
Verma N, Rettenmeier AW, Schmitz-Spanke S (2011) Recent advances in the use of Sus scrofa (pig) as a model system for proteomic studies. Proteomics 11:776–793
Ramayo-Caldas Y, Mach N, Esteve-Codina A, Corominas J, Castelló A, Ballester M, Estellé J, Ibáñez-Escriche N, Fernández AI, Pérez-Enciso M, Folch JM (2012) Liver transcriptome profile in pigs with extreme phenotypes of intramuscular fatty acid composition. BMC Genomics 13:547
Shan L, Wu Q, Li Y, Shang H, Guo K, Wu J, Wei H, Zhao J, Yu J, Li MH (2014) Transcriptome profiling identifies differentially expressed genes in postnatal developing pituitary gland of miniature pig. DNA Res 21:207–216
Puig-Oliveras A, Ramayo-Caldas Y, Corominas J, Estellé J, Pérez-Montarelo D, Hudson NJ, Casellas J, Folch JM, Ballester M (2014) Differences in muscle transcriptome among pigs phenotypically extreme for fatty acid composition. PLoS One 9:e99720
Ropka-Molik K, Zukowski K, Eckert R, Gurgul A, Piórkowska K, Oczkowicz M (2014) Comprehensive analysis of the whole transcriptomes from two different pig breeds using RNA-Seq method. Anim Genet 45:674–684
Sodhi SS, Park WC, Ghosh M, Kim JN, Sharma N, Shin KY, Cho IC, Ryu YC, Oh SJ, Kim SH, Song KD, Hong SP, Cho SA, Kim HB, Jeong DK (2014) Comparative transcriptomic analysis to identify differentially expressed genes in fat tissue of adult Berkshire and Jeju Native Pig using RNA-seq. Mol Biol Rep 41:6305–6315
Xing K, Zhu F, Zhai L, Liu H, Wang Z, Hou Z, Wang C (2014) The liver transcriptome of two full-sibling Songliao black pigs with extreme differences in backfat thickness. J Anim Sci Biotechnol 5:32
Gaur U, Xiong Y, Luo Q, Yuan F, Wu H, Qiao M, Wimmers K, Li K, Mei SQ, Liu GS (2014) Breed-specific transcriptome response of spleen from 6–8 week old piglet after infection with Streptococcus suis type 2. Mol Biol Rep 41:7865–7873
Schroyen M, Tuggle CK (2015) Current transcriptomics in pig immunity research. Mamm Genome 26:1–20
Groebel K, Hoelzle K, Wittenbrink MM, Ziegler U, Hoelzle LE (2009) Mycoplasma suis invades porcine erythrocytes. Infect Immun 77:576–584
Sokoli A, Groebel K, Hoelzle K, Amselgruber WM, Mateos JM, Schneider MK, Ziegler U, Felder KM, Hoelzle LE (2013) Mycoplasma suis infection results endothelial cell damage and activation: new insight into the cell tropism and pathogenicity of hemotrophic mycoplasma. Vet Res 44:6
Stadler J, Jannasch C, Mack SL, Dietz S, Zöls S, Ritzmann M, Hoelzle K, Hoelzle LE (2014) Clinical and haematological characterization of Mycoplasma suis infections in splenectomised and non-splenectomised pigs. Vet Microbiol 172:294–300
Zachary JF, Smith AR (1985) Experimental porcine eperythrozoonosis: T-lymphocyte suppression and misdirected immune responses. Am J Vet Res 46:821–830
Gwaltney SM (1995) Eperythrozoon suis infections in pigs: clinical syndromes and diagnosis. J Swine Health Prod 3:25–27
Guimaraes AM, Vieira RF, Poletto R, Vemulapalli R, Santos AP, de Moraes W, Cubas ZS, Santos LC, Marchant-Forde JN, Timenetsky J, Biondo AW, Messick JB (2011) A quantitative TaqMan PCR assay for the detection of Mycoplasma suis. J Appl Microbiol 111:417–425
Guimaraes AM, Biondo AW, Lara AC, Messick JB (2007) Exploratory study of Mycoplasma suis (Eperythrozoon suis) on four commercial pig farms in southern Brazil. Vet Rec 160:50–53
Ritzmann M, Grimm J, Heinritzi K, Hoelzle K, Hoelzle LE (2009) Prevalence of Mycoplasma suis in slaughter pigs, with correlation of PCR results to hematological findings. Vet Microbiol 133:84–91
Yuan CL, Liang AB, Yao CB, Yang ZB, Zhu JG, Cui L, Yu F, Zhu NY, Yang XW, Hua XG (2009) Prevalence of Mycoplasma suis (Eperythrozoon suis) infection in swine and swine-farm workers in Shanghai, China. Am J Vet Res 70:890–894
Baker HJ, Cassell GH, Lindsey JR (1971) Research complications due to Haemobartonella and Eperythrozoon infections in experimental animals. Am J Pathol 64:625–652
Contamin H, Michel JC (1999) Haemobartonellosis in squirrel monkeys (Saimiri sciureus): antagonism between Haemobartonella sp. and experimental Plasmodium falciparum malaria. Exp Parasitol 91:297–305
Kemming GI, Messick JB, Enders G, Boros M, Lorenz B, Muenzing S, Kisch-Wedel H, Mueller W, Hahmann-Mueller A, Messmer K, Thein E (2004) Mycoplasma haemocanis infection—a kennel disease? Comp Med 54:404–409
Nascimento NC, Santos AP, Guimaraes AMS, Breur GJ, Messick JB (2011) PCR detection of Mycoplasma ovis and ‘Candidatus Mycoplasma haemovis’ in sheep. Vet Clin Pathol 40:598
Dent BT, Stevens KA, Korvick DL, Clymer JW (2013) Mycoplasma suis infection in pigs after splenectomy. Lab Anim 42:125–128
Messick JB, Cooper SK, Huntley M (1999) Development and evaluation of a polymerase chain reaction assay using the 16S rRNA gene for detection of Eperythrozoon suis infection. J Vet Diagn Invest 11:229–236
Guimaraes AM, Santos AP, Timenetsky J, Bower LP, Strait E, Messick JB (2014) Identification of Mycoplasma suis antigens and development of a multiplex microbead immunoassay. J Vet Diagn Invest 26:203–212
Andrews S (2010) A quality control tool for high throughput sequence data. FastQC. http://www.bioinformatics.babraham.ac.uk/projects/fastqc. Accessed 9 Oct 2014
Trapnell C, Pachter L, Salzberg SL (2009) TopHat: discovering splice junctions with RNA-Seq. Bioinformatics 25:1105–1111
Robinson MD, McCarthy DJ, Smyth GK (2010) edgeR: a bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 26:139–140
Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, Salzberg SL, Wold BJ, Pachter L (2010) Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol 28:511–515
Untergasser A, Cutcutache I, Koressaar T, Ye J, Faircloth BC, Remm M, Rozen SG (2012) Primer3—new capabilities and interfaces. Nucleic Acids Res 40:e115
Pfaffl MW (2001) A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res 29:e45
Fu Y, Quan R, Zhang H, Hou J, Tang J, Feng WH (2012) Porcine reproductive and respiratory syndrome virus induces interleukin-15 through the NF-κB signaling pathway. J Virol 86:7625–7636
Jensen LJ, Kuhn M, Stark M, Chaffron S, Creevey C, Muller J, Doerks T, Julien P, Roth A, Simonovic M, Bork P, von Mering C (2009) STRING 8–a global view on proteins and their functional interactions in 630 organisms. Nucleic Acids Res 37:D412–416
Conesa A, Götz S, García-Gómez JM, Terol J, Talón M, Robles M (2005) Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics 21:3674–3676
Felder KM, Hoelzle K, Heinritzi K, Ritzmann M, Hoelzle LE (2010) Antibodies to actin in autoimmune haemolytic anaemia. BMC Vet Res 6:18
Zhang ZH, Jhaveri DJ, Marshall VM, Bauer DC, Edson J, Narayanan RK, Robinson GJ, Lundberg AE, Bartlett PF, Wray NR, Zhao QY (2014) A comparative study of techniques for differential expression analysis on RNA-Seq data. PLoS One 9:e103207
Bin L, Luping D, Bing S, Zhengyu Y, Maojun L, Zhixin F, Yanna W, Haiyan W, Guoqing S, Kongwang H (2014) Transcription analysis of the porcine alveolar macrophage response to Mycoplasma hyopneumoniae. PLoS One 9:e101968
McGowin CL, Radtke AL, Abraham K, Martin DH, Herbst-Kralovetz M (2013) Mycoplasma genitalium infection activates cellular host defense and inflammation pathways in a 3-dimensional human endocervical epithelial cell model. J Infect Dis 207:1857–1868
Hopfe M, Deenen R, Degrandi D, Köhrer K, Henrich B (2013) Host cell responses to persistent mycoplasmas–different stages in infection of HeLa cells with Mycoplasma hominis. PLoS One 8:e54219
Tanji H, Ohto U, Shibata T, Taoka M, Yamauchi Y, Isobe T, Miyake K, Shimizu T (2015) Toll-like receptor 8 senses degradation products of single-stranded RNA. Nat Struct Mol Biol 22:109–115
Cervantes JL, Weinerman B, Basole C, Salazar JC (2012) TLR8: the forgotten relative revindicated. Cell Mol Immunol 9:434–438
Davila S, Hibberd ML, Hari Dass R, Wong HE, Sahiratmadja E, Bonnard C, Alisjahbana B, Szeszko JS, Balabanova Y, Drobniewski F, van Crevel R, van de Vosse E, Nejentsev S, Ottenhoff TH, Seielstad M (2008) Genetic association and expression studies indicate a role of Toll-like receptor 8 in pulmonary tuberculosis. PLoS Genet 4:e1000218
Gantier MP, Irving AT, Kaparakis-Liaskos M, Xu D, Evans VA, Cameron PU, Bourne JA, Ferrero RL, John M, Behlke MA, Williams BR (2010) Genetic modulation of TLR8 response following bacterial phagocytosis. Hum Mutat 31:1069–1079
Cervantes JL, La Vake CJ, Weinerman B, Luu S, O’Connell C, Verardi PH, Salazar JC (2013) Human TLR8 is activated upon recognition of Borrelia burgdorferi RNA in the phagosome of human monocytes. J Leukoc Biol 94:1231–1241
Divanovic S, Trompette A, Atabani SF, Madan R, Golenbock DT, Visintin A, Finberg RW, Tarakhovsky A, Vogel SN, Belkaid Y, Kurt-Jones EA, Karp CL (2005) Negative regulation of Toll-like receptor 4 signaling by the Toll-like receptor homolog RP105. Nat Immunol 6:571–578
El-Benna J, Dang PM, Gougerot-Pocidalo MA (2008) Priming of the neutrophil NADPH oxidase activation: role of p47phox phosphorylation and NOX2 mobilization to the plasma membrane. Semin Immunopathol 30:279–289
Munn DH, Mellor AL (2013) Indoleamine 2,3 dioxygenase and metabolic control of immune responses. Trends Immunol 34:137–143
Schroder K, Tschopp J (2010) The inflammasomes. Cell 140:821–832
Rathinam VA, Vanaja SK, Fitzgerald KA (2012) Regulation of inflammasome signaling. Nat Immunol 13:333–342
Kim J, Ahn H, Woo HM, Lee E, Lee GS (2014) Characterization of porcine NLRP3 inflammasome activation and its upstream mechanism. Vet Res Commun 38:193–200
Yang J, Zhao Y, Shao F (2015) Non-canonical activation of inflammatory caspases by cytosolic LPS in innate immunity. Curr Opin Immunol 32:78–83
Hersh D, Monack DM, Smith MR, Ghori N, Falkow S, Zychlinsky A (1999) The Salmonella invasin SipB induces macrophage apoptosis by binding to caspase-1. Proc Natl Acad Sci U S A 96:2396–2401
Hilbi H, Moss JE, Hersh D, Chen Y, Arondel J, Banerjee S, Flavell RA, Yuan J, Sansonetti PJ, Zychlinsky A (1998) Shigella-induced apoptosis is dependent on caspase-1 which binds to IpaB. J Biol Chem 273:32895–32900
Chipuk JE, Green DR (2006) Dissecting p53-dependent apoptosis. Cell Death Differ 13:994–1002
Mullarky IK, Szaba FM, Kummer LW, Wilhelm LB, Parent MA, Johnson LL, Smiley ST (2007) Gamma interferon suppresses erythropoiesis via interleukin-15. Infect Immun 75:2630–2633
Di Sabatino A, Calarota SA, Vidali F, Macdonald TT, Corazza GR (2011) Role of IL-15 in immune-mediated and infectious diseases. Cytokine Growth Factor Rev 22:19–33
Hicks CA, Willi B, Riond B, Novacco M, Meli ML, Stokes CR, Helps CR, Hofmann-Lehmann R, Tasker S (2015) Protective immunity against infection with Mycoplasma haemofelis. Clin Vaccine Immunol 22:108–118
Kaplan AP, Ghebrehiwet B (2010) The plasma bradykinin-forming pathways and its interrelationships with complement. Mol Immunol 47:2161–2169
Esmon CT (2005) The interactions between inflammation and coagulation. Br J Haematol 131:417–430
Davis AE 3rd, Mejia P, Lu F (2008) Biological activities of C1 inhibitor. Mol Immunol 45:4057–4063
Montes AJ, Wolfe DF, Welles EG, Tyler JW, Tepe E (1994) Infertility associated with Eperythrozoon wenyonii infection in a bull. J Am Vet Med Assoc 204:261–263
Plank G, Heinritzi K (1990) Disseminated intravascular coagulation in eperythrozoonosis of swine. Berl Munch Tierarztl Wochenschr 103:13–18
Cugno M, Cicardi M, Bottasso B, Coppola R, Paonessa R, Mannucci PM, Agostoni A (1997) Activation of the coagulation cascade in C1-inhibitor deficiencies. Blood 89:3213–3218
Kumamoto K, Stewart TA, Johnson AR, Erdös EG (1981) Prolylcarboxypeptidase (angiotensinase C) in human lung and cultured cells. J Clin Invest 67:210–215
Adams GN, Stavrou EX, Fang C, Merkulova A, Alaiti MA, Nakajima K, Morooka T, Merkulov S, Larusch GA, Simon DI, Jain MK, Schmaier AH (2013) Prolylcarboxypeptidase promotes angiogenesis and vascular repair. Blood 122:1522–1531
The authors declare that they have no competing interests.
Conceived and designed the experiments: NCN, AMSG, APS, YC, JBM. Performed the experiments: NCN, APS, YC, JBM. Analyzed the data: NCN, AMSG, APS, YC, LMM, JBM. Contributed reagents/materials/analysis tools: NCN, AMSG, APS, YC, LMM, JBM. Wrote the paper: NCN, AMSG, APS, YC, LMM, JBM. All authors read and approved the final manuscript.
Authors are very grateful to the Purdue Genomics Core Facility for library construction and RNA sequencing, and the Bioinformatics Core at Purdue University, especially Dr Jyothi Thimmapuram and Ketaki Bhide, for the RNA-Seq analysis and bioinformatics support.
Ethics approval and consent to participate
Postdoctoral Fellowship for Naila C. do Nascimento was provided by Brazilian Science Agency CNPq (National Council for Scientific and Technological Development). PhD Fellowship for Ana M. S. Guimaraes was provided by Brazilian Ministry of Education through Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) and Fulbright Commission. CAPES-Fulbright Program, ID 167307/6. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. Support was provided from Hatch Act Formula Grant, Project No. IND020395.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Results of Mycoplasma suis microbead immunoassay (MIA). Median Fluorescence Intensity (MFI) of Mycoplasma suis-infected and sham-inoculated pigs as detected in serum samples collected along the experiment. Serum antibodies were measured using recombinant GrpE (rGrpE) protein in a MIA . Sham-inoculated group (control group): pig #1, pig #2 and pig #3. M. suis-inoculated group (infected group): pig #5, pig #6, pig #7. Day 0 corresponds to inoculation day. MFI cut-off was 74.07. It is important to note that a high cut-off may have been encountered due to the low number of animals to determine a cut-off, which increases the standard deviation. MFI values in bold correspond to positive results based on cut-off values for each group. ND: not done.
Results of RNA-Seq analyses. List of differentially expressed (DE) genes identified in whole blood of Mycoplasma suis-infected pigs (when compared to sham-infected group) separated by all three methods of RNA-Seq analyses (DESeq2, edgeR, Cufflinks). First tab describes all DE genes (p-value ≤ 0.1) that were detected in all three methods, in two methods, followed by each individual method. The following tabs report the statistics of each evaluated gene for each method. Highlighted in grey are statistically significant DE genes (dark grey: p-value ≤ 0.05; light grey p-value ≤ 0.1).
Relative expression (RE) of differentially expressed (DE) in pigs infected with Mycoplasma suis validated by qRT-PCR. Relative expression profile (GAPDH as reference control gene) of DE genes identified in the M. suis-infected pigs compared to the control group (non-infected). Validated genes are: ISG15 ubiquitin-like modifier (ISG15), interleukin 22 receptor, alpha 2 (IL22RA2), CD274 molecule (CD274), BCL2-like 14 (apoptosis facilitator) (BCL2), chemokine (C-C motif) receptor 5 (CCR5), CD180 molecule (CD180), allograft inflammatory factor 1 (AIF1), interleukin 15 (IL15), protein tyrosine phosphatase, receptor type, O (PTPRO), and Toll-like receptor 8 (TLR8). Black dots in each graphic represent the RE of a respective gene for each pig; blue vertical lines represent the standard deviation in each group. * p-value ≤ 0.1, and ** p-value ≤ 0.05.
Gene ontology (GO) categories of all 55 differentially expressed (DE) genes of Mycoplasma suis infected pigs. First tab describes raw results of GO analysis, with all DE genes and their correspondent GO categories (and functions). Second tab shows the identified functions of the “biological process” category of all 55 DE genes and the number of reported genes for each function. Functional categories with the largest number of genes are highlighted in yellow. GO analysis was performed using Blast2Go software .
About this article
Cite this article
do Nascimento, N.C., Guimaraes, A.M.S., dos Santos, A.P. et al. RNA-Seq based transcriptome of whole blood from immunocompetent pigs (Sus scrofa) experimentally infected with Mycoplasma suis strain Illinois. Vet Res 49, 49 (2018). https://doi.org/10.1186/s13567-018-0546-6
- Suis Infection
- Differentially Expressed (DE)
- Infected Pigs
- Hereditary Angioedema
- SERPING1 Gene