Influence of sexual risk behaviour and STI co-infection dynamics on the evolution of HIV set point viral load in MSM

Diana M Hendrickx, Wim Delva, Niel Hens
doi: https://doi.org/10.1101/19011221

Abstract

HIV viral load (VL) is an important predictor of HIV progression and transmission. Anti-retroviral therapy (ART) has been reported to reduce HIV transmission by lowering VL. However, apart from this beneficial effect, increased levels of population mean set-point viral load (SPVL), an estimator for HIV virulence, have been observed in men who have sex with men (MSM) in the decade following the introduction of ART in the Netherlands. Several studies have been devoted to explain these counter-intuitive trends in SPVL. However, to our knowledge, none of these studies has investigated an explanation in which it arises as the result of a sexually transmitted infection (STI) co-factor in detail.

In this study, we adapted an event-based, individual-based model to investigate how STI co-infection and sexual risk behaviour affect the evolution of HIV SPVL in MSM before and after the introduction of ART.

The results suggest that sexual risk behaviour has an effect on SPVL and indicate that more data are needed to test the effect of STI co-factors on SPVL. Furthermore, the observed trends in SPVL cannot be explained by sexual risk behaviour and STI co-factors only.

We recommend to develop mathematical models including also factors related to viral evolution as reported earlier in the literature. However, this requires more complex models, and the collection of more data for parameter estimation than what is currently available.

1 Introduction

HIV viral load (HIV RNA concentration in plasma) is an important and widely used prognostic marker for HIV disease progression and transmission [24, 27, 36]. While immediately after infection rapid HIV replication and high HIV viral load (VL) values are observed, VL declines during the asymptomatic phase and reaches a stable level, called set-point viral load (SPVL) after a few weeks to a few months [13]. Evolution of HIV virulence is often estimated by means of the evolution of the population mean SPVL [22].

It has been reported that anti-retroviral therapy (ART) reduces HIV transmission by decreasing VL [28, 31], and that no HIV transmission happens under successful ART. However, apart from these beneficial effects, increased population mean SPVL levels, corresponding to increased HIV virulence, were observed in MSM in the Netherlands during the period 1995-2007, after the introduction of ART in 1994 [19]. Moreover, an increase in the number of new HIV cases in MSM in the Netherlands has been observed during the period 1996-2008. After 2008, a drop in new HIV diagnoses has been observed [45].

Several studies have been devoted to explain these counter-intuitive trends in HIV SPVL after the introduction of ART. The majority of these studies try to explain these trends from an evolutionary perspective. They explain the observed trends by host genetic factors [43, 52], viral adaptation of HIV to its host population [5, 8, 14, 38] and evolution of HIV SPVL to maintain higher viral fitness [12, 18, 39, 42]. HIV SPVL has been reported to be partially heritable from the infected donor to the susceptible recipient during HIV transmission [2, 12]. This has led to the hypothesis that sexual risk behaviour could influence the evolution of SPVL [11, 17]. A few studies have investigated sexual risk behaviour as a potential factor influencing SPVL [1, 17]. A recent study of Goodreau et al. [17] reports a positive relationship between relational concurrency and mean SPVL.

Mathematical modeling has been shown to be a valuable tool to investigate the evolution of SPVL and HIV virulence. Bezemer et al. [1] used a deterministic compartmental ordinary differential equation (ODE) model to study the influence of sexual risk behaviour on the HIV epidemic. Herbeck et al. [23] used an individual-based model to investigate trends in HIV virulence and community viral load. Their model contained functions for transmission and progression of HIV, SPVL, VL at multiple stages of HIV and SPVL heritability. Sexual relationships were described by a simplified contact network, without including sexual mixing patterns. Roberts et al. [33] designed a deterministic model to assess the impact of ART on viral evolution. Smith et al. [40] developed a compartmental model to study the influence of ART and pre-exposure prophylaxis (PrEP) on the evolution of HIV virulence. Goodreau et al. [17] used an individual- and network-based model to study the relationship between relational concurrency, HIV stages and evolution of HIV SPVL. However, none of these models implemented the effect of a sexually transmitted infection (STI) co-factor whereas the presence of an STI co-factor has been suggested to influence the trend in SPVL [19].

HIV-negative people infected with an STI are reported to have a higher risk for HIV acquisition, while HIV-positive people co-infected with an STI have a higher risk of transmitting HIV to another individual [47]. As a consequence, a higher incidence of an STI would lead to an increased number of new HIV patients. Furthermore, an observational study of Scheer et al [37] in MSM reported that people on ART show an increased risk of developing an STI, and this effect occurs most likely as a consequence of sexual risk behaviour. Moreover, having HIV also increases the risk for acquisition of an STI [26]. Furthermore, Gras et al. [19] observed a positive correlation between SPVL levels and the number of HIV patients. This has led to the hypothesis that the presence of an STI co-factor could influence the evolution of SPVL by increasing the number of new HIV cases.

In this study, we adapted an event-based, individual-based model to investigate the effect of sexual risk behaviour, STI co-infection and their combined effect on the evolution of HIV SPVL in MSM after the introduction of ART.

2 Materials and Methods

2.1 Simpact Cyan 1.0 modeling framework

Simpact Cyan 1.0 is a open-source framework for constructing individual-based models for simulating the transmission, diagnosis and treatment of HIV [29]. The program models each individual in a heterosexual population or a population of MSM, and the sexual relationships between individuals. The formation and dissolution of relationships, as well as birth, mortality, HIV transmission, diagnosis and treatment are represented by events, which have a certain risk of taking place at a certain moment, represented by their hazard function. Models are implemented in continuous time and updated each time an event happens. Simpact Cyan 1.0 also implements a generic sexually transmitted infection (STI) co-factor effect on HIV. Furthermore, it is also possible to simulate interventions by changing certain parameters during the simulation.

More detailed information on the Simpact Cyan 1.0 modeling framework is available in [29] and from http://www.simpact.org/.

2.2 Implementation of viral load and set-point viral load

When an individual becomes infected by HIV, a value for the SPVL is assigned to this individual. SPVL values on a base 10 logaritmic scale were drawn from a binomial distribution using the same mean and standard deviation in both dimensions: Embedded Image The two dimensions represent the persons that are already infected (x) and the persons becoming infected (y). The newly infected individual’s SPVL depends on the SPVL of the infector.

µx, µy: mean of the distribution of SPVL

σx, σy: standard deviation of the distribution of SPVL

ρ: correlation between x and y in the calculation of SPVL

Parameter ρ represents SPVL heritability and was taken from Blanquart et al. [2]. The mean and standard deviation were fitted to the histogram for MSM in Fraser et al. [11]. The values of the parameters are provided in Table S1.

The SPVL is the VL of a person during the chronic stage. In the acute stage, AIDS stage and final AIDS stage, the VL differs from the SPVL such that the HIV transmission hazard increases by a factor 10, 7 and 12 respectively. These multipliers are the default parameters in Simpact Cyan 1.0. They are based on a previous model fit [29] to data generated by Fraser et al. [11].

2.3 Implementation of treatment with anti-retroviral therapy

After being infected with HIV, an individual gets diagnosed as being HIV positive through a diagnosis event. In this study the hazard for diagnosis depends on a single baseline parameter: hazard = exp(baseline).

After being diagnosed, the person’s progression of HIV is monitored by inspecting his CD4 count during an HIV monitoring event. A person is eligible to be treated with ART if his CD4 count is below a given threshold.

In this study, we assume ART is offered to all eligible individuals and that all persons will accept treatment when offered. This is implemented by setting the ART acceptance threshold equal to 1 (see Table S1). Furthermore, we assume that no ART treatment is available before 1994, in correspondence with Palella et al. [32]. To implement that nobody will be treated before 1994, we set the hazard for diagnosis equal to 0 (exp(baseline) = 0). From 1994 to the end of the simulation, availability of ART in increased by

  • increasing the CD4 threshold from 200 in 1994 to 500 in 2014 with steps of 30 in every two years (see Table S2);

  • increasing the baseline parameter for the diagnosis hazard from −1.5 in 1994 to 1 in 2014 with steps of 0.25 in every two years (see Table S2).

When an individual starts with ART, his current VL decreases by 70% (default in Simpact Cyan 1.0).

2.4 Implementation of the STI co-factor effect

The start of the STI, in this study HSV-2, is implemented as an HSV-2 seeding event, where a certain amount of MSM are marked as HSV-2 infected. Transmission of the STI is described by the HSV-2 transmission hazard: Embedded Image In equation (2), the parameter a provides a baseline value for HSV-2 transmission. The parameter d describes the effect of HIV on HSV-2 transmission. Hi is a binary indicator for the HIV status of the HSV-2 infected person: 1 = HIV-positive, 0 = HIV-negative.

The parameters for HSV-2 seeding and transmission were taken from the literature and are provided in Table S1 of the Supplementary Material.

The effects of the STI on HIV are modeled within the HIV transmission hazard: Embedded Image In equation (3), the parameter a provides a baseline value for HIV transmission. Parameters b and c describe the influence of the viral load (V) of the HIV infected person on HIV transmission. Parameters e1 and e2 describe the effect of the STI on HIV transmission. Parameter e1 represents the higher risk of an HIV-positive person co-infected with HSV-2 to transmit HIV to another person. Parameter e2 represents the higher risk of a HIV-negative person infected with HSV-2 to acquire HIV. HSV 2HIV.infected and HSV 2HIV.uninfected are binary indicators for the HSV-2 status of the HIV-infected and HIV-uninfected person respectively (1 = HSV-2-positive, 0 = HSV-2-negative).

The values of the HIV transmission parameters a, b and c are determined during the calibration procedure described in paragraph 2.8. The formula for a hazard depending on viral load is based on a study by Hargrove et al. [21]. Parameters e1 and e2 were taken from the literature and are provided in Table S1 of the Supplementary Material.

2.5 Implementation of sexual behaviour

In Simpact Cyan 1.0, sexual behaviour among MSM is described by MSM formation and dissolution events, which are described by their hazard function. According to previous studies, relationship formation is influenced by

  • the number of relationships each partner has [17] and its assortativity [35];

  • the age of the partners [30, 34] and its assortativity [49].

In this study, the hazard for relationship formation is described by the following formula, where the selection of the parameters is based on the previous studies described above: Embedded Image In equation (4), F is a normalization factor, which roughly divides the formation hazard by the population size. Applying this normalization factor avoids that an increase in population size automatically results in more relationships. Pman1 and Pman2 represent the number of partners the two men already have; ageman1 and ageman2 represent the age of the two men. Parameter α0 represents a baseline value for the formation hazard. α12 describes the influence of the number of relationships the men already have, while α3 captures its assortativity by adding a weight for the effect of the difference in number of partners, between the two men on relationship formation. α12 ranges from minus infinity (-inf) to 0, where -inf corresponds to a population where everyone is monogamous, while 0 corresponds with no influence of the number of partners. The more negative, the less concurrent partnerships are formed. E.g. if this parameter is equal to −5, this means that the hazard is multiplied with exp(−5) (or in other words divided by exp(5)) if one of the men already has a partner. The parameter α3 is 0 in case of no assortativity. A more negative parameter corresponds to more assortativity. E.g. if this parameter is equal to log(0.5), then the formation hazard decreases with 50% if the number of current partners between the two men differs 1. α4 describes the influence of the mean age of the two men, while α5 captures its assortativity by adding a weight for the importance of the difference in age between the two men when forming a relationship. A negative value of α4 means that older persons form less relationships, while a positive value means that more relationships are formed by older persons. Both have been reported in the literature, e.g. in a study of Marcus et al. [30] older persons have less partners, while in another study conducted by Rozenberg et al. [34] less partners were reported for younger people. Therefore, both positive and negative values for parameter α4 were explored in this study. In case of no assortativity, α5 is 0. A more negative parameter corresponds with more assortativity. E.g. α5 = log(0.5) corresponds with a decrease of the formation hazard with 50% if the age between the two men differs 1 year.

The values of the parameters α0, α12, α3, α4, α5 are determined during the calibration procedure described in paragraph 2.8.

In this study, the break-up of relationships is described by the dissolution hazard: Embedded Image The value of the parameter α0 in equation (4) is determined during the calibration procedure described in paragraph 2.8.

2.6 Implementation of a change in risk behaviour with increasing availability of ART

In absence of a change in risk behaviour, the parameter α12 in de formation hazard (equation (4)) is constant over the whole simulation period.

A change in risk behaviour with increasing availability of ART is implemented by increasing the value of α12 with 0.05 every two years between 1994 and 2014, like represented in Table S2 of the Supplementary Material. The initial value of α12 is determined during the calibration procedure described in paragraph 2.8.

2.7 Model scenarios

In this study, we focus on the influence of sexual risk behaviour and STI co-infection dynamics on the evolution of the HIV set point viral load (SPVL) in MSM. Simulations are conducted for the period from 1980 to 2015, including both the period before and after the introduction of ART in 1994. After the introduction of ART, it is also assumed that the CD4 count threshold below which ART is offered increased [9], so that more people are treated. Furthermore, we assume that availability of ART has led to intensified HIV testing [15].

Four model scenarios including ART coverage are considered:

  • no STI co-factor effect and no change in risk behaviour (nSTI-nBC): null scenario;

  • STI co-factor effect (implemented as described in paragraph 2.4.) and no change in risk behaviour (STI-nBC);

  • a change in risk behaviour implemented as described in paragraph 2.5.) and no STI co-factor effect (nSTI-BC);

  • a change in risk behaviour (implemented as described in paragraph 2.5.) and an STI co-factor effect (implemented as described in paragraph 2.4.)(STI-BC).

In case any of the other scenarios has an improved GOF to the SPVL data compared to the null model, and/or better describes the SPVL trends, this suggests that the factors included in this scenario (STI co-factor and/or behavioural change) contribute to explaining the SPVL trends.

Parameters based on prior knowledge used in the model simulations, and intervention events for simulating the increase in ART coverage since 1994, and the increased risk behaviour with ART coverage (if applicable) are described in the Supplementary Material.

2.8 Model calibration

Table 1 presents the available literature data on the HIV epidemic in MSM in the Netherlands, used for model calibration. SPVL and HIV prevalence data are from different, but overlapping time intervals (SPVL: 1985-2007, HIV prevalence: 2003-2012).

View this table:
Table 1.

Available literature data used for model calibration. All data are for the Netherlands. CI: confidence interval; CrI: credible interval. * The values of the mean log10 SPVL in Gras et al. [19] were obtained from data of the mean HIV-1 RNA concentration at set-point of 612 MSM with date of seroconversion between 1985 and 2007. Gras et al. fitted a cubic spline through these data and estimated values for 1985, 1995 and 2007 from this spline.

For the nine parameters in Table 2, there were no literature values available. Therefore, these parameters were fitted to the data in Table 1 by applying an active learning approach [50] using the following steps:

View this table:
Table 2.

Parameters fitted to the data in Table 1. Initial ranges are based on the examples in the RSimpactHelper package https://github.com/wdelva/RSimpactHelp. The parameter names in the second column are directly related to Simpact Cyan 1.0.

  1. Select 10,000 parameter sets by applying Latin Hypercube Sampling (LHS) [41] and using the initial parameter ranges from Table 2.

  2. For each of the 10,000 parameter sets, run a simulation with Simpact Cyan 1.0.

  3. For each simulation, calculate the goodness-of-fit (GOF) based on the sum of squared relative errors [4].

  4. Apply the selection procedure of Castro Sanchez et al [3] based on the GOF measure to narrow the solution space (intervals for the parameters).

  5. Repeat steps 1-4 using the new intervals for the parameters until the GOF does not improve anymore.

When selecting the final parameters, it was also taken into account that the results show realistic point prevalence of partnership concurrency. Unfortunately, no estimates of partnership concurrency were reported for MSM in the Netherlands. Studies in the US report estimates of partnership concurrency equal to 14.8% [16] and 26.3% [48] for MSM. Based on these two studies, we assume that realistic values for the point prevalence of concurrency are between 10% and 30%.

Apart from the overall GOF, also the GOF based on SPVL data only, the GOF based on HIV prevalence data only, and the GOF to the ART coverage of 2013 are reported.

The HIV transmission hazard: hazard = exp(a + bV −c+ additional terms), where V is the current viral load, is based on equation (9) of Hargrove et al. [21]. Its default parameters (a = —1.3997, b = − 12.0220, c = 0.1649) are obtained by fitting parameters a, b and c to the data from Fraser et al. [11]. In this study, additional parameters are added to this hazard in case an STI co-factor is simulated (see paragraph 2.4).

3 Results

3.1 Model calibration

Table 3 shows the fitted parameters for the four model scenarios.

View this table:
Table 3.

Fitted parameters for the four model scenarios. nSTI-nBC: no STI co-factor, no behavioural change; STI-nBC: STI co-factor, no behavioural change; nSTI-BC: no STI co-factor, behavioural change; STI-BC: STI co-factor, behavioural change.

A detailed statistical analysis of the parameter space exploring

  • which parameters are most influenced by the data;

  • which regions contain a large amount of solutions with low relative sum of squared errors;

  • which parameters show associations

is provided in the Supplementary Material.

3.2 Goodness-of-fit

Four goodness-of-fit (GOF) statistics are considered: overall sum of squared relative errors (Figure 1A), sum of squared relative errors to the SPVL data only (Figure 1B), sum of squared relative errors to the HIV prevalence data only (Figure 1C) and squared relative error to the ART coverage in 2013 (Figure 1D). Figure 1A shows that including an STI co-factor effect or/and a behavioural change results in a better overall goodness-of-fit (GOF) to the literature data. However, when considering only GOF to the three data points on SPVL, none of the models showed a better fit to the data than the model not including any of these two effects (see Figure 1B). The improved overall GOF for the models including an STI co-factor (STI-nBC) and/or a behavioural change (nSTI-BC) can be explained by the improved GOF to the HIV prevalence data (see Figure 1C). The model including both an STI co-factor and a behavioural change shows the best GOF to the HIV prevalence data, and also the best overall GOF.

Figure 1.
Figure 1.

Boxplots for top 100 solutions based on overall goodness-of-fit. A) boxplots for overall sum of squared relative errors; B) boxplots for sum of squared relative errors to the SPVL data only; C) boxplots for sum of squared relative errors to the HIV prevalence data only; D) squared relative error to the ART coverage in 2013.

3.3 Trends in set-point viral load, HIV prevalence, ART coverage and HIV incidence

Figures 2 and 3A show the median population mean log10 SPVL and the interquartile range for 100 simulations with the four scenarios. The population mean log10 SPVL was calculated as the mean set-point HIV-RNA concentration at 9-27 months after seroconversion. None of the scenarios generates the increasing trend in population mean HIV SPVL between 1995 and 2007 observed by Gras et al [19]. In the nSTI-nBC (grey in Figure 3A) and nSTI-BC (blue in Figure 3A) scenarios (scenarios without STI co-factor), the mean log10 SPVL decreases from 1995 onwards. In the scenarios including an STI (STI-nBC (red in Figure 3A) and STI-BC (green in Figure 3A)), the mean log10 SPVL is constant between 1995 and 2000, and decreases from 2000 onwards.

Figure 2.
Figure 2.

Mean log10 HIV SPVL (mean set-point HIV-RNA concentration at 9-27 months of untreated HIV infection). Median of 100 simulations with the fitted parameters (solid line) and interquartile range (shaded area) for the period 1980-2015. A: no STI co-factor, no behavioural change (nSTI-nBC); B: STI co-factor, no behavioural change (STI-nBC); C: no STI co-factor, behavioural change (nSTI-BC); D: STI co-factor, behavioural change (STI-BC). See Supplementary Material for more detail on how the figures were generated. The black dots and error bars represent the mean log10 SPVL with 95% confidence interval for 1985, 1995 and 2007 from [19].

Figure 3.
Figure 3.

Mean log10 HIV SPVL (A), prevalence of HIV (B), ART coverage (C) and incidence of HIV (D) - median of 100 simulations with the fitted parameters (solid line) and interquartile range (shaded area) for the period 1980-2015. Grey: no STI co-factor, no behavioural change (nSTI-nBC); red: STI co-factor, no behavioural change (STI-nBC); blue: no STI co-factor, behavioural change (nSTI-BC); darkgreen: STI-cofactor, behavioural change (STI-BC). The black dots and the error bars represent the literature data from Table 1.

Figures 3B and C show how well the simulations of the four scenarios fit to the HIV prevalence and ART coverage data. All scenarios capture the increasing trend in HIV prevalence between 2003 and 2012. However, only the interquartile range for the STI-BC scenario includes the majority of the HIV prevalence data points from the period 2003-2007. In Figure 3C, we observe that the model scenario STI-nBC results in the closest estimate of the ART coverage of 62% in 2013 reported in [20]. Furthermore, for the STI-nBC scenario, also the year where the HIV incidence reaches its peak is the closest to the year of the peak observed in the literature data (2010 for STI-nBC, see Figure 3D; 2008 in the literature [45]). For all other scenarios, the HIV incidence increases until the end of the simulation (2015).

3.4 Relationship between HIV SPVL and point prevalence of concurrency

To investigate the relationship between HIV SPVL and points prevalence of concurrency, the weight for the number of partners (α12 in formula (4)) was varied between −10 and −1, while all other literature and estimated parameters were kept constant at their values in Table S1 and Table 3 respectively. A hundred simulations with values for α12 drawn from a uniform distribution with bounds −10 and −1 were performed. Figure 4 shows the mean log10 HIV SPVL of all individuals with a date of HIV seroconversion between 1980 and 2015 against mean point prevalence of concurrency over the period 1980-2015 for these 100 simulations with the four scenarios. All scenarios confirm the increasing trend of SPVL with relational concurrency described by Goodreau et al [17]. However, this increase in mean log 10 SPVL is small compared to the reported difference in mean log 10 SPVL between 1995 and 2007 (0.67), even when the point prevalence of concurrency changes from 0 to 1 (100%)(change of 0.10-0.15 of mean log 10 SPVL).

Figure 4.
Figure 4.

Mean log10 HIV SPVL of all individuals with a date of HIV seroconversion between 1980 and 2015 against mean point prevalence of concurrency over the period 1980-2015 for 100 simulations obtained by varying the weight for the number of partners between −10 and −1 while keeping all other literature and estimated parameters constant at their values in Table S1 and Table 3 respectively (black dots) and the trendline (blue line) with 95% confidence intervals (shaded area). A) no STI co-factor, no behavioural change; B) STI co-factor, no behavioural change; C) no STI co-factor, behavioural change; D) STI co-factor and behavioural change.

4 Discussion

After the introduction of ART, increasing trends in HIV SPVL and HIV incidence were observed in MSM [19, 45], which are counter-intuitive and not fully understood. In this study, we used an event- and individual-based model to investigate the impact of STI co-infection, changes in sexual risk behaviour and a combination of both on the evolution of HIV SPVL in MSM.

The results suggest that sexual risk behaviour influences population HIV SPVL (see Figure 4). Furthermore, the model including only an STI co-factor and no behavioural change is the only model that simulates the drop in HIV incidence after an initial increase reported by Van Sighem et al. [45](see Figure 3D). However, the HIV prevalence curves are closer to the observed values in [7] when also including behavioural changes. For the ART coverage, the scenario including only an STI co-factor and no behavioural change results in the closest approximation to the reported value in [20].

In this study, increased access to ART was modeled by increasing the CD4 threshold and assuming earlier diagnosis because of increased HIV testing in the presence of ART. In 2011, the WHO recommended treatment in all serodiscordant couples [51]. However, according to the WHO, many people do not disclose their HIV status with their sexual partner(s) (https://www.who.int/hiv/pub/guidelines/9789241501972/en/). The WHO also reports that to be treated with ART when having CD4 counts higher than the CD4 count initiation threshold, it’s necessary that both partners are tested and are willing to disclose their HIV status to each other. Furthermore the estimated proportion of undiagnosed HIV cases in MSM in 2015 in the Netherlands was about 10% [44]. Therefore, it would not be a realistic scenario to assume that for all serodiscordant couples, the HIV positive partner is on ART. In the model scenarios in this study, about 60% of the serodiscordant couples are on ART. To the best of our knowledge, no estimates are available on the proportion of HIV-positive MSM with CD4 counts higher than the CD4 initiation threshold that are willing to undergo a test and share their HIV status with their partner(s). This hampers making an accurate estimation of the percentage of serodiscordant couples on ART.

Although none of the three scenarios including an STI co-factor and/or behavioural changes could explain the increasing trend in population mean SPVL reported in Gras et al. [19], all scenarios can explain the increase in new HIV cases between 1996 and 2008 reported by Van Sighem et al. [45](see Figure 3D). All scenarios could confirm the higher mean population SPVL with higher relational concurrency reported by Goodreau et al [17](see Figure 4).

There are several potential reasons that none of the scenarios can explain the trends in SPVL reported in Gras et al. [19]. First, in addition to the factors explored in this study (treatment, sexual behaviour and STI co-factor effects), also viral evolution has been reported to have an effect on the population mean SPVL. However, to the best of our knowledge, no modeling framework is available that can incorporate treatment, sexual risk behaviour, STI co-infection and viral evolution into a single model. This will be explored in future research. Second, the available data for this case study as reported in the literature is scarce. Making more data available will lead to better estimation of the model parameters, and possibly to improved models that can better explain the trends in SPVL observed in the literature. Furthermore, the data on SPVL from [19] (see Table 1) for 1995 are based on a considerably lower amount of persons than for 1985 and 2007, which creates uncertainty about the observed trends. According to Cumming and Finch [6], methods for significance testing based on confidence intervals (CI) are reliable when both sample sizes are equal or larger than 10, and when the margins of error differ not more than a factor 2. For 1995 the figure with the data of the mean HIV-1 RNA concentration at set-point in Gras et al. [19] displays less than 10 data points. This would mean that although a proportion overlap of 0.387 between the CIs for 1985 and 1995 is obtained, pointing to a significant difference between the two samples according to Cumming and Finch [6], there is uncertainty about the observed decrease between these two time points. However, the values of the mean log10 SPVL for 1985, 1995 and 2007 in Gras et al. [19] were obtained by fitting a cubic spline through the data of the mean HIV-1 RNA concentration at set-point of 612 MSM with date of seroconversion between 1985 and 2007. Third, only a generic STI co-factor effect was included in the model. Including a more specific STI co-factor effect, e.g. a herpes simplex virus (HSV-2) co-factor effect where all stages of HSV-2 are described in detail may lead to improved models for this case study. However, such model has many more parameters than the model used in this study, and as a consequence requires more data than currently available for parameter estimation.

Weiss et al. [49] implemented a model that fully parametrizes the effect of an STI on HIV transmission to study chlamydia and gonorrhea screening and treatment. However, the model doesn’t implement the effect of HIV on STI transmission, while Simpact Cyan 1.0 does. Furthermore, the model of Weiss et al. is implemented in discrete time, while Simpact Cyan 1.0 is implemented in continuous time. Moreover, Simpact Cyan 1.0 also implements CD4 counts, an option that is not available in the model of Weiss et al.

The results of the model simulations in this study identify factors that should be monitored during future research on the evolution of population mean SPVL. These factors include viral evolution, partnership disclosure and natural history of STI’s. In this way, the models can contribute to the study design of future studies. An overview of previous applications of mathematical models to inform study design in epidemiology is provided by Herzog et al. [25].

To implement a model that fully parametrizes the natural history of an STI, and the effect of this STI on HIV transmission and vice versa, we recommend to collect the following data in future studies:

  • monitoring of a cohort of HIV mono-infected, STI mono-infected and HIV-STI co-infected individuals for the study period of interest; collect prevalence and incidence data at least on a yearly basis;

  • collecting characteristics for the cohort such as estimated time of HIV infection, estimated time of STI infection and period of HIV ART treatment if applicable.

In summary, the results of this study suggest that sexual risk behaviour could influence SPVL and indicate that more data are needed to test the effect of STI co-factors on SPVL. Furthermore, the trends described in the literature cannot be explained by sexual risk behaviour and STI co-factors only. Future research to understand SPVL evolution should also consider models that include factors related to viral evolution and describe STI co-factors in more detail. To accomplish these goals, more data has to be collected and made available.

Data Availability

This study concerns a study on mathematical modeling of HIV/AIDS, and uses only information from the literature. It does not include any newly generated data sets or data extracted from repositories.

Supporting Information

See supplementary.pdf

Acknowledgments

The research conducted by DMH and NH in this study was funded by the Fonds Wetenschappelijk Onderzoek - Vlaanderen (Research Foundation – Flanders; FWO, http://www.fwo.be/en/) (Grant agreements G0E8416N and G0B2317N).

WD was supported by a postdoctoral fellowship from FWO (12L5816N).

The computational resources and services used in this work were provided by the VSC (Flemish Supercomputer Center), funded by the Research Foundation - Flanders (FWO) and the Flemish Government – department EWI.

We thank Cécile Kremer (I-Biostat,Data Science Institute, UHasselt) for her support on the literature study to obtain the parameter settings described in Table S1 of the Supplementary Material.

Footnotes

  • ↵* diana.hendrickx{at}uhasselt.be

References

  1. 1.↵
    D. Bezemer, F. de Wolf, M. C. Boerlijst, A. van Sighem, T. D. Hollingsworth, and C. Fraser. 27 years of the HIV epidemic amongst men having sex with men in the Netherlands: an in depth mathematical model-based analysis. Epidemics, 2(2):66–79, 2010.
    OpenUrlCrossRefPubMed
  2. 2.↵
    F. Blanquart, C. Wymant, M. Cornelissen, A. Gall, M. Bakker, D. Bezemer, M. Hall, M. Hillebregt, S. H. Ong, J. Albert, et al. Viral genetic variation accounts for a third of variability in HIV-1 set-point viral load in europe. PLoS biology, 15(6):e2001855, 2017.
    OpenUrlCrossRef
  3. 3.↵
    A. Y. Castro Sanchez, M. Aerts, Z. Shkedy, P. Vickerman, F. Faggiano, G. Salamina, and N. Hens. A mathematical model for HIV and hepatitis C co-infection and its assessment from a statistical perspective. Epidemics, 5(1):56–66, 2013.
    OpenUrlCrossRefPubMedWeb of Science
  4. 4.↵
    K. Chen, S. Guo, Y. Lin, and Z. Ying. Least absolute relative error estimation. Journal of the American Statistical Association, 105(491):1104–1112, 2010.
    OpenUrlCrossRefPubMed
  5. 5.↵
    L. A. Cotton, X. T. Kuang, A. Q. Le, J. M. Carlson, B. Chan, D. R. Chopera, C. J. Brumme, T. J. Markle, E. Martin, A. Shahid, et al. Genotypic and functional impact of HIV-1 adaptation to its host population during the North American epidemic. PLoS genetics, 10(4):e1004295, 2014.
    OpenUrl
  6. 6.↵
    G. Cumming and S. Finch. Inference by eye: confidence intervals and how to read pictures of data. American Psychologist, 60(2):170, 2005.
    OpenUrlCrossRefPubMedWeb of Science
  7. 7.↵
    E. L. O. de Coul, I. Schreuder, S. Conti, A. van Sighem, M. Xiridou, M. G. Van Veen, and J. C. Heijne. Changing patterns of undiagnosed HIV infection in the Netherlands: who benefits most from intensified HIV test and treat policies? PLoS One, 10(7):e0133232, 2015.
    OpenUrlPubMed
  8. 8.↵
    D. Edo-Matas, A. Rachinger, L. C. Setiawan, B. D. Boeser-Nunnink, A. B. D. van’t Wout, P. Lemey, and H. Schuitemaker. The evolution of human immunodeficiency virus type-1 (HIV-1) envelope molecular properties and coreceptor use at all stages of infection in an HIV-1 donor-recipient pair. Virology, 422(1):70–80, 2012.
    OpenUrlCrossRefPubMed
  9. 9.↵
    S. P. Eholié, A. Badje, G. M. Kouame, J.-B. N’takpe, R. Moh, C. Danel, and X. Anglaret. Antiretroviral treatment regardless of CD4 count: the universal answer to a contextual question. AIDS research and therapy, 13(1):27, 2016.
    OpenUrl
  10. 10.
    EuroHIV. HIV/AIDS surveillance in Europe. Mid-year report 2005. saint-maurice, institut de veille sanitaire, 2006.
  11. 11.↵
    C. Fraser, T. D. Hollingsworth, R. Chapman, F. de Wolf, and W. P. Hanage. Variation in HIV-1 set-point viral load: epidemiological analysis and an evolutionary hypothesis. Proceedings of the National Academy of Sciences, 104(44):17441–17446, 2007.
    OpenUrlAbstract/FREE Full Text
  12. 12.↵
    C. Fraser, K. Lythgoe, G. E. Leventhal, G. Shirreff, T. D. Hollingsworth, S. Alizon, and S. Bonhoeffer. Virulence and pathogenesis of HIV-1 infection: an evolutionary perspective. Science, 343(6177):1243727, 2014.
    OpenUrlAbstract/FREE Full Text
  13. 13.↵
    R. B. Geskus, M. Prins, J.-B. Hubert, F. Miedema, B. Berkhout, C. Rouzioux, J.-F. Delfraissy, and L. Meyer. The HIV RNA setpoint theory revisited. Retrovirology, 4(1):65, 2007.
    OpenUrlCrossRefPubMed
  14. 14.↵
    E. F. Gijsbers, H. Schuitemaker, and N. A. Kootstra. HIV-1 transmission and viral adaptation to the host. Future Virology, 7(1):63–71, 2012.
    OpenUrl
  15. 15.↵
    Global HIV Prevention Working Group. HIV prevention in the era of expanded treatment access. Technical report, Washington, DC: Kaiser Family Foundation, 2004.
  16. 16.↵
    S. M. Goodreau, E. S. Rosenberg, S. M. Jenness, N. Luisi, S. E. Stansfield, G. A. Millett, and P. S. Sullivan. Sources of racial disparities in hiv prevalence in men who have sex with men in atlanta, ga, usa: a modelling study. The Lancet HIV, 4(7):e311–e320, 2017.
    OpenUrl
  17. 17.↵
    S. M. Goodreau, S. E. Stansfield, J. T. Murphy, K. C. Peebles, G. S. Gottlieb, N. F. Abernethy, J. T. Herbeck, and J. E. Mittler. Relational concurrency, stages of infection, and the evolution of HIV set point viral load. Virus Evolution, 4(2):vey032, 2018.
    OpenUrl
  18. 18.↵
    L. Gras, R. B. Geskus, S. Jurriaans, M. Bakker, A. van Sighem, D. Bezemer, C. Fraser, J. M. Prins, B. Berkhout, F. de Wolf, et al. Has the rate of CD4 cell count decline before initiation of antiretroviral therapy changed over the course of the Dutch HIV epidemic among MSM? PLoS One, 8(5):e64437, 2013.
    OpenUrl
  19. 19.↵
    L. Gras, S. Jurriaans, M. Bakker, A. van Sighem, D. Bezemer, C. Fraser, J. Lange, J. M. Prins, B. Berkhout, F. de Wolf, et al. Viral load levels measured at set-point have risen over the last decade of the HIV epidemic in the Netherlands. PLoS One, 4(10):e7365, 2009.
    OpenUrlCrossRefPubMed
  20. 20.↵
    S. Gupta and R. Granich. National HIV care continua for key populations: 2010 to 2016. Journal of the International Association of Providers of AIDS Care (JIAPAC), 16(2):125–132, 2017.
    OpenUrl
  21. 21.↵
    J. Hargrove, H. Eastwood, G. Mahiane, and C. van Schalkwyk. How should we best estimate the mean recency duration for the BED method? PloS one, 7(11):e49661, 2012.
    OpenUrlPubMed
  22. 22.↵
    J. T. Herbeck, J. E. Mittler, G. S. Gottlieb, S. M. Goodreau, J. T. Murphy, A. Cori, M. Pickles, and C. Fraser. Evolution of HIV virulence in response to widespread scale up of antiretroviral therapy: a modeling study. Virus Evolution, 2(2):vew028, 2016.
    OpenUrlCrossRef
  23. 23.↵
    J. T. Herbeck, J. E. Mittler, G. S. Gottlieb, and J. I. Mullins. An HIV epidemic model based on viral load dynamics: value in assessing empirical trends in HIV virulence and community viral load. PLoS computational biology, 10(6):e1003673, 2014.
    OpenUrl
  24. 24.↵
    J. T. Herbeck, V. Müller, B. S. Maust, B. Ledergerber, C. Torti, S. Di Giambenedetto, L. Gras, H. F. Günthard, L. P. Jacobson, J. I. Mullins, et al. Is the virulence of HIV changing? A meta-analysis of trends in prognostic markers of HIV disease progression and transmission. AIDS (London, England), 26(2):193, 2012.
    OpenUrlCrossRefPubMedWeb of Science
  25. 25.↵
    S. A. Herzog, S. Blaizot, and N. Hens. Mathematical models used to inform study design or surveillance systems in infectious diseases: a systematic review. BMC infectious diseases, 17(1):775, 2017.
    OpenUrl
  26. 26.↵
    C. Hill, E. McKinney, C. Lowndes, H. Munro, G. Murphy, J. Parry, O. Gill, et al. Epidemiology of herpes simplex virus types 2 and 1 amongst men who have sex with men attending sexual health clinics in England and Wales: implications for HIV prevention and management. Eurosurveillance, 14(47):19418, 2009.
    OpenUrl
  27. 27.↵
    E. L. Korenromp, B. G. Williams, G. P. Schmid, and C. Dye. Clinical prognostic value of RNA viral load and CD4 cell counts during untreated HIV-1 infection—a quantitative review. PloS one, 4(6):e5950, 2009.
    OpenUrlCrossRefPubMed
  28. 28.↵
    T. S. Li, R. Tubiana, C. Katlama, V. Calvez, H. A. Mohand, and B. Autran. Long-lasting recovery in CD4 T-cell function and viral-load reduction after highly active antiretroviral therapy in advanced HIV-1 disease. The Lancet, 351(9117):1682–1686, 1998.
    OpenUrl
  29. 29.↵
    J. Liesenborgs, D. M. Hendrickx, E. Kuylen, D. Niyukuri, N. Hens, and W. Delva. SimpactCyan 1.0: An open-source simulator for individual-based models in HIV epidemiology with R and Python interfaces. Scientific Reports, 9(1):1–13, 2019.
    OpenUrl
  30. 30.↵
    U. Marcus, F. Hickson, P. Weatherburn, A. J. Schmidt, et al. Age biases in a large hiv and sexual behaviour-related internet survey among msm. BMC public health, 13(1):826, 2013.
    OpenUrl
  31. 31.↵
    D. K. Owens, V. Sundaram, L. C. Lazzeroni, L. R. Douglass, P. Tempio, M. Holodniy, G. D. Sanders, V. M. Shadle, V. C. McWhorter, T. Agoncillo, et al. HIV testing of at risk patients in a large integrated health care system. Journal of general internal medicine, 22(3):315–320, 2007.
    OpenUrlCrossRefPubMedWeb of Science
  32. 32.↵
    F. J. Palella Jr, K. M. Delaney, A. C. Moorman, M. O. Loveless, J. Fuhrer, G. A. Satten, D. J. Aschman, S. D. Holmberg, and H. O. S. Investigators. Declining morbidity and mortality among patients with advanced human immunodeficiency virus infection. New England Journal of Medicine, 338(13):853–860, 1998.
    OpenUrlCrossRefPubMedWeb of Science
  33. 33.↵
    H. E. Roberts, P. J. Goulder, and A. R. McLean. The impact of antiretroviral therapy on population-level virulence evolution of HIV-1. Journal of the Royal Society Interface, 12(113):20150888, 2015.
    OpenUrl
  34. 34.↵
    E. S. Rosenberg, P. S. Sullivan, E. A. DiNenno, L. F. Salazar, and T. H. Sanchez. Number of casual male sexual partners and associated factors among men who have sex with men: results from the national hiv behavioral surveillance system. BMC public health, 11(1):189, 2011.
    OpenUrlCrossRefPubMed
  35. 35.↵
    G. Rozhnova, M. F. S. van der Loeff, J. C. Heijne, and M. E. Kretzschmar. Impact of heterogeneity in sexual behavior on effectiveness in reducing hiv transmission with test-and-treat strategy. PLoS computational biology, 12(8):e1005012, 2016.
    OpenUrl
  36. 36.↵
    M. S. Saag, M. Holodniy, D. Kuritzkes, W. O’Brien, R. Coombs, M. Poscher, D. Jacobsen, G. Shaw, D. Richman, and P. Volberding. HIV viral load markers in clinical practice. Nature medicine, 2(6):625, 1996.
    OpenUrlCrossRefPubMedWeb of Science
  37. 37.↵
    S. Scheer, P. L. Chu, J. D. Klausner, M. H. Katz, and S. K. Schwarcz. Effect of highly active antiretroviral therapy on diagnoses of sexually transmitted diseases in people with aids. The Lancet, 357(9254):432–435, 2001.
  38. 38.↵
    I. M. Schellens, M. Navis, H. W. van Deutekom, B. Boeser-Nunnink, B. Berkhout, N. Kootstra, F. Miedema, C. Kesmir, H. Schuitemaker, D. van Baarle, et al. Loss of HIV-1-derived cytotoxic T lymphocyte epitopes restricted by protective HLA-B alleles during the HIV-1 epidemic. Aids, 25(14):1691–1700, 2011.
    OpenUrlCrossRefPubMedWeb of Science
  39. 39.↵
    G. Shirreff, L. Pellis, O. Laeyendecker, and C. Fraser. Transmission selects for HIV-1 strains of intermediate virulence: a modelling approach. PLoS computational biology, 7(10):e1002185, 2011.
    OpenUrl
  40. 40.↵
    D. R. Smith and N. Mideo. Modelling the evolution of HIV-1 virulence in response to imperfect therapy and prophylaxis. Evolutionary applications, 10(3):297–309, 2017.
    OpenUrl
  41. 41.↵
    M. Stein. Large sample properties of simulations using Latin hypercube sampling. Technometrics, 29(2):143–151, 1987.
    OpenUrlCrossRefWeb of Science
  42. 42.↵
    A. C. van der Kuyl, M. Bakker, S. Jurriaans, N. K. Back, A. O. Pasternak, M. Cornelissen, and B. Berkhout. Translational HIV-1 research: from routine diagnostics to new virology insights in Amsterdam, the Netherlands during 1983-2013. Retrovirology, 10(1):93, 2013.
    OpenUrl
  43. 43.↵
    D. van Manen, L. Gras, B. D. Boeser-Nunnink, A. I. van Sighem, I. Maurer, M. M. M. Ruiz, A. M. Harskamp, R. Steingrover, J. M. Prins, F. de Wolf, et al. Rising HIV-1 viral load set point at a population level coincides with a fading impact of host genetic factors on HIV-1 control. Aids, 25(18):2217–2226, 2011.
    OpenUrlCrossRefPubMed
  44. 44.↵
    A. Van Sighem, T. Boender, F. Wit, C. Smit, A. Matser, P. Reiss, et al. Monitoring report 2016. Human immunodeficiency virus (HIV) infection in the Netherlands. Amsterdam: Stichting HIV Monitoring, 2016.
  45. 45.↵
    A. Van Sighem, F. Wit, A. Boyd, C. Smit, A. Matser, and P. Reiss. Monitoring report 2019. Human Immunodeficiency Virus (HIV) infection in the Netherlands. Amsterdam: HIV monitoring foundation, 2019.
  46. 46.
    M. Van Veen, F. Koedijk, I. Van den Broek, E. Op de Coul, I. De Boer, A. Van Sighem, and M. Van der Sande. Sexually transmitted infections in the Netherlands in 2006. RIVM report, 210261003, 2007.
  47. 47.↵
    H. Ward and M. Rönn. The contribution of STIs to the sexual transmission of HIV. Current Opinion in HIV and AIDS, 5(4):305, 2010.
    OpenUrlCrossRef
  48. 48.↵
    K. M. Weiss, S. M. Goodreau, M. Morris, P. Prasad, R. Ramaraju, T. Sanchez, and S. M. Jenness. Egocentric sexual networks of men who have sex with men in the united states: Results from the artnet study. Epidemics, 30:100386, 2020.
    OpenUrlCrossRef
  49. 49.↵
    K. M. Weiss, J. S. Jones, E. J. Anderson, T. Gift, H. Chesson, K. Bernstein, K. Workowski, A. Tuite, E. S. Rosenberg, P. S. Sullivan, et al. Optimizing coverage vs frequency for sexually transmitted infection screening of men who have sex with men. In Open forum infectious diseases, volume 6, page ofz405. Oxford University Press US, 2019.
    OpenUrl
  50. 50.↵
    L. Willem, S. Stijven, E. Vladislavleva, J. Broeckhove, P. Beutels, and N. Hens. Active learning to understand infectious disease models and improve policy making. PLoS computational biology, 10(4):e1003563, 2014.
    OpenUrl
  51. 51.↵
    World Health Organization. Guidance on couples hiv testing and counselling including antiretroviral therapy for treatment and prevention in serodiscordant couples: recommendations for a public health approach. 2012.
  52. 52.↵
    L. Yue, H. A. Prentice, P. Farmer, W. Song, D. He, S. Lakhi, P. Goepfert, J. Gilmour, S. Allen, J. Tang, et al. Cumulative impact of host and viral factors on HIV-1 viral-load control during early infection. Journal of virology, 87(2):708–715, 2013.
    OpenUrlAbstract/FREE Full Text
Back to top
Posted March 18, 2021.
Download PDF

Supplementary Material

Data/Code
Email
Influence of sexual risk behaviour and STI co-infection dynamics on the evolution of HIV set point viral load in MSM
Diana M Hendrickx, Wim Delva, Niel Hens
medRxiv 19011221; doi: https://doi.org/10.1101/19011221
Twitter logo Facebook logo LinkedIn logo Mendeley logo
Citation Tools

Subject Area