Kok Keng Tee

Evolutionary Genetics of Human Enterovirus 71: Origin, Population Dynamics, Natural Selection, and Seasonal Periodicity of the VP1 Gene

ABSTRACT Human enterovirus 71 (EV-71) is one of the major etiologic causes of hand, foot, and mouth disease (HFMD) among young children worldwide, with fatal instances of neurological complications becoming increasingly common. Global VP1 capsid sequences (n = 628) sampled over 4 decades were collected and subjected to comprehensive evolutionary analysis using a suite of phylogenetic and population genetic methods. We estimated that the common ancestor of human EV-71 likely emerged around 1941 (95% confidence interval [CI], 1929 to 1952), subsequently diverging into three genogroups: B, C, and the now extinct genogroup A. Genealogical analysis revealed that diverse lineages of genogroup B and C (subgenogroups B1 to B5 and C1 to C5) have each circulated cryptically in the human population for up to 5 years before causing large HFMD outbreaks, indicating the quiescent persistence of EV-71 in human populations. Estimated phylogenies showed a complex pattern of spatial structure within well-sampled subgenogroups, suggesting endemicity with occasional lineage migration among locations, such that past HFMD epidemics are unlikely to be linked to continuous transmission of a single strain of virus. In addition, rises in genetic diversity are correlated with the onset of epidemics, driven in part by the emergence of novel EV-71 subgenogroups. Using subgenogroup C1 as a model, we observe temporal strain replacement through time, and we investigate the evidence for positive selection at VP1 immunogenic sites. We discuss the consequences of the evolutionary dynamics of EV-71 for vaccine design and compare its phylodynamic behavior with that of influenza virus.

Geographic and Temporal Trends in the Molecular Epidemiology and Genetic Mechanisms of Transmitted HIV-1 Drug Resistance: An Individual-Patient- and Sequence-Level Meta-Analysis

Background Regional and subtype-specific mutational patterns of HIV-1 transmitted drug resistance (TDR) are essential for informing first-line antiretroviral (ARV) therapy guidelines and designing diagnostic assays for use in regions where standard genotypic resistance testing is not affordable. We sought to understand the molecular epidemiology of TDR and to identify the HIV-1 drug-resistance mutations responsible for TDR in different regions and virus subtypes. Methods and Findings We reviewed all GenBank submissions of HIV-1 reverse transcriptase sequences with or without protease and identified 287 studies published between March 1, 2000, and December 31, 2013, with more than 25 recently or chronically infected ARV-naïve individuals. These studies comprised 50,870 individuals from 111 countries. Each set of study sequences was analyzed for phylogenetic clustering and the presence of 93 surveillance drug-resistance mutations (SDRMs). The median overall TDR prevalence in sub-Saharan Africa (SSA), south/southeast Asia (SSEA), upper-income Asian countries, Latin America/Caribbean, Europe, and North America was 2.8%, 2.9%, 5.6%, 7.6%, 9.4%, and 11.5%, respectively. In SSA, there was a yearly 1.09-fold (95% CI: 1.05–1.14) increase in odds of TDR since national ARV scale-up attributable to an increase in non-nucleoside reverse transcriptase inhibitor (NNRTI) resistance. The odds of NNRTI-associated TDR also increased in Latin America/Caribbean (odds ratio [OR] = 1.16; 95% CI: 1.06–1.25), North America (OR = 1.19; 95% CI: 1.12–1.26), Europe (OR = 1.07; 95% CI: 1.01–1.13), and upper-income Asian countries (OR = 1.33; 95% CI: 1.12–1.55). In SSEA, there was no significant change in the odds of TDR since national ARV scale-up (OR = 0.97; 95% CI: 0.92–1.02). An analysis limited to sequences with mixtures at less than 0.5% of their nucleotide positions—a proxy for recent infection—yielded trends comparable to those obtained using the complete dataset. Four NNRTI SDRMs—K101E, K103N, Y181C, and G190A—accounted for >80% of NNRTI-associated TDR in all regions and subtypes. Sixteen nucleoside reverse transcriptase inhibitor (NRTI) SDRMs accounted for >69% of NRTI-associated TDR in all regions and subtypes. In SSA and SSEA, 89% of NNRTI SDRMs were associated with high-level resistance to nevirapine or efavirenz, whereas only 27% of NRTI SDRMs were associated with high-level resistance to zidovudine, lamivudine, tenofovir, or abacavir. Of 763 viruses with TDR in SSA and SSEA, 725 (95%) were genetically dissimilar; 38 (5%) formed 19 sequence pairs. Inherent limitations of this study are that some cohorts may not represent the broader regional population and that studies were heterogeneous with respect to duration of infection prior to sampling. Conclusions Most TDR strains in SSA and SSEA arose independently, suggesting that ARV regimens with a high genetic barrier to resistance combined with improved patient adherence may mitigate TDR increases by reducing the generation of new ARV-resistant strains. A small number of NNRTI-resistance mutations were responsible for most cases of high-level resistance, suggesting that inexpensive point-mutation assays to detect these mutations may be useful for pre-therapy screening in regions with high levels of TDR. In the context of a public health approach to ARV therapy, a reliable point-of-care genotypic resistance test could identify which patients should receive standard first-line therapy and which should receive a protease-inhibitor-containing regimen.

The unique HCV genotype distribution and the discovery of a novel subtype 6u among IDUs co-infected with HIV-1 in Yunnan, China†‡

The Yunnan province is the epicenter of HIV‐1 epidemics in China and a center for drug trafficking to the other parts of the world. In six prefectures of this province, a total of 132 IDUs were recruited to determine the sero‐prevalence of HCV and HIV‐1 and the positive rates were 93.94% and 68.18%, respectively (P < 0.001). Co‐infection with HCV and HIV‐1 was found among 89 IDUs, of whom several HCV fragments were amplified and sequenced. Sequences of the HCV 5′NCR‐C and NS5B region were determined from 82 IDUs. Phylogenetic analyses showed consistent genotyping among 80 IDUs. Among them HCV genotypes 1a, 1b, 3a, 3b, 6a, 6n, and a tentatively assigned novel 6u subtype were found in 1 (1.25%), 16 (20%), 19 (23.75%), 24 (30%), 4 (5%), 9 (11.25%) and 7 (8.75%) individuals, respectively. In two IDUs, genotyping results were discordant, suggesting mixed HCV infections or recombination. The proportion of patients with HCV 1b tended to decrease from the north to south and from the east to west in this province. Genotype 3 and 6 strains were more frequent in the southern prefectures. The novel subtype 6u strains were only detected in Dehong which borders Myanmar. Our findings showed a unique pattern of HCV genotype distribution, which is similar to that in the southeastern Asian countries but distinct from that among the general population in China. Routes of drug trafficking and the resulting high prevalence of HIV‐1 infection may have contributed to this pattern of HCV genotype distribution. J. Med. Virol. 80: 1142–1152, 2008.

Identification of a novel circulating recombinant form (CRF33_01B) disseminating widely among various risk populations in Kuala Lumpur, Malaysia

Summary:A molecular epidemiological investigation was conducted among various risk populations (n = 184) in Kuala Lumpur, Malaysia, in 2003 to 2005, on the basis of nucleotide sequences of protease and reverse transcriptase regions. In addition to circulating HIV-1 strains, including CRF01_AE (57.1%), subtype B (20.1%), and subtype C (0.5%), we detected a candidate with a new circulating recombinant form (CRF). We determined four near-full-length nucleotide sequences with identical subtype structure from epidemiologically unlinked individuals of different risk and ethnic groups. In this chimera, two short subtype B segments were inserted into the gag-RT region in a backbone of CRF01_AE. The recombinant structure was distinct from previously identified CRF15_01B in Thailand. In agreement with the current HIV nomenclature system, this constitutes a novel CRF (CRF33_01B). The overall prevalence of CRF33_01B is 19.0% (35/184). Although the prevalence of CRF33_01B is particularly high among injecting drug users (42.0%, 21/50), it is also detected in a substantial proportion of homo-/bisexual males (18.8%, 3/16) and heterosexuals (9.8%, 9/92). Moreover, unique recombinant forms composed of CRF01_AE and subtype B that have a significant structural relationship with CRF33_01B were detected in 1.6% (3/184) of study subjects, suggesting an ongoing recombination process in Malaysia. This new CRF seems to be bridging viral transmission between different risk populations in this country.

Phylodynamic analysis of the dissemination of HIV-1 CRF01_AE in Vietnam.

To estimate the epidemic history of HIV-1 CRF01_AE in Vietnam and adjacent Guangxi, China, we determined near full-length nucleotide sequences of CRF01_AE from a total of 33 specimens collected in 1997-1998 from different geographic regions and risk populations in Vietnam. Phylogenetic and Bayesian molecular clock analyses were performed to estimate the date of origin of CRF01_AE lineages. Our study reconstructs the timescale of CRF01_AE expansion in Vietnam and neighboring regions and suggests that the series of CRF01_AE epidemics in Vietnam arose by the sequential introduction of founder strains into new locations and risk groups. CRF01_AE appears to have been present among heterosexuals in South-Vietnam for more than a decade prior to its epidemic spread in the early 1990s. In the late 1980s, the virus spread to IDUs in Southern Vietnam and subsequently in the mid-1990s to IDUs further north. Our results indicate the northward dissemination of CRF01_AE during this time.

Emergence of HIV-1 CRF01_AE/B unique recombinant forms in Kuala Lumpur, Malaysia

Objectives:To investigate the molecular epidemiology of HIV-1 and to screen for the emergence of intersubtype recombinants in Kuala Lumpur, Malaysia. Design:A molecular epidemiology study was conducted among HIV-1 seropositive patients attending the University Malaya Medical Center (UMMC) from July 2003 to June 2004. Methods:Protease (PR) and reverse transcriptase (RT) gene sequences were derived from drug resistance genotyping assay of 100 newly diagnosed or antiretroviral-naive patients. These were phylogenetically analysed to determine the subtypes and recombination breakpoint analyses were performed on intersubtype recombinants to estimate the recombination breakpoint(s). Results:CRF01_AE predominated in Kuala Lumpur with 65% in both PR and RT genes. B subtype was detected at 14% and 12% in PR and RT genes, respectively. C subtype was present at 1% in both genes. Overall, the concordance of PR and RT genes in discriminating subtypes/circulating recombinant forms (CRF) was high at 96%. In this study, novel CRF01_AE/B intersubtype recombinants were detected at high prevalence (22%), including those isolates with subtype discordance. Thai variants of CRF01_AE and B subtype were involved in the genesis of these unique recombinant forms (URF). Interestingly, 19 CRF01_AE/B intersubtype recombinant isolates shared similar recombination breakpoints in both PR and RT genes. Several distinct URF were also identified. Conclusion:PR and RT genes can be utilized for subtype/CRF assessment with high degree of agreement, allowing concurrent surveillance of circulating HIV-1 subtypes with antiretroviral drug resistance genotyping tests. The emergence of highly identical CRF01_AE/B intersubtype recombinants suggests the possibility of the appearance of a new circulating recombinant form in Kuala Lumpur.

Identification of a novel second-generation circulating recombinant form (CRF48_01B) in Malaysia: a descendant of the previously identified CRF33_01B.

A molecular epidemiological investigation conducted among injecting drug users in eastern Peninsular Malaysia in 2007 identified a cluster of sequences (n = 3) located outside any known HIV-1 genotype. Analyses of near full-length nucleotide sequences of these strains from individuals with no recognizable linkage revealed that they have an identical subtype structure comprised of CRF01_AE and subtype B′, distinct from any known circulating recombinant forms (CRFs). This novel CRF, designated CRF48_01B, is closely related to CRF33_01B, previously identified in Kuala Lumpur. Phylogenetic analysis of multiple CRF48_01B genome regions showed that CRF48_01B forms a monophyletic cluster within CRF33_01B, suggesting that this new recombinant is very likely a descendant of CRF33_01B. CRF48_01B thus represents one of the first examples of a “second-generation” CRF, generated by additional crossover with pre-existing CRFs. Corroborating these results, Bayesian molecular clock analyses indicated that CRF48_01B emerged in ∼2001, approximately ∼8 years after the emergence of CRF33_01B.

Genetic characterization of the Wyeomyia group of orthobunyaviruses and their phylogenetic relationships

Phylogenetic analyses can give new insights into the evolutionary history of viruses, especially of viruses with segmented genomes. However, sequence information for many viral families or genera is still limited and phylogenies based on single or short genome fragments can be misleading. We report the first genetic analysis of all three genome segments of Wyeomyia group viruses Wyeomyia, Taiassui, Macaua, Sororoca, Anhembi and Cachoeira Porteira (BeAr328208) in the genus Orthobunyavirus of the family Bunyaviridae. In addition, Tucunduba and Iaco viruses were identified as members of the Wyeomyia group. Features of Wyeomyia group members that distinguish them from other viruses in the Bunyamwera serogroup and from other orthobunyaviruses, including truncated NSs sequences that may not counteract the hosts interferon response, were characterized. Our findings also suggest genome reassortment within the Wyeomyia group, identifying Macaua and Tucunduba viruses as M-segment reassortants that, in the case of Tucunduba virus, may have altered pathogenicity, stressing the need for whole-genome sequence information to facilitate characterization of orthobunyaviruses and their phylogenetic relationships.

Genome Sequence of a Novel HIV-1 Circulating Recombinant Form 54_01B from Malaysia

ABSTRACT We report here the first novel HIV-1 circulating recombinant form (CRF) 54_01B (CRF54_01B) isolated from three epidemiologically unlinked subjects of different risk groups in Malaysia. These recently sampled recombinants showed a complex genome organization composed of parental subtype B′ and CRF01_AE, with identical recombination breakpoints observed in the gag, pol, and vif genes. Such a discovery highlights the ongoing active generation and spread of intersubtype recombinants involving the subtype B′ and CRF01_AE lineages and indicates the potential of the new CRF in bridging HIV-1 transmission among different risk groups in Southeast Asia.

Reconstructing the epidemic history of HIV-1 circulating recombinant forms CRF07_BC and CRF08_BC in East Asia: the relevance of genetic diversity and phylodynamics for vaccine strategies.

HIV-1 CRF07_BC and CRF08_BC are closely related circulating recombinant forms (CRFs) with serious public health consequences in China. The temporal and spatial dynamics of these CRFs were determined by estimating their times of divergence, using phylogenetic and Bayesian coalescent methods. Studies of the timelines of CRF07_BC and CRF08_BC trace the expansion of these strains back their origins to Yunnan province. The present study highlights the relevance of incorporating evolutionary and molecular epidemiological analyses into an in-depth understanding of the genesis of HIV epidemic, providing information for determining regional and global public health policies, including future vaccine strategies.