Transmitted drug resistance and molecular transmission network among treatment-naive HIV-1 patients in Wenzhou, China, 2020-2023.

BACKGROUND: Transmitted drug resistance (TDR) increases the risk of antiretroviral therapy (ART) failure in HIV-1 patients. This study investigated the molecular epidemiology of TDR and its transmission networks among newly diagnosed HIV-1 patients in Wenzhou, China. METHODS: We enrolled 1878 ART-naive HIV-1 patients from January 2020 to October 2023. TDR was evaluated using the Stanford University HIV Drug Resistance Database. We performed phylogenetic analysis, genotyping, transmission clustering, and population-based TDR-related factor analysis. RESULTS: Among 1782 patients with successful genotyping, TDR prevalence was 5.7%. Multivariable analysis identified CRF08_BC subtype (adjusted odds ratio [aOR] 18.59, 95% CI 3.79-336.18, p = 0.004), CD4 > 500 cells/mm³ (aOR 2.19, 95% CI 1.16-4.03, p = 0.013), and year 2023 (aOR 1.83, 95% CI 1.11-4.89, p = 0.039) as factors associated with higher TDR risk. The most prevalent NNRTI mutations were K103N, E138A, and V179E. Seven TDR transmission clusters were identified, notably one with V179D that expanded during 2020-2023. CONCLUSIONS: While TDR prevalence in Wenzhou remained lower than in other Chinese regions, an upward trend was observed. Most resistant individuals were in transmission clusters, predominantly middle-aged and elderly. NNRTI resistance was severe and concentrated in efavirenz, nevirapine, and rilpivirine. Enhanced HIV surveillance and wider free antiretroviral options are crucial to control drug-resistant HIV spread in Wenzhou.

Combining molecular transmission network analysis and spatial epidemiology to reveal HIV-1 transmission pattern among the older people in Nanjing, China.

BACKGROUND: In China, the problem of HIV infection among the older people has become increasingly prominent. This study aimed to analyze the pattern and influencing factors of HIV transmission based on a genomic and spatial epidemiological analysis among this population. METHODS: A total of 432 older people who were aged ≥ 50 years, newly diagnosed with HIV-1 between January 2018 and December 2021 and without a history of ART were enrolled. HIV-1 pol gene sequence was obtained by viral RNA extraction and nested PCR. The molecular transmission network was constructed using HIV-TRACE and the spatial distribution analyses were performed in ArcGIS. The multivariate logistic regression analysis was performed to analyze the factors associated with clustering. RESULTS: A total of 382 sequences were successfully sequenced, of which CRF07_BC (52.3%), CRF01_AE (32.5%), and CRF08_BC (6.8%) were the main HIV-1 strains. A total of 176 sequences entered the molecular network, with a clustering rate of 46.1%. Impressively, the clustering rate among older people infected through commercial heterosexual contact was as high as 61.7% and three female sex workers (FSWs) were observed in the network. The individuals who were aged ≥ 60 years and transmitted the virus by commercial heterosexual contact had a higher clustering rate, while those who were retirees or engaged other occupations and with higher education degree were less likely to cluster. There was a positive spatial correlation of clustering rate (Global Moran I = 0.206, P < 0.001) at the town level and the highly aggregated regions were mainly distributed in rural area. We determined three large clusters which mainly spread in the intra-region of certain towns in rural areas. Notably, 54.5% of cases in large clusters were transmitted through commercial heterosexual contact. CONCLUSIONS: Our joint analysis of molecular and spatial epidemiology effectively revealed the spatial aggregation of HIV transmission and highlighted that towns of high aggregation were mainly located in rural area. Also, we found vital role of commercial heterosexual contact in HIV transmission among older people. Therefore, health resources should be directed towards highly aggregated rural areas and prevention strategy should take critical persons as entry points.

Increasing Frequency and Transmission of HIV-1 Non-B Subtypes Among Men Who Have Sex With Men in the Swiss HIV Cohort Study.

BACKGROUND: In Switzerland, HIV-1 transmission among men who have sex with men (MSM) has been dominated by subtype B, whilst non-B subtypes are commonly attributed to infections acquired abroad among heterosexuals. Here, we evaluated the temporal trends of non-B subtypes and the characteristics of molecular transmission clusters (MTCs) among MSM. METHODS: Sociodemographic and clinical data and partial pol sequences were obtained from participants enrolled in the Swiss HIV Cohort Study. For non-B subtypes, maximum likelihood trees were constructed, from which Swiss MTCs were identified and analyzed by transmission group. RESULTS: Non-B subtypes were identified in 8.1% (416/5116) of MSM participants. CRF01_AE was the most prevalent strain (3.5%), followed by subtype A (1.2%), F (1.1%), CRF02_AG (1.1%), C (0.9%), and G (0.3%). Between 1990 and 2019, an increase in the proportion of newly diagnosed individuals (0/123 [0%] to 11/32 [34%]) with non-B subtypes in MSM was found. Across all non-B subtypes, the majority of MSM MTCs were European. Larger MTCs were observed for MSM than heterosexuals. CONCLUSIONS: We found a substantial increase in HIV-1 non-B subtypes among MSM in Switzerland and the occurrence of large MTCs, highlighting the importance of molecular surveillance in guiding public health strategies targeting the HIV-1 epidemic.

HAPHPIPE: Haplotype Reconstruction and Phylodynamics for Deep Sequencing of Intrahost Viral Populations.

Deep sequencing of viral populations using next-generation sequencing (NGS) offers opportunities to understand and investigate evolution, transmission dynamics, and population genetics. Currently, the standard practice for processing NGS data to study viral populations is to summarize all the observed sequences from a sample as a single consensus sequence, thus discarding valuable information about the intrahost viral molecular epidemiology. Furthermore, existing analytical pipelines may only analyze genomic regions involved in drug resistance, thus are not suited for full viral genome analysis. Here, we present HAPHPIPE, a HAplotype and PHylodynamics PIPEline for genome-wide assembly of viral consensus sequences and haplotypes. The HAPHPIPE protocol includes modules for quality trimming, error correction, de novo assembly, alignment, and haplotype reconstruction. The resulting consensus sequences, haplotypes, and alignments can be further analyzed using a variety of phylogenetic and population genetic software. HAPHPIPE is designed to provide users with a single pipeline to rapidly analyze sequences from viral populations generated from NGS platforms and provide quality output properly formatted for downstream evolutionary analyses.

Transmitted drug resistance and transmission clusters among HIV-1 treatment-naïve patients in Guangdong, China: a cross-sectional study.

BACKGROUND: Transmitted drug resistance (TDR) that affects the effectiveness of the first-line antiretroviral therapy (ART) regimen is becoming prevalent worldwide. However, its prevalence and transmission among HIV-1 treatment-naïve patients in Guangdong, China are rarely reported. We aimed to comprehensively analyze the prevalence of TDR and the transmission clusters of HIV-1 infected persons before ART in Guangdong. METHODS: The HIV-1 treatment-naïve patients were recruited between January 2018 and December 2018. The HIV-1 pol region was amplified by reverse transcriptional PCR and sequenced by sanger sequencing. Genotypes, surveillance drug resistance mutations (SDRMs) and TDR were analyzed. Genetic transmission clusters among patients were identified by pairwise Tamura-Nei 93 genetic distance, with a threshold of 0.015. RESULTS: A total of 2368 (97.17%) HIV-1 pol sequences were successfully amplified and sequenced from the enrolled 2437 patients. CRF07_BC (35.90%, 850/2368), CRF01_AE (35.56%, 842/2368) and CRF55_01B (10.30%, 244/2368) were the main HIV-1 genotypes circulating in Guangdong. Twenty-one SDRMs were identified among fifty-two drug-resistant sequences. The overall prevalence of TDR was 2.20% (52/2368). Among the 2368 patients who underwent sequencing, 8 (0.34%) had TDR to protease inhibitors (PIs), 22 (0.93%) to nucleoside reverse transcriptase inhibitors (NRTIs), and 23 (0.97%) to non-nucleoside reverse transcriptase inhibitors (NNRTIs). Two (0.08%) sequences showed dual-class resistance to both NRTIs and NNRTIs, and no sequences showed triple-class resistance. A total of 1066 (45.02%) sequences were segregated into 194 clusters, ranging from 2 to 414 sequences. In total, 15 (28.85%) of patients with TDR were included in 9 clusters; one cluster contained two TDR sequences with the K103N mutation was observed. CONCLUSIONS: There is high HIV-1 genetic heterogeneity among patients in Guangdong. Although the overall prevalence of TDR is low, it is still necessary to remain vigilant regarding some important SDRMs.

Reconstruction of the Genetic History and the Current Spread of HIV-1 Subtype A in Germany.

HIV-1 non-B infections have been increasing in Europe for several years. In Germany, subtype A belongs to the most abundant non-B subtypes showing an increasing prevalence of 8.3% among new infections in 2016. Here we trace the origin and examine the current spread of the German HIV-1 subtype A epidemic. Bayesian coalescence and birth-death analyses were performed with 180 German HIV-1 pol sequences and 528 related and publicly available sequences to reconstruct the population dynamics and fluctuations for each of the transmission groups. Our reconstructions indicate two distinct sources of the German subtype A epidemic, with an Eastern European and an Eastern African lineage both cocirculating in the country. A total of 13 German-origin clusters were identified; among these, 6 clusters showed recent activity. Introductions leading to further countrywide spread originated predominantly from Eastern Africa when introduced before 2005. Since 2005, however, spreading introductions have occurred exclusively within the Eastern European clade. Moreover, we observed changes in the main route of subtype A transmission. The beginning of the German epidemic (1985 to 1995) was dominated by heterosexual transmission of the Eastern African lineage. Since 2005, transmissions among German men who have sex with men (MSM) have been increasing and have been associated with the Eastern European lineage. Infections among people who inject drugs dominated between 1998 and 2005. Our findings on HIV-1 subtype A infections provide new insights into the spread of this virus and extend the understanding of the HIV epidemic in Germany.IMPORTANCE HIV-1 subtype A is the second most prevalent subtype worldwide, with a high prevalence in Eastern Africa and Eastern Europe. However, an increase of non-B infections, including subtype A infections, has been observed in Germany and other European countries. There has simultaneously been an increased flow of refugees into Europe and especially into Germany, raising the question of whether the surge in non-B infections resulted from this increased immigration or whether German transmission chains are mainly involved. This study is the first comprehensive subtype A study from a western European country analyzing in detail its phylogenetic origin, the impact of various transmission routes, and its current spread. The results and conclusions presented provide new and substantial insights for virologists, epidemiologists, and the general public health sector. In this regard, they should be useful to those authorities responsible for developing public health intervention strategies to combat the further spread of HIV/AIDS.

Transmission dynamics and evolutionary history of 2019-nCoV.

To investigate the time origin, genetic diversity, and transmission dynamics of the recent 2019-nCoV outbreak in China and beyond, a total of 32 genomes of virus strains sampled from China, Thailand, and the USA with sampling dates between 24 December 2019 and 23 January 2020 were analyzed. Phylogenetic, transmission network, and likelihood-mapping analyses of the genome sequences were performed. On the basis of the likelihood-mapping analysis, the increasing tree-like signals (from 0% to 8.2%, 18.2%, and 25.4%) over time may be indicative of increasing genetic diversity of 2019-nCoV in human hosts. We identified three phylogenetic clusters using the Bayesian inference framework and three transmission clusters using transmission network analysis, with only one cluster identified by both methods using the above genome sequences of 2019-nCoV strains. The estimated mean evolutionary rate for 2019-nCoV ranged from 1.7926 × 10-3 to 1.8266 × 10-3 substitutions per site per year. On the basis of our study, undertaking epidemiological investigations and genomic data surveillance could positively impact public health in terms of guiding prevention efforts to reduce 2019-nCOV transmission in real-time.

Determining the Origins of Human Immunodeficiency Virus Type 1 Drug-resistant Minority Variants in People Who Are Recently Infected Using Phylogenetic Reconstruction.

BACKGROUND: Drug-resistant minority variants (DRMinVs) detected in patients who recently acquired human immunodeficiency virus type 1 (HIV-1) can be transmitted, generated de novo through virus replication, or technical errors. The first form is likely to persist and result in treatment failure, while the latter two could be stochastic and transient. METHODS: Ultradeep sequencing of plasma samples from 835 individuals with recent HIV-1 infection in the United Kingdom was performed to detect DRMinVs at a mutation frequency between 2% and 20%. Sequence alignments including >110 000 HIV-1 partial pol consensus sequences from the UK HIV Drug Resistance Database (UK-HDRD), linked to epidemiological and clinical data from the HIV and AIDS Reporting System, were used for transmission cluster analysis. Transmission clusters were identified using Cluster Picker with a clade support of >90% and maximum genetic distances of 4.5% or 1.5%, the latter to limit detection to likely direct transmission events. RESULTS: Drug-resistant majority variants (DRMajVs) were detected in 66 (7.9%) and DRMinVs in 84 (10.1%) of the recently infected individuals. High levels of clustering to sequences in UK-HDRD were observed for both DRMajV (n = 48; 72.7%) and DRMinV (n = 63; 75.0%) sequences. Of these, 43 (65.2%) with DRMajVs were in a transmission cluster with sequences that harbored the same DR mutation compared to only 3 (3.6%) sequences with DRMinVs (P < .00001, Fisher exact test). Evidence of likely direct transmission of DRMajVs was observed for 25/66 (37.9%), whereas none were observed for the DRMinVs (P < .00001). CONCLUSIONS: Using a densely sampled HIV-infected population, we show no evidence of DRMinV transmission among recently infected individuals.

HIV co-infection is associated with increased transmission risk in patients with chronic hepatitis C virus.

Molecular epidemiological analysis of viral pathogens can identify factors associated with increased transmission risk. We investigated the frequency of genetic clustering in a large data set of NS34A, NS5A, and NS5B viral sequences from patients with chronic hepatitis C virus (HCV). Within a subset of patients with longitudinal samples, Receiver Operator Characteristic (ROC) analysis was applied which identified a threshold of 0.02 substitutions/site as most appropriate for clustering. From the 7457 patients with chronic HCV infection included in this analysis, we inferred 256 clusters comprising 541 patients (7.3%). We found that HCV/HIV co-infection, young age, and high HCV viral load were all associated with increased clustering frequency, an indicator of increased transmission risk. In light of previous work on HCV/HIV co-infection in acute HCV cohorts, our results suggest that patients with HCV/HIV co-infection may disproportionately be the source of new HCV infections and treatment efforts should be geared towards viral elimination in this vulnerable population.

AUTO-TUNE: SELECTING THE DISTANCE THRESHOLD FOR INFERRING HIV TRANSMISSION CLUSTERS.

Molecular surveillance of viral pathogens and inference of transmission networks from genomic data play an increasingly important role in public health efforts, especially for HIV-1. For many methods, the genetic distance threshold used to connect sequences in the transmission network is a key parameter informing the properties of inferred networks. Using a distance threshold that is too high can result in a network with many spurious links, making it difficult to interpret. Conversely, a distance threshold that is too low can result in a network with too few links, which may not capture key insights into clusters of public health concern. Published research using the HIV-TRACE software package frequently uses the default threshold of 0.015 substitutions/site for HIV pol gene sequences, but in many cases, investigators heuristically select other threshold parameters to better capture the underlying dynamics of the epidemic they are studying. Here, we present a general heuristic scoring approach for tuning a distance threshold adaptively, which seeks to prevent the formation of giant clusters. We prioritize the ratio of the sizes of the largest and the second largest cluster, maximizing the number of clusters present in the network. We apply our scoring heuristic to outbreaks with different characteristics, such as regional or temporal variability, and demonstrate the utility of using the scoring mechanism's suggested distance threshold to identify clusters exhibiting risk factors that would have otherwise been more difficult to identify. For example, while we found that a 0.015 substitutions/site distance threshold is typical for US-like epidemics, recent outbreaks like the CRF07_BC subtype among men who have sex with men (MSM) in China have been found to have a lower optimal threshold of 0.005 to better capture the transition from injected drug use (IDU) to MSM as the primary risk factor. Alternatively, in communities surrounding Lake Victoria in Uganda, where there has been sustained hetero-sexual transmission for many years, we found that a larger distance threshold is necessary to capture a more risk factor-diverse population with sparse sampling over a longer period of time. Such identification may allow for more informed intervention action by respective public health officials.

SARS-CoV-2 population dynamics in immunocompetent individuals in a closed transmission chain shows genomic diversity over the course of infection.

BACKGROUND: SARS-CoV-2 remains rapidly evolving, and many biologically important genomic substitutions/indels have characterised novel SARS-CoV-2 lineages, which have emerged during successive global waves of the pandemic. Worldwide genomic sequencing has been able to monitor these waves, track transmission clusters, and examine viral evolution in real time to help inform healthcare policy. One school of thought is that an apparent greater than average divergence in an emerging lineage from contemporary variants may require persistent infection, for example in an immunocompromised host. Due to the nature of the COVID-19 pandemic and sampling, there were few studies that examined the evolutionary trajectory of SARS-CoV-2 in healthy individuals. METHODS: We investigated viral evolutionary trends and participant symptomatology within a cluster of 16 SARS-CoV-2 infected, immunocompetent individuals with no co-morbidities in a closed transmission chain. Longitudinal nasopharyngeal swab sampling allowed characterisation of SARS-CoV-2 intra-host variation over time at both the dominant and minor genomic variant levels through Nimagen-Illumina sequencing. RESULTS: A change in viral lineage assignment was observed in individual infections; however, there was only one indel and no evidence of recombination over the period of an acute infection. Minor and dominant genomic modifications varied between participants, with some minor genomic modifications increasing in abundance to become the dominant viral sequence during infection. CONCLUSIONS: Data from this cohort of SARS-CoV-2-infected participants demonstrated that long-term persistent infection in an immunocompromised host was not necessarily a prerequisite for generating a greater than average frequency of amino acid substitutions. Amino acid substitutions at both the dominant and minor genomic sequence level were observed in immunocompetent individuals during infection showing that viral lineage changes can occur generating viral diversity.

AUTO-TUNE: selecting the distance threshold for inferring HIV transmission clusters.

Molecular surveillance of viral pathogens and inference of transmission networks from genomic data play an increasingly important role in public health efforts, especially for HIV-1. For many methods, the genetic distance threshold used to connect sequences in the transmission network is a key parameter informing the properties of inferred networks. Using a distance threshold that is too high can result in a network with many spurious links, making it difficult to interpret. Conversely, a distance threshold that is too low can result in a network with too few links, which may not capture key insights into clusters of public health concern. Published research using the HIV-TRACE software package frequently uses the default threshold of 0.015 substitutions/site for HIV pol gene sequences, but in many cases, investigators heuristically select other threshold parameters to better capture the underlying dynamics of the epidemic they are studying. Here, we present a general heuristic scoring approach for tuning a distance threshold adaptively, which seeks to prevent the formation of giant clusters. We prioritize the ratio of the sizes of the largest and the second largest cluster, maximizing the number of clusters present in the network. We apply our scoring heuristic to outbreaks with different characteristics, such as regional or temporal variability, and demonstrate the utility of using the scoring mechanism's suggested distance threshold to identify clusters exhibiting risk factors that would have otherwise been more difficult to identify. For example, while we found that a 0.015 substitutions/site distance threshold is typical for US-like epidemics, recent outbreaks like the CRF07_BC subtype among men who have sex with men (MSM) in China have been found to have a lower optimal threshold of 0.005 to better capture the transition from injected drug use (IDU) to MSM as the primary risk factor. Alternatively, in communities surrounding Lake Victoria in Uganda, where there has been sustained heterosexual transmission for many years, we found that a larger distance threshold is necessary to capture a more risk factor-diverse population with sparse sampling over a longer period of time. Such identification may allow for more informed intervention action by respective public health officials.

Molecular cluster, transmission characteristics, origin and dynamics analysis of HIV-1 CRF59_01B in China: A molecular epidemiology study.

PURPOSE: This study investigated for the HIV-1 CRF59_01B epidemic's spatiotemporal dynamics and its transmission networks in China. METHODS: Between 2007 and 2020, a total of 250 partial pol gene sequences of HIV-1 CRF59_01B were collected from four regions (10 Chinese provinces). Phylogenetic tree construction and cluster identification were then performed. The Bayesian skyline and birth-death susceptible-infected-removed models were employed for the phylodynamic analyses of subtypes and large clusters, respectively. Phylogenetic analyses and trait diffusion of these sequences were performed using Bayesian phylogenetic methods (beast-classic package). Distance-based molecular network analyses were performed to identify putative relationships. RESULTS: Using a genetic distance threshold of 1.3 %, We identified 45 clusters that included 62.40 % (156/250) of the sequences. Three clusters (6.67 %, 3/45) had 10 or more sequences, and were considered "large clusters". Six clusters (13.33 %) included sequences from different regions (Southeast, Northeast, Southeast, and Central China). Thirteen clusters (28.89 %) included sequences of men who had sex with men only, three clusters (6.67 %) included sequences of heterosexuals only, and 12 clusters (26.67 %) included sequences of both groups. The substitution rate of CRF59_01B was 1.91 × 10-3 substitutions per site per year [95 % highest posterior density (HPD) interval: 1.39 × 10-3-2.49 × 10-3)], the time to the most recent common ancestor of CRF59_01B was to be 1992.83 (95 % HPD: 1977.97-2002.81). A Bayesian skyline plot revealed that the effective population size of CRF59_01B increased from 2000 to 2015 and remained stable after 2015. The large clusters showed continuous growth from 2013 to 2020. Phylogeographic analysis showed that CRF59_01B B most likely originated in Southeast China, with a posterior probability of 97.44 %, and then spread to other regions. CONCLUSIONS: Our study revealed the temporal and geographical origins of HIV-1 CRF59_01B as well as the process of transmission among various regions and risk groups in China, which can help develop targeted HIV prevention strategies.

Distinct Clusters of HIV-1 CRF01_AE in Zhejiang, China: High-Risk Transmission Cluster 4 Requires Heightened Surveillance.

Background: HIV-1 CRF01_AE is becoming the predominant HIV-1 subtype among patients in China. The distribution and characteristics of transmission clusters of HIV-1 CRF01_AE in Zhejiang, Eastern China remains unclear. This study analyzed the epidemiologic characteristics and transmission clusters of HIV-1 CRF01_AE in Zhejiang. Methods: Plasma samples obtained from 152 patients of HIV-1 CRF01_AE not undergoing ART were used to amplify HIV-1 pol and env gene. CRF01_AE drug resistance mutations (DRM) prevalence was analysed using Stanford University's HIV Drug Resistance Database. A phylogenetic tree was constructed using FastTree (version 2.1.11) based on the GTR nucleotide substitution model and visualized using Figtree (version 1.4.4) and The Interactive Tree of Life; the Chinese HIV Gene Sequence Data Platform was used to construct genetic transmission networks. Results: Majority samples could be grouped into CRF01_AE transmission Clusters 1 (11.2%), 4 (64.5%), and 5 (7.2%). The CD4+ T-cell counts in Cluster 1, 4a, 4b are lower than 5 were 15, 38, 30, and 248 cells/mm3, respectively (P < 0.05). The high X4 tropism rates were 13.2%, 11.8%, 20.0%, and 0.0% in Clusters 1, 4a, 4b, and 5, respectively. DRM rates in Clusters 4a and 4b were 17.6%, and 25.45% respectively (P < 0.05), whereas they were 17.6% and 18.2% in Clusters 1 and 5, respectively. In total, 24 transmission genetic networks, comprising 72 sequences and 61 links, were discovered; of them, 61.2%, 11.7%, and 18.2% were from Clusters 4, 1, and 5, respectively (P < 0.05). Conclusion: In Zhejiang, different CRF01_AE clusters displayed unique clinic features. Cluster 4, particularly Cluster 4b, was considered a high-risk transmission cluster. The surveillance of epidemiology of HIV-1 should be enhanced to minimize its transmission.

Vaccine-induced reduction of COVID-19 clusters in school settings in Japan during the epidemic wave caused by B.1.1.529 (Omicron) BA.2, 2022.

Clusters of COVID-19 in high-risk settings, such as schools, have been deemed a critical driving force of the major epidemic waves at the societal level. In Japan, the vaccination coverage among students remained low up to early 2022, especially for 5-11-year-olds. The vaccination of the student population only started in February 2022. Given this background and considering that vaccine effectiveness against school transmission has not been intensively studied, this paper proposes a mathematical model that links the occurrence of clustering to the case count among populations aged 0-19, 20-59, and 60+ years of age. We first estimated the protected (immune) fraction of each age group either by infection or vaccination and then linked the case count in each age group to the number of clusters via a time series regression model that accounts for the time-varying hazard of clustering per infector. From January 3 to May 30, 2022, there were 4,722 reported clusters in school settings. Our model suggests that the immunity offered by vaccination averted 226 (95% credible interval: 219-232) school clusters. Counterfactual scenarios assuming elevated vaccination coverage with faster roll-out reveal that additional school clusters could have been averted. Our study indicates that even relatively low vaccination coverage among students could substantially lower the risk of clustering through vaccine-induced immunity. Our results also suggest that antigenically updated vaccines that are more effective against the variant responsible for the ongoing epidemic may greatly help decrease not only the incidence but also the unnecessary loss of learning opportunities among school-age students.

Transmitted drug resistance and transmission clusters among ART-naïve HIV-1-infected individuals from 2019 to 2021 in Nanjing, China.

Background: Transmitted drug resistance (TDR) is an increasingly prevalent problem worldwide, which will significantly compromise the effectiveness of HIV treatments. However, in Nanjing, China, there is still a dearth of research on the prevalence and transmission of TDR among ART-naïve HIV-1-infected individuals. This study aimed to understand the prevalence and transmission of TDR in Nanjing. Methods: A total of 1,393 participants who were newly diagnosed with HIV-1 and had not received ART between January 2019 and December 2021 were enrolled in this study. HIV-1 pol gene sequence was obtained by viral RNA extraction and nested PCR amplification. Genotypes, TDR and transmission cluster analyses were conducted using phylogenetic tree, Stanford HIV database algorithm and HIV-TRACE, respectively. Univariate and multivariate logistic regression analyses were performed to identify the factors associated with TDR. Results: A total of 1,161 sequences were successfully sequenced, of which CRF07_BC (40.6%), CRF01_AE (38.4%) and CRF105_0107 (6.3%) were the main HIV-1 genotypes. The overall prevalence of TDR was 7.8%, with 2.0% to PIs, 1.0% to NRTIs, and 4.8% to NNRTIs. No sequence showed double-class resistance. Multivariate logistic regression analysis revealed that compared with CRF01_AE, subtype B (OR = 2.869, 95%CI: 1.093-7.420) and female (OR = 2.359, 95%CI: 1.182-4.707) were risk factors for TDR. Q58E was the most prevalent detected protease inhibitor (PI) -associated mutation, and V179E was the most frequently detected non-nucleoside reverse transcriptase inhibitor (NNRTI) -associated mutation. A total of 613 (52.8%) sequences were segregated into 137 clusters, ranging from 2 to 74 sequences. Among 44 individuals with TDR (48.4%) within 21 clusters, K103N/KN was the most frequent TDR-associated mutation (31.8%), followed by Q58E/QE (20.5%) and G190A (15.9%). Individuals with the same TDR-associated mutations were usually cross-linked in transmission clusters. Moreover, we identified 9 clusters in which there was a transmission relationship between drug-resistant individuals, and 4 clusters in which drug-resistant cases increased during the study period. Conclusion: The overall prevalence of TDR in Nanjing was at a moderate level during the past 3 years. However, nearly half of TDR individuals were included in the transmission clusters, and some drug-resistant individuals have transmitted in the clusters. Therefore, HIV drug-resistance prevention, monitoring and response efforts should be sustained and expanded to reduce the prevalence and transmission of TDR in Nanjing.

HIV-1 Transmission Cluster in Injection Drug Users in Nairobi City, Kenya.

Background: While there is a striking increase in the prevalence of HIV in injection drug users, information on envelope-gene subtypes and transmission clusters in injection drug users is scarce. Method: In a cross-sectional study, 247 injection drug users were recruited via out-rich method. Deoxyribonucleic acid was extracted from dry blood spot samples, amplified by Polymerase Chain Reaction and sequenced. Subtyping was performed using COntext-based Modeling for Expeditious Typing (COMET) and Recombinant Identification Program (RIP) tools. Phylogenetic diversity and Transmission clusters were identified using MEGA version 6.0 and TreeLink, respectively. Results: Overall, 42 (17.0%) injection drug users were sero-positive for HIV-1. Of the 37 samples successfully sequenced, 29 (78.4%) sequences were identified as A1, 6 (16.2%) as AG while 1 (2.7%) as A1/G/AE and A1/C recombinants. The HIV subtypes formed clusters with little genetic diversity. Conclusion: The high HIV prevalence was associated with transmission clusters and diversity in subtypes indicating ongoing local transmission. Therefore, there is need for comprehensive HIV care tailored to this population.

Exploring the Evolutionary History and Phylodynamics of Human Immunodeficiency Virus Type 1 Outbreak From Unnao, India Using Phylogenetic Approach.

Certain rural and semiurban settings in the Unnao district, Uttar Pradesh, India observed an unprecedented increase in the detection of HIV cases during July 2017. Subsequent investigations through health camps and a follow-up case-control study attributed the outbreak to the unsafe injection exposures during treatment. In this study, we have undertaken a secondary analysis to understand the phylogenetic aspects of the outbreak-associated HIV-1 sequences along with the origin and phylodynamics of these sequences. The initial phylogenetic analysis indicated separate monophyletic grouping and there was no mixing of outbreak-associated sequences with sequences from other parts of India. Transmission network analysis using distance-based and non-distance-based methods revealed the existence of transmission clusters within the monophyletic Unnao clade. The median time to the most recent common ancestor (tMRCA) for sequences from Unnao using the pol gene region was observed to be 2011.87 [95% highest posterior density (HPD): 2010.09-2013.53], while the estimates using envelope (env) gene region sequences traced the tMRCA to 2010.33 (95% HPD: 2007.76-2012.99). Phylodynamics estimates demonstrated that the pace of this local epidemic has slowed down in recent times before the time of sampling, but was certainly on an upward track since its inception till 2014.

Transmission cluster characteristics of global, regional, and lineage-specific SARS-CoV-2 phylogenies.

The SARS-CoV-2 pandemic has been presenting in periodic waves and multiple variants, of which some dominated over time with increased transmissibility. SARS-CoV-2 is still adapting in the human population, thus it is crucial to understand its evolutionary patterns and dynamics ahead of time. In this work, we analyzed transmission clusters and topology of SARS-CoV-2 phylogenies at the global, regional (North America) and clade-specific (Delta and Omicron) epidemic scales. We used the Nextstrain's nCov open global all-time phylogeny (September 2022, 2,698 strains, 2,243 for North America, 499 for Delta21A, and 543 for Omicron20M), with Nextstrain's clade annotation and Pango lineages. Transmission clusters were identified using Phylopart, DYNAMITE, and several tree imbalance measures were calculated, including staircase-ness, Sackin and Colless index. We found that the phylogenetic clustering profiles of the global epidemic have highest diversification at a distance threshold of 3% (divergence of 10, where the tree sampled median is 49). Phylopart and DYNAMITE clusters moderately-to-highly agree with the Pango nomenclature and the Nextstrain's clade. At the regional and clade-specific scale, transmission clustering profiles tend to flatten and similar clusters are found at distance thresholds between 0.05% and 25%. All the considered phylogenies exhibit high tree imbalance with respect to what expected in random phylogenies, suggesting short infection times and antigenic drift, perhaps due to progressive transition from innate to adaptive immunity in the population.

Plasma Virome Reveals Blooms and Transmission of Anellovirus in Intravenous Drug Users with HIV-1, HCV, and/or HBV Infections.

Intravenous drug users (IDUs) are a high-risk group for HIV-1, hepatitis C virus (HCV), and hepatitis B virus (HBV) infections, which are the leading causes of death in IDUs. However, the plasma virome of IDUs and how it is influenced by above viral infections remain unclear. Using viral metagenomics, we determined the plasma virome of IDUs and its association with HIV-1, HCV, and/or HBV infections. Compared with healthy individuals, IDUs especially those with major viral infections had higher viral abundance and diversity. Anelloviridae dominated plasma virome. Coinfections of multiple anelloviruses were common, and anelloviruses from the same genus tended to coexist together. In this study, 4,487 anellovirus ORF1 sequences were identified, including 1,620 (36.1%) with less than 69% identity to any known sequences, which tripled the current number. Compared with healthy controls (HC), more anellovirus sequences were observed in neg-IDUs, and HIV-1, HCV, and/or HBV infections further expanded the sequence number in IDUs, which was characterized by the emergence of novel divergent taxons and blooms of resident anelloviruses. Pegivirus was mainly identified in infected IDUs. Five main pegivirus transmission clusters (TCs) were identified by phylogenetic analysis, suggesting a transmission link. Similar anellovirus profiles were observed in IDUs within the same TC, suggesting transmission of anellome among IDUs. Our data suggested that IDUs suffered higher plasma viral burden especially anelloviruses, which was associated with HIV-1, HCV, and/or HBV infections. Blooms in abundance and unprecedented diversity of anellovirus highlighted active evolution and replication of this virus in blood circulation, and an uncharacterized role it may engage with the host. IMPORTANCE Virome is associated with immune status and determines or influences disease progression through both pathogenic and resident viruses. Increased viral burden in IDUs especially those with major viral infections indicated the suboptimal immune status and high infection risks of these population. Blooms in abundance and unprecedented diversity of anellovirus highlighted its active evolution and replication in the blood circulation, and sensitive response to other viral infections. In addition, transmission cluster analysis revealed the transmission link of pegivirus among IDUs, and the individuals with transmission links shared similar anellome profiles. In-depth monitoring of the plasma virome in high-risk populations is not only needed for surveillance for emerging viruses and transmission networks of major and neglected bloodborne viruses, but also important for a better understanding of commensal viruses and their role it may engage with immune system.