- Research article
- Open Access
Temporal analysis of mRNA expression profiles in Orientia infected C3HeB/FeJ mouse
BMC Microbiology volume 20, Article number: 3 (2020)
Scrub typhus causes up to 35% mortality if left untreated. One billion people living in the endemic regions are at risk. In spite of its heavy disease burden in some of the most populated areas in the world, there is no vaccine available. Although the disease can be effectively treated by proper antibiotics, timely and accurate diagnosis remains a challenge. Orientia tsutsugamushi infects a variety of mammalian cells in vitro and replicates in the cytoplasm of the infected cells. Microarray analysis has been used extensively to study host-pathogen interactions in in vitro models to understand pathogenesis. However there is a lack of in vivo studies.
In this study, C3HeB/FeJ (C3H) mice were infected by O. tsutsugamushi via the intraperitoneal route and monitored gene expression at 10 different time points post infection. We observed two distinct types of expression profiles in the genes that we analyzed. There are two valleys (4–18 h and 2–4 days) with low number of differentially expressed genes (DEG) with three peaks with high number of DEG at 2 h, 1-day and 7-day post infection. Further analysis revealed that pathways like complement and coagulation cascade, and blood clotting cascade pathways showed significant global changes throughout entire time course. Real time quantitative Polymerase Chain Reaction (RT-qPCR) confirmed the change of expression for genes involved in complement and coagulation cascade. These results suggested dynamic regulation of the complement and coagulation cascades throughout most of the time post infection while some other specific pathways, such as fatty acid metabolism and tryptophan metabolism, are turned on or off at certain times post infection.
The findings highlight the complex interconnection among all different biological pathways. It is conceivable that specific pathways such as cell growth control and cell development in the host are affected by Orientia in the initial phase of infection for Orientia to grow intracellularly. Once Orientia is replicating successfully inside the host as infection progresses, the infection could activate pathways involved in cellular immune responses to defend for host cell survival and try to eliminate the pathogen.
Scrub typhus is a febrile illness caused by an obligate intracellular bacterium, Orientia tsutaugamushi. The disease is transmitted by the bite of the larvae of various species of trombiculid mites. Humans are the accidental hosts while the natural hosts for these mites are rodents. Scrub typhus is endemic in Asia-Pacific region and up to 23% of all febrile illnesses can be attributed to scrub typhus [1, 2]. Scrub typhus is endemic within a “tsutsugamushi” triangle area that is about 13 million square kilometers which includes Pakistan, India and Nepal to the west, Japan to the east, southeastern Siberia, China, and Korea to the north and Indonesia, Philippines, northern Australia and the intervening Pacific islands to the south [3,4,5,6,7,8,9,10,11,12]. This area is one of the most populated areas and an estimated 1 billion people living in this area are at risk. Recent evidence has suggested that the disease is expanding beyond the traditional Asia-Pacific area to include the Middle East, South America and Africa [2, 13,14,15]. The disease is characterized by fever, rash, eschar, pneumonitis, meningitis, and in some severe cases, patients can have disseminated intravascular coagulation that may lead to circulatory failure . In most scrub typhus patients, the symptoms are relatively mild, a timely, accurate and differential diagnosis can be challenging . Therefore, although the disease can be effectively treated by doxycycline, oftentimes the patients are left untreated due to the lack of accurate diagnosis. The mortality in this untreated population can be as high as 35% [1, 17]. Furthermore, cases of ineffective doxycycline treatment have been documented in Thailand and India, suggesting the potential rise of antibiotic resistant strains of O. tsutsugamushi in endemic areas [18, 19].
Orientia is extremely labile and can infect a variety of mammalian cells ex vivo, including endothelial cells , dendritic cells , monocytes  and polymorphonuclear leukocytes . Once they are in mammalian cells, it is believed to replicate in the cytoplasm of the infected cells. This is supported by previous studies showing that O. tsutsugamushi induces phagocytosis of host cells and then escapes from the phagosome in 30 min by lysing the phagosomal membrane . Nevertheless, the exact mechanisms by which O. tsutsugamushi enter the host cells and escape from the host endocytic pathway have not yet been elucidated. Human endothelial cells are believed to be the target cells mainly because the results of immunohistochemistry using autopsy tissues of scrub typhus patients . In contrast to autopsy results, dendritic cells and monocytes rather than endothelial cells are the target cells in eschars of scrub typhus patients , suggesting that the cells first encounter (i.e., dendritic cells and monocytes) O. tsutsugamushi can be different from those cells where O. tsutsugamushi are subsequently disseminated (endothelium) and accumulated. The involvement of local host immune responses and the mechanism for following systemic dissemination of Orientia may be very important to elucidate infection mechanism. A lot of efforts have focused on the immune responses upon infection of O. tsutsugamushi in various animal models to better understand the pathological observations in patients. Traditionally, the intraperitoneal (IP) inoculation is the route of choice to infect mice and has been used to evaluate the efficacy of several vaccine candidates [25,26,27,28]. Other routes for inoculation, such as intradermal (ID) [29, 30], intravenous (IV) [31,32,33,34], and more recently footpad , have been used in mouse models to study immunological responses. From all these studies, it appears that parameters like variations in route of inoculation, the mouse strains used, and strains of O. tsutsugamushi inoculum, can often lead to different challenge outcomes (e.g., lethal vs. non-lethal challenge models). However, it is difficult to conclude whether certain route or strain of mouse or strain of O. tsutsugamushi inoculum is really superior to the other as these studies had different emphases. Some models observed whether mice succumb to O. tsutsugamushi infection to evaluate the efficacy of vaccine candidates (i.e., IP), others suggested that a good mouse model is the one that shows similar pathology and target cells as observed in human scrub typhus (i.e., ID), and still others suggested that a good model should use an inoculation route that is similar to the natural infection route (i.e., ID or footpad). In addition to mouse models, a non-human primate monkey model has recently been used to evaluate the leading vaccine candidates [36, 37]. The study of the immune responses in non-human primate model is not as comprehensive as that in mouse model.
Microarray technology has been used extensively to investigate the host transcriptomic responses in various infectious diseases. Results from experiments conducted ex vivo or in mouse model have revealed insights into mechanisms of pathogenesis [38,39,40,41]. Similarly, microarray technology has been used to study host responses upon infections in animal models [42,43,44,45].
In this study, we inoculated C3H mice with O. tsutsugamushi via the IP route and characterized the changes in gene expression profiles at 10 different time points post inoculation. We observed two valleys (4–18 h and 2–4 days) of low number of differential expressed genes (DEG) with three peaks of high number of DEG at 2 h, 1 day and 7 days post inoculation. A spike pathway network analysis was performed and showed, in addition to a more general change of expression in pathways like the complement and coagulation cascades throughout most of the time course, there were also pathways that appear to be turning on or off at certain times post-inoculation. Further RT-qPCR was performed to confirm the expression of genes involved in complement and coagulation cascade. The findings highlight the complex connection of all different biological pathways. It is conceivable that during the initial phase of infection, Orientia is actively affecting pathways such as cell growth control, extracellular matrix adhesion and cell development in the host in order to survive. At later time post-infection, genes involved in cellular immune responses are triggered to defend infection and try to eliminate the pathogen.
Five mice per time point were observed for apparent clinical symptoms. All mice appeared healthy during the first 9 days with no signs of disease or illness. There was no eschar formation or any observable reaction at the injection site in any of the mice by 9 days post inoculation (dpi). Ruffled fur and loss of mobility were observed starting 9 dpi. All the 5 mice monitored for survival succumbed to the infection by 12 dpi, with an average survival of 11.2 dpi.
Antibody response following inoculation
Serum samples from each mouse were used to monitor the antibody responses elicited by O. tsutsugamushi infection. No Orientia-specific antibodies were detected in infected mice before or at 7 dpi. At 10 dpi, two mice were found non-responsive to outside stimuli and were euthanized to collect blood and tissues. The sera from these two mice had detectable IgG and IgM antibodies against the Orientia 56 kDa immunodominant protein, but the antibodies against the 47 kDa protein antigen were still undetectable. Antibodies were detected by 10 dpi, indicating that antibodies rose to detectable levels between 7 dpi and 10 dpi. This observation coincided with the appearance of ruffled fur and loss of mobility.
Detection of circulating soluble cell adhesion molecules
In order to investigate surrogate markers of endothelial and leukocyte activation following IP inoculation of O. tsutsugamushi, serum levels of circulating L-selectin, ICAM-1, and VCAM-1 were measured at selected time points. As shown in Fig. 1a, the amount of ICAM decreased as early as 8-h post inoculation (hpi) and continued to be lower than sham samples until 7 dpi when infected samples had more ICAM. Similar patterns were found for L-Selectin (Fig. 1b). The trend for VCAM (Fig. 1c) was different. The infected samples started to show increased amounts as early as 1 dpi and were higher than sham samples from 1 dpi to 7 dpi. Statistically significant differences in the amounts of these cell adhesion molecules between sham and infected samples were observed with ICAM at 8 hpi and 1 dpi. (Fig. 2a), for L-Selectin at 4 dpi and 7 dpi. (Fig. 2b), and for VCAM at 7 dpi. (Fig. 2c).
Tissue tropism of O. tsutsugamushi in IP inoculated mice
To determine the distribution of O. tsutsugamushi in mouse tissues after IP inoculation, lung, liver, and spleen were collected for bacterial quantitation by the Orientia-specific 47 kDa qPCR assay. Negligible amounts of O. tsutsugamushi were detected in all three tissues at 4 dpi (Table 1). Appreciable amounts of Orientia DNA were detected at 7 dpi and 10 dpi. Among these tissues, the liver had the highest amount of O. tsutsugamushi at all time points (4 dpi, 7 dpi, and 10 dpi), followed by the spleen, and then the lung. The fold increase of Orientia between 4 dpi and 7 dpi was similar to that between 7 dpi and 10 dpi for both liver and lung, leading to very high Orientia load in these tissues. The fold increase was much more moderate for the spleen.
Correlation of fold change
We applied Pearson’s test to correlate the fold change of expression at each time point. The higher the correlation indicates the higher similarity in gene responses between the two time points (Fig. 2). Noticeably, the gene response at 2 hpi was negatively correlated to 8 hpi (− 0.307) and 4 dpi (− 0.410), suggesting the expression between 2 hpi was opposite to that of 8 hpi and 4 dpi. In contrast, 12 hpi vs 2 dpi (0.349) and 8 hpi vs 4 dpi (0.317) have positive correlations, indicating that the gene responses between these time points are sharing similar patterns.
Individual differentially expressed probes analysis
A rank-sum test was used to identify the Differentially Expressed Probes (DEPs) (defined by p-value< 0.05 and |log(fold change)| > 0.5) for individual time points. We also applied moderated t-test (limma, R package). The number of DEPs from t-test is generally 1.5 times more than that of the rank-sum test, but the trend is the same. The stacked bar chart of Fig. 3 highlights the number of up-regulated (red) and down-regulated (green) DEPs. Significantly larger numbers of DEPs were observed at 2 hpi, 1 dpi, and 7 dpi than those of the adjacent times. After an initial response to IP inoculation at 2 hpi, it seemed that at 6 hpi and 8 hpi. Fewer probes showed significant changes, while from 8 hpi there is a sharp jump in DEPs (~ 2000 DEPs). However, it suddenly dropped to ~ 400 DEPs at 2 dpi and then gradually increased to more than 800 DEPs. Therefore, the three peaks at 2 hpi, 1 dpi, and 7 dpi may suggest a three-phase response, with an initial response to infection (i.e. 2 hpi) then 2 cycles (i.e., early (4 hpi – 1 dpi), and late stage (2 dpi – 7 dpi)) of increasing number of significantly expressed genes. The initial phase at 2 hpi might be different from the later two phases as they took longer to develop.
Gene cluster analysis
We further dissected the DEGs into spike patterns based on number of significantly up- or down-regulated probes at different time points. Figure 4 panel A shows top 20 spike patterns with the number of up- and down-regulated probes in descending order. These results revealed that three different spike patterns showed the highest number of probes with significant changes occurred at 2 hpi (pattern #8, #9 and #20), 18 hpi (pattern #5, #6 and #11), and 7 dpi (pattern #3, #4 and #18). On 1 dpi, 5 different patterns showed most up- and down-regulated genes. The majority of these 20 spike patterns had only single time point with most significantly up- and down-regulated probes. The exceptions were pattern #11, #18 and #20. Figure 4 panel B shows the number of probes /genes that are significantly up- or down-regulated at each time point. We further studied the enriched pathways of the top 20 clusters. As shown in Fig. 5, Rac1 pathway (up), TNFα signaling via NFκB (down), TGFβ signaling (down), interleukins signaling (up), and cytokine signaling (down) pathways are in general the major enriched pathways from 2 hpi to 1 dpi. Various signaling pathways were commonly expressed from 18 hpi, 1 dpi, and 7 dpi. Many metabolic pathways emerged at 2 dpi but disappeared later. Inflammation signaling pathways were activated between 18 hpi and 1 dpi. Various metabolic pathways were up regulated at 2 dpi and 4 dpi. At 7 dpi many pathways encompassing various biological functions such as WNT, TNF, mTOR, interferon, cytokines/chemokines, and toll-like receptor pathways were shown to be significantly up-regulated. The detailed enriched functions are listed in Additional file 1: Table S1.
Global pathway analysis
Pathway enrichment analysis identified the top 30 pathways with significant changes throughout the entire time course (Fig. 6). In general, around 40% of all genes involved in these pathways were significantly regulated. Several clusters pointed to pathways related to platelet degranulation and adhesion, coagulation, immune system (interferon alpha and immune-regulatory interaction), complement pathway, and iron metabolism.
Confirmation of selected pathways using qRT-PCR
The differential expression profiles of genes involved in blood coagulation cascades and classic complement pathways were confirmed using qRT-PCR. The fold-change of gene expression resulting after Orientia infection was confirmed. Four genes, including genes involved in complement pathway (C2 and cfb) and coagulation pathways (F5 and SerpinC1), were shown in Fig. 7 as examples to demonstrate that similar trends were observed by microarray probes (dotted lines) and qRT-PCR (solid lines).
Significant variations of gene expression in many key biological pathways are observed in C3H mice infected by Orientia. The C3H mouse model is a lethal model and all mice that received the 250 mLD50 of Orientia were sick and succumbed to the infection within 21 dpi. The expression profiles, while varying significantly for specific pathways, can be grouped into two valleys (4–18 hpi and 2–4 dpi) of differentially expressed genes and three peaks of differentially expressed genes at 2 hpi, 1 dpi, and 7 dpi. There were no genes associated with certain pathways showing a stable and continuous up- or down-regulation throughout the entire time course.
The IP challenge route resulted in dissemination of Orientia to different tissues in a time-dependent fashion with the maximum copy number of Orientia detected in the liver. Other challenge routes have shown different tissue tropisms, suggesting that different challenge routes may lead to different immune responses [29,30,31,32,33,34,35]. Orientia was detected in the liver as early as 7dpi, which is slightly earlier than the visible clinical signs, i.e., ruffled fur. It is intriguing that the lung, which is one of the target tissues associated with severe ST, appeared to have accumulated more Orientia in comparison to the spleen. We have also observed that higher numbers of Orientia accumulated in the lung were associated with delayed recovery in mice partially protected by immunization using recombinant protein antigens (data not shown). Leukocytes and endothelial cells were activated based on the level of surrogate markers. All three of these markers showed highest levels in infected mice at 7 dpi with a general upward trend suggesting that the activation occurs at around the same time when signs of sickness appeared and dissemination of Orientia to different tissues were detected. Both IL-1 and TNFα pathways were significantly up-regulated at 7 dpi and this may contribute to the upward trend of the three markers as the infection progressed.
There are pathways in which 40% of all genes involved were significantly regulated throughout the entire time course. Among these pathways, the complement, platelet and coagulation pathways are particularly interesting as they are among the top 30 pathways. These include intrinsic and extrinsic coagulation pathways, complement pathway, and platelet degranulation and adhesion pathways. While these 7 pathways have their respective unique pathophysiological roles, they do share some commonalities, particularly with respect to their innate defences against external threats (i.e., Orientia invasion) . There are several interesting observations to support our results. First, it was observed in an ectoparasite infection fish model that a non-lethal dosage of infection resulted in significant down-regulation of many genes in these pathways  similar to what we observed in the early phase of infection. On the contrary, when a lethal dose of infection was administered in the fish model, many of these same genes were up-regulated, just as what was observed in the late phase of Orientia infection. Second, our results are similar to a genome-wide transcriptome profiling of a murine acute melioidosis model in which early activation of complement system was observed at 24 h post-infection to maintain host cellular homeostasis and eliminate intracellular bacteria . Furthermore, a recent prospective study from Laos has shown that markers of coagulation activation and all inflammatory cytokines are significantly elevated, and in vivo coagulation activation ears prominent and seems to be related to a strong proinflammatory response  in scrub typhus pappatients. Finally, the activation of the complement system can also directly induce platelet activation and aggregation , and trigger the induction of tissue factor expression and activity in endothelial cells and neutrophils. The activation of endothelial cells may explain the procoagulant properties in scrub typhus patients with acute respiratory distress syndrome [51,52,53].
Among the 20 clusters of gene expression patterns, unique and shared pathway enrichments were observed at different time points. Especially noteworthy pathways in the early phase (2 hpi) of infection included Rac1, TNFα signaling via NFκB and G_α signaling events. These genes are mainly responsible for cell growth control, cytokine secretion, cell motility, extracellular matrix adhesion, cell transformation and invasion, and cell development. In the case of TNFα signalling via NFκB, these genes signify that NFκB-dependent signalling pathway could be crucial for the initial cellular responses to Orientia via the IP challenge route. As early as 2 hpi, genes responsible for the reorganization of cell-cell interaction and cell surface were differentially expressed, possibly due to the invasion of Orientia. As the infection continues, genes involved in additional pathways were differentially expressed. For example, the TGFβ signalling pathway and extracellular matrix–receptor interaction pathway were activated between 4 hpi and 8 hpi, cytokine signalling, FAS signalling, different interleukin signalling, Fcγ receptor mediated phagocytosis and Fcγ receptor signalling were all significantly regulated during 12 hpi to 18 hpi. These are consistent with the idea that early gene regulation involved Orientia invasion and host cellular immune responses.
During the late phase, from 1 dpi and 4 dpi, pathways related to many cellular functions were affected. In addition to pathways involved in cellular immune responses and MAPK/JAK_STAT pathways, several metabolic pathways were significantly affected, such as pyruvate metabolism and TCA cycle and respiratory electron transfer, bile acid metabolism, metabolism of lipids and lipoproteins, fatty acid β-oxidation, metabolism of amino acids, and derivatives, arachidonic acid metabolism, and tryptophan metabolism. These observations suggest that as the infection progressed, host cellular responses shifted from establishment of infection to affecting host metabolic pathways to potentially accommodate the need of Orientia growth. Specifically, we observed that Ido2 was one of the up-regulated genes at 2 dpi and its gene product is known to catabolize tryptophan to kynurenine. This is of interest because indoleamine 2,3-dioxygenase (Ido1), another gene product involved in catabolizing of tryptophan, is up-regulated in scrub typhus patients . Although Ido1 was not observed in this study, it is possible that a different route of infection for the two hosts (i.e., human and mouse) may explain the difference in observation of either Ido1 or Ido2 up-regulation in Orientia infected hosts. Many genes participating in fatty acid β-oxidation were also up-regulated, including Hmgcs2, that is involved in the ketogenesis pathway leading up to the generation of 3-hydroxybutyrate. Interestingly enough, the accumulation of 3-hydroxybutyrate has been observed in two metabolic studies in Orientia infected mouse models [55, 56]. These results demonstrated a significant perturbation of energy production.
Antigen presenting and processing, and Class I MHC mediated antigen presenting and processingwere significantly regulated at 7 dpi. These pathways are particularly interesting since antibody responses were also detectable after 7 dpi. In addition, genes such as Il15ra (1.5 folds), Il12rb1 (4.1 folds), Il18 (2.7 folds) were significantly up-regulated only at 7 dpi. Since these cytokines are all IFNγ-induced cytokines, it is not surprising that these cytokines are released into serum samples in high concentrations in scrub typhus patients . Furthermore, the transcriptomic study of monocytes infected by Orientia had shown many genes up-regulated at 8 hpi and 1 dpi, including, Il-1b, Cxcl10, Oas1, and Mx1 . These genes were also observed significantly up-regulated by 3.6, 16.0, 7.5 and 16.0 folds at 7 dpi, respectively. It is not clear whether any of these genes were specifically responsible for the overall outcome. Nevertheless, the combination of these genes, when they are all significantly up- or down-regulated simultaneously, is likely to be the major contributing factor for the final outcome of Orientia infection.
It is noted that IP is not a natural route for infection nor is the amount of live Orientia used be comparable with the chigger-fed mouse model . Regardless, there are similarities between the two models that could still provide information regarding cellular responses to Orientia infection. For example, in both cases, the time for onset of symptoms is around 9 dpi, and both routes resulted in animal death within 21 dpi. On the contrary, multiple variables are different between the two models. For example, the infection by chiggers may also involve an exchange of proteins between chigger and animal host, while IP injection does not. Consequently, it is important to explore the regulatory pathways in a chigger challenge model in order to better understand the cellular regulation upon Orientia infection.
Many pathways are dynamically regulated at different time points post Orientia IP infection in a mouse model. During the initial phase of infection, genes involved in pathways related to cell growth control, cytokine secretion, cell motility, extracellular matrix adhesion, cell transformation and invasion, and cell development are significantly regulated, supporting the idea that the Orientia is actively affecting the host in order to survive. At later time post-inoculation, genes involved in cellular immune responses are triggered to defend infection and try to eliminate the pathogen. We observed good correlation between regulation of various metabolic pathways and the change of amounts of various transcripts and proteins involved in them, suggesting that the growth of Orientia after initial invasion affected many different metabolic activities. Furthermore, the regulation of complement and coagulation pathways may play a role in the overall homeostasis of pathogens and host immune responses.
The procedures for animal study were approved by the Walter Reed Army Institute of Research (WRAIR) Naval Medical Research Center (NMRC) Institutional Animal Care and Use Committee (IACUC). The facility’s animal care and use program is accredited by Association for Assessment and Accreditation of Laboratory Animal Care International (AAALAC). Female eight-week old C3HeB/FeJ (C3H) inbred mice were purchased from Jackson Laboratories (Bar Harbor, Maine). Mice were divided into two groups and housed in an ABSL3 laboratory in cages (5 mice per cage) in accordance with the most current SOPs established by Veterinary Service Program (VSP) of NMRC. Feeding and water were provided ad libidum. One group of mice was inoculated with Snyder’s buffer (sham) to serve as naïve controls and those in the other group (infected) were infected via IP route  with 250 mLD50 of liver/spleen homogenate from mice infected by the Karp strain O. tsutsugamushi from New Guinea as described previously . We used Snyder’s buffer as a negative control instead of normal non-Orienita infected mouse liver-spleen homogenate based on results obtained from previously performed experiments (data not shown). The results from these experiments showed that sera from normal liver-spleen homogenate inoculated mice exhibited similar amounts of Th1 (TNF-α, IFN-γ, IL-2, and IL-12) and Th2 (IL-4, IL-5, and IL-10) cytokines as those from sera collected from sham inoculated mice. Therefore, the results suggest that only limited if any immune responses are activated by liver-spleen homogenate from uninfected mice. After inoculation, mice were observed twice daily for up to 12 days. Clinical signs, such as ruffled fur, mobility and morbidity were recorded. Each group of mice was divided into 5 mice per time point for a total of 10 different time points (2-h (hr), 4-h, 6-h, 8-h, 12-h, 18-h, 24-h, 2-day, 4-day and 7-day, Fig. 8). At the indicated times, mice from each group were euthanized by CO2, blood was drawn (detail described below) followed by collecting lung, liver, and spleen. Five additional mice were infected with the same infectious dose and monitored daily for up to 21 days or until they succumb to infection to ensure that 250 mLD50 is a lethal.
Blood collection, RNA extraction and data acquisition at indicated times
Before euthanasia whole blood from each mouse was collected by cardiac puncture. A total of 0.6 mL blood was immediately transferred to PAXgene™ blood RNA Tubes (Qiagen, Valencia, CA, USA). The samples were mixed and then flash frozen in liquid nitrogen until use. Additional blood was centrifuged to obtain serum for serology. RNA extraction was performed on a Qiacube™ (Qiagen, Valencia, CA, USA) and the quantity of extracted RNA was measured by NanoDrop™ 2000 spectrometer (Thermo Fisher Scientific, Waltham, MA, USA) as per manufacturer’s guideline. Agilent® SurePrint G3 Mouse GE 8x60K Microarray slides with associated Feature Extraction software v18.104.22.168 (Agilent Technologies, Santa Clara, CA, USA) were used following recommended protocol. Each SurePrint G3 Mouse GE 8x60K Microarray slide contains 27,122 RefSeq (Entrez) gene count and 4578 lncRNA transcripts. All genes have multiple probes.
Recombinant protein antigen based ELISA for antibody detection
Detection of soluble cell adhesion molecules (sCAMs)
Serum samples from all time points were assessed for sL-selectin, sICAM-1 and sVCAM-1 by using ELISA kits (R & D Systems, Minneapolis, MN, USA) following the manufacturer’s instructions. Samples were assayed in duplicate. The optical density at 465 nm (Vmax/Kinetic Microplate Reader, Molecular Devices, San Jose, CA, USA) was measured. The quantity of each molecule was calculated based on a standard curve established with known amount of the said molecule. The average quantity of each molecule in mouse sera from sham inoculated mice and from Orientia inoculated mice was compared by t-test at each time point.
Quantitation of O. tsutsugamushi
The quantity of O. tsutsugamushi in liver, lung, and spleen after inoculation was determined by a qPCR assay targeting the Orientia-specific 47 kDa gene . Mouse-specific primers for the single copy gene mouse complement factor D (cfd) were used to quantitate the amount of mouse gene equivalents in the extracted DNA as described by Sunyakumthorn et al. .
A real-time quantitative PCR (RT-qPCR) analysis was carried out using the RT2 Profiler PCR Array System (SA Biosciences Corp, Frederick, MD) focusing on the blood clotting cascade and classical complement pathway according to the manufacturer’s protocol. Five hundred ng of total RNA was used to produce cDNA using the RT2 First Strand Kit (Qiagen, Germantown, MD, USA), followed by qPCR assays using the RT2 SYBR Green/Rox Mastermix Kit (Qiagen, Germantown, MD, USA) in an Applied Biosystems 7900HT (Applied Biosystems, Foster City, CA, USA) under the recommended conditions. The RT-qPCR results were analyzed with SDS 2.3 software (Applied Biosystems, Foster City, CA, USA).
Data pre-processing to identify differentially expressed genes
For those samples that failed in RNA extraction or had a low RIN (RIN < 6.5), they were removed from microarray experiment. All data from microarray experiment was preprocessed and quality control was performed using within-chip Lowess and between-chip quantile normalization. Two outliers were identified by Principal Component Analysis and were removed from downstream analysis as they could bias the analysis. Batch correction of extraction time was done using Bioconductor limma 3.7 package  in R 3.2 . Differentially expressed genes (DEGs) were identified using Wilcoxon rank-sum test.
Clustering and enrichment analysis
The clustering and pathway analysis were performed using Functional Heatmap (https://bioinfo-abcc.ncifcrf.gov/Heatmap/index.php) and custom R code. Functional Heatmap first transfers gene expression levels into average fold change curves, and then transfers the changes above 0.5 as symbol “+”, below − 0.5 as symbol “-“, and the rest as “0”. Thus we generalized the curves into strings. Finally, the strings with similar appearance were grouped into patterns and the genes of selected patterns were sent to pathway enrichment analysis using the hypergeometric test. The pathway view showed the temporal pathway expression for known pathways collected from five pathway databases KEGG version 80 , Wiki pathway , Biocarta , Reactome , and GSEA . The dark green, green, red, and dark red stand for the percentage of the genes with log2 fold change of (−∞, − 1), (− 1,-0.5), (0.5, 1), and (1, +∞), respectively. Pathways with 95% common components were merged.
Availability of data and materials
The datasets used and analyzed during the current study are available from the corresponding author upon reasonable request. The results from analysis of the datasets are presented in this published article as tables, figures and supplementary information files.
Association for Assessment and Accreditation of Laboratory Animal Care International
Animal biosafety level 3
complement factor D
Differential expressed probes
days post inoculation
Gene set enrichment analysis
Institutional Animal Care and Use Committee
Kyoto Encyclopedia of Genes and Genomes
lethal dosage causing 50% death in a mouse
Naval Medical Research Center
hours post inoculation
soluble cell adhesion molecules
soluble intercellular adhesion molecule-1
Standard operation procedures
soluble vascular cell adhesion molecule-1
Veterinary services program
Brown GW, et al. Scrub typhus: a common cause of illness in indigenous populations. Trans R Soc Trop Med Hyg. 1976;70(5–6):444–8.
Izzard L, et al. Isolation of a novel Orientia species (O. chuto sp. nov.) from a patient infected in Dubai. J Clin Microbiol. 2010;48(12):4404–9.
Jiang J, et al. Laboratory diagnosis of two scrub typhus outbreaks at Camp Fuji, Japan in 2000 and 2001 by enzyme-linked immunosorbent assay, rapid flow assay, and Western blot assay using outer membrane 56-kD recombinant proteins. Am J Trop Med Hyg. 2003;69(1):60–6.
Kumar K, et al. Outbreak investigation of scrub typhus in Himachal Pradesh (India). J Commun Dis. 2004;36(4):277–83.
Chang WH. Current status of tsutsugamushi disease in Korea. J Korean Med Sci. 1995;10(4):227–38.
Sharma A, et al. Investigation of an outbreak of scrub typhus in the himalayan region of India. Jpn J Infect Dis. 2005;58(4):208–10.
Wang S, et al. Natural foci of tsutsugamushi disease in the Nan Peng lie islands in China. Chin Med J. 2002;115(2):272–5.
Lewis MD, et al. Scrub typhus reemergence in the Maldives. Emerg Infect Dis. 2003;9(12):1638–41.
Durand AM, et al. Scrub typhus in the Republic of Palau, Micronesia. Emerg Infect Dis. 2004;10(10):1838–40.
Tay ST, Kamalanathan M, Rohani MY. Antibody prevalence of Orientia tsutsugamushi, rickettsia typhi and TT118 spotted fever group rickettsiae among Malaysian blood donors and febrile patients in the urban areas. Southeast Asian J Trop Med Public Health. 2003;34(1):165–70.
Odorico DM, et al. New Orientia tsutsugamushi strain from scrub typhus in Australia. Emerg Infect Dis. 1998;4(4):641–4.
Lee YS, et al. Epidemiology of scrub typhus in eastern Taiwan, 2000-2004. Jpn J Infect Dis. 2006;59(4):235–8.
Weitzel T, et al. Endemic scrub typhus in South America. N Engl J Med. 2016;375(10):954–61.
Balcells ME, et al. Endemic scrub typhus-like illness, Chile. Emerg Infect Dis. 2011;17(9):1659–63.
Ghorbani R, et al. A case of scrub typhus probably acquired in Africa. Clin Infect Dis. 1997;25:1473–4.
Kawamura A, Tanaka H, Tamura A. Tsutsugamushi Disease. Tokyo: University of Tokyo Press; 1995.
Brown GW, et al. Single dose doxycycline therapy for scrub typhus. Trans R Soc Trop Med Hyg. 1978;72(4):412–6.
Watt G, et al. Scrub typhus infections poorly responsive to antibiotics in northern Thailand. Lancet. 1996;348(9020):86–9.
Mathai E, et al. Outbreak of scrub typhus in southern India during the cooler months. Ann N Y Acad Sci. 2003;990:359–64.
Murata M, et al. Proliferating sites of rickettsia tsutsugamushi in mice by different routes of inoculation evidenced with immunofluorescence. Jpn J Exp Med. 1985;55(5):193–9.
Paris DH, et al. Orientia tsutsugamushi in human scrub typhus eschars shows tropism for dendritic cells and monocytes rather than endothelium. PLoS Negl Trop Dis. 2012;6(1):e1466.
Rikihisa Y, Ito S. Entry of rickettsia tsutsugamushi into polymorphonuclear leukocytes. Infect Immun. 1982;38(1):343–50.
Urakami H, Tsuruhara T, Tamura A. Penetration of rickettsia tsutsugamushi into cultured mouse fibroblasts (L cells): an electron microscopic observation. Microbiol Immunol. 1983;27(3):251–63.
Moron CG, et al. Identification of the target cells of Orientia tsutsugamushi in human cases of scrub typhus. Mod Pathol. 2001;14(8):752–9.
Yu Y, et al. Induction of protective immunity against scrub typhus with a 56-kilodalton recombinant antigen fused with a 47-kilodalton antigen of Orientia tsutsugamushi Karp. Am J Trop Med Hyg. 2005;72(4):458–64.
Seong SY, et al. Neutralization epitopes on the antigenic domain II of the Orientia tsutsugamushi 56-kDa protein revealed by monoclonal antibodies. Vaccine. 2000;19(1):2–9.
Seong SY, et al. Induction of homologous immune response to rickettsia tsutsugamushi Boryong with a partial 56-kilodalton recombinant antigen fused with the maltose-binding protein MBP-Bor56. Infect Immun. 1997;65(4):1541–5.
Ni YS, et al. Protection against scrub typhus by a plasmid vaccine encoding the 56-KD outer membrane protein antigen gene. Am J Trop Med Hyg. 2005;73(5):936–41.
Sunyakumthorn P, et al. An intradermal inoculation model of scrub typhus in Swiss CD-1 mice demonstrates more rapid dissemination of virulent strains of Orientia tsutsugamushi. PLoS One. 2013;8(1):e54570.
Soong L, et al. An intradermal inoculation mouse model for immunological investigations of acute scrub typhus and persistent infection. PLoS Negl Trop Dis. 2016;10(8):e0004884.
Shelite TR, et al. A hematogenously disseminated Orientia tsutsugamsushi-infected murine model of scrub typhus. PLoS Negl Trop Dis. 2014;8(7):e2966.
Shelite TR, et al. IL-33-dependent endothelial activation contributes to apoptosis and renal injury in Orientia tsutsugamushi-infected mice. PLoS Negl Trop Dis. 2016;10(3):e0004467.
Xu G, et al. CD8+ T cells provide immune protection against murine disseminated endotheliotropic Orientia tsutsugamushi infection. PLoS Negl Trop Dis. 2017;11(7):e0005763.
Soong L, et al. Strong type 1, but impaired type 2, immune responses contribute to Orientia tsutsugamushi-induced pathology in mice. PLoS Negl Trop Dis. 2014;8(9):e3191.
Keller CA, et al. Dissemination of Orientia tsutsugamushi and inflammatory responses in a murine model of scrub typhus. PLoS Negl Trop Dis. 2014;8(8):e3064.
Chattopadhyay S, et al. Scrub typhus vaccine candidate Kp r56 induces humoral and cellular immune responses in cynomolgus monkeys. Infect Immun. 2005;73(8):5039–47.
Paris DH, et al. A nonhuman primate scrub typhus model: protective immune responses induced by pKarp47 DNA vaccination in cynomolgus macaques. J Immunol. 2015;194(4):1702–16.
Gupta P, et al. Differential host gene responses from infection with neurovirulent and partially-neurovirulent strains of Venezuelan equine encephalitis virus. BMC Infect Dis. 2017;17(1):309.
Priya GB, et al. Genome wide host gene expression analysis in mice experimentally infected with Pasteurella multocida. PLoS One. 2017;12(7):e0179420.
Sahlberg Bang C, et al. Global gene expression profiling and antibiotic susceptibility after repeated exposure to the carbon monoxide-releasing molecule-2 (CORM-2) in multidrug-resistant ESBL-producing uropathogenic Escherichia coli. PLoS One. 2017;12(6):e0178541.
VanLeuven JT, et al. Lung epithelial cells have virus-specific and shared gene expression responses to infection by diverse respiratory viruses. PLoS One. 2017;12(6):e0178408.
Hammamieh R, et al. Blood genomic profiles of exposures to Venezuelan equine encephalitis in Cynomolgus macaques (Macaca fascicularis). Virol J. 2007;4:82.
Giotis ES, et al. Transcriptomic profiling of virus-host cell interactions following chicken Anaemia virus (CAV) infection in an in vivo model. PLoS One. 2015;10(8):e0134866.
Baskin CR, et al. Integration of clinical data, pathology, and cDNA microarrays in influenza virus-infected pigtailed macaques (Macaca nemestrina). J Virol. 2004;78(19):10420–32.
Cheng Q, et al. Quantitative trait loci and candidate genes for neutrophil recruitment in sterile inflammation mapped in AXB-BXA recombinant inbred mice. PLoS One. 2015;10(5):e0124117.
Oikonomopoulou K, et al. Interactions between coagulation and complement--their role in inflammation. Semin Immunopathol. 2012;34(1):151–65.
Yin F, et al. Transcriptome and analysis on the complement and coagulation cascades pathway of large yellow croaker (Larimichthys crocea) to ciliate ectoparasite Cryptocaryon irritans infection. Fish Shellfish Immunol. 2016;50:127–41.
Chin CY, Monack DM, Nathan S. Genome wide transcriptome profiling of a murine acute melioidosis model reveals new insights into how Burkholderia pseudomallei overcomes host innate immunity. BMC Genomics. 2010;11:672.
Paris DH, et al. Coagulation and inflammation in scrub typhus and murine typhus--a prospective comparative study from Laos. Clin Microbiol Infect. 2012;18(12):1221–8.
Polley MJ, Nachman RL. Human platelet activation by C3a and C3a des-arg. J Exp Med. 1983;158(2):603–15.
Kurup A, et al. Scrub typhus with sepsis and acute respiratory distress syndrome. J Clin Microbiol. 2013;51(8):2787–90.
Angelakis E, et al. Orientia tsutsugamushi in lung of patient with acute respiratory distress syndrome, France, 2013. Emerg Infect Dis. 2015;21(2):373–5.
Kambas K, et al. C5a and TNF-alpha up-regulate the expression of tissue factor in intra-alveolar neutrophils of patients with the acute respiratory distress syndrome. J Immunol. 2008;180(11):7368–75.
Prachason T, et al. Activation of indoleamine 2,3-dioxygenase in patients with scrub typhus and its role in growth restriction of Orientia tsutsugamushi. PLoS Negl Trop Dis. 2012;6(7):e1731.
Jung J, et al. Metabolic responses to Orientia tsutsugamushi infection in a mouse model. PLoS Negl Trop Dis. 2015;9(1):e3427.
Chao CC, et al. Metabolic characterization of serum from mice challenged with Orientia tsutsugamushi –infected mites. New Microbes New Infect. 2018;23:70–6.
Chierakul W, et al. Differential expression of interferon-gamma and interferon-gamma-inducing cytokines in Thai patients with scrub typhus or leptospirosis. Clin Immunol. 2004;113(2):140–4.
Tantibhedhyangkul W, et al. Orientia tsutsugamushi stimulates an original gene expression program in monocytes: relationship with gene expression in patients with scrub typhus. PLoS Negl Trop Dis. 2011;5(5):e1028.
Ching WM, Lurchachaiwong W, Zhang Z, Awoyomi T, Chao CC, Schuster A. Evaluation of a Recombinant Vaccine Candidate r56Lc-1 in a Chigger Challenge Mouse Model. J Vaccines Vaccination. 2014;05(06). https://doi.org/10.4172/2157-7560.1000257.
Chan TC, et al. Development of a rapid method for determining the infectious dose (ID)50 of Orientia tsutsugamushi in a scrub typhus mouse model for the evaluation of vaccine candidates. Vaccine. 2003;21(31):4550–4.
Chen HW, et al. Kinetics and magnitude of antibody responses against the conserved 47-kilodalton antigen and the variable 56-kilodalton antigen in scrub typhus patients. Clin Vaccine Immunol. 2011;18(6):1021–7.
Frey A, Di Canzio J, Zurakowski D. A statistically defined endpoint titer determination method for immunoassays. J Immunol Methods. 1998;221(1–2):35–41.
Huber E, et al. Loop-Mediated Isothermal Amplification Assay Targeting the 47-Kda Gene of Orientia tsutsugamushi: A Rapid and Sensitive Alternative to Real-Time PCR. J Med Microbiol Diagn. 2012;01(04):2161–0703.
Ritchie ME, et al. Limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43(7):e47.
Team RC. R: a language and environment for statistical computing; 2013.
Kanehisa M, Goto S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27–30.
Slenter DN, et al. WikiPathways: a multifaceted pathway database bridging metabolomics to other omics research. Nucleic Acids Res. 2017;46(D1):D661–7.
Nishimura D, BioCarta. Biotech Software & Internet Report: The Computer Software. J Scient. 2001;2(3):117–20.
Joshi-Tope G, et al. Reactome: a knowledgebase of biological pathways. Nucleic Acids Res. 2005;33(suppl_1):D428–32.
Subramanian A, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci. 2005;102(43):15545–50.
The authors wish to thank the NMRC Vet Med staff for all of their help in conducting the animal study.
This work was supported by work unit number 6000.RAD1.J.A0311. The funder had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Ethics approval and consent to participate
The study protocol was reviewed and approved by the Walter Reed Army Institute of Research and Naval Medical Research Institutional Animal Care and Use Committee in compliance with all applicable Federal regulations governing the protection of animals in research.
Consent for publication
The authors declare that they have no competing interests. CCC, RH, MJ, and WMC are US Government employees. The work of these individuals was prepared as part of official government duties. Title 17 U.S.C. §105 provides that ‘Copyright protection under this title is not available for any work of the United States Government’. Title 17 U.S.C. §101 defines a U.S. Government work as a work prepared by a military service member or employee of the U.S. Government as part of that person’s official duties. The views expressed are those of the authors and do not necessarily reflect the official policy or position of the Department of the Navy, Department of the Army, Department of Defense, nor the U.S. Government.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Chao, C., Yang, R., Zhang, Z. et al. Temporal analysis of mRNA expression profiles in Orientia infected C3HeB/FeJ mouse. BMC Microbiol 20, 3 (2020). https://doi.org/10.1186/s12866-019-1684-3
- Scrub typhus
- cDNA microarray
- Orientia tsutsugamushi infection
- Semi-quantitative PCR
- Gene expression profiles