Comparative Transcriptomic Analyses of Haemophilus parasuis Reveal Differently Expressed Genes among Strains with Different Virulence Degrees

Haemophilus parasuis is commonly found in the upper respiratory tract of the pigs. Some isolates of H. parasuis can lead to both pneumonia and Glässer's disease of pigs with severe clinical symptoms. The virulence-associated genes for the various degrees of virulence observed in H. parasuis remains poorly understood. In the present study, we identified the differentially expressed genes between YK1603 (non-virulent strain) and XM1602 (moderately virulent strain) or CY1201 (highly virulent strain) of H. parasuis using Illumina sequencing technique. In comparison to YK1603, a total of 195 genes were significantly changed in CY1201, of which 71 genes were up-regulated and 124 genes were down-regulated, whereas 705 genes were significantly changed in XM1602, of which 415 genes were up-regulated and 290 genes were down-regulated. The enriched analysis of Gene Ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways on the differentially expressed genes showed that both enriched main GO terms and KEGG pathways appear to be different between the two kinds of comparision: CY1201 versus YK1603, and XM1602 versus YK1603. Based on real-time PCR technique, on the whole, it was confirmed that the expression of ten genes: lpxL, tbpB, kdtA, waaQ, oapA, napA, ptsH, mmsA, lpxM, and lpxB were agreement with the findings in Illumina sequencing analysis. These identified genes might participate in the regulation of a wide range of biological process involved in virulence of H. parasuis, such as phosphotransferase system and ABC transporters. Our results from this study provide a new way to gain insight into the virulent mechanisms of H. parasuis.