Eco-climate drivers shape virome diversity in a globally invasive tick species.

Spillovers of viruses into human occur more frequently under warmer conditions, particularly arboviruses. The invasive tick species Haemaphysalis longicornis poses a significant public health threat due to its global expansion and its potential to carry a wide range of pathogens. We analyzed meta-transcriptomic data from 3595 adult H. longicornis ticks collected between 2016 and 2019 in 22 provinces across China, encompassing diverse ecological conditions. Generalized additive modelling revealed that climate factors exerted a stronger influence on the virome of H. longicornis compared to other ecological factors, such as ecotypes, distance to coastline, animal host, tick gender, and anti-viral immunity. We investigated the mechanistic understanding of how climate changes drive the tick virome using causality inference and emphasized its significance for public health. Our findings demonstrated that higher temperatures and lower relative humidity/precipitation contribute to variations in animal host diversity, leading to an increased diversity of tick virome, particularly the evenness of vertebrate associated viruses. This finding may explain the evolution of tick-borne viruses into generalists across multiple hosts, thereby increasing the probability of spillover events involving tick-borne pathogens. Deep learning projections indicate that the diversity of H. longicornis virome is expected to increase in 81.9% of regions under the SSP8.5 scenario from 2019-2030. Extension of surveillance should be implemented to avert the spread of tick-borne diseases.

ARGNet: using deep neural networks for robust identification and classification of antibiotic resistance genes from sequences.

BACKGROUND: Emergence of antibiotic resistance in bacteria is an important threat to global health. Antibiotic resistance genes (ARGs) are some of the key components to define bacterial resistance and their spread in different environments. Identification of ARGs, particularly from high-throughput sequencing data of the specimens, is the state-of-the-art method for comprehensively monitoring their spread and evolution. Current computational methods to identify ARGs mainly rely on alignment-based sequence similarities with known ARGs. Such approaches are limited by choice of reference databases and may potentially miss novel ARGs. The similarity thresholds are usually simple and could not accommodate variations across different gene families and regions. It is also difficult to scale up when sequence data are increasing. RESULTS: In this study, we developed ARGNet, a deep neural network that incorporates an unsupervised learning autoencoder model to identify ARGs and a multiclass classification convolutional neural network to classify ARGs that do not depend on sequence alignment. This approach enables a more efficient discovery of both known and novel ARGs. ARGNet accepts both amino acid and nucleotide sequences of variable lengths, from partial (30-50 aa; 100-150 nt) sequences to full-length protein or genes, allowing its application in both target sequencing and metagenomic sequencing. Our performance evaluation showed that ARGNet outperformed other deep learning models including DeepARG and HMD-ARG in most of the application scenarios especially quasi-negative test and the analysis of prediction consistency with phylogenetic tree. ARGNet has a reduced inference runtime by up to 57% relative to DeepARG. CONCLUSIONS: ARGNet is flexible, efficient, and accurate at predicting a broad range of ARGs from the sequencing data. ARGNet is freely available at https://github.com/id-bioinfo/ARGNet , with an online service provided at https://ARGNet.hku.hk . Video Abstract.

Association between Gut Microbiota Composition and Long-Term Vaccine Immunogenicity following Three Doses of CoronaVac.

BACKGROUND: Neutralizing antibody level wanes with time after COVID-19 vaccination. We aimed to study the relationship between baseline gut microbiota and immunogenicity after three doses of CoronaVac. METHODS: This was a prospective cohort study recruiting three-dose CoronaVac recipients from two centers in Hong Kong. Blood samples were collected at baseline and one year post-first dose for virus microneutralization (vMN) assays to determine neutralization titers. The primary outcome was high immune response (defined as with vMN titer ≥ 40). Shotgun DNA metagenomic sequencing of baseline fecal samples identified potential bacterial species and metabolic pathways using Linear Discriminant Analysis Effect Size (LEfSe) analysis. Univariate and multivariable logistic regression models were used to identify high response predictors. RESULTS: In total, 36 subjects were recruited (median age: 52.7 years [IQR: 47.9-56.4]; male: 14 [38.9%]), and 18 had low immune response at one year post-first dose vaccination. Eubacterium rectale (log10LDA score = 4.15, p = 0.001; relative abundance of 1.4% vs. 0, p = 0.002), Collinsella aerofaciens (log10LDA score = 3.31, p = 0.037; 0.39% vs. 0.18%, p = 0.038), and Streptococcus salivarius (log10LDA score = 2.79, p = 0.021; 0.05% vs. 0.02%, p = 0.022) were enriched in low responders. The aOR of high immune response with E. rectale, C. aerofaciens, and S. salivarius was 0.03 (95% CI: 9.56 × 10-4-0.32), 0.03 (95% CI: 4.47 × 10-4-0.59), and 10.19 (95% CI: 0.81-323.88), respectively. S. salivarius had a positive correlation with pathways enriched in high responders like incomplete reductive TCA cycle (log10LDA score = 2.23). C. aerofaciens similarly correlated with amino acid biosynthesis-related pathways. These pathways all showed anti-inflammation functions. CONCLUSION: E. rectale,C. aerofaciens, and S. salivarius correlated with poorer long-term immunogenicity following three doses of CoronaVac.

Robust expansion of phylogeny for fast-growing genome sequence data.

Massive sequencing of SARS-CoV-2 genomes has urged novel methods that employ existing phylogenies to add new samples efficiently instead of de novo inference. 'TIPars' was developed for such challenge integrating parsimony analysis with pre-computed ancestral sequences. It took about 21 seconds to insert 100 SARS-CoV-2 genomes into a 100k-taxa reference tree using 1.4 gigabytes. Benchmarking on four datasets, TIPars achieved the highest accuracy for phylogenies of moderately similar sequences. For highly similar and divergent scenarios, fully parsimony-based and likelihood-based phylogenetic placement methods performed the best respectively while TIPars was the second best. TIPars accomplished efficient and accurate expansion of phylogenies of both similar and divergent sequences, which would have broad biological applications beyond SARS-CoV-2. TIPars is accessible from https://tipars.hku.hk/ and source codes are available at https://github.com/id-bioinfo/TIPars.

Binding affinity between coronavirus spike protein and human ACE2 receptor.

Coronaviruses (CoVs) pose a major risk to global public health due to their ability to infect diverse animal species and potential for emergence in humans. The CoV spike protein mediates viral entry into the cell and plays a crucial role in determining the binding affinity to host cell receptors. With particular emphasis on α- and ß-coronaviruses that infect humans and domestic animals, current research on CoV receptor use suggests that the exploitation of the angiotensin-converting enzyme 2 (ACE2) receptor poses a significant threat for viral emergence with pandemic potential. This review summarizes the approaches used to study binding interactions between CoV spike proteins and the human ACE2 (hACE2) receptor. Solid-phase enzyme immunoassays and cell binding assays allow qualitative assessment of binding but lack quantitative evaluation of affinity. Surface plasmon resonance, Bio-layer interferometry, and Microscale Thermophoresis on the other hand, provide accurate affinity measurement through equilibrium dissociation constants (KD). In silico modeling predicts affinity through binding structure modeling, protein-protein docking simulations, and binding energy calculations but reveals inconsistent results due to the lack of a standardized approach. Machine learning and deep learning models utilize simulated and experimental protein-protein interaction data to elucidate the critical residues associated with CoV binding affinity to hACE2. Further optimization and standardization of existing approaches for studying binding affinity could aid pandemic preparedness. Specifically, prioritizing surveillance of CoVs that can bind to human receptors stands to mitigate the risk of zoonotic spillover.

Characterization of an Emergent Chicken H3N8 Influenza Virus in Southern China: a Potential Threat to Public Health.

Although influenza A viruses of several subtypes have occasionally infected humans, to date only those of the H1, H2, and H3 subtypes have led to pandemics and become established in humans. The detection of two human infections by avian H3N8 viruses in April and May of 2022 raised pandemic concerns. Recent studies have shown the H3N8 viruses were introduced into humans from poultry, although their genesis, prevalence, and transmissibility in mammals have not been fully elucidated. Findings generated from our systematic influenza surveillance showed that this H3N8 influenza virus was first detected in chickens in July 2021 and then disseminated and became established in chickens over wider regions of China. Phylogenetic analyses revealed that the H3 HA and N8 NA were derived from avian viruses prevalent in domestic ducks in the Guangxi-Guangdong region, while all internal genes were from enzootic poultry H9N2 viruses. The novel H3N8 viruses form independent lineages in the glycoprotein gene trees, but their internal genes are mixed with those of H9N2 viruses, indicating continuous gene exchange among these viruses. Experimental infection of ferrets with three chicken H3N8 viruses showed transmission through direct contact and inefficient transmission by airborne exposure. Examination of contemporary human sera detected only very limited antibody cross-reaction to these viruses. The continuing evolution of these viruses in poultry could pose an ongoing pandemic threat. IMPORTANCE A novel H3N8 virus with demonstrated zoonotic potential has emerged and disseminated in chickens in China. It was generated by reassortment between avian H3 and N8 virus(es) and long-term enzootic H9N2 viruses present in southern China. This H3N8 virus has maintained independent H3 and N8 gene lineages but continues to exchange internal genes with other H9N2 viruses to form novel variants. Our experimental studies showed that these H3N8 viruses were transmissible in ferrets, and serological data suggest that the human population lacks effective immunological protection against it. With its wide geographical distribution and continuing evolution in chickens, other spillovers to humans can be expected and might lead to more efficient transmission in humans.

Probable airborne transmission of Burkholderia pseudomallei causing an urban outbreak of melioidosis during typhoon season in Hong Kong, China.

Between January 2015 and October 2022, 38 patients with culture-confirmed melioidosis were identified in the Kowloon West (KW) Region, Hong Kong. Notably, 30 of them were clustered in the Sham Shui Po (SSP) district, which covers an estimated area of 2.5 km2. Between August and October 2022, 18 patients were identified in this district after heavy rainfall and typhoons. The sudden upsurge in cases prompted an environmental investigation, which involved collecting 20 air samples and 72 soil samples from residential areas near the patients. A viable isolate of Burkholderia pseudomallei was obtained from an air sample collected at a building site five days after a typhoon. B. pseudomallei DNA was also detected in 21 soil samples collected from the building site and adjacent gardening areas using full-length 16S rRNA gene sequencing, suggesting that B. psuedomallei is widely distributed in the soil environment surrounding the district. Core genome-multilocus sequence typing showed that the air sample isolate was phylogenetically clustered with the outbreak isolates in KW Region. Multispectral satellite imagery revealed a continuous reduction in vegetation region in SSP district by 162,255 m2 from 2016 to 2022, supporting the hypothesis of inhalation of aerosols from the contaminated soil as the transmission route of melioidosis during extreme weather events. This is because the bacteria in unvegetated soil are more easily spread by winds. In consistent with inhalational melioidosis, 24 (63.2%) patients had pneumonia. Clinicians should be aware of melioidosis during typhoon season and initiate appropriate investigation and treatment for patients with compatible symptoms.

Skin infectome of patients with a tick bite history.

Introduction: Ticks are the most important obligate blood-feeding vectors of human pathogens. With the advance of high-throughput sequencing, more and more bacterial community and virome in tick has been reported, which seems to pose a great threat to people. Methods: A total of 14 skin specimens collected from tick-bite patients with mild to severe symptoms were analyzed through meta-transcriptomic sequencings. Results: Four bacteria genera were both detected in the skins and ticks, including Pseudomonas, Acinetobacter, Corynebacterium and Propionibacterium, and three tick-associated viruses, Jingmen tick virus (JMTV), Bole tick virus 4 (BLTV4) and Deer tick mononegavirales-like virus (DTMV) were identified in the skin samples. Except of known pathogens such as pathogenic rickettsia, Coxiella burnetii and JMTV, we suggest Roseomonas cervicalis and BLTV4 as potential new agents amplified in the skins and then disseminated into the blood. As early as 1 day after a tick-bite, these pathogens can transmit to skins and at most four ones can co-infect in skins. Discussion: Advances in sequencing technologies have revealed that the diversity of tick microbiome and virome goes far beyond our previous understanding. This report not only identifies three new potential pathogens in humans but also shows that the skin barrier is vital in preventing horizontal transmissions of tick-associated bacteria or virus communities to the host. It is the first research on patients' skin infectome after a tick bite and demonstrates that more attention should be paid to the cutaneous response to prevent tick-borne illness.

Global environmental resistome: Distinction and connectivity across diverse habitats benchmarked by metagenomic analyses.

The widely distributed antibiotic resistance genes (ARGs) were unevenly proliferated in various habitats. Great endeavors are needed to resolve the resistome features that can differentiate or connect different habitats. This study retrieved a broad spectrum of resistome profiles from 1723 metagenomes categorized into 13 habitats, encompassing industrial, urban, agricultural, and natural environments, and spanning most continents and oceans. The resistome features (ARG types, subtypes, indicator ARGs, and emerging mobilizable ARGs: mcr and tet(X)) in these habitats were benchmarked via a standardized workflow. We found that wastewater and wastewater treatment works were characterized to be reservoirs of more diverse genotypes of ARGs than any other habitats including human and livestock fecal samples, while fecal samples were with higher ARG abundance. Bacterial taxonomy composition was significantly correlated with resistome composition across most habitats. Moreover, the source-sink connectivities were disentangled by developing the resistome-based microbial attribution prediction model. Environmental surveys with standardized bioinformatic workflow proposed in this study will help comprehensively understand the transfer of ARGs in the environment, thus prioritizing the critical environments with high risks for intervention to tackle the problem of ARGs.

Evaluation of Mycobacterium tuberculosis enrichment in metagenomic samples using ONT adaptive sequencing and amplicon sequencing for identification and variant calling.

Sensitive detection of Mycobacterium tuberculosis (TB) in small percentages in metagenomic samples is essential for microbial classification and drug resistance prediction. However, traditional methods, such as bacterial culture and microscopy, are time-consuming and sometimes have limited TB detection sensitivity. Oxford nanopore technologies (ONT) MinION sequencing allows rapid and simple sample preparation for sequencing. Its recently developed adaptive sequencing selects reads from targets while allowing real-time base-calling to achieve sequence enrichment or depletion during sequencing. Another common enrichment method is PCR amplification of the target TB genes. In this study, we compared both methods using ONT MinION sequencing for TB detection and variant calling in metagenomic samples using both simulation runs and those with synthetic and patient samples. We found that both methods effectively enrich TB reads from a high percentage of human (95%) and other microbial DNA. Adaptive sequencing with readfish and UNCALLDE achieved a 3.9-fold and 2.2-fold enrichment compared to the control run. We provide a simple automatic analysis framework to support the detection of TB for clinical use, openly available at https://github.com/HKU-BAL/ONT-TB-NF . Depending on the patient's medical condition and sample type, we recommend users evaluate and optimize their workflow for different clinical specimens to improve the detection limit.

New insights into the impact of microbiome on horizontal and vertical transmission of a tick-borne pathogen.

BACKGROUND: The impact of host skin microbiome on horizontal transmission of tick-borne pathogens , and of pathogen associated transstadial and transovarial changes in tick microbiome are largely unknown, but are important to control increasingly emerging tick-borne diseases worldwide. METHODS: Focusing on a rickettsiosis pathogen, Rickettsia raoultii, we used R. raoultii-positive and R. raoultii-negative Dermacentor spp. tick colonies to study the involvement of skin microbiota in cutaneous infection with rickettsiae in laboratory mice, and the function of the tick microbiome on maintenance of rickettsiae through all tick developmental stages (eggs, larvae, nymphs, adults) over two generations. RESULTS: We observed changes in the skin bacteria community, such as Chlamydia, not only associated with rickettsial colonization but also with tick feeding on skin. The diversity of skin microbiome differed between paired tick-bitten and un-bitten sites. For vertical transmission, significant differences in the tick microbiota between pathogenic rickettsia-positive and -negative tick chorts was observed across all developmental stages at least over two generations, which appeared to be a common pattern not only for R. raoultii but also for another pathogenic species, Candidatus Rickettsia tarasevichiae. More importantly, bacterial differences were complemented by functional shifts primed for genetic information processing during blood feeding. Specifically, the differences in tick microbiome gene repertoire between pathogenic Rickettsia-positive and -negative progenies were enriched in pathways associated with metabolism and hormone signals during vertical transmission. CONCLUSIONS: We demonstrate that host skin microbiome might be a new factor determining the transmission of rickettsial pathogens through ticks. While pathogenic rickettsiae infect vertebrate hosts during blood-feeding by the tick, they may also manipulate the maturation of the tick through changing the functional potential of its microbiota over the tick's life stages. The findings here might spur the development of new-generation control methods for ticks and tick-borne pathogens. Video Abstract.

Metavirome of 31 tick species provides a compendium of 1,801 RNA virus genomes.

The increasing prevalence and expanding distribution of tick-borne viruses globally have raised health concerns, but the full repertoire of the tick virome has not been assessed. We sequenced the meta-transcriptomes of 31 different tick species in the Ixodidae and Argasidae families from across mainland China, and identified 724 RNA viruses with distinctive virome compositions among genera. A total of 1,801 assembled and complete or nearly complete viral genomes revealed an extensive diversity of genome architectures of tick-associated viruses, highlighting ticks as a reservoir of RNA viruses. We examined the phylogenies of different virus families to investigate virome evolution and found that the most diverse tick-associated viruses are positive-strand RNA virus families that demonstrate more ancient divergence than other arboviruses. Tick-specific viruses are often associated with only a few tick species, whereas virus clades that can infect vertebrates are found in a wider range of tick species. We hypothesize that tick viruses can exhibit both 'specialist' and 'generalist' evolutionary trends. We hope that our virome dataset will enable much-needed research on vertebrate-pathogenic tick-associated viruses.

Genotype and organ effect on the occupancy of phyllosphere prokaryotes in different rice landraces.

Bacteria community provides essential ecological services to rice plants. The bacterial diversity of rice varies across host plant genotype and organs. This study employed 16S rDNA amplicon sequencing to characterise the bacterial community associated with three rice landraces using leaf blade and stem samples. The prokaryotic community found in these rice landraces comprised of two kingdoms, 12 phyla, 25 classes, 40 orders, 80 families, and 118 genera. Proteobacteria (53.9%) was the most abundant phylum. The most abundant genus was an undefined genus under Cyanobacteria (33.0%). Homogeneity of prokaryotic community was observed across the three rice landraces, which may suggest a high similarity in biological and genetical properties of the rice landraces. The difference in prokaryotic composition between leaf blade and stem was depicted based on principal coordinate analysis. This study observed that the prokaryotic inhabitants in rice plants is predominantly determined by rice plant organs.

Association between inflammatory cytokines and anti-SARS-CoV-2 antibodies in hospitalized patients with COVID-19.

BACKGROUND: COVID-19 patients may experience "cytokine storm" when human immune system produces excessive cytokines/chemokines. However, it remains unclear whether early responses of inflammatory cytokines would lead to high or low titers of anti-SARS-CoV-2 antibodies. METHODS: This retrospective study enrolled a cohort of 272 hospitalized patients with laboratory-confirmed SARS-CoV-2. Laboratory assessments of serum cytokines (IL-2R, IL-6, IL-8, IL-10, TNF-α), anti-SARS-CoV-2 IgG/IgM antibodies, and peripheral blood biomarkers were conducted during hospitalization. RESULTS: At hospital admission, 36.4% patients were severely ill, 51.5% patients were ≥ 65 years, and 60.3% patients had comorbidities. Higher levels of IL-2R and IL-6 were observed in older patients (≥65 years). Significant differences of IL-2R (week 2 to week ≥5 from symptom onset), IL-6 (week 1 to week ≥5), IL-8 (week 2 to week ≥5), and IL-10 (week 1 to week 3) were observed between moderately-ill and severely ill patients. Anti-SARS-CoV-2 IgG titers were significantly higher in severely ill patients than in moderately ill patients, but such difference was not observed for IgM. High titers of early-stage IL-6, IL-8, and TNF-α (≤2 weeks after symptom onset) were positively correlated with high titers of late-stage IgG (≥5 weeks after symptom onset). Deaths were mostly observed in severely ill older patients (45.9%). Survival analyses revealed risk factors of patient age, baseline COVID-19 severity, and baseline IL-6 that affected survival time, especially in severely ill older patients. CONCLUSION: Early responses of elevated cytokines such as IL-6 reflect the active immune responses, leading to high titers of IgG antibodies against COVID-19.

Unlocking capacities of viral genomics for the COVID-19 pandemic response.

More than any other infectious disease epidemic, the COVID-19 pandemic has been characterized by the generation of large volumes of viral genomic data at an incredible pace due to recent advances in high-throughput sequencing technologies, the rapid global spread of SARS-CoV-2, and its persistent threat to public health. However, distinguishing the most epidemiologically relevant information encoded in these vast amounts of data requires substantial effort across the research and public health communities. Studies of SARS-CoV-2 genomes have been critical in tracking the spread of variants and understanding its epidemic dynamics, and may prove crucial for controlling future epidemics and alleviating significant public health burdens. Together, genomic data and bioinformatics methods enable broad-scale investigations of the spread of SARS-CoV-2 at the local, national, and global scales and allow researchers the ability to efficiently track the emergence of novel variants, reconstruct epidemic dynamics, and provide important insights into drug and vaccine development and disease control. Here, we discuss the tremendous opportunities that genomics offers to unlock the effective use of SARS-CoV-2 genomic data for efficient public health surveillance and guiding timely responses to COVID-19.

Genomic characterization of a new CRF01_AE/CRF07_BC case from a MSM patient in Guangdong, China.

The CRF01_AE and CRF07_BC clades dominate the human immunodeficiency virus (HIV) epidemics in China. Both clades have been identified in the men who have sex with men (MSM) population in Guangdong province, raising a serious concern of possible complex recombination events ahead. Here, we report the first case of CRF01_AE/CRF07_BC recombinant sampled from a MSM patient in southern China. The genomic structure of this case is a mosaic with some regions resembling the CRF01_AE and CRF07_BC clades. Our phylogenetic analyses show that the two parental lineages of this recombinant virus were mainly found in the MSM population. This case has a different genomic composition compared with other recombinants descended from the same parental clades CRF01_AE and CRF07_BC. Our finding suggests that the MSM populations have become a hotspot for expanding viral diversity through the viral recombination mechanism. Therefore, further epidemiologic surveillance and monitoring should be conducted within the MSM populations to help advance our knowledge of viral transmission mechanisms. Additionally, these measures will serve to enhance the control and prevention of HIV/acquired immunodeficiency syndrome in China.

Genetic Variation of Multiple Serotypes of Enteroviruses Associated with Hand, Foot and Mouth Disease in Southern China.

Enteroviruses (EVs) species A are a major public health issue in the Asia-Pacific region and cause frequent epidemics of hand, foot and mouth disease (HFMD) in China. Mild infections are common in children; however, HFMD can also cause severe illness that affects the central nervous system. To molecularly characterize EVs, a prospective HFMD virological surveillance program was performed in China between 2013 and 2016. Throat swabs, rectal swabs and stool samples were collected from suspected HFMD patients at participating hospitals. EVs were detected using generic real-time and nested reverse transcription-polymerase chain reactions (RT-PCRs). Then, the complete VP1 regions of enterovirus A71 (EV-A71), coxsackievirus A16 (CVA16) and CVA6 were sequenced to analyze amino acid changes and construct a viral molecular phylogeny. Of the 2836 enrolled HFMD patients, 2,517 (89%) were EV positive. The most frequently detected EVs were CVA16 (32.5%, 819), CVA6 (31.2%, 785), and EV-A71 (20.4%, 514). The subgenogroups CVA16_B1b, CVA6_D3a and EV-A71_C4a were predominant in China and recombination was not observed in the VP1 region. Sequence analysis revealed amino acid variations at the 30, 29 and 44 positions in the VP1 region of EV-A71, CVA16 and CVA6 (compared to the respective prototype strains BrCr, G10 and Gdula), respectively. Furthermore, in 21 of 24 (87.5%) identified EV-A71 samples, a known amino acid substitution (D31N) that may enhance neurovirulence was detected. Our study provides insights about the genetic characteristics of common HFMD-associated EVs. However, the emergence and virulence of the described mutations require further investigation.

Sequencing and characterisation of complete mitogenome DNA for Rasbora sarawakensis (Cypriniformes: Cyprinidae: Rasbora) with phylogenetic consideration.

The Blueline Rasbora (Rasbora sarawakensis) is a small ray-finned fish categorized under the genus Rasbora in the Cyprinidae family. In this study, the complete mitogenome sequence of R. sarawakensis was sequenced using four primers targeting overlapping regions. The mitogenome is 16,709 bp in size, accommodating 22 transfer RNA genes, 13 protein-coding genes, two ribosomal RNA genes and a putative control region. Identical gene organisation was detected between this species and other genus counterparts. The heavy strand houses 28 genes while the light strand stores the other nine genes. Most protein-coding genes employ ATG as start codon, excluding COI gene, which utilizes GTG instead. The central conserved sequence blocks (CSB-F, CSB-E and CSB-D), variable sequence blocks (CSB-3, CSB-2 and CSB-1) as well as the terminal associated sequence (TAS) are conserved in the control region. The maximum likelihood phylogenetic tree revealed the divergence of R. sarawakensis from the basal region of the Rasbora clade, where its evolutionary relationships with R. maculatus and R. pauciperforata are poorly resolved as indicated by the low bootstrap values. This work acts as steppingstone towards further molecular evolution and population genetics studies of Rasbora genus in future.