Introduction

In the context of Human Immunodeficiency Virus (HIV) treatment, adherence to antiretroviral medication is an important health behavior that facilitates viral suppression, the driver of better health outcomes and prevention of onward HIV transmission. The sexual partnership within HIV-serodifferent couples, where one partner is living with HIV and the other is HIV-negative, puts the HIV-negative partner at high risk for HIV acquisition when adequate strategies for HIV prevention and treatment are not in place [1, 2]. Programs for HIV-serodifferent couples focus on supporting HIV treatment for the partner living with HIV and time-limited pre-exposure prophylaxis (PrEP) using antiretroviral medication for HIV-negative partners, when HIV viremia may not be suppressed in the partner living with HIV [3,4,5]. When PrEP is integrated into antiretroviral treatment (ART) programs for HIV-serodifferent couples, it bridges the immediate HIV prevention gap prior to viral suppression in the partner living with HIV and has been demonstrated to have 95% effectiveness in reducing HIV incidence [6]. A randomized trial of oral antiretroviral therapy for HIV-serodiscordant heterosexual couples in Kenya and Uganda found 67% and 75% reduction in the incidence of HIV with tenofovir and combination tenofovir–emtricitabine respectively [7]. The HIV-negative partners chose to take PrEP to reduce the HIV transmission anxiety, improve the relationship, and together build shared protection against HIV as a couple [8].

Behavioral modeling, or observational learning, refers to the process of individuals learning and acquiring new behaviors from watching others [9]. Observational learning has four processes, including attention, retention, production, and motivation. Individuals first direct their attention to the modeled behaviors or interactions (attention) and remember them (retention), which is dependent on individuals’ intellectual capacities. Individuals then perform the modeled behavior when the circumstances and physical skills allow them to (production). Motivation is the internal factor that drives the production of the modeled behavior, depending on the expectations about costs and benefits [7]. Within HIV-serodifferent couples using antiretrovirals for prevention and treatment, behavioral modeling is one mechanism that could explain each partner’s medication adherence behavior. Specifically, partners living with HIV could build self-efficacy for taking HIV medications by observing their HIV-negative partners manage medication side effects, adherence, and HIV-related stigma. The perceived benefits on their health, relationship quality, and shared future with their partners may also motivate ART adherence [10]. Observing the partners taking ART could reinforce PrEP adherence among HIV-negative partners in return, creating a positive feedback loop within the partners.

To our knowledge, no measure has been developed yet to capture whether heterosexual HIV-serodifferent couples exhibit behavioral modeling in their antiretroviral use (PrEP for prevention, ART for treatment). The current study aimed to (1) use factor analysis to identify the underlying factors that the items on a questionnaire are measuring and examine the internal consistency of a new behavioral modeling scale; (2) examine the association between behavioral modeling with medication-taking behaviors among HIV-serodifferent couples; and (3) explore the association between behavioral modeling with viral load levels among partners with HIV and Tenofovir (TFV) drug levels among HIV-negative partners.

Methods

Study Design and Participants

This study is a longitudinal observational study that was nested within the Partners PrEP Program, a stepped wedge cluster randomized trial evaluating the integration of PrEP delivery into ART programs (#NCT03586128) 11] in Kampala and Wakiso, Uganda from June 2018 to December 2021. Heterosexual sero-different couples were recruited into the program through couples-based HIV testing and received access to standard of care for HIV care after the intervention was launched at their clinic.

Participant inclusion criteria were as follows. For both members of the couple: (1) Age ≥ 18 years; (2) Ability and willingness to provide informed consent; (3) Being sexually active with each other; (4) Willingness to engage with the clinic system collaboratively. For the HIV-positive members (index participants): (1) Diagnosis as HIV-positive per the national HIV testing algorithm; (2) Recent identification as a member of an HIV serodiscordant couple; (3) Absence of current enrollment in any HIV treatment clinical trial. For the HIV-negative members (partner participants): (1) Diagnosis as HIV-negative using the national HIV testing algorithm; (2) Recent identification within an HIV serodiscordant relationship; (3) Absence of current enrollment in any HIV treatment clinical trial; (4) Non-usage of PrEP at present; (5) Eligibility for PrEP as defined by the WHO or Ugandan national guidelines [12].

Programmatic data were abstracted from patient medical charts including dates when PrEP and ART were dispensed and HIV viral load results from testing that occurred approximately 3–6 months after ART initiation in alignment with national guidelines. Additionally, couples were invited to enroll for supplemental research procedures that included interviewer-administered questionnaires on every quarter for about 24 months to ascertain data on demographics, sexual behavior, modeled behavior, and self-reported PrEP and ART adherence. A dried blood spot was also collected from HIV-negative partners in the research cohort and a random sample of 10% were selected for quantification of tenofovir (TFV)-diphosphate levels via validated liquid chromatography/tandem mass spectrometry methods [13].

PrEP Delivery Program

The 12 ART clinics received training on PrEP delivery, ongoing technical assistance visits from the training team, and PrEP commodities before the program was launched. The PrEP training protocol was adapted from the national PrEP training curriculum in Uganda, with adjustments to tailor the service to HIV-serodifferent couples. The technical assistance visits were conducted by the training team who were experienced in PrEP delivery and clinical management. Once clinics entered the PrEP delivery stage, trained staff offered PrEP to the HIV-negative members of recently diagnosed HIV-serodifferent couples [14]. HIV-negative partners were counseled to use PrEP for at least 6 months after their partners started ART and supported to discontinue when the partner sustained ART use and there were no other potential sources of exposure (e.g., other partners).

Measurements

Behavioral Modeling Scale

Through discussions and initial pilot testing with a group of clinicians, ART clients and their partners, and research staff, we developed a 28-item questionnaire with a four-point Likert-type response scale to measure behavioral modeling. The items were created based on the observational learning concepts of the Social Cognitive Theory [7]. Nineteen items mapped to the four processes of observational learning: attention (e.g., “I often see my partner taking PrEP/ART”), retention (e.g., “I remind my partner to take PrEP”), production (e.g., “I take my ART at the same time my partner takes PrEP/ART”), and motivation (e.g., “When my partner takes PrEP/ART I have hope for our future together”) [9]. Additionally, we added 9 items to capture the relationship quality between the partners (e.g., “I feel close to my partner”). These items were based on previous research that suggested a strong association between relationship satisfaction and PrEP adherence, and that role modeling is more effective when it comes from someone the observer feels close to and values [15]. Four items were removed due to poor face validity. Face validity refers to whether a question appears to measure the construct it aims to measure [16]. For example, one item that was removed was “I have people I can talk to about my relationship when I need support”, which was more about social support than behavioral modeling.

The final scale contained 24 items. The items were identical for both partners within a couple, with “PrEP” or “ART” used respectively for the HIV-negative partners and partners living with HIV. The response options ranged from (0) “Strongly disagree” to (3) “Strongly agree”. The scale was scored by computing a total raw score ranging from 0 to 72. Separate total scores were computed once subscales were identified. We computed total scores for both HIV-negative partners and the partners living with HIV to capture each partner’s self-report of the other partners’ behaviors that could have been observed.

Self-Reported ART/PrEP-Taking Behavior

Medication -aking behaviors were measured using three self-reported items adapted from Wilson et al. (1) whether the partner has taken ART/PrEP since the last clinic visit (yes/no); (2) how well the client took their ART/PrEP medication as directed in the past month (0 = very poor to 5 = excellent), and (3) how often the client took their ART/PrEP medication in the past month (0 = not at all to 3 = every day). We performed linear transformation on the item responses for the second and the third item to compose the adherence score (with a range 0-100), with higher scores indicating better adherence [17].

Statistical Analysis

All analyses were conducted in R [18]. Demographic characteristics of the sample and medication-taking behaviors across three time points (months 1, 3, and 6 after enrollment) were presented and compared between HIV-negative and -positive partners using Chi-square tests for categorical variables and ANOVA for continuous variables. We conducted an exploratory factor analysis (EFA) to examine the structure of the behavioral modeling scale, using data from the first follow-up visit (e.g., 1 month after enrollment). The rationale for using factor analysis was based on this being a novel scale without previous exploration of its structure in this specific cultural and study setting. We used the promax factor rotation method that allows factor correlations and results in simple and interpretable factor structure [19]. The final factor structure was determined after an evaluation of eigenvalues, residual variances, chi-square, Bayesian Information Criterion (BIC), Comparative Fit Index (CFI), Root Mean Square Error of Approximation (RMSEA), and variance accounted for by different factor structures [20]. The internal consistency of the scale is indicated by Cronbach’s α. The data from both HIV-positive and negative partners were combined to conduct EFA and examine the internal consistency of the scale.

We examined the association between a participant’s behavioral modeling scale and subscales and their own self-reported PrEP/ART adherence score across three time points using generalized estimating equations (GEE) to account for repeated measurements. The regression coefficients and 95% confidence interval (CI) from the GEE model indicated the strength of the associations and the uncertainty of the estimates. We used logistic regression to examine the associations between the behavioral modeling scale and (a) outcome of viral suppression (HIV RNA < 1000 copies/ml) at 6 months post ART initiation and, in separate models, (b) outcome of TFV levels associated with high adherence (TFV > 700 fmol/punch) [21]. The odds ratio and 95%CI indicated the strength and precision of observed associations. Since we had only one time point with viral load and TVF outcome, behavioral model scale scores were averaged across month 1, 3, and 6 for this analysis. We also examined associations between behavioral modeling scale score and continuous TFV levels and viral load outcomes using linear regression. For all analyses, adjusted models were generated including participant age (continuous) and sex (male, female) as a priori determined confounding factors.

Ethics

The protocol was approved by ethics committees at the Uganda National HIV/AIDS Research Committee (ARC 194), the Uganda National Council for Science and Technology (HS 2381), and the University of Washington Human Subjects Division (STUDY00000320). The Kampala Capital City Authority and Wakiso District gave permission for data to be abstracted from clinic medical records. All participants for research procedures underwent an informed consent process and provided written consent in their preferred language.

Results

Descriptive Characteristics

In total, both members of 149 serodifferent couples enrolled for the research cohort at baseline, with 141 retained at month 1, 3 and 139 at month 6. More than half of the partners living with HIV were female (64%). The median age of the participants was 28 years (IQR 9.8). Based on self-report adherence score, the partners living with HIV were more adherent to ART than their partners were to PrEP, across all three follow-up visits (Table 1).

Table 1 Baseline characteristics of HIV-serodifferent couples

Exploratory Factor Analysis

The Kaiser-Meyer-Olkin index of sampling adequacy for the behavioral modeling scale was 0.92, suggesting the data were suitable for factor analysis. Eigenvalues from factor analysis with promax rotation suggested a 5-factor model would be the most explanatory, with the first five components having an eigenvalue larger than one and the five factors accounting for 66% of the total variance (Table 2). As the number of factors increased from three to seven, Chi-square statistics and RMSEA both decreased, while CFI increased, suggesting improvement in model performance. Model fit index BIC decreased initially as the number of factors increased from three to five, suggesting improvement in model performance. BIC started to increase after the factor number exceeded five, suggesting that the model performance started to decrease. The fit statistics together supported the five-factor model. Additionally, the five-factor model provided a meaningful interpretation of the factor structure based on the factor loadings (Table 3). The five factors were: (1) attention to partner behavior, (2) collective action, (3) partner as a role model, (4) motivation, and (5) relationship quality. The mean and SD of factor correlations were presented in Table 4. When we split the dataset by HIV status, the factor structure was replicated in both subgroups (Supplement Tables 1, 2 and 3).

Table 2 Fit statistics of explorative factor analysis models
Table 3 Factor loadings of five-factor model
Table 4 Factor correlations

Internal Consistency

The scale demonstrated a high level of internal consistency (Cronbach’s α = 0.93, 95% confidence interval [CI]: 0.92–0.95). The five subscales also demonstrated high internal consistency: relationship quality (α = 0.93, 95% CI: 0.92–0.94), attention to partner behavior (α = 0.92, 95% CI: 0.91–0.93), partner as a role model (α = 0.85, 95% CI: 0.81–0.89), motivation (α = 0.83, 95% CI: 0.80–0.87), and collective action (α = 0.80, 95% CI: 0.76–0.84). The internal consistency of the scale was similar between partners living with and without HIV (Supplement Table 4).

Behavioral Modeling, Adherence Score, and Viral Suppression

The mean of the total score for our 24-item behavioral modeling measure was 57.2 (SD = 11.0) at month 1, 56.3 (SD = 12.6) at month 3, and 56.4 (SD = 2.5) at month 6, with no statistically significant difference between partners of the same couple. Compared to partners living with HIV, HIV-negative partners reported higher scores on motivation and role modeling subscales at 1-month (t(280) = 3.34, p < 0.001; t(180) = 2.74, p = 0.006) and 3-month follow-up visits (t(278) = 4.69, p < 0.001; t(278) = 3.26, p = 0.001). Compared to partners living with HIV, HIV-negative partners reported lower scores on relationship subscale (t(280) = -2.99, p = 0.005) and collective action (t(280) = -2.27, p = 0.026) at 6-month follow-up visit (Table 5).

Table 5 Behavioral modeling score at all visits by HIV status, mean (SD)

HIV-negative partners’ behavior modeling of PrEP taking, indicated by the composite behavioral modeling scale score reported by HIV-positive partners, was positively associated with their HIV-positive partners’ ART adherence (β = 0.06, 95% CI: 0.01–0.10) among partners living with HIV, after accounting for age, and sex. In separate GEE models, subscales of motivation (β = 0.18, 95% CI: 0.02–0.38), relationship quality (β = 0.19, 95% CI: 0.02–0.33), attention (β = 0.18, 95% CI: 0.02–0.35), and role modeling (β = 0.50, 95% CI: 0.07–0.93) were positively associated with self-reported ART adherence score (Table 6).

Table 6 Associations between behavioral modeling scores and self-reported ART adherence score (partners living with HIV)

HIV-positive partners’ behavioral modeling of ART taking, indicated by the composite behavioral modeling scale score reported by HIV-negative partners, was positively associated with their self-reported PrEP adherence (β = 0.25, 95% CI: 0.17–0.33) after controlling for age and sex. In separate GEE models, subscales for motivation (β = 1.09, 95% CI: 0.71–1.47), relationship quality (β = 0.56, 95% CI: 0.34–0.78), attention (β = 0.69, 95% CI: 0.39–0.99), collective action (β = 0.64, 95% CI: 0.35–0.93), and role modeling (β = 2.76, 95% CI: 1.87–3.66) were positively associated with ART adherence score (Table 7).

Table 7 GEE models: behavioral modeling and self-reported PrEP adherence score (HIV-negative partners)

When averaged across three follow-up visits, neither the behavioral modeling scale composite score, nor subscale scores were associated with either dichotomized or continuous parameterizations of viral load or TFV levels at the month 6 follow-up visit (Tables 8, 9, 10 and 11). The results were similar when limited to data from the 6-month time point. Self-reported ART adherence mapped to viral load data with 91.0% of people who took ART every day having suppressed viremia, a greater frequency than those who reported gaps in taking ART. Similarly, self-reported PrEP use was consistent with TFV level data, with 66.7% of people who took PrEP every day having detectable TFV levels, which was more frequent than people who self-reported low adherence to PrEP (Tables 12 and 13).

Table 8 Associations between modeling behavior and HIV viral load (continuous)
Table 9 Association between modeling behavior and viral suppression (binary, HIV RNA < 1000 copies/ml)
Table 10 Association between modeling behavior and TFV diphosphate levels (continuous)
Table 11 Association between modeling behavior and TFV detection (binary, TFV > 700 fmol/punch)
Table 12 Concordance between ART adherence and viral suppression (HIV RNA < 1000 copies/ml)
Table 13 Concordance between PrEP adherence and TFV detection (TFV > 700 fmol/punch)

Discussion

To our knowledge, this is the first study that developed a scale to measure dyadic behavioral modeling within HIV-serodifferent couples. We found that modeled ART-taking behavior was associated with PrEP use and that modeled PrEP-taking behavior was associated with ART use. This association was bidirectional irrespective of sex. The role model subscale was positively associated with self-reported adherence to ART and PrEP.

The final factor structure of the scale (attention, collective action, motivation, relationship, and role model) mapped well onto the four processes of observational learning. Notably, we grouped the items that were intended to measure retention (item 13–15) and production (item 16–17) and named this factor “collective action”, representing the joint efforts of both partners to support each other in medication adherence. Collective action is a construct of unique value in the context of measuring modeled behavior within a couple. While social cognitive theory recognizes individuals’ ability to actively construct their environment to achieve their individual goals, it also emphasizes collective action, where a group of individuals work together to enact change that can benefit the group [9]. The collective action within the couple, such as reminding each other to take medication, serves the goals of maintaining optimal health for both partners and preserving the long-term relationship. Prior work also suggested that the shared experience of taking antiretroviral medication together which often involves mutual reminders, and providing emotional and material support to each other, has been documented to draw members of the couples closer to each other and prevent their separation, contributing to the success of integrated delivery of ART and PrEP to HIV serodifferent couples [22].

Additionally, the measure included a “partner as a role model” factor that reflected a partner’s medication adherence behavior. Although this factor is composed of only two items, the role model subscale consistently predicted the adherence score of ART/PrEP and had the largest coefficient among all subscales. This is consistent with the literature on role modeling and behavior change. A study using peers as role models who were similar to participants along multiple dimensions (e.g., sex, age, occupation) found that the peer modeling intervention increased participants’ physical fitness, compared to the group that only received general health information [23]. A review of peer-based interventions on HIV-related outcomes found that peer-based interventions increased linkage and retention to HIV care [24]. The improved health behaviors might be due to the positive effect of exposure to peer models on participant confidence and knowledge. The behaviors of the partner might contribute to the establishment of a norm within the couple that is mutually reinforcing the adherent behavior [25].

When applied in clinical settings and research, our novel behavioral modeling scale can be used to outline intervention targets across five dimensions (motivation, relationship, attention, collective action, role modeling) for members of serodifferent couples to improve medication adherence and HIV prevention and treatment outcomes. This is consistent with the perspectives of sero-different couples where the participants reported factors such as relationship quality and share protection that motivated them to adopt HIV prevention strategies coupled with viral control [8]. HIV-serodifferent couples are uniquely positioned where the sustained HIV viremic control through ART and maximizing prevention benefit of PrEP can be attained through collective medication-taking efforts from both partners. Intervention programs that strategically target the different dimensions of behavioral modeling between partners can create a sense of collective action and promote optimal adherence to antiretrovirals of both partners. This measure was developed for the specific context of serodifferent couples to capture dimensions of dyadic behavioral modeling that are associated with adherence behavior in HIV prevention and treatment. However, this measure has the potential to be adapted and generalized to the measurement of dyadic behavioral modeling effect in other health-related research. Further studies are needed to examine the psychometric properties of this measure in other health behavior settings and cultural contexts.

This study has several limitations. First, we were not able to utilize all three self-reported adherence items when calculating the adherence scores due to the adaption of the first item in this study. Future studies may use the original adherence items to test the predictive effects of this behavioral modeling. Second, despite the concordance between self-reported adherence data and biological outcomes (viral suppression and TFV level), our data are insufficiently powered to determine whether the scale and subscales are associated with objective adherence measures, viral suppression and TFV. Future studies can collect a larger sample of biological outcomes and examine whether behavioral modeling has differential predictive effects for self-reported versus biological outcomes.

Conclusions

Our novel behavioral modeling measure successfully assessed five dimensions of modeled behavior (i.e., attention, motivation, relationship quality, collective action, and role modeling), which are components of observational learning. Furthermore, our results suggested the behavioral modeling measure was associated with self-reported medication-adherence behavior among members of HIV-serodifferent couples. With further validation, this measure could be applied in clinical and research settings to outline modifiable behavioral change within partners using the five dimensions of behavioral modeling of medication taking. Furthermore, the results suggest that further interventions with HIV-serodifferent couples could think beyond role modeling, and include dimensions such as collective action and relationship quality between couples, to maximize adherence to antiretrovirals and protection against HIV transmission.