Individuals living with HIV experience a much higher risk of progression from latent M. tuberculosis infection to active tuberculosis (TB) disease relative to individuals with intact immune systems. A several-month daily course of a single drug during latent infection (i.e. isoniazid preventive therapy (IPT)) has proved in clinical trials to substantially reduce an HIV-infected individual's risk of TB disease. As a result of these findings and ongoing studies, the World Health Organization has produced strong guidelines for implementing IPT on a community-wide scale for individuals with HIV at risk of TB disease. To date, there has been limited use of IPT at a community-wide level. In this paper, we present a new co-network model for HIV and TB co-epidemics to address questions about how the population-level impact of community-wide IPT may differ from the individual-level impact of IPT offered to selected individuals. In particular, we examine how the effect of clustering of contacts within high-TB incidence communities may affect the rates of re-infection with TB and how this clustering modifies the expected population-level effects of IPT. We find that populations with clustering of respiratory contacts experience aggregation of TB cases and high numbers of re-infection events. While, encouragingly, the overall population-level effects of community-wide IPT appear to be sustained regardless of network structure, we find that in populations where these contacts are highly clustered, there is dramatic heterogeneity in the impact of IPT: in some sub-regions of these populations, TB is nearly eliminated, while in others, repeated re-infection almost completely undermines the effect of IPT. Our findings imply that as IPT programmes are brought to scale, we should expect local heterogeneity of effectiveness as a result of the complex patterns of disease transmission within communities.
The emergence of HIV over the past three decades has complicated the global control of tuberculosis (TB) . In some sub-Saharan African countries with generalized (more than 1% of the population infected) HIV epidemics, TB incidence rates have increased over threefold in less than 20 years (for example, Lesotho and South Africa [2,3]).
Despite the widespread availability of effective anti-TB drug treatment, TB remains a leading cause of death for HIV-infected individuals. The latest estimates indicate that nearly a third of all HIV/AIDS-related deaths might be attributable to TB . This grim statistic has motivated the adoption of new approaches for reducing the burden of HIV-associated TB through the use of intensive case identification, infection control to reduce nosocomial transmission and isoniazid preventive therapy (IPT) (the three I's ).
IPT reduces the risk of progression from latent infection with M. tuberculosis to active TB disease through the use of a single drug given for six to 12 months ; some recent trials include treatment arms of up to 36 months of therapy . This beneficial effect has been demonstrated for both individuals with HIV  and without HIV/TB co-infection . Since the primary mechanism by which HIV epidemics result in a rising incidence of TB is by increasing the rate at which co-infected individuals progress from latent M. tuberculosis infection to active TB disease, community-wide IPT (i.e. IPT delivered to all eligible HIV co-infected individuals in a community) promises to reduce the incidence of HIV-associated TB. Several mathematical models have identified IPT interventions as a important tool for long-term reductions in rates of TB incidence [10–12].
However, because of the complex dynamics underlying transmission of TB and HIV, it is not straightforward to extrapolate the observed effect of IPT delivered to individuals to the combined effect of IPT delivered to communities. Effects of interventions at the community level may not mirror those observed at the individual level. For example, the overall impact of IPT may be underestimated if IPT reduces onward TB transmission; in this case, even individuals not eligible for or not given IPT treatment may benefit. In contrast, the overall impact may be overestimated if individuals remain at high risk of re-infection after their IPT treatment is concluded.
The overall outcome after IPT is delivered to an individual will depend not only on the effectiveness of IPT in preventing disease in that individual, but also on the continuing risk of M. tuberculosis infection that persists beyond the time that IPT is administered. Assuming that prevalent cases of TB are treated, the continuing risk to an individual after IPT is a function of two factors: (i) the prevalence of latent infection among their respiratory contacts and (ii) the risk that these latent infections will progress to TB and result in re-infection.
The prevalence of latent infection among an individual's respiratory contacts in turn depends on the prevalence of latent M. tuberculosis infection in the population, and the extent to which these infections are aggregated or clustered in a community. The risk that such infections will progress to disease depends on the prevalence and aggregation of HIV infections (since HIV/TB co-infection is the largest known risk factor for progression) and on the probability that latent M. tuberculosis among these HIV co-infected contacts are also given IPT treatment. Thus, the overall effectiveness of IPT depends on the structure of the networks underlying disease transmission (which drives the clustering of HIV and TB on these networks), the prevalence of TB and HIV in the communities and the coverage levels and patterns of distribution of the IPT intervention.
Here we use a model to examine patterns of HIV and TB dynamics to relate the individual-level impact of IPT to its potential impact at the community level. The model uses two distinct contact networks representing respiratory and sexual transmission. Studies recording different types of contact over time, for example [13,14], can inform us about possible respiratory networks; however, not only are these studies in settings different from those we wish to represent (Europe instead of Africa), it is also difficult to extract from the results the contact patterns that would be sustained and close enough to allow TB transmission. In general, it is assumed that respiratory networks will have more local ‘clustering’ of contacts than sexual networks , since respiratory contacts will be shared within family, school and social groups to a greater extent than sexual contacts.
There is some evidence that TB cases are clustered, for example, a study documented TB disease patterns in an region of Cape Town, South Africa (an area of high TB incidence), and using spatial analysis techniques found there was clustering of active TB cases. They linked these ‘hot-spots’ of disease to other factors like overcrowding and Shebeen (informal bar) locations . A more recent study of TB cases in The Gambia found that there was significant clustering of cases and suggested that using cluster detection techniques (such as contact tracing and case detection in the community) to identify clusters and their sources may have a positive impact on treatment programmes . An analysis of the social network of the initial cases in a TB outbreak in British Columbia revealed that all the patients frequented common locations in the community and had shared contacts, even if they had no direct contact , again indicating a clustering of TB transmission. It is only recently that studies such as these have revealed the impact that network structures have on TB transmission; to investigate this further, we use varying types of respiratory networks to explore the effect on the epidemic.
We use a dynamic model of HIV and TB transmission over these coupled networks and investigate the effects of local network structure on the effectiveness of IPT interventions at the population and the individual level. At the population level, we examine the impact of IPT on TB incidence and prevalence over time, while at the individual level, we compare cohorts of individuals receiving no IPT with those receiving varying-length IPT courses. To further elucidate the impact of network structure, we examine patterns of spatial clustering in TB infection events and explore whether spatial heterogeneity can explain variability in the impact of IPT.
2.1. Network generation
The respiratory contact networks were generated using the method described by Read & Keeling , with an average degree of 12 and a low variance in the degree distribution. In this algorithm, every individual is assigned a location and connected to other nodes with a probability proportional to their distance apart. On more local networks, nodes form connections preferentially with those closer to them, thereby forming a more clustered structure, while on global graphs, connections are, on average, longer and there are fewer cliques. The ‘location’ represents both geographical and social location: we wish to represent that some people are at higher risk than others for exposure to TB, both because of their shared environment and because of whom they are in contact with. This network generation method allows us to investigate how changing these risks affects the epidemic: we investigated dynamics on local, intermediate and more global networks.
In contrast, in the sexual contact network, most individuals have a small degree and there are a few highly connected individuals. The network was generated following the configuration model method  using a power law degree distribution . The network is static, so the degree represents the total number of sexual contacts over the duration of the simulation. The mean degree of an individual is around 1.8, with 5 per cent of individuals having no sexual contacts. Examples of the networks are shown in figure 1.
2.2. Disease transmission
From the respiratory and sexual contact networks, we create the co-network by combining them in a multi-graph, so every individual may have contacts in both networks. Each individual is in one of seven states: susceptible to both diseases, infected with latent TB (exposed), infected with active TB, HIV-positive, infected with latent TB and HIV-positive, infected with active TB and HIV-positive (co-infected), dead. Each state (except if dead) has corresponding parameters reflecting mortality, TB progression and recovery from active TB (table 1). At each (monthly) time step, the status change for each individual is determined stochastically, based on the probability of infection with HIV and/or TB by infectious contacts on the networks and the relevant mortality, progression and recovery rates. The possible transitions between states in the network model are illustrated in figure 2 and described in detail in appendix A.
While the natural history of TB is modelled in some detail, we include only a single state of HIV infection for simplicity, recognizing that in reality, the state of HIV infection is heterogeneous .
2.3. IPT cohort trials
We generated populations of 10 000 individuals with either more or less clustered respiratory networks. We scaled the TB and HIV transmission parameters to allow for stable and comparable epidemics on each network structure. We then simulated IPT trials by running three simulations, ‘enrolling’ a cohort of 200 HIV-infected, IPT-eligible individuals in each: the control cohort received no IPT (and served as our reference group), the first intervention cohort received six months of IPT and the second intervention cohort received 36 months of IPT. We then followed these cohorts for 10 years to examine the simulated effects of these interventions. Each simulation was repeated 10 times on each network structure. Because we knew precise timing of infection and disease events in these cohorts, we could construct Kaplan–Meier curves for displaying event-free survival. To examine whether there was spatial variability in the effects of IPT, for each simulation, we split the population of individuals into four regions and plotted separate Kaplan–Meier curves for each sub-area.
2.4. Spatial analysis
Heat maps and spatial statistical analysis were used to find significant spatial clusters in the number of TB re-infection events per individual. We used SatScan1 to identify statistically significant circular clusters of individuals re-infected three or more times. Finally, we calculated Ripley's K-statistic and its derivative, the L-statistic. The K-statistic identifies clustering in point patterns  and the L-statistic is used to identify deviations from complete spatial randomness; these were examined only for individuals re-infected three or more times. We then used a Monte Carlo approach to determine the significance of the results.
The model simulations recreate realistic HIV and TB trends from typical settings in sub-Saharan Africa. The initial conditions were chosen to approximately match those of high-TB settings such as Botswana in 1990 [2,3]. With no intervention, the prevalence of HIV and latent TB in the population grows steadily, slowing after 20 years. As is seen in regions of the world with high TB and HIV incidence, for all networks in our model, the co-infection rate is very high; up to 75 per cent of those with active TB are co-infected with HIV.
IPT has a dramatic impact on TB in the simulations, an effect that increases with IPT coverage (figure 3). These effects at the population level do not exhibit a strong dependence on the respiratory contact structure of the population. The effect of increased coverage is not linear; while there is a substantial difference in TB incidence between 10 and 25 per cent coverage with IPT, there is little difference between 25 and 50 per cent coverage. Each line in the figure represents an average over 50 simulations, and the percentage by which the individual runs varied was on average 6 per cent for locally and globally connected networks (results not shown).
Figure 4 illustrates the result of our individual-level (cohort-based) comparisons of the performance of IPT in the different networks. In figure 4a,b, the Kaplan–Meier survival plots show the per cent of the cohort that has not developed active TB since IPT began. While the endpoints at 120 months are similar on both network topographies, the difference between the IPT and non-IPT cohorts emerges more rapidly in the more global network. Indeed, in the more local networks, the overall disease dynamics are slower , but after 10 years, the performance of IPT is better than in the more global networks. Unsurprisingly, we find that a strategy that uses a 36 month course of IPT has a more pronounced effect than a six month course of IPT.
When we now examine the more local effects of IPT by splitting the population into spatial regions (figure 4c), we see considerably more variation in the effect of IPT between areas on the more local networks. On the local networks, a 36 month course of IPT had a more consistent effect in the four regions than the six month course. However, on more global networks, there was little variability across the regions for both the six and 36 month courses.
The spatial analysis reveals the origin of the high heterogeneity between regions in local networks. We find that the total number of re-infection events throughout the simulation is far lower for more global networks, consistent with the idea that TB infection is contained within tightly linked clusters of contacts and within these more infection events are likely to be re-infections. However, similar numbers of individuals are re-infected in each network. This means that there are more individuals suffering multiple re-infections on local networks; these individuals were re-infected more often (3.3 times on average) than those re-infected on global networks (2.1 times on average) (table 2).
On the local network (figure 5a), there are obvious clusters of high numbers of re-infection events (lightly shaded areas), and there is a low level of re-infection events elsewhere (dark shaded areas) (online version in colour). The circular clusters found by SatScan correlate with the areas highlighted in the heat map, implying that in local networks, individuals with multiple re-infection events are grouped with other re-infected individuals. In the more global networks (figure 5b,c), there are fewer clusters, with a lower number of re-infection events per individual in these clusters, and the distribution of re-infection events is more uniform. Results from the L-statistic confirm this interpretation (figure 5d–f).
Despite the anticipated population benefits of IPT, this intervention has only been implemented in a limited number of settings and in 2007, a mere 0.1 per cent of eligible individuals actually began a course of IPT . Possible reasons for the slow pace of adoption of IPT interventions include the challenges associated with administering daily single drug preventive therapy for a very large number of individuals without symptomatic TB disease  and the difficulty in ruling out the presence of active TB disease in HIV co-infected individuals . Increasing the uptake of IPT in communities with high HIV prevalence may require projections of the cost-effectiveness of such interventions and further demonstration that available diagnostic tools can effectively rule out subclinical TB disease before monotherapy is used. To encourage access and usage of IPT, the World Health Organization (WHO) recently revised IPT guidelines to suggest that the absence of symptoms (current cough, fever, weight loss or night sweats) is sufficient for ruling out active disease and initiating preventive therapy . They strongly recommend that eligible individuals complete a six month course of IPT, and in high-transmission areas, a 36 month course is advised.
However, recent studies aiming to assess the effectiveness and optimal duration of IPT have come to varying conclusions. A study in Botswana of HIV-positive individuals compared six and 36 months of IPT. Among those who tested positive for latent TB, in the six month cohort, the TB rate was 2.22 per 100 person years compared with 0.57 per 100 person years in the 36 month cohort. In individuals who tested negative for latent TB, there was little difference between the six and 36 month courses [7,35]. By contrast, a trial in India of HIV-positive patients without active TB concluded that six and 36 month courses of IPT were equally effective [35,36].
Our work suggests that differences between respiratory contact patterns and TB incidence between these settings may help explain why IPT strategies may have variable impact in different settings. While IPT treatment given to individuals prevents progression to disease, the impact of IPT at the community level depends on (i) the prevalence of TB infection among an individual's respiratory contacts, and (ii) the incidence of disease among those contacts, resulting in re-infection occurring after the IPT course is complete. The first of these is related to the extent of clustering of TB among shared groups of respiratory contacts, together with the overall prevalence of TB infection in the population. The second is mediated by the level of HIV/TB co-infection and the level of IPT coverage.
We have shown that the effect of IPT is highly variable in locally clustered networks compared with more globally connected ones. This is because in locally clustered networks, latent and incident TB infections occur in clusters, so that some areas suffer very high rates of incidence and re-infection, while other clusters are effectively protected from disease. In these protected clusters, treatment of a sufficient number of TB cases results in a reduction or elimination of onward TB transmission in the area. This occurs as a result of IPT when (i) a sufficient portion of the local TB burden is among HIV-infected individuals, and (ii) IPT coverage of eligible individuals is sufficiently high. We would therefore expect IPT to be most undermined by network and dynamical effects when TB infections are clustered among local groups of shared respiratory contacts, when HIV/TB co-infection is low among prevalent TB cases and/or when IPT coverage levels are too low to prevent onward transmission. This effect is exacerbated by local clustering of respiratory contacts, because it is in these circumstances that individuals who have been infected with TB once (whether or not they received IPT) are likely to be re-infected rapidly after their IPT treatment ends.
In terms of policy guidelines, our results imply that IPT will perform best if entire clusters of infected individuals at high risk of progression can be treated simultaneously. The benefits of IPT may be eroded if infection is not cleared from these clusters. This provides support to the conclusions from studies on clustering of TB cases in Africa , advising that cluster detection techniques (of both high and low case numbers) may help TB control programmes use their resources efficiently and effectively. Contact patterns and clustering are difficult to measure, particularly because relationships are dynamic and contacts may occur in casual settings where individuals are not known. However, techniques such as contact tracing from known TB cases have proved useful for identifying possible clusters (covering both past and current contacts) ; these efforts to identify clustered risk may be supported by social network analysis  and targeted investigation within settings of possible transmission such as bars .
It has previously been observed that clustering of contacts reduces disease spread over a network [39,40]; we also saw this effect in our simulations. Indeed, a high degree of clustering reduces disease spread compared with what would be expected from the initial growth rate of an epidemic or outbreak . Because epidemics spread first within a community and then (typically later) between them, networks with community structure can maintain epidemics at quite low incidence levels overall, but with a large number of generations of infectives over a relatively long period of time. By contrast, globally connected networks tend to have relatively fast, high-incidence, but short epidemics . This can make estimating the efficiency of control measures very difficult.
Our original hypothesis was that high levels of re-infection in locally clustered networks would undermine IPT programmes by allowing those individuals who had recently received IPT to be at risk of TB disease despite their treatment. However, instead we saw that the same mechanism that exposed some individuals to high levels of re-infection (clustering) also worked to protect other individuals from re-infection, via some clusters enjoying the near-elimination of TB after the introduction of IPT. As the protective effect of clustering is much weaker on globally connected networks, IPT performance on local networks can actually be better than on global ones at the population level, despite substantial regional variability.
The model is limited by the network structures we assumed. It is difficult to obtain comprehensive network data for populations, so our networks are caricatures of the structures of contact networks. In addition, for simplicity, we assumed that the networks were static; this assumption probably alters the transmission dynamics across the population. However, these limitations do not alter the overall qualitative insights from the model.
In this model, we have investigated single, fixed-length courses of IPT, which is consistent with current policy guidelines. We showed that the network structure affects both the frequency and distribution of TB re-infection, and thereby may modify the expected effectiveness of the impact of these types of IPT interventions. To mitigate the effect of re-infection, it may be useful to consider continuous or multiple courses of IPT, as suggested by some studies . These policies could potentially improve the impact of IPT by providing a longer period of protection from re-infection and progression, but there are concerns that longer or multiple courses of IPT may provide increased selective pressure and may facilitate the emergence of drug-resistant strains of TB (despite there being no evidence that IPT causes acquired resistance in individual patients ). In the recent publication of TB guidelines, the WHO concluded that IPT does not increase the risk of developing drug-resistant TB, so this should not be a barrier to its provision. Drug-resistant strains have been identified in most high-TB incidence areas [43–45], and IPT would be expected to exert selective pressure at the population level if given in the context of co-circulating resistance and sensitive strains [46,47]. Our model assumes only one, drug-sensitive strain of TB and does not consider the evolution of new drug-resistant strains or the effect on existing levels of resistance.
In conclusion, we have developed a new multi-network model for HIV and TB and find that local network structure can induce high levels of repeated re-infection that may undermine the projected effectiveness of IPT. At the same time, strong local clustering may provide a protective effect wherein TB is essentially eliminated from some areas and re-introduction of TB is limited by a scarcity of long-range connections. The combination of these dynamics means that IPT's effectiveness in clustered networks may be variable; this insight will be useful as we attempt to interpret and reconcile results from ongoing IPT trials.
T.C. was funded by NIH grants U54 GM088558-01 (MIDAS), R01 AI083036 and DP2OD006663 from the Office of the Director, US National Institutes of Health. The content is solely the responsibility of the authors and does not necessarily represent the official views of the Office of the Director of the US NIH, or the NIH. H.L.M. was supported by the Bristol Centre for Complexity Sciences and EPSRC grant EP/5011214.
Appendix A. Supplementary methods
The probability per time step of an individual with KHIV HIV-positive contacts becoming infected with HIV is equal to pHIV = 1 − (1 − αHIV)KHIV, where αHIV is the probability of transmission of HIV per infected contact, per month. The probability per time step of an individual becoming latently infected with TB is equal to pTB = 1 − (1−αTB)KTB(1 − εαTB)KTB,HIV, where αTB is the probability of transmission of TB and where the individual has KTB respiratory contacts with active TB and KTB,HIV sexual contacts with active TB. Sexual contacts are also respiratory contacts, so TB may be transmitted over the sexual contact network. The probability of transmission is reduced by ε over the sexual network to reflect the fact that these contacts are generally likely to be seen less often than those in the respiratory network.
An individual can progress from latent to active TB infection with a higher probability of progressing in the first 5 years after the initial infection . Individuals may be re-infected by infectious contacts, which is modelled as follows: if the re-infected individual is in the fast progression stage (within 5 years of initial infection), there is no change to its progression probability; however, if they are in the slow progression stage, they will have an elevated progression probability. Owing to some partial immunity acquired from the initial infection, this probability is lower than the fast progression probability. HIV-positive individuals have 25 per cent partial immunity and individuals not infected with HIV have 65 per cent partial immunity . The elevated progression probability caused by re-infection lasts for 5 years.
In the network model, transmission is more efficient on global networks, causing more active TB and more re-infection. To allow comparison of IPT on different networks in the same epidemiological setting, the transmission of TB was scaled between the more local and global networks so that the incidence of active TB was the same in each network.
Finally, to allow for the possibility of relapse and represent the immunity acquired after recovering from an active infection, recovery from active TB is modelled as moving to the slow progressing latent class.
If an individual dies, it is replaced in the following time step by a susceptible individual. This individual has the same respiratory contacts as the dead individual, but forms new contacts at random in the sexual network while retaining the sexual degree of the dead individual. The rearrangement of sexual contacts means old contacts do not immediately form a new sexual contact and the new individual is not automatically placed in an area of high HIV prevalence.
A.1. IPT and antiretroviral therapy
An IPT programme can be introduced in the simulation at any time step. The probability of an individual receiving IPT at any time step is such that over a year of the programme, a set percentage of those eligible will receive IPT (an eligible individual is HIV positive without active TB). The IPT course lasts for a set time period (here, nine months) during which the individual loses any latent TB infection and will not become infected again, although after the course of IPT is complete and they are susceptible to TB, they may become latently infected. An individual only receives IPT once in their lifetime, reflecting the current policy guidelines.
In the model, every individual who receives IPT also receives antiretroviral therapy (ART). It is assumed that ART reduces the infectivity of the individual by 80 per cent per contact (trial results vary between 67 and 92% [48–51]), reduces their mortality rate by 50 per cent and continues for the rest of their life.
↵1 SatScan is freely available software that was developed under the joint auspices of (i) Martin Kulldorff, (ii) the National Cancer Institute, and (iii) Farzad Mostashari of the New York City Department of Health and Mental Hygiene.
- Received March 22, 2011.
- Accepted March 30, 2011.
- This journal is © 2011 The Royal Society