Universal Testing, Expanded Treatment, and Incidence of HIV Infection in Botswana.

BACKGROUND: The feasibility of reducing the population-level incidence of human immunodeficiency virus (HIV) infection by increasing community coverage of antiretroviral therapy (ART) and male circumcision is unknown. METHODS: We conducted a pair-matched, community-randomized trial in 30 rural or periurban communities in Botswana from 2013 to 2018. Participants in 15 villages in the intervention group received HIV testing and counseling, linkage to care, ART (started at a higher CD4 count than in standard care), and increased access to male circumcision services. The standard-care group also consisted of 15 villages. Universal ART became available in both groups in mid-2016. We enrolled a random sample of participants from approximately 20% of households in each community and measured the incidence of HIV infection through testing performed approximately once per year. The prespecified primary analysis was a permutation test of HIV incidence ratios. Pair-stratified Cox models were used to calculate 95% confidence intervals. RESULTS: Of 12,610 enrollees (81% of eligible household members), 29% were HIV-positive. Of the 8974 HIV-negative persons (4487 per group), 95% were retested for HIV infection over a median of 29 months. A total of 57 participants in the intervention group and 90 participants in the standard-care group acquired HIV infection (annualized HIV incidence, 0.59% and 0.92%, respectively). The unadjusted HIV incidence ratio in the intervention group as compared with the standard-care group was 0.69 (P = 0.09) by permutation test (95% confidence interval [CI], 0.46 to 0.90 by pair-stratified Cox model). An end-of-trial survey in six communities (three per group) showed a significantly greater increase in the percentage of HIV-positive participants with an HIV-1 RNA level of 400 copies per milliliter or less in the intervention group (18 percentage points, from 70% to 88%) than in the standard-care group (8 percentage points, from 75% to 83%) (relative risk, 1.12; 95% CI, 1.09 to 1.16). The percentage of men who underwent circumcision increased by 10 percentage points in the intervention group and 2 percentage points in the standard-care group (relative risk, 1.26; 95% CI, 1.17 to 1.35). CONCLUSIONS: Expanded HIV testing, linkage to care, and ART coverage were associated with increased population viral suppression. (Funded by the President's Emergency Plan for AIDS Relief and others; Ya Tsie ClinicalTrials.gov number, NCT01965470.).

Sample size considerations in the design of cluster randomized trials of combination HIV prevention.

Background Cluster randomized trials have been utilized to evaluate the effectiveness of HIV prevention strategies on reducing incidence. Design of such studies must take into account possible correlation of outcomes within randomized units. Purpose To discuss power and sample size considerations for cluster randomized trials of combination HIV prevention, using an HIV prevention study in Botswana as an illustration. Methods We introduce a new agent-based model to simulate the community-level impact of a combination prevention strategy and investigate how correlation structure within a community affects the coefficient of variation - an essential parameter in designing a cluster randomized trial. Results We construct collections of sexual networks and then propagate HIV on them to simulate the disease epidemic. Increasing level of sexual mixing between intervention and standard-of-care (SOC) communities reduces the difference in cumulative incidence in the two sets of communities. Fifteen clusters per arm and 500 incidence cohort members per community provide 95% power to detect the projected difference in cumulative HIV incidence between SOC and intervention communities (3.93% and 2.34%) at the end of the third study year, using a coefficient of variation 0.25. Although available formulas for calculating sample size for cluster randomized trials can be derived by assuming an exchangeable correlation structure within clusters, we show that deviations from this assumption do not generally affect the validity of such formulas. Limitations We construct sexual networks based on data from Likoma Island, Malawi, and base disease progression on longitudinal estimates from an incidence cohort in Botswana and in Durban as well as a household survey in Mochudi, Botswana. Network data from Botswana and larger sample sizes to estimate rates of disease progression would be useful in assessing the robustness of our model results. Conclusion Epidemic modeling plays a critical role in planning and evaluating interventions for prevention. Simulation studies allow us to take into consideration available information on sexual network characteristics, such as mixing within and between communities as well as coverage levels for different prevention modalities in the combination prevention package.

The vitamin D for COVID-19 (VIVID) trial: A pragmatic cluster-randomized design.

OBJECTIVES: To determine the effect of vitamin D supplementation on disease progression and post-exposure prophylaxis for COVID-19 infection. We hypothesize that high-dose vitamin D3 supplementation will reduce risk of hospitalization/death among those with recently diagnosed COVID-19 infection and will reduce risk of COVID-19 infection among their close household contacts. METHODS: We report the rationale and design of a planned pragmatic, cluster randomized, double-blinded trial (N = 2700 in total nationwide), with 1500 newly diagnosed individuals with COVID-19 infection, together with up to one close household contact each (~1200 contacts), randomized to either vitamin D3 (loading dose, then 3200 IU/day) or placebo in a 1:1 ratio and a household cluster design. The study duration is 4 weeks. The primary outcome for newly diagnosed individuals is the occurrence of hospitalization and/or mortality. Key secondary outcomes include symptom severity scores among cases and changes in the infection (seroconversion) status for their close household contacts. Changes in vitamin D 25(OH)D levels will be assessed and their relation to study outcomes will be explored. CONCLUSIONS: The proposed pragmatic trial will allow parallel testing of vitamin D3 supplementation for early treatment and post-exposure prophylaxis of COVID-19. The household cluster design provides a cost-efficient approach to testing an intervention for reducing rates of hospitalization and/or mortality in newly diagnosed cases and preventing infection among their close household contacts.

Population uptake of HIV testing, treatment, viral suppression, and male circumcision following a community-based intervention in Botswana (Ya Tsie/BCPP): a cluster-randomised trial.

BACKGROUND: In settings with high HIV prevalence and treatment coverage, such as Botswana, it is unknown whether uptake of HIV prevention and treatment interventions can be increased further. We sought to determine whether a community-based intervention to identify and rapidly treat people living with HIV, and support male circumcision could increase population levels of HIV diagnosis, treatment, viral suppression, and male circumcision in Botswana. METHODS: The Ya Tsie Botswana Combination Prevention Project study was a pair-matched cluster-randomised trial done in 30 communities across Botswana done from Oct 30, 2013, to June 30, 2018. 15 communities were randomly assigned to receive HIV prevention and treatment interventions, including enhanced HIV testing, earlier antiretroviral therapy (ART), and strengthened male circumcision services, and 15 received standard of care. The first primary endpoint of HIV incidence has already been reported. In this Article, we report findings for the second primary endpoint of population uptake of HIV prevention services, as measured by proportion of people known to be HIV-positive or tested HIV-negative in the preceding 12 months; proportion of people living with HIV diagnosed and on ART; proportion of people living with HIV on ART with viral suppression; and proportion of HIV-negative men circumcised. A longitudinal cohort of residents aged 16-64 years from a random, approximately 20% sample of households across the 15 communities was enrolled to assess baseline uptake of study outcomes; we also administered an end-of-study survey to all residents not previously enrolled in the longitudinal cohort to provide study end coverage estimates. Differences in intervention uptake over time by randomisation group were tested via paired Student's t test. The study has been completed and is registered with ClinicalTrials.gov (NCT01965470). FINDINGS: In the six communities participating in the end-of-study survey, 2625 residents (n=1304 from standard-of-care communities, n=1321 from intervention communities) were enrolled into the 20% longitudinal cohort at baseline from Oct 30, 2013, to Nov 24, 2015. In the same communities, 10 791 (86%) of 12 489 eligible enumerated residents not previously enrolled in the longitudinal cohort participated in the end-of-study survey from March 30, 2017, to Feb 25, 2018 (5896 in intervention and 4895 in standard-of-care communities). At study end, in intervention communities, 1228 people living with HIV (91% of 1353) were on ART; 1166 people living with HIV (88% of 1321 with available viral load) were virally suppressed, and 673 HIV-negative men (40% of 1673) were circumcised in intervention communities. After accounting for baseline differences, at study end the proportion of people living with HIV who were diagnosed was significantly higher in intervention communities (absolute increase of 9% to 93%) compared with standard-of-care communities (absolute increase of 2% to 88%; prevalence ratio [PR] 1·08 [95% CI 1·02-1·14], p=0·032). Population levels of ART, viral suppression, and male circumcision increased from baseline in both groups, with greater increases in intervention communities (ART PR 1·12 [95% CI 1·07-1·17], p=0·018; viral suppression 1·13 [1·09-1·17], p=0·017; male circumcision 1·26 [1·17-1·35], p=0·029). INTERPRETATION: It is possible to achieve very high population levels of HIV testing and treatment in a high-prevalence setting. Maintaining these coverage levels over the next decade could substantially reduce HIV transmission and potentially eliminate the epidemic in these areas. FUNDING: US President's Emergency Plan for AIDS Relief through the Centers for Disease Control and Prevention.

Linear Growth Faltering Among HIV-Exposed Uninfected Children.

BACKGROUND: HIV-exposed uninfected (HEU) children experience increased mortality compared with their HIV-unexposed uninfected (HUU) peers. It is unclear whether HEU children are also at increased risk for undernutrition, a modifiable risk factor for mortality. METHODS: We conducted a cross-sectional, population-based survey of children <5 years of age in 5 health districts in Botswana. Linear mixed-effects models were used to assess continuous outcomes, and generalized estimating equations were used to estimate relative risks of stunting, wasting, and underweight between HEU (n = 396) and HUU (n = 1109) children. Secondary analyses examined potential mediation by low birth weight. RESULTS: The association between maternal HIV exposure and child stunting varied significantly by child age (P < 0.01). HEU children <1 and ≥2 years of age had 1.85 [95% confidence interval (CI): 1.03 to 3.31; P = 0.04] and 1.41 (95% CI: 1.06 to 1.88; P = 0.02) times the risk of stunting compared with HUU children after multivariate adjustment, respectively. During the period of 1-2 years of age, when breastfeeding cessation occurred among HUU children, HUU children had increased risk of stunting compared with HEU children who were predominantly formula fed (relative risk: 1.56; 95% CI: 1.05 to 2.32; P = 0.03). A mediation analysis estimated that 67% of the excess risk of stunting among HEU children ≥2 years was attributable to low birth weight (P = 0.02). There was no difference in risk of wasting or underweight. CONCLUSION: HEU children are at increased risk of stunting compared with their HUU peers; however, interventions to increase birth weight may significantly ameliorate this excess risk. Interventions to support optimal growth during weaning are needed for all breast-fed children.

Botswana's progress toward achieving the 2020 UNAIDS 90-90-90 antiretroviral therapy and virological suppression goals: a population-based survey.

BACKGROUND: HIV programmes face challenges achieving high rates of HIV testing and treatment needed to optimise health and to reduce transmission. We used data from the Botswana Combination Prevention Project study survey to assess Botswana's progress toward achieving UNAIDS targets for 2020: 90% of all people living with HIV knowing their status, 90% of these receiving sustained antiretroviral therapy (ART), and 90% of those having virological suppression (90-90-90). METHODS: A population-based sample of individuals was recruited and interviewed in 30 rural and periurban communities from Oct 30, 2013, to Nov 24, 2015, as part of a large, ongoing community-randomised trial designed to assess the effect of a combination prevention package on HIV incidence. A random sample of about 20% of households in each community was selected. Consenting household residents aged 16-64 years who were Botswana citizens or spouses of citizens responded to a questionnaire and had blood drawn for HIV testing in the absence of documentation of positive HIV status. Viral load testing was done in all HIV-infected participants, irrespective of treatment status. We used modified Poisson generalised estimating equations to obtain prevalence ratios, corresponding Huber robust SEs, and 95% Wald CIs to examine associations between individual sociodemographic factors and a binary outcome indicating achievement of the three individual and combined overall 90-90-90 targets. The study is registered at ClinicalTrials.gov, number NCT01965470. FINDINGS: 81% of enumerated eligible household members took part in the survey (10% refused and 9% were absent). Among 12 610 participants surveyed, 3596 (29%) were infected with HIV, and 2995 (83·3%, 95% CI 81·4-85·2) of these individuals already knew their HIV status. Among those who knew their HIV status, 2617 (87·4%, 95% CI 85·8-89·0) were receiving ART (95% of those eligible by national guidelines, and 73% of all infected people). Of the 2609 individuals receiving ART with a viral load measurement, 2517 (96·5%, 95% CI 96·0-97·0) had viral load of 400 copies per mL or less. Overall, 70·2% (95% CI 67·5-73·0) of HIV-infected people had virological suppression, close to the UNAIDS target of 73%. INTERPRETATION: UNAIDS 90-90-90 targets are achievable even in resource-constrained settings with high HIV burden. FUNDING: US President's Emergency Plan for AIDS Relief through the Centers for Disease Control and Prevention.

Importance of Viral Sequence Length and Number of Variable and Informative Sites in Analysis of HIV Clustering.

To improve the methodology of HIV cluster analysis, we addressed how analysis of HIV clustering is associated with parameters that can affect the outcome of viral clustering. The extent of HIV clustering and tree certainty was compared between 401 HIV-1C near full-length genome sequences and subgenomic regions retrieved from the LANL HIV Database. Sliding window analysis was based on 99 windows of 1,000 bp and 45 windows of 2,000 bp. Potential associations between the extent of HIV clustering and sequence length and the number of variable and informative sites were evaluated. The near full-length genome HIV sequences showed the highest extent of HIV clustering and the highest tree certainty. At the bootstrap threshold of 0.80 in maximum likelihood (ML) analysis, 58.9% of near full-length HIV-1C sequences but only 15.5% of partial pol sequences (ViroSeq) were found in clusters. Among HIV-1 structural genes, pol showed the highest extent of clustering (38.9% at a bootstrap threshold of 0.80), although it was significantly lower than in the near full-length genome sequences. The extent of HIV clustering was significantly higher for sliding windows of 2,000 bp than 1,000 bp. We found a strong association between the sequence length and proportion of HIV sequences in clusters, and a moderate association between the number of variable and informative sites and the proportion of HIV sequences in clusters. In HIV cluster analysis, the extent of detectable HIV clustering is directly associated with the length of viral sequences used, as well as the number of variable and informative sites. Near full-length genome sequences could provide the most informative HIV cluster analysis. Selected subgenomic regions with a high extent of HIV clustering and high tree certainty could also be considered as a second choice.

Estimated age and gender profile of individuals missed by a home-based HIV testing and counselling campaign in a Botswana community.

INTRODUCTION: It would be useful to understand which populations are not reached by home-based HIV-1 testing and counselling (HTC) to improve strategies aimed at linking these individuals to care and reducing rates of onward HIV transmission. METHODS: We present the results of a baseline home-based HTC (HBHTC) campaign aimed at counselling and testing residents aged 16 to 64 for HIV in the north-eastern sector of Mochudi, a community in Botswana with about 44,000 inhabitants. Collected data were compared with population references for Botswana, the United Nations (UN) estimates based on the National Census data and the Botswana AIDS Impact Survey IV (BAIS-IV). Analyzed data and references were stratified by age and gender. RESULTS: A total of 6238 age-eligible residents were tested for HIV-1; 1247 (20.0%; 95% CI 19.0 to 21.0%) were found to be HIV positive (23.7% of women vs. 13.4% of men). HIV-1 prevalence peaked at 44% in 35- to 39-year-old women and 32% in 40- to 44-year-old men. A lower HIV prevalence rate, 10.9% (95% CI 9.5 to 12.5%), was found among individuals tested for the first time. A significant gender gap was evident in all analyzed subsets. The existing HIV transmission network was analyzed by combining phylogenetic mapping and household structure. Between 62.4 and 71.8% of all HIV-positive individuals had detectable virus. When compared with the UN and BAIS-IV estimates, the proportion of men missed by the testing campaign (48.5%; 95% CI 47.0 to 50.0%) was significantly higher than the proportion of missed women (14.2%; 95% CI 13.2 to 15.3%; p<0.0001). The estimated proportion of missed men peaked at about 60% in the age group 30 to 39 years old. The proportions of missed women were substantially smaller, at approximately 28% within the age groups 30 to 34 and 45 to 49 years old. CONCLUSIONS: The HBHTC campaign seems to be an efficient tool for reaching individuals who have never been tested previously in southern African communities. However, about half of men from 16 to 64 years old were not reached by the HBHTC, including about 60% of men between 30 and 40 years old. Alternative HTC strategies should be developed to bring these men to care, which will contribute to reduction of HIV incidence in communities.

Impact of sampling density on the extent of HIV clustering.

Identifying and monitoring HIV clusters could be useful in tracking the leading edge of HIV transmission in epidemics. Currently, greater specificity in the definition of HIV clusters is needed to reduce confusion in the interpretation of HIV clustering results. We address sampling density as one of the key aspects of HIV cluster analysis. The proportion of viral sequences in clusters was estimated at sampling densities from 1.0% to 70%. A set of 1,248 HIV-1C env gp120 V1C5 sequences from a single community in Botswana was utilized in simulation studies. Matching numbers of HIV-1C V1C5 sequences from the LANL HIV Database were used as comparators. HIV clusters were identified by phylogenetic inference under bootstrapped maximum likelihood and pairwise distance cut-offs. Sampling density below 10% was associated with stochastic HIV clustering with broad confidence intervals. HIV clustering increased linearly at sampling density >10%, and was accompanied by narrowing confidence intervals. Patterns of HIV clustering were similar at bootstrap thresholds 0.7 to 1.0, but the extent of HIV clustering decreased with higher bootstrap thresholds. The origin of sampling (local concentrated vs. scattered global) had a substantial impact on HIV clustering at sampling densities ≥10%. Pairwise distances at 10% were estimated as a threshold for cluster analysis of HIV-1 V1C5 sequences. The node bootstrap support distribution provided additional evidence for 10% sampling density as the threshold for HIV cluster analysis. The detectability of HIV clusters is substantially affected by sampling density. A minimal genotyping density of 10% and sampling density of 50-70% are suggested for HIV-1 V1C5 cluster analysis.

Phylogenetic relatedness of circulating HIV-1C variants in Mochudi, Botswana.

BACKGROUND: Determining patterns of HIV transmission is increasingly important for the most efficient use of modern prevention interventions. HIV phylogeny can provide a better understanding of the mechanisms underlying HIV transmission networks in communities. METHODS: To reconstruct the structure and dynamics of a local HIV/AIDS epidemic, the phylogenetic relatedness of HIV-1 subtype C env sequences obtained from 785 HIV-infected community residents in the northeastern sector of Mochudi, Botswana, during 2010-2013 was estimated. The genotyping coverage was estimated at 44%. Clusters were defined based on relatedness of HIV-1C env sequences and bootstrap support of splits. RESULTS: The overall proportion of clustered HIV-1C env sequences was 19.1% (95% CI 17.5% to 20.8%). The proportion of clustered sequences from Mochudi was significantly higher than the proportion of non-Mochudi sequences that clustered, 27.0% vs. 14.7% (p = 5.8E-12; Fisher exact test). The majority of clustered Mochudi sequences (90.1%; 95% CI 85.1% to 93.6%) were found in the Mochudi-unique clusters. None of the sequences from Mochudi clustered with any of the 1,244 non-Botswana HIV-1C sequences. At least 83 distinct HIV-1C variants, or chains of HIV transmission, in Mochudi were enumerated, and their sequence signatures were reconstructed. Seven of 20 genotyped seroconverters were found in 7 distinct clusters. CONCLUSIONS: The study provides essential characteristics of the HIV transmission network in a community in Botswana, suggests the importance of high sampling coverage, and highlights the need for broad HIV genotyping to determine the spread of community-unique and community-mixed viral variants circulating in local epidemics. The proposed methodology of cluster analysis enumerates circulating HIV variants and can work well for surveillance of HIV transmission networks. HIV genotyping at the community level can help to optimize and balance HIV prevention strategies in trials and combined intervention packages.

Improved double-robust estimation in missing data and causal inference models.

Recently proposed double-robust estimators for a population mean from incomplete data and for a finite number of counterfactual means can have much higher efficiency than the usual double-robust estimators under misspecification of the outcome model. In this paper, we derive a new class of double-robust estimators for the parameters of regression models with incomplete cross-sectional or longitudinal data, and of marginal structural mean models for cross-sectional data with similar efficiency properties. Unlike the recent proposals, our estimators solve outcome regression estimating equations. In a simulation study, the new estimator shows improvements in variance relative to the standard double-robust estimator that are in agreement with those suggested by asymptotic theory.