- Research article
- Open Access
RNA-Seq based transcriptome of whole blood from immunocompetent pigs (Sus scrofa) experimentally infected with Mycoplasma suis strain Illinois
- Naíla C. do Nascimento†1Email author,
- Ana M. S. Guimaraes†1, 4View ORCID ID profile,
- Andrea P. dos Santos†1,
- Yuefeng Chu1, 5,
- Lucas M. Marques2, 3 and
- Joanne B. Messick1Email author
© The Author(s) 2018
- Received: 7 December 2017
- Accepted: 12 March 2018
- Published: 18 June 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–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–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–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.
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
Primers used for validation of differentially expressed (DE) genes in Mycoplasma suis -infected pigs by quantitative reverse-transcriptase polymerase chain reaction (qRT-PCR)
Primer forward (5′–3′)
Primer reverse (5′–3′)
Amplicon length (bp)
ISG15 ubiquitin-like modifier
Interleukin 22 receptor, alpha 2
BCL2-like 14 (apoptosis facilitator)
Chemokine (C-C motif) receptor 5
Allograft inflammatory factor 1
Protein tyrosine phosphatase, receptor type, O
Toll-like receptor 8
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
RNA-Seq results and validation by qRT-PCR
Differentially expressed (DE) genes identified in Mycoplasma suis -infected pigs tested by quantitative reverse-transcriptase polymerase chain reaction (qRT-PCR)
DE up-regulated genes
ISG15 ubiquitin-like modifier
DE down-regulated genes
Interleukin 22 receptor, alpha 2
BCL2-like 14 (apoptosis facilitator)
chemokine (C–C motif) receptor 5
Allograft inflammatory factor 1
Protein tyrosine phosphatase, receptor type, O
Toll-like receptor 8
Network analysis of the relationships between proteins encoded by DE genes
Gene ontology analysis of differentially expressed (DE) genes of M. suis infected pigs
Four main gene ontology (GO) functions, identified as biological processes, of differentially expressed (DE) genes of Mycoplasma suis infected pigs
Number of DE genes
CCR5, CYBB, IDO1, TLR8, CASP4, CASP1, CXCL10, AIF1, CD180
CCR5, TLR8, CCL8, CXCL10, PDL1
CASP4, BCL2L14, CASP1, AIF1, MERTK
C2, C1-INH, CASP4, CASP1, LAP3, PRCP
Cellular response to lipopolysaccharide
CCR5, HAVCR2, CXCL10, CD180
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–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–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–56]. Intracellular pathogens and LPS have been described to directly interfere with these caspases [56–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.
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.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Lunney JK (2007) Advances in swine biomedical model genomics. Int J Biol Sci 3:179–184View ArticlePubMedPubMed CentralGoogle Scholar
- Fan N, Lai L (2013) Genetically modified pig models for human diseases. J Genet Genomics 40:67–73View ArticlePubMedGoogle Scholar
- Gün G, Kues WA (2014) Current progress of genetically engineered pig models for biomedical research. Bioresour Open Access 3:255–264View ArticleGoogle Scholar
- 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–398View ArticlePubMedPubMed CentralGoogle Scholar
- 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–793View ArticlePubMedGoogle Scholar
- 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:547View ArticlePubMedPubMed CentralGoogle Scholar
- 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–216View ArticlePubMedGoogle Scholar
- 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:e99720View ArticlePubMedPubMed CentralGoogle Scholar
- 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–684View ArticlePubMedGoogle Scholar
- 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–6315View ArticlePubMedGoogle Scholar
- 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:32View ArticlePubMedPubMed CentralGoogle Scholar
- 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–7873View ArticlePubMedGoogle Scholar
- Schroyen M, Tuggle CK (2015) Current transcriptomics in pig immunity research. Mamm Genome 26:1–20View ArticlePubMedGoogle Scholar
- Groebel K, Hoelzle K, Wittenbrink MM, Ziegler U, Hoelzle LE (2009) Mycoplasma suis invades porcine erythrocytes. Infect Immun 77:576–584View ArticlePubMedGoogle Scholar
- 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:6View ArticlePubMedPubMed CentralGoogle Scholar
- 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–300View ArticlePubMedGoogle Scholar
- Zachary JF, Smith AR (1985) Experimental porcine eperythrozoonosis: T-lymphocyte suppression and misdirected immune responses. Am J Vet Res 46:821–830PubMedGoogle Scholar
- Gwaltney SM (1995) Eperythrozoon suis infections in pigs: clinical syndromes and diagnosis. J Swine Health Prod 3:25–27Google Scholar
- 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–425View ArticlePubMedGoogle Scholar
- 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–53View ArticlePubMedGoogle Scholar
- 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–91View ArticlePubMedGoogle Scholar
- 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–894View ArticlePubMedGoogle Scholar
- Baker HJ, Cassell GH, Lindsey JR (1971) Research complications due to Haemobartonella and Eperythrozoon infections in experimental animals. Am J Pathol 64:625–652PubMedPubMed CentralGoogle Scholar
- Contamin H, Michel JC (1999) Haemobartonellosis in squirrel monkeys (Saimiri sciureus): antagonism between Haemobartonella sp. and experimental Plasmodium falciparum malaria. Exp Parasitol 91:297–305View ArticlePubMedGoogle Scholar
- 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–409PubMedGoogle Scholar
- 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:598Google Scholar
- Dent BT, Stevens KA, Korvick DL, Clymer JW (2013) Mycoplasma suis infection in pigs after splenectomy. Lab Anim 42:125–128View ArticleGoogle Scholar
- 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–236View ArticlePubMedGoogle Scholar
- 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–212View ArticlePubMedGoogle Scholar
- 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–1111View ArticlePubMedPubMed CentralGoogle Scholar
- Robinson MD, McCarthy DJ, Smyth GK (2010) edgeR: a bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 26:139–140View ArticlePubMedGoogle Scholar
- 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–515View ArticlePubMedPubMed CentralGoogle Scholar
- Untergasser A, Cutcutache I, Koressaar T, Ye J, Faircloth BC, Remm M, Rozen SG (2012) Primer3—new capabilities and interfaces. Nucleic Acids Res 40:e115View ArticlePubMedPubMed CentralGoogle Scholar
- Pfaffl MW (2001) A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res 29:e45View ArticlePubMedPubMed CentralGoogle Scholar
- 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–7636View ArticlePubMedPubMed CentralGoogle Scholar
- 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–416View ArticlePubMedGoogle Scholar
- 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–3676View ArticlePubMedGoogle Scholar
- Felder KM, Hoelzle K, Heinritzi K, Ritzmann M, Hoelzle LE (2010) Antibodies to actin in autoimmune haemolytic anaemia. BMC Vet Res 6:18View ArticlePubMedPubMed CentralGoogle Scholar
- 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:e103207View ArticlePubMedPubMed CentralGoogle Scholar
- 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:e101968View ArticlePubMedPubMed CentralGoogle Scholar
- 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–1868View ArticlePubMedPubMed CentralGoogle Scholar
- 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:e54219View ArticlePubMedPubMed CentralGoogle Scholar
- 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–115View ArticlePubMedGoogle Scholar
- Cervantes JL, Weinerman B, Basole C, Salazar JC (2012) TLR8: the forgotten relative revindicated. Cell Mol Immunol 9:434–438View ArticlePubMedPubMed CentralGoogle Scholar
- 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:e1000218View ArticlePubMedPubMed CentralGoogle Scholar
- 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–1079View ArticlePubMedGoogle Scholar
- 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–1241View ArticlePubMedPubMed CentralGoogle Scholar
- 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–578View ArticlePubMedPubMed CentralGoogle Scholar
- 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–289View ArticlePubMedGoogle Scholar
- Munn DH, Mellor AL (2013) Indoleamine 2,3 dioxygenase and metabolic control of immune responses. Trends Immunol 34:137–143View ArticlePubMedGoogle Scholar
- Schroder K, Tschopp J (2010) The inflammasomes. Cell 140:821–832View ArticlePubMedGoogle Scholar
- Rathinam VA, Vanaja SK, Fitzgerald KA (2012) Regulation of inflammasome signaling. Nat Immunol 13:333–342View ArticlePubMedPubMed CentralGoogle Scholar
- 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–200View ArticlePubMedGoogle Scholar
- Yang J, Zhao Y, Shao F (2015) Non-canonical activation of inflammatory caspases by cytosolic LPS in innate immunity. Curr Opin Immunol 32:78–83View ArticlePubMedGoogle Scholar
- 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–2401View ArticlePubMedPubMed CentralGoogle Scholar
- 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–32900View ArticlePubMedGoogle Scholar
- Chipuk JE, Green DR (2006) Dissecting p53-dependent apoptosis. Cell Death Differ 13:994–1002View ArticlePubMedGoogle Scholar
- 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–2633View ArticlePubMedPubMed CentralGoogle Scholar
- 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–33View ArticlePubMedGoogle Scholar
- 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–118View ArticlePubMedGoogle Scholar
- Kaplan AP, Ghebrehiwet B (2010) The plasma bradykinin-forming pathways and its interrelationships with complement. Mol Immunol 47:2161–2169View ArticlePubMedGoogle Scholar
- Esmon CT (2005) The interactions between inflammation and coagulation. Br J Haematol 131:417–430View ArticlePubMedGoogle Scholar
- Davis AE 3rd, Mejia P, Lu F (2008) Biological activities of C1 inhibitor. Mol Immunol 45:4057–4063View ArticlePubMedPubMed CentralGoogle Scholar
- 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–263PubMedGoogle Scholar
- Plank G, Heinritzi K (1990) Disseminated intravascular coagulation in eperythrozoonosis of swine. Berl Munch Tierarztl Wochenschr 103:13–18PubMedGoogle Scholar
- 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–3218PubMedGoogle Scholar
- Kumamoto K, Stewart TA, Johnson AR, Erdös EG (1981) Prolylcarboxypeptidase (angiotensinase C) in human lung and cultured cells. J Clin Invest 67:210–215View ArticlePubMedPubMed CentralGoogle Scholar
- 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–1531View ArticlePubMedPubMed CentralGoogle Scholar