Genetic diversity and geographic distribution of genetically distinct rabies viruses in the Philippines

Molecular analysis and geographic distribution of the recent Indonesian rabies virus

Background and Aim

Some Indonesian islands, including Sumatra, Kalimantan, Sulawesi, Java, and East Nusa Tenggara, have endemic rabies. Rabies outbreaks in Bali began from 2008 to 2011 and continue to occur sporadically. This study aimed to study the molecular analysis and geographical distribution of Indonesian rabies virus (RABV) from 2016 to 2021 and compare to previous periods.

Materials and Methods

Virus isolates from 2016 to 2021 were extracted from dog brains and sequenced at the nucleoprotein gene locus. They were compared with data sequences available in the GenBank database. Indonesian RABV from the previous three periods (before 1989, 1997-2003, and 2008-2010) was extracted from the GenBank database. The genetic diversity in this study was based on the N gene of Indonesian RABV.

Results

Asian RABV, which is genetically close to the Indonesian virus, is a virus from China (ASIA-3 cluster) and from the Southeast Asia region, namely, virus isolates from Sarawak and Malaysia and some Cambodian isolates. Rabies virus, which was isolated from the Bali islands, was the new cluster first detected and published in Bali, Indonesia, in 2008, while RABV from West Sumatra Province, which was isolated from 2016 to 2021, was also considered a new cluster that is genetically distant from other clusters in Indonesia.

Conclusion

The RABV in Indonesia is divided into five clusters. The isolates from West Sumatra Province from 2016 to 2021 were a new cluster genetically distant from other Indonesian viruses.

Impact and mitigation of sampling bias to determine viral spread: Evaluating discrete phylogeography through CTMC modeling and structured coalescent model approximations

Bayesian phylogeographic inference is a powerful tool in molecular epidemiological studies, which enables reconstruction of the origin and subsequent geographic spread of pathogens. Such inference is, however, potentially affected by geographic sampling bias. Here, we investigated the impact of sampling bias on the spatiotemporal reconstruction of viral epidemics using Bayesian discrete phylogeographic models and explored different operational strategies to mitigate this impact. We considered the continuous-time Markov chain (CTMC) model and two structured coalescent approximations (Bayesian structured coalescent approximation [BASTA] and marginal approximation of the structured coalescent [MASCOT]). For each approach, we compared the estimated and simulated spatiotemporal histories in biased and unbiased conditions based on the simulated epidemics of rabies virus (RABV) in dogs in Morocco. While the reconstructed spatiotemporal histories were impacted by sampling bias for the three approaches, BASTA and MASCOT reconstructions were also biased when employing unbiased samples. Increasing the number of analyzed genomes led to more robust estimates at low sampling bias for the CTMC model. Alternative sampling strategies that maximize the spatiotemporal coverage greatly improved the inference at intermediate sampling bias for the CTMC model, and to a lesser extent, for BASTA and MASCOT. In contrast, allowing for time-varying population sizes in MASCOT resulted in robust inference. We further applied these approaches to two empirical datasets: a RABV dataset from the Philippines and a SARS-CoV-2 dataset describing its early spread across the world. In conclusion, sampling biases are ubiquitous in phylogeographic analyses but may be accommodated by increasing the sample size, balancing spatial and temporal composition in the samples, and informing structured coalescent models with reliable case count data.

Clinical, epidemiological, and spatial features of human rabies cases in Metro Manila, the Philippines from 2006 to 2015

Rabies remains a public health problem in the Philippines despite the widespread provision of rabies vaccines and rabies immunoglobulin (RIG) as post-exposure prophylaxis (PEP). Detailed descriptions of recent human rabies cases in the Philippines are scarce. This study aimed to describe the clinical, epidemiological, and spatial features of human rabies cases between January 1, 2006, and December 31, 2015. We conducted a retrospective hospital-based case record review of all patients admitted to one referral hospital in Manila who received a clinical diagnosis of rabies. During the 10-year study period there were 575 patients (average 57.5 cases per year, range 57 to 119) with a final diagnosis of rabies. Most patients were male (n = 404, 70.3%) and aged ≥ 20 years (n = 433, 75.3%). Patients mostly came from the National Capital Region (n = 160, 28.0%) and the adjacent Regions III (n = 197, 34.4%) and IV-A (n = 168, 29.4%). Case mapping and heatmaps showed that human rabies cases were continuously observed in similar areas throughout the study period. Most patients had hydrophobia (n = 444, 95.5%) and/or aerophobia (n = 432, 93.3%). The leading causative animals were dogs (n = 421, 96.3%) and cats (n = 16, 3.7%). Among 437 patients with animal exposure history, only 42 (9.6%) had been administered at least one rabies vaccine. Two patients (0.5%), young children bitten on their face, had received and a full course of rabies vaccine. Human rabies patients were continuously admitted to the hospital, with no notable decline over the study period. The geographical area in which human rabies cases commonly occurred also did not change. Few patients received PEP and there were two suspected cases of PEP failure. The retrospective design of this study was a limitation; thus, prospective studies are required.

Rabies remains a public health problem in the Philippines despite improvements in the availability of rabies vaccines and rabies immunoglobulin (RIG) as post-exposure prophylaxis (PEP). The incidence of rabies is highest in Metro Manila and surrounding areas. We reviewed the records of all human rabies patients admitted to the national infectious disease hospital in Manila between 2006 and 2015. This hospital treats most cases in this area. During the 10-year study period, human rabies cases were continuously admitted to the hospital, with no notable decline in numbers by year. Most patients were adult men bitten by domestic dogs. The geographical areas in which cases commonly occurred during the 10-year period also did not change over time. Only 9.6% of patients had received at least one dose of a rabies vaccine as PEP. Although the risk of PEP failure is reported to be almost zero, we identified two suspected cases of PEP failure. The retrospective design of this study was a limitation, and the exact details of PEP were not reliably available. As human rabies death is a significant public health concern, the circumstances of each case should be prospectively investigated. Further research is required to understand how to reduce the number of rabies cases.

Molecular Epidemiology of Rabies in Wild Canidae in Tunisia

Rabies is a viral zoonosis that is transmissible to humans via domestic and wild animals. There are two epidemiological cycles for rabies, the urban and the sylvatic cycles. In an attempt to study the epidemiological role of wild canidae in rabies transmission, the present study aimed to analyze the genetic characteristics of virus isolates and confirm prior suggestions that rabies is maintained through a dog reservoir in Tunisia. Virus strains isolated from wild canidae were subject to viral sequencing, and Bayesian phylogenetic analysis was performed using Beast2 software. Essentially, the virus strains isolated from wild canidae belonged to the Africa-1 clade, which clearly diverges from fox-related strains. Our study also demonstrated that genetic characteristics of the virus isolates were not as distinct as could be expected if a wild reservoir had already existed. On the contrary, the geographic landscape is responsible for the genetic diversity of the virus. The landscape itself could have also acted as a natural barrier to the spread of the virus.

On the Use of Phylogeographic Inference to Infer the Dispersal History of Rabies Virus: A Review Study

Rabies is a neglected zoonotic disease which is caused by negative strand RNA-viruses belonging to the genus Lyssavirus. Within this genus, rabies viruses circulate in a diverse set of mammalian reservoir hosts, is present worldwide, and is almost always fatal in non-vaccinated humans. Approximately 59,000 people are still estimated to die from rabies each year, leading to a global initiative to work towards the goal of zero human deaths from dog-mediated rabies by 2030, requiring scientific efforts from different research fields. The past decade has seen a much increased use of phylogeographic and phylodynamic analyses to study the evolution and spread of rabies virus. We here review published studies in these research areas, making a distinction between the geographic resolution associated with the available sequence data. We pay special attention to environmental factors that these studies found to be relevant to the spread of rabies virus. Importantly, we highlight a knowledge gap in terms of applying these methods when all required data were available but not fully exploited. We conclude with an overview of recent methodological developments that have yet to be applied in phylogeographic and phylodynamic analyses of rabies virus.

Genome-based local dynamics of canine rabies virus epidemiology, transmission, and evolution in Davao City, Philippines, 2018-2019

Rabies is a fatal zoonotic and neglected tropical disease caused by the rabies virus (RABV) and is associated with neuronal dysfunction and death, with dogs as the predominant carrier. The Philippines plans to eradicate rabies by 2022, but this is challenged with sub-optimal coverage of vaccination programs coupled with sustained transmission chains, making it unable to eradicate the disease. We investigated the dynamics of canine rabies in the highly urbanized Davao City of the Philippines and its neighboring localities by assessing genetic relationships, transmission patterns, selection pressure, and recombination events using the whole genome sequence of 49 RABV cases from June 2018 to May 2019, majority of which (46%) were from the district of Talomo, Davao City. Although phylogeographic clustering was observed, local variants also exhibited genetic sub-lineages. Phylogenetic and spatial transmission analysis provided evidence for intra- and inter-city transmission predominantly through the Talomo district of Davao City. Around 84% of the cases were owned dogs, but the genetic similiarity of RABVs from stray and owned dogs further alluded to the role of the former as transmission vectors. The high rate of improper vaccination among the affected dogs (80%) was also a likely contributor to transmission. The RABV population under Investigation is generally under strong purifying selection with no evidence of vaccine evasion due to the genetic homogeneity of viruses from vaccinated and improperly vaccinated dogs. However, some homologous recombination (HR) events were identified along the G and L genes, also predominantly associated with viruses from Talomo. The complementary findings on epidemiology, transmission, and recombination for Talomo suggest that high incidence areas can be seeds for virus dispersal and evolution. We recommend further Investigations on the possibility of HR in future large-scale genome studies. Finally, districts associated with these phenomena can be targeted for evidence-based local strategies that can help break RABV transmission chains and prevent emergence of novel strains in Davao City.

Evaluation of the diagnostic accuracy of lateral flow devices as a tool to diagnose rabies in post-mortem animals

Implementation of lateral flow devices (LFDs) for rabies antigen detection is expected to improve surveillance through the efficient detection of rabid animals in resource-limited settings; however, the use of LFDs for diagnosis remains controversial because some commercially available kits show low sensitivity. Therefore, we compared the diagnostic efficacy of three LFDs (ADTEC, Bionote, and Elabscience kits) paralleled with the direct fluorescent antibody test (dFAT) using fresh samples and investigated the diagnostic accuracies. To do so, we evaluated rabies-suspected samples submitted to the Regional Animal Disease Diagnostic Laboratory III, Philippines. Furthermore, we conducted real-time RT-PCR and sequencing to measure the accuracy of field laboratory diagnosis. The total number of animals submitted during this study period was 184 cases, including negative control samples. Of these, 53.9% (84 cases) were positive in the dFAT. Dogs were the most common rabies-suspected animal (n = 135). The sensitivities of the ADTEC and Bionote kits were 0.88 (74 cases) and 0.95 (80 cases), respectively. The specificity of both kits was 1.00 (100 cases). Furthermore, the sensitivity and specificity of the ADTEC kit after directly homogenizing the samples in assay buffer without dilution in phosphate-buffered saline (ADTEC kit DM) were 0.94 (79 cases) and 1.00 (100 cases), respectively. By contrast, there were no positive results using the Elabscience kit among all dFAT-positive samples. The sensitivity and specificity of LFDs make these tests highly feasible if properly used. Therefore, LFD tests can be used to strengthen the surveillance of rabies-infected animals in endemic and resource-limited settings.

Spatio-temporal evolution of canine rabies in Tunisia, 2011-2016

Tunisia is an endemic country for dog mediated rabies. An increase in canine rabies cases during the last decade has been suspected. Since no studies have been conducted on rabies spatial distribution, the present work was focused on spatiotemporal evolution of rabies in Tunisia during the 2011-2016 period with a special focus on the reservoir species. Data collected concerned suspected dogs that originate from the whole country. Surveillance indicators such as positive fractions and number of suspected dogs received at the laboratory have been calculated. Spatiotemporal hotspots were then mapped, spatial and spatio-temporal analysis were carried out using discrete Poisson spatial model and space-time permutation models available in SaTScan9 software. The study revealed that an actual increase in canine rabies incidence occurred in Tunisia since 2012. Spatial and spatio-temporal analysis identified clusters centered in the North and in the Center East of the country. Spatio-temporal clusters were non overlapping, indicating that this spatial distribution is not fixed through time. A large heterogeneity in surveillance indicators such as number of suspected dogs was associated to the distance to the laboratory or to insufficient coordination between governorates.

Rapid in-country sequencing of whole virus genomes to inform rabies elimination programmes

Genomic surveillance is an important aspect of contemporary disease management but has yet to be used routinely to monitor endemic disease transmission and control in low- and middle-income countries. Rabies is an almost invariably fatal viral disease that causes a large public health and economic burden in Asia and Africa, despite being entirely vaccine preventable. With policy efforts now directed towards achieving a global goal of zero dog-mediated human rabies deaths by 2030, establishing effective surveillance tools is critical. Genomic data can provide important and unique insights into rabies spread and persistence that can direct control efforts. However, capacity for genomic research in low- and middle-income countries is held back by limited laboratory infrastructure, cost, supply chains and other logistical challenges. Here we present and validate an end-to-end workflow to facilitate affordable whole genome sequencing for rabies surveillance utilising nanopore technology. We used this workflow in Kenya, Tanzania and the Philippines to generate rabies virus genomes in two to three days, reducing costs to approximately £60 per genome. This is over half the cost of metagenomic sequencing previously conducted for Tanzanian samples, which involved exporting samples to the UK and a three- to six-month lag time. Ongoing optimization of workflows are likely to reduce these costs further. We also present tools to support routine whole genome sequencing and interpretation for genomic surveillance. Moreover, combined with training workshops to empower scientists in-country, we show that local sequencing capacity can be readily established and sustainable, negating the common misperception that cutting-edge genomic research can only be conducted in high resource laboratories. More generally, we argue that the capacity to harness genomic data is a game-changer for endemic disease surveillance and should precipitate a new wave of researchers from low- and middle-income countries.

Rapid in-country sequencing of whole virus genomes to inform rabies elimination programmes

Genomic surveillance is an important aspect of contemporary disease management but has yet to be used routinely to monitor endemic disease transmission and control in low- and middle-income countries. Rabies is an almost invariably fatal viral disease that causes a large public health and economic burden in Asia and Africa, despite being entirely vaccine preventable. With policy efforts now directed towards achieving a global goal of zero dog-mediated human rabies deaths by 2030, establishing effective surveillance tools is critical. Genomic data can provide important and unique insights into rabies spread and persistence that can direct control efforts. However, capacity for genomic research in low- and middle-income countries is held back by limited laboratory infrastructure, cost, supply chains and other logistical challenges. Here we present and validate an end-to-end workflow to facilitate affordable whole genome sequencing for rabies surveillance utilising nanopore technology. We used this workflow in Kenya, Tanzania and the Philippines to generate rabies virus genomes in two to three days, reducing costs to approximately £60 per genome. This is over half the cost of metagenomic sequencing previously conducted for Tanzanian samples, which involved exporting samples to the UK and a three- to six-month lag time. Ongoing optimization of workflows are likely to reduce these costs further. We also present tools to support routine whole genome sequencing and interpretation for genomic surveillance. Moreover, combined with training workshops to empower scientists in-country, we show that local sequencing capacity can be readily established and sustainable, negating the common misperception that cutting-edge genomic research can only be conducted in high resource laboratories. More generally, we argue that the capacity to harness genomic data is a game-changer for endemic disease surveillance and should precipitate a new wave of researchers from low- and middle-income countries.

The phylodynamics of the rabies virus in the Russian Federation

Near complete rabies virus N gene sequences (1,110 nt) were determined for 82 isolates obtained from different regions of Russia between 2008 and 2016. These sequences were analyzed together with 108 representative GenBank sequences from 1977-2016 using the Bayesian coalescent approach. The timing of the major evolutionary events was estimated. Most of the isolates represented the steppe rabies virus group C, which was found over a vast geographic region from Central Russia to Mongolia and split into three groups (C0-C2) with discrete geographic prevalence. A single strain of the steppe rabies virus lineage was isolated in the far eastern part of Russia (Primorsky Krai), likely as a result of a recent anthropogenic introduction. For the first time the polar rabies virus group A2, previously reported in Alaska, was described in the northern part of European Russia and at the Franz Josef Land. Phylogenetic analysis suggested that all currently circulating rabies virus groups in the Russian Federation were introduced within the few last centuries, with most of the groups spreading in the 20th century. The dating of evolutionary events was highly concordant with the historical epidemiological data.

Epidemiology and transmission characteristics of human adenovirus type 7 caused acute respiratory disease outbreak in military trainees in East China

BACKGROUND

Human adenovirus type 7 (HAdV7) is globally attracting great concern as its high morbidity and severity in respiratory diseases, especially in Asia.

OBJECTIVE

To investigate the clinical and epidemiologic characteristics of HAdV7 infection outbreak in East China.

METHODS

The clinical samples were collected from the patients of an ARD outbreak in East Chinafor the detection of causative pathogens by multiplex PCR. The molecular type of human adenovirus isolates were identified by sequencing and homologous comparison based on their hexon genes. The spatiotemporal dynamics of global HAdV7 was investigated using the phylogenetic and phylogeographic analyses. Total 67 referenced HAdV7 hexon sequences (>800 bp) from GenBank were selected for constructing the maximum likelihood tree by MEGA 5.1.0, grouped according to the tree topology for the further migration analysis by PAUP* 4.0 and MigraPhyla 1.0 b to understand the transmission patterns of HAdV7 in global epidemics.

RESULTS

The results showed HAdV7 as the causative pathogen in this outbreak, and the outbreak strains had the hexon sequences highly identical with the isolates in Shaanxi (2012). The origin of HAdV7 was inferred as California, meanwhile a total of 21 migration routes were acquired. HAdV7 in this outbreak was statistically proven dispersed from Shaanxi province (2012).

CONCLUSIONS

The analyses of epidemiology and transmission pattern of HAdV7 would not only enrich the molecular biological basic database but also provide theoretical basis for HAdV7 prevention and control strategy.

Molecular and mathematical modeling analyses of inter-island transmission of rabies into a previously rabies-free island in the Philippines

Rabies is endemic in the Philippines and dog bites are a major cause of rabies cases in humans. The rabies control program has not been successful in eliminating rabies because of low vaccination coverage among dogs. Therefore, more effective and feasible strategies for rabies control are urgently required in the country. To control rabies, it is very important to know if inter-island transmission can occur because rabies can become endemic once the virus is introduced in areas that previously had no reported cases. Our molecular epidemiological study suggests that inter-island transmission events can occur; therefore, we further investigated these inter-island transmission using phylogenetic and modeling approaches. We investigate inter-island transmission between Luzon and Tablas Islands in the Philippines. Phylogenetic analysis and mathematical modeling demonstrate that there was a time lag of several months to a year from rabies introduction to initial case detection, indicating the difficulties in recognizing the initial rabies introductory event. There had been no rabies cases reported in Tablas Island; however, transmission chain was sustained on this island after the introduction of rabies virus because of low vaccination coverage among dogs. Across the islands, a rabies control program should include control of inter-island dog transportation and rabies vaccination to avoid viral introduction from the outside and to break transmission chains after viral introduction. However, this program has not yet been completely implemented and transmission chains following inter-island virus transmission are still observed. Local government units try to control dog transport; however, it should be more strictly controlled, and a continuous rabies control program should be implemented to prevent rabies spread even in rabies-free areas.

Emergence of a sylvatic enzootic formosan ferret badger-associated rabies in Taiwan and the geographical separation of two phylogenetic groups of rabies viruses

Taiwan had been declared rabies-free in humans and domestic animals for five decades until July 2013, when surprisingly, three Formosan ferret badgers (FB) were diagnosed with rabies. Since then, a variety of wild carnivores and other wildlife species have been found dead, neurologically ill, or exhibiting aggressive behaviors around the island. To determine the affected animal species, geographic areas, and environments, animal bodies were examined for rabies by direct fluorescent antibody test (FAT). The viral genomes from the brains of selected rabid animals were sequenced for the phylogeny of rabies viruses (RABV). Out of a total of 1016 wild carnivores, 276/831 (33.2%) Formosan FBs were FAT positive, with occasional biting incidents in 1 dog and suspected spillover in 1 house shrew. All other animals tested, including dogs, cats, bats, mice, house shrews, and squirrels, were rabies-negative. The rabies was badger-associated and confined to nine counties/cities in sylvatic environments. Phylogeny of nucleoprotein and glycoprotein genes from 59 Formosan FB-associated RABV revealed them to be clustered in two distinct groups, TWI and TWII, consistent with the geographic segregation into western and eastern Taiwan provided by the Central Mountain Range and into northern rabies-free and central-southern rabies-affected regions by a river bisecting western Taiwan. The unique features of geographic and genetic segregation, sylvatic enzooticity, and FB-association of RABV suggest a logical strategy for the control of rabies in this nation.

Rabies Virus Infection in Ferret Badgers (Melogale moschata subaurantiaca) in Taiwan: A Retrospective Study

Fifteen ferret badgers (Melogale moschata subaurantiaca), collected 2010-13 and stored frozen, were submitted for rabies diagnosis by direct fluorescent antibody test and reverse transcription PCR. We detected seven positive animal samples, including some from 2010, which indicated that the ferret badger population in Taiwan had been affected by rabies prior to 2010.

[Researches on virology at the Tohoku University Research Center in the Philippines]

Tohoku University Graduate School of Medicine has established the Tohoku-RITM Collaborative Research Center on Emerging and Re-emerging Diseases at Research Institute for Tropical Medicine (RITM) in the Philippines in 2008. Our aim of the center is to conduct operational researches, which can contribute to control of infectious diseases in the Philippines. Therefore most of our researches in the Philippines are being conducted in the fields. Main research themes include severe acute respiratory infections in children, influenza disease burden study, molecular epidemiology of rabies, and viral etiology of acute diarrhea. The study on severe acute respiratory infections in children in Leyte Island has recruited hospitalized cases with severe pneumonia. We showed that enterovirus 68 was one of important causative agents in severe pneumonia cases. We also conducted other analyses including molecular epidemiology of respiratory syncytial virus (RSV) and pathogenesis of human rhinoviruses (HRV). Based on these studies, we initiated more comprehensive researches in the Philippines since 2010.

Molecular characterization of cryptically circulating rabies virus from ferret badgers, Taiwan

After the last reported cases of rabies in a human in 1959 and a nonhuman animal in 1961, Taiwan was considered free from rabies. However, during 2012-2013, an outbreak occurred among ferret badgers in Taiwan. To examine the origin of this virus strain, we sequenced 3 complete genomes and acquired multiple rabies virus (RABV) nucleoprotein and glycoprotein sequences. Phylogeographic analyses demonstrated that the RABV affecting the Taiwan ferret badgers (RABV-TWFB) is a distinct lineage within the group of lineages from Asia and that it has been differentiated from its closest lineages, China I (including isolates from Chinese ferret badgers) and the Philippines, 158-210 years ago. The most recent common ancestor of RABV-TWFB originated 91-113 years ago. Our findings indicate that RABV could be cryptically circulating in the environment. An understanding of the underlying mechanism might shed light on the complex interaction between RABV and its host.

Phylogeographic analysis of rabies viruses in the Philippines

Rabies still remains a public health threat in the Philippines. A significant number of human rabies cases, about 200-300 cases annually, have been reported, and the country needs an effective strategy for rabies control. To develop an effective control strategy, it is important to understand the transmission patterns of the rabies viruses. We conducted phylogenetic analyses by considering the temporal and spatial evolution of rabies viruses to reveal the transmission dynamics in the Philippines. After evaluating the molecular clock and phylogeographic analysis, we estimated that the Philippine strains were introduced from China around the beginning of 20th century. Upon this introduction, the rabies viruses evolved within the Philippines to form three major clades, and there was no indication of introduction of other rabies viruses from any other country. However, within the Philippines, island-to-island migrations were observed. Since then, the rabies viruses have diffused and only evolved within each island group. The evolutionary pattern of these viruses was strongly shaped by geographical boundaries. The association index statistics demonstrated a strong spatial structure within the island group, indicating that the seas were a significant geographical barrier for viral dispersal. Strong spatial structure was also observed even at a regional level, and most of the viral migrations (79.7% of the total median number) in Luzon were observed between neighboring regions. Rabies viruses were genetically clustered at a regional level, and this strong spatial structure suggests a geographical clustering of transmission chains and the potential effectiveness of rabies control that targets geographical clustering. Dog vaccination campaigns have been conducted independently by local governments in the Philippines, but it could be more effective to implement a coordinated vaccination campaign among neighboring areas to eliminate geographically-clustered rabies transmission chains.