Abstract:
BACKGROUND : Babesia rossi is a leading cause of morbidity and mortality among the canine population of sub-
Saharan Africa, but pathogenesis remains poorly understood. Previous studies of B. rossi infection were derived from
clinical cases, in which neither the onset of infection nor the infectious inoculum was known. Here, we performed
controlled B. rossi inoculations in canines and evaluated disease progression through clinical tests and whole blood
transcriptomic profiling.
RESULTS : Two subjects were administered a low inoculum (104 parasites) while three received a high (108 parasites).
Subjects were monitored for 8 consecutive days; anti-parasite treatment with diminazene aceturate was
administered on day 4. Blood was drawn prior to inoculation as well as every experimental day for assessment of
clinical parameters and transcriptomic profiles. The model recapitulated natural disease manifestations including
anemia, acidosis, inflammation and behavioral changes. Rate of disease onset and clinical severity were proportional
to the inoculum. To analyze the temporal dynamics of the transcriptomic host response, we sequenced mRNA
extracted from whole blood drawn on days 0, 1, 3, 4, 6, and 8. Differential gene expression, hierarchical clustering,
and pathway enrichment analyses identified genes and pathways involved in response to hemolysis, metabolic
changes, and several arms of the immune response including innate immunity, adaptive immunity, and response to
viral infection.
CONCLUSIONS : This work comprehensively characterizes the clinical and transcriptomic progression of B. rossi
infection in canines, thus establishing a large mammalian model of severe hemoprotozoal disease to facilitate the
study of host-parasite biology and in which to test novel anti-disease therapeutics. The knowledge gained from the
study of B. rossi in canines will not only improve our understanding of this emerging infectious disease threat in
domestic dogs, but also provide insight into the pathobiology of human diseases caused by Babesia and
Plasmodium species.
Description:
Additional file 1: Table S1. A summary of MultiQC data. Shown are
the number of reads mapped (in millions), percent of these reads aligned
to the canine reference genome, and percent GC content of these reads.
Samples corresponding to low inoculum day 4 and high inoculum days
3 and 4 were not included in the mean.
Additional file 2: Table S2. DEGs on every day in both inoculum
cohorts. Genes that were differentially expressed on every experimental
day in both cohorts (except day 1 in the low inoculum) are shown, along
with their corresponding cluster.
Additional file 3: Table S3. Pearson correlation values of erythroblast
transcript trajectory with reticulocyte count. The expression level of each
gene (in CPM) was correlated with the reticulocyte count (109/L) on each
day. This allowed for correlation based on expression trajectory through
time.
Additional file 4: Fig. S1. Top 10 Gene Ontology (GO)- and Reactomeannoted
pathways by FDR on each day in each inoculum cohort (A –
low; B – high). x-axis: percentage of genes in pathway that were identified
as differentially expressed on the respective day; y-axis: pathway
name; point color: absolute number of genes associated with the pathway
that were identified as differentially expressed on the corresponding
day. Note that color does not indiciate up- or down-reulgation of genes
in pathway.