Skip to main content
  • Research article
  • Open access
  • Published:

Genotypic characterization of Orientia tsutsugamushi from patients in two geographical locations in Sri Lanka

Abstract

Background

To date more than 20 antigenically distinct strains of Orientia tsutsugamushi (OT) reported within the tsutsugamushi triangle that cause an undifferentiated acute febrile illness in humans. Genotypic characterization of OT in different geographic regions or within the same country, is important in order to establish effective diagnostics, clinical management and to develop effective vaccines. Genetic and antigenic characterization of OT causing human disease in OT-endemic regions is not known for Sri Lanka.

Methods

Adult patients and children who were admitted with an acute febrile illness and presumed to having acute scrub typhus based on presence of an eschar and other supporting clinical features were recruited. Eschar biopsies and buffy coat samples collected from patients who were confirmed having OT by IFA were further studied by real time PCR (Orientia 47 kD) and nested PCR (Orientia 56 kD) amplification. DNA sequences were obtained for 56 kD gene amplicons and phylogenetic comparisons were analyzed using currently available data in GenBank [Neucleotide substitution per 100 residues, 1000 Bootstrap Trials].

Results

Twenty eschar biopsies (Location1,19, Location 2,1) and eight buffy coat samples (Location1,6, Location2,2) examined by real time PCR revealed Orientia amplicons in 16 samples. DNA sequences were obtained for the 56 kD gene amplicons in 12 eschars and 4 buffy coat samples. The genotypes of the Location1 samples revealed that, 7 exhibiting close homology with JP1 [distantly related to UT177 Thai (Karp related)], five had close homology with Kato strain, two had close homology with JGv and JG AF [Distantly related to Kawasaki M63383] and one had close homology with Gilliam strain. The Location 2 strain was closely related to Kuroki-Boryong L04956, the genotype which is distributed in far eastern Asia. Similar to other patients in the cohort this patient also had never travelled out of Sri Lanka.

Conclusions

We observed all three main OT genotypes in Sri Lanka, and the majority fell into Thai Karp related clade. These results demonstrate great antigenic diversity of OT in the studied areas of Sri Lanka.

Peer Review reports

Background

Scrub typhus (ST), or tsutsugamushi disease, is an acute febrile illness in humans caused by infection with Orientia tsutsugamushi (OT) following a bite of an infected mite vector of the genus Leptotrombidium [1]. Scrub typhus is endemic in the Asia-Pacific region, extending from Afghanistan to China, Korea, the islands of the western Pacific and Indian Oceans, and northern Australia [2, 3]. This endemic region is often referred to as the tsutsugamushi triangle, and hosts approximately 1 billion people [4]. The vectors can be found in a variety of ecological conditions, from the mountainous regions of northern India to the tropical climates of the Malay Peninsula and Indonesia [2]. Transovarial transmission of OT within vectors appears to be essential to maintenance of the agent in nature; thus, the mite serves as both the vector and the reservoir [2]. Transmission of the etiologic agent to vertebrate hosts occurs during feeding of the larval or “chigger” stage of mites [2]. While Orientia is vertically maintained in Leptotrombidium mite populations, it may be transmitted horizontally from mites to vertebrate hosts [2]. The transmission to humans is incidental. Currently, there is no vaccine against ST [5].

The disease mimics several other tropical febrile illnesses, and can vary from mild to fatal disease, with reported mortality rates of 35%–50% during the pre-antibiotic era [6, 7] to 7–9% currently, in different geographical regions [5]. Although detection of eschars assists early diagnosis, the occurrence of this sign is highly variable and depends on the host, geographical region, and possibly the bacterial strain [8]. ST became more familiar during World War II, when soldiers deployed to endemic regions were affected in great numbers [2]. Research since World War II has highlighted dramatic antigenic variation among strains of OT [2] Although great inter-strain variability in virulence has been shown in mouse models, it is not clear whether virulence for mice can be directly applicable to humans [2]. However, the antigenic variation may result in a spectrum of illness from very mild to disease that is often fatal when untreated [2]. To date, 20 antigenically distinct strains of OT have been reported, including the initially serologically distinguished prototypic strains Karp, Gilliam, and Kato [9]. The introduction of chloramphenicol and, later, the tetracyclines dramatically eliminated mortality among cases where treatment was started early [1012]. Although relapses and reinfections occur [12], infection is generally responsive to treatment with antibiotics such as doxycycline and azithromycin, even when antigenically diverse strains are involved [2].

The documented history of scrub typhus in Sri Lanka dates back to the Second World War [5]. Similar to other countries in the tsutsugamushi triangle where scrub typhus is endemic, Sri Lanka also experiences a variable burden due to the illness. The incidence of the disease seems to be based on human activity, and ecological and climatic factors. Although OT in Sri Lanka is confirmed with IFA-based diagnostic techniques, the genotypic characterization of infecting OT strains is unknown. Such genotypic characterization of OT is important to establish effective diagnostics and to develop effective vaccines.

The current study was designed to identify the strains of OT causing clinical illness in two areas along the western coastal belt of Sri Lanka. Both areas (Fig. 1), which are approximately 100 km apart, are endemic for OT. They are likely to have similar ecological characteristics, both are in the wet zone and 5–10 km from the coast.

Fig. 1
figure 1

A map of Sri Lanka to highlight the geographical locations of the two study sites

Methods

Paediatric and adult patients who presented with an acute febrile illness to Balapitiya Base Hospital (Location 1) and the Professorial Medical Unit of Colombo North Teaching Hospital (Location 2) in Sri Lanka were examined clinically in order to identify an eschar; the essential inclusion criteria for the study. Consecutive patients who fulfilled the essential criteria (acute fever and an eschar) were recruited for the study. For patients who were more than 18 years of age a punch biopsy of the eschar was performed after obtaining informed written consent. For those who were less than 18 years old, as they were considered not competent for consenting according to the Sri Lankan law, a 3 cc sample of blood was collected after obtaining informed written consent from their guardians. Patients who were considered severely ill, declined to consent, or had received anti-rickettsial antibiotics for more than 24 h were excluded from the study. Relevant demographic, clinical, and laboratory data of all patients were collected into a predesigned data collectionform.

Biopsy tissue and buffy coat samples were stored in RNAlater (Life Technologies, Grand Island, NY) according to the manufacturer’s instructions. The samples were stored at−80°C until they were sent to the University of Texas Medical Branch, Galveston, USA, for further analysis. The patients were further tested by serology to detect IgG antibodies against Karp, Kato or Gilliam antigens by IFA in order to confirm OT infection. The IFA-IgG test was done using antigen coated slides procured from Fuller Laboratories (Fullerton, California, USA) adhering to manufacturer’s instructions.

Prior to processing, biopsy tissues stored in RNAlater were rinsed in phosphate buffered saline (PBS). A sterile scalpel blade was used to cut a one fourth sized wedge out of each tissue sample for DNA extraction. The biopsy wedges were placed in 2 mL microcentrifuge tubes with 200 μL of PBS and two 4 mm stainless steel grinding balls. Tissues were homogenized with a Retsch MM300 mixer mill (Bio-Rad, Hercules, CA) for 2 min at 30 Hz. DNA was extracted from the homogenate using the DNeasy Blood and Tissue Kit (QIAGEN, Valencia, CA).

Buffy coat samples stored in RNAlater were centrifuged in 2.0 mL microcentrifuge tubes at 5000 × g for 5 min. The supernatant was removed, and the remaining cell pellet was resuspended in 200 μL of PBS with 200 μL of lysis buffer (QIAGEN). DNA was extracted using the DNeasy Blood and Tissue Kit.

Real time PCR was performed on extracted DNA using primers designed to amplify a segment of the Orientia 47-kD antigen gene (Additional file 1: Table S1). [13] Of the samples positive for OT DNA, nested PCR to amplify a portion of the 56-kD protein gene was performed as previously described (Additional file 1: Table S1). [14] PCR amplicons were cloned using chemically competent E. coli cells and the TOPO TA Cloning Kit (Life Technologies). Plasmid DNA was purified with the PureLink Plasmid Miniprep Kit (Life Technologies) and sequenced using a 3130xl Genetic Analyzer (Life Technologies).

The partial 56-kD sequences were trimmed of primer sequences using EditSeq of the Lasergene Core Suite software version 12.0 (DNASTAR, Madison, WI) and then identified by comparing them with sequences available in GenBank via BLAST database search (http,//blast.ncbi.nlm.nih.gov). Phylogenetic analyses of the study and reference sequences were performed using MegAlign of the Lasergene Core Suite (DNASTAR). Sequence alignment was performed using CLUSTAL W, and a phylogenetic tree was constructed with 1,000 bootstrap replications using the neighbor-joining method. For residue conservation analysis, nucleotide sequences were translated from open reading frames of DNA sequences using EditSeq. These nucleotide sequences were then aligned and analyzed using MegAlign.

Results

During the period of six months from January - August 2013, 19 eschar biopsies and 6 buffy coat samples were collected from study site 1 (Balapitiya), and 1 eschar biopsy and 2 buffy coat samples were collected from study site 2 (Ragama). Real time PCR for the Orientia 47-kD protein gene showed amplicons in 16 samples (15/25 at study site 1 and 1/3 at study site 2). Of those positive by real time PCR, nested PCR for the 56-kD gene amplified DNA from 12 eschars and 4 buffy coat samples (Additional file 2: Figure S1).

The genotypes of the Location 1 samples showed that seven exhibited close homology with JP1 [distantly related to UT177 Thai (Karp-related)], five had close homology with Kato strain, two had close homology with JGv and JG AF [distantly related to Kawasaki M63383], one had close homology with Gilliam strain. The Location 2 strain was closely related to Kuroki-Boryong L04956 (Fig. 2).

Fig. 2
figure 2

Phylogenetic tree to demonstrate the genetic relationship of O. tsutsugamushi isolates from the two study sites

Alignment and residue conservation analysis of the nucleotide sequences demonstrated that genotypic differences between samples conferred nucleotide heterogeneity, which could result in antigenic variation in the studied region (Additional file 3: Figure S2 and Additional file 4: Figure S3).

Discussion

Originally, three distinctive antigenic prototypes of OT (Karp, Kato, and Gilliam) were described using serologic testing (via complement fixation) [15]. The serologic diversity is based largely on the 56-kDa type specific antigen (TSA), [16]. Later, over 20 additional strains of OT were found [9] by indirect immunofluorescence assays (IFA) using strain- or type-specific monoclonal antibodies or hyperimmune sera which recognize the 56-kDa TSA and by the genetic sequences of this genetically diverse protein. Recent efforts to catalog the strains present in various locations have identified novel subtypes. Multi-locus sequence typing (MLST) analysis of Thai isolates has shown a very high rate of recombination, which along with gene duplication and horizontal gene transfer accounts for the genetic diversity among OT isolates in various locations of Thailand [17, 18].

An increased prevalence of Karp and Karp-like strains has been reported in the Southeast Asian countries, including Malaysia (56%), Thailand (65%), Myanmar (46%), and Philippines (78%) [1921], Boryong is the predominant strain reported from Korea [22]. However, several studies in India have reported a high prevalence of Kato-like organisms (65%) and the majority have sequence similarity with OT isolates from Korea, Thaiwan, Vietnam and Cambodia [23]. Around 35% were Karp-like, and the rest were related to Kawasaki and Gilliam types.

This is the first study from Sri Lanka to describe the genotypes of OT, in two geographical areas that are about 100 km apart along the western coastal belt within the wet zone of Sri Lanka with probably similar ecological conditions. We were able to recruit 25 patients from study site 1 and 3 patients in study site 2 during the 6 months of the study. However, there were only 15/25 at study site 1 and 1/3 at study site 2 were positive by real time PCR, nested PCR for the 56-kD gene. The patients who were negative by PCR fulfilled criteria for recruitment to the study and had very high IFA-IgG titres against Orientia antigens. However, it is very likely that they had either received anti-rickettsial antibiotics prior to admission to hospital by the general practitioners or were recruited at a very late stage of the infection and therefore had already cleared the agent from their blood or tissues.

Although the study included only 16 positive samples, it is interesting to note that we have observed great genetic diversity of OT within the relatively small geographical areas that were studied. In contrast to what has been reported from India, but similar to reports from southeastern parts of Asia, we noted a slightly higher prevalence of Karp [Thai-related] genotypes followed by Kato-related, Kawasaki-related and Gilliam-related genotypes. However, it is interesting to note the close homology with the Kuroki-Boryong organism in one patient from study site 2 (Ragama) who, like the others, had never travelled out of Sri Lanka. The main limitation of this study is that it was carried out over a short period of 6 months and included only 16 positive samples (15 from study site 1 and only 1 from study site 2). Although we observed a close temporal distribution of genotypes within the 6 months, a longer study recruiting more epidemiologically representative samples, is likely to give more insight into the pattern of circulation of the different strains of O. tsutsugamushi in the two selected geographical regions. Design of such a study requires an epidemiological survey of the two study sites in order to understand the exact prevalence of this re-emerging infection in the two areas – such data is currently lacking.

A 2011 study in Korea comparing the clinical severity between Boryong and Karp genotypes suggested that eschars and rashes were found in 97% and 94% of the patients infected with the Boryong cluster compared to 73.7% and 68.4% of the patients infected with the Karp cluster, respectively, suggesting clinical variation among infections with different serotypes of OT [24]. In a previous study conducted in the same geographical region, we found a high eschar detection rate (89%) that is in keeping with the above findings [25].

Since this study analyzed only 16 positive samples, it is difficult to arrive at conclusions on the complete composition of OT genotypes in Sri Lanka. However, since Sri Lanka has diverse ecological and environmental conditions within the matter of 2–3 h travel, a larger scale study that covers all ecological zones could contribute to understanding of the complex diversity of OT infection and its transmission within the country.

Conclusions

We observed five genotypes of O. tsutsugamushi demonstrating great genetic diversity in relatively small geographic areas of Sri Lanka, with the majority falling into the Karp-related clade. We believe that such knowledge will contribute to potential prediction of clinical illness, production of future diagnostic tests, and development of vaccines.

Abbreviations

CA:

Canada

DNA:

Deoxy ribonucleic acis

IFA:

Indirect immunofluorescence antibody

IgG:

Immunoglobin G

MLST:

Multilocus sequence typing

NY:

New York

OT:

Orientia tsutsugamushi

PBA:

Phosphate buffered salina

PCR:

Polymerase chain reaction

RNA:

Ribonucleic acid

ST:

Scrub typhus

TSA:

Type specific antigen

USA:

United States of America

References

  1. Traub R, Wisseman Jr CL. The ecology of chigger-borne rickettsiosis (scrub typhus). J Med Entomol. 1974;11:237–303.

    Article  CAS  PubMed  Google Scholar 

  2. Kelly DJ, Fuerst PA, Ching WM, Richards AL. Scrub typhus, the geographic distribution of phenotypic and genotypic variants of Orientia tsutsugamushi. Clin Infect Dis. 2009;15(48 Suppl 3):S203–30. doi:10.1086/596576.

    Article  Google Scholar 

  3. Oaks Jr SC, Ridgway RL, Shirai A, Twartz JC. Scrub typhus. Inst Med Res Malays Bull. 1983;21:1–98.

    Google Scholar 

  4. Chattopadhyay S, Richards AL. Scrub typhus vaccines, past history and recent developments. Hum Vaccines. 2007;3:73–80.

    Article  CAS  Google Scholar 

  5. WHO. Frequently asked questions; Scrub typhus. http,//www.searo.who.int/entity/emerging_diseases/CDS_faq_Scrub_Typhus.pdf. Accessed on May 1, 2016.

  6. Smadel JE, Elisberg BL. Scrub typhus rickettsia. In: Horsfall Jr FL, Tamm I, editors. Viral and rickettsial infections of man. 4th ed. Philadelphia: Lippincott; 1965. p. 1130–42.

    Google Scholar 

  7. Kawamura A, Tanaka H. Rickettsiosis in Japan. Jpn J Exp Med. 1988;58:169–84.

    PubMed  Google Scholar 

  8. Kim DM, Won KJ, Park CY, Yu KD, Kim HS, Yang TY, et al. Distribution of eschars on the body of scrub typhus patients, a prospective study. Am J Trop Med Hyg. 2007;76:806–9.

    PubMed  Google Scholar 

  9. Kelly DJ, Richards AL, Temenak J, Strickman D, Dasch GA. The past and present threat of rickettsial diseases to military medicine and international public health. Clin Infect Dis. 2002;34 Suppl 4:S145–69.

    Article  PubMed  Google Scholar 

  10. Smadel JE, Woodward TE, Ley Jr HL, Green R, Mankikar DS. Chloromycetin in the treatment of scrub typhus. Science. 1948;108:160–1.

    Article  CAS  PubMed  Google Scholar 

  11. Smadel JE, Traub R, Ley Jr HL, Philip CB, Woodward TE, Lewthwaite R. Chloramphenicol (Chloromycetin) in the chemoprophylaxis of scrub typhus (tsutsugamushi disease). II. Results with volunteers exposed in hyperendemic areas of scrub typhus. Am J Hyg. 1949;50:75–91.

    CAS  PubMed  Google Scholar 

  12. Smadel JE, Woodward TE, Ley Jr HL, Lewthwaite R. Chloramphenicol (Chloromycetin) in the treatment of tsutsugamushi disease (scrub typhus). J Clin Invest. 1949;28:1196–215.

    Article  PubMed Central  Google Scholar 

  13. Jiang J, Chan TC, Temenak JJ, Dasch GA, Ching WM, Richards AL. Development of a quantitative real-time polymerase chain reaction assay specific for Orientia tsutsugamushi. Am J Trop Med Hyg. 2004;40:351–6.

    CAS  Google Scholar 

  14. Kim DM, Kim HL, Park CY, Yang TY, Lee JH, Yang JT, et al. Clinical usefulness of eschar polymerase chain reaction for the diagnosis of scrub typhus, A prospective study. Clin Infect Dis. 2006;43:1296–300.

    Article  PubMed  Google Scholar 

  15. Shishido A. Identification and serological classification of the causative agent of scrub typhus in Japan. Jpn J Med Sci Biol. 1962;15:308–21.

    Google Scholar 

  16. Shishido A. Strain variation of Rickettsia orientalis in the complement fixation test. Jpn J Med Sci Biol. 1964;17:59–72.

    Article  CAS  PubMed  Google Scholar 

  17. Nakayama K, Yamashita A, Kurokawa K, Morimoto T, Ogawa M, Fukuhara M. The whole-genome sequencing of the obligate intracellular bacterium Orientia tsutsugamushi revealed massive gene amplification during reductive genome evolution. DNA Res. 2008;15:185–99.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Sonthayanon P, Peacock SJ, Chierakul W, Wuthiekanun V, Blacksell SD, Holden MT. High rates of homologous recombination in the mite endosymbiont and opportunistic human pathogen Orientia tsutsugamushi. PLoS Negl Trop Dis. 2010;4:752.

    Article  Google Scholar 

  19. Shirai A, Robinson DM, Brown GW, Gan E, Huxsoll DL. Antigenic analysis by direct immunofluorescence of 114 isolates of Rickettsia tsutsugamushi recovered from febrile patients in rural Malaysia. Jpn J Med Sci Biol. 1979;32:337–44.

    Article  CAS  PubMed  Google Scholar 

  20. Blacksell SD, Luksameetanasan R, Kalambaheti T, Aukkanit N, Paris DH, McGready R. Genetic typing of the 56-kDa type-specific antigen gene of contemporary Orientia tsutsugamushi isolates causing human scrub typhus at two sites in north-eastern and western Thailand. FEMS Immunol Med Microbiol. 2008;52:335–42.

    Article  CAS  PubMed  Google Scholar 

  21. Bengtson IA. A serological study of 37 cases of tsutsugamushi disease (scrub typhus) occurring in Burma and the Philippine Islands. Public Health Rep. 1946;61:887–94.

    Article  Google Scholar 

  22. Jeong HW, Choi YK, Baek YH, Seong MH. Phylogenetic analysis of the 56-kDa type-specific protein genes of Orientia tsutsugamushi in central Korea. J Korean Med Sci. 2012;27:1315–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Varghese GM, Janardhanan J, Trowbridge P, Peter JV, Prakash JA, Sathyendra S, et al. Scrub typhus in South India, clinical and laboratory manifestations, genetic variability, and outcome. Int J Infect Dis. 2013;7:e981–7. doi:10.1016/j.ijid.2013.05.017. Epub 2013 Jul 26.

    Article  Google Scholar 

  24. Kim DM, Yun NR, Neupane GP, Shin SH, Ryu SY, Yoon HJ. Differences in clinical features according to Boryoung and Karp genotypes of Orientia tsutsugamushi. PLoS One. 2011;6:e22731. doi:10.1371/journal.pone.0022731. Epub 2011 Aug 15.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Premaratna R, Loftis AD, Chandrasena TG, Dasch GA, de Silva HJ. Rickettsial infections and their clinical presentations in the Western Province of Sri Lanka, a hospital-based study. Int J Infect Dis. 2008;12:198–202. Epub 2007 Sep 27.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We thank Dr. Kumara TADN, Dr. Liyanarachchi MS, Balapitiya Base hospital, Sri Lanka

Mr. Bandara, NKBKRGW, Department of Microbiology, Faculty of Medicine, University of Kelaniya, Sri Lanka

Dr. Weerasinghe M., Professorial Medical Unit and

Dr. Gunasekara C, Dr. Mendis S, Department of Dermatology and the staff Colombo North Teaching Hospital, Ragama, Sri Lanka for the support given for this work.

Availability of material and data

The first and the second authors possess relevant clinical and experimental data.

Authors’ contributions

RP, LB, SDN, NdeS, TGANC, DHW, HJdeS, Designed the study. RP, SDN, NDeS, Collected relevant clinical materials. TGANC, Carried out relevant tests locally. LB, DHW, carried out advanced laboratory testing. RP, LB, DHW, TGANC and HJdeS, Wrote the manuscript. All authors read and approved the final version of the manuscript.

Competing interests

Authors declare no competing interests.

Consent for publication

All participants above 18 years of age and guardians of those who are less than 18 years consented for publication of relevant data without revealing their identity.

Ethics approval and consent to participate

The Ethics Review Committee of the Faculty of Medicine, University of Kelaniya, Sri Lanka granted ethical approval for the study.

All participants above 18 years of age and guardians of those who are less than 18 years consented for the participation in the study.

Financial Support

This work was supported by the Carmage and Martha Walls Distinguished Chair in Tropical Diseases Endowment.

Author information

Authors and Affiliations

Authors

Corresponding authors

Correspondence to Ranjan Premaratna or Lucas S. Blanton.

Additional files

Additional file 1:

Table S1. Primers designed to amplify a segment of the Orientia 47-kD antigen gene and a portion of the 56-kD protein gene. (DOCX 14 kb)

Additional file 2:

Figure S1. Example of representative PCR amplicons on a 1% DNA agarose gel. Lane 1, 100 bp DNA ladder; lane 2, eschar #8 (Karp-related); lane 3, eschar #2 (Kato-related); lane 4, eschar #12 (Gilliam-related); lane 5, negative control; lane 6, Karp strain positive control; lane 7, Kato strain positive control; lane 8, Gilliam strain positive control. (JPG 9 kb)

Additional file 3:

Figure S2. Analysis of nucleotide sequence alignment showing the percent identity and divergence as each eschar and buffy coat nucleotide sequence is compared to one another. (JPG 121 kb)

Additional file 4:

Figure S3. Graphical report showing matching versus differing residues within the nucleotide sequence alignment generated from eschars and buffy coats. (JPG 210 kb)

Rights and permissions

Open Access This 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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Premaratna, R., Blanton, L.S., Samaraweera, D.N. et al. Genotypic characterization of Orientia tsutsugamushi from patients in two geographical locations in Sri Lanka. BMC Infect Dis 17, 67 (2017). https://doi.org/10.1186/s12879-016-2165-z

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12879-016-2165-z

Keywords