Chronology of the HIV-1 CRF07_BC expansion in East Asia

Genetic history of hepatitis C virus genotype 6 in Taiwan

Unlike hepatitis C virus (HCV) genotype (GT) 6, which is widely circulated in Southeast Asia and South China, GT 6 was not reported in Taiwan until 2006. GT 1b and 2a, also known as global HCV subtypes, have been reported as major GTs circulating in Taiwan. Because of improvement in genotyping kits and sequencing techniques for the subtyping of HCV, an increasing number of GT 6 subtypes have been reported, especially subtype 6a among intravenous drug users with human immunodeficiency virus infection after an outbreak since 2003. Thus, HCV GT 6 infection is regarded to be closely associated with injection drug use. However, recently, we found an unexpectedly high GT 6 prevalence in the general population in Tainan, southern Taiwan. Most of these GT 6 samples belonged to a putative novel subtype closely related to 6g and 6w instead of 6a. Phylogenetic analyses indicated that this putative 6g-related novel subtype and 6w could be indigenous in southern Taiwan for centuries. Southern Taiwan could be the origin of HCV subtype 6w. This finding might change the perspective of HCV epidemiology in Taiwan.

Geographic Transmission and Epidemic History of HIV-1 CRF01_AE, CRF07_BC, and HCV Subtype-6w among Taiwanese Persons Who Inject Drugs

Persons who inject drugs (PWID) and their risk-related behaviors (e.g., unprotected sex and sharing needles/syringes/other injection equipment) have caused severe public health problems, especially in the rapid spread of HIV-1 and HCV. Here, we reconstructed the epidemic history of HIV-1 circulating recombinant form (CRF) 01_AE, CRF07_BC, and HCV subtype-6w among Taiwanese PWID. The timescales were estimated using phylogenetic and Bayesian coalescent analyses. The results revealed that CRF01_AE started to circulate in the Taiwanese PWID population in central Taiwan at 1992.5 (95% credible region: 1988.8–1995.9) and spread to other regions of Taiwan, while CRF07_BC was first identified in southern Taiwan at 2000.0 (95% CR: 1997.8–2002.2) and then spread northward to central-northern Taiwan. All HCV-6 strains were from Asia (that is, China, Myanmar, Taiwan, and Vietnam) and originated in 1928.1 (95% CR: 1890.2–1966.0). Furthermore, subtype-6w isolates from different regions of Taiwan appeared to share a common source that existed in the mid-1990s (95% CR: 1985.0–2001.8) or thereabouts. The routes of drug trafficking and the resulting high prevalence of HIV-1/HCV co-infections among PWID might have contributed to the virus transmission and promoted its spread worldwide. Long-term monitoring and policy implementation in at-risk populations would be useful for disease control.

Hepatitis C Virus Subtypes Novel 6g-Related Subtype and 6w Could Be Indigenous in Southern Taiwan with Characteristic Geographic Distribution

Hepatitis C virus (HCV) genotype (GT) 6 is the most genetically diverse GT and mainly distributed in Southeast Asia and south China but not Taiwan. Earlier studies showed the major HCV GTs in Taiwan were GT 1b and 2 with very rare GT 6 except in injection drug users (IDUs), and subtype 6a is the main GT 6 subtype among IDUs. Recently, we reported a much higher prevalence (18.3%) of GT 6 in Tainan City, southern Taiwan. This study was designed to clarify the subtypes of GT 6 in this endemic area. A total of 3022 (1343 men and 1679 women) HCV viremic patients were enrolled. Subtypes of GT 6 were determined by sequencing of core/E1 and nonstructural protein 5B in 322 of 518 GT 6 patients. The overall GT 6 prevalence rate was 17.1% (518/3022), with higher prevalence districts (>25%) located in northern Tainan. A novel 6g-related subtype is the most prevalent subtype (81.0%), followed by 6w (10.8%), 6a (7.5%), and 6n (0.7%). The high GT 6 prevalence in Tainan was mainly due to a novel 6g-related subtype and 6w. These two subtypes could be indigenous in Tainan with characteristic geographic distribution.

Phylogenetic and temporal dynamics of human immunodeficiency virus type 1 CRF01_AE and CRF07_BC among recently infected antiretroviral therapy-naïve men who have sex with men in Jiangsu province, China, 2012 to 2015: A molecular epidemiology-based study

The prevalence and incidence of human immunodeficiency virus type 1 (HIV-1) among men who have sex with men (MSM) are on the rise throughout China. With a large population of MSM, Jiangsu Province is facing an escalating HIV-1 epidemic.The aim of this study was to explore the phylogenetic and temporal dynamics of HIV-1 CRF01_AE and CRF07_BC among antiretroviral therapy (ART)-naïve MSM recently infected with HIV-1 in Jiangsu Province.We recruited MSM in Jiangsu Province (Suzhou, Wuxi, Nantong, Taizhou and Yancheng) 2012 to 2015. We collected information on demographics and sexual behaviors and a blood sample for HIV genome RNA extraction, RT-PCR amplification, and DNA sequencing. Multiple alignments were made using Gene Cutter, with the selected reference sequences of various subtypes/recombinants from the Los Alamos HIV-1 database. Phylogenetic and Bayesian evolutionary analysis was performed by MEGA version 6.0, Fasttree v2.1.7. and BEAST v1.6.2. Categorical variables were analyzed using χ test (or Fisher exact test where necessary). χ test with trend was used to assess the evolution of HIV-1 subtype distribution over time. All data were analyzed using SPSS20.0 software package (IBM Company, New York, NY).HIV-1 phylogenetic analysis revealed a broad viral diversity including CRF01_AE (60.06%), CRF07_BC (22.29%), subtype B (5.88%), CRF67_01B (5.26%), CRF68_01B (2.79%), CRF55_01B (1.55%), CRF59_01B (0.93%), and CRF08_BC (0.62%). Two unique recombination forms (URFs) (0.62%) were also detected. Four epidemic clusters and 1 major cluster in CRF01_AE and CRF07_BC were identified. The introduction of CRF01_AE strain (2001) was earlier than CRF07_BC strain (2004) into MSM resided in Jiangsu based on the time of the most recent common ancestor.Our study demonstrated HIV-1 subtype diversity among ART-naïve MSM recently infected with HIV-1 in Jiangsu. We first depicted the spatiotemporal dynamics, traced the dates of origin for the HIV-1 CRF01_AE/07_BC strains and made inference for the effective population size among newly infected ART-naïve MSM in Jiangsu from 2012 to 2015. A real-time surveillance of HIV-1 viral diversity and phylodynamics of epidemic cluster would be of great value to the monitoring of the epidemic and control of transmission, improvement of antiretroviral therapy strategies, and design of vaccines.

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.

HIV-1 Genetic Diversity and Its Impact on Baseline CD4+T Cells and Viral Loads among Recently Infected Men Who Have Sex with Men in Shanghai, China

The HIV-1 epidemic among men who have sex with men (MSM) has been spreading throughout China. Shanghai, a central gathering place for MSM, is facing a continuously increasing incidence of HIV-1 infection. In order to better understand the dynamics of HIV-1 diversity and its influence on patient's immune status at baseline on diagnosis, 1265 newly HIV-1-infected MSM collected from January 2009 to December 2013 in Shanghai were retrospectively analyzed for genetic subtyping, CD4+T cell counts, and viral loads. HIV-1 phylogenetic analysis revealed a broad viral diversity including CRF01_AE (62.13%), CRF07_BC (24.51%), subtype B (8.06%), CRF55_01B (3.24%), CER67_01B (0.95%), CRF68_01B (0.4%), CRF08_BC (0.08%) and CRF59_01B (0.08%). Twenty-four unique recombination forms (URFs) (1.98%) were identified as well. Bayesian inference analysis indicated that the introduction of CRF01_AE strain (1997) was earlier than CRF07_BC strain (2001) into MSM population in Shanghai based on the time of the most recent common ancestor (tMRCA). Three epidemic clusters and five sub-clusters were found in CRF01_AE. Significantly lower CD4+T cell count was found in individuals infected with CRF01_AE than in those infected with CRF07_BC infection (P<0.01), whereas viral load was significantly higher those infected with CRF01_AE than with CRF07_BC (P<0.01). In addition, the patients with >45 years of age were found to have lower CD4+T cell counts and higher viral loads than the patients with <25 years of age (P<0.05). This study reveals the presence of HIV-1 subtype diversity in Shanghai and its remarkable influence on clinical outcome. A real-time surveillance of HIV-1 viral diversity and phylodynamics of epidemic cluster, patient's baseline CD4+T cell count and viral load would be of great value to monitoring of disease progression, intervention for transmission, improvement of antiretroviral therapy strategy and design of vaccines.

Cross-border sexual transmission of the newly emerging HIV-1 clade CRF51_01B

A novel HIV-1 recombinant clade (CRF51_01B) was recently identified among men who have sex with men (MSM) in Singapore. As cases of sexually transmitted HIV-1 infection increase concurrently in two socioeconomically intimate countries such as Malaysia and Singapore, cross transmission of HIV-1 between said countries is highly probable. In order to investigate the timeline for the emergence of HIV-1 CRF51_01B in Singapore and its possible introduction into Malaysia, 595 HIV-positive subjects recruited in Kuala Lumpur from 2008 to 2012 were screened. Phylogenetic relationship of 485 amplified polymerase gene sequences was determined through neighbour-joining method. Next, near-full length sequences were amplified for genomic sequences inferred to be CRF51_01B and subjected to further analysis implemented through Bayesian Markov chain Monte Carlo (MCMC) sampling and maximum likelihood methods. Based on the near full length genomes, two isolates formed a phylogenetic cluster with CRF51_01B sequences of Singapore origin, sharing identical recombination structure. Spatial and temporal information from Bayesian MCMC coalescent and maximum likelihood analysis of the protease, gp120 and gp41 genes suggest that Singapore is probably the country of origin of CRF51_01B (as early as in the mid-1990s) and featured a Malaysian who acquired the infection through heterosexual contact as host for its ancestral lineages. CRF51_01B then spread rapidly among the MSM in Singapore and Malaysia. Although the importation of CRF51_01B from Singapore to Malaysia is supported by coalescence analysis, the narrow timeframe of the transmission event indicates a closely linked epidemic. Discrepancies in the estimated divergence times suggest that CRF51_01B may have arisen through multiple recombination events from more than one parental lineage. We report the cross transmission of a novel CRF51_01B lineage between countries that involved different sexual risk groups. Understanding the cross-border transmission of HIV-1 involving sexual networks is crucial for effective intervention strategies in the region.

A new migration map of HIV-1 CRF07_BC in China: analysis of sequences from 12 provinces over a decade

BACKGROUND

As one of prevalence HIV-1 strains in IDUs in Asia, the origination and full transmission map of CRF07_BC is of great interested and remains unclear. In the study, we collected 769 CRF07_BC derived sequences (including 45 sequences generated in our laboratory) from 12 provinces in China for reconstructing transmission map. Meanwhile, ample historic epidemic evidences were also reviewed to assist sequences analysis.

METHODOLOGY/PRINCIPAL FINDINGS

In the study, we collected 769 CRF07_BC derived sequences and identified 138 independent sequences from 12 provinces in China for subsequent phylogeographic tree analysis, Bayes Factor test and the estimation of state tMRCA. The analyses demonstrated that CRF07_BC was originated in 1993 in IDU in Yunnan province and then initially spread to Guangxi (eastern neighbor to Yunnan) in 1994, to Xinjiang (northwest) in 1995 and to Sichuan (northern neighbor to Yunnan) in 1996. The subsequent transmissions occurred from Yunnan to Liaoning (northeast) in 1997 and to Jiangsu in 1998. Interestingly, after the early introduction of CRF07_BC into Guangxi, Xinjiang and Sichuan, these three regions served as secondary epicenters for further spreading into Gansu, Ningxia, Qinghai, Beijing and Hunan during 1999-2001. These analyzed results are in accordance with early epidemic investigations.

CONCLUSIONS/SIGNIFICANCE

Our data not just reconstructed the migration map of CRF07_BC, but also firstly revealed the active roles of these secondary epicenters in the dynamic migration of CRF07_BC in China.

Phylodynamic inference of infectious diseases caused by HIV, enterovirus 71 and the 2009 swine-origin human influenza virus

The rapidly evolving nature of viruses, particularly RNA viruses, exhibit complex phylodynamic behaviors, accumulating mutations across the genome on a timescale observable by humans. Phylodynamic investigations provide unique and quantitative evolutionary data on the origin, progression in terms of time and space, and epidemic and transmission history of pathogens. This review describes recent phylodynamic conjectures on three distinct categories of human viruses, including HIV, enterovirus 71 and the swine-origin pandemic influenza virus H1N1, which are of enormous importance to public health.

The human immunodeficiency virus-1 genotype diversity and drug resistance mutations profile of volunteer blood donors from Chinese blood centers

BACKGROUND

The global human immunodeficiency virus (HIV)-1 epidemic is becoming increasingly diverse and complex. Molecular epidemiologic characteristics were studied for HIV-1-infected blood donors from five Chinese regions to determine genotype diversity and drug resistance mutations (DRMs) profile.

STUDY DESIGN AND METHODS

HIV-1 confirmed-reactive serum samples were collected from 172 blood donors from five blood centers during 2007 to 2010. HIV-1 Pol including whole protease and partial reverse transcriptase genes was amplified, sequenced, and analyzed for the subtype determination and drug resistance profile description.

RESULT

A total of 113 amplified sequences including 82 from Kunming blood center and 31 from four other blood centers had the following genotype characteristics: G (0.9%), B (2.7%), circulating recombinant form (CRF) 01_AE (32.7%), CRF07_BC (22.1%), and CRF08_BC (41.6%). Female donors represent 45.1% of all cases and 63.9% cases with DRMs. The prevalence of samples with potential low or higher resistance among Chinese blood donors is 4.4%.

CONCLUSION

HIV-1 infection in Chinese blood donors is genetically diverse and the subtype distribution reflects that from the high-risk populations. Our results support continuous molecular epidemiologic surveillance for HIV-1 in blood donors as a part of a comprehensive HIV control program.

Phylogeographic analyses reveal a crucial role of Xinjiang in HIV-1 CRF07_BC and HCV 3a transmissions in Asia

Background

China faces an increasing prevalence of two HIV-1 circulating recombinant forms (CRFs) 07_BC and 08_BC. Both CRFs_BC were previously demonstrated to originate in Yunnan and spread to Liaoning from Yunnan via injection drug use (IDU) in China. Supposing it is true, we are unable to answer why only CRF07_BC, rather than both CRFs_BC together, was transmitted to Xinjiang.

Methodology/Principal Findings

We investigated the phylogeography of CRF07_BC and CRF08_BC using multiple HIV-1 genomic regions with Bayesian phylogeography method. Phylogenetic reconstructions showed that all CRF07_BC sequences were divided into two clades, Yunnan and Xinjiang, and all strains from other regions of mainland China clustered within the Xinjiang clade. Significant geographic diffusion links of Xinjiang with other regions (including Liaoning, Beijing, Jiangsu and Guangdong) were supported by Bayes factor tests. The temporal dynamics analyses showed that CRF07_BC spread from Xinjiang to Liaoning in 1996.10, and to Jiangsu in 2000.9. The analyses of CRF08_BC not only confirmed the previous conclusion on temporal and spatial dynamics of CRF08_BC, but also indicated that the CRF08_BC strains from Guangdong and Shanghai originated from Yunnan. The analyses of HCV 3a showed that it was introduced into Xinjiang in the early 1980s, and spread from Xinjiang to Yunnan in 1990.10 and to Jiangsu in 1999.2, and further from Yunnan to Guangxi in 1995.3. The temporal and spatial dynamics of HCV 3a were similar to some extent to that of HIV-1 CRF07_BC and/or CRF08_BC, suggesting a possible association in migration patterns between HCV and HIV-1 through IDU. In addition, HCV 3a spread from Xinjiang to Pakistan, implying a drug trafficking route linking them.

Conclusions/Significance

Xinjiang, as the most important transfer station for drug trafficking from Golden Crescent to other regions of China, plays a very crucial role in the transmission of viruses (e.g., HIV-1 and HCV) through IDU in Asia.

New trends of HCV infection in China revealed by genetic analysis of viral sequences determined from first-time volunteer blood donors

Recently, we studied hepatitis C virus (HCV) sera-prevalence among 559 890 first-time volunteer blood donors in China. From randomly selected 450 anti-HCV positive donors, we detected HCV RNA in 270 donors. In this study, we amplified HCV E1 and/or NS5B sequences from 236 of these donors followed by DNA sequencing and phylogenetic analysis. The results indicate new trends of HCV infection in China. The HCV genotype distribution differed according to the donors' region of origin. Among donors from Guangdong province, we detected subtypes 6a, 1b, 3a, 3b, 2a, and 1a at frequencies of 49.7%, 31.0%, 7.6%, 5.5%, 4.1%, and 2.1%, respectively. Among donors from outside Guangdong, we detected 1b, 2a, 6a, 3b, 3a, 6e, and 6n at frequencies 57.1%, 13.2%, 11.0%, 9.9%, 4.4%, 2.2%, and 2.2%, respectively. Although we found no significant differences among regions in age or gender, subtype 6a was more common (P < 0.001) in donors from Guangdong than those from elsewhere, whilst subtypes 1b (P < 0.02) and 2a (P < 0.001) were more frequent outside Guangdong. Disregarding origins, the male/female ratio was higher for subtype 6a-infected donors (P < 0.05) than for subtype 1b donors, whilst the mean age of subtype 2a donors was 8-10 years older (P < 0.05) than that for all other subtypes. Detailed phylogenetic analysis of our sequence data provides further insight into the transmission of HCV within China, and between China and other countries. The predominance of HCV 6a among blood donors in Guangdong is striking and mandates studies into risk factors for its acquisition.

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.

Molecular epidemiology of HCV genotypes among injection drug users in Taiwan: Full-length sequences of two new subtype 6w strains and a recombinant form_2b6w

Human immunodeficiency virus type 1 (HIV-1) circulating recombinant form (CRF) 07_BC strain has caused serious outbreaks among injection drug users in Taiwan since 2004. The objective of this study was to conduct a molecular epidemiological study of HCV genotypes in intravenous drug users in Taiwan. Blood samples and questionnaires from 591 intravenous drug users infected with HIV-1 were collected nationwide. In total, 180 samples were selected for HCV genotyping using multiplex PCR and phylogenetic analysis of the core, E1 and NS5B regions. The Inno-Lipa assay was used to confirm multiple infections with different genotypes. Eighty percent had a single infection with subtype 1b being the most common subtype (24%), 12% had double infections and two had triple infections. In addition, three recombinant forms (RFs)-2a1a, 3a1b, and 2b6w were identified. Phylogenetic analyses showed that the 3a, 6a, and 6n strains were clustered with strains present in Thailand and mainland China. Full-length sequence analysis showed that two 6w strains shared 89.4-90.2% sequence homology with the 6(r) strain from the Guangdong Province, China. Bootscan analysis revealed that the recombination breakpoint of RF_2b6w was located at the NS2-NS3 junction. In summary, the distribution of HCV genotypes among Taiwanese intravenous drug users was complex and more than 12% of the drug users were infected with more than one genotype of HCV.

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.

Temporal and spatial dynamics of human immunodeficiency virus type 1 circulating recombinant forms 08_BC and 07_BC in Asia

Human immunodeficiency virus type 1 (HIV-1) CRF08_BC and CRF07_BC are two major recombinants descended from subtypes B' and C. Despite their massive epidemic impact in China, their migration patterns and divergence times remain unknown. Phylogenetic and population genetic analyses were performed on 228 HIV-1 sequences representing CRF08_BC, CRF07_BC, and subtype C strains from different locations across China, India, and Myanmar. Genome-specific rates of evolution and divergence times were estimated using a Bayesian Markov chain Monte Carlo framework under various evolutionary models. CRF08_BC originated in 1990.3 (95% credible region [CR], 1988.6 to 1991.9) in Yunnan province before spreading to Guangxi (south) and Liaoning (northeast) around 1995. Inside Guangxi region, the eastward expansion of CRF08_BC continued from Baise city (west) to Binyang (central) between 1997 and 1998 and later spread into Pingxiang around 1999 in the south, mainly through injecting drug users. Additionally, CRF07_BC diverged from its common ancestor in 1993.3 (95% CR, 1991.2 to 1995.2) before crossing the border into southern Taiwan in late 1990s. Phylogenetic analysis indicates that both CRF08_BC and CRF07_BC can trace their origins to Yunnan. The parental Indian subtype C lineage likely entered China around 1981.2 (95% CR, 1976.7 to 1985.9). Using a multiple unlinked locus model, we also showed that the dates of divergence calculated in this study may not be significantly affected by intrasubtype recombination among different lineages. This is the first phylodynamic study depicting the spatiotemporal dynamics of HIV/AIDS in East Asia.