Comparative Analysis of Immune Response Genes Induced by a Virulent or Attenuated Strain of .

RNA-seq technology has been widely used for the characterization of the transcriptome profile induced by several diseases in both humans and animals. In the present study, RNA-seq was used to identify the differential expression of genes associated with the immune response in cattle infected with two different strains of , both derived from the same Mexican field isolate, which exhibit distinct phenotypic characteristics: the virulent strain, capable of producing acute clinical signs, and the attenuated strain, capable of stimulating a protective immune response when used as an immunogen with an efficacy greater than 80%. The differential gene expression analysis performed revealed a total of 620 differentially expressed genes (DEGs). However, the intersection of the edgeR and DESeq2 programs used in the bioinformatics analysis only identified 247 DEGs, of which 108 genes were enriched to be closely correlated with the bovine immune response based on gene ontology terms; most of the DEGs obtained encode proteins associated with the major histocompatibility complex, immunoglobulins, and T-cell surface receptors. The infection caused by the attenuated strain induced higher transcription of immune response genes compared to the infection caused by the virulent strain; nonetheless, in both infections, a greater down-regulation than up-regulation was observed. Different immunoglobulin-associated genes were found to be up-regulated in the group inoculated with the attenuated strain, whereas these were down-regulated in the virulent strain-inoculated group. In addition, an up-regulation of the HSPA6, CD163, and SLC11a1 genes was observed in the group inoculated with the virulent strain, previously reported in other Apicomplexan infections. The findings provide relevant information that could contribute to clarifying the immune response associated with an acute bovine babesiosis infection by .