Phylogenetic Analysis of Hepatitis C Virus Infections in a Large Belgian Cohort Using Next-Generation Sequencing of Full-Length Genomes.

The hepatitis C virus (HCV) epidemic in Western countries is primarily perpetuated by the sub-populations of men who have sex with men (MSM) and people who inject drugs (PWID). Understanding the dynamics of transmission in these communities is crucial for removing the remaining hurdles towards HCV elimination. We sequenced 269 annotated HCV plasma samples using probe enrichment and next-generation sequencing, obtaining 224 open reading frames of HCV (OR497849-OR498072). Maximum likelihood phylogenies were generated on the four most prevalent subtypes in this study (HCV1a, 1b, 3a, 4d) with a subsequent transmission cluster analysis. The highest rate of clustering was observed for HCV4d samples (13/17 (76.47%)). The second highest rate of clustering was observed in HCV1a samples (42/78 (53.85%)) with significant association with HIV-positive MSM. HCV1b and HCV3a had very low rates of clustering (2/83 (2.41%) and (0/29)). The spread of the prevalent subtype HCV1b appears to have been largely curtailed, and we demonstrate the onwards transmission of HCV1a and HCV4d in the HIV-positive MSM population across municipal borders. More systematic data collection and sequencing is needed to allow a better understanding of the HCV transmission among the community of PWID and overcome the remaining barriers for HCV elimination in Belgium.

High viral suppression and detection of dolutegravir-resistance associated mutations in treatment-experienced Tanzanian adults living with HIV-1 in Dar es Salaam.

To curb HIV infection rate in Tanzania, antiretroviral therapy (ART) has been scaled up since 2006, and in 2019, the country shifted to regimen including dolutegravir as a default first line. We assessed the success of ART and the contribution of HIV drug resistance (HIVDR) to unsuppressed viral loads. Between February and May 2023 a cross-sectional survey with random sampling was conducted in the six clinics in an urban cohort in Dar es Salaam. Patients with unsuppresed viral loads (local criteria viral load (VL) ≥ 1000 copies/mL) were tested for HIVDR mutations using the WHO adapted protocol for plasma samples. Mutations were interpreted using the Stanford HIVDR database. In total 600 individuals participated in this survey, the majority were female (76.83%), mean age ([Formula: see text] standard deviation) was 44.0 ([Formula: see text] 11.6) years. The median duration on ART (interquartile range) was 6.5 (3.9-10.2) years. Approximately 99% were receiving tenofovir + lamivudine + dolutegravir as a fixed dose combination. VL testing was successful in 99.67% (598/600) of survey patients and only 33 had VL ≥ 1000 copies/mL, resulting in a viral suppression level of 94.48% (565/598, 95% CI 92.34-96.17%). For 23 samples, protease and reverse transcriptase (RT) genotyping were successful, with 13 sequences containing RT inhibitor surveillance drug resistance mutations (SDRMs) (56.5%). No SDRM against protease inhibitors were detected. Thirty samples were successfully genotyped for integrase with 3 sequences (10.08%) containing integrase strand transfer inhibitor (INSTI) SDRMs. In samples successfully genotyped in the three genetic regions, 68.18% (16/22) had a genotypic susceptibility score (GSS) ≥ 2.5 for the concurrent regimen, implying factors beyond drug resistance caused the unsuppressed viral load. For five patients, GSS indicated that HIVDR may have caused the unsuppressed viral load. All three patients with INSTI resistance mutations were highly resistant to dolutegravir and accumulated nucleoside and non-nucleoside RT inhibitor HIVDR mutations. Although in this cohort the last 95 UNAIDS target was almost achieved, HIVDR mutations, including INSTIs resistance mutations were detected in HIV-positive individuals taking ART for at least one year. We recommend the design and implementation of high-impact interventions to prevent the increase of HIVDR, failure of dolutegravir and address the non-resistance factors in the study area.

Recommendations on data sharing in HIV drug resistance research.

• Human immunodeficiency virus (HIV) drug resistance has implications for antiretroviral treatment strategies and for containing the HIV pandemic because the development of HIV drug resistance leads to the requirement for antiretroviral drugs that may be less effective, less well-tolerated, and more expensive than those used in first-line regimens. • HIV drug resistance studies are designed to determine which HIV mutations are selected by antiretroviral drugs and, in turn, how these mutations affect antiretroviral drug susceptibility and response to future antiretroviral treatment regimens. • Such studies collectively form a vital knowledge base essential for monitoring global HIV drug resistance trends, interpreting HIV genotypic tests, and updating HIV treatment guidelines. • Although HIV drug resistance data are collected in many studies, such data are often not publicly shared, prompting the need to recommend best practices to encourage and standardize HIV drug resistance data sharing. • In contrast to other viruses, sharing HIV sequences from phylogenetic studies of transmission dynamics requires additional precautions as HIV transmission is criminalized in many countries and regions. • Our recommendations are designed to ensure that the data that contribute to HIV drug resistance knowledge will be available without undue hardship to those publishing HIV drug resistance studies and without risk to people living with HIV.

HIV-1 subtype diversity and phylogenetic insight into non-B subtype transmission in Slovenia, 1989-2013.

INTRODUCTION: Disease progression, drug resistance mutations, and treatment strategies may vary by HIV-1 subtype. This study determined HIV-1 subtypes circulating in Slovenia, a Central European country with an HIV-1 epidemic driven by men who have sex with men, focusing on molecular epidemiology of non-B subtypes. METHODS: A total of 367 HIV-1 sequences were included. Subtype was assigned by employing eight different HIV subtyping tools coupled with maximum likelihood phylogenetic analyses. RESULTS: The subtyping tools COMET, jpHMM, and REGA 3.0 exhibited the best performance on the dataset studied. Phylogenetic analyses showed a 14.7% prevalence of non-B subtypes, with subtype A detected most frequently (4.9%), followed by CRF02_AG (2.4%), subtype C (1.1%), subtypes D, G, and CRF01_AE (0.8% each), and subtypes F and CRF22_01A1 (0.3% each). A subtype could not be assigned to 12 sequences (3.3%), indicating potential unique recombinant forms. Non-B subtypes were significantly associated with a heterosexual route of transmission and infection acquired in Eastern Europe, Africa, or Asia. CONCLUSION: In a country where subtype B is predominant, non-B subtypes were observed in one out of seven patients, a non-negligible proportion, which underlines the importance of systematic surveillance of HIV subtype diversity and the corresponding molecular epidemiology.

High seroprevalence of Leishmania infantum is linked to immune activation in people with HIV: a two-stage cross-sectional study in Bahia, Brazil.

Visceral leishmaniasis is an opportunistic disease in HIV-1 infected individuals, unrecognized as a determining factor for AIDS diagnosis. The growing geographical overlap of HIV-1 and Leishmania infections is an emerging challenge worldwide, as co-infection increases morbidity and mortality for both infections. Here, we determined the prevalence of people living with HIV (PWH) with a previous or ongoing infection by Leishmania infantum and investigated the virological and immunological factors associated with co-infection. We adopted a two-stage cross-sectional cohort (CSC) design (CSC-I, n = 5,346 and CSC-II, n = 317) of treatment-naïve HIV-1-infected individuals in Bahia, Brazil. In CSC-I, samples collected between 1998 and 2013 were used for serological screening for leishmaniasis by an in-house Enzyme-Linked Immunosorbent Assay (ELISA) with SLA (Soluble Leishmania infantum Antigen), resulting in a prevalence of previous or ongoing infection of 16.27%. Next, 317 PWH were prospectively recruited from July 2014 to December 2015 with the collection of sociodemographic and clinical data. Serological validation by two different immunoassays confirmed a prevalence of 15.46 and 8.20% by anti-SLA, and anti-HSP70 serology, respectively, whereas 4.73% were double-positive (DP). Stratification of these 317 individuals in DP and double-negative (DN) revealed a significant reduction of CD4+ counts and CD4+/CD8+ ratios and a tendency of increased viral load in the DP group, as compared to DN. No statistical differences in HIV-1 subtype distribution were observed between the two groups. However, we found a significant increase of CXCL10 (p = 0.0076) and a tendency of increased CXCL9 (p = 0.061) in individuals with DP serology, demonstrating intensified immune activation in this group. These findings were corroborated at the transcriptome level in independent Leishmania- and HIV-1-infected cohorts (Swiss HIV Cohort and Piaui Northeast Brazil Cohort), indicating that CXCL10 transcripts are shared by the IFN-dominated immune activation gene signatures of both pathogens and positively correlated to viral load in untreated PWH. This study demonstrated a high prevalence of PWH with L. infantum seropositivity in Bahia, Brazil, linked to IFN-mediated immune activation and a significant decrease in CD4+ levels. Our results highlight the urgent need to increase awareness and define public health strategies for the management and prevention of HIV-1 and L. infantum co-infection.

Changes to virus taxonomy and the ICTV Statutes ratified by the International Committee on Taxonomy of Viruses (2023).

This article reports changes to virus taxonomy and taxon nomenclature that were approved and ratified by the International Committee on Taxonomy of Viruses (ICTV) in April 2023. The entire ICTV membership was invited to vote on 174 taxonomic proposals that had been approved by the ICTV Executive Committee in July 2022, as well as a proposed revision of the ICTV Statutes. All proposals and the revised ICTV Statutes were approved by a majority of the voting membership. Of note, the ICTV continued the process of renaming existing species in accordance with the recently mandated binomial format and included gene transfer agents (GTAs) in the classification framework by classifying them as viriforms. In total, one class, seven orders, 31 families, 214 genera, and 858 species were created.

Analyses of Early ZIKV Genomes Are Consistent with Viral Spread from Northeast Brazil to the Americas.

The Americas, particularly Brazil, were greatly impacted by the widespread Zika virus (ZIKV) outbreak in 2015 and 2016. Efforts were made to implement genomic surveillance of ZIKV as part of the public health responses. The accuracy of spatiotemporal reconstructions of the epidemic spread relies on the unbiased sampling of the transmission process. In the early stages of the outbreak, we recruited patients exhibiting clinical symptoms of arbovirus-like infection from Salvador and Campo Formoso, Bahia, in Northeast Brazil. Between May 2015 and June 2016, we identified 21 cases of acute ZIKV infection and subsequently recovered 14 near full-length sequences using the amplicon tiling multiplex approach with nanopore sequencing. We performed a time-calibrated discrete phylogeographic analysis to trace the spread and migration history of the ZIKV. Our phylogenetic analysis supports a consistent relationship between ZIKV migration from Northeast to Southeast Brazil and its subsequent dissemination beyond Brazil. Additionally, our analysis provides insights into the migration of ZIKV from Brazil to Haiti and the role Brazil played in the spread of ZIKV to other countries, such as Singapore, the USA, and the Dominican Republic. The data generated by this study enhances our understanding of ZIKV dynamics and supports the existing knowledge, which can aid in future surveillance efforts against the virus.

Virus taxonomy and the role of the International Committee on Taxonomy of Viruses (ICTV).

The taxonomy of viruses is developed and overseen by the International Committee on Taxonomy of Viruses (ICTV), which scrutinizes, approves and ratifies taxonomic proposals, and maintains a list of virus taxa with approved names (https://ictv.global). The ICTV has approximately 180 members who vote by simple majority. Taxon-specific Study Groups established by the ICTV have a combined membership of over 600 scientists from the wider virology community; they provide comprehensive expertise across the range of known viruses and are major contributors to the creation and evaluation of taxonomic proposals. Proposals can be submitted by anyone and will be considered by the ICTV irrespective of Study Group support. Thus, virus taxonomy is developed from within the virology community and realized by a democratic decision-making process. The ICTV upholds the distinction between a virus or replicating genetic element as a physical entity and the taxon category to which it is assigned. This is reflected by the nomenclature of the virus species taxon, which is now mandated by the ICTV to be in a binomial format (genus + species epithet) and is typographically distinct from the names of viruses. Classification of viruses below the rank of species (such as, genotypes or strains) is not within the remit of the ICTV. This article, authored by the ICTV Executive Committee, explains the principles of virus taxonomy and the organization, function, processes and resources of the ICTV, with the aim of encouraging greater understanding and interaction among the wider virology community.

Four principles to establish a universal virus taxonomy.

A universal taxonomy of viruses is essential for a comprehensive view of the virus world and for communicating the complicated evolutionary relationships among viruses. However, there are major differences in the conceptualisation and approaches to virus classification and nomenclature among virologists, clinicians, agronomists, and other interested parties. Here, we provide recommendations to guide the construction of a coherent and comprehensive virus taxonomy, based on expert scientific consensus. Firstly, assignments of viruses should be congruent with the best attainable reconstruction of their evolutionary histories, i.e., taxa should be monophyletic. This fundamental principle for classification of viruses is currently included in the International Committee on Taxonomy of Viruses (ICTV) code only for the rank of species. Secondly, phenotypic and ecological properties of viruses may inform, but not override, evolutionary relatedness in the placement of ranks. Thirdly, alternative classifications that consider phenotypic attributes, such as being vector-borne (e.g., "arboviruses"), infecting a certain type of host (e.g., "mycoviruses," "bacteriophages") or displaying specific pathogenicity (e.g., "human immunodeficiency viruses"), may serve important clinical and regulatory purposes but often create polyphyletic categories that do not reflect evolutionary relationships. Nevertheless, such classifications ought to be maintained if they serve the needs of specific communities or play a practical clinical or regulatory role. However, they should not be considered or called taxonomies. Finally, while an evolution-based framework enables viruses discovered by metagenomics to be incorporated into the ICTV taxonomy, there are essential requirements for quality control of the sequence data used for these assignments. Combined, these four principles will enable future development and expansion of virus taxonomy as the true evolutionary diversity of viruses becomes apparent.

Application of Phylodynamic Tools to Inform the Public Health Response to COVID-19: Qualitative Analysis of Expert Opinions.

BACKGROUND: In the wake of the SARS-CoV-2 pandemic, scientists have scrambled to collect and analyze SARS-CoV-2 genomic data to inform public health responses to COVID-19 in real time. Open source phylogenetic and data visualization platforms for monitoring SARS-CoV-2 genomic epidemiology have rapidly gained popularity for their ability to illuminate spatial-temporal transmission patterns worldwide. However, the utility of such tools to inform public health decision-making for COVID-19 in real time remains to be explored. OBJECTIVE: The aim of this study is to convene experts in public health, infectious diseases, virology, and bioinformatics-many of whom were actively engaged in the COVID-19 response-to discuss and report on the application of phylodynamic tools to inform pandemic responses. METHODS: In total, 4 focus groups (FGs) occurred between June 2020 and June 2021, covering both the pre- and postvariant strain emergence and vaccination eras of the ongoing COVID-19 crisis. Participants included national and international academic and government researchers, clinicians, public health practitioners, and other stakeholders recruited through purposive and convenience sampling by the study team. Open-ended questions were developed to prompt discussion. FGs I and II concentrated on phylodynamics for the public health practitioner, while FGs III and IV discussed the methodological nuances of phylodynamic inference. Two FGs per topic area to increase data saturation. An iterative, thematic qualitative framework was used for data analysis. RESULTS: We invited 41 experts to the FGs, and 23 (56%) agreed to participate. Across all the FG sessions, 15 (65%) of the participants were female, 17 (74%) were White, and 5 (22%) were Black. Participants were described as molecular epidemiologists (MEs; n=9, 39%), clinician-researchers (n=3, 13%), infectious disease experts (IDs; n=4, 17%), and public health professionals at the local (PHs; n=4, 17%), state (n=2, 9%), and federal (n=1, 4%) levels. They represented multiple countries in Europe, the United States, and the Caribbean. Nine major themes arose from the discussions: (1) translational/implementation science, (2) precision public health, (3) fundamental unknowns, (4) proper scientific communication, (5) methods of epidemiological investigation, (6) sampling bias, (7) interoperability standards, (8) academic/public health partnerships, and (9) resources. Collectively, participants felt that successful uptake of phylodynamic tools to inform the public health response relies on the strength of academic and public health partnerships. They called for interoperability standards in sequence data sharing, urged careful reporting to prevent misinterpretations, imagined that public health responses could be tailored to specific variants, and cited resource issues that would need to be addressed by policy makers in future outbreaks. CONCLUSIONS: This study is the first to detail the viewpoints of public health practitioners and molecular epidemiology experts on the use of viral genomic data to inform the response to the COVID-19 pandemic. The data gathered during this study provide important information from experts to help streamline the functionality and use of phylodynamic tools for pandemic responses.

Analysis of mutational history of multidrug-resistant genotypes with a mutagenetic tree model.

Human immunodeficiency virus (HIV) can develop resistance to all antiretroviral drugs. Multidrug resistance, however, is a rare event in modern HIV treatment, but can be life-threatening, particular in patients with very long therapy histories and in areas with limited access to novel drugs. To understand the evolution of multidrug resistance, we analyzed the EuResist database to uncover the accumulation of mutations over time. We hypothesize that the accumulation of resistance mutations is not acquired simultaneously and randomly across viral genotypes but rather tends to follow a predetermined order. The knowledge of this order might help to elucidate potential mechanisms of multidrug resistance. Our evolutionary model shows an almost monotonic increase of resistance with each acquired mutation, including less well-known nucleoside reverse transcriptase (RT) inhibitor-related mutations like K223Q, L228H, and Q242H. Mutations within the integrase (IN) (T97A, E138A/K G140S, Q148H, N155H) indicate high probability of multidrug resistance. Hence, these IN mutations also tend to be observed together with mutations in the protease (PR) and RT. We followed up with an analysis of the mutation-specific error rates of our model given the data. We identified several mutations with unusual rates (PR: M41L, L33F, IN: G140S). This could imply the existence of previously unknown virus variants in the viral quasispecies. In conclusion, our bioinformatics model supports the analysis and understanding of multidrug resistance.

Recent changes to virus taxonomy ratified by the International Committee on Taxonomy of Viruses (2022).

This article reports the changes to virus taxonomy approved and ratified by the International Committee on Taxonomy of Viruses (ICTV) in March 2022. The entire ICTV was invited to vote on 174 taxonomic proposals approved by the ICTV Executive Committee at its annual meeting in July 2021. All proposals were ratified by an absolute majority of the ICTV members. Of note, the Study Groups have started to implement the new rule for uniform virus species naming that became effective in 2021 and mandates the binomial 'Genus_name species_epithet' format with or without Latinization. As a result of this ratification, the names of 6,481 virus species (more than 60 percent of all species names currently recognized by ICTV) now follow this format.

Full-genome next-generation sequencing of hepatitis C virus to assess the accuracy of genotyping by the commercial assay LiPA and the prevalence of resistance-associated substitutions in a Belgian cohort.

BACKGROUND: Although most currently used regimens for Hepatitis C virus (HCV) infections can be initiated without prior knowledge of genotype and subtype, genotyping is still useful to identify patients who might benefit from a personalized treatment due to resistance to direct-acting antivirals (DAA). OBJECTIVES: To assess the utility of full-genome next-generation sequencing (FG-NGS) for HCV genotyping. STUDY DESIGN: 138 HCV plasma samples previously genotyped by VERSANT HCV Genotype Assay (LiPA) were subjected to FG-NGS and phylogenetically genotyped Genome Detective. Consensuses were analysed by HCV-GLUE for resistance-associated substitutions (RASs) and their impact on treatment response was investigated. RESULTS: 102/138 (73.9%) samples were sequenced to a genome coverage and depth of >90% of the HCV open reading frame covered by >100 reads/site. Concordant genotype and subtype results were assigned in 97.1% and 79.4% of samples, respectively. FG-NGS resolved the subtype of 13.7% samples that had ambiguous calls by LiPA and identified one dual infection and one recombinant strain. At least one RAS was found for the HCV genes NS3, NS5A, and NS5B in 2.91%, 36.98% and 27.3% samples, respectively. Irrespective of the observed RAS, all patients responded well to DAA treatment, except for HCV1b-infected patients treated with Zepatier (33.3% failure rate (5/15)). CONCLUSION: While LiPA and FG-NGS showed overall good concordance, FG-NGS improved specificity for subtypes, recombinant and mixed infections. FG-NGS enabled the detection of RAS, but its predictive value for treatment outcome in DAA-naïve patients remains uncertain. With additional refinements, FG-NGS may be the way forward for HCV genotyping.

Women in the European Virus Bioinformatics Center.

Viruses are the cause of a considerable burden to human, animal and plant health, while on the other hand playing an important role in regulating entire ecosystems. The power of new sequencing technologies combined with new tools for processing "Big Data" offers unprecedented opportunities to answer fundamental questions in virology. Virologists have an urgent need for virus-specific bioinformatics tools. These developments have led to the formation of the European Virus Bioinformatics Center, a network of experts in virology and bioinformatics who are joining forces to enable extensive exchange and collaboration between these research areas. The EVBC strives to provide talented researchers with a supportive environment free of gender bias, but the gender gap in science, especially in math-intensive fields such as computer science, persists. To bring more talented women into research and keep them there, we need to highlight role models to spark their interest, and we need to ensure that female scientists are not kept at lower levels but are given the opportunity to lead the field. Here we showcase the work of the EVBC and highlight the achievements of some outstanding women experts in virology and viral bioinformatics.

Genome-wide diversity of Zika virus: Exploring spatio-temporal dynamics to guide a new nomenclature proposal.

The Zika virus (ZIKV) disease caused a public health emergency of international concern that started in February 2016. The overall number of ZIKV-related cases increased until November 2016, after which it declined sharply. While the evaluation of the potential risk and impact of future arbovirus epidemics remains challenging, intensified surveillance efforts along with a scale-up of ZIKV whole-genome sequencing provide an opportunity to understand the patterns of genetic diversity, evolution, and spread of ZIKV. However, a classification system that reflects the true extent of ZIKV genetic variation is lacking. Our objective was to characterize ZIKV genetic diversity and phylodynamics, identify genomic footprints of differentiation patterns, and propose a dynamic classification system that reflects its divergence levels. We analysed a curated dataset of 762 publicly available sequences spanning the full-length coding region of ZIKV from across its geographical span and collected between 1947 and 2021. The definition of genetic groups was based on comprehensive evolutionary dynamics analyses, which included recombination and phylogenetic analyses, within- and between-group pairwise genetic distances comparison, detection of selective pressure, and clustering analyses. Evidence for potential recombination events was detected in a few sequences. However, we argue that these events are likely due to sequencing errors as proposed in previous studies. There was evidence of strong purifying selection, widespread across the genome, as also detected for other arboviruses. A total of 50 sites showed evidence of positive selection, and for a few of these sites, there was amino acid (AA) differentiation between genetic clusters. Two main genetic clusters were defined, ZA and ZB, which correspond to the already characterized 'African' and 'Asian' genotypes, respectively. Within ZB, two subgroups, ZB.1 and ZB.2, represent the Asiatic and the American (and Oceania) lineages, respectively. ZB.1 is further subdivided into ZB.1.0 (a basal Malaysia sequence sampled in the 1960s and a recent Indian sequence), ZB.1.1 (South-Eastern Asia, Southern Asia, and Micronesia sequences), and ZB.1.2 (very similar sequences from the outbreak in Singapore). ZB.2 is subdivided into ZB.2.0 (basal American sequences and the sequences from French Polynesia, the putative origin of South America introduction), ZB.2.1 (Central America), and ZB.2.2 (Caribbean and North America). This classification system does not use geographical references and is flexible to accommodate potential future lineages. It will be a helpful tool for studies that involve analyses of ZIKV genomic variation and its association with pathogenicity and serve as a starting point for the public health surveillance and response to on-going and future epidemics and to outbreaks that lead to the emergence of new variants.

The Impact of Genital Ulcers on HIV Transmission Has Been Underestimated-A Critical Review.

In the early 1990s, several observational studies determined that genital ulcer disease (GUD), in either the index or the exposed person, facilitates HIV transmission. Several meta-analyses have since presented associated risk ratios (RR) over the baseline per-act transmission probability (PATP) usually in the range of 2-5. Here we review all relevant observational studies and meta-analyses, and show that the estimation of RRs was, in most cases, biased by assuming the presence of GUD at any time during long follow-up periods, while active genital ulcers were present in a small proportion of the time. Only two studies measured the GUD co-factor effect in PATPs focusing on acts in which ulcers were present, and both found much higher RRs (in the range 11-112). We demonstrate that these high RRs can be reconciled with the studies on which currently accepted low RRs were based, if the calculations are restricted to the actual GUD episodes. Our results indicate that the effect of genital ulcers on the PATP of HIV might be much greater than currently accepted. We conclude that the medical community should work on the assumption that HIV risk is very high during active genital ulcers.

Exploring the mechanisms behind HIV drug resistance in sub-Saharan Africa: conceptual mapping of a complex adaptive system based on multi-disciplinary expert insights.

BACKGROUND: HIV drug resistance (HIVDR) continues to threaten the effectiveness of worldwide antiretroviral therapy (ART). Emergence and transmission of HIVDR are driven by several interconnected factors. Though much has been done to uncover factors influencing HIVDR, overall interconnectedness between these factors remains unclear and African policy makers encounter difficulties setting priorities combating HIVDR. By viewing HIVDR as a complex adaptive system, through the eyes of multi-disciplinary HIVDR experts, we aimed to make a first attempt to linking different influencing factors and gaining a deeper understanding of the complexity of the system. METHODS: We designed a detailed systems map of factors influencing HIVDR based on semi-structured interviews with 15 international HIVDR experts from or with experience in sub-Saharan Africa, from different disciplinary backgrounds and affiliated with different types of institutions. The resulting detailed system map was conceptualized into three main HIVDR feedback loops and further strengthened with literature evidence. RESULTS: Factors influencing HIVDR in sub-Saharan Africa and their interactions were sorted in five categories: biology, individual, social context, healthcare system and 'overarching'. We identified three causal loops cross-cutting these layers, which relate to three interconnected subsystems of mechanisms influencing HIVDR. The 'adherence motivation' subsystem concerns the interplay of factors influencing people living with HIV to alternate between adherence and non-adherence. The 'healthcare burden' subsystem is a reinforcing loop leading to an increase in HIVDR at local population level. The 'ART overreliance' subsystem is a balancing feedback loop leading to complacency among program managers when there is overreliance on ART with a perceived low risk to drug resistance. The three subsystems are interconnected at different levels. CONCLUSIONS: Interconnectedness of the three subsystems underlines the need to act on the entire system of factors surrounding HIVDR in sub-Saharan Africa in order to target interventions and to prevent unwanted effects on other parts of the system. The three theories that emerged while studying HIVDR as a complex adaptive system form a starting point for further qualitative and quantitative investigation.