Analyzing Patient-Sharing Network Using an Administrative Claim Database in Japan.

We analyzed the behavior of patient with a focus on patient-sharing based on the methodology of network analysis. We used an administrative healthcare claims database from September of the years 2008-2020 to identify shared patients with hypertension. The patients' behavior of visiting multiple medical facilities was extracted as graphical data, and we calculated density and centrality as indicators to evaluate the structure of the patient sharing network. Our findings indicate that density, reciprocity, and transitivity increased over time, and that centrality and PageRank were correlated.

Changes in Demand Volume and Patient/Health Care Provider Characteristics of First-Time Telehealth Users: A Comparative Analysis Before and After the COVID-19 Policy Response Using the Administrative Claims Database.

Introduction: The COVID-19 pandemic has led to a decrease in demand for medical services in Japan, but the utilization of telehealth, which the Japanese government has recently promoted, has seen a temporary increase. This study aims to analyze the trend of telehealth utilization and changes in patient characteristics following the policy response to COVID-19. Methods: This retrospective study analyzed data from 26,152 adult patients who used telehealth for the first time between April 2019 and April 2021 in Mie Prefecture, Japan. An interrupted time series analysis was conducted to evaluate changes in the number of first-time patients before and after April 2020. Results: The number of telehealth users increased by 111.87% after April 2020, but the trend showed a declining slope thereafter. Patient characteristics and disease types showed different trends. The percentage of patients choosing a hospital over a clinic increased for the first time. Conclusions: After the policy response to COVID-19, the number of first-time telehealth users overall increased immediately, but gradually showed a declining trend. However, some diseases have shown both an immediate increase and a continued upward trend in telehealth utilization. Patients with these diseases may be candidates for adopting telehealth services in clinical settings.

Exploring Current Practice, Knowledge, and Challenges of Sexually Transmitted Infection/HIV Management and Pre-Exposure Prophylaxis Among Japanese Health Care Professionals: A Cross-Sectional Web Survey.

We conducted a web-based survey targeting physicians in specialties of treating sexually transmitted infection (STI) and/or human immunodeficiency virus (HIV) patients to understand the current STI/HIV care practices and their acceptability of and barriers to the prescription of pre-exposure prophylaxis (PrEP) in Japan. A descriptive analysis was used to summarize survey responses. Univariate and multivariable logistic regression were performed to identify factors associated with willingness to prescribe PrEP. Of 316 survey respondents, 57 were specialized in HIV, 90 STI/Urology/Proctology, 55 Obstetrics/Gynecology, and 114 General Practice/Internal Medicine/Dermatology. Proportion of HIV-specialized physicians who interview the patients about risk behaviors tended to be higher than other physician groups (84.2% vs. 54.8%, 47.3%, and 50.9%, respectively), and 53 - 75% of non-HIV-specialized physicians reported that they were incapable of making decisions on HIV medications. Higher PrEP knowledge enhanced the willingness to recommend and prescribe PrEP drugs (odds ratio: 2.31, 95% confidence interval: 1.30-4.10, p = 0.0044), and 45.4% physicians with no PrEP knowledge raised the concern of incapability to respond and manage when an individual is infected with HIV. Educational opportunities on management and prevention measures for both STI and HIV may encourage non-HIV-specialized physicians to be involved in HIV care and to enhance initiation of HIV tests and adoption of PrEP.

The Association Between Telehealth Utilization and Policy Responses on COVID-19 in Japan: Interrupted Time-Series Analysis.

BACKGROUND: Telehealth using telephones or online communication is being promoted as a policy initiative in several countries. However, there is a lack of research on telehealth utilization in a country such as Japan that offers free access to medical care and regulates telehealth provision-particularly with respect to COVID-19. OBJECTIVE: The present study aimed to clarify telehealth utilization, the characteristics of patients and medical institutions using telehealth, and the changes to telehealth in Japan in order to support the formulation of policy strategies for telehealth provision. METHODS: Using a medical administrative claim database of the National Health Insurance and Advanced Elderly Medical Service System in Mie Prefecture, we investigated patients who used telehealth from January 2017 to September 2021. We examined telehealth utilization with respect to both patients and medical institutions, and we determined their characteristics. Using April 2020 as the reference time point for COVID-19, we conducted an interrupted time-series analysis (ITSA) to assess changes in the monthly proportion of telehealth users to beneficiaries. RESULTS: The number of telehealth users before the reference time point was 13,618, and after the reference time point, it was 28,853. Several diseases and conditions were associated with an increase in telehealth utilization. Telehealth consultations were mostly conducted by telephone and for prescriptions. The ITSA results showed a sharp increase in the proportion of telehealth use to beneficiaries after the reference time point (rate ratio 2.97; 95% CI 2.14-2.31). However, no apparent change in the trend of increasing or decreasing telehealth use was evident after the reference time point (rate ratio 1.00; 95% CI 1.00-1.01). CONCLUSIONS: We observed a sharp increase in telehealth utilization after April 2020, but no change in the trend of telehealth use was evident. We identified changes in the characteristics of patients and providers using telehealth.

Short Communication: HIV-DRLink: A Tool for Reporting Linked HIV-1 Drug Resistance Mutations in Large Single-Genome Data Sets Using the Stanford HIV Database.

The prevalence of HIV-1 drug resistance is increasing worldwide and monitoring its emergence is important for the successful management of populations receiving combination antiretroviral therapy. It is likely that pre-existing drug resistance mutations linked on the same viral genomes are predictive of treatment failure. Because of the large number of sequences generated by ultrasensitive single-genome sequencing (uSGS) and other similar next-generation sequencing methods, it is difficult to assess each sequence individually for linked drug resistance mutations. Several software/programs exist to report the frequencies of individual mutations in large data sets, but they provide no information on linkage of resistance mutations. In this study, we report the HIV-DRLink program, a research tool that provides resistance mutation frequencies as well as their genetic linkage by parsing and summarizing the Sierra output from the Stanford HIV Database. The HIV-DRLink program should only be used on data sets generated by methods that eliminate artifacts due to polymerase chain reaction recombination, for example, standard single-genome sequencing or uSGS. HIV-DRLink is exclusively a research tool and is not intended to inform clinical decisions.

Nation-Wide Viral Sequence Analysis of HIV-1 Subtype B Epidemic in 2003-2012 Revealed a Contribution of Men Who Have Sex With Men to the Transmission Cluster Formation and Growth in Japan.

Background: To better understand the epidemiology of human immunodeficiency virus type 1 (HIV-1) subtype B transmission in Japan, phylodynamic analysis of viral pol sequences was conducted on individuals newly diagnosed as HIV-1 seropositive. Methodology: A total of 5,018 patients newly diagnosed with HIV-1 infection and registered in the Japanese Drug Resistance HIV Surveillance Network from 2003 to 2012 were enrolled in the analysis. Using the protease-reverse transcriptase nucleotide sequences, their subtypes were determined, and phylogenetic relationships among subtype B sequences were inferred using three different methods: distance-matrix, maximum likelihood, and Bayesian Markov chain Monte Carlo. Domestically spread transmission clusters (dTCs) were identified based on the following criteria: >95% in interior branch test, >95% in Bayesian posterior probability and <10% in depth-first searches for sub-tree partitions. The association between dTC affiliation and individuals' demographics was analyzed using univariate and multivariate analyses. Results: Among the cases enrolled in the analysis, 4,398 (87.6%) were classified as subtype B. Many of them were Japanese men who had sex with men (MSM), and 3,708 (84.3%) belonged to any of 312 dTCs. Among these dTCs, 243 (77.9%) were small clusters with <10 individuals, and the largest cluster consisted of 256 individuals. Most dTCs had median time of the most recent common ancestor between 1995 and 2005, suggesting that subtype B infection was spread among MSMs in the second half of the 1990s. Interestingly, many dTCs occurred within geographical regions. Comparing with singleton cases, TCs included more MSM, young person, and individuals with high CD4+ T-cell count at the first consultation. Furthermore, dTC size was significantly correlated with gender, age, transmission risks, recent diagnosis and relative population size of the region mainly distributed. Conclusions: Our study clarified that major key population of HIV-1 subtype B epidemic in Japan is local MSM groups. The study suggests that HIV-1 subtype B spread via episodic introductions into the local MSM groups, some of the viruses spread to multiple regions. Many cases in dTC were diagnosed during the early phase of infection, suggesting their awareness to HIV risks.

Ultrasensitive single-genome sequencing: accurate, targeted, next generation sequencing of HIV-1 RNA.

BACKGROUND: Although next generation sequencing (NGS) offers the potential for studying virus populations in unprecedented depth, PCR error, amplification bias and recombination during library construction have limited its use to population sequencing and measurements of unlinked allele frequencies. Here we report a method, termed ultrasensitive Single-Genome Sequencing (uSGS), for NGS library construction and analysis that eliminates PCR errors and recombinants, and generates single-genome sequences of the same quality as the "gold-standard" of HIV-1 single-genome sequencing assay but with more than 100-fold greater depth. RESULTS: Primer ID tagged cDNA was synthesized from mixtures of cloned BH10 wild-type and mutant HIV-1 transcripts containing ten drug resistance mutations. First, the resultant cDNA was divided and NGS libraries were generated in parallel using two methods: uSGS and a method applying long PCR primers to attach the NGS adaptors (LP-PCR-1). Second, cDNA was divided and NGS libraries were generated in parallel comparing 3 methods: uSGS and 2 methods adapted from more recent reports using variations of the long PCR primers to attach the adaptors (LP-PCR-2 and LP-PCR-3). Consistently, the uSGS method amplified a greater proportion of cDNAs, averaging 30% compared to 13% for LP-PCR-1, 21% for LP-PCR-2 and 14% for LP-PCR-3. Most importantly, when the uSGS sequences were binned according to their primer IDs, 94% of the bins did not contain PCR recombinant sequences versus only 55, 75 and 65% for LP-PCR-1, 2 and 3, respectively. Finally, when uSGS was applied to plasma samples from HIV-1 infected donors, both frequent and rare variants were detected in each sample and neighbor-joining trees revealed clusters of genomes driven by the linkage of these mutations, showing the lack of PCR recombinants in the datasets. CONCLUSIONS: The uSGS assay can be used for accurate detection of rare variants and for identifying linkage of rare alleles associated with HIV-1 drug resistance. In addition, the method allows accurate in-depth analyses of the complex genetic relationships of viral populations in vivo.

Clonally expanded CD4+ T cells can produce infectious HIV-1 in vivo.

Reservoirs of infectious HIV-1 persist despite years of combination antiretroviral therapy and make curing HIV-1 infections a major challenge. Most of the proviral DNA resides in CD4(+)T cells. Some of these CD4(+)T cells are clonally expanded; most of the proviruses are defective. It is not known if any of the clonally expanded cells carry replication-competent proviruses. We report that a highly expanded CD4(+) T-cell clone contains an intact provirus. The highly expanded clone produced infectious virus that was detected as persistent plasma viremia during cART in an HIV-1-infected patient who had squamous cell cancer. Cells containing the intact provirus were widely distributed and significantly enriched in cancer metastases. These results show that clonally expanded CD4(+)T cells can be a reservoir of infectious HIV-1.

Unique Flap Conformation in an HIV-1 Protease with High-Level Darunavir Resistance.

Darunavir (DRV) is one of the most powerful protease inhibitors (PIs) for treating human immunodeficiency virus type-1 (HIV-1) infection and presents a high genetic barrier to the generation of resistant viruses. However, DRV-resistant HIV-1 infrequently emerges from viruses exhibiting resistance to other protease inhibitors. To address this resistance, researchers have gathered genetic information on DRV resistance. In contrast, few structural insights into the mechanism underlying DRV resistance are available. To elucidate this mechanism, we determined the crystal structure of the ligand-free state of a protease with high-level DRV resistance and six DRV resistance-associated mutations (including I47V and I50V), which we generated by in vitro selection. This crystal structure showed a unique curling conformation at the flap regions that was not found in the previously reported ligand-free protease structures. Molecular dynamics simulations indicated that the curled flap conformation altered the flap dynamics. These results suggest that the preference for a unique flap conformation influences DRV binding. These results provide new structural insights into elucidating the molecular mechanism of DRV resistance and aid to develop PIs effective against DRV-resistant viruses.

Characteristics of Transmitted Drug-Resistant HIV-1 in Recently Infected Treatment-Naive Patients in Japan.

OBJECTIVES: Progress in antiretroviral treatment has led to fewer virological failure cases, but 10%-20% of treatment-naive HIV/AIDS cases are reported to harbor drug-resistant strains, suggesting transmission of drug-resistant HIV. We aimed to determine the trend in prevalence of transmitted drug-resistant (TDR) HIV in Japan, particularly in recently infected patients. METHODS: Drug-resistance test was performed on 3904 HIV-1-infected cases newly diagnosed between 2007 and 2012. The number of cases infected within 6 months [recent seroconverters (RS)] was estimated by BED assay of 2700 plasma samples. Characteristics of RS cases were further analyzed. RESULTS: The overall prevalence of TDR was 9.1%, ranging from 7.3% in 2008% to 12.5% in 2010. Among 1403 subtype B/E/D cases with >50 CD4 T cell counts and >1000 HIV copies per milliliter, 468 (33.4%) were estimated to be RS. The prevalence of RS was significantly higher among cases who were male, Japanese, and men who have sex with men. The prevalence of TDR did not differ significantly between recent and long-term seroconverters (8.5% vs. 9.2%, respectively, P = 0.68). Common mutations in both groups were M46I/L and T215 revertants. Furthermore, sequences with these mutations, K103N and D30N/N88D formed clusters on phylogenetic trees. CONCLUSION: Our study clarified an increase in prevalence of TDR in Japan from 2007 to 2012. The phylogenetic clustering of cases with M46I/L or T215 revertants suggests that HIV with these mutations have become circulating strains. Furthermore, detailed analyses showed that Japanese men who have sex with men are more aware of their risk of HIV infection.

HIV-1 CRF01_AE and Subtype B Transmission Networks Crossover: A New AE/B Recombinant Identified in Japan.

The major circulating HIV-1 strains in Japan have been subtype B (B) followed by CRF01_AE (AE) in newly diagnosed HIV/AIDS cases. These two subtypes have distinct epidemiological characteristics; B predominates in men who have sex with men, while AE is observed mostly in heterosexuals engaging in high-risk sex. However, transmission networks of these two high-risk populations appear to be crossing over and diffusing. Here we report the emergence of previously unidentified HIV-1 AE/B recombinants in Japan. We initially identified 13 cases with discordant subtyping results with AE (gag MA)/B (pol PR-RT)/AE (env C2V3) by molecular phylogenetic analysis of 1,070 cases who visited Nagoya Medical Center from 1997 to 2012. Genetic characterization of full-length sequences demonstrated that they shared an identical recombinant structure, and was designated as CRF69_01B by the Los Alamos HIV National Laboratory. By reviewing gag, pol, and env sequences collected in the Japanese Drug Resistance HIV-1 Surveillance Network, we found five other CRF69_01B probable cases from different areas in Japan, suggesting that the strain is transmitted widely throughout the country. The time of the most recent common ancestor analyses estimated that CRF69_01B emerged between 1991 and 1995, soon after AE was introduced from neighboring countries in the mid-1990s. Understanding the current epidemic strains is important for the diagnosis and treatment of HIV/AIDS, as well as for the development of globally effective HIV vaccines.

Quasispecies Analyses of the HIV-1 Near-full-length Genome With Illumina MiSeq.

Human immunodeficiency virus type-1 (HIV-1) exhibits high between-host genetic diversity and within-host heterogeneity, recognized as quasispecies. Because HIV-1 quasispecies fluctuate in terms of multiple factors, such as antiretroviral exposure and host immunity, analyzing the HIV-1 genome is critical for selecting effective antiretroviral therapy and understanding within-host viral coevolution mechanisms. Here, to obtain HIV-1 genome sequence information that includes minority variants, we sought to develop a method for evaluating quasispecies throughout the HIV-1 near-full-length genome using the Illumina MiSeq benchtop deep sequencer. To ensure the reliability of minority mutation detection, we applied an analysis method of sequence read mapping onto a consensus sequence derived from de novo assembly followed by iterative mapping and subsequent unique error correction. Deep sequencing analyses of aHIV-1 clone showed that the analysis method reduced erroneous base prevalence below 1% in each sequence position and discarded only < 1% of all collected nucleotides, maximizing the usage of the collected genome sequences. Further, we designed primer sets to amplify the HIV-1 near-full-length genome from clinical plasma samples. Deep sequencing of 92 samples in combination with the primer sets and our analysis method provided sufficient coverage to identify >1%-frequency sequences throughout the genome. When we evaluated sequences of pol genes from 18 treatment-naïve patients' samples, the deep sequencing results were in agreement with Sanger sequencing and identified numerous additional minority mutations. The results suggest that our deep sequencing method would be suitable for identifying within-host viral population dynamics throughout the genome.

Longitudinal Detection and Persistence of Minority Drug-Resistant Populations and Their Effect on Salvage Therapy.

BACKGROUND: Drug-resistant HIV are more prevalent and persist longer than previously demonstrated by bulk sequencing due to the ability to detect low-frequency variants. To clarify a clinical benefit to monitoring minority-level drug resistance populations as a guide to select active drugs for salvage therapy, we retrospectively analyzed the dynamics of low-frequency drug-resistant population in antiretroviral (ARV)-exposed drug resistant individuals. MATERIALS AND METHODS: Six HIV-infected individuals treated with ARV for more than five years were analyzed. These individuals had difficulty in controlling viremia, and treatment regimens were switched multiple times guided by standard drug resistance testing using bulk sequencing. To detect minority variant populations with drug resistance, we used a highly sensitive allele-specific PCR (AS-PCR) with detection thresholds of 0.3-2%. According to ARV used in these individuals, we focused on the following seven reverse transcriptase inhibitor-resistant mutations: M41L, K65R, K70R, K103N, Y181C, M184V, and T215F/Y. Results of AS-PCR were compared with bulk sequencing data for concordance and presence of additional mutations. To clarify the genetic relationship between low-frequency and high-frequency populations, AS-PCR amplicon sequences were compared with bulk sequences in phylogenetic analysis. RESULTS: The use of AS-PCR enabled detection of the drug-resistant mutations, M41L, K103N, Y181C, M184V and T215Y, present as low-frequency populations in five of the six individuals. These drug resistant variants persisted for several years without ARV pressure. Phylogenetic analysis indicated that pre-existing K103N and T215I variants had close genetic relationships with high-frequency K103N and T215I observed during treatment. DISCUSSION AND CONCLUSION: Our results demonstrate the long-term persistence of drug-resistant viruses in the absence of drug pressure. The rapid virologic failures with pre-existing mutant viruses detectable by AS-PCR highlight the clinical importance of low-frequency drug-resistant viruses. Thus, our results highlight the usefulness of AS-PCR and support its expanded evaluation in ART clinical management.

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.

Contact sensitivity in patients with recalcitrant atopic dermatitis.

Patients with atopic dermatitis are usually responsive to conventional treatment such as topical steroids; however, they are sometimes refractory to the treatment. The influence of contact sensitivities on the course of patients with recalcitrant atopic dermatitis is not known. The aim of this study was to investigate whether contact sensitivities affect the course of patients with recalcitrant atopic dermatitis. We evaluated 45 patients with atopic dermatitis who had failed conventional therapy. Patch testing was performed with the Japanese standard series, metal series and/or suspected items. A total of 15 patients had a positive patch test reaction to at least one allergen. The most common allergens were nickel, topical drugs and rubber accelerators. Avoidance of products or food containing allergic substances greatly or partially improved skin symptoms in nine patients. These results suggest that contact allergens and metals may be critical factors causing eczematous lesions in patients with recalcitrant atopic dermatitis.

Japanese external quality assessment program to standardize HIV-1 drug-resistance testing (JEQS2010 program) using in vitro transcribed RNA as reference material.

To design appropriate antiretroviral therapy regimens and avoid the emergence of human immunodeficiency virus (HIV)-1 variants with reduced susceptibility to antiretroviral drugs, genotypic drug-resistance testing (HIV genotyping) is strongly recommended. To monitor the quality of HIV genotyping in Japan, we performed an external quality assessment (EQA), named the Japanese external quality assessment program, to standardize HIV genotyping (JEQS). To accurately evaluate the quality of HIV genotyping, we employed as reference material (RM) a well-characterized sample, in vitro transcribed RNA (trRNA) that includes the HIV gag-pol sequence, and created a JEQS2010 panel consisting of three single variant and three mixed trRNA samples. All 11 participating laboratories showed high concordance rates (>96%) for the single variant samples. Eight laboratories also showed good rates of detecting minor variants, but three laboratories failed to detect the variants comprising one-half of the sample. These three laboratories used a common primer that had four internal mismatches to the minor trRNA clone. This program showed the usefulness of trRNA as RM, the high quality of HIV genotyping, and extensive interlaboratory variation in the ability to detect minor variants. These results suggest that improving the quality of HIV genotyping in Japan requires regularly implementing the EQA program and improving the HIV genotyping protocol in each laboratory.

Phylodynamic analysis reveals CRF01_AE dissemination between Japan and neighboring Asian countries and the role of intravenous drug use in transmission.

BACKGROUND: One major circulating HIV-1 subtype in Southeast Asian countries is CRF01_AE, but little is known about its epidemiology in Japan. We conducted a molecular phylodynamic study of patients newly diagnosed with CRF01_AE from 2003 to 2010. METHODS: Plasma samples from patients registered in Japanese Drug Resistance HIV-1 Surveillance Network were analyzed for protease-reverse transcriptase sequences; all sequences undergo subtyping and phylogenetic analysis using distance-matrix-based, maximum likelihood and Bayesian coalescent Markov Chain Monte Carlo (MCMC) phylogenetic inferences. Transmission clusters were identified using interior branch test and depth-first searches for sub-tree partitions. Times of most recent common ancestor (tMRCAs) of significant clusters were estimated using Bayesian MCMC analysis. RESULTS: Among 3618 patient registered in our network, 243 were infected with CRF01_AE. The majority of individuals with CRF01_AE were Japanese, predominantly male, and reported heterosexual contact as their risk factor. We found 5 large clusters with ≥5 members and 25 small clusters consisting of pairs of individuals with highly related CRF01_AE strains. The earliest cluster showed a tMRCA of 1996, and consisted of individuals with their known risk as heterosexual contacts. The other four large clusters showed later tMRCAs between 2000 and 2002 with members including intravenous drug users (IVDU) and non-Japanese, but not men who have sex with men (MSM). In contrast, small clusters included a high frequency of individuals reporting MSM risk factors. Phylogenetic analysis also showed that some individuals infected with HIV strains spread in East and South-eastern Asian countries. CONCLUSIONS: Introduction of CRF01_AE viruses into Japan is estimated to have occurred in the 1990s. CFR01_AE spread via heterosexual behavior, then among persons connected with non-Japanese, IVDU, and MSM. Phylogenetic analysis demonstrated that some viral variants are largely restricted to Japan, while others have a broad geographic distribution.