Amplicon-Based Next Generation Sequencing for Rapid Identification of Rickettsia and Ectoparasite Species from Entomological Surveillance in Thailand.

BACKGROUND: Next generation sequencing (NGS) technology has been used for a wide range of epidemiological and surveillance studies. Here, we used amplicon-based NGS to species identify Rickettsia and their arthropod hosts from entomological surveillance. METHODS: During 2015-2016, we screened 1825 samples of rodents and ectoparasites collected from rodents and domestic mammals (dog, cat, and cattle) across Thailand for Rickettsia. The citrate synthase gene was amplified to identify Rickettsia to species, while the Cytochrome Oxidase subunit I (COI) and subunit II (COII) genes were used as target genes for ectoparasite identification. All target gene amplicons were pooled for library preparation and sequenced with Illumina MiSeq platform. RESULT: The highest percentage of Rickettsia DNA was observed in fleas collected from domestic animals (56%) predominantly dogs. Only a few samples of ticks from domestic animals, rodent fleas, and rodent tissue were positive for Rickettisia DNA. NGS based characterization of Rickettsia by host identified Rickettsia asembonensis as the most common bacteria in positive fleas collected from dogs (83.2%) while "Candidatus Rickettsia senegalensis" was detected in only 16.8% of Rickettsia positive dog fleas. Sequence analysis of COI and COII revealed that almost all fleas collected from dogs were Ctenocephalides felis orientis. Other Rickettsia species were detected by NGS including Rickettsia heilongjiangensis from two Haemaphysalis hystricis ticks, and Rickettsia typhi in two rodent tissue samples. CONCLUSION: This study demonstrates the utility of NGS for high-throughput sequencing in the species characterization/identification of bacteria and ectoparasite for entomological surveillance of rickettsiae. A high percentage of C. f. orientis are positive for R. asembonensis. In addition, our findings indicate there is a risk of tick-borne Spotted Fever Group rickettsiosis, and flea-borne murine typhus transmission in Tak and Phangnga provinces of Thailand.

Metagenomic Approach to Characterizing Disease Epidemiology in a Disease-Endemic Environment in Northern Thailand.

In this study, we used a metagenomic approach to analyze bacterial communities from diverse populations (humans, animals, and vectors) to investigate the role of these microorganisms as causative agents of disease in human and animal populations. Wild rodents and ectoparasites were collected from 2014 to 2018 in Nan province, Thailand where scrub typhus is highly endemic. Samples from undifferentiated febrile illness (UFI) patients were obtained from a local hospital. A total of 200 UFI patient samples were obtained and 309 rodents and 420 pools of ectoparasites were collected from rodents (n = 285) and domestic animals (n = 135). The bacterial 16S rRNA gene was amplified and sequenced with the Illumina. Real-time PCR and Sanger sequencing were used to confirm the next-generation sequencing (NGS) results and to characterize pathogen species. Several pathogens were detected by NGS in all populations studied and the most common pathogens identified included Bartonella spp., Rickettsia spp., Leptospira spp., and Orientia tsutsugamushi. Interestingly, Anaplasma spp. was detected in patient, rodent and tick populations, although they were not previously known to cause human disease from this region. Candidatus Neoehrlichia, Neorickettsia spp., Borrelia spp., and Ehrlichia spp. were detected in rodents and their associated ectoparasites. The same O. tsutsugamushi genotypes were shared among UFI patients, rodents, and chiggers in a single district indicating that the chiggers found on rodents were also likely responsible for transmitting to people. Serological testing using immunofluorescence assays in UFI samples showed high prevalence (IgM/IgG) of Rickettsia and Orientia pathogens, most notably among samples collected during September-November. Additionally, a higher number of seropositive samples belonged to patients in the working age population (20-60 years old). The results presented in this study demonstrate that the increased risk of human infection or exposure to chiggers and their associated pathogen (O. tsutsugamushi) resulted in part from two important factors; working age group and seasons for rice cultivation and harvesting. Evidence of pathogen exposure was shown to occur as there was seropositivity (IgG) in UFI patients for bartonellosis as well as for anaplasmosis. Using a metagenomic approach, this study demonstrated the circulation and transmission of several pathogens in the environment, some of which are known causative agents of illness in human populations.

Longitudinal study on enterovirus A71 and coxsackievirus A16 genotype/subgenotype replacements in hand, foot and mouth disease patients in Thailand, 2000-2017.

BACKGROUND: Enterovirus A71 (EV-A71) and coxsackievirus A16 (CV-A16) are the major causative agents of hand, foot and mouth disease (HFMD) worldwide, particularly in the Asia-Pacific region. Several strains have emerged, circulated, and faded out over time in recent decades. This study investigated the EV-A71 and CV-A16 circulating strains and replacement of genotypes/subgenotypes in Thailand during the years 2000-2017. METHODS: The complete VP1 regions of 92 enteroviruses obtained from 90 HFMD patients, one asymptomatic adult contact case, and one encephalitic case were sequenced and investigated for serotypes, genotypes, and subgenotypes using a phylogenetic analysis. RESULTS: The 92 enterovirus isolates were identified as 67 (72.8%) EV-A71 strains comprising subgenotypes B4, B5, C1, C2, C4a, C4b and C5, and 25 (27.2%) CV-A16 strains comprising subgenotypes B1a and B1b. Genotypic/subgenotypic replacements were evidenced during the study period. EV-A71 B5 and C4a have been the major circulating strains in Thailand for more than a decade, and CV-A16 B1a has been circulating for almost two decades. CONCLUSIONS: This study provides chronological data on the molecular epidemiology of EV-A71 and CV-A16 subgenotypes in Thailand. Subgenotypic replacement frequently occurred with EV-A71, but not CV-A16. Monitoring for viral genetic and subgenotypic changes is important for molecular diagnosis, vaccine selection, and vaccine development.

Heterogeneity of Orientia tsutsugamushi genotypes in field-collected trombiculid mites from wild-caught small mammals in Thailand.

Trombiculid mites are the vectors of scrub typhus, with infected larval mites (chiggers) transmitting the causative agent, Orientia tsutsugamushi, during feeding. Co-existence of multiple O. tsutsugamushi strains within infected mites has previously been reported in naturally infected, laboratory-reared mite lines using molecular methods to characterize the 56-kDa type-specific antigen (TSA) gene. In the current study, more advanced next-generation sequencing technology was used to reveal the heterogeneity of O. tsutsugamushi genotypes in field-collected trombiculid mites from rodents and small mammals in scrub typhus-endemic areas of Thailand. Twenty-eight trombiculid mites collected from 10 small mammals were positive for O. tsutsugamushi, corresponding to a prevalence rate of 0.7% within the mite population. Twenty-four of the infected mites were Leptotrombidium spp., indicating that this genus is the main vector for O. tsutsugamushi transmission in Thailand. In addition, O. tsutsugamushi was detected in the mite genera Ascoschoengastia, Blankaartia, Gahrliepia, and Lorillatum. Of the 10 infested small animal hosts, six had 2-10 infected mites feeding at the time of collection. Deep sequencing was used to characterize mixed infections (two to three O. tsutsugamushi genotypes within an individual mite), and 5 of the 28 infected mites (17.9%) contained mixed infections. Additionally, 56-kDa TSA gene sequence analysis revealed identical bacterial genotypes among co-feeding mites with single or mixed infections. These results suggest that co-feeding transmission may occur during the feeding process, and could explain the occurrence of mixed infections in individual mites, as well as the recovery of multiple infected mites from the same host. This study also revealed highly diverse within-host O. tsutsugamushi genotypes. The occurrence of multiple O. tsutsugamushi genotypes within individual mites has important implications, and could provide a mechanism for pathogen evolution/diversification in the mite vector.

Francisella-Like Endosymbiont Detected in Haemaphysalis Ticks (Acari: Ixodidae) From the Republic of Korea.

A total of 6,255 ticks belonging to three genera and six species (Haemaphysalis flava Neumann, Haemaphysalis longicornis Neumann, Haemaphysalis phasiana Saito, Ixodes nipponensis Kitaoka & Saito, Ixodes persulcatus Schulze, and Amblyomma testudinarium Koch) collected from May-August, 2013, at four southwestern provinces in the Republic of Korea (ROK) were submitted to the Armed Forces Research Institute of Medical Sciences and assayed for selected tick-borne pathogens. One pool each of H. flava and H. phasiana was positive by PCR and sequencing for a Francisella-like endosymbiont, while all pools were negative for Francisella tularensis, the causative agent of tularemia.

A novel flavivirus detected in two Aedes spp. collected near the demilitarized zone of the Republic of Korea.

Flaviviruses comprise a large and diverse group of positive-stranded RNA viruses, including tick-, mosquito- and unknown-vector-borne flaviviruses. A novel flavivirus was detected in pools of Aedes vexans nipponii (n=1) and Aedes esoensis (n=3) collected in 2012 and 2013 near the demilitarized zone (DMZ), Republic of Korea (ROK). Phylogenetic analyses of the NS5, E gene and complete polyprotein coding sequence (CDS) showed that the novel virus fell within the Aedes-borne flaviviruses (ABFVs), with nucleotide identity ranging from 57.8-75.1 %, 46.1-74.2 % and 51.1-76.2 %, respectively. While the novel ABFV was distant from other flaviviruses within the group, it formed a clade with Ilomantsi virus (ILOV). Sequence alignments of the partial NS5 gene, full-length E gene and polyprotein CDS between the novel virus and ILOV showed approximately 76.2 % nucleotide identity and 90 % amino acid identity, respectively. The ABFV identified in Aedes mosquitoes from the ROK is a novel ABFV based on the sequence analyses and is designated as Panmunjeom flavivirus (PANFV).

Transovarial Transmission of Co-Existing Orientia tsutsugamushi Genotypes in Laboratory-Reared Leptotrombidium imphalum.

The co-evolution of Orientia tsutsugamushi and its vector/host Leptotrombidium mites is important for this bacterium to survive and exist in its environment. The data in this study demonstrated that O. tsutsugamushi has adapted to take advantage of the parasitic nature of the host's larval stage and thus increase its chance of transmission to a vertebrate host and potentially to other vector mites by increasing its density at the time of transmission. Our data demonstrated that during the larval stage the density of O. tsutsugamushi was at its highest level compared to other life stages (Kruskal-Wallis, p < 0.0001). We further revealed that the different O. tsutsugamushi 56-kDa type-specific antigen (TSA) genotypes within the mite were maintained and preserved during transovarial transmission from the Leptotrombidium imphalum, lines Li-3 and Li-5. No sequence difference of 56-kDa TSA gene (variable domain I-III, 765 bp) was observed between the UT302-like genotype found in mothers and their offspring (100% identity). However, one or two nonsynonymous mutations in the 56-kDa TSA gene were observed in the Karp-like genotypes found in the F1 offspring with a percent difference ranging from 0.13 to 0.26 for nucleotide sequences and from 0.39 to 0.78 for amino acid sequences. Additionally, the composition of co-existing O. tsutsugamushi genotypes was maintained in L. imphalum lines through transsovarial and transstadial transmission processes; however, the proportion of these genotypes in each stage varied (larva, nymph, adult). These results show some of the key characteristics of O. tsutsugamushi maintenance within and transmission among its vector/host L. imphalum.