- 1Nuffield Department of Population Health, Infectious Disease Epidemiology Unit, National Perinatal Epidemiology Unit, University of Oxford, Oxford, United Kingdom
- 2Centre for Statistics in Medicine, Nuffield Department of Orthopaedics, Rheumatology and Musculoskeletal Sciences, University of Oxford, Botnar Research Centre, Oxford, United Kingdom
Introduction: Global HIV infections due to HIV-1 recombinants are increasing and impede prevention and treatment efforts. Key populations suffer most new HIV infections, but their role in the spread of HIV-1 recombinants is unknown. We conducted a global analysis of the associations between key populations and HIV-1 recombinants.
Methods: We searched PubMed, EMBASE, CINAHL, and Global Health for HIV-1 subtyping studies published from 1/1/1990 to 31/12/2015. Unpublished data was collected through a global survey. We included studies with HIV-1 subtyping data of key populations collected during 1990-2015. Key populations assessed were heterosexual people (HET), men who have sex with men (MSM), people who inject drugs (PWID), vertical transmissions (VERT), commercial sex workers (CSW), and transfusion-associated infections (BLOOD). Logistic regression was used to determine associations of key populations with HIV-1 recombinants. Subgroup analyses were performed for circulating recombinant forms (CRFs), unique recombinant forms (URFs), regions, and time periods.
Results: Eight hundred and eighty five datasets including 77,284 participants from 83 countries were included. Globally, PWID were associated with the greatest odds of recombinants and CRFs (OR 2.6 [95% CI 2.46–2.74] and 2.99 [2.83–3.16]), compared to HET. CSW were associated with increased odds of recombinants and URFs (1.59 [1.44–1.75] and 3.61 [3.15–4.13]). VERT and BLOOD were associated with decreased odds of recombinants (0.58 [0.54–0.63] and 0.43 [0.33–0.56]). MSM were associated with increased odds of recombinants in 2010–2015 (1.43 [1.35–1.51]). Subgroup analyses supported our main findings.
Discussion: As PWID, CSW, and MSM are associated with HIV-1 recombinants, increased preventative measures and HIV-1 molecular surveillance are crucial within these key populations.
Systematic review registration: PROSPERO [CRD42017067164].
1. Introduction
In 2021, 38.4 million people were living with HIV worldwide and 1.5 million people became newly infected (1). The HIV pandemic remains a major global health challenge, and its extreme global genetic diversity impedes treatment and prevention efforts (2). Global temporal analysis indicates that the HIV-1 pandemic is diversifying, with increases in both the numbers of distinct HIV-1 variants and proportions of recombinant strains (3–5). Increasing diversity impacts HIV diagnosis and treatment, drug resistance, viral load measurement, transmission, disease progression, immune responses, and vaccine development (2, 6–10).
After zoonotic transmission from chimpanzees to humans in Central Africa around 1900, the HIV-1 group M epidemic rapidly diversified into distinct subtypes, designated by the letters A–D, F–H, and J–L (11, 12). HIV-1 subtypes spread across the globe throughout the 20th century, resulting in HIV-1 subtype distributions that greatly vary by region (3, 13). The genetic complexity of the HIV pandemic continues to increase over time, largely driven by the high mutation and recombination rates of the error-prone reverse transcriptase enzyme (14). Recombination occurs when an individual is co-infected with multiple strains which combine into a new variant (15). The resulting variants are designated as circulating recombinant forms (CRFs) or unique recombinant forms (URFs). CRFs, which are characterized by community spread, must be fully sequenced and found in at least three epidemiologically unlinked individuals. More than 120 distinct CRFs have been described to date, and more CRFs continue to be identified (16). URFs are unique recombinant sequences without evidence of onward transmission. The proportion of recombinants has been increasing over time, both globally and in most regions, and recombinants now constitute nearly a quarter of all HIV-1 infections globally (4). In addition to increasing the genetic complexity of the HIV pandemic, recombination may confer an evolutionary advantage, leading to altered transmission and/or virulence (17, 18).
In 2021, 70% of new HIV infections occurred within key populations and their sexual partners, though these populations account for <5% of the global population (1). It is estimated that men who have sex with men (MSM) have 28 times the risk of HIV infection relative to heterosexual (HET) adult men, female commercial sex workers (CSW) have 30 times the risk relative to other adult women, and people who inject drugs (PWID) have 35 times the risk compared to those who do not inject drugs (1). Additionally, people in areas without comprehensive blood screening are particularly vulnerable to HIV infection through transfusions with infected blood (BLOOD) (19), and children born to mothers with HIV can become infected via vertical transmission (VERT) during pregnancy, labor, delivery, or breastfeeding (20). Prior work indicates that HIV can follow a chain of transmission among these groups, spreading from PWID to CSW who transmit the virus to their HET clients. The virus can then be transmitted to the client's female sexual partner before VERT transmission of HIV infection to children (20, 21). Transmission among MSM and during blood transfusions has also played a major historical role in the spread of HIV, particularly across Asia, Europe, and North America (19, 21, 22).
Though these key populations are known to play a role in HIV transmission, it is unclear what role they play in the spread of HIV-1 recombinants. Since these populations often face difficulties accessing HIV services and have an increased risk of infection (1), potentially by multiple strains, they may be more likely to develop novel HIV strains. These recombinant strains may cross from key populations into the general population, making the overall HIV epidemic more complex.
The global proportion of HIV infections with recombinants is increasing and key populations globally account for most new HIV infections. However, there is an evidence gap regarding the global association of key populations with HIV-1 recombinants. To address this gap, we conducted a global analysis of the association between multiple key populations and HIV-1 recombinants using the largest global HIV-1 molecular epidemiology database assembled to date.
2. Materials and methods
2.1. Data collection
Data on the global distribution of HIV-1 subtypes and recombinants among key populations were obtained through a systematic literature review (PROSPERO: CRD42017067164), review of specialist journals and reports, and global survey of experts (3). We searched PubMed (29,825 citations retrieved), Embase (Ovid) (25,914 citations), CINAHL (Ebscohost) (451 citations), and Global Health (Ovid) (9,707 citations) for studies reporting HIV-1 subtyping data published from Jan 1, 1990 to Dec 31, 2015. This time period covers the period for which reliable estimates of national HIV prevalence were available. Search terms were Medical Subject Headings (MeSH) and Emtree terms, free text words, and synonyms, including “HIV,” “Subtype,” “recombinant,” “CRF,” and “URF” (Appendix pp2–5). No language or methodology filters were used. All references retrieved were combined in Endnote reference manager, and duplicates removed (Endnote X9; Clarivate Analytics, Philadelphia, PA). Authors RE, JY, LD-T and JH screened titles and abstracts, retrieved relevant full text articles, and assessed articles against the eligibility criteria. Additional published data were derived from the WHO HIV Drug Resistance Report 2012 (23), published reviews and reports on HIV diversity, and papers indexed on Scopus that referenced previous publications on global HIV-1 diversity (Appendix pp6–8). Additionally, four specialist journals (AIDS, Journal of AIDS, Journal of Virology, AIDS Research and Human Retroviruses) were screened for relevant articles published between January 1990 and February 2016. Using a data collection template, unpublished original HIV-1 subtyping data was collected through a global survey of members in the WHO-UNAIDS Network for HIV Isolation and Characterisation.
2.2. Eligibility criteria and data extraction
Studies were eligible for inclusion if they were prevalence studies of key populations living with HIV with original HIV-1 subtyping data, known country and year of sample collection (1990–2015), and a minimum of 20 participants. Studies that only contained incident infections or untyped samples were excluded. Full-length genomes or any genome segment could be used for subtyping, no minimum sequence length was specified, all online subtyping tools were accepted, and subtyping data from each included dataset was assumed to be correct.
Authors RE, JY, LD-T, and JH extracted the following information for each data set: country, city or region, sample collection year(s), study type, key population, HIV-1 subtyping method(s), and genome segment(s) analyzed. The primary outcome was the number of each key population designated by the original authors as each HIV-1 subtype (A, B, C, D, F, G, H, J, K), CRFs, and URFs. Country designation was based on where samples were taken. One subtype/CRF/URF was assigned to each participant. Subtyping methods included sequencing, heteroduplex mobility assay, and serotyping. The vast majority of data was acquired by sequencing (100% in 2010–2015), mostly of partial genome sequences, mainly pol (94.4% in 2010–2015) (3). Contributing researchers were assumed to have obtained consent from participants, and no personal identifiable information was retrieved. Formal assessment of individual study quality was not performed. Discrepancies were resolved by the senior reviewer (JH).
2.3. Key populations
Based on the populations specified by each study, participants were categorized as heterosexual (HET), men who have sex with men (MSM), people who inject drugs (PWID), vertical transmissions (VERT), commercial sex workers (CSW), and transfusion-associated infections (BLOOD) by author NN and confirmed by JH. Studies involving multiple key populations were assigned to the key population comprising at least 95% of data or excluded if no single key population met the 95% threshold. Studies with unspecified or indeterminate key populations were excluded. Any discrepancies or ambiguities were resolved by JH.
2.4. Meta-analysis
As most studies provided data on a single key population, one-stage meta-analysis of individual-participant data of different studies was performed. For logistic regression, HIV-1 variants were categorized as “Subtype” or “Recombinant” (CRF/URF). A univariate binomial logistic regression model was constructed to analyse the global association of each key population with HIV-1 recombinants. To assess the global association of key populations with CRFs and URFs separately, logistic regression was repeated using a multinomial model.
Countries were grouped into 14 regions (Appendix p9) and data were assigned to four periods: 1990–1999, 2000–2004, 2005–2009, and 2010–2015 (Appendix p10). All participants in each dataset were assigned to periods based on the midpoint year of the reported sample collection period. Datasets of which sampling years were evenly split between two periods (e.g., 2003–2006) were excluded from time-stratified analyses. To assess temporal and geographic differences in the associations with recombinants, the binomial logistic model was separately stratified into subgroups by time period and region.
For all logistic models, “Subtype” was used as the reference group. Odds Ratios (ORs) were reported with 95% confidence intervals. In the Appendix, pairwise ORs are reported globally and for each region (pp11–13) and period subgroup analyses were repeated for the multinomial logistic regression model (pp14). Statistical analyses were performed using STATA 17.0 (StataCorp LLC, College Station, TX). This systematic review is reported according to the PRISMA guidelines, as applicable.
3. Results
3.1. Data collection
A total of 885 datasets including 77,284 participants from key populations from 83 countries were included (Figure 1). The systematic literature search yielded 208 datasets comprising 23,988 participants. Six hundred and seventeen datasets with 48,984 participants were collected from the global survey, and 60 datasets with 4,312 participants were obtained from other published sources.
Figure 1. Data collection flowchart. *For example, HIV-positive immigrants only. †For example, data only provided for subtype B and non-B participants. ‡For example, subtypes referred to disease states, not HIV subtypes. §For example, HIV subtyping data could not be assigned to a specific key population, as multiple key populations were present in the study.
Most included participants were heterosexual people (58.2%) and MSM (25.7%), with smaller proportions of participants representing PWID (8.1%), VERT (5.3%), CSW (2.2%), and BLOOD (0.6%) (Table 1). Data from the HET population was the largest in each time period, while there was no data for CSW in the most recent period. Most participants were derived from Western and central Europe, and North America (WCENA) followed by East, West, and Southern Africa. HET data was available in every geographic region while all BLOOD participants were derived from East Asia and WCENA. HET was selected as the reference group as it was the only population represented in all regions and time periods.
3.2. Global association of key populations with recombinants
The global distribution of HIV-1 subtypes, CRFs and URFs among key populations in 1990–2015 is shown in Table 2. During 1990–2015, the largest proportion of recombinant infections was found among PWID (52.8%) and CSW (40.5%), followed by HET (30.1%), MSM (29.4%), VERT (19.9%), and BLOOD (15.6%). PWID had the highest proportion of CRFs (49.1%) and CSW had the highest proportion of URFs (17.6%). The proportion of recombinant infections grew consistently across periods for HET, MSM, and BLOOD and across the first three periods for PWID. Global proportions of individual CRFs for each key population are included in the Appendix p16.
PWID were associated with the greatest odds of recombinants relative to all other key populations (Figure 2A; Table 3). Relative to HET, PWID and CSW were associated with increased odds of recombinants [OR 2.6 (95% CI 2.46–2.74) and 1.59 (1.44–1.75)], while VERT and BLOOD were associated with decreased odds [0.58 (0.54–0.63) and 0.43 (0.33–0.56)]. MSM did not have a significant association with recombinants relative to HET [0.97 (0.93–1.01)] (Figure 2A; Table 3).
Figure 2. Global associations of key populations with HIV-1 recombinants, CRFs, and URFs relative to heterosexual people (1990–2015). (A) Odds ratios for HIV-1 recombinants, CRFs, and URFs, compared to HIV-1 subtypes, of key populations relative to heterosexual people (1990–2015). (B) Odds ratios for HIV-1 recombinants, compared to HIV-1 subtypes, of key populations relative to heterosexual people in each time period (1990–1999, 2000–2004, 2005–2009, 2010–2015). No data on recombinants was available for BLOOD in 1990–1999, and no data was available for CSW in 2010–2015. Error bars represent the 95% confidence intervals. Square areas are proportional to the number of participants in each key population analyzed. Odds ratios and 95% CI are provided in the Appendix p15 (*P < 0.05, **P < 0.01). BLOOD, blood/plasma transfusion associated infections; CRF, circulating recombinant form; CSW, commercial sex workers; HET, Heterosexual; MSM, men who have sex with men; PWID, people who inject drugs; URF, Unique Recombinant Form; VERT, vertical transmission (mother to child).
Table 3. Regional associations of HIV-1 recombinants, CRFs, and URFs with key populations relative to heterosexual people, 1990–2015.
Relative to HET, independent associations of each key population with CRFs and URFs varied substantially (Figure 2A; Appendix p15). PWID were significantly associated with increased odds of CRFs [2.99 (2.83–3.16)], while CSW were significantly associated with increased odds of URFs [3.61 (3.15–4.13)]. VERT was associated with decreased odds of CRFs [0.49 (0.44–0.53)], and BLOOD was associated with decreased odds of both CRFs and URFs [0.43 (0.32–0.56) and 0.45 (0.26–0.77)]. MSM were associated with slightly increased odds of CRFs [1.09 (1.05–1.14)], but decreased odds of URFs [0.44 (0.40–0.48)] (Figure 2A; Appendix p15).
Relative to HET, the strength of associations with recombinants across time differed by key population (Figure 2B; Appendix pp14, 15). PWID were associated with increased odds of recombinants across all periods. CSW were initially associated with increased odds of recombinants, but the strength of the association decreased with time before leading to decreased odds in the 2005–2009 period. VERT alternated from increased to decreased odds of recombinants across time. BLOOD was associated with decreased odds of recombinants during 2000–2009 but was not significant in the most recent period. MSM had decreased odds of recombinants during 1990–2009 but trended upwards and were associated with increased odds in 2010–2015 [1.43 (1.35–1.51)].
3.3. Regional association of key populations with recombinants
The regional distribution of HIV-1 subtypes, CRFs, and URFs for each key population is included in the Appendix pp17–19. The association of key populations with recombinants varied by region (Table 3). Compared to HET, PWID had the greatest odds of recombinants in Eastern Europe and central Asia [EECA; OR 19.98 (95% CI 6.30–63.34)], followed by Latin America [4.16 (3.02–5.74)] and East Asia [3.20 (2.50–4.09)]. PWID were significantly associated with CRFs in EECA [54.89 (7.62–395.26)] and East Asia [3.26 (2.55–4.18)], and both CRFs and URFs in Latin America [2.72 (1.45–5.11) and 4.80 (3.36–6.85)]. In contrast, PWID had a strong negative association with recombinants in Southeast Asia [0.06 (0.05–0.07)].
CSW had the largest odds of recombinants in Latin America [15.58 (10.63–22.85)] and East Africa [3.36 (2.85–3.97)]. In both regions, CSW were significantly associated with URFs [22.46 (15.17–33.23) and 3.81 (3.22–4.51)]. However, in West Africa, CSW had decreased odds of recombinants [0.68 (0.56–0.82)], CRFs [0.77 (0.64–0.93)], and URFs [0.06 (0.01–0.23)].
VERT was associated with decreased odds of recombinants in Southern Africa [0.14 (0.03–0.56)] and WCENA [0.69 (0.55–0.88)]. VERT was not independently associated with CRFs or URFs in Southern Africa, but had decreased odds of both CRFs [0.74 (0.58–0.94)] and URFs [0.51 (0.28–0.91)] in WCENA. In Latin America [3.10 (2.50–3.85)], West Africa [3.39 (2.04–5.63)], and SE Asia [12.12 (1.69–86.98)] VERT was associated with increased odds of recombinants. In Latin America, VERT was significantly associated with increased odds of both CRFs [1.85 (1.20–2.86)] and URFs [3.65 (2.86–4.66)], but it was only associated with increased odds of CRFs in West Africa [3.40 (2.03–5.68)] and SE Asia [12.28 (1.71–88.12)].
BLOOD was associated with decreased odds of recombinants in East Asia [0.04 (0.02–0.06)] and WCENA [0.70 (0.51–0.96)]. In both regions, BLOOD was also associated with decreased odds of CRFs [0.04 (0.02–0.06) and 0.65 (0.45–0.92)].
In East Asia, MSM were associated with increased odds of recombinants [1.94 (1.60–2.34)] and CRFs [1.92 (1.59–2,32)]. In West Africa, MSM were associated with increased odds of recombinants [1.44 (1.14–1.82)] and URFs [2.44 (1.80–3.30)]. In Latin America and WCENA, MSM had decreased odds of recombinants [0.46 (0.31–0.67) and 0.18 (0.16–0.20)] and URFs [0.66 (0.45–0.97) and 0.18 (0.15–0.23)].
4. Discussion
A strong association between PWID and recombinants and CRFs was observed globally across all periods and in most regions. Only in SE Asia, where CRF01_AE has a prevalence of ~70–80% (3), were PWID associated with decreased odds of recombinant strains. The prevalence of recombinant epidemics among PWID in most regions, where “pure” HIV-1 subtypes are typically the most prevalent overall (3), and subtype-based epidemics among PWID in SE Asia, where recombinant strains are highly prevalent, indicates that HIV-1 circulates among PWID via transmission networks distinct from the HET population. This finding extends previous studies suggesting that HIV is transmitted among independent PWID networks across multiple continents (24, 25). Furthermore, the association with recombinants across all periods indicates that PWID play a major role in the global diversification of HIV-1.
CSW were associated with increased odds of recombinants and URFs, particularly in the periods 1990–1999 and 2000–2004. Across East Africa and Latin America, CSW were significantly more likely to be infected with URFs than the HET population. These findings highlight that novel HIV strains frequently arise within the CSW population. URFs arise independently and lack evidence of transmission, minimizing the likelihood that the observed association is due to reverse causation. The decreased odds of CRFs in West Africa may be related to the high prevalence of CRF02_AG (3), similarly to the case of PWID and CRF01_AE in SE Asia. Additional data is required to identify factors contributing to the diminishing global association of CSWs with recombinants across time.
Though VERT was associated with decreased odds of recombinants and CRFs, results greatly varied across times and regions. While biological differences in recombinant strains may cause increased rates of vertical transmission relative to “pure” subtypes (26), high levels of heterogeneity indicate that VERT is not a major driver of increasing HIV-1 diversity.
BLOOD was associated with decreased odds of recombinants and CRFs in both East Asia and WCENA. Particularly in East Asia, where CRF01_AE is highly prevalent, blood transfusion recipients were significantly less likely to have a recombinant strain of HIV than the heterosexual population, which may reflect the geographical origins of the blood donor base. However, the small number of datasets means that the observed association is subject to limitations of power and representativeness. Additional data is required to clarify the association between BLOOD and recombinants.
MSM did not have a significant global association with recombinants overall, likely due to a positive association with CRFs and negative association with URFs. The positive association between MSM and CRFs was strongest in East Asia where the prevalence of CRFs has grown from 25.9 to 75.5% during 1990–2015 (3). Within this region, MSM had nearly double the odds of CRFs as HET, indicating that MSM may be at the forefront of the growing epidemic across East Asia. Similar results were seen in West Africa where the proportion of URFs grew from 3.4 to 15.5% over the same period (3), and findings indicated that MSM had 2.44 times the odds of being infected with URFs. These associations suggest that MSM likely play a major role in the spread of new strains in some regions. Despite an overall association with recombinants that was not significant and historical associations with HIV-1 subtype B (27), the significant association in 2010–2015 and increasing trend across time indicate that MSM may be associated with an increased risk of recombinants.
A key strength of this study is its unprecedented large size, including 77,284 participants from 83 countries, collected from key populations globally during 1990–2015. To our knowledge, this is the first comprehensive analysis of the association between key populations and HIV-1 subtypes and recombinants at a global and regional level. Additionally, data was collected through both a literature search and a global survey, with the inclusion of unpublished data enabling increased regional coverage and improved coverage of recent time periods.
The study also had some limitations. Estimates of associations of key populations with HIV-1 variants are dependent on the underlying data. There was notable variation in coverage by key population, geographic region, and time period. Although the numbers of participants were generally high, the limited number of datasets included from BLOOD means that results must be interpreted with caution. Conclusions could not be independently drawn for any key populations in the Caribbean, Ethiopia, Oceania, and Middle East and North Africa (MENA) due to persistent data gaps that have been previously noted (1, 28). Similarly, the absence of data for certain key populations in some regions (e.g., MSM data in Central Africa, East Africa, and MENA) may reflect limited access to healthcare due to sociolegal restrictions (29). Most data were not drawn from nationally representative surveys and we were unable to weigh country-level data according to relative numbers of people of key populations living with HIV in each country, as comprehensive global data on key populations is not available. Hence, reported distributions of HIV-1 variants should not be interpreted as representative of key populations in each region or globally. As HIV subtyping data for most studies was primarily based on pol sequencing rather than the whole genome (3), recombination outside of this genome region was likely missed, leading to an underestimation of recombinants. Seventy four CRFs were described at the time of data collection (up until 2015), contributing to the discrepancy between the 48 CRFs identified within the datasets contributing to this study and the >120 CRFs that have been described to date (16). Findings could be subject to bias due to heterogeneity in study design, inclusion/exclusion criteria, subtyping methods, and rates of treatment and migration across regions. In particular, differences in participant recruitment and the definition of key populations between studies could affect observed associations with recombinants. Lastly, insufficient data was available to conduct analysis for transgender women.
Among key populations, increased risk of HIV infection, potentially by multiple strains, and difficulty accessing treatment, potentially leading to increased viral loads, may contribute to the formation and onward spread of HIV-1 recombinants (18, 30). The increasing diversity of the HIV pandemic has implications across diagnosis, treatment, and prevention (2, 6–10). Efforts to prevent the spread of novel HIV strains should consider approaches for key populations such as PWID, CSW, and MSM that are at increased risk of developing and transmitting recombinants. In the case of PWID, this may require prevention-based approaches such as distribution of sterile injection equipment (31, 32), opiate substitution treatment (33), and increased access to antiretroviral therapies (34). For CSW and MSM, prevention efforts should focus on increasing availability of the dapivirine vaginal ring for cisgender women, oral TDF-based pre-exposure prophylaxis (PrEP), and long-acting injectable cabotegravir (1, 35, 36). Increased HIV testing among key populations will help detect and treat new HIV infections early. These efforts can help limit the spread of traits from newly-emergent, highly virulent strains (18). Structural reform may also be necessary as the criminalization of these three key populations is associated with worse HIV outcomes and inadequate viral suppression (37, 38), potentially accelerating HIV-1 diversification.
In summary, this is the first study to comprehensively analyse the global association of key populations with HIV-1 recombinants. PWID, CSW, and MSM were significantly associated with recombinants globally and across multiple regions. As key populations and their partners account for 70% of new HIV infections (1), it is apparent that key populations are driving the genetic diversification of the global HIV-1 pandemic, posing a challenge to diagnostics, treatments, and vaccines against HIV. Therefore, additional surveillance of HIV-1 molecular epidemiology and increased preventative measures should be targeted toward these key populations.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding author.
Ethics statement
Ethical review and approval was not required for the study on human participants in accordance with the local legislation and institutional requirements. Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.
Author contributions
NN assessed eligibility of manuscripts, conducted the analyses, designed figures and tables, interpreted the data, and wrote the manuscript. RE, JY, and LD-T screened the electronic literature search results for relevant manuscripts, assessed their eligibility, extracted data, and collected additional published data. SK designed and did the electronic literature search. JH conceived, designed, coordinated the study, wrote the systematic review protocol, assisted with the literature search, assessed eligibility of manuscripts, collected additional published data, conducted the global survey, extracted data, designed the analysis plan, interpreted the data, and wrote the manuscript. All authors read and approved the final version of the manuscript.
WHO-UNAIDS Network for HIV Isolation and Characterisation
Alash'le G. Abimiku, Simon Agwale, Chris Archibald, Boaz Avidor, María Gabriela Barbás, Francoise Barre-Sinoussi, Banson Barugahare, El Hadj Belabbes, Silvia Bertagnolio, Deborah Birx, Aleksei F. Bobkov, James Brandful, Helba Bredell, Catherine A. Brennan, James Brooks, Marie Bruckova, Luigi Buonaguro, Franco Buonaguro, Stefano Buttò, Anne Buve, Mary Campbell, Jean Carr, Alex Carrera, Manuel Gómez Carrillo, Connie Celum, Beth Chaplin, Macarthur Charles, Dimitrios Chatzidimitriou, Zhiwei Chen, Katsumi Chijiwa, David Cooper, Philip Cunningham, Anoumou Dagnra, Cillian F. de Gascun, Julia Del Amo, Elena Delgado, Ursula Dietrich, Dominic Dwyer, Dennis Ellenberger, Barbara Ensoli, Max Essex, Herve Fleury, Peter N. Fonjungo, Vincent Foulongne, Deepak A. Gadkari, Feng Gao, Federico García, Roger Garsia, Guy Michel Gershy-Damet, Judith R. Glynn, Ruth Goodall, Zehava Grossman, Monick Lindenmeyer Guimarães, Beatrice Hahn, Raph L. Hamers, Osamah Hamouda, Ray Handema, Xiang He, Joshua Herbeck, David D. Ho, Africa Holguin, Mina Hosseinipour, Gillian Hunt, Masahiko Ito, Mohamed Ali Bel Hadj Kacem, Erin Kahle, Pontiano Kaleebu, Marcia Kalish, Adeeba Kamarulzaman, Chun Kang, Phyllis Kanki, Edward Karamov, Jean-Claude Karasi, Kayitesi Kayitenkore, Tony Kelleher, Dwip Kitayaporn, Leondios G. Kostrikis, Claudia Kucherer, Claudia Lara, Thomas Leitner, Kirsi Liitsola, Jai Lingappa, Marek Linka, Ivette Lorenzana de Rivera, Vladimir Lukashov, Shlomo Maayan, Luzia Mayr, Francine McCutchan, Nicolas Meda, Elisabeth Menu, Fred Mhalu, Doreen Mloka, John L. Mokili, Brigitte Montes, Orna Mor, Mariza Morgado, Fausta Mosha, Awatef Moussi, James Mullins, Rafael Najera, Mejda Nasr, Nicaise Ndembi, Joel R. Neilson, Vivek R. Nerurkar, Florian Neuhann, Claudine Nolte, Vlad Novitsky, Philippe Nyambi, Marianna Ofner, Fem J. Paladin, Anna Papa, Jean Pape, Neil Parkin, Chris Parry, Martine Peeters, Alexandra Pelletier, Lucía Pérez-Álvarez, Deenan Pillay, Angie Pinto, Trinh Duy Quang, Cecilia Rademeyer, Filimone Raikanikoda, Mark A. Rayfield, Jean-Marc Reynes, Tobias Rinke de Wit, Kenneth E. Robbins, Morgane Rolland, Christine Rousseau, Jesus Salazar-Gonzales, Hanan Salem, Mika Salminen, Horacio Salomon, Paul Sandstrom, Mario L. Santiago, Abdoulaye D. Sarr, Bryan Schroeder, Michel Segondy, Philippe Selhorst, Sylvester Sempala, Jean Servais, Ansari Shaik, Yiming Shao, Amine Slim, Marcelo A. Soares, Elijah Songok, Debbie Stewart, Julie Stokes, Shambavi Subbarao, Ruengpung Sutthent, Jun Takehisa, Amilcar Tanuri, Kok Keng Tee, Kiran Thapa, Michael Thomson, Tyna Tran, Willy Urassa, Hiroshi Ushijima, Philippe van de Perre, Guido van der Groen, Kristel van Laethem, Joep van Oosterhout, Ard van Sighem, Eric van Wijngaerden, Anne-Mieke Vandamme, Jurgen Vercauteren, Nicole Vidal, Lesley Wallace, Carolyn Williamson, Dawit Wolday, Jianqing Xu, Chunfu Yang, Linqi Zhang, Rong Zhang.
Affiliations
Institute of Human Virology, University of Maryland, Baltimore, MD, USA (A. G. Abimiku, J. Carr); Gede Foundation, Abuja, Nigeria (S. Agwale); Public Health Agency of Canada, Ottawa, ON, Canada (C. Archibald, J. Brooks, M. Ofner, P. Sandstrom, J. Stokes); Tel-Aviv Sourasky Medical Center, Tel-Aviv, Israel (B. Avidor); Ministerio de Salud, Córdoba, Argentina (M. G. Barbás); Institut Pasteur, Paris, France (F. Barre-Sinoussi, E. Menu); Ministry of Health, Entebbe, Uganda (B. Barugahare); National Reference Laboratory on HIV/AIDS, Institut Pasteur d'Algérie, Algiers, Algeria (E. Belabbes); World Health Organization, Geneva, Switzerland (S. Bertagnolio); Office of the Global AIDS Coordinator, Washington, DC, USA (D. Birx); The D I Ivanovsky Institute of Virology, Moscow, Russia (A. F. Bobkov); Noguchi Memorial Institute for Medical Research, University of Ghana, Accra, Ghana (J. Brandful); University of Cape Town, Cape Town, South Africa (H. Bredell, C. Rademeyer, P. Selhorst, D. Stewart, C. Williamson); Abbott Laboratories, Chicago, IL, USA (C. A. Brennan); National Institute of Public Health, Prague, Czech Republic (M. Bruckova, M. Linka); AIDS Reference Center, National Cancer Institute “Fond. G. Pascale”, Naples, Italy (L. Buonaguro, F. Buonaguro); National AIDS Center, Istituto Superiore di Sanità, Rome, Italy (S. Buttò, B. Ensoli); Institute of Tropical Medicine, Antwerp, Belgium (A. Buve, G. van der Groen); University of Washington School of Medicine, Seattle, WA, USA (M. Campbell, C. Celum, J. Herbeck, E. Kahle, J. Lingappa, J. Mullins, M. Rolland, C. Rousseau); St Vincent's Hospital, Sydney, Australia (A. Carrera, P. Cunningham); University of Buenos Aires, Buenos Aires, Argentina (M. Carrillo, H. Salomon); Harvard T H Chan School of Public Health, Boston, MA, USA (B. Chaplin, P. Kanki); Gheskio Center, Port-au-Prince, Haiti (M. Charles, C. Nolte, J. Pape); Aristotle University of Thessaloniki, Thessaloniki, Greece (D. Chatzidimitriou, A. Papa); Chinese Academy of Medical Sciences, Peking Union Medical School, Beijing, China (Z. Chen, L. Zhang, R. Zhang); Fukuoka Institute of Health and Environmental Sciences, Kyushu University Hospital, Dazaifu, Japan (K. Chijiwa); The Kirby Institute, Sydney, NSW, Australia (D. Cooper, T. Kelleher, A. Pinto, A. Shaik); Faculté des Sciences de la Santé, Université de Lomé, Lomé, Togo (A. Dagnra); University College Dublin, Dublin, Ireland (C. de Gascun); Instituto de Salud Carlos III, Madrid, Spain (J. Del Amo, E. Delgado, R. Najera, L. Pérez-Álvarez, M. Thomson); Chemotherapeutisches Forschungsinstitut, Georg-Speyer-Haus, Frankfurt, Germany (U. Dietrich); Pathology West, Westmead Hospital, Westmead, NSW, Australia (D. Dwyer, K. Thapa, T. Tran); Centers for Disease Control and Prevention, Atlanta, GA, USA (D. Ellenberger, P. N. Fonjungo, M. A. Rayfield, K. E. Robbins, S. Subbarao, C. Yang); Harvard School of Public Health, Boston, MA, USA (M. Essex, V. Novitsky, A. D. Sarr); Duke University Medical Center, Durham, NC, USA (F. Gao); University of Bordeaux, Bordeaux, France (H. Fleury); Montpellier University Hospital, Montpellier, France (V. Foulongne, P. van de Perre); National AIDS Research Institute, Pune, India (D. A. Gadkari); Complejo Hospitalario Universitario de Granada, Granada, Spain (F. García); Royal Prince Alfred Hospital, Sydney, Australia (R. Garsia, H. Salem); HIV Laboratory Programme on AIDS/AFRO, World Health Organisation, Ouagadougou, Burkina Faso (G. M. Gershy-Damet); London School of Hygiene and Tropical Medicine, London, UK (J. R. Glynn); University College London, London, UK (R. Goodall); National HIV Reference Laboratory, Ministry of Health, Tel Aviv, Israel (Z. Grossman, O. Mor); Instituto Oswaldo Cruz, FIOCRUZ, Rio de Janeiro, Brazil (M. L. Guimarães, M. Morgado); University of Pennsylvania, Philadelphia, PN, USA (B. Hahn); Amsterdam Institute for Global Health and Development, Amsterdam, Netherlands (R. L. Hamers, T. Rinke de Wit); Robert Koch Institute, Berlin, Germany (O. Hamouda, C. Kucherer); Yamanashi Medical University, Yamanashi, Japan (R. Handema, M. Ito); National Center for AIDS/STD Control and Prevention, China CDC, Beijing, China (X. He, Y. Shao, J. Xu); Aaron Diamond AIDS Research Center, The Rockefeller University, New York, NY, USA (D. D. Ho, L. G. Kostrikis); Ramón y Cajal Research Institute, Hospital Universitario Ramón y Cajal de Madrid, Madrid, Spain (A. Holguin); University of North Carolina, Chapel Hill, NC, USA (M. Hosseinipour); National Institute for Communicable Diseases, Johannesburg, South Africa (G. Hunt); Charles Nicolle Hospital, Tunis, Tunisia (M. Kacem, A. Moussi, M. Nasr, A. Slim); Medical Research Council, Entebbe, Uganda (P. Kaleebu, C. Parry); Vanderbilt Institute for Global Health, Vanderbilt University School of Medicine, Nashville, TN, USA (M. Kalish); University of Malaya, Kuala Lumpur, Malaysia (A. Kamarulzaman, K-K Tee); Institute for Molecular Biology and Genetics and Medical College, Seoul National University, Seoul, Korea (C. Kang); Gamaleya Center for Epidemiology and Microbiology, Moscow, Russian (E. Karamov); National Reference Laboratory, Kigali, Rwanda (J-C Karasi); Emory University School of Medicine, Atlanta, GA, USA (K. Kayitenkore); HIV/AIDS Collaboration, Nonthaburi, Thailand (D. Kitayaporn); Karolinska Institute, Huddinge University Hospital, Stockholm, Sweden (C. Lara); Los Alamos National Laboratory, Los Alamos, NM, USA (T. Leitner); National Institute for Health and Wellfare, Helsinki, Finland (K. Liitsola, M. Salminen); National Autonomous University of Honduras, Tegucigalpa, Honduras (I. Lorenzana de Rivera); Academic Medical Center, University of Amsterdam, Amsterdam, Netherlands (V. Lukashov); Hadassah U Hospital, Jerusalem, Israel (S. Maayan); New York University School of Medicine, New York, NY, USA (L. Mayr, P. Nyambi); Henry M. Jackson Foundation for the Advancement of Military Medicine, Bethesda, MD, USA (F. McCutchan); Centre Muraz, Bobo-Dioulasso, Burkina Faso (N. Meda); Muhimbili University of Health Sciences, Dar-es-salaam, Tanzania (F. Mhalu, D. Mloka, F. Mosha, W. Urassa); University of Edinburgh, Edinburgh, UK (J. L. Mokili); Montpellier University Hospital, Montpellier, France (B. Montes, M. Segondy); Institute of Human Virology, Abuja, Nigeria (N. Ndembi); University of Washington, Seattle, WA, USA (J. R. Neilson); University of Hawaii, Honolulu, HI, USA (V. R. Nerurkar); University Clinic Heidelberg, Heidelberg, Germany & Lighthouse Trust, Lilongwe, Malawi (F. Neuhann); Research Institute for Tropical Medicine, Muntinlupa City, Manila, Philippines (F. J. Paladin, M. L. Santiago); Data First Consulting, Inc, Belmont, CA, USA (N. Parkin); University of Montpellier, Montpellier, France (M. Peeters, N. Vidal); Centre de Recherche Public-Santé, Luxembourg, Luxembourg (A. Pelletier, J. Servais); Africa Health Research Institute, Durban, KwaZulu-Natal, South Africa & Division of Infection and Immunity, University College London, London, UK (D. Pillay); Institute of International Health, University of Tokyo, Tokyo, Japan (T. D. Quang); University of Sydney, Sydney, NSW, Australia (F. Raikanikoda); Institut Pasteur du Cambodge, Phnom Penh, Cambodia (J-M Reynes); University of Alabama at Birmingham, Birmingham, AL, USA (J. Salazar-Gonzales); Auckland City Hospital, Auckland, New Zealand (B. Schroeder); Uganda Virus Research Institute, Entebbe, Uganda (S. Sempala); Instituto Nacional de Câncer, Rio de Janeiro, Brazil (M. A. Soares); Kenya Medical Research Institute, Nairobi, Kenya (E. Songok); National HIV Repository and Bioinformatic Center, Siriraj Hospital, Mahidol University, Thailand (R. Sutthent); Laboratory of Viral Pathogenesis, Kyoto University, Kyoto, Japan (J. Takehisa); Federal University of Rio de Janeiro, Rio de Janeiro, Brazil (A. Tanuri); Aino Health Science Center and Aino University, Tokyo, Japan (H. Ushijima); Rega Institute for Medical Research, KU Leuven, Belgium (K. van Laethem, E. van Wijngaerden, A.-M. Vandamme, J. Vercauteren); Department of Medicine, Blantyre, Malawi (J. van Oosterhout); Stichting HIV Monitoring, Amsterdam, Netherlands (A. van Sighem); Health Protection Scotland, Glasgow, UK (L. Wallace); and Ethiopian Health & Nutrition Research Institute, Addis Ababa, Ethiopia (D. Wolday).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpubh.2023.1153638/full#supplementary-material
References
2. Hemelaar J. The origin and diversity of the HIV-1 pandemic. Trends Mol Med. (2012) 18:182–92. doi: 10.1016/j.molmed.2011.12.001
3. Hemelaar J, Elangovan R, Yun J, Dickson-Tetteh L, Fleminger I, Kirtley S, et al. Global and regional molecular epidemiology of HIV-1, 1990–2015: a systematic review, global survey, and trend analysis. Lancet Infect Dis. (2019) 19:143–55. doi: 10.1016/S1473-3099(18)30647-9
4. Hemelaar J, Loganathan S, Elangovan R, Yun J, Dickson-Tetteh L, Kirtley S, et al. Country Level diversity of the HIV-1 pandemic between 1990 and 2015. J Virol. (2020) 95:1–15. doi: 10.1128/JVI.01580-20
5. Hemelaar J, Elangovan R, Yun J, Dickson-Tetteh L, Kirtley S, Gouws-Williams E, et al. Global and regional epidemiology of HIV-1 recombinants in 1990–2015: a systematic review and global survey. Lancet HIV. (2020) 7:e772–81. doi: 10.1016/S2352-3018(20)30252-6
6. Taylor BS, Hammer SM. The challenge of HIV-1 subtype diversity. N Engl J Med. (2008) 358:1590–602. doi: 10.1056/NEJMra0706737
7. Barouch DHD, Korber B. HIV-1 vaccine development after STEP. Annu Rev Med. (2010) 61:153–67. doi: 10.1146/annurev.med.042508.093728
8. Hemelaar J. Implications of HIV diversity for the HIV-1 pandemic. J Infect. (2013) 66:391–400. doi: 10.1016/j.jinf.2012.10.026
9. Fauci AS. An HIV vaccine is essential for ending the HIV/AIDS pandemic. J Am Med Assoc. (2017) 318:1535–6. doi: 10.1001/jama.2017.13505
10. Elangovan R, Jenks M, Yun J, Dickson-Tetteh L, Kirtley S, Hemelaar J, et al. Global and regional estimates for subtype-specific therapeutic and prophylactic HIV-1 vaccines: a modeling study. Front Microbiol. (2021) 12:1636. doi: 10.3389/fmicb.2021.690647
11. Robertson DL, Anderson JP, Bradac JA, Carr JK, Foley B, Funkhouser RK, et al. HIV-1 nomenclature proposal. Science. (2000) 288:55. doi: 10.1126/science.288.5463.55d
12. Worobey M, Gemmel M, Teuwen DE, Haselkorn T, Kunstman K, Bunce M, et al. Direct evidence of extensive diversity of HIV-1 in Kinshasa by 1960. Nature. (2008) 455:661–4. doi: 10.1038/nature07390
13. Tebit DDMD, Arts EJE. Tracking a century of global expansion and evolution of HIV to drive understanding and to combat disease. Lancet Infect Dis. (2011) 11:45–56. doi: 10.1016/S1473-3099(10)70186-9
14. Roberts JD, Bebenek K, Kunkel TA. The accuracy of reverse transcriptase from HIV-1. Science. (1988) 242:1171–3. doi: 10.1126/science.2460925
15. Fang G, Weiser B, Kuiken C, Philpott SM, Rowland-Jones S, Plummer F, et al. Recombination following superinfection by HIV-1. AIDS. (2004) 18:153–9. doi: 10.1097/00002030-200401230-00003
16. Los Alamos National Laboratory (2020). HIV Circulating Recombinant Forms (CRFs). Available online at: https://www.hiv.lanl.gov/components/sequence/HIV/crfdb/crfs.comp (accessed July 09, 2023).
17. Ariën KK, Vanham G, Arts EJ. Is HIV-1 evolving to a less virulent form in humans?. Nat Rev Microbiol. (2007) 5:141–51. doi: 10.1038/nrmicro1594
18. Wymant C, Bezemer D, Blanquart F, Ferretti L, Gall A, Hall M, et al. A highly virulent variant of HIV-1 circulating in the Netherlands. Science. (2022) 375:540–5. doi: 10.1126/science.abk1688
19. Scott SR, Wu Z. Risks and challenges of HIV infection transmitted via blood transfusion. Biosaf Health. (2019) 1:124–8. doi: 10.1016/j.bsheal.2019.12.001
20. Kourtis AP, Bulterys M, Nesheim SR, Lee FK. Understanding the timing of HIV transmission from mother to infant. JAMA. (2001) 285:709–12. doi: 10.1001/jama.285.6.709
21. Ruxrungtham K, Brown T, Phanuphak P. HIV/AIDS in Asia. Lancet. (2004) 364:69–82. doi: 10.1016/S0140-6736(04)16593-8
22. Qian ZH, Vermund SH, Wang N. Risk of HIV/AIDS in China: subpopulations of special importance. Sex Transm Infect. (2005) 81:442–7. doi: 10.1136/sti.2004.014258
23. Bertagnolio S World Health Organization World World Health Organization Department of HIV/AIDS. WHO HIV Drug Resistance Report 2012' (Geneva) (2012).
24. Tuan NA, Fylkesnes K, Thang BD, Hien NT, Long NT, Van Kinh N, et al. Human immunodeficiency virus (HIV) infection patterns and risk behaviours in different population groups and provinces in Viet Nam. Bull World Health Organ. (2007) 85:35–41. doi: 10.2471/BLT.05.026153
25. Jarlais DCD, Sypsa V, Feelemyer J, Abagiu AO, Arendt V, Broz D, et al. HIV outbreaks among people who inject drugs in Europe, North America, and Israel. Lancet HIV. (2020) 7:e434–42. doi: 10.1016/S2352-3018(20)30082-5
26. Koulinska IN, Villamor E, Msamanga G, Fawzi W, Blackard J, Renjifo B, et al. Risk of HIV-1 transmission by breastfeeding among mothers infected with recombinant and non-recombinant HIV-1 genotypes. Virus Res. (2006) 120:191–8. doi: 10.1016/j.virusres.2006.03.007
27. Lunar MM, Vandamme A-M, TomaŽič J, Karner P, Vovko TD, Pečavar B, et al. Bridging epidemiology with population genetics in a low incidence MSM-driven HIV-1 subtype B epidemic in Central Europe. BMC Infect Dis. (2015) 15:1–12. doi: 10.1186/s12879-015-0802-6
28. Mumtaz GR, Chemaitelly H, AlMukdad S, Osman A, Fahme S, Rizk NA, et al. Status of the HIV epidemic in key populations in the Middle East and north Africa: knowns and unknowns. Lancet HIV. (2022) 9:e506–16. doi: 10.1016/S2352-3018(22)00093-5
29. UNAIDS. Evidence for Eliminating HIV-related Stigma and Discrimination — Guidance for Countries to Implement Effective Programmes to Eliminate HIV-Related Stigma and Discrimination in Six Settings'. Geneva (2020).
30. Vuilleumier S, Bonhoeffer S. Contribution of recombination to the evolutionary history of HIV. Curr Opin HIV AIDS. (2015) 10:84–9. doi: 10.1097/COH.0000000000000137
31. Abdul-Quader AS, Feelemyer J, Modi S, Stein ES, Briceno A, Semaan S, et al. Effectiveness of structural-level needle/syringe programs to reduce HCV and HIV infection among people who inject drugs: a systematic review. AIDS Behav. (2013) 17:2878–92. doi: 10.1007/s10461-013-0593-y
32. Aspinall EJ, Nambiar D, Goldberg DJ, Hickman M, Weir A, Van Velzen E, et al. Are needle and syringe programmes associated with a reduction in HIV transmission among people who inject drugs: a systematic review and meta-analysis. Int J Epidemiol. (2014) 43:235–48. doi: 10.1093/ije/dyt243
33. MacArthur GJ, van Velzen E, Palmateer N, Kimber J, Pharris A, Hope V, et al. Interventions to prevent HIV and Hepatitis C in people who inject drugs: a review of reviews to assess evidence of effectiveness. Int J Drug Policy. (2014) 25:34–52. doi: 10.1016/j.drugpo.2013.07.001
34. Braitstein P, Brinkhof MWG, Dabis F, Schechter M, Boulle A, Miotti P, et al. Mortality of HIV-1-infected patients in the first year of antiretroviral therapy: comparison between low-income and high-income countries. Lancet. (2006) 367:817–24. doi: 10.1016/S0140-6736(06)68337-2
35. Bigna JJ, Temfack E. PrEP for men who have sex with men in Africa. Lancet HIV. (2021) 8:e388–9. doi: 10.1016/S2352-3018(21)00024-2
36. WHO. Consolidated Guidelines on HIV, Viral Hepatitis and STI Prevention, Diagnosis, Treatment and Care for Key Populations' (Geneva) (2022).
37. Kavanagh MM, Agbla SC, Joy M, Aneja K, Pillinger M, Case A, et al. Law, criminalisation and HIV in the world: have countries that criminalise achieved more or less successful pandemic response? BMJ Global Health. (2021) 6:e006315. doi: 10.1136/bmjgh-2021-006315
38. Ward Z, Stone J, Bishop C, Ivakin V, Eritsyan K, Deryabina A, et al. Costs and impact on HIV transmission of a switch from a criminalisation to a public health approach to injecting drug use in eastern Europe and central Asia: a modelling analysis. Lancet HIV. (2022) 9:e42–53. doi: 10.1016/S2352-3018(21)00274-5
Keywords: HIV, key populations, recombinant, CRF, URF, molecular epidemiology
Citation: Nchinda N, Elangovan R, Yun J, Dickson-Tetteh L, Kirtley S, Hemelaar J and WHO-UNAIDS Network for HIV Isolation and Characterisation (2023) Global associations of key populations with HIV-1 recombinants: a systematic review, global survey, and individual participant data meta-analysis. Front. Public Health 11:1153638. doi: 10.3389/fpubh.2023.1153638
Received: 29 January 2023; Accepted: 28 June 2023;
Published: 27 July 2023.
Edited by:
Malik Sallam, The University of Jordan, JordanReviewed by:
Christopher Williams, State University of New York, United StatesAugust Guang, Brown University, United States
Copyright © 2023 Nchinda, Elangovan, Yun, Dickson-Tetteh, Kirtley, Hemelaar and WHO-UNAIDS Network for HIV Isolation and Characterisation. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Joris Hemelaar, joris.hemelaar@ndph.ox.ac.uk