Pol Sequences Research Articles

Minimal Cross-resistance to Tenofovir in Children and Adolescents Failing ART Makes Them Eligible for Tenofovir-Lamivudine-Dolutegravir Treatment.

Background:Fixed-dose combination of dolutegravir (DTG) with tenofovir disoproxil fumarate (TDF) and lamivudine (3TC) likely improves adherence and has a favorable resistance profile. We evaluated predicted efficacy of TLD (TDF-3TC-DTG) in children and adolescents failing abacavir (ABC), zidovudine (AZT), or TDF containing regimens.Methods:Drug resistance mutations were analyzed in a retrospective dataset of individuals <19 years of age, failing ABC (n = 293) AZT (n = 288) or TDF (n = 69) based treatment. Pol sequences were submitted to Stanford HIVdb v8.9. Genotypic susceptibility scores were calculated for various DTG-containing regimens.Results:Genotypes were assessed for 650 individuals with a median age of 14 years (IQR 10-17 years). More individuals failed a protease inhibitor (PI)-based (78.3%) than a non-nucleoside reverse transcriptase inhibitors (NNRTI)-based (21.7%) regimen. Most individuals in the AZT group (n = 288; 94.4%) failed a PI-based regimen, compared with 71.0% and 64.2% in the TDF (n = 69) and ABC group (n = 293). Genotypic sensitivity scores <2 to TLD were observed in 8.5% and 9.4% of ABC- and AZT-exposed individuals, compared with 23.2% in the TDF group. The M184V mutation was most often detected in the ABC group (70.6%) versus 60.0% and 52.4% in TDF and AZT groups. The presence of K65R was rare (n = 13, 2.0%) and reduced TLD susceptibility was commonly caused by accumulation of nucleoside reverse transcriptase inhibitor (NRTI) mutations.Conclusions:Cross-resistance to TDF was limited, further reducing concerns about use of transition to TLD in children and adolescents. The NADIA trial has subsequently shown that patients failing a TDF/3TC/EFV regimen can safely be transitioned to a TLD regimen but we do not have data for patients failing an ABC/3TC/NNRTI or PI regimens. Frequent virological monitoring is recommended after switch to DTG, especially in children continuing ABC in the backbone. Clinical studies correlating predicted resistance with clinical outcomes, especially in settings without access to genotyping, are required.

The Molecular Epidemiological and Immunological Characteristics of HIV-1 CRF01_AE/B Recombinants in Nanjing, China.

Human immunodeficiency virus-type 1 (HIV-1) CRF01_AE/B recombinants are newly emerging strains that are spreading rapidly in Southern and Eastern China. This study aimed to elucidate the molecular epidemiological characteristics of HIV-1 CRF01_AE/B recombinants in Nanjing and to explore the impact of these novel strains on the immunological status. A total of 1,013 blood samples from newly diagnosed HIV-1-infected patients were collected in Nanjing from 2015 to 2019, among which 958 partial Pol sequences were sequenced successfully. We depicted the molecular epidemiological characteristics of CRF01_AE/B recombinants by the molecular evolutionary analysis, Bayesian system evolution analysis, and transmission network analysis. The generalized additive mixed model was applied to evaluate the CD4+ T-cell count change of CRF01_AE/B recombinants. The Kaplan–Meier analysis was performed to assess the time from combined antiretroviral therapy (cART) initiation to immune reconstruction. We have identified 102 CRF01_AE/B recombinants (102/958, 10.65%) in Nanjing, including CRF67_01B (45/102, 44.12%), CRF68_01B (35/102, 34.31%), and CRF55_01B (22/102, 12.57%). According to the Bayesian phylogenetic inference, CRF55_01B had a rapid decline stage during 2017–2019, while CRF67_01B and CRF68_01B have experienced a fast growth phase during 2014–2015 and then remained stable. We have constructed 83 transmission networks, in which three larger clusters were composed of CRF67_01B and CRF68_01B. CRF01_AE/B recombinants manifested a faster decrease rate of CD4+ T-cell count than CRF_07BC but similar to CRF01_AE. The probability of achieving immune reconstruction in CRF01_AE/B recombinants was lower than CRF07_BC in the subgroup of baseline CD4+ T-cell count at cART initiation <300 cells/μl. In summary, CRF67_01B and CRF68_01B were the major strains of CRF01_AE/B recombinants in Nanjing, which have formed large transmission clusters between Nanjing and other provinces. CRF01_AE/B recombinants might be associated with rapid disease progression and poor immune reconstruction. The continuous epidemiological monitoring of CRF01_AE/B recombinants should be highly emphasized.

Genetic Diversity and Low Therapeutic Impact of Variant-Specific Markers in HIV-1 Pol Proteins.

The emergence and spread of new HIV-1 variants pose a challenge for the effectiveness of antiretrovirals (ARV) targeting Pol proteins. During viral evolution, non-synonymous mutations have fixed along the viral genome, leading to amino acid (aa) changes that can be variant-specific (V-markers). Those V-markers fixed in positions associated with drug resistance mutations (DRM), or R-markers, can impact drug susceptibility and resistance pathways. All available HIV-1 Pol sequences from ARV-naïve subjects were downloaded from the United States Los Alamos HIV Sequence Database, selecting 59,733 protease (PR), 6,437 retrotranscriptase (RT), and 6,059 integrase (IN) complete sequences ascribed to the four HIV-1 groups and group M subtypes and circulating recombinant forms (CRFs). Using a bioinformatics tool developed in our laboratory (EpiMolBio), we inferred the consensus sequences for each Pol protein and HIV-1 variant to analyze the aa conservation in Pol. We analyzed the Wu–Kabat protein variability coefficient (WK) in PR, RT, and IN group M to study the susceptibility of each site to evolutionary replacements. We identified as V-markers the variant-specific aa changes present in >75% of the sequences in variants with >5 available sequences, considering R-markers those V-markers that corresponded to DRM according to the IAS-USA2019 and Stanford-Database 9.0. The mean aa conservation of HIV-1 and group M consensus was 82.60%/93.11% in PR, 88.81%/94.07% in RT, and 90.98%/96.02% in IN. The median group M WK was 10 in PR, 4 in RT, and 5 in IN. The residues involved in binding or catalytic sites showed a variability <0.5%. We identified 106 V-markers: 31 in PR, 28 in RT, and 47 in IN, present in 11, 12, and 13 variants, respectively. Among them, eight (7.5%) were R-markers, present in five variants, being minor DRM with little potential effect on ARV susceptibility. We present a thorough analysis of Pol variability among all HIV-1 variants circulating to date. The relatively high aa conservation observed in Pol proteins across HIV-1 variants highlights their critical role in the viral cycle. However, further studies are needed to understand the V-markers’ impact on the Pol proteins structure, viral cycle, or treatment strategies, and periodic variability surveillance studies are also required to understand PR, RT, and IN evolution.

Evidence of bovine immunodeficiency virus: A molecular survey in water buffalo populations of Iran.

BackgroundBovine immunodeficiency virus (BIV) is a member of the Retroviridae family causing a progressive lifelong infection in cattle and buffaloes.ObjectiveDespite the worldwide distribution of the virus, the studies concerning the prevalence of BIV in buffalo populations have not been conducted in Iran as yet.MethodsThe BIV proviral DNA was surveyed in 120 whole blood samples of water buffaloes in southwestern Iran. Nested PCR was employed to amplify a 298‐bp fragment of the pol gene. The BIV Pol sequence was detected in 9.1% of the samples. Among PCR‐positive samples, two amplified fragments were confirmed by nucleotide sequencing.Results and conclusionsThe studied sequences were completely identical to each other and had more than 98%–99% nucleotide homology to R‐29 and HXB3 sequences previously deposited in GenBank. Some point mutations that caused coding substitutions were observed in the studied isolates, compared to other strains. A phylogenetic tree was generated based on the BIV Pol nucleotide sequences reported from other countries. All the BIV strains originated from a unique main cluster and then separated from each other over time. This is the first report on the molecular detection of BIV infections in water buffalo populations in Iran. The wide distribution of BIV in different countries including Iran indicates the importance of the infection as it relates to animal health. Although buffaloes show greater resistance to diseases, they should be considered a health risk to cattle. Furthermore, BIV has negative effects on buffalo milk production and can predispose them to secondary infections. Hence, the findings of this study can advance our understanding of the occurrence of BIV infection in Iran, which can play an important role in the distribution of the disease worldwide.

Phylodynamics of HIV in the Mexico City Metropolitan Region.

Evolutionary analyses of viral sequences can provide insights into transmission dynamics, which in turn can optimize prevention interventions. Here, we characterized the dynamics of HIV transmission within the Mexico City metropolitan area. HIV pol sequences from persons recently diagnosed at the largest HIV clinic in Mexico City (between 2016 and 2021) were annotated with demographic/geographic metadata. A multistep phylogenetic approach was applied to identify putative transmission clades. A data set of publicly available sequences was used to assess international introductions. Clades were analyzed with a discrete phylogeographic model to evaluate the timing and intensity of HIV introductions and transmission dynamics among municipalities in the region. A total of 6,802 sequences across 96 municipalities (5,192 from Mexico City and 1,610 from the neighboring State of Mexico) were included (93.6% cisgender men, 5.0% cisgender women, and 1.3% transgender women); 3,971 of these sequences formed 1,206 clusters, involving 78 municipalities, including 89 clusters of ≥10 sequences. Discrete phylogeographic analysis revealed (i) 1,032 viral introductions into the region, over one-half of which were from the United States, and (ii) 354 migration events between municipalities with high support (adjusted Bayes factor of ≥3). The most frequent viral migrations occurred between northern municipalities within Mexico City, i.e., Cuauhtémoc to Iztapalapa (5.2% of events), Iztapalapa to Gustavo A. Madero (5.4%), and Gustavo A. Madero to Cuauhtémoc (6.5%). Our analysis illustrates the complexity of HIV transmission within the Mexico City metropolitan area but also identifies a spatially active transmission area involving a few municipalities in the north of the city, where targeted interventions could have a more pronounced effect on the entire regional epidemic. IMPORTANCE Phylogeographic investigation of the Mexico City HIV epidemic illustrates the complexity of HIV transmission in the region. An active transmission area involving a few municipalities in the north of the city, with transmission links throughout the region, is identified and could be a location where targeted interventions could have a more pronounced effect on the entire regional epidemic, compared with those dispersed in other manners.

Using molecular network analysis to explore the characteristics of HIV-1 transmission in a China-Myanmar border area.

BackgroundThe China-Myanmar border area is considered a hot spot of active HIV-1 recombination in Southeast Asia. To better understand the characteristics of HIV-1 transmission in this area, a cross-sectional HIV-1 molecular epidemiological survey was conducted in Baoshan Prefecture of Yunnan Province.MethodsIn total, 708 newly reported HIV-1 cases in Baoshan Prefecture from 2019 to 2020 were included in this study. HIV-1 gag, pol and env genes were sequenced, and the spatial and demographic distributions of HIV-1 genotypes were analyzed. The characteristics of HIV-1 transmission were investigated using the HIV-1 molecular network method.ResultsIn the 497 samples with genotyping results, 19 HIV-1 genotypes were found, with URFs being the predominant strains (30.2%, 150/497). The main circulating HIV-1 strains were mostly distributed in the northern area of Baoshan. URFs were more likely identified in Burmese individuals, intravenous drug users and those younger than 50 years old. CRF08_BC was more likely detected in farmers and those of Han ethnicity, CRF01_AE in the young and those of Han ethnicity, and CRF07_BC in the subpopulation with junior middle school education and higher. Moreover, CRF118_BC and CRF64_BC were more likely found in the subpopulation aged ≥40 years and ≥50 years, respectively. Among 480 individuals with pol sequence detection, 179 (37.3%) were grouped into 78 clusters, with Baoshan natives being more likely to be in the network. The proportion of the linked individuals showed significant differences when stratified by the regional origin, marital status, age and county of case reporting. In the molecular network, recent infections were more likely to occur among nonfarmers and individuals aged below 30 years.ConclusionsHIV-1 genetics has become complex in Baoshan. HIV-1 molecular network analysis provided transmission characteristics in the local area, and these findings provided information to prioritize transmission-reduction interventions.

Recent HIV infections: evaluation of a simple identification score for newly diagnosed patients.

ABSTRACTOBJECTIVE Recognize incident infection to better characterize the groups that fuel HIV epidemic. We propose a simple score to identify recent infections among newly diagnosed patients as a HIV surveillance tool.METHODS Newly diagnosed patients were defined as recent infections when a negative serological test in the previous year was available. Laboratory tests, such as the avidity index (Bio-Rad, according to the CEPHIA protocol), chemiluminescent intensity (CMIA, architect, Abbott), and the nucleotide ambiguity index of partial pol sequences were used as proxies of recency. A simple score based on clinical symptoms of acute retroviral syndrome during the previous year, CD4+ T cell count, and viral load at admission was tested to assess the predictive power, using receiver operating characteristic (ROC) curves, to identify recent cases of infection.RESULTS We evaluated 204 recently diagnosed patients who were admitted to the Ambulatório de Referência em Moléstias Infecciosas de Santo André (Santo André Reference Infectious Diseases Outpatient Clinic), in the metropolitan region of São Paulo, Brazil, recruited between 2011 and 2018. An HIV-negative test in the year prior to enrollment was documented in 37% of participants. The proportion of cases classified as recent infections (less than one year), according to the laboratory proxies were: 37% (67/181) for an avidity index < 40%, 22% (30/137) for a CMIA < 200, and 68% (124/181) for an ambiguity index < 0.5%. Using different combinations of recency definitions, our score showed an area under the ROC curve from 0.66 to 0.87 to predict recency.CONCLUSIONS Using data from patients’ interviews and routine laboratory tests at admission, a simple score may provide information on HIV recency and thus, a proxy for HIV incidence to guide public policies. This simple for the Brazilian public health system and other low- and middle-income countries.

Frequently Transmission and Close Relationship Among Immigrants in the China-Myanmar Border Region Indicated by Molecular Transmission Analysis From a Cross-Sectional Data.

BackgroundAccurate identification of molecular transmission clusters (MTCs) and understanding the dynamics of human immunodeficiency virus (HIV) transmission are necessary to develop targeted interventions to prevent HIV transmission. We evaluated the characteristics of antiretroviral therapy-naïve individuals who belonged to HIV-1 MTCs in the China–Myanmar border region to inform targeted effective HIV intervention.MethodsPhylogenetic analyses were undertaken on HIV-1 pol sequences to characterize subtypes or circulating recombinant forms and identify MTCs. MTCs were defined as those with 2 or more sequences having bootstrap support > 80% and a pairwise gene distance less than or equal to 0.03. Factors correlated with MTCs were evaluated using logistic regression analysis. The chi-square test was used to compare differences between Chinese and Burmese participants belonging to MTCs.ResultsA total of 900 people had their pol gene successfully sequenced. Twenty-one MTCs were identified and included 110 individuals (12.2%). Individuals in MTCs were more likely to be Burmese [aOR = 2.24 (95% CI: 1.33, 3.79), P = 0.003], be younger [aOR = 0.34 (95% CI: 0.20, 0.58), P < 0.001 for age 26–50 vs. 25 years or younger], have a lower CD4 T cell count [aOR = 2.86 (95% CI: 1.34, 6.11), P = 0.007 for < 200 vs. 350 or greater], and have subtypes CRF07_BC or C [CRF07_BC: aOR = 7.88 (95% CI: 3.55, 17.52), P < 0.001; C: aOR = 2.38 (95% CI: 1.23, 4.62), P = 0.010 compared to CRF01_AE]. In MTCs, Burmese were younger (89.7 vs. 57.7% for age 25 years or younger), had a lower education level (41.0 vs. 8.5% for illiterate), were more likely to be infected through injection drug use (35.9 vs. 12.7%), and had a higher proportion of subtype BC (33.3 vs. 15.5%) and CRF01_AE (20.5 vs. 8.5%) compared to Chinese (P < 0.05 for all).ConclusionBurmese participants were more likely to belong to MTCs, and most MTCs had both Burmese and Chinese participants. These data highlight the bidirectional transmission of HIV-1 frequently transmission and close relationship among immigrants in the China–Myanmar border region. Local health departments should pay more attention to HIV screening and intervention to immigrants Burmese with the characteristics of younger age, having lower CD4 T cell count and infected with HIV subtypes CRF07_ BC or C.

Transmission Clusters, Predominantly Associated With Men Who Have Sex With Men, Play a Main Role in the Propagation of HIV-1 in Northern Spain (2013-2018).

Viruses of HIV-1-infected individuals whose transmission is related group phylogenetically in transmission clusters (TCs). The study of the phylogenetic relations of these viruses and the factors associated with these individuals is essential to analyze the HIV-1 epidemic. In this study, we examine the role of TCs in the epidemiology of HIV-1 infection in Galicia and the Basque County, two regions of northern Spain. A total of 1,158 HIV-1-infected patients from both regions with new diagnoses (NDs) in 2013–2018 were included in the study. Partial HIV-1 pol sequences were analyzed phylogenetically by approximately maximum-likelihood with FastTree 2. In this analysis, 10,687 additional sequences from samples from HIV-1-infected individuals collected in Spain in 1999–2019 were also included to assign TC membership and to determine TCs’ sizes. TCs were defined as those which included viruses from ≥4 individuals, at least 50% of them Spaniards, and with ≥0.95 Shimodaira-Hasegawa-like node support in the phylogenetic tree. Factors associated to TCs were evaluated using odds ratios (OR) and their 95% CI. Fifty-one percent of NDs grouped in 162 TCs. Male patients (OR: 2.6; 95% CI: 1.5–4.7) and men having sex with men (MSM; OR: 2.1; 95% CI: 1.4–3.2) had higher odds of belonging to a TC compared to female and heterosexual patients, respectively. Individuals from Latin America (OR: 0.3; 95% CI: 0.2–0.4), North Africa (OR: 0.4; 95% CI: 0.2–1.0), and especially Sub-Saharan Africa (OR: 0.02; 95% CI: 0.003–0.2) were inversely associated to belonging to TCs compared to native Spaniards. Our results show that TCs are important components of the HIV-1 epidemics in the two Spanish regions studied, where transmission between MSM is predominant. The majority of migrants were infected with viruses not belonging to TCs that expand in Spain. Molecular epidemiology is essential to identify local peculiarities of HIV-1 propagation. The early detection of TCs and prevention of their expansion, implementing effective control measures, could reduce HIV-1 infections.

First report of fibropapillomatosis (FP) and Chelonid alphaherpesvirus 5 (ChHV5) in a green sea turtle (Chelonia mydas) from the historically fibropapillomatosis-free Fernando de Noronha Archipelago, Northeastern Brazil

Fibropapillomatosis (FP) is an infectious disease caused by Chelonid alphaherpesvirus 5 (ChHV5). Nevertheless, its clinical manifestations are considered multifactorial. Due to its relevance, FP is currently monitored in sea turtle populations in the United States, Australia, Caribbean, and Brazil. Between 2000 and 2020, the TAMAR Project/ TAMAR Project Foundation analyzed the prevalence of FP in nine states and oceanic islands along the Brazilian coast, including Fernando de Noronha Archipelago (FNA), a historically FP-free area. A total of 4,435 green sea turtles (Chelonia mydas) were monitored from 2010 to 2016. Additionally, in 2012 and 2014, 43 FP-free skin samples were analyzed for ChHV5 using a qualitative PCR for the UL30 polymerase (pol) sequence. In 2015, a bilateral ocular nodule characterized as an FP tumor was reported in one of the monitored individuals undergoing rehabilitation. Tissue samples were collected following surgical removal of the tumor. Characterization of a 454 bp UL30 polymerase gene revealed a ChHV5 sequence previously reported in other areas of the Atlantic Brazilian coast. In the years following this finding from January 2017 to March 2020, a total of 360 C. mydas were monitored in the same area and no FP tumors were detected. This is the first report of FP and the first detection of ChHV5 in FNA, a finding of great concern considering this site’s historical absence of FP occurrence. This study highlights the importance of monitoring this disease in historically FP-free areas of the Brazilian Atlantic coast.

Significant Upregulation of HERV-K (HML-2) Transcription Levels in Human Lung Cancer and Cancer Cells

Lung cancer is the second most common cancer worldwide and the leading cause of cancer death in the world. Therefore, there is an urgent need to develop new and effective biomarkers for diagnosis and treatment. Under this circumstance, human endogenous retroviruses (HERVs) were recently introduced as novel biomarkers for cancer diagnosis. This study focused on the correlation between lung cancer and HERV-K (HML-2) transcription levels. At the cellular level, different types of lung cancer cells and human normal lung epithelial cells were used to analyze the transcription levels of the HERV-K (HML-2) gag, pol, and env genes by RT–qPCR. At the level of lung cancer patients, blood samples with background information from 734 lung cancer patients and 96 healthy persons were collected to analyze the transcription levels of HERV-K (HML-2) gag, pol, and env genes. The results showed that the transcriptional levels of the HERV-K (HML-2) gag, pol, and env genes in lung cancer cells and lung cancer patient blood samples were significantly higher than those in the healthy controls, which was also verified by RNAScope ISH technology. In addition, we also found that there was a correlation between the abnormal transcription levels of HERV-K (HML-2) genes in lung cancer patients and the clinicopathological parameters of lung cancer. We also identified the distribution locations of the gag, pol, and env primer sequences on each chromosome and analyzed the function of these loci. In conclusion, HERV-K (HML-2) genes may be a potential biomarker for the diagnosis of lung cancer.

Phylogeographic Assessment Reveals Geographic Sources of HIV-1 Dissemination Among Men Who Have Sex With Men in Kenya.

HIV-1 transmission dynamics involving men who have sex with men (MSM) in Africa are not well understood. We investigated the rates of HIV-1 transmission between MSM across three regions in Kenya: Coast, Nairobi, and Nyanza. We analyzed 372 HIV-1 partial pol sequences sampled during 2006–2019 from MSM in Coast (N = 178, 47.9%), Nairobi (N = 137, 36.8%), and Nyanza (N = 57, 15.3%) provinces in Kenya. Maximum-likelihood (ML) phylogenetics and Bayesian inference were used to determine HIV-1 clusters, evolutionary dynamics, and virus migration rates between geographic regions. HIV-1 sub-subtype A1 (72.0%) was most common followed by subtype D (11.0%), unique recombinant forms (8.9%), subtype C (5.9%), CRF 21A2D (0.8%), subtype G (0.8%), CRF 16A2D (0.3%), and subtype B (0.3%). Forty-six clusters (size range 2–20 sequences) were found—half (50.0%) of which had evidence of extensive HIV-1 mixing among different provinces. Data revealed an exponential increase in infections among MSM during the early-to-mid 2000s and stable or decreasing transmission dynamics in recent years (2017–2019). Phylogeographic inference showed significant (Bayes factor, BF > 3) HIV-1 dissemination from Coast to Nairobi and Nyanza provinces, and from Nairobi to Nyanza province. Strengthening HIV-1 prevention programs to MSM in geographic locations with higher HIV-1 prevalence among MSM (such as Coast and Nairobi) may reduce HIV-1 incidence among MSM in Kenya.

CRF07_BC is associated with slow HIV disease progression in Chinese patients

HIV subtypes convey important epidemiological information and possibly influence the rate of disease progression. In this study, HIV disease progression in patients infected with CRF01_AE, CRF07_BC, and subtype B was compared in the largest HIV molecular epidemiology study ever done in China. A national data set of HIV pol sequences was assembled by pooling sequences from public databases and the Beijing HIV laboratory network. Logistic regression was used to assess factors associated with the risk of AIDS at diagnosis ([AIDSAD], defined as a CD4 count < 200 cells/µL) in patients with HIV subtype B, CRF01_AE, and CRF07_BC. Of the 20,663 sequences, 9,156 (44.3%) were CRF01_AE. CRF07_BC was responsible for 28.3% of infections, followed by B (13.9%). In multivariable analysis, the risk of AIDSAD differed significantly according to HIV subtype (OR for CRF07_BC vs. B: 0.46, 95% CI 0.39─0.53), age (OR for ≥ 65 years vs. < 18 years: 4.3 95% CI 1.81─11.8), and transmission risk groups (OR for men who have sex with men vs. heterosexuals: 0.67 95% CI 0.6─0.75). These findings suggest that HIV diversity in China is constantly evolving and gaining in complexity. CRF07_BC is less pathogenic than subtype B, while CRF01_AE is as pathogenic as B.

Dissemination Dynamics of HIV-1 Subtype B Pandemic and Non-pandemic Lineages Circulating in Amazonas, Brazil.

The HIV-1 epidemic in the Amazonas state, as in most of Brazil, is dominated by subtype B. The state, nonetheless, is singular for its significant co-circulation of the variants BCAR, which can mostly be found in the Caribbean region, and BPAN, a clade that emerged in the United States and aggregates almost the totality of subtype B infections world-wide. The Amazonian HIV-1 epidemic provides a unique scenario to compare the epidemic potential of BPAN and BCAR clades spreading in the same population. To reconstruct the spatiotemporal dynamic and demographic history of both subtype B lineages circulating in Amazonas, we analyzed 1,272 HIV-1 pol sequences sampled in that state between 2009 and 2018. Our phylogeographic analyses revealed that while most BCAR infections resulted from a single successful founder event that took place in the Amazonas state around the late 1970s, most BPAN infections resulted from the expansion of multiple clusters seeded in the state since the late 1980s. Our data support the existence of at least four large clusters of the pandemic form in Amazonas, two of them nested in Brazil’s largest known subtype B cluster (BBR–I), and two others resulting from new introductions detected here. The reconstruction of the demographic history of the most prevalent BPAN (n = 4) and BCAR (n = 1) clades identified in Amazonas revealed that all clades displayed a continuous expansion [effective reproductive number (Re) > 1] until most recent times. During the period of co-circulation from the late 1990s onward, the Re of Amazonian BPAN and BCAR clusters behaved quite alike, fluctuating between 2.0 and 3.0. These findings support that the BCAR and BPAN variants circulating in the Brazilian state of Amazonas displayed different evolutionary histories, but similar epidemic trajectories and transmissibility over the last two decades, which is consistent with the notion that both subtype B variants display comparable epidemic potential. Our findings also revealed that despite significant advances in the treatment of HIV infections in the Amazonas state, BCAR and BPAN variants continue to expand and show no signs of the epidemic stabilization observed in other parts of the country.

Quantifying rates of HIV-1 flow between risk groups and geographic locations in Kenya: A country-wide phylogenetic study.

In Kenya, HIV-1 key populations including men having sex with men (MSM), people who inject drugs (PWID) and female sex workers (FSW) are thought to significantly contribute to HIV-1 transmission in the wider, mostly heterosexual (HET) HIV-1 transmission network. However, clear data on HIV-1 transmission dynamics within and between these groups are limited. We aimed to empirically quantify rates of HIV-1 flow between key populations and the HET population, as well as between different geographic regions to determine HIV-1 ‘hotspots’ and their contribution to HIV-1 transmission in Kenya. We used maximum-likelihood phylogenetic and Bayesian inference to analyse 4058 HIV-1 pol sequences (representing 0.3 per cent of the epidemic in Kenya) sampled 1986–2019 from individuals of different risk groups and regions in Kenya. We found 89 per cent within-risk group transmission and 11 per cent mixing between risk groups, cyclic HIV-1 exchange between adjoining geographic provinces and strong evidence of HIV-1 dissemination from (i) West-to-East (i.e. higher-to-lower HIV-1 prevalence regions), and (ii) heterosexual-to-key populations. Low HIV-1 prevalence regions and key populations are sinks rather than major sources of HIV-1 transmission in Kenya. Targeting key populations in Kenya needs to occur concurrently with strengthening interventions in the general epidemic.

HIV-1 drug resistance mutations among individuals with low-level viraemia while taking combination ART in Botswana.

To assess whether a single instance of low-level viraemia (LLV) is associated with the presence of drug resistance mutations (DRMs) and predicts subsequent virological failure (VF) in adults receiving ART in 30 communities participating in the Botswana Combination Prevention Project. A total of 6078 HIV-1 C pol sequences were generated and analysed using the Stanford HIV drug resistance database. LLV was defined as plasma VL = 51-999 copies/mL and VF was defined as plasma VL ≥ 1000 copies/mL. Among 6078 people with HIV (PWH), 4443 (73%) were on ART for at least 6 months. Of the 332 persons on ART with VL > 50 copies/mL, 175 (4%) had VL ≥ 1000 copies/mL and 157 (4%) had LLV at baseline. The prevalence of any DRM was 57 (36%) and 78 (45%) in persons with LLV and VL ≥ 1000 copies/mL, respectively. Major DRMs were found in 31 (20%) with LLV and 53 (30%) with VL ≥ 1000 copies/mL (P = 0.04). Among the 135 PWH with at least one DRM, 17% had NRTI-, 35% NNRTI-, 6% PI- and 3% INSTI-associated mutations. Among the 3596 participants who were followed up, 1709 (48%) were on ART for ≥6 months at entry and had at least one subsequent VL measurement (median 29 months), 43 (3%) of whom had LLV. The OR of experiencing VF in persons with LLV at entry was 36-fold higher than in the virally suppressed group. A single LLV measurement while on ART strongly predicted the risk of future VF, suggesting the use of VL > 50 copies/mL as an indication for more intensive adherence support with more frequent VL monitoring.

Pretreatment drug resistance in people living with HIV: A large retrospective cohort study in Chongqing, China.

ObjectivesThe emergence of pretreatment drug resistance (PDR) caused by increased usage of antiretroviral therapy (ART) represents a significant challenge to HIV management. In this study, we evaluated the prevalence of PDR in people living with HIV (PLWH) in Chongqing, China.MethodsWe retrospectively collected the data of 1110 ART‐naïve PLWH in Chongqing from January 1, 2018 to June 30, 2021. HIV‐1 genotypes and drug resistance were analyzed using the HIV‐1 pol sequence. Risk factors associated with PDR were evaluated via the logistic regression model.ResultsNine genotypes were detected among 1110 participants, with CRF07_BC (55.68%) being the dominant genotype, followed by CRF01_AE (21.44%), CRF08_BC (14.14%), and other genotypes (8.74%). Of all the participants, 24.14% exhibited drug resistance mutations (DRMs). The predominant DRMs for non‐nucleoside reverse transcriptase inhibitors (NNRTIs) and nucleoside reverse transcriptase inhibitors (NRTIs) were V179D/E/A/DIN (13.60%) and M184V/I (1.44%), respectively, whereas only two major DRMs (M46L and I54L) were identified for protease inhibitors (PIs). The total prevalence of PDR was 10.54%, with 2.43%, 7.66%, and 1.71% participants exhibiting PDR to NRTIs, NNRTIs, and PIs, respectively. Furthermore, female PLWH, delays in ART initiation, and the CRF08_BC genotype were associated with a higher risk of PDR.ConclusionsOur study provides the first large cohort data on the prevalence of PDR in Chongqing, China. HIV‐1 genotypes are diverse and complex, with a moderate level of PDR, which does not reach the threshold for the initiation of a public health response. Nevertheless, continuous surveillance of PDR is both useful and advisable.

Drug Resistance to HIV-1 Integrase Inhibitors Among Treatment-Naive Patients in Beijing, China.

IntroductionIntegrase strand transfer inhibitors (INSTIs) are important drugs that are currently used as the first line treatment for HIV-1 patients. The aim of this study was to characterize HIV-1 INSTI mutations among ART-naive patients in Beijing from 2019–2021.Methods865 ART-naive patients were enrolled in this study between January 2019 and June 2021 in Beijing. The amplification of the entire pol gene containing the reverse transcriptase, protease and integrase regions was performed using a validated In-house SBS method. HIV-1 subtypes and circulating recombinant forms (CRFs) were determined using the COMET online tool (http://comet.retrovirology.lu). Stanford HIV-1 drug resistance database (HIVdb version 8.9) was used to analyze the mutations.Results865 HIV-1 pol sequences were successfully amplified and sequenced. Among them, no major INSTI-related mutations were identified, but 12 polymorphic accessory mutations were found. Two patients have E138A and G163R mutations respectively and both could cause low-level resistance to RAL and EVG. Furthermore, one patient having S230R mutation resulted in low-level resistance to RAL, EVG, DTG and BIC.ConclusionThe prevalence of INSTIs mutations remains low, which demonstrated that INSTIs have good applicability currently in our city. Nevertheless, it is very important to monitor the INSTI-related mutations in Beijing.

Molecular Epidemiology and Transmission Dynamics of the HIV-1 Epidemic in Ethiopia: Epidemic Decline Coincided With Behavioral Interventions Before ART Scale-Up.

BackgroundEthiopia is one of the sub-Saharan countries hit hard by the HIV epidemic. Previous studies have shown that subtype C dominates the Ethiopian HIV-1 epidemic, but the evolutionary and temporal dynamics of HIV-1 in Ethiopia have not been closely scrutinized. Understanding the evolutionary and epidemiological pattern of HIV is vital to monitor the spread, evaluate and implement HIV prevention strategies.MethodsWe analyzed 1,276 Ethiopian HIV-1 subtype C polymerase (pol sequences), including 144 newly generated sequences, collected from different parts of the country from 1986 to 2017. We employed state-of-art maximum likelihood and Bayesian phylodynamic analyses to comprehensively describe the evolutionary dynamics of the HIV-1 epidemic in Ethiopia. We used Bayesian phylodynamic models to estimate the dynamics of the effective population size (Ne) and reproductive numbers (Re) through time for the HIV epidemic in Ethiopia.ResultsOur analysis revealed that the Ethiopian HIV-1 epidemic originated from two independent introductions at the beginning of the 1970s and 1980s from eastern and southern African countries, respectively, followed by epidemic growth reaching its maximum in the early 1990s. We identified three large clusters with a majority of Ethiopian sequences. Phylodynamic analyses revealed that all three clusters were characterized by high transmission rates during the early epidemic, followed by a decline in HIV-1 transmissions after 1990. Re was high (4–6) during the earlier time of the epidemic but dropped significantly and remained low (Re < 1) after the mid-1990. Similarly, with an expected shift in time, the effective population size (Ne) steadily increased until the beginning of 2000, followed by a decline and stabilization until recent years. The phylodynamic analyses corroborated the modeled UNAIDS incidence and prevalence estimates.ConclusionThe rapid decline in the HIV epidemic took place a decade before introducing antiretroviral therapy in Ethiopia and coincided with early behavioral, preventive, and awareness interventions implemented in the country. Our findings highlight the importance of behavioral interventions and antiretroviral therapy scale-up to halt and maintain HIV transmissions at low levels (Re < 1). The phylodynamic analyses provide epidemiological insights not directly available using standard surveillance and may inform the adjustment of public health strategies in HIV prevention in Ethiopia.

DRUG RESISTANCE MUTATIONS AND TRANSMISSION CLUSTERS OF THE HIV-1 CRF01_AE SUB-EPIDEMIC IN BULGARIA

Background. In Bulgaria the second most predominant HIV-1 strain after subtype B is CRF01_AE. Material and Methods. 270 HIV-1 polymerase (pol) CRF01_AE sequences collected between 1995–2019 were analyzed with the phylogenetic approach. To identify transmission clusters, we used two different genetic distance thresholds (d), i. e. 1.5% and 0.5%. Results. Genetic distance d=1.5% defined more distant infections in a huge cluster of 154 sequences composed mostly of people, who inject drugs (PWID), whereas when (d) was reduced to 0.5%, we determined more recent transmissions and the large cluster disintegrated into smaller ones. Different drug resistance mutations (DRM) were detected in patient’s HIV-1 pol sequences and were most common in male heterosexual (HET) single sequences. Conclusions. Our data showed repeated introduction of CRF01_AE in Bulgaria and rapid spread of the infection among PWID groups. Molecular monitoring of the epidemic among PWID communities could help reduce the spread of HIV-1 infection.