Spatial patterns of tuberculosis and HIV co-infection in Ethiopia.

BACKGROUND: Tuberculosis (TB) and human immunodeficiency virus (HIV) are the leading causes of infectious-disease-related deaths in Ethiopia, but little is known about their spatial distribution across the country. The aim of this study was to investigate the spatial patterns of TB and HIV co-infection in Ethiopia at the district level.
METHODS: We conducted an ecological study using TB and HIV data reported from all regions of Ethiopia through the national Health Management Information System (HMIS), between June 2015 and June 2017. Spatial clustering was assessed using Moran's I statistic and Getis-Ord statistic. Spatial binomial regression models were constructed separately for the prevalence of TB among people living with HIV and for the prevalence of HIV among TB patients, with and without spatial components using a Bayesian approach.
RESULTS: A total of 1,830,880 HIV and 192,359 TB patients were included in the analysis. The prevalence of HIV among TB patients was 7.34%; hotspots were observed in districts located in Amhara, Afar, and Gambela regions, and cold spots were observed in Oromiya and Southern Nations, Nationalities, and People (SNNP) regions. The prevalence of TB among people living with HIV varied from 0.7% in Oromia region to 14.5% in Afar region. Hotspots of TB prevalence among people living with HIV were observed in districts located in Gambela, Afar, Somali, and Oromiya regions; whereas the cold spots were observed in districts located in Amhara and Tigray regions. The ecological-level factors associated with the prevalence of TB among people living with HIV were low wealth index (OR: 1.49; 95% CrI: 1.05, 2.05), low adult literacy rate (OR: 0.67; 95% CrI: 0.46, 0.94), and distance to an international border (OR: 0.61; 95% CrI: 0.40, 0.91). The factors associated with the prevalence of HIV among TB patients were poor health care access (OR: 0.76; 95% CrI: 0.59, 0.95), low wealth index (OR: 1.31; 95% CrI: 1.01, 1.67), and low adult literacy rate (OR: 1.37; 95% CrI: 1.03, 1.78).
CONCLUSION: Our study provides evidence for geographic clustering of TB/HIV co-infection in Ethiopia. Health care access, proximity to international borders, and demographic factors such as low wealth index and adult literacy were significantly associated with the prevalence of TB/HIV co-infection.

Introduction

Tuberculosis (TB) and human immunodeficiency virus (HIV) are the first and the second leading cause of death attributable to infectious diseases worldwide, respectively. The synergy between TB and HIV infection is strong. While TB is the leading cause of death among people living with HIV [1], HIV infection is the most important risk factor for developing active TB [2, 3]. Tuberculosis and HIV co-infection has been a major public health challenge to the global efforts in achieving elimination goals for both diseases. There were an estimated 36 million people living with HIV globally in 2017, with approximately one third of these people being co-infected with TB [4]. The highest burden of TB/HIV co-infection is reported from sub-Saharan Africa [4].

Ethiopia is one of the 30 high TB- and HIV-burden countries in the world [5-7], with incidence of 0.17 cases per 1000 population for HIV and 1.64 per 1000 for TB [8-10]. Studies conducted in Ethiopia have shown that TB/HIV co-infected individuals have a greater risk of psychosocial problems, low quality of life, and poorer physical health than HIV infected individuals without active TB [11-13]. The prevalence of TB/HIV co-infection has been found to vary widely in Ethiopia [14, 15] and other countries [16-18], partially because of differences in health care access and other socio-demographic factors such as wealth index and literacy rate.

Although several small scale studies have reported the epidemiology of TB/HIV co-infection in many parts of Ethiopia [19, 20], to our knowledge, there are no reports on the spatial distribution of TB-HIV co-infection nationally. Identifying areas where the burden of TB/HIV co-infection is concentrated might help identify populations at higher risk of co-infection, and knowledge of high and low burden areas is required for successful surveillance programs and optimal resource allocation [21], including identifying areas where existing TB and HIV collaborative programs and integrated services need to be strengthened. Therefore, this study aimed to assess the spatial distribution of TB/HIV co-infection in Ethiopia at the district level and investigate ecological-level factors associated with the prevalence of TB among people living with HIV and the prevalence of HIV among TB patients.

Methods

Study design and setting

An ecological study was conducted in Ethiopia. Our study estimated the prevalence of TB/HIV co-infection to the district level, using the national TB/HIV case notification data, reported from June 2015 to June 2017. The district is the third administrative level in Ethiopia. Districts comprise a well-defined population within a clearly delineated administrative and geographical area. Districts contain networks of primary health care units such as health centers, health posts and a district hospital. Health related data are compiled at the district level and are reported to the Zone Health Department, Region Health Bureau and then to the Federal Ministry of Health. Both TB and HIV are notifiable diseases in Ethiopia. These diseases are reported by the district Health Office to the Federal Ministry of Health through the Health Management Information System (HMIS) on a quarterly basis [22, 23]. The data include the number of: TB patients enrolled into the national TB program, TB cases who were tested for HIV, TB patients with HIV positive test results; HIV positive clients who received HIV care during the reporting period; HIV positive clients who were screened for TB during their visit; and HIV positive clients with active TB. TB and HIV diagnosis and treatment are integrated and are provided free of charge in Ethiopia across the country.

Data variables

The dependent variable in this study was numbers of TB and HIV co-infection in two study populations: people living with HIV and TB patients. This was determined as the prevalence of TB among people living with HIV and the prevalence of HIV among TB patients. Individual data on infection status with TB and HIV for the entire population or a representative sample were not available, precluding the use of a single (multinomial) model incorporating all possible combinations of TB and HIV mono- and co-infections.

All TB/HIV co-infected individuals reported to the national HMIS during the study period were included in the study. Study participants whose current address of residence was not registered, or unknown were excluded from the analysis because their data was unable to be georeferenced to the district level. The independent variables were socio-economic factors (such as low wealth index, average number of persons per room, the unemployment rate, and the adult literacy rate), behavioural factors (such as chat chewing and alcohol drinking), poor health care access (i.e. difficulty of getting advice or treatment due to lack of money, or distance to a health facility), distance to an international border (i.e. straight-line distance to the nearest international border), and climatic factors (such as mean temperature, precipitation, and aridity). A table showing the summary of independent variables, their sources of data and units of measurement is presented as supplementary information (S1 Table). These independent variables were extracted from the Ethiopian Demographic and Health Survey (EDHS 2011 and EDHS 2016), the Ethiopia Atlas of Population Density, and EDHS Spatial Data. The variables were selected a priori due to their plausible association with the prevalence of HIV/TB co-infection.

A Polygon shapefile for Ethiopia’s administrative boundaries at the district level was obtained from the Open Africa website (https://africaopendata.org/group/ethiopia). The TB and HIV datasets and the study covariates were merged with the district level polygons using ArcGIS (ESRI, Redlands, CA) geographical information system (GIS) software [24].

Data analysis

Descriptive statistics and geospatial analyses were conducted. The prevalence of TB among people living with HIV was calculated as the number of TB patients with an HIV positive test result in a district divided by the total number of clients enrolled in HIV care who were screened for TB during their visit in the same district. Similarly, the prevalence of HIV among TB patients was calculated as the number of TB patients with an HIV positive test result in a district divided by the total number of TB patients who were tested for HIV.

Spatial analysis

Spatial autocorrelation (based on the feature locations and attribute values) in prevalence of co-infection was measured using the Global Moran's I statistic in order to assess the overall pattern and trend of the data. Local spatial clustering was explored using the Getis-Ord statistic. An initial candidate set of demographics, behavioural, climatic and environmental covariates were considered for inclusion in the spatial models to investigate the ecological-level factors associated with the clustering of TB and HIV co-infection prevalence. Variable selection was made using fixed-effects univariable regression models and all variables with P<0.2 were selected for inclusion in final multivariable models. Variables were also checked for the presence of multi-collinearity, and those variables with a variance inflation factor greater than six were excluded from the final spatial model. The residuals of the final models were examined for the presence of spatial autocorrelation.

Spatial binomial regression models were constructed separately for the prevalence of TB among people living with HIV and the prevalence of HIV among TB patients in WinBUGS version 1.4 (MRC Biostatistics Unit, Cambridge, and Imperial College London, UK). Models were developed with and without spatial components using a Bayesian approach.

We constructed the model with a spatially structured and unstructured random effects. We assumed that the observed prevalence of TB among people living with HIV (r) at zone (i) had a binomial distribution with a total number of TB patients who were tested for HIV at the zone (ni) and the modelled prevalence of TB among people living with HIV at the zone (pi): ri ~ Binomial (ni, pi). For this model, pi was specified as follows: logit (p) = α+∑*X+U+V where α is the intercept; ∑*X is the matrix of m independent zone-specific variables X (i.e. proportion of the population with a low wealth index, alcohol drinking, poor health care access, adult literacy rate, and distance to an international border) measured at each zone i, multiplied by their coefficients (β). U are unstructured random effects; and V are the spatially structured random effects, modelled using a conditional autoregressive (CAR) approach. The CAR component was defined using an adjacency matrix to determine the spatial relationships between zones. The adjacency matrix for each zone was generated using ArcGIS based on the queen definition, whereby two areas are considered neighbours if they share a common boundary or vertex. A weight of 1 was given if two zones were neighbours and 0 if otherwise. The posterior parameters were estimated using a Bayesian Markov Chain Monte Carlo (MCMC) simulation. Non-informative priors were used for α (uniform prior with bounds—∞ and +∞) and the coefficients (normal prior with mean = 0 and precision 1× 10–6). The priors for the precision of the unstructured and spatially structured random effects were given non-informative gamma distributions with shape and scale parameters set at 0.001. The deviance information criterion (DIC) statistic was calculated to select the best-fitting model (models with a lower DIC statistic are considered to show a better compromise between model fit and parsimony). Convergence of the parameter values was observed using visual examination of history and density plots of the model chains. Convergence of the models was successfully achieved after 400,000 iterations. The models were run for a further 100,000 iterations and the posterior values were stored after 500,000 iterations. Using these procedures and steps, similar but separate models were constructed for the prevalence of HIV among TB patients. Models were also developed for each outcome excluding either the unstructured or structured random effects.

Ethics clearance

This study was approved by the Australian National University Human Research Ethics Committee (protocol number 2016/218) and a letter of endorsement was obtained from the Ethiopian Ministry of Health. As this study used district-level aggregated data, informed consent was not obtained from the study participants.

Results

A total of 1,830,880 HIV and 192,359 TB patients reported to the national HMIS from June 2015 to June 2017 were included in this study. These data covered all zones and regions of Ethiopia, except for 23 districts that had no TB or HIV data available. From the total 1,830,880 HIV positive clients who received HIV care in Ethiopia during the study period, 1,685,303 (92%) were screened for TB. A higher proportion of clients enrolled in HIV care who were screened for TB during the study period were female (61.6%; n = 1,128,917) than male (38.4%; n = 701,963), whereas a higher proportion of TB cases enrolled to Directly Observed Therapy Short Course (DOTS) and who were tested for HIV were male (54.8%; n = 105,387) than female (45.2%; n = 86,972). Table 1 presents the total number of TB and HIV patients in Ethiopia, stratified by sex. The overall prevalence of TB among people living with HIV was 0.84%; and varied greatly by sex, with a higher TB prevalence found among males (1.2%) than females (0.64%). The prevalence of HIV among TB patients was 7.34% and was slightly lower (7.09%) for males than females (7.72%). The highest prevalence of TB among people living with HIV was reported in Somalia (14.5%), Gambela (9.6%), and Afar (8.5%) regions (Table 2; Fig 1). The highest prevalence of HIV among TB patients was reported in Addis Ababa (22.5%), and the lowest prevalence was reported in Southern Nations, Nationalities, and People (SNNP) (3.8%) and Oromia (4.8%) regions (Table 2; Fig 2). The prevalence of TB among people living with HIV and the prevalence of HIV among TB patients at country, region, and zone levels are provided in the supplementary tables (S2 Table).

Table 1

The number of TB, HIV and TB/HIV patients reported in Ethiopia, 2015–2017.

Variable	Male	Female	Total
Total number of HIV positive clients who received HIV care during the reporting period	701,963	1,128,917	1,830,880
Number of clients enrolled in HIV care who were screened for TB during their visit	647,807	1,037,496	1,685,303
Total number of HIV positive clients with active TB	7,531	6,621	14,152
Number of TB cases enrolled to DOTS and who were tested for HIV	105,387	86,972	192,359
Number of TB patients with HIV positive test result	7,472	6,716	14,188
Number of TB patients with HIV Negative test result	96,917	79,589	176,506
Number of TB patients with unknown or undocumented HIV test result	998	667	1,665

HIV: Human Immunodeficiency Virus; TB: Tuberculosis; DOTS: Directly Observed Therapy, Short Course.

Table 2

The prevalence of TB among HIV patients, and the prevalence of HIV among TB patients in all regions of Ethiopia.

	TB prevalence among HIV patients (%)	HIV prevalence among TB patients (%)
Addis Ababa	1.0	22.5
Afar	8.5	8.1
Amhara	0.6	8.2
Beneshangul Gumuz	7.7	6.3
Dire Dawa	1.0	12.6
Gambela	9.6	20.8
Harari	0.9	8.9
Oromiya	0.7	4.8
SNNP	1.6	3.8
Somali	14.5	5.5
Tigray	0.9	8.4
Ethiopia	0.8	7.4

Fig 1

Geographical distribution of the prevalence of tuberculosis among people living with HIV (PLHIV) in Ethiopia, June 2016 and June 2017.

Fig 2

Geographical distribution of the prevalence of HIV among tuberculosis patients in Ethiopia, June 2016 and June 2017.

Spatial clustering of TB/HIV co-infection

There was evidence of spatial clustering of TB/HIV co-infection, with a global Moran’s I = 0.13 (p value < 0.001) for the prevalence of TB among people living with HIV and 0.17 (p value <0.001) for the prevalence of HIV among TB patients. The Getis-Ord statistics showed that the prevalence of HIV among TB patients in Ethiopia was clustered in districts located in Amhara region and the prevalence of HIV among TB patients in the same region is relatively widespread. Hotspots of TB prevalence among people living with HIV were observed in districts located in Gambela, Afar, Somali, and Oromiya regions; whereas the cold spots were observed in districts located in Amhara and Tigray regions (Fig 3). For the prevalence of HIV among TB patients, hotspots were observed in districts located in Amhara, Afar, and Gambela regions, and cold spots were observed in Oromiya and SNNP regions (Fig 4).

Fig 3

Spatial clustering of TB prevalence among people living with HIV (PLHIV) at district-level in Ethiopia, 2015 to 2017, based on the Getis-Ord Gi* statistic.

Fig 4

Spatial clustering of HIV prevalence among TB patients at district-level in Ethiopia, 2015 to 2017, based on the Getis-Ord Gi* statistic.

Factors associated with spatial clustering of TB/HIV co-infection

Table 3 shows the univariate models fitted for both the prevalence of TB among people living with HIV and the prevalence of HIV among TB patients. For both outcomes, the best-fitting model was the one with spatially structured random effects, but no spatially unstructured random effects. In the best-fitting model for the prevalence of TB among people living with HIV, low wealth index (odds ratio (OR): 1.49; 95% CrI: 1.05, 2.05), low adult literacy rate (OR: 0.67; 95% CrI: 0.46, 0.94), and proximity (distance) to international borders (OR: 0.61; 95% CrI: 0.40, 0.91) were significantly associated with TB (Table 4). When the model was fitted with the prevalence of HIV among TB patients, poor health care access (OR: 0.76; 95% CrI: 0.59, 0.95), low wealth index (OR: 1.31; 95% CrI: 1.01, 1.67), and low adult literacy rate (OR: 1.37; 95% CrI: 1.03, 1.78), were significantly associated with HIV prevalence (Table 4). When the residuals of the final model were examined for spatial autocorrelation, after incorporating the socio-economic and health care access variables into the model, spatial autocorrelation was apparent (Global Moran’s index = 0.09; p-value = 0.04), indicating that the covariates were not able to fully explain the spatial clustering of TB/HIV co-infections. Maps showing the residuals of structured random effects for HIV prevalence among TB patients and for TB prevalence among people living with HIV are presented in supplementary figures (S1 Fig and S2 Fig).

Table 3

Univariate model for prevalence of TB among HIV patients, and prevalence of HIV among TB patients in Ethiopia, 2015–2017.

	TB prevalence among HIV patients		HIV prevalence among TB patients
Independent variables	Coefficient (95% CI)	P-value	Coefficient (95% CI)	P-value
Low wealth index	0.12 (0.06, 0.18)	0.00	0.01 (-0.06, 0.09)	0.19
Average number of persons per room	2.32 (0.70, 3.95)	0.006	0.73 (-0.81, 2.29)	0.34
Unemployed population	-0.123 (-0.46, 0.21)	0.463	0.57 (0.22, 0.93)	0.002
Adult literacy rate	-0.14 (-0.24, -0.04)	0.004	0.04 (-0.05, 0.15)	0.15
Chat chewing	-0.007 (-0.09, 0.08)	0.868	-0.02 (-0.12, 0.08)	0.69
Alcohol drinking	-0.06 (-0.11, -0.007)	0.027	-0.02 (-0.08, 0.04)	0.18
Health care access problem	0.044 (-0.043, 0.13)	0.19	-0.09 (-0.19, 0.003)	0.058
Poor knowledge towards TB	0.13 (0.007, 0.26)	0.038	0.05 (-0.09, 0.19)	0.49
Poor attitude towards TB	0.104 (0.005, 0.203)	0.039	-0.05 (-0.17, 0.06)	0.34
Enhanced vegetation index	0.007 (-0.01, 0.01)	0.932	-0.008 (-0.03, 0.009)	0.34
Rainfall	-0.0005 (-0.002, 0.001)	0.54	0.01 (-0.03, 0.01)	0.24
Aridity	-0.005 (-0.01, 0.0009)	0.098	-0.005 (-0.01, 0.001)	0.10
Mean temperature	1.095 (0.63, 1.55)	0.001	0.53 (0.008, 1.05)	0.04
Proximity to international border	-0.02 (-0.03, -0.09)	0.001	-0.012 (-0.028, 0.06)	0.12

Table 4

Bayesian logistic regression model for the association of ecological level covariates and spatially structured and unstructured random effect with the prevalence of TB among HIV patients and prevalence of HIV among TB patients in Ethiopia, 2015 to 2017.

	TB prevalence among HIV patients			HIV prevalence among TB patients
	Spatial unstructured	Spatially structured	Both spatially structured and unstructured	Spatial unstructured	Spatially structured	Both spatially structured and unstructured
Independent variables	OR (95% CI)	OR (95% CI)	OR (95% CI)	OR (95% CI)	OR (95% CI)	OR (95% CI)
Low wealth index	1.53 (1.08, 2.10)	1.49 (1.05, 2.05)	1.53 (1.08, 2.10)	1.24 (0.94, 1.59)	1.31 (1.01, 1.67)	1.31 (1.01, 1.67)
Adult literacy rate	0.866 (0.62, 1.18)	0.67 (0.46, 0.94)	0.80 (0.55, 1.11)	1.27 (0.98, 1.62)	1.37 (1.03, 1.78)	1.36 (1.03, 1.78)
Alcohol drinking	0.69 (0.55, 1.85)	0.74 (0.55, 1.98)	0.71 (0.53, 0.92)	0.99 (0.83, 1.19)	1.10 (0.87, 1.37)	1.08 (0.84, 1.36)
Poor health care access	0.86 (0.67, 1.10)	0.80 (0.59, 1.07)	0.78 (0.58, 1.02)	0.78 (0.64, 0.95)	0.76 (0.59, 0.95)	0.74 (0.58, 0.93)
Proximity to international border	0.71 (0.55, 0.90)	0.61 (0.40, 0.91)	0.74 (0.53, 1.00)	0.90 (0.74, 1.09)	0.77 (0.55 (1.03)	0.82 (0.61, 1.08)
Intercept	-4.02 (-4.24, -3.80)	-4.06 (-4.16, -3.96)	-4.03 (-4.22, -4.03)	-2.72 (-2.89, -2.55)	-2.72 (-2.78, -2.67)	-2.72 (-2.82, -2.62)
σ2 (variance of spatially unstructured random effect)	1.21 (0.82, 1.68)	-	2.74 (0.96, 10.33)	1.69 (1.16, 2.33)	-	24.75 (1.80, 134.6)
σ2 (variance of spatial structured random effect)	-	0.32 (0.21, 0.45)	36.29 (0.34, 381.1)	-	0.50 (0.34, 0.69)	9.54 (0.39, 16.08)
DIC value	618.1	298.7	1230.0	673.2	494.9	495.2

Independent variables were standardised to have mean = 0 and standard deviation = 1; CI, Bayesian credible interval; OR, odds ratio; DIC, deviance information criterion.

Discussion

TB and HIV co-infection is a major public health concern in Ethiopia. The present study examined the spatial distribution of TB/HIV co-infection in the whole of Ethiopia at the district level, by fitting two separate models for the prevalence of TB among people living with HIV and the prevalence of HIV among TB patients. This study found spatial clustering of HIV/TB co-infection in Ethiopia was associated with poor health care access, distance to international borders and socio-economic factors.

The prevalence of TB among people living with HIV varied strongly between districts and zones in Ethiopia, with prevalence as high as 35.9% in Etang Special zone and as low as 0.1% in Burji Special zone. The spatial heterogeneity in the prevalence of TB among people living with HIV was significantly associated with underlying variability in area-level socio-economic variables such as low wealth index, low adult literacy rate, and close distance to the international border. The observed association between the prevalence of TB and socio-economic factors has already been described in the existing literature [25-27]. A high prevalence of TB has previously been reported among people with low wealth quintile and low educational status in several countries such as South Africa [28], Kenya [29], China [30], and Brazil [31].

An interesting finding in our study was the association between the prevalence of TB among people with HIV and distance to the international border; specifically that a higher prevalence of TB was reported in the districts near to the international border as compared to districts located far from the international border. The map in our local spatial clustering analysis also showed a similar finding, in which hot-spot districts are located in Gambela, Somali, Afar, and Oromiya regions near the Ethiopia-Sudan, Ethiopia-Eritrea and Ethiopia-Kenya borders. These districts are characterised by high numbers of migrants and poor access to health care services. Our previous studies conducted in northwest Ethiopia to assess the spatial clustering of paediatric TB and MDR-TB in the general population have also detected similar hot-spot areas near the Ethiopia-Sudan and Ethiopia-Eritrea borders [32, 33]. These findings suggest that effective public health interventions are urgently needed in these remote settings to control TB transmission in the border regions.

Our study also showed that the overall prevalence of HIV infection among TB patients was 7.4%; and while this figure is considerably higher than the national HIV prevalence among the general population (0.9%) [34], it is in line with the WHO estimates of HIV prevalence among TB patients for Ethiopia (7%) [7]. A high prevalence of HIV has been consistently reported among TB patients in previous studies [35, 36]. A recent systematic review and meta-analysis conducted in sub-Saharan Africa showed a high prevalence (31.8%) of HIV among TB patients, with a gradual decline in prevalence from 33.7% in the period before 2000 to 25.7% in the period after 2010 [37].

Our study also identified wide variation in the prevalence of HIV among TB patients at the sub-national level, ranging from 0.8% in East Hararge zone to 57.1% in Liben zone. A hotspot of HIV prevalence among TB patients was observed in districts located in Amhara, Afar, and Gambela regions. This indicates that the overall national prevalence of HIV may mask the local spatial heterogeneity of the HIV epidemic and this highlights a clustered HIV transmission pattern at the district level within the country. Previous studies have also found spatial clustering of HIV prevalence among TB patients at the local level [38, 39]. Such findings suggest that the burden of HIV among TB patients should be monitored at the sub-national level in order to effectively utilize prevention strategies and strengthen integrating services. This is particularly important in a country such as Ethiopia, where the burden of HIV and TB is high.

The spatial clustering of HIV prevalence among TB patients was significantly associated with health care access and socio-economic factors such low wealth index and low adult literacy rate. This finding is consistent with previous studies conducted in sub-Saharan Africa where a high prevalence of HIV infection was reported in areas with a low socio-economic index [40, 41]. The negative association between poor health care access and prevalence of HIV among TB patients likely relates to access to diagnostic facilities. Although HIV counselling and testing are provided in Ethiopia at all hospitals and health centres, their coverage is poor in rural and remote areas [42]. This suggests that identifying geographic areas with low wealth index and providing adequate health care services would be vital to improve local HIV identification, prevention and control programs.

The observed spatial clustering of TB and HIV co-infection in Ethiopia generally suggests that although country-level estimates are important for international comparison, they may mask important local geographic variations in the prevalence of both TB and HIV. This indicates that TB and HIV are poorly controlled in some districts. Unless effective interventions are designed and put in place, these hotspot districts could continue to be a source of TB and HIV transmission and create operational challenges for the design and implementation of collaborative TB-HIV interventions. Thus, a better understanding of the subnational variations of TB/HIV co-infection will help with the scale up of WHO-recommended TB/HIV collaborative activities to address the impact of TB/HIV co-infection and to achieve the national and global targets [31].

Our analysis had a large sample size and wide geographical coverage across Ethiopia. However, it has some limitations. First, the study was conducted based on routinely collected data that are available from the national HMIS. These data may not reflect the exact national figures as under-detection and under-reporting of cases is a potential limitation of the data. Second, since the unit of spatial analysis was districts, we were unable to estimate small scale spatial variation in HIV and TB prevalence at the community level. However, we believe that our results provide a useful operational framework for planning and setting priorities as districts in Ethiopia are the most important administrative level where local resources are allocated for each sector. Third, it would be more meaningful if we assessed the spatial distribution of TB, HIV and TB/HIV co-infection, using a multinomial model with general population as the denominator. The multinomial outcomes would be not infected, HIV mono-infection, TB mono-infection and HIV/TB co-infection. However, we were unable to do this due to the lack of population level screening data; our cohort was restricted to individuals who received secondary screening after receiving a positive diagnosis for TB or HIV. Lastly, although we used several ecological level datasets, some important covariates such as drug resistant patterns were not included in the analysis. These data were not available in the HMIS and further analysis should consider their inclusion when they become available.

Conclusion

Our study provides evidence for geographic clustering of TB/HIV co-infection in Ethiopia and identified priority regions, zones and districts where TB/HIV collaborative activities should be strengthened to achieve greater impact with minimal additional investment.

Summary of independent variables, sources of data and definition of variables.

(DOCX)

Click here for additional data file.

The prevalence of TB, HIV and TB/HIV co-infection in Ethiopia at country, region, and zone levels, between June 2016 and June 2017.

(DOCX)

Click here for additional data file.

Resdiduals of structured random effects for the prevalence of HIV among tuberculosis patients in Ethiopia, June 2016 and June 2017.

(TIF)

Click here for additional data file.

Resdiduals of structured random effects for the prevalence of tuberculosis among people lving with HIV (PLHIV) in Ethiopia, June 2016 and June 2017.

(TIF)

Click here for additional data file.

Transfer Alert

This paper was transferred from another journal. As a result, its full editorial history (including decision letters, peer reviews and author responses) may not be present.

3 Nov 2019

PONE-D-19-28080

Spatial patterns of tuberculosis and HIV co-infection in Ethiopia

PLOS ONE

Dear Dr. Alene,

Thank you for submitting your manuscript to PLOS ONE. After careful consideration, we feel that it has merit but does not fully meet PLOS ONE’s publication criteria as it currently stands. Therefore, we invite you to submit a revised version of the manuscript that addresses the points raised during the review process.

We would appreciate receiving your revised manuscript by Dec 18 2019 11:59PM. When you are ready to submit your revision, log on to https://www.editorialmanager.com/pone/ and select the 'Submissions Needing Revision' folder to locate your manuscript file.

If you would like to make changes to your financial disclosure, please include your updated statement in your cover letter.

To enhance the reproducibility of your results, we recommend that if applicable you deposit your laboratory protocols in protocols.io, where a protocol can be assigned its own identifier (DOI) such that it can be cited independently in the future. For instructions see: http://journals.plos.org/plosone/s/submission-guidelines#loc-laboratory-protocols

Please include the following items when submitting your revised manuscript:

A rebuttal letter that responds to each point raised by the academic editor and reviewer(s). This letter should be uploaded as separate file and labeled 'Response to Reviewers'.

A marked-up copy of your manuscript that highlights changes made to the original version. This file should be uploaded as separate file and labeled 'Revised Manuscript with Track Changes'.

An unmarked version of your revised paper without tracked changes. This file should be uploaded as separate file and labeled 'Manuscript'.

Please note while forming your response, if your article is accepted, you may have the opportunity to make the peer review history publicly available. The record will include editor decision letters (with reviews) and your responses to reviewer comments. If eligible, we will contact you to opt in or out.

We look forward to receiving your revised manuscript.

Kind regards,

HASNAIN SEYED EHTESHAM

Academic Editor

PLOS ONE

Journal Requirements:

1. When submitting your revision, we need you to address these additional requirements. Please ensure that your manuscript meets PLOS ONE's style requirements, including those for file naming. The PLOS ONE style templates can be found at http://www.journals.plos.org/plosone/s/file?id=wjVg/PLOSOne_formatting_sample_main_body.pdf and http://www.journals.plos.org/plosone/s/file?id=ba62/PLOSOne_formatting_sample_title_authors_affiliations.pdf

2. In ethics statement in the manuscript and in the online submission form, please provide additional information about the patient records used in your retrospective study. Specifically, please ensure that you have discussed whether all data were fully anonymized before you accessed them and/or whether the IRB or ethics committee waived the requirement for informed consent. If patients provided informed written consent to have data from their medical records used in research, please include this information.

3. We note you have included a table to which you do not refer in the text of your manuscript. Please ensure that you refer to Tables 3 and 4 in your text; if accepted, production will need this reference to link the reader to the Table.

Additional Editor Comments:

Minor Revision

[Note: HTML markup is below. Please do not edit.]

Reviewers' comments:

Reviewer's Responses to Questions

Comments to the Author

1. Is the manuscript technically sound, and do the data support the conclusions?

The manuscript must describe a technically sound piece of scientific research with data that supports the conclusions. Experiments must have been conducted rigorously, with appropriate controls, replication, and sample sizes. The conclusions must be drawn appropriately based on the data presented.

Reviewer #1: Yes

Reviewer #2: Yes

**********

2. Has the statistical analysis been performed appropriately and rigorously?

Reviewer #1: Yes

Reviewer #2: Yes

**********

3. Have the authors made all data underlying the findings in their manuscript fully available?

The PLOS Data policy requires authors to make all data underlying the findings described in their manuscript fully available without restriction, with rare exception (please refer to the Data Availability Statement in the manuscript PDF file). The data should be provided as part of the manuscript or its supporting information, or deposited to a public repository. For example, in addition to summary statistics, the data points behind means, medians and variance measures should be available. If there are restrictions on publicly sharing data—e.g. participant privacy or use of data from a third party—those must be specified.

Reviewer #1: Yes

Reviewer #2: Yes

**********

4. Is the manuscript presented in an intelligible fashion and written in standard English?

PLOS ONE does not copyedit accepted manuscripts, so the language in submitted articles must be clear, correct, and unambiguous. Any typographical or grammatical errors should be corrected at revision, so please note any specific errors here.

Reviewer #1: Yes

Reviewer #2: Yes

**********

5. Review Comments to the Author

Please use the space provided to explain your answers to the questions above. You may also include additional comments for the author, including concerns about dual publication, research ethics, or publication ethics. (Please upload your review as an attachment if it exceeds 20,000 characters)

Reviewer #1: Comments to the authors:

Tuberculosis (TB) is the leading cause of death for individuals infected with Human immunodeficiency virus (HIV). Conversely, HIV is the most important risk factor in the progression of TB from the latent to the active status. In order to manage this double epidemic situation, an integrated approach that includes HIV management in TB patients was proposed by the World Health Organization and was implemented in many countries world-wide. To enable targeted intervention using the integrated approach, areas with high disease prevalence rates for TB and HIV need to be identified. The authors of the manuscript identified the hot and cold spot regions of HIV/TB co-infection in the districts of Ethiopia. They used a large data set to get inference on HIV/TB co-infection. In this regard, particularly in Ethiopia, the data regarding the HIV/TB co-infection are scare. This study successfully identified co-infection hot spots and associated ecological variables. The study could be used in the future to effectively manage HIV/TB co-infection in the Ethiopia. The present study is well designed, executed and the inferences drawn are well supported by the data. The manuscript is well written and understandable. Although, the concept of the study is not new, but could provide a reference for the management of TB/HIV co-infection in Ethiopia. The manuscript could be accepted for publication after incorporating the given suggestions.

Abstract:

1. The authors only mentioned hot spots for TB/HIV or HIV/TB con-infection in the districts of Ethiopia. They are advised to include cold spot regions in the abstract to distinguish between highly prevalent regions of co-infection to the least one.

Introduction:

2. The authors should also include in the introduction section, the earlier studies reporting the spatial distribution of co-infection of HIV/TB in other countries/regions. For eg. Aturinde A et al. 2019, Peres DA et al. 2019, Wei W et al. 2016.

Results:

3. The authors should include Figure S1 and Figure S2 in the main manuscript.

4. Supplementary Table S2 was incorrectly labeled as supplementary table S1. The authors should correct the labeling mistake.

Reviewer #2: The paper entitled ‘Spatial patterns of tuberculosis and HIV co-infection in Ethiopia’ by Clements et al is a comparative study of the geospatial distribution of prevalence of TB/HIV co-infection in Ethiopia. Additionally, the authors have investigated the ecological factors like health care access, distance to international borders and socioeconomic factors associated with the prevalence of HIV among people living with TB and vice versa at the district level. Furthermore, the authors found gender bias in the prevalence of TB in people living with HIV being prevalent in males than in females, whereas prevalence of HIV among TB patients was slightly lower for males than females.

The study is elaborative and very well designed with large sample size and wide geographical coverage across Ethiopia. The paper is well-written with minor syntax mistakes (which require minor corrections). Font throughout the article including the references should remain same all across. The technical procedures are correct and the results are clear-cut. Limitations of the study have been put forward by the authors very clearly. On the whole, the data is convincing and the results significantly add up to the already existing information on the particular area of research.

**********

6. PLOS authors have the option to publish the peer review history of their article (what does this mean?). If published, this will include your full peer review and any attached files.

If you choose “no”, your identity will remain anonymous but your review may still be made public.

Do you want your identity to be public for this peer review? For information about this choice, including consent withdrawal, please see our Privacy Policy.

Reviewer #1: No

Reviewer #2: No

[NOTE: If reviewer comments were submitted as an attachment file, they will be attached to this email and accessible via the submission site. Please log into your account, locate the manuscript record, and check for the action link "View Attachments". If this link does not appear, there are no attachment files to be viewed.]

While revising your submission, please upload your figure files to the Preflight Analysis and Conversion Engine (PACE) digital diagnostic tool, https://pacev2.apexcovantage.com/. PACE helps ensure that figures meet PLOS requirements. To use PACE, you must first register as a user. Registration is free. Then, login and navigate to the UPLOAD tab, where you will find detailed instructions on how to use the tool. If you encounter any issues or have any questions when using PACE, please email us at figures@plos.org. Please note that Supporting Information files do not need this step.

Submitted filename: PONE-D-19-28080_reviewer.docx

Click here for additional data file.

11 Nov 2019

Manuscript Number: PONE-D-19-28080

Article Title: Spatial patterns of tuberculosis and HIV co-infection in Ethiopia

Authors: Kefyalew Addis Alene, Kerri Viney, Hannah C Moore, Maereg Wagaw, Archie CA Clements

Journal: PLOS ONE

Response to Editor’s Comments:

1. When submitting your revision, we need you to address these additional requirements. Please ensure that your manuscript meets PLOS ONE's style requirements, including those for file naming.

Response: We followed the PLOS ONE style templates when submitting the revised version of the manuscript.

2. In ethics statement in the manuscript and in the online submission form, please provide additional information about the patient records used in your retrospective study. Specifically, please ensure that you have discussed whether all data were fully anonymized before you accessed them and/or whether the IRB or ethics committee waived the requirement for informed consent. If patients provided informed written consent to have data from their medical records used in research, please include this information.

Response: We have now mentioned (on page 7 line 177-178) that “This study was approved by the Australian National University Human Research Ethics Committee (protocol number 2016/218) and a letter of endorsement was obtained from the Ethiopian Ministry of Health. As this study used district-level aggregated data, informed consent was not obtained from the study participants.”

3. We note you have included a table to which you do not refer in the text of your manuscript. Please ensure that you refer to Tables 3 and 4 in your text; if accepted, production will need this reference to link the reader to the Table.

Response: Tables are now appropriately referred in the text.

Response to Reviews’ Comments

Reviewer #1: Tuberculosis (TB) is the leading cause of death for individuals infected with Human immunodeficiency virus (HIV). Conversely, HIV is the most important risk factor in the progression of TB from the latent to the active status. In order to manage this double epidemic situation, an integrated approach that includes HIV management in TB patients was proposed by the World Health Organization and was implemented in many countries world-wide. To enable targeted intervention using the integrated approach, areas with high disease prevalence rates for TB and HIV need to be identified. The authors of the manuscript identified the hot and cold spot regions of HIV/TB co-infection in the districts of Ethiopia. They used a large data set to get inference on HIV/TB co-infection. In this regard, particularly in Ethiopia, the data regarding the HIV/TB co-infection are scare. This study successfully identified co-infection hot spots and associated ecological variables. The study could be used in the future to effectively manage HIV/TB co-infection in the Ethiopia. The present study is well designed, executed and the inferences drawn are well supported by the data. The manuscript is well written and understandable. Although, the concept of the study is not new, but could provide a reference for the management of TB/HIV co-infection in Ethiopia. The manuscript could be accepted for publication after incorporating the given suggestions.

Abstract:

1. The authors only mentioned hot spots for TB/HIV or HIV/TB con-infection in the districts of Ethiopia. They are advised to include cold spot regions in the abstract to distinguish between highly prevalent regions of co-infection to the least one.

Response: It is now included in the revised version of the manuscript as “…hotspots were observed in districts located in Amhara, Afar, and Gambela regions, and cold spots were observed in Oromiya and Southern Nations, Nationalities, and People (SNNP) regions. …. Hotspots of TB prevalence among people living with HIV were observed in districts located in Gambela, Afar, Somali, and Oromiya regions; whereas the cold spots were observed in districts located in Amhara and Tigray regions.”

Introduction:

2. The authors should also include in the introduction section, the earlier studies reporting the spatial distribution of co-infection of HIV/TB in other countries/regions. For eg. Aturinde A et al. 2019, Peres DA et al. 2019, Wei W et al. 2016.

Response: These recent studies are now cited in the revised version of the manuscript on page 3 line 74.

Results:

3. The authors should include Figure S1 and Figure S2 in the main manuscript.

Response: Figure S1 and Figure S2 are now included in the main manuscript on page 16 and 17.

4. Supplementary Table S2 was incorrectly labeled as supplementary table S1. The authors should correct the labeling mistake.

Response: This is now corrected.

Reviewer #2: The paper entitled ‘Spatial patterns of tuberculosis and HIV co-infection in Ethiopia’ by Clements et al is a comparative study of the geospatial distribution of prevalence of TB/HIV co-infection in Ethiopia. Additionally, the authors have investigated the ecological factors like health care access, distance to international borders and socioeconomic factors associated with the prevalence of HIV among people living with TB and vice versa at the district level. Furthermore, the authors found gender bias in the prevalence of TB in people living with HIV being prevalent in males than in females, whereas prevalence of HIV among TB patients was slightly lower for males than females. The study is elaborative and very well designed with large sample size and wide geographical coverage across Ethiopia. The paper is well-written with minor syntax mistakes (which require minor corrections). Font throughout the article including the references should remain same all across. The technical procedures are correct and the results are clear-cut. Limitations of the study have been put forward by the authors very clearly. On the whole, the data is convincing and the results significantly add up to the already existing information on the particular area of research.

Response: We greatly appreciate the reviewer’s efforts to carefully review the paper. All the errors and font styles are now corrected in the revised version of the manuscript.

21 Nov 2019

Spatial patterns of tuberculosis and HIV co-infection in Ethiopia

PONE-D-19-28080R1

Dear Dr. Alene,

We are pleased to inform you that your manuscript has been judged scientifically suitable for publication and will be formally accepted for publication once it complies with all outstanding technical requirements.

Within one week, you will receive an e-mail containing information on the amendments required prior to publication. When all required modifications have been addressed, you will receive a formal acceptance letter and your manuscript will proceed to our production department and be scheduled for publication.

Shortly after the formal acceptance letter is sent, an invoice for payment will follow. To ensure an efficient production and billing process, please log into Editorial Manager at https://www.editorialmanager.com/pone/, click the "Update My Information" link at the top of the page, and update your user information. If you have any billing related questions, please contact our Author Billing department directly at authorbilling@plos.org.

If your institution or institutions have a press office, please notify them about your upcoming paper to enable them to help maximize its impact. If they will be preparing press materials for this manuscript, you must inform our press team as soon as possible and no later than 48 hours after receiving the formal acceptance. Your manuscript will remain under strict press embargo until 2 pm Eastern Time on the date of publication. For more information, please contact onepress@plos.org.

With kind regards,

HASNAIN SEYED EHTESHAM

Academic Editor

PLOS ONE

Additional Editor Comments (optional):

The Authors have revised the manuscript by addressing all the comments and questions raised by the Reviewers. New references have been added. Supplementary figures 1 and 2 have now been added to the main manuscript. Corrections to table numbering have been made.

I recommend this manuscript for publication.

Reviewers' comments:

26 Nov 2019

PONE-D-19-28080R1

Spatial patterns of tuberculosis and HIV co-infection in Ethiopia

Dear Dr. Alene:

I am pleased to inform you that your manuscript has been deemed suitable for publication in PLOS ONE. Congratulations! Your manuscript is now with our production department.

If your institution or institutions have a press office, please notify them about your upcoming paper at this point, to enable them to help maximize its impact. If they will be preparing press materials for this manuscript, please inform our press team within the next 48 hours. Your manuscript will remain under strict press embargo until 2 pm Eastern Time on the date of publication. For more information please contact onepress@plos.org.

For any other questions or concerns, please email plosone@plos.org.

Thank you for submitting your work to PLOS ONE.

With kind regards,

PLOS ONE Editorial Office Staff

on behalf of

Prof HASNAIN SEYED EHTESHAM

Academic Editor

PLOS ONE