Skip to main content

Spatiotemporal dynamics of HIV-1 transmission in France (1999–2014) and impact of targeted prevention strategies



Characterizing HIV-1 transmission networks can be important in understanding the evolutionary patterns and geospatial spread of the epidemic. We reconstructed the broad molecular epidemiology of HIV from individuals with primary HIV-1 infection (PHI) enrolled in France in the ANRS PRIMO C06 cohort over 15 years.


Sociodemographic, geographic, clinical, biological and pol sequence data from 1356 patients were collected between 1999 and 2014. Network analysis was performed to infer genetic relationships, i.e. clusters of transmission, between HIV-1 sequences. Bayesian coalescent-based methods were used to examine the temporal and spatial dynamics of identified clusters from different regions in France. We also evaluated the use of network information to target prevention efforts. Participants were mostly Caucasian (85.9%) and men (86.7%) who reported sex with men (MSM, 71.4%). Overall, 387 individuals (28.5%) were involved in clusters: 156 patients (11.5%) in 78 dyads and 231 participants (17%) in 42 larger clusters (median size: 4, range 3–41). Compared to individuals with single PHI (n = 969), those in clusters were more frequently men (95.9 vs 83%, p < 0.01), MSM (85.8 vs 65.6%, p < 0.01) and infected with CRF02_AG (20.4 vs 13.4%, p < 0.01). Reconstruction of viral migrations across time suggests that Paris area was the major hub of dissemination of both subtype B and CRF02_AG epidemics. By targeting clustering individuals belonging to the identified active transmission network before 2010, 60 of the 143 onward transmissions could have been prevented.


These analyses support the hypothesis of a recent and rapid rise of CRF02_AG within the French HIV-1 epidemic among MSM. Combined with a short turnaround time for sample processing, targeting prevention efforts based on phylogenetic monitoring may be an efficient way to deliver prevention interventions but would require near real time targeted interventions on the identified index cases and their partners.


Identifying and monitoring HIV-1 transmission networks can be important in understanding the evolutionary patterns and geospatial spread of the epidemic. Recent advances in molecular epidemiology have greatly enhanced our ability to characterize the dynamics and the structure of HIV transmission networks over space and time using HIV-1 pol sequences generated for routine [13]. A better understanding of the dynamics of the HIV-1 epidemic can assist preventive measures [25]. In the past decade, there has been an increase in the circulation of non-B strains and Circulating Recombinant Forms (CRFs) of HIV-1 in Europe and North America. In France and other European countries, HIV-1 subtype B is still predominant but the proportion of non-B infected individuals has progressively increased [68]. This increase of non-B viral infections has been reported in both newly diagnosed chronic HIV-1 infections [9, 10] and individuals with primary or recent HIV-1 infection (PHI) [1114]. Among HIV-1 non-B subtypes, CRF02_AG is one of the most prevalent recombinant forms in the world, responsible for at least 8% of total infections [15]. Though CRF02_AG is predominantly transmitted within heterosexuals in Sub-Saharan Africa, it has been increasingly reported among men who have sex with men (MSM) [16]. In a recent study conducted in newly diagnosed patients living in Europe, the proportion of circulating recombinant form CRF02_AG increased significantly between 2002 and 2010 [17]. A similar trend was also observed in Western and Central Europe, with an increased proportion of CRF02_AG from 5% in 2000–2003 to 8% in 2004–2007 [15]. In France, recent data showed a spread of non-B subtypes in individuals of French origin and that the MSM group are particularly involved in this dynamic [14]. Altogether, these reports emphasize the need for a better understanding of the spread of various HIV-1 subtypes through these transmission networks. The recent advances in molecular epidemiology have greatly enhanced our ability to evaluate the dynamic of these transmission networks [2, 1821]. Several recent studies have also used clustering approaches to characterize potential correlates of HIV transmission [3, 4, 22] and there is a growing interest in using these approaches to implement and evaluate prevention interventions [1, 2, 23, 24].

In this study, we analyzed HIV-1 pol sequences generated over a period of more than 15 years from individuals enrolled during PHI in France in the National Agency for Research on AIDS and hepatitis (ANRS) PRIMO Cohort to reconstruct the broad molecular epidemiology of the HIV epidemic in France. We then determined if clustering analyses could be used efficiently to target prevention interventions in newly diagnosed PHI individuals belonging to an active cluster of transmission (index cases).


Study population

The study protocol was approved by the Paris Cochin Ethics Committee, and all patients gave their written informed consent.

The multicenter ANRS CO6 PRIMO cohort has enrolled in France more than 1900 participants with PHI since November 1996, which were diagnosed on the basis of a negative or incomplete Western blot (no anti-p68 or anti-p34) with detectable HIV-1 RNA for 96% of cases or on the basis of an interval of <3–6 months between a negative and a positive enzyme-linked immunosorbent assay (ELISA) for the remaining cases [25]. All were treatment naïve at enrollment. Demographic, behavioral, biological (CD4, HIV-RNA) and clinical data of the participants were collected and organized anonymously in a common electronic archive.

HIV-1 pol sequencing

For all participants, plasma samples were collected at enrollment in the ANRS-PRIMO cohort, centralized in the Virology Laboratory of Necker Hospital and stored for genotypic resistance testing. HIV RNA was extracted and pol-amplified products were sequenced using published primers (, HXB2 coordinates 2530–3334) [6, 26].

Transmission network

Sequence curation, alignment, and network inference were performed using the freely available software (, After quality control procedures no contaminant sequences were identified [27], and the partial transmission network was inferred based on the nucleotide genetic distances between bulk HIV-1 pol sequences from each participant [28].

Similar to previous studies investigating the structure and dynamics of HIV transmission networks [2, 1821], we linked two individuals (nodes) in networks whenever their pol sequences were ≤1.5% distant (TN93 distance measure). This relatively conservative genetic distance cutoff was determined based on previous within host evolutionary rate estimates [29] where HIV sequences from mono-infected participants collected after almost a decade had less than 1% divergence from baseline.

The degree (connectivity) of each individual was defined as the number of links (edges in the transmission network) to other individuals [30]. We also have explored the potential impact of drug resistance mutations (DRM) on network inferences by excluding codons associated with DRM [28, 31]. To prevent spurious sequence linkage due to nucleotide ambiguities, genetic distances between ambiguous nucleotides and known nucleobases were averaged (i.e., R was considered as 50% A and 50% G) [20]. Dyads and clusters were defined as connected components of the network comprising 2 nodes and >2 nodes, respectively. Singletons were defined as individuals without an identified phylogenetic connection.

Phylogenetic analysis and demographic history

HIV-1 subtypes were determined uploading sequences individually into the REGA HIV-1 automated Subtyping Tool version 2.0 ( and confirmed by in-house phylogenetic analysis.

The age of the most recent common ancestor (TMRCA, years) and the ancestral geographic movements were jointly estimated using a Markov Chain Monte Carlo (MCMC) framework as implemented in BEAST v1.8.1 [32]. We used a discretized gamma distribution (GTR + 4Γ) to account for among-site rate variation. Time scales of the trees were calibrated with the sampling dates available. The temporal scale of evolutionary process was estimated from the sampling dates of the sequences using a relaxed uncorrelated lognormal molecular clock model and a gamma prior on clock rate. Different parametric demographic models (a constant population size, exponential and logistic growth) and a nonparametric Bayesian skyline plot (BSP) were compared, and the best models were selected as confirmed by a higher Bayes Factor (BF) support implemented in BEAST [32]. MCMC simulations were run for 50–200 × 106 chain steps, sub-sampling parameters every 20,000 steps. After removing 10% of burn-in and combining evolutionary parameters and trees using LogCombiner. Convergence of the chains was inspected using Tracer.v.1.5. The TMRCA estimates were expressed as mean and 95% highest posterior density (HPD) years before the most recent sampling dates.

The identification of significant migration pathways was performed using discrete non-reversible diffusion models and a Bayesian stochastic search variable selection (BSSVS) approach [33]. We first applied a discrete diffusion model and geographic locations were recorded at the tips of pol phylogenies. To quantify the dissemination process, we estimated the number of viral migrations among locations using ‘Markov Jump’ counts [34] of location-state transitions along the posterior tree distribution [35]. In an attempt to maximize spatial information and put spatiotemporal dynamics in a demographic context, we compiled location information into 10 equally populated areas (regions 1–10). We also included sequences originating from French overseas departments (region 11).

Statistical analyses

All available data including demographics, HIV risk factor, baseline CD4 count, viral load and HIV-1 subtype or CRF were compared between individuals who clustered and those who did not and analyzed to determine if these factors changed over time. Categorical variables were compared using the Fisher exact test, and continuous variables were evaluated using the Wilcoxon rank test. p values <0.05 were considered statistically significant. The calculations of all statistical tests were performed by using Graph-Pad Prism 6.0c software (GraphPad Software, Inc., San Diego, CA) and the computing R environment.

Impact of targeted prevention

We evaluated whether phylogenetic monitoring can be used to efficiently target prevention interventions [5]. Here, we evaluated the impact of a prevention intervention targeting clustering individuals (index cases) belonging to an active cluster before 2010 on onward HIV transmission within this network. Briefly, we estimated the number of infections that would have been prevented from 2010 to 2014 if identified clustering individuals (index cases) enrolled in the PRIMO Cohort had received the prevention intervention based on the network density (defined as \(\frac{{{\sum }{\text{degrees}}\left( {\text{cluster}} \right)}}{{{\text{size}}\left( {\text{cluster}} \right)}}\)) prior to 2010. We then considered that an onward transmission would have been prevented if a participant (contact of index cases) was diagnosed after the date of intervention (i.e. from 2010 to 2014) and belonged to an active transmission cluster targeted by the intervention.


Population characteristics

A total of 1356 individuals with PHI enrolled in the ANRS PRIMO cohort CO6 between 1999 and 2014 were included. Participants were preferentially Caucasian (85.9%), male (86.7%) and MSM (71.4%). More than one-third (37%) were diagnosed in Paris (n = 348) and its surrounding areas (n = 167). The median age was 35 years (range from 17 to 79 years), and a vast majority (69%) was between 20 and 40 years old. Over this period, the proportion of HIV-infected individuals reporting MSM risk increased over the study period from 66% (n = 309) in 1999–2005 to 77.7% (n = 310) in 2011–2014 (Chi square test for trend, p < 0.01). The socio-demographic characteristics of 1356 individuals are summarized in Table 1.

Table 1 Population characteristics at primary infection

Subtype epidemics

Overall, the HIV-1 subtype B (72.5%, n = 983) was largely predominant. Among the non-B viral strains, CRF02_AG was the most prevalent (56% of non-B, n = 209) (Table 1). MSM individuals were significantly more likely to be infected with subtype B (80%) than CRF02_AG (12.2%) and other non-B lineages (p < 0.001). Among heterosexuals, subtype B clade was also the most prevalent (42%) but prevalence of CRF02_AG (27%) and other non-B subtypes (30.7%) increased overtime. We also showed a significant increase of non-B subtypes among MSM over the study period (p < 0.001).

Transmission network characteristics

The HIV-1 reverse transcriptase sequences generated from each participant were used to infer the transmission network. Overall, the mean genetic distance between pairs of sequences was 6.5% (s.d. ± 2%) and 5.9% (s.d. ± 2%) among subtype B and CRF02_AG infections, respectively. We used a pairwise distance below the threshold of 1.5% to define a link between individuals [2]. Using this threshold, 387 individuals (28.5%) were connected to at least one other study participant. Exclusion of DRM codons did not impact the observed HIV transmission network. Connected nodes (Fig. 1) were arranged in 44 large clusters (i.e., connected to more than one other individual) ranging in size from 3 to 41 individuals (median = 4) and 78 dyads (i.e., connected with a unique participant). Then, we evaluated the factors associated with clustering. Comparison between people in clusters (i.e. connected) versus people not in cluster (i.e. singletons) revealed that individuals who are linked in a cluster were more likely men (95.9 vs 83%, p < 0.01) and significantly younger (median age = 32.5 vs 36.0 years, p < 0.01) than singletons. MSM (85.8% clustering vs 65.6% non-clustering, p < 0.01) and white participants (91.5 vs 83.7%, p < 0.01) were significantly more likely to cluster. There were no significant associations between clustering and: baseline CD4 T cell count (p = 0.25), baseline viral load (p = 0.15) or being diagnosed in Paris and its suburban area (37.0 vs 38.4%, p = 0.66) (Table 1). Individuals infected with CRF02_AG viruses were also more likely to be connected to another participant than other subtype infections (20.4 vs 13.4%) (Additional file 1: Figure S1). As might be expected, the probability that an individual was connected with another individual increased significantly during the study period (from 1999 to 2014) (p < 0.01).

Fig. 1

Inferred HIV transmission clusters. HIV-1 transmission cluster diagrams illustrating the structure and demographics of the putative transmission clusters identified in the PRIMO ANRS CO6 cohort. A total of 387 of the 1356 (28.5%) individuals were connected with at least one other individual. Color indicates the reported transmission risk [red MSM; green heterosexual (HTS), purple others]; and shape denotes gender (ellipse male, square women). All edges represent a genetic distance of ≤1.5% separating nodes. All shapes are labeled according to the HIV-1 subtype. NA not available. White and unfilled dots correspond to missing informations

Within the transmission network, three distinctly large viral clusters B1 (subtype B), AG1 and AG2 (CRF02_AG) were identified with a respective size of 9, 14 and 41. Cluster B1 included nine white MSM (median age: 32 years old [21–37]). It was first identified in 2001 with no additional cases since 2010, mostly originating from the Nantes area (West coast) (Fig. 2). Cluster AG1 was composed of fourteen white MSM (median age: 31 years old [20–48]) infected with CRF02_AG-related variants. It was first identified in 2011 and had a peak identification of individuals in 2013. Individuals in this cluster were mostly originating from the Marseille area (South-East) but more widely spread throughout the country (Fig. 2). Finally, cluster AG2 was composed of forty-one white men (median age: 32 years old [21–45]) mostly originating from Paris and its suburban area (n = 34, 83%). It was first noted in 2000 and accrued individuals relatively steadily throughout the observation period (Fig. 2). All but one reported MSM sexual risk exposure. To further evaluate the timing of these clusters, we found that the cluster B1 had a TMRCA in 1997 (95% HPD: 1995–2000) (Fig. 3), and the CRF02_AG clusters AG1 and AG2 had a TMRCA in 2003 (95% HPD: 2000–2008) and 2000 (95% HPD: 1998–2001) respectively (Fig. 3).

Fig. 2

Characteristics of the 3 larger clusters B1, AG1 and AG2. a Transmission network of the three larger cluster AG1 (n = 14), AG2 (N = 41) and B1 (n = 9) and evolution of the main clusters over the study period. b Map representing the number of clustering individuals by location of residence. c Ancestral root state probabilities. The root state probabilities are presented with the color codes corresponding to the 11 equally populated regions

Fig. 3

Bayesian time-scaled tree of the HIV transmission network of subtype B (a) and CRF02_AG (b) pol sequences in clusters from the participants enrolled in the PRIMO ANRS cohort between 1999 and 2014. Time scaled in year. Nodes and branches are colored according to the most probable location state of their descendent nodes. Tips are colored according to the recorded location of sampling

Spatiotemporal patterns of the subtype B and CRF02_AG epidemics

Discrete diffusion models can offer insights into the origins and epidemiological links within the set of locations from which infections were sampled [33, 36].

To analyze the spatial spread of subtype B and CRF02_AG epidemic among individuals with PHI in France between 1999 and 2014 at a national scale, we compiled all subtype B (n = 983) and CRF02_AG (n = 209) sequences. Given the demographic characteristics of France with wide range of population density, we performed our analyzes after dividing the country in 10 equally populated area and an 11th area accounting for individuals originating from overseas departments (n = 86) accounting for 6.3% of the participants.

Geographic locations throughout the phylogenetic histories were estimated by applying a discrete asymmetric Bayesian phylogeographic approach, which allows for possibly different dispersal rates between two locations depending on the directionality of diffusion [37]. Reconstruction of viral migrations across time suggests that Paris area was the major hub of dissemination of both subtype B and CRF02_AG epidemics.

Similarly, we evaluated the dynamics of the three main clusters of subtype B and CRF02_AG lineages. The most probable root location for phylogeny of cluster B1 was the Western region R3 with posterior state probability of 0.72. The most probable root state probabilities were South Eastern regions with posterior state probability and Paris region R1 with posterior state probability of 0.95 for the AG1 and AG2 clusters, respectively (Fig. 2). These results showed that the diffusion of the larger CRF02_AG cluster (AG2) ignited in the most densely populated region (Paris, R1) and thereafter spread to other areas (Fig. 2).

Targeted ART for prevention

Given the depth of sampling and the characteristics of the ANRS PRIMO transmission network, we evaluated the number of onward transmissions that could have been prevented by an enhanced prevention intervention targeting newly diagnosed clustering HIV-infected individuals (index cases).

We identified that 223 of the 872 (25.6%) newly diagnosed individuals before 2010 in the PRIMO cohort belonged to a transmission network. We estimated that an enhanced prevention intervention on these clustering index individuals could have potentially prevented 60 of the 143 (42%) onward new infections linked to these index cases. The 320 individuals diagnosed between 2010 and 2014 that did not cluster with any individuals infected before 2010 would not have been protected by any prevention intervention given to prior cohort participants. We next evaluated the impact of intervention targeted to participants enrolled before 2010 by cluster densities defined as \(\frac{{{\sum }{\text{degrees}}\left( {\text{cluster}} \right)}}{{{\text{size}}\left( {\text{cluster}} \right)}}.\) We found that intensified prevention intervention provided to all individuals belonging to clusters with densities ≥2 (n = 79 individuals) diagnosed before 2010 would have potentially prevented 33 transmissions linked to these index cases.


Here, we investigated the transmission network and spatiotemporal dynamics of HIV-1 epidemic in 1356 primary infected individuals enrolled in the French ANRS PRIMO CO6 cohort between 1999 and 2014. We particularly focused on the major circulating subtypes B and CRF02_AG to provide a better understanding on the dynamics of these epidemics in France and their relative evolution [6, 12, 26, 3840].

As the epidemic has matured, patterns of HIV transmission have changed, albeit MSM represent a large proportion of the affected populations in many countries [41]. In Western Europe, the incidence of HIV-1 among MSM has increased or at least has remained constant over the last decade [4244]. The increased proportion of MSM is also observed among all new HIV diagnoses in France (up to 43% in 2013) [45]. Consistent with these observations, we found that the proportion of MSM among all PHI enrolled in the ANRS-PRIMO cohort significantly increased over the last 15 years (from 66% in 1999–2005 to 77.7% in 2011–2014). However, the high frequency of MSM among recently HIV-infected patients can also be partially explained by targeted and repeated screening of MSM populations. For example in the United Kingdom, the proportion of individuals tested for HIV among people attending to STI clinics was the highest among MSM and increased overtime (86% in 2013 vs 78% in 2009) [46]. More recently, a multicenter preventive trial in France (ANRS IPERGAY) also showed a high incidence rate of HIV among MSM, up to 9 per 100 person-years in Paris [47].

By combining methods from classical and molecular epidemiology, we were able to infer and characterize the HIV-1 transmission networks among identified PHI individuals in France over 15 years. Similar to previous data showing that the HIV-1 networks derived mostly from populations of MSM [2, 39, 48], we observed that MSM were significantly more prevalent in connected individuals than singletons. We also showed that non-B subtypes were frequent in primary infected individuals (28.5%), although slightly less than among all new HIV diagnoses in France [10, 45]. We observed that 57% of the individuals infected with CRF02_AG variants reported MSM sexual risk exposure. These results are in line with reports from other European countries showing an increase in proportion of non-B clades in the MSM population [49] suggesting that the sociodemographic boundaries between HIV-1 subtypes globally are diminishing in Western Europe.

By applying Bayesian phylogeographic inference using discrete non-reversible models to pol geo-referenced sequences, we also investigated the spatial patterns of subtype B and CRF02_AG clades in France among PHI individuals. We found that these two distinct epidemic lineages have ignited in the most urbanized region of Paris (R1, “Ile de France”) as illustrated by the most probable location state of the descendant nodes with over 90% of all viral lineages movements originating from Paris area (i.e. viral dispersal from Paris) for both lineages (96.5% [95% CI 95.1–97.5%] and 94.1% [92.2–95.6%] for subtype B and CRF02_AG respectively). These results suggest the key role of Paris as a hub for new HIV infections and the potential spatial expansion from this region to the rest of the country. This might be explained by the high prevalence of MSM in the ANRS-PRIMO cohort, who are more likely originating from urbanized and highly-densified area [50].

Interestingly, the three larger clusters of HIV-1 subtype B (cluster B1, n = 9) and CRF02_AG (clusters AG1, n = 14 and AG2, n = 41) identified displayed also very distinct evolutionary patterns. While cluster B1 did not increase since 2010, both clusters AG1 and AG2 continued to grow up to end of the follow up period in late 2014. Altogether, these observations suggest that CRF02_AG clade was introduced and disseminated within highly connected networks of MSM, which may explain the successful rapid dissemination and increased prevalence of CRF02_AG subtype [51].

Understanding the dynamic of HIV transmission is crucial in the design of effective interventions and recently Individuals contributes disproportionately to the spread of the HIV epidemic [5254]. Considering the limited time frame of HIV transmission, targeted prevention strategies focusing on PHI may have a significant impact on the HIV epidemic [5255]. Here, the extensive collection of samples along with demographic and clinical data of the 1356 PHI participants from the PRIMO ANRS cohort allowed us to better time the observed transmission clusters. It also helped to better evaluate the direct effect of theoretical prevention interventions at an individual level rather than relying on population effect of the intervention (i.e. indirect effect).

Given that the French 2013 guidelines for ART of HIV-1 infection in adults [56] recommend that ART should be initiated in any HIV-positive person, whatever his/her CD4 T cell count, even when >500/mm3, we first considered a prevention intervention strategy based on immediate ART introduction [57]. We found that an intensified intervention targeting all clustering participants with PHI diagnosed before 2010 (n = 223 index cases) would have prevented 60 out of the identified 143 onward transmissions identified between 2010 and 2014. Given the depth of sampling and structure of the transmission network, we evaluated intervention strategies based on network connectivity. We found that prevention targeted to newly infected individuals belonging to intermediate and high density clusters (density ≥ 2, n = 79 index cases) would have potentially prevented 33 onward new infections linked to these index cases. In light with a recent study showing that near real-time phylogenetic monitoring of routinely collected HIV genotypes are a promising resource for public health intervention on localized outbreaks of HIV transmission [23], our results also emphasize the need of an efficient sample processing, a rapid turnaround time on HIV sequence generation and a near real time based monitoring of clustering analysis. Beyond the rapid identification of the index cases, an effective and efficient prevention strategy would require a combination of interventions targeting these indices (i.e. immediate provision of ART and enhanced adherence support services), and the identification and immediate provision of pre-exposure prophylaxis (PrEP) to the contacts of these index cases [58]. Identifying and monitoring HIV clusters should be also an invaluable method to (1) determine the leading edge of local HIV transmission, (3) characterize potential correlates of HIV transmission (i.e. recreational drug use, high risk venues) [3, 4, 52], (3) investigate hotspots of transmission [23] and (4) further intensify HIV screening, providing PrEP and partner tracing [59] in population with high rates of transmission.

This study has a number of limitations. First, HIV clustering and network inferences are directly affected by the sampling density [60]: studies with low sampling density showed minimal HIV clustering [61], while in depth sampling allows more accurate characterization of HIV transmission network [48, 49, 62]. In these deeper sampling density studies, the proportion of clustering sequences varied between 28 and 41%, though these studies were performed in heterogeneous population, in other countries and with various methods. While the number of individuals in our dataset is smaller than other studies, we reported an overall clustering proportion of 28.5% and up to 52.5% (332 out of 636) among MSM, consistent with previous reports. This overall clustering rate is also consistent with previous reports where interventions have been deemed useful [5] and allowed the characterization of the dynamics of HIV transmission [19]. We have also estimated that the PRIMO ANRS Cohort is well representative of new HIV infections occurring in France. Indeed, based on the data recorded through the mandatory notification of new HIV infections, 12% of new HIV diagnoses are done at time of primary infection, corresponding to 600–700 cases each year [63]. Therefore, the cases included in the PRIMO ANRS cohort represent approximately 15% of all primary infections diagnosed in the country. In addition, the proportion of MSM included in our study (71.4%) is similar to that observed at the national level (73.6%) [64]. Altogether, while our dataset is not ideal and that we could have missed a number of HIV infected individuals in the local network due to the lack of diagnosis of all PHI, it is sufficient to use the presented techniques to understand the transmission dynamics underlying the sampled epidemic.

Second, recent HIV infection accounts for 39% of new HIV diagnoses in France [45]. Here, we did not take into account the contribution of transmissions from chronically infected individuals to partners diagnosed after the PHI stage, which is also likely to have an important impact on epidemic spread [65] and rely on convenient sampling. Hence, individuals who are not tested or disengaged from care are not represented in this cohort and may also contribute to the epidemic. Though the number of early diagnoses has increased since 2011 in France, the number of late diagnoses has remained stable among MSM, and 25% of the estimated 6220 new HIV diagnoses in 2013 in France were still diagnosed at a late stage (<200 CD4 or AIDS) [45]. Another limitation of all convenience based sampling methods is that the presence of a directed link between two individuals simply reflects recent relatedness of the virus, possibly through a series of unobserved intermediaries [2]. Finally, the use of HIV pol region for HIV subtyping may have led to potential misclassification of CRF forms. However, it should not impact the HIV network inferences and the dynamic of CRF02_AG clade within the HIV transmission network presented in this study.


In summary, this study showed the rapid evolutionary dynamics of subtype B and CRF02_AG among PHI in France since the late 1990s. While subtype B remains the most prevalent lineage, the rapid diffusion of CRF02_AG in highly connected networks of MSM could lead to a substantial and rapid reshaping of the HIV epidemic in Western Countries. Combined with a short turnaround time for sample processing, our findings also show that identifying hotspots of HIV transmission and near real-time monitoring based on phylogenetic analyses can be an effective prevention intervention in combination with Public Health Recommendations for early treatment introduction, enhanced adherence support and partner tracing for immediate diagnosis, provision of PrEP and counseling.


  1. 1.

    Smith DM, May SJ, Tweeten S, Drumright L, Pacold ME, Kosakovsky Pond SL, et al. A public health model for the molecular surveillance of HIV transmission in San Diego, California. AIDS. 2009;23(2):225–32.

    Article  PubMed  PubMed Central  Google Scholar 

  2. 2.

    Little SJ, Kosakovsky Pond SL, Anderson CM, Young JA, Wertheim JO, Mehta SR, et al. Using HIV networks to inform real time prevention interventions. PLoS ONE. 2014;9(6):e98443.

    Article  PubMed  PubMed Central  Google Scholar 

  3. 3.

    Wertheim JO, Leigh Brown AJ, Hepler NL, Mehta SR, Richman DD, Smith DM, et al. The global transmission network of HIV-1. J Infect Dis. 2014;209(2):304–13.

    Article  PubMed  Google Scholar 

  4. 4.

    Poon AF, Joy JB, Woods CK, Shurgold S, Colley G, Brumme CJ, et al. The impact of clinical, demographic and risk factors on rates of HIV transmission: a population-based phylogenetic analysis in British Columbia, Canada. J Infect Dis. 2015;211(6):926–35.

    Article  PubMed  Google Scholar 

  5. 5.

    Wang X, Wu Y, Mao L, Xia W, Zhang W, Dai L, et al. Targeting HIV prevention based on molecular epidemiology among deeply sampled subnetworks of men who have sex with men. Clin Infect Dis. 2015;61(9):1462–8.

    Article  PubMed  PubMed Central  Google Scholar 

  6. 6.

    Chaix M-L, Descamps D, Harzic M, Schneider V, Deveau C, Tamalet C, et al. Stable prevalence of genotypic drug resistance mutations but increase in non-B virus among patients with primary HIV-1 infection in France. AIDS. 2003;17(18):2635–43.

    Article  PubMed  Google Scholar 

  7. 7.

    Descamps D, Chaix M-L, André P, Brodard V, Cottalorda J, Deveau C, et al. French national sentinel survey of antiretroviral drug resistance in patients with HIV-1 primary infection and in antiretroviral-naive chronically infected patients in 2001–2002. J Acquir Immune Defic Syndr. 2005;38(5):545–52.

    Article  PubMed  Google Scholar 

  8. 8.

    Barin F, Courouce AM, Pillonel J, Buzelay L. Increasing diversity of HIV-1 M serotypes in French blood donors over a 10-year period (1985–1995). Retrovirus Study Group of the French Society of Blood Transfusion. AIDS. 1997;11(12):1503–8.

    CAS  Article  PubMed  Google Scholar 

  9. 9.

    Semaille C, Barin F, Cazein F, Pillonel J, Lot F, Brand D, et al. Monitoring the dynamics of the HIV epidemic using assays for recent infection and serotyping among new HIV diagnoses: experience after 2 years in France. J Infect Dis. 2007;196(3):377–83.

    Article  PubMed  Google Scholar 

  10. 10.

    Descamps D, Assoumou L, Chaix ML, Chaillon A, Pakianather S, de Rougemont A, et al. National sentinel surveillance of transmitted drug resistance in antiretroviral-naive chronically HIV-infected patients in France over a decade: 2001–2011. J Antimicrob Chemother. 2013;68(11):2626–31.

    CAS  Article  PubMed  Google Scholar 

  11. 11.

    Chaix ML, Descamps D, Wirden M, Bocket L, Delaugerre C, Tamalet C, et al. Stable frequency of HIV-1 transmitted drug resistance in patients at the time of primary infection over 1996–2006 in France. Aids. 2009;23(6):717–24.

    PubMed  Google Scholar 

  12. 12.

    Chaix ML, Seng R, Frange P, Tran L, Avettand-Fenoel V, Ghosn J, et al. Increasing HIV-1 non-B subtype primary infections in patients in France and effect of HIV subtypes on virological and immunological responses to combined antiretroviral therapy. Clin Infect Dis. 2013;56(6):880–7.

    CAS  Article  PubMed  Google Scholar 

  13. 13.

    Frange P, Galimand J, Vidal N, Goujard C, Deveau C, Souala F, et al. New and old complex recombinant HIV-1 strains among patients with primary infection in 1996–2006 in France: the French ANRS CO06 primo cohort study. Retrovirology. 2008;5:69.

    Article  PubMed  PubMed Central  Google Scholar 

  14. 14.

    Brand D, Moreau A, Cazein F, Lot F, Pillonel J, Brunet S, et al. Characteristics of patients recently infected with HIV-1 non-B subtypes in France: a nested study within the mandatory notification of new HIV diagnoses. J Clin Microbiol. 2014;52(11):4010–6.

    Article  PubMed  PubMed Central  Google Scholar 

  15. 15.

    Hemelaar J, Gouws E, Ghys PD, Osmanov S. Global trends in molecular epidemiology of HIV-1 during 2000–2007. AIDS. 2011;25(5):679–89.

    Article  PubMed  PubMed Central  Google Scholar 

  16. 16.

    Ndiaye HD, Tchiakpe E, Vidal N, Ndiaye O, Diop AK, Peeters M, et al. HIV type 1 subtype C remains the predominant subtype in men having sex with men in Senegal. AIDS Res Hum Retroviruses. 2013;29(9):1265–72.

    PubMed  Google Scholar 

  17. 17.

    Hofstra LM, Sauvageot N, Albert J, Alexiev I, Garcia F, Struck D, et al. Transmission of HIV drug resistance and the predicted effect on current first-line regimens in Europe. Clin Infect Dis. 2016;62(5):655–63.

    Article  PubMed  Google Scholar 

  18. 18.

    Hoenigl M, Chaillon A, Kessler HH, Haas B, Stelzl E, Weninger K, et al. Characterization of HIV transmission in South-East Austria. PLoS ONE. 2016;11(3):e0151478.

    Article  PubMed  PubMed Central  Google Scholar 

  19. 19.

    Mehta SR, Wertheim JO, Brouwer KC, Wagner K, Chaillon A, Strathdee SA, et al. HIV transmission networks in the San Diego-Tijuana border region. EBioMedicine. 2015;2(10):1456–63.

    Article  PubMed  PubMed Central  Google Scholar 

  20. 20.

    Wertheim JO, Oster AM, Hernandez AL, Saduvala N, Banez Ocfemia MC, Hall HI. The international dimension of the U.S. HIV transmission network and onward transmission of HIV recently imported into the United States. AIDS Res Hum Retroviruses. 2016;32(10–11):1046–53.

    Article  PubMed  Google Scholar 

  21. 21.

    Whiteside YO, Song R, Wertheim JO, Oster AM. Molecular analysis allows inference into HIV transmission among young men who have sex with men in the United States. AIDS. 2015;29(18):2517–22.

    Article  PubMed  PubMed Central  Google Scholar 

  22. 22.

    Brenner BG, Roger M, Stephens D, Moisi D, Hardy I, Weinberg J, et al. Transmission clustering drives the onward spread of the HIV epidemic among men who have sex with men in Quebec. J Infect Dis. 2011;204(7):1115–9.

    Article  PubMed  PubMed Central  Google Scholar 

  23. 23.

    Poon AF, Gustafson R, Daly P, Zerr L, Demlow SE, Wong J, et al. Near real-time monitoring of HIV transmission hotspots from routine HIV genotyping: an implementation case study. The lancet HIV. 2016;3(5):e231–8.

    Article  PubMed  Google Scholar 

  24. 24.

    Ragonnet-Cronin M, Ofner-Agostini M, Merks H, Pilon R, Rekart M, Archibald CP, et al. Longitudinal phylogenetic surveillance identifies distinct patterns of cluster dynamics. J Acquir Immune Defic Syndr. 2010;55:102–8.

    Article  PubMed  Google Scholar 

  25. 25.

    Goujard C, Bonarek M, Meyer L, Bonnet F, Chaix M-L, Deveau C, et al. CD4 cell count and HIV DNA level are independent predictors of disease progression after primary HIV type 1 infection in untreated patients. Clin Infect Dis. 2006;42(5):709–15.

    CAS  Article  PubMed  Google Scholar 

  26. 26.

    Ghosn J, Pellegrin I, Goujard C, Deveau C, Viard JP, Galimand J, et al. HIV-1 resistant strains acquired at the time of primary infection massively fuel the cellular reservoir and persist for lengthy periods of time. AIDS. 2006;20(2):159–70.

    Article  PubMed  Google Scholar 

  27. 27.

    Smith D, Delport W, Butler D, Little S, Richman D, Pond SK. Response to comment on “The origins of sexually transmitted HIV among men who have sex with men”. Sci Transl Med. 2010;2(50):501r1.

    Article  PubMed  PubMed Central  Google Scholar 

  28. 28.

    Lewis F, Hughes GJ, Rambaut A, Pozniak A, Leigh Brown AJ. Episodic sexual transmission of HIV revealed by molecular phylodynamics. PLoS Med. 2008;5(3):e50.

    Article  PubMed  PubMed Central  Google Scholar 

  29. 29.

    Hightower GK, May SJ, Pérez-Santiago J, Pacold ME, Wagner GA, Little SJ, et al. HIV-1 clade B pol evolution following primary infection. PLoS ONE. 2013;8(6):e68188.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  30. 30.

    Helleringer S, Kohler HP. Sexual network structure and the spread of HIV in Africa: evidence from Likoma Island, Malawi. AIDS. 2007;21(17):2323–32.

    Article  PubMed  Google Scholar 

  31. 31.

    Wheeler WH, Ziebell RA, Zabina H, Pieniazek D, Prejean J, Bodnar UR, et al. Prevalence of transmitted drug resistance associated mutations and HIV-1 subtypes in new HIV-1 diagnoses, U.S.-2006. AIDS. 2010;24(8):1203–12.

    Article  PubMed  Google Scholar 

  32. 32.

    Drummond AJ, Suchard MA, Xie D, Rambaut A. Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol. 2012;29(8):1969–73.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  33. 33.

    Lemey P, Rambaut A, Drummond AJ, Suchard MA. Bayesian phylogeography finds its roots. PLoS Comput Biol. 2009;5(9):e1000520.

    Article  PubMed  PubMed Central  Google Scholar 

  34. 34.

    Minin VN, Suchard MA. Counting labeled transitions in continuous-time Markov models of evolution. J Math Biol. 2008;56(3):391–412.

    Article  PubMed  Google Scholar 

  35. 35.

    Nunes MRT, Faria NR, Vasconcelos HB, Medeiros D, Silva de Lima CP, Carvalho VL, et al. Phylogeography of dengue virus serotype 4, Brazil, 2010-2011. Emerg Infect Dis. 2012;18(11):1858–64.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  36. 36.

    Faria NR, Suchard MA, Rambaut A, Lemey P. Toward a quantitative understanding of viral phylogeography. Curr Opin Virol. 2011;1(5):423–9.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  37. 37.

    Drummond AJ, Rambaut A. BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol Biol. 2007;7:214.

    Article  PubMed  PubMed Central  Google Scholar 

  38. 38.

    Frange P, Assoumou L, Descamps D, Cheret A, Goujard C, Tran L, et al. HIV-1 subtype B-infected MSM may have driven the spread of transmitted resistant strains in France in 2007–12: impact on susceptibility to first-line strategies. J Antimicrob Chemother. 2015;70(7):2084–9.

    CAS  PubMed  Google Scholar 

  39. 39.

    Frange P, Meyer L, Deveau C, Tran L, Goujard C, Ghosn J, et al. Recent HIV-1 infection contributes to the viral diffusion over the French territory with a recent increasing frequency. PLoS ONE. 2012;7(2):e31695.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  40. 40.

    Meyer L, Seng R, Allegre T, Timsit J, Talamali A, Reynes J, et al. Increasing frequency of self-reported orogenital HIV-1 transmission among men having sex with men: the ANRS PRIMO Cohort. J Acquir Immune Defic Syndr. 2013;63(5):e164–6.

    Article  PubMed  Google Scholar 

  41. 41.

    Hamers FF, Downs AM. The changing face of the HIV epidemic in Western Europe: what are the implications for public health policies? Lancet. 2004;364(9428):83–94.

    Article  PubMed  Google Scholar 

  42. 42.

    Le Vu S, Le Strat Y, Barin F, Pillonel J, Cazein F, Bousquet V, et al. Population-based HIV-1 incidence in France, 2003–08: a modelling analysis. Lancet Infect Dis. 2010;10(10):682–7.

    Article  PubMed  Google Scholar 

  43. 43.

    Semaille C, Cazein F, Lot F, Pillonel J, Le Vu S, Le Strat Y, et al. Recently acquired HIV infection in men who have sex with men (MSM) in France, 2003–2008. Euro Surveill. 2009;14(48):1424–7.

    Google Scholar 

  44. 44.

    Sullivan PS, Hamouda O, Delpech V, Geduld JE, Prejean J, Semaille C, et al. Reemergence of the HIV epidemic among men who have sex with men in North America, Western Europe, and Australia, 1996–2005. Ann Epidemiol. 2009;19(6):423–31.

    Article  PubMed  Google Scholar 

  45. 45.

    Cazein F PJ, Le Strat Y, Pinget R, Le Vu S, Brunet S, et al. New HIV and AIDS diagnoses, France, 2003–2013. Bull Epidémiol Hebd. 2015;9–10:152–61.

    Google Scholar 

  46. 46.

    Yin Z BA, Hughes G, Nardone A, Gill ON, Delpech VC, et al. HIV in the United Kingdom 2014 report: data to end 2013. Public Health, London. 2014.

  47. 47.

    McCormack SM, Veronica N, Michel MJ. PrEP in Europe—expectations, opportunities and barriers. JIAS. 2016.

  48. 48.

    Leigh Brown AJ, Lycett SJ, Weinert L, Hughes GJ, Fearnhill E, Dunn DT, et al. Transmission network parameters estimated from HIV sequences for a nationwide epidemic. J Infect Dis. 2011;204(9):1463–9.

    Article  PubMed  PubMed Central  Google Scholar 

  49. 49.

    Ragonnet-Cronin M, Lycett SJ, Hodcroft EB, Hue S, Fearnhill E, Brown AE, et al. Transmission of Non-B HIV subtypes in the United Kingdom is increasingly driven by large non-heterosexual transmission clusters. J Infect Dis. 2016;213(9):1410–8.

    Article  PubMed  Google Scholar 

  50. 50.

    Giraud C. Residential spaces: looking at where homosexual men in Paris choose to live. Soc Contemp. 2011;81:151–76.

    Article  Google Scholar 

  51. 51.

    Yebra G, Kalish ML, Leigh Brown AJ. Reconstructing the HIV-1 CRF02_AG and CRF06_cpx epidemics in Burkina Faso and West Africa using early samples. Infect Genet Evol. 2016;46:209–18.

    Article  PubMed  Google Scholar 

  52. 52.

    Brenner BG, Roger M, Routy JP, Moisi D, Ntemgwa M, Matte C, et al. High rates of forward transmission events after acute/early HIV-1 infection. J Infect Dis. 2007;195:951–9.

    CAS  Article  PubMed  Google Scholar 

  53. 53.

    Hollingsworth TD, Anderson RM, Fraser C. HIV-1 transmission, by stage of infection. J Infect Dis. 2008;198(5):687–93.

    Article  PubMed  Google Scholar 

  54. 54.

    Volz EM, Ionides E, Romero-Severson EO, Brandt M-G, Mokotoff E, Koopman JS. HIV-1 transmission during early infection in men who have sex with men: a phylodynamic analysis. PLoS Med. 2013;10(12):e1001568.

    Article  PubMed  PubMed Central  Google Scholar 

  55. 55.

    Mehta SR, Murrell B, Anderson CM, Kosakovsky Pond SL, Wertheim JO, Young JA, et al. Using HIV sequence and epidemiologic data to assess the effect of self-referral testing for acute HIV infection on incident diagnoses in San Diego, California. Clin Infect Dis. 2016;63(1):101–7.

    Article  PubMed  Google Scholar 

  56. 56.

    Hoen B, Bonnet F, Delaugerre C, Delobel P, Goujard C, L’Henaff M, et al. French 2013 guidelines for antiretroviral therapy of HIV-1 infection in adults. J Int AIDS Soc. 2014;17:19034.

    Article  PubMed  PubMed Central  Google Scholar 

  57. 57.

    Granich RM, Gilks CF, Dye C, De Cock KM, Williams BG. Universal voluntary HIV testing with immediate antiretroviral therapy as a strategy for elimination of HIV transmission: a mathematical model. Lancet. 2009;373(9657):48–57.

    Article  PubMed  Google Scholar 

  58. 58.

    Molina JM, Capitant C, Spire B, Pialoux G, Cotte L, Charreau I, et al. On-demand preexposure prophylaxis in men at high risk for HIV-1 infection. N Engl J Med. 2015;373(23):2237–46.

    CAS  Article  PubMed  Google Scholar 

  59. 59.

    Green N, Hoenigl M, Chaillon A, Anderson CM, Pond SL, Smith DM, et al. Partner services in adults with acute and early HIV infection, San Diego, 1996–2014. AIDS. 2016;31(2):287–93.

    Article  Google Scholar 

  60. 60.

    Novitsky V, Moyo S, Lei Q, DeGruttola V, Essex M. Impact of sampling density on the extent of HIV clustering. AIDS Res Hum Retroviruses. 2014;30(12):1226–35.

    Article  PubMed  PubMed Central  Google Scholar 

  61. 61.

    Yirrell DL, Pickering H, Palmarini G, Hamilton L, Rutemberwa A, Biryahwaho B, et al. Molecular epidemiological analysis of HIV in sexual networks in Uganda. AIDS. 1998;12(3):285–90.

    CAS  Article  PubMed  Google Scholar 

  62. 62.

    Bezemer D, Van SA, Lukashov VV, Van Der Hoek L, Back N, Schuurman R, et al. Transmission networks of HIV-1 among men having sex with men in the Netherlands. AIDS. 2010;24:271–82.

    Article  PubMed  Google Scholar 

  63. 63.

    Cazein FLF, Pillonel J, Le Strat Y, Sommen C, Pinget R, et al. Découvertes de séropositivité VIH et sida - France, 2003–2012. Bull Epidémiol Hebd. 2014;9–10:154–62.

    Google Scholar 

  64. 64.

    Brand D, Capsec J, Chaillon A, Cazein F, Le Vu S, Moreau A, et al. HIV surveillance combining an assay for identification of very recent infection and phylogenetic analyses on dried spots. AIDS. 2016;31(3):407–16.

    Google Scholar 

  65. 65.

    Sigaloff KC, Lange JM, Montaner J. Global response to HIV: treatment as prevention, or treatment for treatment? Clin Infect Dis. 2014;59(Suppl 1):S7–s11.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Authors’ contributions

Conceived and designed the experiments: MLC, AE, LM, CR. Performed the experiments: LM, MLC, CR. Analyzed the data: AC, MLC. Contributed reagents/materials/analysis tools: MLC, LM, CG. Wrote the paper: AC, MLC, FB, AE, CD, CR, CG, LM, PF, DMS, JG, GP, OR. All authors read and approved the final manuscript.


We thank all the patients who are participating in the ANRS PRIMO cohort and the ANRS PRIMO Cohort scientific committee: V Avettand-Fenoël, F Barin, C Bourgeois, ML Chaix, A Cheret, S Couffin-Cadiergues, JF Delfraissy, A Essat, H Fischer, C Goujard, C Lascoux-Combe, C Lecuroux, L Marchand, L Meyer, C Rouzioux, A Saez-Cirion, R Seng. We also thank all physicians participating to the ANRS Primo Cohort study cohort; the full list can be found at

Competing interests

The authors declare that they have no competing interests.

Availability of data and material

The datasets during and/or analysed during the current study available from the corresponding author on reasonable request.

Ethics approval and consent to participate

The study protocol was approved by the Paris Cochin Ethics Committee, and all patients gave their written informed consent.


This study was supported by the ANRS (Agence Nationale de Recherches sur le SIDA et les Hépatites virales), Paris, France; the Department of Veterans Affairs and Grants from the National Institutes of Health AI093163, AI036214, AI027763, AI120009, MH101012, TR000098-04AI100665, MH097520, DA034978, DA041007, AI068636, MH081482, AmFAR Grant 108537, the James B. Pendleton Charitable Trust and the Bettencourt-Schueller Foundation. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Author information


Author notes

  1. Asma Essat and Pierre Frange contributed equally to this work



    Corresponding author

    Correspondence to Antoine Chaillon.

    Additional file

    Additional file 1: Figure S1.

    Contingency Plots by Risk (A) and by HIV subtype (B) among clustering and not clustering individuals. A Mosaic plot of the clustering and not clustering individuals according to the risk of HIV acquisition; B Mosaic plot of the clustering and not clustering individuals according to the HIV subtype. MSM: Man who have sex with Men; HTS: Heterosexual individual; IDU: Injection Drug User.

    Rights and permissions

    Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

    Reprints and Permissions

    About this article

    Verify currency and authenticity via CrossMark

    Cite this article

    Chaillon, A., Essat, A., Frange, P. et al. Spatiotemporal dynamics of HIV-1 transmission in France (1999–2014) and impact of targeted prevention strategies. Retrovirology 14, 15 (2017).

    Download citation


    • HIV-1
    • Transmission network
    • Phylogeography
    • Primary infection
    • Treatment as prevention