Tropical Infectious Diseases Research & Education Centre

Human arboviral diseases transmitted by Aedes aegypti such as dengue, Zika, chikungunya and yellow fever remain global public health threats to date. Of these diseases, dengue fever is particularly prevalent in Southeast Asia. Relentless vector control efforts are performed to curtail disease transmissions through which pyrethroid insecticides are broadly used as the first line of defense to control Ae. aegypti, especially in the course of disease outbreaks. Here, we compile the largest contemporary database for susceptibility profiles and underlying mechanisms involved in Ae. aegypti resistant to pyrethroids in Southeast Asia. The extensive use of pyrethroids inevitably elicit different levels of resistance to numerous populations despite the presence of geographical isolation. The most common mechanisms of pyrethroid resistance that have been identified in Ae. aegypti includes mutations in the voltage sensitive sodium channel gene (Vssc gene) and metabolic-mediated insecticide resistance. Aedes aegypti develops resistance to pyrethroids by acquisition of one or several amino acid substitution(s) in this Vssc gene. Enzymes involved in metabolic-mediated detoxification (i.e. monooxygenases, glutathione-S-transferases and esterases) have been reported to be related to pyrethroid resistance but many specific contributory enzymes are not completely studied. An inadequate amount of data from some countries indicates an urgent need for further study to fill the knowledge gaps. Perspectives and future research needs are also discussed.

Human metapneumovirus (HMPV) is an important viral respiratory pathogen worldwide. Current knowledge regarding the genetic diversity, seasonality and transmission dynamics of HMPV among adults and children living in tropical climate remains limited. HMPV prevailed at 2.2% (n = 86/3,935) among individuals presented with acute respiratory tract infections in Kuala Lumpur, Malaysia between 2012 and 2014. Seasonal peaks were observed during the northeast monsoon season (November-April) and correlated with higher relative humidity and number of rainy days (P < 0.05). Phylogenetic analysis of the fusion and attachment genes identified the co-circulation of three known HMPV sub-lineages, A2b and B1 (30.2% each, 26/86) and B2 (20.9%, 18/86), with genotype shift from sub-lineage B1 to A2b observed in 2013. Interestingly, a previously unrecognized sub-lineage of A2 was identified in 18.6% (16/86) of the population. Using a custom script for network construction based on the TN93 pairwise genetic distance, we identified up to nine HMPV transmission clusters circulating as multiple sub-epidemics. Although no apparent major outbreak was observed, the increased frequency of transmission clusters (dyads) during seasonal peaks suggests the potential roles of transmission clusters in driving the spread of HMPV. Our findings provide essential information for therapeutic research, prevention strategies, and disease outbreak monitoring of HMPV.

Carbapenem resistant Enterobacteriaceae is a growing concern worldwide including Malaysia. The emergence of this pathogen is worrying because carbapenem is one of the 'last-line' antibiotics. The main objective of this study was to determine the prevalence of genetic mechanisms and clinical risk factors of carbapenem resistant Klebsiella pneumoniae (K. pneumoniae) in Malaysia. In this study, seventeen carbapenem resistant K. pneumoniae strains isolated from a tertiary teaching hospital in 2013 were studied. Minimal inhibitory concentration (MIC) of the bacterial strains was determined and genes associated with carbapenemases and extended-spectrum-beta-lactamases (ESBLs) were sequenced and compared with the closest representatives published in public domains. All strains were also sub-typed using pulsed-field gel electrophoresis (PFGE) and multilocus sequence typing (MLST). Statistical analyses were performed to determine the correlation between risk factors for acquiring carbapenem resistant K. pneumoniae and in-hospital mortality. The predominant carbapenemase was bla OXA-48, detected in 12 strains (70.59%). Other carbapenemases detected in this study were bla KPC-2, bla IMP-8, bla NMC-A and bla NDM-1. Nine different pulsotypes were identified and nine strains which were affiliated with ST101, the predominant sequence type had similar PFGE patterns (similarity index of 85%). Based on univariate statistical analysis, resistance to imipenem and usage of mechanical ventilation showed a statistically significant effect separately to in-hospital mortality. The diverse genetic mechanisms harbored by these carbapenem resistant K. pneumoniae facilitates its spread and complicates its detection. Thus, correlation between microbiological trends with host characteristics and clinical factors will provide a better insight of rational treatment strategies and pathogen control.

The current war in Yemen has displaced millions of people from their homes into living in cramped shelters where the healthcare is limited. The breakdown of Yemen's healthcare and sanitation systems has facilitated the spread of infectious diseases including mosquito-borne diseases. The present study aimed to describe the prevalence of dengue virus (DENV) infection among the febrile patients of the Taiz governorate, Yemen as well as their knowledge, attitude and preventive practices (KAPs) regarding dengue fever (DF), and to investigate the factors associated with dengue preventive practices during the war. A total of 384 clinically dengue-suspected patients who sought health care in Taiz, Yemen during the period from July 2016 until October 2016 were recruited for the study. Serum samples were obtained and screened for the presence of DENV RNA and anti-DENV antibodies by reverse transcription-recombinase polymerase amplification (RT-RPA) and dengue IgM/IgG-capture ELISA, respectively. KAP questionnaires were obtained from all participants too. In the study, dengue was laboratory confirmed in approximately 49.3% (189/384) of the clinically suspected dengue patients. In general, 67.1% of the patients had low knowledge scores regarding DF. Low scores for knowledge about DF was significantly associated with those in the age groups of ≤20 years and 21-30 years, illiterates and patients with non-skilled jobs or jobless. The most common preventive practices reported by participants were covering stored water (78.6%) and putting a screen on the house's windows (65.3%). A low proportion of participants (6.7%) had 51-100% of good DF preventive practices. Low scores of positive attitudes toward DF was identified as a risk factor. The study participants showed poor knowledge about DF and their ways of dealing with the various aspects of DF prevention was quite limited, hence, preventive measures against the disease were less likely to be undertaken. Findings from the study highlight the peril of dengue in Taiz, Yemen, which is now comparable to that of endemic regions. The ongoing civil war with disruption in regular health services compounded by the low knowledge about DF as well as the limited DF preventive practices could result in entrenchment of dengue in Yemen.

BACKGROUND: Phenotypic identification systems are established methods for laboratory identification of bacteria causing human infections. Here, the utility of phenotypic identification systems was compared against 16S rDNA identification method on clinical isolates obtained during a 5-year study period, with special emphasis on isolates that gave unsatisfactory identification. METHODS: One hundred and eighty-seven clinical bacteria isolates were tested with commercial phenotypic identification systems and 16S rDNA sequencing. Isolate identities determined using phenotypic identification systems and 16S rDNA sequencing were compared for similarity at genus and species level, with 16S rDNA sequencing as the reference method. RESULTS: Phenotypic identification systems identified ~46% (86/187) of the isolates with identity similar to that identified using 16S rDNA sequencing. Approximately 39% (73/187) and ~15% (28/187) of the isolates showed different genus identity and could not be identified using the phenotypic identification systems, respectively. Both methods succeeded in determining the species identities of 55 isolates; however, only ~69% (38/55) of the isolates matched at species level. 16S rDNA sequencing could not determine the species of ~20% (37/187) of the isolates. CONCLUSION: The 16S rDNA sequencing is a useful method over the phenotypic identification systems for the identification of rare and difficult to identify bacteria species. The 16S rDNA sequencing method, however, does have limitation for species-level identification of some bacteria highlighting the need for better bacterial pathogen identification tools.

While fleas are often perceived simply as a biting nuisance and a cause of allergic dermatitis, they represent important disease vectors worldwide, especially for bacterial zoonoses such as plague (transmitted by rodent fleas) and some of the rickettsioses and bartonelloses. The cosmopolitan cat ( Ctenocephalides felis ) and dog ( Ctenocephalides canis ) fleas, as well as Ctenocephalides orientis (restricted to tropical and subtropical Asia), breed in human dwellings and are vectors of cat-scratch fever (caused by Bartonella spp.) and Rickettsia spp., including Rickettsia felis (agent of flea-borne spotted fever) and Rickettsia asembonensis , a suspected pathogen. These Rickettsia spp. are members of a phylogenetic clade known as the ‘transitional group’, which includes both human pathogens and arthropod-specific endosymbionts. The relatively depauperate flea microbiome can also contain other endosymbionts, including a diverse range of Wolbachia strains. Here, we present circularized genome assemblies for two C. orientis -derived pathogens ( Bartonella clarridgeiae and R. asembonensis ) from Malaysia, a novel Wolbachia strain ( w Cori), and the C. orientis mitochondrion; all were obtained by direct metagenomic sequencing of flea tissues. Moreover, we isolated two Wolbachia strains from Malaysian C. felis into tick cell culture and recovered circularized genome assemblies for both, one of which ( w CfeF) is newly sequenced. We demonstrate that the three Wolbachia strains are representatives of different major clades (‘supergroups’), two of which appear to be flea-specific. These Wolbachia genomes exhibit unique combinations of features associated with reproductive parasitism or mutualism, including prophage WO, cytoplasmic incompatibility factors and the biotin operon of obligate intracellular microbes. The first circularized assembly for R. asembonensis includes a plasmid with a markedly different structure and gene content compared to the published plasmid; moreover, this novel plasmid was also detected in cat flea metagenomes from the USA. Analysis of loci under positive selection in the transitional group revealed genes involved in host–pathogen interactions that may facilitate host switching. Finally, the first B. clarridgeiae genome from Asia exhibited large-scale genome stability compared to isolates from other continents, except for SNPs in regions predicted to mediate interactions with the vertebrate host. These findings highlight the paucity of data on the genomic diversity of Ctenocephalides -associated bacteria and raise questions regarding how interactions between members of the flea microbiome might influence vector competence.

The recommended test guidelines for Zika virus (ZIKV) include using both molecular and serological tools. While the molecular tools are useful for detecting acute infection, the serological tools are useful for the detection of previous infections. Nevertheless, detection of ZIKV-specific antibodies remains a challenge due to the high cross-reactivity between ZIKV and other flaviviruses such as dengue virus (DENV) and Japanese encephalitis virus (JEV). The objective of this study is to evaluate the commercially available enzyme-linked immunosorbent assay (ELISA) for the detection of ZIKV IgG. In this study, we evaluated 6 commercially available anti-ZIKV IgG ELISA kits. Pre-characterized serum panels consisting of 70 sera were selected for the evaluation. The diagnostic accuracy of each ELISA kits was determined and compared to the gold standard, Foci Reduction Neutralization Test (FRNT). The present study established that the performance of the NS1-based anti-ZIKV IgG ELISA kit was superior to that which uses of the E protein as antigen. Overall, commercial ZIKV IgG ELISA showed varying test performances, with some achieving moderate to high test sensitivities and specificities. When compared against the FRNT, the test sensitivities ranged from 7.1% to 78.6%, whereas, the test specificities ranged from 40.0% to 100%. Limitation to the study includes the cross reactivity between flavivirus and specificity of the kit in addressing the cross reactivity.

OBJECTIVES: This study reported SARS-CoV-2 whole genome sequencing results from June 2021 to January 2022 from seven genome sequencing centers in Malaysia as part of the national surveillance program. METHODS: COVID-19 samples that tested positive by reverse transcription polymerase chain reaction and with cycle threshold values <30 were obtained throughout Malaysia. Sequencing of SARS-CoV-2 complete genomes was performed using Illumina, Oxford Nanopore, or Ion Torrent platforms. A total of 6163 SARS-CoV-2 complete genome sequences were generated over the surveillance period. All sequences were submitted to the Global Initiative on Sharing All Influenza Data database. RESULTS: From June 2021 to January 2022, Malaysia experienced the fourth wave of COVID-19 dominated by the Delta variant of concern, including the original B.1.617.2 lineage and descendant AY lineages. The B.1.617.2 lineage was identified as the early dominant circulating strain throughout the country but over time, was displaced by AY.59 and AY.79 lineages in Peninsular (west) Malaysia, and the AY.23 lineage in east Malaysia. In December 2021, pilgrims returning from Saudi Arabia facilitated the introduction and spread of the BA.1 lineage (Omicron variant of concern) in the country. CONCLUSION: The changing trends of circulating SARS-CoV-2 lineages were identified, with differences observed between west and east Malaysia. This initiative highlighted the importance of leveraging research expertise in the country to facilitate pandemic response and preparedness.

Despite the fact that coronavirus disease 2019 (COVID-19) treatment and management are now considerably regulated, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is still one of the leading causes of death in 2022. The availability of COVID-19 vaccines, FDA-approved antivirals, and monoclonal antibodies in low-income countries still poses an issue to be addressed. Natural products, particularly traditional Chinese medicines (TCMs) and medicinal plant extracts (or their active component), have challenged the dominance of drug repurposing and synthetic compound libraries in COVID-19 therapeutics. Their abundant resources and excellent antiviral performance make natural products a relatively cheap and readily available alternative for COVID-19 therapeutics. Here, we deliberately review the anti-SARS-CoV-2 mechanisms of the natural products, their potency (pharmacological profiles), and application strategies for COVID-19 intervention. In light of their advantages, this review is intended to acknowledge the potential of natural products as COVID-19 therapeutic candidates.

Resistance to pyrethroid insecticides is widespread in Indonesian Aedes aegypti (Linnaeus), the primary vector of dengue viruses. This study aims to investigate the mutations in the voltage-gated sodium channel (Vgsc) conferring pyrethroid resistance against Ae. aegypti populations from Indonesia. Molecular genotyping of mutations using polymerase chain reaction assay and direct DNA sequencing were performed at positions 989 and 1,016 in IIS6 region, and 1,534 in IIIS6 region of the voltage-gated sodium channel (Vgsc) in nine populations of Indonesian Ae. aegypti. The V1016G and S989P genotyping identified the RR genotype to be predominant in six out of nine populations of Ae. aegypti, whereas the SS genotype occurred only in minority. Interestingly, co-occurrence of the V1016G and S989P mutations was detected in the aforementioned six populations with high frequency. Genotyping of F1534C showed all nine populations exhibited the SS genotype, with merely two individuals from a population were heterozygous (RS). Significant correlations were demonstrated between the allele frequencies of the V1016G mutation and the survivability rates as well as resistance ratios in pyrethroid adult bioassays. This signifies the V1016G can contribute more to the insensitivity of Vgsc than the F1534C. Homozygous 1016G mosquitoes were likelier to survive pyrethroid exposure. Identification of underlying mechanisms resulting in insecticide resistance is advantageous in developing effective mosquito control programs in Indonesia.

BACKGROUND: Plasmodium knowlesi, a simian malaria parasite infection, increases as Plasmodium falciparum and Plasmodium vivax infections decrease in Johor, Malaysia. Therefore, this study aimed to identify the distribution of vectors involved in knowlesi malaria transmission in Johor. This finding is vital in estimating hotspot areas for targeted control strategies. METHODS: Anopheles mosquitoes were collected from the location where P. knowlesi cases were reported. Cases of knowlesi malaria from 2011 to 2019 in Johor were analyzed. Internal transcribed spacers 2 (ITS2) and cytochrome c oxidase subunit I (COI) genes were used to identify the Leucosphyrus Group of Anopheles mosquitoes. In addition, spatial analysis was carried out on the knowlesi cases and vectors in Johor. RESULTS: One hundred and eighty-nine cases of P. knowlesi were reported in Johor over 10 years. Young adults between the ages of 20-39 years comprised 65% of the cases. Most infected individuals were involved in agriculture and army-related occupations (22% and 32%, respectively). Four hundred and eighteen Leucosphyrus Group Anopheles mosquitoes were captured during the study. Anopheles introlatus was the predominant species, followed by Anopheles latens. Spatial analysis by Kriging interpolation found that hotspot regions of P. knowlesi overlapped or were close to the areas where An. introlatus and An. latens were found. A significantly high number of vectors and P. knowlesi cases were found near the road within 0-5 km. CONCLUSIONS: This study describes the distribution of P. knowlesi cases and Anopheles species in malaria-endemic transmission areas in Johor. Geospatial analysis is a valuable tool for studying the relationship between vectors and P. knowlesi cases. This study further supports that the Leucosphyrus Group of mosquitoes might be involved in transmitting knowlesi malaria cases in Johor. These findings may provide initial evidence to prioritize diseases and vector surveillance.

The COVID-19 pandemic, caused by the novel coronavirus SARS-CoV-2, has profoundly impacted global healthcare systems and spurred extensive research efforts over the past three years. One critical aspect of the disease is the intricate interplay between the virus and the host immune response, particularly the role of inflammatory gene expression in severe COVID-19. While numerous previous studies have explored the role of genetic polymorphisms in COVID-19, research specifically focusing on inflammatory genes and their associations with disease severity remains limited. This review explores the relationship between severe COVID-19 outcomes and genetic polymorphisms within key inflammatory genes. By investigating the impact of genetic variations on immune responses, which include cytokine production and downstream signalling pathways, we aim to provide a comprehensive overview of how genetic polymorphisms contribute to the variability in disease presentation. Through an in-depth analysis of existing literature, we shed light on potential therapeutic targets and personalized approaches that may enhance our understanding of disease pathogenesis and treatment strategies.

The Leucosphyrus Group of mosquitoes are the major simian malaria vectors in Malaysia. Accurate species identification is required to help in curbing the spread of simian malaria. The aim of the study is to provide an accurate molecular method for identifying the four important Anopheles vector species found in Malaysia. Mosquito specimens were collected from various localities in Malaysia, where simian malaria cases were reported. DNA from 122 mosquito specimens was tested to develop a multiplex polymerase chain reaction (PCR) assay. The specificity of this assay was tested against other mosquito species. Molecular identification of the species was further confirmed by analysing the internal transcribed spacer 2 (ITS2) DNA region of the specimens. Anopheles balabacensis and An. latens showed two distinct clades in the phylogenetic tree. The multiplex PCR assay was developed based on the ITS2 region for the identification of Anopheles introlatus (298–299 bp), Anopheles latens (197–198 bp), Anopheles cracens (421–426 bp), and Anopheles balabacensis (224–228 bp). This method will be useful to accurately identify the major Anopheles Leucosphyrus Group species in Malaysia, which are difficult to identify morphologically, to determine the correct vector as well as its geographical distribution.

The primary objective of this paper is to serve as a valuable resource for policymakers who are confronted with the evolving landscape of the coronavirus disease 2019 (COVID-19), considering both free and cost-based vaccination approaches. The potential consequences of shifting from free to cost-based vaccination are explored, encompassing its impact on global vaccine equity and prioritization, economic well-being, healthcare systems and delivery, public health policies, and vaccine distribution strategies. Examining past studies on willingness to pay for the initial COVID-19 vaccine dose and booster shots provides insights into how individuals value COVID-19 vaccinations and underscores the significance of addressing issues related to affordability. If COVID-19 vaccinations incur expenses, using effective communication strategies that emphasize the importance of vaccination and personal health benefits can increase willingness to pay. Making COVID-19 vaccines accessible through public health programs or health insurance can help alleviate financial barriers and increase vaccination rates.

Many human clinical cases attributed to vector-borne pathogens are underreported in Malaysia, especially in rural localities where healthcare infrastructures are lacking. Here, 217 small mammals, consisting of rodents and tree shrews, were trapped in oil palm plantations in the Peninsular Malaysia states of Johor and Perak. Species identification was performed using morphological and DNA barcoding analyses, and 203 small mammals were included in the detection of selected vector-borne bacteria. The DNA extracted from the spleens was examined for Orientia tsutsugamushi, Borrelia spp., Bartonella spp. and Rickettsia spp. using established PCR assays. The small mammals collected in this study included Rattus tanezumi R3 mitotype (n = 113), Rattus argentiventer (n = 24), Rattus tiomanicus (n = 22), Rattus exulans (n = 17), Rattus tanezumi sensu stricto (n = 1) and Tupaia glis (n = 40). Orientia tsutsugamushi, Borrelia spp. and Bartonella phoceensis were detected in the small mammals with the respective detection rates of 12.3%, 5.9% and 4.9%. Rickettsia spp., however, was not detected. This study encountered the presence of both Lyme disease and relapsing fever-related borreliae in small mammals collected from the oil palm plantation study sites. All three microorganisms (Orientia tsutsugamushi, Borrelia spp. and Bartonella phoceensis) were detected in the R. tanezumi R3 mitotype, suggesting that the species is a competent host for multiple microorganisms. Further investigations are warranted to elucidate the relationships between the ectoparasites, the small mammals and the respective pathogens.

Dengue fever (DF) is a mosquito-borne disease still with no effective treatment or vaccine available. A randomized, placebo-controlled, double-blinded, parallel-group trial was undertaken to evaluate the efficacy of oral intake of Lactococcus lactis strain plasma (LC-Plasma) on the presentation and severity of DF-like symptoms among healthy volunteers. Study participants (320) were assigned into two groups, and consumed either placebo or LC-Plasma tablets (approximately 100 billion cells/day) for 8 weeks. The clinical symptoms of DF were self-recorded through questionnaires, and exposure to DENV was determined by serum antibody and/or DENV antigen tests. No significant differences between groups were observed for exposure to DENV, or the symptomatic ratio. Results obtained showed that participants from the LC-Plasma group reported a significant reduction in the cumulative incidence days of DF-like symptoms, which include fever (p < 0.001), muscle pain (p < 0.005), joint pain (p < 0.001), and pain behind the eyes (p < 0.001), compared to that of the placebo group. Subgroup analysis revealed a significantly (p < 0.05) reduced severity score in the LC-Plasma group when study sites were separately analyzed. Overall, our findings suggest that LC-Plasma supplementation reduces the cumulative days with DF-like symptoms, and the severity of the symptoms. Daily oral intake of LC-Plasma, hence, is shown to mitigate the DF-like symptoms.

In this study, we present the draft genome sequence of B. pseudohinzii BH370 recovered from the trachea and lung tissues of an ICR mouse in Malaysia. The genome consists of 4,474,040 bp with a GC content of 66.4%. Annotation using RAST algorithm displayed 5119 protein encoding and 52 RNA genes. The CRISPR- cas genomic sequences previously reported in B. pseudohinzii were identified. The nucleotide sequences of BH370 was deposited into the European Nucleotide Archive under the genome assembly accession number FPJN01000000.

The COVID-19 pandemic has caused extensive loss of lives and economic hardship. In response, infectious disease experts and vaccine developers promptly responded by bringing forth candidate vaccines, some of which have been listed in the World Health Organization's Emergency Use Listing. Notwithstanding the diverse worldwide population genetics, the vaccines thus far developed are generic in nature for use worldwide. Differences in the human leukocyte antigen (HLA) in different populations, variation of the T cell epitopes, and the propensity of SARS-CoV-2 genetic mutations left room for improvement of the vaccines. Here, we discussed the implications of COVID-19 vaccination and SARS-CoV-2 infection by taking into consideration SARS-CoV-2 mutations, T cell epitopes, risk factors, and current platforms of candidate vaccines based on the HLA types that are commonly present in Peninsular Malaysia Chinese, Indian, and Malay populations. The HLA types associated with protection against and susceptibility to severe SARS-CoV-2 infection were identified based on reported case-control and cohort studies. The relevance of including the non-spike SARS-CoV-2 proteins in the future COVID-19 vaccines is also highlighted. This review is meant to trigger researchers to acknowledge the importance of investigating the possible relationships between the HLA haplotype and the SARS-CoV-2 strains circulating in different populations.

Despite the fact that COVID-19 treatment and management are now considerably regulated, SARS-CoV-2 is still one of the leading causes of death in 2022. The availability of COVID-19 vaccines, FDA-approved antivirals, and monoclonal antibodies in low-income countries still poses an issue to be addressed. Natural products, particularly traditional Chinese medicines (TCMs) and medicinal plant extracts (or their active component) have challenged the dominance of drug repurposing and synthetic compound libraries in COVID-19 therapeutics. Abundant resources and excellent antiviral performance make natural products a relatively cheap and easily available alternative for COVID-19 therapeutics. Here, we deliberately review the anti-SARS-CoV-2 mechanisms of the natural products, their potency (pharmacological profiles), and application strategies for COVID-19 intervention. In light of their advantages, this review is intended to acknowledge the potential of natural products as COVID-19 therapeutic candidates.

Aedes albopictus is one of the main mosquito vectors responsible for transmitting arboviruses to humans and animals. The ability of this mosquito to support virus transmission has been linked to vector competence, which is partly attributed to the genetic disparities in Ae. albopictus population. At present, little is known about the biologically important traits of Ae. albopictus in Malaysia. Thus, the study aims to determine the genetic variation of Ae. albopictus based on the mitochondria-encoded sequences of cytochrome oxidase subunit I (COI). A statistical parsimony network of 253 taxa aligned as 321 characters of the COI gene revealed 42 haplotypes (H1-H42), of which H1 was the most widespread haplotype in Peninsular Malaysia. Three highly divergent haplotypes (H21, H30, and H31) were detected from the northern population. Overall, haplotype and nucleotide diversities were 0.576 and 0.003, respectively, with low genetic differentiation (FST = 0.039) and high gene flow (Nm = 12.21) across all populations.