Surveillance, isolation and genomic characterization of Pteropine orthoreovirus of probable bat origin among patients with acute respiratory infection in Malaysia.

Pteropine orthoreovirus (PRV), an emerging bat-borne virus, has been linked to cases of acute respiratory infections (ARI) in humans. The prevalence, epidemiology and genomic diversity of PRV among ARI of unknown origin were studied. Among 632 urban outpatients tested negative for all known respiratory viruses, 2.2% were PRV-positive. Patients mainly presented with moderate to severe forms of cough, sore throat and muscle ache, but rarely with fever. Phylogenetic analysis revealed that over 90% of patients infected with the Melaka virus (MelV)-like PRV, while one patient infected with the Pulau virus previously found only in fruit bats. Human oral keratinocytes and nasopharyngeal epithelial cells were susceptible to clinical isolates of PRV, including the newly isolated MelV-like 12MYKLU1034. Whole genome sequence of 12MYKLU1034 using Nanopore technique revealed a novel reassortant strain. Evolutionary analysis of the global PRV strains suggests the continuous evolution of PRV through genetic reassortment among PRV strains circulating in human, bats and non-human primate hosts, creating a spectrum of reassortant lineages with complex evolutionary characteristics. In summary, the role of PRV as a common etiologic agent of ARI is evident. Continuous monitoring of PRV prevalence, pathogenicity and diversity among human and animal hosts is important to trace the emergence of novel reassortants.

The IDentif.AI-x pandemic readiness platform: Rapid prioritization of optimized COVID-19 combination therapy regimens.

IDentif.AI-x, a clinically actionable artificial intelligence platform, was used to rapidly pinpoint and prioritize optimal combination therapies against COVID-19 by pairing a prospective, experimental validation of multi-drug efficacy on a SARS-CoV-2 live virus and Vero E6 assay with a quadratic optimization workflow. A starting pool of 12 candidate drugs developed in collaboration with a community of infectious disease clinicians was first narrowed down to a six-drug pool and then interrogated in 50 combination regimens at three dosing levels per drug, representing 729 possible combinations. IDentif.AI-x revealed EIDD-1931 to be a strong candidate upon which multiple drug combinations can be derived, and pinpointed a number of clinically actionable drug interactions, which were further reconfirmed in SARS-CoV-2 variants B.1.351 (Beta) and B.1.617.2 (Delta). IDentif.AI-x prioritized promising drug combinations for clinical translation and can be immediately adjusted and re-executed with a new pool of promising therapies in an actionable path towards rapidly optimizing combination therapy following pandemic emergence.

Application of a VP4/VP2-inferred transmission clusters in estimating the impact of interventions on rhinovirus transmission.

BACKGROUND: Despite the clinical burden attributable to rhinovirus (RV) infections, the RV transmission dynamics and the impact of interventions on viral transmission remain elusive. METHODS: A total of 3,935 nasopharyngeal specimens were examined, from which the VP4/VP2 gene was sequenced and genotyped. RV transmission clusters were reconstructed using the genetic threshold of 0.005 substitutions/site, estimated from the global VP4/VP2 sequences. A transmission cluster is characterized by the presence of at least two individuals (represent by nodes), whose viral sequences are genetically linked (represent by undirected edges) at the estimated genetic distance threshold supported by bootstrap value of ≥ 90%. To assess the impact of facemask, pleconaril and social distancing on RV transmission clusters, trials were simulated for interventions with varying efficacy and were evaluated based on the reduction in the number of infected patients (nodes) and the reduction in the number of nodes-connecting edges. The putative impact of intervention strategies on RV transmission clusters was evaluated through 10,000 simulations. RESULTS: A substantial clustering of 168 RV transmission clusters of varying sizes were observed. This suggests that RV disease burden observed in the population was largely due to multiple sub-epidemics, predominantly driven by RV-A, followed by RV-C and -B. No misclassification of RV species and types were observed, suggesting the specificity and sensitivity of the analysis. Through 10,000 simulations, it was shown that social distancing may be effective in decelerating RV transmission, by removing more than 95% of nodes and edges within the RV transmission clusters. However, facemask removed less than 8% and 66% of nodes and edges, respectively, conferring moderate advantage in limiting RV transmission. CONCLUSION: Here, we presented a network-based approach of which the degree of RV spread that fuel disease transmission in the region was mapped for the first time. The utilization of RV transmission clusters in assessing the putative impact of interventions on disease transmission at the population level was demonstrated.

Polymersomes as Stable Nanocarriers for a Highly Immunogenic and Durable SARS-CoV-2 Spike Protein Subunit Vaccine.

Multiple successful vaccines against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) are urgently needed to address the ongoing coronavirus disease 2019 (Covid-19) pandemic. In the present work, we describe a subunit vaccine based on the SARS-CoV-2 spike protein coadministered with CpG adjuvant. To enhance the immunogenicity of our formulation, both antigen and adjuvant were encapsulated with our proprietary artificial cell membrane (ACM) polymersome technology. Structurally, ACM polymersomes are self-assembling nanoscale vesicles made up of an amphiphilic block copolymer comprising poly(butadiene)-b-poly(ethylene glycol) and a cationic lipid, 1,2-dioleoyl-3-trimethylammonium-propane. Functionally, ACM polymersomes serve as delivery vehicles that are efficiently taken up by dendritic cells (DC1 and DC2), which are key initiators of the adaptive immune response. Two doses of our formulation elicit robust neutralizing antibody titers in C57BL/6 mice that persist at least 40 days. Furthermore, we confirm the presence of functional memory CD4+ and CD8+ T cells that produce T helper type 1 cytokines. This study is an important step toward the development of an efficacious vaccine in humans.

Genetic diversity of Coxsackievirus A21 associated with sporadic cases of acute respiratory infections in Malaysia.

BACKGROUND: Coxsackievirus A21 (CVA21), a member of Enterovirus C from the Picornaviridae family, has been associated with respiratory illnesses in humans. METHODS: A molecular epidemiological investigation of CVA21 was conducted among patients presenting with acute upper respiratory illnesses in the ambulatory settings between 2012 and 2014 in Kuala Lumpur, Malaysia. RESULTS: Epidemiological surveillance of acute respiratory infections (n = 3935) showed low-level detection of CVA21 (0.08%, 1.4 cases/year) in Kuala Lumpur, with no clear seasonal distribution. Phylogenetic analysis of the new complete genomes showed close relationship with CVA21 strains from China and the United States. Spatio-temporal mapping of the VP1 gene determined 2 major clusters circulating worldwide, with inter-country lineage migration and strain replacement occurring over time. CONCLUSIONS: The study highlights the emerging role of CVA21 in causing sporadic acute respiratory outbreaks.

Spike S2 Subunit: The Dark Horse in the Race for Prophylactic and Therapeutic Interventions against SARS-CoV-2.

In the midst of the unceasing COVID-19 pandemic, the identification of immunogenic epitopes in the SARS-CoV-2 spike (S) glycoprotein plays a vital role in the advancement and development of intervention strategies. S is expressed on the exterior of the SARS-CoV-2 virion and contains two subunits, namely the N-terminal S1 and C-terminal S2. It is the key element for mediating viral entry as well as a crucial antigenic determinant capable of stimulating protective immune response through elicitation of anti-SARS-CoV-2 antibodies and activation of CD4+ and CD8+ cells in COVID-19 patients. Given that S2 is highly conserved in comparison to the S1, here, we provide a review of the latest findings on the SARS-CoV-2 S2 subunit and further discuss its potential as an attractive and promising target for the development of prophylactic vaccines and therapeutic agents against COVID-19.

Genome sequencing and phylogenetic reconstruction reveal a potential fourth rhinovirus species and its worldwide distribution.

Genome sequences of members of a potential fourth rhinovirus (RV) species, provisionally denoted as rhinovirus A clade D, from patients with acute respiratory infection were determined. Bayesian coalescent analysis estimated that clade D emerged around the 1940s and diverged further around 2006-2007 into two distinctive sublineages (RV-A8-like and RV-A45-like) that harbored unique "clade-defining" substitutions. Similarity plots and bootscan mapping revealed a recombination breakpoint located in the 5'-UTR region of members of the RV-A8-like sublineage. Phylogenetic reconstruction revealed the distribution of clade D viruses in the Asia Pacific region and in Europe, underlining its worldwide distribution.

The effects of age on clinical characteristics, hospitalization and mortality of patients with influenza-related illness at a tertiary care centre in Malaysia.

BACKGROUND: Age is an established risk factor for poor outcomes in individuals with influenza-related illness, and data on its influence on clinical presentations and outcomes in the South-East Asian settings are scarce. The aim of this study was to determine the above among adults with influenza-related upper respiratory tract infection at a teaching hospital in Malaysia. METHODS: A retrospective case-note analysis was conducted on a cohort of 3935 patients attending primary care at the University Malaya Medical Centre, Malaysia from February 2012 till May 2014 with URTI symptoms. Demographics, clinical characteristics, medical and vaccination history were obtained from electronic medical records, and mortality data from the National Registration Department. Comparisons were made between those aged <25, ≥25 to <65 and ≥65 years. RESULTS: 470 (11.9%) had PCR-confirmed influenza virus infection. Six (1.3%) received prior influenza vaccination. Those aged ≥65 years were more likely to have ≥2 comorbidities (P < .001) and were less likely to present with fever (P = .004). One-third of those aged ≥65 years experienced hospitalization, intensive care admission or death within a year compared to 10% in the ≥25 to <65 years. Age ≥65 years was an independent predictor of hospitalization and death (OR = 9.97; 95% CI = 3.11-31.93) compared to those aged <25 years. CONCLUSION: Older patients in our cohort were more likely to have comorbidities and present with atypical features, with older age being an independent predictor of poor health outcomes. Our findings will now inform future health policies on older persons and economic modelling of adult vaccination programmes.

The role of human Metapneumovirus genetic diversity and nasopharyngeal viral load on symptom severity in adults.

BACKGROUND: Human metapneumovirus (HMPV) is established as one of the causative agents of respiratory tract infections. To date, there are limited reports that describe the effect of HMPV genotypes and/or viral load on disease pathogenesis in adults. This study aims to determine the role of HMPV genetic diversity and nasopharyngeal viral load on symptom severity in outpatient adults with acute respiratory tract infections. METHODS: Severity of common cold symptoms of patients from a teaching hospital was assessed by a four-category scale and summed to obtain the total symptom severity score (TSSS). Association between the fusion and glycoprotein genes diversity, viral load (quantified using an improved RT-qPCR assay), and symptom severity were analyzed using bivariate and linear regression analyses. RESULTS: Among 81/3706 HMPV-positive patients, there were no significant differences in terms of demographics, number of days elapsed between symptom onset and clinic visit, respiratory symptoms manifestation and severity between different HMPV genotypes/sub-lineages. Surprisingly, elderly patients (≥65 years old) had lower severity of symptoms (indicated by TSSS) than young and middle age adults (p = 0.008). Nasopharyngeal viral load did not correlate with nor predict symptom severity of HMPV infection. Interestingly, at 3-5 days after symptom onset, genotype A-infected patients had higher viral load compared to genotype B (4.4 vs. 3.3 log10 RNA copies/µl) (p = 0.003). CONCLUSIONS: Overall, HMPV genetic diversity and viral load did not impact symptom severity in adults with acute respiratory tract infections. Differences in viral load dynamics over time between genotypes may have important implications on viral transmission.

Viral Load and Sequence Analysis Reveal the Symptom Severity, Diversity, and Transmission Clusters of Rhinovirus Infections.

Background: Rhinovirus (RV) is one of the main viral etiologic agents of acute respiratory illnesses. Despite the heightened disease burden caused by RV, the viral factors that increase the severity of RV infection, the transmission pattern, and seasonality of RV infections remain unclear. Methods: An observational study was conducted among 3935 patients presenting with acute upper respiratory illnesses in the ambulatory settings between 2012 and 2014. Results: The VP4/VP2 gene was genotyped from all 976 RV-positive specimens, where the predominance of RV-A (49%) was observed, followed by RV-C (38%) and RV-B (13%). A significant regression in median nasopharyngeal viral load (VL) (P < .001) was observed, from 883 viral copies/µL at 1-2 days after symptom onset to 312 viral copies/µL at 3-4 days and 158 viral copies/µL at 5-7 days, before declining to 35 viral copies/µL at ≥8 days. In comparison with RV-A (median VL, 217 copies/µL) and RV-B (median VL, 275 copies/µL), RV-C-infected subjects produced higher VL (505 copies/µL; P < .001). Importantly, higher RV VL (median, 348 copies/µL) was associated with more severe respiratory symptoms (Total Symptom Severity Score ≥17, P = .017). A total of 83 phylogenetic-based transmission clusters were identified in the population. It was observed that the relative humidity was the strongest environmental predictor of RV seasonality in the tropical climate. Conclusions: Our findings underline the role of VL in increasing disease severity attributed to RV-C infection, and unravel the factors that fuel the population transmission dynamics of RV.

Genome Sequence of a Novel HIV-1 Circulating Recombinant Form (CRF77_cpx) Identified among Blood Donors in Malaysia.

We report here the first HIV-1 circulating recombinant form (CRF) complex identified among the blood donors in Malaysia. The CRF77_cpx mosaic genome consists of parental subtypes B', C, and CRF01_AE and is structurally related to CRF07_BC. The identification of CRF77_cpx underlines the genetic complexity and mobility of HIV-1 among the blood donors.

Identification and evolutionary dynamics of two novel human coronavirus OC43 genotypes associated with acute respiratory infections: phylogenetic, spatiotemporal and transmission network analyses.

Human coronavirus OC43 (HCoV-OC43) is commonly associated with respiratory tract infections in humans, with five genetically distinct genotypes (A to E) described so far. In this study, we obtained the full-length genomes of HCoV-OC43 strains from two previously unrecognized lineages identified among patients presenting with severe upper respiratory tract symptoms in a cross-sectional molecular surveillance study in Kuala Lumpur, Malaysia, between 2012 and 2013. Phylogenetic, recombination and comparative genomic analyses revealed two distinct clusters diverging from a genotype D-like common ancestor through recombination with a putative genotype A-like lineage in the non-structural protein (nsp) 10 gene. Signature amino acid substitutions and a glycine residue insertion at the N-terminal domain of the S1 subunit of the spike gene, among others, exhibited further distinction in a recombination pattern, to which these clusters were classified as genotypes F and G. The phylogeographic mapping of the global spike gene indicated that the genetically similar HCoV-OC43 genotypes F and G strains were potentially circulating in China, Japan, Thailand and Europe as early as the late 2000s. The transmission network construction based on the TN93 pairwise genetic distance revealed the emergence and persistence of multiple sub-epidemic clusters of the highly prevalent genotype D and its descendant genotypes F and G, which contributed to the spread of HCoV-OC43 in the region. Finally, a more consistent nomenclature system for non-recombinant and recombinant HCoV-OC43 lineages is proposed, taking into account genetic recombination as an important feature in HCoV evolution and classification.

Whole-Genome Phylogenetic Analysis of Influenza B/Phuket/3073/2013-Like Viruses and Unique Reassortants Detected in Malaysia between 2012 and 2014.

Reassortment of genetic segments between and within influenza B lineages (Victoria and Yamagata) has been shown to generate novel reassortants with unique genetic characteristics. Based on hemagglutinin (HA) and neuraminidase (NA) genes, recent surveillance study has identified reassortment properties in B/Phuket/3073/2013-like virus, which is currently used in the WHO-recommended influenza vaccine. To understand the potential reassortment patterns for all gene segments, four B/Phuket/3073/2013-like viruses and two unique reassortants (one each from Yamagata and Victoria) detected in Malaysia from 2012-2014 were subjected to whole-genome sequencing. Each gene was phylogenetically classified into lineages, clades and sub-clades. Three B/Phuket/3073/2013-like viruses from Yamagata lineage were found to be intra-clade reassortants, possessing PA and NA genes derived from Stockholm/12-like sub-clade, while the remaining genes from Wisconsin/01-like sub-clade (both sub-clades were within Yamagata Clade 3/Yam-3). However, the other B/Phuket/3073/2013-like virus had NS gene that derived from Stockholm/12-like sub-clade instead of Wisconsin/01-like sub-clade. One inter-clade reassortant had Yamagata Clade 2/Yam-2-derived HA and NP, and its remaining genes were Yam-3-derived. Within Victoria Clade 1/Vic-1 in Victoria lineage, one virus had intra-clade reassortment properties: HA and PB2 from Vic-1B sub-clade, MP and NS from a unique sub-clade "Vic-1C", and the remaining genes from Vic-1A sub-clade. Although random reassortment event may generate unique reassortants, detailed phylogenetic classification of gene segments showed possible genetic linkage between PA and NA genes in B/Phuket/3073/2013-like viruses, which requires further investigation. Understanding on reassortment patterns in influenza B evolution may contribute to future vaccine design.

Performance of a Taqman Assay for Improved Detection and Quantification of Human Rhinovirus Viral Load.

Human rhinovirus (HRV) is the major aetiology of respiratory tract infections. HRV viral load assays are available but limitations that affect accurate quantification exist. We developed a one-step Taqman assay using oligonucleotides designed based on a comprehensive list of global HRV sequences. The new oligonucleotides targeting the 5'-UTR region showed high PCR efficiency (E = 99.6%, R2 = 0.996), with quantifiable viral load as low as 2 viral copies/µl. Assay evaluation using an External Quality Assessment (EQA) panel yielded a detection rate of 90%. When tested on 315 human enterovirus-positive specimens comprising at least 84 genetically distinct HRV types/serotypes (determined by the VP4/VP2 gene phylogenetic analysis), the assay detected all HRV species and types, as well as other non-polio enteroviruses. A commercial quantification kit, which failed to detect any of the EQA specimens, produced a detection rate of 13.3% (42/315) among the clinical specimens. Using the improved assay, we showed that HRV sheds in the upper respiratory tract for more than a week following acute infection. We also showed that HRV-C had a significantly higher viral load at 2-7 days after the onset of symptoms (p = 0.001). The availability of such assay is important to facilitate disease management, antiviral development, and infection control.

Genetic diversity, seasonality and transmission network of human metapneumovirus: identification of a unique sub-lineage of the fusion and attachment genes.

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.

Diversity and Evolutionary Histories of Human Coronaviruses NL63 and 229E Associated with Acute Upper Respiratory Tract Symptoms in Kuala Lumpur, Malaysia.

The human alphacoronaviruses HCoV-NL63 and HCoV-229E are commonly associated with upper respiratory tract infections (URTI). Information on their molecular epidemiology and evolutionary dynamics in the tropical region of southeast Asia however is limited. Here, we analyzed the phylogenetic, temporal distribution, population history, and clinical manifestations among patients infected with HCoV-NL63 and HCoV-229E. Nasopharyngeal swabs were collected from 2,060 consenting adults presented with acute URTI symptoms in Kuala Lumpur, Malaysia, between 2012 and 2013. The presence of HCoV-NL63 and HCoV-229E was detected using multiplex polymerase chain reaction (PCR). The spike glycoprotein, nucleocapsid, and 1a genes were sequenced for phylogenetic reconstruction and Bayesian coalescent inference. A total of 68/2,060 (3.3%) subjects were positive for human alphacoronavirus; HCoV-NL63 and HCoV-229E were detected in 45 (2.2%) and 23 (1.1%) patients, respectively. A peak in the number of HCoV-NL63 infections was recorded between June and October 2012. Phylogenetic inference revealed that 62.8% of HCoV-NL63 infections belonged to genotype B, 37.2% was genotype C, while all HCoV-229E sequences were clustered within group 4. Molecular dating analysis indicated that the origin of HCoV-NL63 was dated to 1921, before it diverged into genotype A (1975), genotype B (1996), and genotype C (2003). The root of the HCoV-229E tree was dated to 1955, before it diverged into groups 1-4 between the 1970s and 1990s. The study described the seasonality, molecular diversity, and evolutionary dynamics of human alphacoronavirus infections in a tropical region.

Molecular epidemiology and evolutionary histories of human coronavirus OC43 and HKU1 among patients with upper respiratory tract infections in Kuala Lumpur, Malaysia.

BACKGROUND: Despite the worldwide circulation of human coronavirus OC43 (HCoV-OC43) and HKU1 (HCoV-HKU1), data on their molecular epidemiology and evolutionary dynamics in the tropical Southeast Asia region is lacking. METHODS: The study aimed to investigate the genetic diversity, temporal distribution, population history and clinical symptoms of betacoronavirus infections in Kuala Lumpur, Malaysia between 2012 and 2013. A total of 2,060 adults presented with acute respiratory symptoms were screened for the presence of betacoronaviruses using multiplex PCR. The spike glycoprotein, nucleocapsid and 1a genes were sequenced for phylogenetic reconstruction and Bayesian coalescent inference. RESULTS: A total of 48/2060 (2.4 %) specimens were tested positive for HCoV-OC43 (1.3 %) and HCoV-HKU1 (1.1 %). Both HCoV-OC43 and HCoV-HKU1 were co-circulating throughout the year, with the lowest detection rates reported in the October-January period. Phylogenetic analysis of the spike gene showed that the majority of HCoV-OC43 isolates were grouped into two previously undefined genotypes, provisionally assigned as novel lineage 1 and novel lineage 2. Sign of natural recombination was observed in these potentially novel lineages. Location mapping showed that the novel lineage 1 is currently circulating in Malaysia, Thailand, Japan and China, while novel lineage 2 can be found in Malaysia and China. Molecular dating showed the origin of HCoV-OC43 around late 1950s, before it diverged into genotypes A (1960s), B (1990s), and other genotypes (2000s). Phylogenetic analysis revealed that 27.3 % of the HCoV-HKU1 strains belong to genotype A while 72.7 % belongs to genotype B. The tree root of HCoV-HKU1 was similar to that of HCoV-OC43, with the tMRCA of genotypes A and B estimated around the 1990s and 2000s, respectively. Correlation of HCoV-OC43 and HCoV-HKU1 with the severity of respiratory symptoms was not observed. CONCLUSIONS: The present study reported the molecular complexity and evolutionary dynamics of human betacoronaviruses among adults with acute respiratory symptoms in a tropical country. Two novel HCoV-OC43 genetic lineages were identified, warranting further investigation on their genotypic and phenotypic characteristics.

Genome Sequence of Complex HIV-1 Unique Recombinant Forms Sharing a Common Recombination Breakpoint Identified in Malaysia.

Three strains of HIV-1 unique recombinant forms (URFs) descended from subtypes B, B', and CRF01_AE were identified among people who inject drugs in Kuala Lumpur, Malaysia. These three URFs shared a common recombination breakpoint in the reverse transcriptase region, indicating frequent linkage within the drug-injecting networks in Malaysia.

Co-infections and transmission networks of HCV, HIV-1 and HPgV among people who inject drugs.

Co-infections with human immunodeficiency virus type 1 (HIV-1) and human pegivirus (HPgV) are common in hepatitis C virus (HCV)-infected individuals. However, analysis on the evolutionary dynamics and transmission network profiles of these viruses among individuals with multiple infections remains limited. A total of 228 injecting drug users (IDUs), either HCV- and/or HIV-1-infected, were recruited in Kuala Lumpur, Malaysia. HCV, HIV-1 and HPgV genes were sequenced, with epidemic growth rates assessed by the Bayesian coalescent method. Based on the sequence data, mono-, dual- and triple-infection were detected in 38.8%, 40.6% and 20.6% of the subjects, respectively. Fifteen transmission networks involving HCV (subtype 1a, 1b, 3a and 3b), HIV-1 (CRF33_01B) and HPgV (genotype 2) were identified and characterized. Genealogical estimates indicated that the predominant HCV, HIV-1 and HPgV genotypes were introduced into the IDUs population through multiple sub-epidemics that emerged as early as 1950s (HCV), 1980s (HIV-1) and 1990s (HPgV). By determining the difference in divergence times between viral lineages (ΔtMRCA), we also showed that the frequency of viral co-transmission is low among these IDUs. Despite increased access to therapy and other harm reduction interventions, the continuous emergence and coexistence of new transmission networks suggest persistent multiple viral transmissions among IDUs.