Abstract

Background

A substantial proportion of persons on antiretroviral therapy (ART) considered lost to follow-up have actually transferred their human immunodeficiency virus (HIV) care to other facilities. However, the relationship between facility switching and virologic outcomes, including viral rebound, is poorly understood.

Methods

We used data from 40 communities (2015–2020) in the Rakai Community Cohort Study to estimate incidence of facility switching and viral rebound. Persons aged 15–49 years with serologically confirmed HIV who self-reported ART use and contributed ≥1 follow-up visit were included. Facility switching and virologic outcomes were assessed between 2 consecutive study visits (ie, index and follow-up visits, interval of approximately 18 months). Those who reported different HIV treatment facilities between index and follow-up study visits were classified as having switched facilities. Virologic outcomes included viral rebound among individuals initially suppressed (<200 copies/mL). Multivariable Poisson regression was used to estimate associations between facility switching and viral rebound.

Results

Overall, 2257 persons who self-reported ART use (median age, 35 years; 65% female, 92% initially suppressed) contributed 3335 visit-pairs and 5959 person-years to the analysis. Facility switching was common (4.8 per 100 person-years; 95% confidence interval [CI], 4.2–5.5) and most pronounced in persons aged <30 years and fishing community residents. Among persons suppressed at their index visit (n = 2076), incidence of viral rebound was more than twice as high in persons who switched facilities (adjusted incidence rate ratio = 2.27; 95% CI, 1.16–4.45).

Conclusions

Facility switching was common and associated with viral rebound among persons initially suppressed. Investments in more agile, person-centered models for mobile clients are needed to address system inefficiencies and bottlenecks that can disrupt HIV care continuity.

Universal antiretroviral therapy (ART) provision through a “Treat All” approach has improved human immunodeficiency virus (HIV) outcomes throughout Africa. In the “Treat All” era, ART coverage in settings with generalized HIV epidemics such as Uganda has exceeded 90%, accelerating gains in population-level viral load suppression (VLS) [1–3]. Given the success of mass HIV treatment scale-up efforts, focus has been redirected toward identifying and addressing the unsuppressed fraction of persons with HIV, including those experiencing HIV treatment interruptions or defaulting from care altogether. A growing body of literature has attempted to characterize care attrition in persons with HIV, deemed lost to follow-up from HIV services [4]. This broadly includes individuals who miss appointments or ART refills and are not successfully traced by the health facility.

Nevertheless, some studies have shown that sizeable numbers of these untraced persons have not defaulted from HIV treatment but rather transferred care to other health facilities [5]. In 2 South African studies, more than a third of clients who transferred clinics were misclassified as lost to follow-up [6, 7]. Available evidence suggests that facility switching, both planned and unplanned, is motivated by numerous factors, including stigma/privacy concerns, care quality, travel distances, and migration [8–10]. Whether formal (through clinic-facilitated transfer processes) or informal (through client-initiated self-transfers unregistered by the health system, also called “silent transfers”), facility switching has the potential to disrupt HIV care continuity, with some studies attributing prolonged treatment gaps to clinic transfers [11–13].

Despite the frequency of facility transfers in many contexts and their potential to destabilize care, few studies have examined HIV clinical outcomes in persons who switch facilities, and none have leveraged population-level data. The available literature presents mixed evidence on whether facility switching improves or worsens virologic outcomes. In Malawi, 3-year survival among patients on ART who transferred their care to other facilities was comparable to that of the nontransferring population [14]. Conversely, other studies of South African clinical cohorts have reported suboptimal VLS outcomes following facility transfers [15, 16]. Importantly, these studies reported on facility switches documented by the health system and include time periods that predate “Treat All” scale-up. Contemporary population-level data, which capture both formal and silent facility transfers, can better illuminate switching dynamics across independent facility networks and how these dynamics relate to virologic outcomes in the “Treat All” era.

We nested a study of health facility switching in a prospective HIV surveillance cohort in southern Uganda to examine the longitudinal relationship between facility switching and viral load status among persons who self-reported ART use. Estimating the incidence of facility switching and its relationship to longitudinal virologic outcomes can guide differentiation of HIV services to clients who move between facilities, optimizing their treatment outcomes.

METHODS

Design

Data were derived from the Rakai Community Cohort Study (RCCS), an open, population-based HIV surveillance cohort in 40 communities in and around Rakai District, Uganda. Located in south-central Uganda, the RCCS surveillance area consists of rural (agrarian) and semiurban (trading) inland communities with moderate HIV prevalence (9%–26%) as well as hyperendemic fishing communities (HIV prevalence, 37%–43%) along Lake Victoria [17]. RCCS communities are also characterized by substantial in- and out-migration, including mobility within the study area as well as migration to and from communities outside Rakai [18].

Individuals aged ≥15 years across censused communities complete structured interviews approximately every 18 months for assessment of household characteristics, sexual behaviors, and HIV service utilization. HIV status is ascertained among those aged 15–49 years using a validated algorithm of 3 rapid tests and confirmatory enzyme immunoassay testing [19]. Viral load testing is then performed on stored plasma of persons with serologically confirmed HIV using the Abbott RealTime HIV-1 assay (Abbott Molecular, Inc, Des Plaines, IL).

For this study, inclusion was restricted to persons with HIV who self-reported ART use over ≥2 consecutive visits during the following 3 survey intervals: round 17 (February 2015 to September 2016), round 18 (October 2016 to May 2018), and round 19 (June 2018 to October 2020).

Measures

The unit of analysis was a visit-pair, defined as 2 consecutive study visits during the observation period. Each visit-pair consisted of an index visit and a follow-up visit. Participants could, therefore, contribute up to 2 visit-pairs to the analysis, where an individual’s follow-up visit in the first visit-pair served as the index visit in their second visit-pair.

Longitudinal Virologic Outcomes

Four categorical viral load outcomes were identified within visit-pairs, each measuring the following longitudinal virologic patterns over visit intervals: durable suppression (VLS across visits), new suppression (VLS at follow-up visit only), viral rebound (VLS at index visit only), and persistent viremia (VLS at neither visit). A VLS cut point of <200 HIV RNA copies/mL was used to maximize detection of persons who exhibited lower-level viremia and minimize misclassification of persons who experienced episodic, inconsequential viremic “blips” as unsuppressed [20]. ART use (current, previous, never) is self-reported in the RCCS and is a valid proxy for treatment coverage relative to assay-based drug detection methods (sensitivity, 77%; specificity, 99%) [21].

Health Facility Switching

In Uganda, HIV care is provided at level 2 facilities (ie, parish level) and higher. Clinics or hospitals attended for HIV care, both within and outside the study area, were enumerated at each study visit for persons currently on ART. Participants were considered to have switched facilities if they reported HIV care enrollment at 2 distinct health facilities over the visit-pair. By comparison, individuals were classified as not having switched facilities (“stable attendance”) if their designated health facility remained constant over the visit-pair.

Distance From Households to Health Facilities

Global positioning system (GPS) coordinates (latitude and longitude) were captured for censused households and enumerated health facilities within the RCCS surveillance area using handheld GPS devices. Health facilities outside the study area were identified and geocoded using a combination of sources, including the Uganda Ministry of Health’s National Health Facility Master List [22], a spatial database of public health facilities in sub-Saharan Africa [23], and Google Earth Pro version 9 (Google LLC, Mountain View, CA). The distance (in kilometers) from an individual’s household to their designated HIV care facility was measured using a geodesic distance calculation, which estimates the shortest spherical distance between 2 points along the earth’s surface. Although road network distance calculations measure travel distances more accurately than geodesic distance, numerous households in the study area were not connected to existing road networks and, thus, could not be mapped onto existing transit infrastructure [24]. Geodesic distances do not vary substantially from travel distances in this context and are, thus, appropriate proxies for travel distances from households to facilities [25].

To estimate the relative increase or decrease in distances to health facilities that resulted from facility switching, distances (in kilometers) between households and health facilities estimated at follow-up (d2) were subtracted from distances between households and facilities at the index visit (d1). The final travel distance function (d1d2) corresponded to the change in travel distance attributed to switching, with positive values (d1d2 > 0) representing switches to facilities closer to households and negative values (d1d2 < 0) representing switches to facilities further away from households.

Sociodemographics

Other measured independent variables included age, sex, marital status, educational attainment, primary occupation, religion, household wealth (ascertained from a locally validated index [26] of asset ownership), in-migration (immigration into the study area since the last census), community type, and intercommunity mobility (residence in a different community at follow-up). Among females, pregnancy was ascertained at the index and follow-up visits through self-report and urine dipstick screening. Excluding intercommunity mobility and pregnancy, sociodemographic variables were assigned time-invariant values from the participant’s index visit.

Statistical Analyses

Data were managed and analyzed using Stata/IC 15.1 (StataCorp LLC, College Station, TX). Facility switching incidence rates were estimated, overall and stratified by age, sex, and community type, by dividing the total number of switching events by the sum of total person-years contributed within each visit-pair (accrued person-time between index and follow-up visit, represented by visit dates). Events were assumed to occur at the midpoint of visit-pairs.

Next, Poisson regression with generalized estimating equations, independent covariance matrices, and robust standard errors (accounting for nonindependence of visit-pairs contributed by the same individual) estimated associations of facility switching with new VLS among initially viremic persons and viral rebound among initially suppressed persons, respectively, reported as adjusted incidence rate ratios (adjIRR) with 95% confidence intervals (CIs). Visit-pairs that exhibited persistent viremia and durable VLS served as the reference groups for new suppression and viral rebound, respectively. A best-fitting working correlation structure was selected using the quasi-likelihood under the independence model criterion [27] (Supplementary Table 1). Statistical models were first fit for the overall sample, then stratified by sex and community type. Multivariable models adjusted for sociodemographics were associated with facility switching and VLS at the P < .1 level. To examine the potential effect of the coronavirus disease 2019 (COVID-19) pandemic on results, multivariable models were reanalyzed, excluding all follow-up visits after March 2020.

To address selection biases induced by excluding participants with <2 consecutive visits on ART, inverse probability of selection weights were calculated using marginal probabilities of study inclusion, derived from logistic regression (Supplementary Table 2, Supplementary Figure 1) and included in multivariable analysis.

Ethics

The Johns Hopkins University School of Medicine Institutional Review Board and the Uganda Virus Research Institute Research and Ethics Committee approved the study protocol. Adults (aged ≥18 years) provided written informed consent prior to study procedures. Written assent and parental consent were obtained for unemancipated minors aged 15–17 years.

RESULTS

Overall, 33 219 individuals had ≥1 RCCS visit(s) from February 2015 to October 2020, of which 5814 (17.5%) were serologically confirmed to be diagnosed HIV (Figure 1). Of these, 913 (15.7%) persons without a follow-up visit during the observation period were excluded, and the remaining 3080 (62.8%) had ≥2 visits during the observation period. An additional 823 individuals were excluded due to ART nonuse (n = 257) or nonsequential study visits on ART (n = 566), yielding a final analytic sample of 2257 individuals contributing 3335 visit-pairs and 5959 person-years to the analysis.

Flow chart RCCS inclusion into the analytic sample. Abbreviations: ART, antiretroviral therapy; HIV, human immunodeficiency virus; RCCS, Rakai Community Cohort Study.
Figure 1.

Flow chart RCCS inclusion into the analytic sample. Abbreviations: ART, antiretroviral therapy; HIV, human immunodeficiency virus; RCCS, Rakai Community Cohort Study.

Table 1 presents descriptive participant characteristics at their first visit during the observation period (n = 2257), stratified by sex. The median age was 35 years (interquartile range [IQR], 30–40), and most participants were female (64.7%) and exhibited VLS (92.0%) at their first visit. Relative to the analytic cohort, individuals lost to follow-up (<2 consecutive visits) were significantly (P < .001) younger (<30 years, 36.0% vs 23.2%) and more likely to recently in-migrate into the study area (31.2% vs 14.6%) and exhibit HIV viremia (14.5% vs 8.0%) at their initial visit (Supplementary Table 2).

Table 1.

Participant-level Sample Characteristics at the First Study Visit in the Observation Period by Sex

CharacteristicOverall
N = 2257
Male
n = 797 (35.3%)
Female
n = 1460 (64.7%)
χ²
P Value
Index survey visit (calendar period)<.001
 Feb. 2015–Sept. 20161599 (70.8)519 (65.1)1080 (74.0)
 Oct. 2016–May 2018658 (29.2)278 (34.9)380 (26.0)
Age, median (IQR), ya35 (30–40)36 (31–41)34 (29–39)<.001
Age group, y<.001
 15–29522 (23.1)125 (15.7)397 (27.2)
 30–391117 (49.5)417 (52.3)700 (47.9)
 40–49618 (27.4)255 (32.0)363 (24.9)
Current marital status<.001
 Never married137 (6.1)37 (4.6)100 (6.9)
 Currently married1369 (60.7)554 (69.5)815 (55.8)
 Previously married751 (33.3)206 (35.9)545 (37.3)
Educational attainment.093
 No formal education196 (8.7)69 (8.7)127 (8.7)
 Primary1648 (73.0)604 (75.8)1044 (71.5)
 Secondary356 (15.8)106 (13.3)250 (17.1)
 Technical/University57 (2.5)18 (2.3)39 (2.7)
Primary occupation<.001
 Agriculture or housework836 (37.0)167 (20.9)669 (45.8)
 Trading or shopkeeping408 (21.3)107 (13.4)373 (25.6)
 Bar work, waitressing, or sex work189 (8.4)2 (0.3)187 (12.8)
 Fishing-related occupation354 (15.7)349 (43.8)5 (0.3)
 Other398 (17.6)172 (21.6)226 (15.5)
Religion.007
 Catholic/Christian1995 (88.4)727 (91.2)1268 (86.9)
 Muslim244 (10.8)64 (8.0)180 (12.3)
 Other or none18 (0.8)6 (0.8)12 (0.8)
Household wealth<.001
 Lowest965 (42.8)376 (47.2)589 (40.3)
 Low-middle453 (20.1)164 (20.6)289 (19.8)
 High-middle511 (22.6)155 (19.4)356 (24.4)
 Highest324 (14.3)98 (12.3)226 (15.5)
 Missing4 (0.2)4 (0.5)n/a
Recent migration.001
 Long-term resident1935 (85.7)709 (89.0)1226 (84.0)
 In-migrant322 (14.3)88 (11.0)234 (16.0)
Community type<.001
 Agrarian687 (30.4)210 (26.4)477 (32.7)
 Trading477 (21.1)122 (15.3)355 (24.3)
 Fishing1093 (48.4)465 (58.3)628 (43.0)
Distance to health facility, median (IQR), kma24.5 (14.4–33.4)28.5 (14.3–33.5)23.6 (14.4–33.4).007
Viral load, median (IQR), log10 copies/mLa7.0 (5.0–9.1)6.8 (5.0–9.3)7.3 (5.1–9.0).466
Viral load status.072
 Unsuppressed (≥200 copies/mL)181 (8.0)75 (9.4)106 (7.3)
 Suppressed (<200 copies/mL)2076 (92.0)722 (90.6)1354 (92.7)
CharacteristicOverall
N = 2257
Male
n = 797 (35.3%)
Female
n = 1460 (64.7%)
χ²
P Value
Index survey visit (calendar period)<.001
 Feb. 2015–Sept. 20161599 (70.8)519 (65.1)1080 (74.0)
 Oct. 2016–May 2018658 (29.2)278 (34.9)380 (26.0)
Age, median (IQR), ya35 (30–40)36 (31–41)34 (29–39)<.001
Age group, y<.001
 15–29522 (23.1)125 (15.7)397 (27.2)
 30–391117 (49.5)417 (52.3)700 (47.9)
 40–49618 (27.4)255 (32.0)363 (24.9)
Current marital status<.001
 Never married137 (6.1)37 (4.6)100 (6.9)
 Currently married1369 (60.7)554 (69.5)815 (55.8)
 Previously married751 (33.3)206 (35.9)545 (37.3)
Educational attainment.093
 No formal education196 (8.7)69 (8.7)127 (8.7)
 Primary1648 (73.0)604 (75.8)1044 (71.5)
 Secondary356 (15.8)106 (13.3)250 (17.1)
 Technical/University57 (2.5)18 (2.3)39 (2.7)
Primary occupation<.001
 Agriculture or housework836 (37.0)167 (20.9)669 (45.8)
 Trading or shopkeeping408 (21.3)107 (13.4)373 (25.6)
 Bar work, waitressing, or sex work189 (8.4)2 (0.3)187 (12.8)
 Fishing-related occupation354 (15.7)349 (43.8)5 (0.3)
 Other398 (17.6)172 (21.6)226 (15.5)
Religion.007
 Catholic/Christian1995 (88.4)727 (91.2)1268 (86.9)
 Muslim244 (10.8)64 (8.0)180 (12.3)
 Other or none18 (0.8)6 (0.8)12 (0.8)
Household wealth<.001
 Lowest965 (42.8)376 (47.2)589 (40.3)
 Low-middle453 (20.1)164 (20.6)289 (19.8)
 High-middle511 (22.6)155 (19.4)356 (24.4)
 Highest324 (14.3)98 (12.3)226 (15.5)
 Missing4 (0.2)4 (0.5)n/a
Recent migration.001
 Long-term resident1935 (85.7)709 (89.0)1226 (84.0)
 In-migrant322 (14.3)88 (11.0)234 (16.0)
Community type<.001
 Agrarian687 (30.4)210 (26.4)477 (32.7)
 Trading477 (21.1)122 (15.3)355 (24.3)
 Fishing1093 (48.4)465 (58.3)628 (43.0)
Distance to health facility, median (IQR), kma24.5 (14.4–33.4)28.5 (14.3–33.5)23.6 (14.4–33.4).007
Viral load, median (IQR), log10 copies/mLa7.0 (5.0–9.1)6.8 (5.0–9.3)7.3 (5.1–9.0).466
Viral load status.072
 Unsuppressed (≥200 copies/mL)181 (8.0)75 (9.4)106 (7.3)
 Suppressed (<200 copies/mL)2076 (92.0)722 (90.6)1354 (92.7)

Abbreviation: IQR, interquartile range.

aMedians and IQRs were estimated and compared by gender using Wilcoxon rank sum tests. Continuous viral load was only estimated among participants with detectable viral loads (>50 human immunodeficiency virus RNA copies/mL). Bolded values represent P-values that are statistically significant at the P <.05 level.

Table 1.

Participant-level Sample Characteristics at the First Study Visit in the Observation Period by Sex

CharacteristicOverall
N = 2257
Male
n = 797 (35.3%)
Female
n = 1460 (64.7%)
χ²
P Value
Index survey visit (calendar period)<.001
 Feb. 2015–Sept. 20161599 (70.8)519 (65.1)1080 (74.0)
 Oct. 2016–May 2018658 (29.2)278 (34.9)380 (26.0)
Age, median (IQR), ya35 (30–40)36 (31–41)34 (29–39)<.001
Age group, y<.001
 15–29522 (23.1)125 (15.7)397 (27.2)
 30–391117 (49.5)417 (52.3)700 (47.9)
 40–49618 (27.4)255 (32.0)363 (24.9)
Current marital status<.001
 Never married137 (6.1)37 (4.6)100 (6.9)
 Currently married1369 (60.7)554 (69.5)815 (55.8)
 Previously married751 (33.3)206 (35.9)545 (37.3)
Educational attainment.093
 No formal education196 (8.7)69 (8.7)127 (8.7)
 Primary1648 (73.0)604 (75.8)1044 (71.5)
 Secondary356 (15.8)106 (13.3)250 (17.1)
 Technical/University57 (2.5)18 (2.3)39 (2.7)
Primary occupation<.001
 Agriculture or housework836 (37.0)167 (20.9)669 (45.8)
 Trading or shopkeeping408 (21.3)107 (13.4)373 (25.6)
 Bar work, waitressing, or sex work189 (8.4)2 (0.3)187 (12.8)
 Fishing-related occupation354 (15.7)349 (43.8)5 (0.3)
 Other398 (17.6)172 (21.6)226 (15.5)
Religion.007
 Catholic/Christian1995 (88.4)727 (91.2)1268 (86.9)
 Muslim244 (10.8)64 (8.0)180 (12.3)
 Other or none18 (0.8)6 (0.8)12 (0.8)
Household wealth<.001
 Lowest965 (42.8)376 (47.2)589 (40.3)
 Low-middle453 (20.1)164 (20.6)289 (19.8)
 High-middle511 (22.6)155 (19.4)356 (24.4)
 Highest324 (14.3)98 (12.3)226 (15.5)
 Missing4 (0.2)4 (0.5)n/a
Recent migration.001
 Long-term resident1935 (85.7)709 (89.0)1226 (84.0)
 In-migrant322 (14.3)88 (11.0)234 (16.0)
Community type<.001
 Agrarian687 (30.4)210 (26.4)477 (32.7)
 Trading477 (21.1)122 (15.3)355 (24.3)
 Fishing1093 (48.4)465 (58.3)628 (43.0)
Distance to health facility, median (IQR), kma24.5 (14.4–33.4)28.5 (14.3–33.5)23.6 (14.4–33.4).007
Viral load, median (IQR), log10 copies/mLa7.0 (5.0–9.1)6.8 (5.0–9.3)7.3 (5.1–9.0).466
Viral load status.072
 Unsuppressed (≥200 copies/mL)181 (8.0)75 (9.4)106 (7.3)
 Suppressed (<200 copies/mL)2076 (92.0)722 (90.6)1354 (92.7)
CharacteristicOverall
N = 2257
Male
n = 797 (35.3%)
Female
n = 1460 (64.7%)
χ²
P Value
Index survey visit (calendar period)<.001
 Feb. 2015–Sept. 20161599 (70.8)519 (65.1)1080 (74.0)
 Oct. 2016–May 2018658 (29.2)278 (34.9)380 (26.0)
Age, median (IQR), ya35 (30–40)36 (31–41)34 (29–39)<.001
Age group, y<.001
 15–29522 (23.1)125 (15.7)397 (27.2)
 30–391117 (49.5)417 (52.3)700 (47.9)
 40–49618 (27.4)255 (32.0)363 (24.9)
Current marital status<.001
 Never married137 (6.1)37 (4.6)100 (6.9)
 Currently married1369 (60.7)554 (69.5)815 (55.8)
 Previously married751 (33.3)206 (35.9)545 (37.3)
Educational attainment.093
 No formal education196 (8.7)69 (8.7)127 (8.7)
 Primary1648 (73.0)604 (75.8)1044 (71.5)
 Secondary356 (15.8)106 (13.3)250 (17.1)
 Technical/University57 (2.5)18 (2.3)39 (2.7)
Primary occupation<.001
 Agriculture or housework836 (37.0)167 (20.9)669 (45.8)
 Trading or shopkeeping408 (21.3)107 (13.4)373 (25.6)
 Bar work, waitressing, or sex work189 (8.4)2 (0.3)187 (12.8)
 Fishing-related occupation354 (15.7)349 (43.8)5 (0.3)
 Other398 (17.6)172 (21.6)226 (15.5)
Religion.007
 Catholic/Christian1995 (88.4)727 (91.2)1268 (86.9)
 Muslim244 (10.8)64 (8.0)180 (12.3)
 Other or none18 (0.8)6 (0.8)12 (0.8)
Household wealth<.001
 Lowest965 (42.8)376 (47.2)589 (40.3)
 Low-middle453 (20.1)164 (20.6)289 (19.8)
 High-middle511 (22.6)155 (19.4)356 (24.4)
 Highest324 (14.3)98 (12.3)226 (15.5)
 Missing4 (0.2)4 (0.5)n/a
Recent migration.001
 Long-term resident1935 (85.7)709 (89.0)1226 (84.0)
 In-migrant322 (14.3)88 (11.0)234 (16.0)
Community type<.001
 Agrarian687 (30.4)210 (26.4)477 (32.7)
 Trading477 (21.1)122 (15.3)355 (24.3)
 Fishing1093 (48.4)465 (58.3)628 (43.0)
Distance to health facility, median (IQR), kma24.5 (14.4–33.4)28.5 (14.3–33.5)23.6 (14.4–33.4).007
Viral load, median (IQR), log10 copies/mLa7.0 (5.0–9.1)6.8 (5.0–9.3)7.3 (5.1–9.0).466
Viral load status.072
 Unsuppressed (≥200 copies/mL)181 (8.0)75 (9.4)106 (7.3)
 Suppressed (<200 copies/mL)2076 (92.0)722 (90.6)1354 (92.7)

Abbreviation: IQR, interquartile range.

aMedians and IQRs were estimated and compared by gender using Wilcoxon rank sum tests. Continuous viral load was only estimated among participants with detectable viral loads (>50 human immunodeficiency virus RNA copies/mL). Bolded values represent P-values that are statistically significant at the P <.05 level.

Incidence and Predictors of Facility Switching

Overall, 275 facility switching events were identified in 254 participants (11.3%) over the observation period, yielding an incidence rate of 4.8 per 100 person-years (95% CI, 4.2–5.5; Supplementary Table 3). Figure 2 displays age-stratified incidence rates of facility switching by sex and community type. Facility switching rates were comparable among males and females. Across community types, facility switching incidence was highest among individuals aged 25–29 years (inland, 5.1 per 100 person-years; 95% CI, 3.2–8.1 and fishing, 10.4 per 100 person-years; 95% CI, 7.7–14.0); however, facility switching was higher in fishing than inland communities across age strata.

Age-specific incidence rates of facility switching (per 100 person-years) by sex and community type. Error bars represent 95% confidence intervals for incidence rates of facility switching. Person-years are defined as accrued person-time (in years) between study visits rather than cumulative person-time accrued during the entirety of the observation period.
Figure 2.

Age-specific incidence rates of facility switching (per 100 person-years) by sex and community type. Error bars represent 95% confidence intervals for incidence rates of facility switching. Person-years are defined as accrued person-time (in years) between study visits rather than cumulative person-time accrued during the entirety of the observation period.

Table 2 summarizes sociodemographic and clinical predictors of facility switching at the visit-pair level (n = 3335). Younger people (<30 years, 29.4% vs 19.1%), households in the poorest wealth quartile (56.0% vs 40.5%), fishing community residents (68.7% vs 44.6%), in-migrants (22.2% vs 9.2%), and individuals who moved to new communities within the study area (7.3% vs 2.0%) were significantly (P < .001) more likely to switch facilities. Facility switching was also significantly more common among individuals who exhibited viremia at their index study visit than in suppressed persons (20.4% vs 5.8%, P < .001). Among females, facility switching was unassociated with pregnancy at the index or follow-up visit (15.3% vs 13.5%, P = .491).

Table 2.

Sociodemographic and Human Immunodeficiency Virus–related Outcomes at the Visit-Pair Level by Facility Switching Status

CharacteristicOverall
N = 3335
Stable Attendance
n = 3060 (91.7%)
Facility Switching
n = 275 (8.3%)
χ²
P Value
Index survey visit (calendar period).018
 Feb. 2015–Sept. 20161599 (47.9)1486 (48.6)113 (41.1)
 Oct. 2016–May 20181736 (52.1)1574 (51.4)162 (58.9)
Age, median (IQR), ya36 (31–41)36 (31–41)35 (29–40).001
Age group, y<.001
 15–29665 (20.0)584 (19.1)81 (29.4)
 30–391625 (48.7)1506 (49.2)119 (43.3)
 40–491045 (31.3)970 (31.7)75 (27.3)
Gender.072
 Male1122 (33.6)1043 (34.1)79 (28.7)
 Female2213 (66.4)2017 (65.9)196 (71.3)
Current marital status.080
 Never married189 (5.7)173 (5.6)16 (5.8)
 Currently married2035 (61.0)1884 (61.6)151 (54.9)
 Previously married1111 (33.3)1003 (32.8)108 (39.3)
Educational attainment.349
 No formal education290 (8.7)259 (8.5)31 (11.3)
 Primary2436 (73.0)2236 (73.1)200 (72.7)
 Secondary527 (15.8)488 (15.9)39 (14.2)
 Technical/University82 (2.5)77 (2.5)5 (1.8)
Primary occupation.003
 Agriculture or housework1299 (38.9)1219 (39.8)80 (29.1)
 Trading or shopkeeping680 (20.4)611 (20.0)69 (25.1)
 Bar work, waitressing, or sex work278 (8.3)248 (8.1)30 (10.9)
 Fishing-related occupation479 (14.4)430 (14.1)49 (17.8)
 Other599 (18.0)552 (18.0)47 (17.1)
Religion.715
 Catholic/Christian2947 (88.4)2703 (88.3)244 (88.7)
 Muslim362 (10.8)332 (10.8)30 (10.9)
 Other or none26 (0.8)25 (0.9)1 (0.4)
Household wealth<.001
 Lowest1394 (41.8)1240 (40.5)154 (56.0)
 Low-middle639 (19.2)594 (19.4)45 (16.4)
 High-middle786 (23.6)740 (24.2)46 (16.7)
 Highest512 (15.3)483 (15.8)29 (10.5)
 Missing4 (0.1)3 (0.1)1 (0.4)
Recent migration<.001
 Long-term resident2993 (89.7)2779 (90.8)214 (77.8)
 In-migrant342 (10.3)281 (9.2)61 (22.2)
Community type<.001
 Agrarian1075 (32.2)1033 (33.7)42 (15.3)
 Trading707 (21.2)663 (21.7)44 (16.0)
 Fishing1553 (46.6)1364 (44.6)189 (68.7)
Intercommunity mobility<.001
 No movement within communities3254 (97.6)2999 (98.0)255 (92.7)
 Moved to a new community81 (2.4)61 (2.0)20 (7.3)
Current pregnancy (females only), n = 2213303 (13.7)273 (13.5)30 (15.3).491
Distance to health facility, median (IQR), kma24.4 (14.2–33.4)24.4 (14.3–33.4)21.3 (12.6–34.9).146
Viral load status at index visit<.001
 Unsuppressed (≥200 copies/mL)234 (7.0)178 (5.8)56 (20.4)
 Suppressed (<200 copies/mL)3101 (93.0)2882 (94.2)219 (79.6)
CharacteristicOverall
N = 3335
Stable Attendance
n = 3060 (91.7%)
Facility Switching
n = 275 (8.3%)
χ²
P Value
Index survey visit (calendar period).018
 Feb. 2015–Sept. 20161599 (47.9)1486 (48.6)113 (41.1)
 Oct. 2016–May 20181736 (52.1)1574 (51.4)162 (58.9)
Age, median (IQR), ya36 (31–41)36 (31–41)35 (29–40).001
Age group, y<.001
 15–29665 (20.0)584 (19.1)81 (29.4)
 30–391625 (48.7)1506 (49.2)119 (43.3)
 40–491045 (31.3)970 (31.7)75 (27.3)
Gender.072
 Male1122 (33.6)1043 (34.1)79 (28.7)
 Female2213 (66.4)2017 (65.9)196 (71.3)
Current marital status.080
 Never married189 (5.7)173 (5.6)16 (5.8)
 Currently married2035 (61.0)1884 (61.6)151 (54.9)
 Previously married1111 (33.3)1003 (32.8)108 (39.3)
Educational attainment.349
 No formal education290 (8.7)259 (8.5)31 (11.3)
 Primary2436 (73.0)2236 (73.1)200 (72.7)
 Secondary527 (15.8)488 (15.9)39 (14.2)
 Technical/University82 (2.5)77 (2.5)5 (1.8)
Primary occupation.003
 Agriculture or housework1299 (38.9)1219 (39.8)80 (29.1)
 Trading or shopkeeping680 (20.4)611 (20.0)69 (25.1)
 Bar work, waitressing, or sex work278 (8.3)248 (8.1)30 (10.9)
 Fishing-related occupation479 (14.4)430 (14.1)49 (17.8)
 Other599 (18.0)552 (18.0)47 (17.1)
Religion.715
 Catholic/Christian2947 (88.4)2703 (88.3)244 (88.7)
 Muslim362 (10.8)332 (10.8)30 (10.9)
 Other or none26 (0.8)25 (0.9)1 (0.4)
Household wealth<.001
 Lowest1394 (41.8)1240 (40.5)154 (56.0)
 Low-middle639 (19.2)594 (19.4)45 (16.4)
 High-middle786 (23.6)740 (24.2)46 (16.7)
 Highest512 (15.3)483 (15.8)29 (10.5)
 Missing4 (0.1)3 (0.1)1 (0.4)
Recent migration<.001
 Long-term resident2993 (89.7)2779 (90.8)214 (77.8)
 In-migrant342 (10.3)281 (9.2)61 (22.2)
Community type<.001
 Agrarian1075 (32.2)1033 (33.7)42 (15.3)
 Trading707 (21.2)663 (21.7)44 (16.0)
 Fishing1553 (46.6)1364 (44.6)189 (68.7)
Intercommunity mobility<.001
 No movement within communities3254 (97.6)2999 (98.0)255 (92.7)
 Moved to a new community81 (2.4)61 (2.0)20 (7.3)
Current pregnancy (females only), n = 2213303 (13.7)273 (13.5)30 (15.3).491
Distance to health facility, median (IQR), kma24.4 (14.2–33.4)24.4 (14.3–33.4)21.3 (12.6–34.9).146
Viral load status at index visit<.001
 Unsuppressed (≥200 copies/mL)234 (7.0)178 (5.8)56 (20.4)
 Suppressed (<200 copies/mL)3101 (93.0)2882 (94.2)219 (79.6)

Abbreviation: IQR, interquartile range.

aMedians and IQRs were estimated and compared by facility switching status using Wilcoxon rank-sum tests.

Bolded values represent P-values that are statistically significant at the P <.05 level.

Table 2.

Sociodemographic and Human Immunodeficiency Virus–related Outcomes at the Visit-Pair Level by Facility Switching Status

CharacteristicOverall
N = 3335
Stable Attendance
n = 3060 (91.7%)
Facility Switching
n = 275 (8.3%)
χ²
P Value
Index survey visit (calendar period).018
 Feb. 2015–Sept. 20161599 (47.9)1486 (48.6)113 (41.1)
 Oct. 2016–May 20181736 (52.1)1574 (51.4)162 (58.9)
Age, median (IQR), ya36 (31–41)36 (31–41)35 (29–40).001
Age group, y<.001
 15–29665 (20.0)584 (19.1)81 (29.4)
 30–391625 (48.7)1506 (49.2)119 (43.3)
 40–491045 (31.3)970 (31.7)75 (27.3)
Gender.072
 Male1122 (33.6)1043 (34.1)79 (28.7)
 Female2213 (66.4)2017 (65.9)196 (71.3)
Current marital status.080
 Never married189 (5.7)173 (5.6)16 (5.8)
 Currently married2035 (61.0)1884 (61.6)151 (54.9)
 Previously married1111 (33.3)1003 (32.8)108 (39.3)
Educational attainment.349
 No formal education290 (8.7)259 (8.5)31 (11.3)
 Primary2436 (73.0)2236 (73.1)200 (72.7)
 Secondary527 (15.8)488 (15.9)39 (14.2)
 Technical/University82 (2.5)77 (2.5)5 (1.8)
Primary occupation.003
 Agriculture or housework1299 (38.9)1219 (39.8)80 (29.1)
 Trading or shopkeeping680 (20.4)611 (20.0)69 (25.1)
 Bar work, waitressing, or sex work278 (8.3)248 (8.1)30 (10.9)
 Fishing-related occupation479 (14.4)430 (14.1)49 (17.8)
 Other599 (18.0)552 (18.0)47 (17.1)
Religion.715
 Catholic/Christian2947 (88.4)2703 (88.3)244 (88.7)
 Muslim362 (10.8)332 (10.8)30 (10.9)
 Other or none26 (0.8)25 (0.9)1 (0.4)
Household wealth<.001
 Lowest1394 (41.8)1240 (40.5)154 (56.0)
 Low-middle639 (19.2)594 (19.4)45 (16.4)
 High-middle786 (23.6)740 (24.2)46 (16.7)
 Highest512 (15.3)483 (15.8)29 (10.5)
 Missing4 (0.1)3 (0.1)1 (0.4)
Recent migration<.001
 Long-term resident2993 (89.7)2779 (90.8)214 (77.8)
 In-migrant342 (10.3)281 (9.2)61 (22.2)
Community type<.001
 Agrarian1075 (32.2)1033 (33.7)42 (15.3)
 Trading707 (21.2)663 (21.7)44 (16.0)
 Fishing1553 (46.6)1364 (44.6)189 (68.7)
Intercommunity mobility<.001
 No movement within communities3254 (97.6)2999 (98.0)255 (92.7)
 Moved to a new community81 (2.4)61 (2.0)20 (7.3)
Current pregnancy (females only), n = 2213303 (13.7)273 (13.5)30 (15.3).491
Distance to health facility, median (IQR), kma24.4 (14.2–33.4)24.4 (14.3–33.4)21.3 (12.6–34.9).146
Viral load status at index visit<.001
 Unsuppressed (≥200 copies/mL)234 (7.0)178 (5.8)56 (20.4)
 Suppressed (<200 copies/mL)3101 (93.0)2882 (94.2)219 (79.6)
CharacteristicOverall
N = 3335
Stable Attendance
n = 3060 (91.7%)
Facility Switching
n = 275 (8.3%)
χ²
P Value
Index survey visit (calendar period).018
 Feb. 2015–Sept. 20161599 (47.9)1486 (48.6)113 (41.1)
 Oct. 2016–May 20181736 (52.1)1574 (51.4)162 (58.9)
Age, median (IQR), ya36 (31–41)36 (31–41)35 (29–40).001
Age group, y<.001
 15–29665 (20.0)584 (19.1)81 (29.4)
 30–391625 (48.7)1506 (49.2)119 (43.3)
 40–491045 (31.3)970 (31.7)75 (27.3)
Gender.072
 Male1122 (33.6)1043 (34.1)79 (28.7)
 Female2213 (66.4)2017 (65.9)196 (71.3)
Current marital status.080
 Never married189 (5.7)173 (5.6)16 (5.8)
 Currently married2035 (61.0)1884 (61.6)151 (54.9)
 Previously married1111 (33.3)1003 (32.8)108 (39.3)
Educational attainment.349
 No formal education290 (8.7)259 (8.5)31 (11.3)
 Primary2436 (73.0)2236 (73.1)200 (72.7)
 Secondary527 (15.8)488 (15.9)39 (14.2)
 Technical/University82 (2.5)77 (2.5)5 (1.8)
Primary occupation.003
 Agriculture or housework1299 (38.9)1219 (39.8)80 (29.1)
 Trading or shopkeeping680 (20.4)611 (20.0)69 (25.1)
 Bar work, waitressing, or sex work278 (8.3)248 (8.1)30 (10.9)
 Fishing-related occupation479 (14.4)430 (14.1)49 (17.8)
 Other599 (18.0)552 (18.0)47 (17.1)
Religion.715
 Catholic/Christian2947 (88.4)2703 (88.3)244 (88.7)
 Muslim362 (10.8)332 (10.8)30 (10.9)
 Other or none26 (0.8)25 (0.9)1 (0.4)
Household wealth<.001
 Lowest1394 (41.8)1240 (40.5)154 (56.0)
 Low-middle639 (19.2)594 (19.4)45 (16.4)
 High-middle786 (23.6)740 (24.2)46 (16.7)
 Highest512 (15.3)483 (15.8)29 (10.5)
 Missing4 (0.1)3 (0.1)1 (0.4)
Recent migration<.001
 Long-term resident2993 (89.7)2779 (90.8)214 (77.8)
 In-migrant342 (10.3)281 (9.2)61 (22.2)
Community type<.001
 Agrarian1075 (32.2)1033 (33.7)42 (15.3)
 Trading707 (21.2)663 (21.7)44 (16.0)
 Fishing1553 (46.6)1364 (44.6)189 (68.7)
Intercommunity mobility<.001
 No movement within communities3254 (97.6)2999 (98.0)255 (92.7)
 Moved to a new community81 (2.4)61 (2.0)20 (7.3)
Current pregnancy (females only), n = 2213303 (13.7)273 (13.5)30 (15.3).491
Distance to health facility, median (IQR), kma24.4 (14.2–33.4)24.4 (14.3–33.4)21.3 (12.6–34.9).146
Viral load status at index visit<.001
 Unsuppressed (≥200 copies/mL)234 (7.0)178 (5.8)56 (20.4)
 Suppressed (<200 copies/mL)3101 (93.0)2882 (94.2)219 (79.6)

Abbreviation: IQR, interquartile range.

aMedians and IQRs were estimated and compared by facility switching status using Wilcoxon rank-sum tests.

Bolded values represent P-values that are statistically significant at the P <.05 level.

Distance From Households to Health Facilities

At index visits, the median distance from households to health facilities was 24.5 km (IQR, 14.4–33.4; Table 1). Distance from households to health facilities, however, was unassociated with subsequent facility switching (21.3 vs 24.4 km; P = .146; Table 2). The largest fraction of facility switches were within 5 km of participants’ households (37.3%), followed by switches to facilities located ≥5 km from their households (33.2%; Supplementary Figure 2). The median travel distance difference (d1d2) between origin and destination clinics was 3.8 km (IQR, −11.5  to 9.5), indicating that persons tended to switch to facilities closer to their homes.

Facility Switching, New HIV Suppression, and Viral Rebound

Most visit-pairs exhibited durable VLS (90.5%) during the observation period (Supplementary Table 4); however, durable VLS was substantially less common in visit-pairs where facility switching was observed (75.1% vs 91.9%).

Table 3 presents results from multivariable analysis of facility switching and new suppression and viral rebound, respectively. While new VLS was unassociated with facility switching in the overall sample (adjIRR = 1.17; 95% CI, .76–1.78), individuals in inland communities who exhibited viremia initially were twice as likely to achieve VLS at follow-up if they switched facilities (adjIRR = 2.15; 95% CI, 1.27–3.62).

Table 3.

Incidence Rate Ratios of New/Renewed Suppression and Viral Rebound by Facility Switching Status Stratified by Sex and Community Type

CharacteristicNo. of EventsPerson-YearsIR Per 100
Person-Years (95% CI)
IRR (95% CI)adjIRR (95% CI)IPW adjIRR
(95% CI)
New suppression (ref. persistent viremia)
 Overall
  Stable attendance8124333.3 (26.8–41.3)Ref.Ref.Ref.
  Facility switching297439.1 (27.6–55.5)1.18 (.78–1.78)1.17 (.76–1.78)1.16 (.75–1.79)
 Male
  Stable attendance369637.3 (27.1–51.3)Ref.Ref.Ref.
  Facility switching113234.3 (18.9–62.0)0.92 (.47–1.78)0.89 (.47–1.69)0.89 (.50–1.66)
 Female
  Stable attendance4514730.6 (22.8–41.1)Ref.Ref.Ref.
  Facility switching184242.9 (27.8–66.2)1.40 (.83–2.37)1.45 (.87–2.44)1.40 (.79–2.45)
 Inland communities
  Stable attendance4013330.0 (21.9–41.0)Ref.Ref.Ref.
  Facility switching121770.6 (45.7–95.5)2.35 (1.38–4.02)2.15 (1.27–3.62)1.94 (1.10–3.42)
 Fishing communities
  Stable attendance4110037.2 (27.6–50.2)Ref.Ref.Ref.
  Facility switching175729.8 (18.7–47.6)0.80 (.46–1.38)0.72 (.42–1.27)0.75 (.43–1.29)
Viral rebound (ref. durable viral load suppression)
 Overall
  Stable attendance7150781.4 (1.1–1.8)Ref.Ref.Ref.
  Facility switching123923.1 (1.7–5.4)2.19 (1.19–4.05)2.27 (1.16–4.45)2.05 (1.02–4.13)
 Male
  Stable attendance3716852.2 (1.6–3.0)Ref.Ref.Ref.
  Facility switching5995.0 (2.1–12.2)2.29 (.90–5.85)2.02 (.71–5.74)1.83 (.62–5.41)
 Female
  Stable attendance3433921.0 (.7–1.4)Ref.Ref.Ref.
  Facility switching72922.4 (1.1–5.1)2.39 (1.06–5.41)2.40 (.98–5.89)2.13 (.83–5.47)
 Inland communities
  Stable attendance3528671.2 (.9–1.7)Ref.Ref.Ref.
  Facility switching61274.7 (2.1–10.7)3.88 (1.62–9.28)4.35 (1.76–10.73)3.98 (1.57–10.07)
 Fishing communities
  Stable attendance3622101.6 (1.2–2.3)Ref.Ref.Ref.
  Facility switching62652.3 (1.0–5.1)1.39 (.58–3.31)1.31 (.54–3.15)1.19 (.47–3.01)
CharacteristicNo. of EventsPerson-YearsIR Per 100
Person-Years (95% CI)
IRR (95% CI)adjIRR (95% CI)IPW adjIRR
(95% CI)
New suppression (ref. persistent viremia)
 Overall
  Stable attendance8124333.3 (26.8–41.3)Ref.Ref.Ref.
  Facility switching297439.1 (27.6–55.5)1.18 (.78–1.78)1.17 (.76–1.78)1.16 (.75–1.79)
 Male
  Stable attendance369637.3 (27.1–51.3)Ref.Ref.Ref.
  Facility switching113234.3 (18.9–62.0)0.92 (.47–1.78)0.89 (.47–1.69)0.89 (.50–1.66)
 Female
  Stable attendance4514730.6 (22.8–41.1)Ref.Ref.Ref.
  Facility switching184242.9 (27.8–66.2)1.40 (.83–2.37)1.45 (.87–2.44)1.40 (.79–2.45)
 Inland communities
  Stable attendance4013330.0 (21.9–41.0)Ref.Ref.Ref.
  Facility switching121770.6 (45.7–95.5)2.35 (1.38–4.02)2.15 (1.27–3.62)1.94 (1.10–3.42)
 Fishing communities
  Stable attendance4110037.2 (27.6–50.2)Ref.Ref.Ref.
  Facility switching175729.8 (18.7–47.6)0.80 (.46–1.38)0.72 (.42–1.27)0.75 (.43–1.29)
Viral rebound (ref. durable viral load suppression)
 Overall
  Stable attendance7150781.4 (1.1–1.8)Ref.Ref.Ref.
  Facility switching123923.1 (1.7–5.4)2.19 (1.19–4.05)2.27 (1.16–4.45)2.05 (1.02–4.13)
 Male
  Stable attendance3716852.2 (1.6–3.0)Ref.Ref.Ref.
  Facility switching5995.0 (2.1–12.2)2.29 (.90–5.85)2.02 (.71–5.74)1.83 (.62–5.41)
 Female
  Stable attendance3433921.0 (.7–1.4)Ref.Ref.Ref.
  Facility switching72922.4 (1.1–5.1)2.39 (1.06–5.41)2.40 (.98–5.89)2.13 (.83–5.47)
 Inland communities
  Stable attendance3528671.2 (.9–1.7)Ref.Ref.Ref.
  Facility switching61274.7 (2.1–10.7)3.88 (1.62–9.28)4.35 (1.76–10.73)3.98 (1.57–10.07)
 Fishing communities
  Stable attendance3622101.6 (1.2–2.3)Ref.Ref.Ref.
  Facility switching62652.3 (1.0–5.1)1.39 (.58–3.31)1.31 (.54–3.15)1.19 (.47–3.01)

Abbreviations: adjIRR, adjusted incidence rate ratio; CI, confidence interval; IPW, inverse probability of selection weights; IR, incidence rate; IRR, incidence rate ratio; Ref., reference.

Bolded values represent p-values that are statistically significant at the P < .05 level. Person-years are defined as accrued person-time (in years) between study visits rather than cumulative person-time accrued during the entirety of the observation period. IRRs were generated from Poisson regression with generalized estimating equations, independent correlation structures, and robust standard errors. Visit-pairs that exhibited persistent viremia were the reference category for new suppression. Durable viral load suppression was the reference category for viral rebound. Multivariable models adjusted for index survey visit, age, sex, current marital status, primary occupation, household wealth, community type, and migration history. Stabilized IPWs were implemented in sensitivity analyses to correct for potential selection biases induced by exclusion of individuals with fewer than 2 consecutive visits during the observation period in the analytic cohort.

Table 3.

Incidence Rate Ratios of New/Renewed Suppression and Viral Rebound by Facility Switching Status Stratified by Sex and Community Type

CharacteristicNo. of EventsPerson-YearsIR Per 100
Person-Years (95% CI)
IRR (95% CI)adjIRR (95% CI)IPW adjIRR
(95% CI)
New suppression (ref. persistent viremia)
 Overall
  Stable attendance8124333.3 (26.8–41.3)Ref.Ref.Ref.
  Facility switching297439.1 (27.6–55.5)1.18 (.78–1.78)1.17 (.76–1.78)1.16 (.75–1.79)
 Male
  Stable attendance369637.3 (27.1–51.3)Ref.Ref.Ref.
  Facility switching113234.3 (18.9–62.0)0.92 (.47–1.78)0.89 (.47–1.69)0.89 (.50–1.66)
 Female
  Stable attendance4514730.6 (22.8–41.1)Ref.Ref.Ref.
  Facility switching184242.9 (27.8–66.2)1.40 (.83–2.37)1.45 (.87–2.44)1.40 (.79–2.45)
 Inland communities
  Stable attendance4013330.0 (21.9–41.0)Ref.Ref.Ref.
  Facility switching121770.6 (45.7–95.5)2.35 (1.38–4.02)2.15 (1.27–3.62)1.94 (1.10–3.42)
 Fishing communities
  Stable attendance4110037.2 (27.6–50.2)Ref.Ref.Ref.
  Facility switching175729.8 (18.7–47.6)0.80 (.46–1.38)0.72 (.42–1.27)0.75 (.43–1.29)
Viral rebound (ref. durable viral load suppression)
 Overall
  Stable attendance7150781.4 (1.1–1.8)Ref.Ref.Ref.
  Facility switching123923.1 (1.7–5.4)2.19 (1.19–4.05)2.27 (1.16–4.45)2.05 (1.02–4.13)
 Male
  Stable attendance3716852.2 (1.6–3.0)Ref.Ref.Ref.
  Facility switching5995.0 (2.1–12.2)2.29 (.90–5.85)2.02 (.71–5.74)1.83 (.62–5.41)
 Female
  Stable attendance3433921.0 (.7–1.4)Ref.Ref.Ref.
  Facility switching72922.4 (1.1–5.1)2.39 (1.06–5.41)2.40 (.98–5.89)2.13 (.83–5.47)
 Inland communities
  Stable attendance3528671.2 (.9–1.7)Ref.Ref.Ref.
  Facility switching61274.7 (2.1–10.7)3.88 (1.62–9.28)4.35 (1.76–10.73)3.98 (1.57–10.07)
 Fishing communities
  Stable attendance3622101.6 (1.2–2.3)Ref.Ref.Ref.
  Facility switching62652.3 (1.0–5.1)1.39 (.58–3.31)1.31 (.54–3.15)1.19 (.47–3.01)
CharacteristicNo. of EventsPerson-YearsIR Per 100
Person-Years (95% CI)
IRR (95% CI)adjIRR (95% CI)IPW adjIRR
(95% CI)
New suppression (ref. persistent viremia)
 Overall
  Stable attendance8124333.3 (26.8–41.3)Ref.Ref.Ref.
  Facility switching297439.1 (27.6–55.5)1.18 (.78–1.78)1.17 (.76–1.78)1.16 (.75–1.79)
 Male
  Stable attendance369637.3 (27.1–51.3)Ref.Ref.Ref.
  Facility switching113234.3 (18.9–62.0)0.92 (.47–1.78)0.89 (.47–1.69)0.89 (.50–1.66)
 Female
  Stable attendance4514730.6 (22.8–41.1)Ref.Ref.Ref.
  Facility switching184242.9 (27.8–66.2)1.40 (.83–2.37)1.45 (.87–2.44)1.40 (.79–2.45)
 Inland communities
  Stable attendance4013330.0 (21.9–41.0)Ref.Ref.Ref.
  Facility switching121770.6 (45.7–95.5)2.35 (1.38–4.02)2.15 (1.27–3.62)1.94 (1.10–3.42)
 Fishing communities
  Stable attendance4110037.2 (27.6–50.2)Ref.Ref.Ref.
  Facility switching175729.8 (18.7–47.6)0.80 (.46–1.38)0.72 (.42–1.27)0.75 (.43–1.29)
Viral rebound (ref. durable viral load suppression)
 Overall
  Stable attendance7150781.4 (1.1–1.8)Ref.Ref.Ref.
  Facility switching123923.1 (1.7–5.4)2.19 (1.19–4.05)2.27 (1.16–4.45)2.05 (1.02–4.13)
 Male
  Stable attendance3716852.2 (1.6–3.0)Ref.Ref.Ref.
  Facility switching5995.0 (2.1–12.2)2.29 (.90–5.85)2.02 (.71–5.74)1.83 (.62–5.41)
 Female
  Stable attendance3433921.0 (.7–1.4)Ref.Ref.Ref.
  Facility switching72922.4 (1.1–5.1)2.39 (1.06–5.41)2.40 (.98–5.89)2.13 (.83–5.47)
 Inland communities
  Stable attendance3528671.2 (.9–1.7)Ref.Ref.Ref.
  Facility switching61274.7 (2.1–10.7)3.88 (1.62–9.28)4.35 (1.76–10.73)3.98 (1.57–10.07)
 Fishing communities
  Stable attendance3622101.6 (1.2–2.3)Ref.Ref.Ref.
  Facility switching62652.3 (1.0–5.1)1.39 (.58–3.31)1.31 (.54–3.15)1.19 (.47–3.01)

Abbreviations: adjIRR, adjusted incidence rate ratio; CI, confidence interval; IPW, inverse probability of selection weights; IR, incidence rate; IRR, incidence rate ratio; Ref., reference.

Bolded values represent p-values that are statistically significant at the P < .05 level. Person-years are defined as accrued person-time (in years) between study visits rather than cumulative person-time accrued during the entirety of the observation period. IRRs were generated from Poisson regression with generalized estimating equations, independent correlation structures, and robust standard errors. Visit-pairs that exhibited persistent viremia were the reference category for new suppression. Durable viral load suppression was the reference category for viral rebound. Multivariable models adjusted for index survey visit, age, sex, current marital status, primary occupation, household wealth, community type, and migration history. Stabilized IPWs were implemented in sensitivity analyses to correct for potential selection biases induced by exclusion of individuals with fewer than 2 consecutive visits during the observation period in the analytic cohort.

Conversely, individuals who exhibited VLS at their index visit were twice as likely to exhibit viremia at follow-up if they switched facilities (adjIRR = 2.27; 95% CI, 1.16–4.45). A similar, albeit more pronounced, effect size was observed among inland community residents (adjIRR = 4.35; 95% CI, 1.76–10.73). Findings were similar in weighted analyses (Table 3) and sensitivity analyses restricted to visits through March 2020 (Supplementary Table 5), indicating limited effects of the COVID-19 pandemic restrictions on study findings. Travel distance differences from households to health facilities were also unassociated with virologic outcomes in visit-pairs where facility switching was observed (Supplementary Figure 3).

DISCUSSION

In this population-based study, facility switching was common and significantly associated with HIV viremia, specifically among individuals who exhibited VLS prior to facility switching. We found that facility switching was associated with a 2-fold rate increase in rebounding viremia among initially suppressed persons, reaffirming how facility switching can destabilize HIV care continuity [11–13]. Loss of transfer paperwork, delays in processing transfers, and denial of services to “silently transferring” clients are documented mechanisms through which HIV treatment interruptions can occur during facility switching [8, 12]. Furthermore, the relationship between facility switching and viral rebound was most pronounced in inland communities, where facility switching is less common and, thus, a potentially more disruptive event than in fishing communities. Fishing community residents, whose mobility within and outside the study area can be intermittent and episodic [28, 29], may be more experienced than inland community residents at navigating the HIV care system, particularly during periods of mobility. Addressing systemic bottlenecks to issuing and administering client transfers, including retroactive processing of clients who self-transferred to facilities, is imperative to maintaining care continuity for persons who switch clinics.

Facility switching prevalence (approximately 8%) and incidence (approximately 5 per 100 person-years) in this study population were comparable to estimates from other studies in high-prevalence settings [6, 7, 15, 30]. Consistent with other studies in sub-Saharan Africa [8, 9, 31], mobility was a critical driver of facility switching in the present study. Facility switching was significantly concentrated in highly mobile groups (ie, young people, in-migrants, fishing community residents) who already shoulder a disproportionate burden of HIV viremia in the overall population [32]. Likewise, durable VLS was substantially less common in the context of facility switching (approximately 75%), revealing suboptimal virologic outcomes among persons on ART with facility switching histories. Expanding access to multimonth ART dispensing for viremic persons (who are ineligible for extended prescription durations in most settings) could mitigate care retention challenges for individuals who switched facilities.

Among individuals who switched health facilities, there was substantial heterogeneity in household proximity to health facilities after switching. Although persons switched predominantly to facilities closer to their households, comparable fractions of switching occurred when individuals transferred their care to facilities further from their households, suggesting facility switching may be motivated by factors beyond travel distance. Importantly, distances from households to health facilities were unassociated with facility switching, reaffirming that individuals may prefer receiving HIV services from clinics located further from their homes [25, 33, 34]. The proliferation of differentiated service delivery models (DSDMs) has also expanded treatment access by reducing travel time and clinic visit frequency [35]. Sustained DSDM scale-up could continue rendering travel distances to clinics negligible, minimizing the burden of facility transfers for clients and health systems alike.

This study has some key limitations. First, identified associations of facility switching with viral rebound remain temporally ambiguous. Viral rebound is a plausible consequence of facility switching, but facility switching may also be a response to rebounding viremia. Poor-quality care and mistreatment by providers, for instance, are known predictors of poor virologic outcomes motivating facility switching [36–38]. Future studies should seek to temporally disentangle the relationship between VLS and facility switching. Second, only RCCS participants observed over ≥2 rounds in the cohort were eligible for inclusion; thus, highly mobile persons with greater propensity for facility switching were likely excluded, inducing selection bias, though this was potentially attenuated through inverse probability weighting. Third, facility switching and viral load status were ascertained only at the start and end of visit-pairs, likely underestimating switching events or changes in viremia that occurred during the interval. Likewise, pregnancy ascertainment between approximately 18-month study intervals likely underestimates pregnancies in the study population and their potential relationship with facility switching. Pregnant persons may be diagnosed with HIV during antenatal care and initiate ART at one facility, then switch to another facility postpartum, which has been documented in other settings such as South Africa [6, 39]. Fourth, facility switching was ascertained by self-report and is, thus, susceptible to response and recall biases. Fifth, despite widespread DSDM implementation in Uganda [40], this study did not measure the proportion of persons enrolled in community-based DSDMs who may not necessarily need to visit their designated health facilities for ART refills. Sixth, unlike other clinic-based studies, this study could not discern whether observed facility switching events were formal transfers or self-transfers, which could further clarify the relationship between facility switching and viremia. Last, unmeasured clinical factors likely associated with virologic outcomes (eg, ART use duration, ART regimen, HIV drug resistance) were unavailable.

CONCLUSIONS

In this population-based study of persons on ART in southern Uganda, the frequency of facility switching over a 5-year period was high and significantly associated with loss of VLS. Although facility switching may benefit persons on ART when desired by clients themselves, it can place strain on underresourced health systems with limited capacity to surveil client movements between clinics. Investments in more agile, person-centered models for mobile and complex clients, including retroactive transfer processing and expedited reenrollments for self-transfers, could address system inefficiencies and bottlenecks that disrupt HIV care continuity.

Supplementary Data

Supplementary materials are available at Clinical Infectious Diseases online. Consisting of data provided by the authors to benefit the reader, the posted materials are not copyedited and are the sole responsibility of the authors, so questions or comments should be addressed to the corresponding author.

Notes

Acknowledgments. The authors thank all Rakai Community Cohort Study (RCCS) participants for their time and contributions to the present study. They also acknowledge the dedication and efforts of all RCCS data collection, laboratory, service linkage, and data management staff.

Disclaimer. The contents presented here are solely the responsibility of the authors and do not necessarily reflect the views or policies of the funding agencies, nor do mentions of trade names, commercial products, or organizations imply endorsement by the US government.

Financial support. This work was supported by the National Institute of Allergy and Infectious Diseases (U01AI100031, U01AI075115, R01AI110324, R01AI128779, R01AI123002, R01AI143333, R01AI155080, T32AI102623), the National Institute of Mental Health (R01MH105313), the National Cancer Institute (75N91019D00024), and the US President’s Emergency Plan for AIDS Relief through the Centers for Disease Control and Prevention under the terms of NU2GGH00081. J. G. R. and K. B. R. were supported by the National Institute of Mental Health (F31MH126796, K01MH129226).

References

1

Grabowski
 
MK
,
Serwadda
 
DM
,
Gray
 
RH
, et al.  
Combination HIV prevention and HIV incidence in Uganda
.
N Engl J Med
 
2017
;
377
:
2154
66
.

2

Kagaayi
 
J
,
Chang
 
LW
,
Ssempijja
 
V
, et al.  
Impact of combination HIV interventions on HIV incidence in hyperendemic fishing communities in Uganda: a prospective cohort study
.
Lancet HIV
 
2019
;
6
:
e680
7
.

3

Uganda Ministry of Health
. Uganda Population-Based HIV Impact Assessment (UPHIA), 2016–2017: Final Report. ICAP; 2019. Available at: https://phia.icap.columbia.edu/wp-content/uploads/2019/07/UPHIA_Final_Report_Revise_07.11.2019_Final_for-web.pdf. Accessed 16 December 2022.

4

Kebede
 
HK
,
Mwanri
 
L
,
Ward
 
P
,
Gesesew
 
HA
.
Predictors of lost to follow up from antiretroviral therapy among adults in sub-Saharan Africa: a systematic review and meta-analysis
.
Infect Dis Poverty
 
2021
;
10
:
33
.

5

Wilkinson
 
LS
,
Skordis-Worrall
 
J
,
Ajose
 
O
,
Ford
 
N
.
Self-transfer and mortality amongst adults lost to follow-up in ART programmes in low- and middle-income countries: systematic review and meta-analysis
.
Trop Med Int Health
 
2015
;
20
:
365
79
.

6

Clouse
 
K
,
Vermund
 
SH
,
Maskew
 
M
, et al.  
Mobility and clinic switching among postpartum women considered lost to HIV care in South Africa
.
J Acquir Immune Defic Syndr
 
2017
;
74
:
383
9
.

7

Etoori
 
D
,
Wringe
 
A
,
Kabudula
 
CW
, et al.  
Misreporting of patient outcomes in the South African national HIV treatment database: consequences for programme planning, monitoring, and evaluation
.
Front Public Health
 
2020
;
8
:
100
.

8

Bengtson
 
AM
,
Go
 
V
,
Kumwenda
 
W
, et al.  
“A way of escaping”: a qualitative study exploring reasons for clinic transferring and its impact on engagement in care among women in option B+
.
AIDS Care
 
2020
;
32
:
72
5
.

9

Sikazwe
 
I
,
Eshun-Wilson
 
I
,
Sikombe
 
K
, et al.  
Patient-reported reasons for stopping care or switching clinics in Zambia: a multisite, regionally representative estimate using a multistage sampling-based approach in Zambia
.
Clin Infect Dis
 
2021
;
73
:
e2294
302
.

10

Beeman
 
A
,
Bengtson
 
AM
,
Swartz
 
A
,
Colvin
 
CJ
,
Lurie
 
MN
.
Cyclical engagement in HIV care: a qualitative study of clinic transfers to re-enter HIV care in Cape Town, South Africa
.
AIDS Behav
 
2022
;
26
:
2387
96
.

11

Teasdale
 
CA
,
Sogaula
 
N
,
Yuengling
 
KA
, et al.  
High risk of loss to follow-up among South African children on ART during transfer, a retrospective cohort analysis with community tracing
.
J Int AIDS Soc
 
2017
;
20
:
21748
.

12

Sikombe
 
K
,
Mody
 
A
,
Kadota
 
J
, et al.  
Understanding patient transfers across multiple clinics in Zambia among HIV infected adults
.
PLoS One
 
2020
;
15
:
e0241477
.

13

Bengtson
 
AM
,
Espinosa Dice
 
AL
,
Kirwa
 
K
,
Cornell
 
M
,
Colvin
 
CJ
,
Lurie
 
MN
.
Patient transfers and their impact on gaps in clinical care: differences by gender in a large cohort of adults living with HIV on antiretroviral therapy in South Africa
.
AIDS Behav
 
2021
;
25
:
3337
46
.

14

Yu
 
JKL
,
Tok
 
TS
,
Tsai
 
JJ
, et al.  
What happens to patients on antiretroviral therapy who transfer out to another facility?
 
PLoS One
 
2008
;
3
:
e2065
.

15

Nglazi
 
MD
,
Kaplan
 
R
,
Orrell
 
C
, et al.  
Increasing transfers-out from an antiretroviral treatment service in South Africa: patient characteristics and rates of virological non-suppression
.
PLoS One
 
2013
;
8
:
e57907
.

16

Odayar
 
J
,
Chi
 
BH
,
Phillips
 
TK
, et al.  
Transfer of patients on antiretroviral therapy attending primary health care services in South Africa
.
J Acquir Immune Defic Syndr
 
2022
;
90
:
309
15
.

17

Chang
 
LW
,
Grabowski
 
MK
,
Ssekubugu
 
R
, et al.  
Heterogeneity of the HIV epidemic in agrarian, trading, and fishing communities in Rakai, Uganda: an observational epidemiological study
.
Lancet HIV
 
2016
;
3
:
e388
96
.

18

Grabowski
 
MK
,
Lessler
 
J
,
Bazaale
 
J
, et al.  
Migration, hotspots, and dispersal of HIV infection in Rakai, Uganda
.
Nature Communications
 
2020
;
11
:
976
.

19

Galiwango
 
RM
,
Musoke
 
R
,
Lubyayi
 
L
, et al.  
Evaluation of current rapid HIV test algorithms in Rakai, Uganda
.
J Virol Methods
 
2013
;
192
:
25
7
.

20

Ellman
 
TM
,
Alemayehu
 
B
,
Abrams
 
EJ
,
Arpadi
 
S
,
Howard
 
AA
,
El-Sadr
 
WM
.
Selecting a viral load threshold for routine monitoring in resource-limited settings: optimizing individual health and population impact
.
J Int AIDS Soc
 
2017
;
20
:
e25007
.

21

Grabowski
 
MK
,
Reynolds
 
SJ
,
Kagaayi
 
J
, et al.  
The validity of self-reported antiretroviral use in persons living with HIV: a population-based study
.
AIDS
 
2018
;
32
:
363
9
.

22

Uganda Ministry of Health
. National Health Facility Master List 2018: A Complete List of All Health Facilities in Uganda. Government of the Republic of Uganda; 2018. Available at: http://library.health.go.ug/health-infrastructure/health-facility-inventory/national-health-facility-master-facility-list-2018. Accessed 22 November 2022.

23

Maina
 
J
,
Ouma
 
PO
,
Macharia
 
PM
, et al.  
A spatial database of health facilities managed by the public health sector in sub Saharan Africa
.
Sci Data
 
2019
;
6
:
134
.

24

Love
 
RF
,
Morris
 
JG
.
Modelling inter-city road distances by mathematical functions
.
J Oper Res Soc
 
1972
;
23
:
61
71
.

25

Billioux
 
VG
,
Grabowski
 
MK
,
Ssekasanvu
 
J
, et al.  
HIV viral suppression and geospatial patterns of HIV antiretroviral therapy treatment facility use in Rakai, Uganda
.
AIDS
 
2018
;
32
:
819
24
.

26

Santelli
 
JS
,
Chen
 
I
,
Makumbi
 
F
, et al.  
Household wealth and HIV incidence over time, rural Uganda, 1994–2018
.
AIDS
 
2021
;
35
:
1835
43
.

27

Pan
 
W
.
Akaike’s information criterion in generalized estimating equations
.
Biometrics
 
2001
;
57
:
120
5
.

28

Rosen
 
JG
,
Nakyanjo
 
N
,
Isabirye
 
D
, et al.  
Antiretroviral treatment sharing among highly mobile Ugandan fisherfolk living with HIV: a qualitative study
.
AIDS Care
 
2020
;
32
:
912
5
.

29

Schuyler
 
AC
,
Edelstein
 
ZR
,
Mathur
 
S
, et al.  
Mobility among youth in Rakai, Uganda: trends, characteristics, and associations with behavioural risk factors for HIV
.
Glob Public Health
 
2017
;
12
:
1033
50
.

30

Baldé
 
A
,
Lièvre
 
L
,
Maiga
 
AI
, et al.  
Risk factors for loss to follow-up, transfer or death among people living with HIV on their first antiretroviral therapy regimen in Mali
.
HIV Med
 
2019
;
20
:
47
53
.

31

Leslie
 
HH
,
Mooney
 
AC
,
Gilmore
 
HJ
, et al.  
Prevalence, motivation, and outcomes of clinic transfer in a clinical cohort of people living with HIV in North West Province, South Africa
.
BMC Health Serv Res
 
2022
;
22
:
1
10
.

32

Billioux
 
VG
,
Chang
 
LW
,
Reynolds
 
SJ
, et al.  
Human immunodeficiency virus care cascade among sub-populations in Rakai, Uganda: an observational study
.
J Int AIDS Soc
 
2017
;
20
:
21590
.

33

Akullian
 
AN
,
Mukose
 
A
,
Levine
 
GA
,
Babigumira
 
JB
.
People living with HIV travel farther to access healthcare: a population-based geographic analysis from rural Uganda
.
J Int AIDS Soc
 
2016
;
19
:
20171
.

34

Belay
 
YA
,
Yitayal
 
M
,
Atnafu
 
A
,
Taye
 
FA
.
Patient experiences and preferences for antiretroviral therapy service provision: implications for differentiated service delivery in northwest Ethiopia
.
AIDS Res Ther
 
2022
;
19
:
30
.

35

Grimsrud
 
A
,
Bygrave
 
H
,
Doherty
 
M
, et al.  
Reimagining HIV service delivery: the role of differentiated care from prevention to suppression
.
J Int AIDS Soc
 
2016
;
19
:
21484
.

36

Layer
 
EH
,
Brahmbhatt
 
H
,
Beckham
 
SW
, et al.  
“I pray that they accept me without scolding:” experiences with disengagement and re-engagement in HIV care and treatment services in Tanzania
.
AIDS Patient Care STDS
 
2014
;
28
:
483
8
.

37

Bor
 
J
,
Gage
 
A
,
Onoya
 
D
, et al.  
Variation in HIV care and treatment outcomes by facility in South Africa, 2011–2015: a cohort study
.
PLoS Medicine
 
2021
;
18
:
e1003479
.

38

Jaibat
 
L
,
Cui
 
Z
,
Mitra
 
S
, et al.  
Perceived mistreatment in health care settings and its relationship with HIV clinical outcomes in HIV-positive people who use drugs in Vancouver, Canada
.
AIDS Behav
 
2022
;
27
:
1636
46
.

39

Phillips
 
TK
,
Mogoba
 
P
,
Brittain
 
K
, et al.  
Long-term outcomes of HIV-infected women receiving antiretroviral therapy after transferring out of an integrated maternal and child health service in South Africa
.
J Acquir Immune Defic Syndr
 
2020
;
83
:
202
9
.

40

Zakumumpa
 
H
,
Rujumba
 
J
,
Kwiringira
 
J
,
Katureebe
 
C
,
Spicer
 
N
.
Understanding implementation barriers in the national scale-up of differentiated ART delivery in Uganda
.
BMC Health Serv Res
 
2020
;
20
:
222
.

Author notes

Potential conflicts of interest. The authors: No reported conflicts of interest. All authors have submitted the ICMJE Form for Disclosure of Potential Conflicts of Interest.

This article is published and distributed under the terms of the Oxford University Press, Standard Journals Publication Model (https://dbpia.nl.go.kr/pages/standard-publication-reuse-rights)

Supplementary data