Emerging Genomic Trends on Rabies Virus in Davao Region, Philippines, 2018-2021.

Rabies, caused by the rabies virus (RABV), remains a significant public health issue in the Philippines despite efforts to control it. To eliminate rabies by 2030, effective surveillance strategies are crucial. In this study, we examined RABV evolution and phylodynamics in the Davao Region using genome sequences from Davao City and nearby provinces. We adapted the RABV ARTIC Protocol for Oxford Nanopore High-Throughput Sequencing to optimize workflow efficiency under limited resources. Comparing new virus samples collected from June 2019 to June 2021 (n = 38) with baseline samples from June 2018 to May 2019 (n = 49), new sub-clades were observed in the phylogenetic tree, suggesting divergence from older variants that were previously undetected. Most of the new viruses belonged to the Asian SEA4_A1.1.1 lineage, but new (SEA4_B1 and SEA4_B1.1) and emerging (SEA4_B1.1_E1) lineages that have never been reported in the Philippines were also identified. The baseline study reported phylogeographic clustering of RABV isolates from the same areas. However, this pattern was disrupted in the current biosurveillance, with variants detected in areas outside the original cluster. Furthermore, our findings revealed significant transmission routes between Davao City and neighboring provinces, contrasting with the predominantly intra-city transmission observed in the baseline study. These results underscore the need for ongoing and timely genomic surveillance to monitor genetic diversity changes and the emergence of novel strains, as well as to track alterations in transmission pathways. Implementing cost-effective next-generation sequencing workflows will facilitate the integration of genomic surveillance into rabies control programs, particularly in resource-limited settings. Collaborations between different sectors can empower local laboratories and experts in genomic technologies and analysis.

In silico evaluation of Philippine Natural Products against SARS-CoV-2 Main Protease.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the causative agent of COVID-19, is a novel strain of coronavirus first reported in December 2019 which rapidly spread throughout the world and was subsequently declared a pandemic by the World Health Organization (WHO) in March 2020. Although vaccines, as well as treatments, have been rapidly developed and deployed, these are still spread thin, especially in the developing world. There is also a continuing threat of the emergence of mutated variants which may not be as responsive to available vaccines and drugs. Accessible and affordable sources of antiviral drugs against SARS-CoV-2 offer wider options for the clinical treatment of populations at risk for severe COVID-19. Using in silico methods, this study identified potential inhibitors against the SARS-CoV-2 main protease (Mpro), the protease directly responsible for the activation of the viral replication enzyme, from a consolidated database of 1516 Philippine natural products. Molecular docking experiments, along with in silico ADME predictions, determined top ligands from this database with the highest potential inhibitory effects against Mpro. Molecular dynamic trajectories of the apo and diosmetin-7-O-b-D-glucopyranoside (DG) in complex with the protein predicted potential mechanisms of action for the ligand-by separating the Cys145-His41 catalytic dyad and by influencing the protein network through key intra-signaling residues within the Mpro binding site. These findings show the inhibitory potential of DG against the SARS-CoV-2 Mpro, and further validation is recommended through in vitro or in vivo experimentation.

Multifaceted Assessment of Wastewater-Based Epidemiology for SARS-CoV-2 in Selected Urban Communities in Davao City, Philippines: A Pilot Study.

Over 60 countries have integrated wastewater-based epidemiology (WBE) in their COVID-19 surveillance programs, focusing on wastewater treatment plants (WWTP). In this paper, we piloted the assessment of SARS-CoV-2 WBE as a complementary public health surveillance method in susceptible communities in a highly urbanized city without WWTP in the Philippines by exploring the extraction and detection methods, evaluating the contribution of physico-chemical-anthropogenic factors, and attempting whole-genome sequencing (WGS). Weekly wastewater samples were collected from sewer pipes or creeks in six communities with moderate-to-high risk of COVID-19 transmission, as categorized by the City Government of Davao from November to December 2020. Physico-chemical properties of the wastewater and anthropogenic conditions of the sites were noted. Samples were concentrated using a PEG-NaCl precipitation method and analyzed by RT-PCR to detect the SARS-CoV-2 N, RdRP, and E genes. A subset of nine samples were subjected to WGS using the Minion sequencing platform. SARS-CoV-2 RNA was detected in twenty-two samples (91.7%) regardless of the presence of new cases. Cycle threshold values correlated with RNA concentration and attack rate. The lack of a sewershed map in the sampled areas highlights the need to integrate this in the WBE planning. A combined analysis of wastewater physico-chemical parameters such as flow rate, surface water temperature, salinity, dissolved oxygen, and total dissolved solids provided insights on the ideal sampling location, time, and method for WBE, and their impact on RNA recovery. The contribution of fecal matter in the wastewater may also be assessed through the coliform count and in the context of anthropogenic conditions in the area. Finally, our attempt on WGS detected single-nucleotide polymorphisms (SNPs) in wastewater which included clinically reported and newly identified mutations in the Philippines. This exploratory report provides a contextualized framework for applying WBE surveillance in low-sanitation areas.

The first complete mitochondrial genome sequence of Cynopterusbrachyotis (Chiroptera, Pteropodidae) from the Philippines.

The technical limitations of capillary sequencing in providing insights on phylogeny have been greatly aided in recent years by the implementation of next generation sequencing platforms which can generate whole mitochondrial genome (mitogenome) sequences. In this study, enriched mitochondrial DNA of Cynopterusbrachyotis from Mindanao, Philippines was sequenced using the Illumina MiSeq platform. A total of 653,967 clean paired-end reads was assembled using a MIRA-MITObim pipeline, resulting in a consensus mitogenome sequence length of 17,382 bases and a GC content of 41.48%, which is consistent with other published mitogenomes in fruit bats. The assembled C.brachyotis mitogenome was annotated using the MITOS online server and was able to resolve all mitochondrial genes, except for one transfer RNA gene (trnT) which may be further resolved by additional capillary sequencing of the region. Sequence analysis showed that the Philippine C.brachyotis is only 90%-91% homologous with other Cynopterus spp., based on its full mitogenome sequence. Phylogenetic analysis of fruit bat mitogenomes, deposited in online repositories, revealed that the Philippine C.brachyotis in this study has diverged from Asian Cynopterus, namely Cynopterusbrachyotis and Cynopterussphinx from other parts of Asia (100% bootstrap support) with the latter two forming a separate clade. This divergence at the species level was consistent with phylogentic inference using cytochrome oxidase 1 (CO1) and cytochrome B (cytb) gene markers. Our results strengthen the previously reported hypothesis that the Cynopteruscf.brachyotis in the Philippines is distinct from its Asian counterparts and should be, therefore, treated as a new species.

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.

A cointegration analysis of rabies cases and weather components in Davao City, Philippines from 2006 to 2017.

Rabies is a lethal viral disease and dogs are the major disease reservoir in the Philippines. Spatio-temporal variations in environmental factors are known to affect disease dynamics. Some rabies-affected countries considered investigating the role of weather components in driving rabies cases and it has helped them to strategize their control efforts. In this study, cointegration analysis was conducted between the monthly reported rabies cases and the weather components, such as temperature and precipitation, to verify the effect of weather components on rabies incidence in Davao City, Philippines. With the Engle-Granger cointegration tests, we found that rabies cases are cointegrated into each of the weather components. It was further validated, using the Granger causality test, that each weather component predicts the rabies cases and not vice versa. Moreover, we performed the Johansen cointegration test to show that the weather components simultaneously affect the number of rabies cases, which allowed us to estimate a vector-error correction model for rabies incidence as a function of temperature and precipitation. Our analyses showed that canine rabies in Davao City was weather-sensitive, which implies that rabies incidence could be projected using established long-run relationship among reported rabies cases, temperature, and precipitation. This study also provides empirical evidence that can guide local health officials in formulating preventive strategies for rabies control and eradication based on weather patterns.

Spatiotemporal dynamics and risk factors of rotavirus A circulation in backyard pig farms in a Philippine setting.

Rotavirus A is one of the leading etiological agents of porcine gastroenteritis, a condition which results to stunted growth among piglets. Moreover, there is increasing evidence for zoonosis of rotavirus A (RVA), which is also the principal cause of diarrhea in children. In the absence of rigorous animal health monitoring in Philippine backyard farms, there is therefore a strong need for RVA surveillance. In this study, 30 randomly selected backyard farms were subjected to surveillance for RVA for 12 months. Results show that RVA detection at a monthly farm-level rate ranged from 0 to 52%, with an overall average of 23%. RVA had higher detection rates in adult pigs compared to young piglets and was most prevalent in non-diarrheic stools, indicating asymptomatic circulation of the virus. Spatiotemporal analysis demonstrated that the viral circulation exhibits a seasonal pattern that peaks and forms geographical clusters during the cooler months of the year, suggesting farm-to-farm transmission. Risk factor analysis identified specific farm conditions that increase the likelihood of RVA circulation: presence of gilts, larger herd size, presence of other animals, and abiotic factors such as low relative humidity and low altitude. The same analysis also revealed three major management practices that can help reduce the pressure of infection in these farms: sanitation and waste disposal, animal grouping, and diet. This new perspective on porcine RVA circulation will benefit the underprivileged backyard farmers and help empower them to protect both animal and public health.