Effects of Music and Counseling Programs on Mental Health and Antiretroviral Treatment Adherence of Adolescents and Young Adults Living With HIV in Northern Thailand,a Randomized Trial

Abstract

Background

Mental health significantly impacts antiretroviral treatment (ART) adherence in adolescents and young adults with perinatal HIV (AYA-PHIV), posing challenges for healthcare providers. This study investigated the effects of music and counseling programs on mental health and treatment adherence of AYA-PHIV.

Methods

In a randomized controlled study, the participants were assigned to the Music and Counseling Program (Arm 1), the Music Program (Arm 2), the Counseling Program (Arm 3), and a no-intervention group (Arm 4). Before and after 10 weekly sessions of a 90-minute study intervention, ART adherence and mental health questionnaires were completed.

Results

After the intervention, Arm 1 participants had a significantly increased median self-esteem score (17-25, P = .05). A significantly decrease in mean Thai-IHSS scores were observed in Arm 1 and 3 participants.

Conclusions

The study demonstrated that music and counseling programs improved self-esteem and reduced internalized HIV-related stigma in AYA-PHIV. Most participants reported better ART adherence and improved virologic outcomes.

Keywords

Trotzer's concept of group counseling music therapy mental health adherence HIV

Introduction

In 2023, there were ∼1.01 million adolescents aged 15 to 19 years living with HIV globally.¹ Meanwhile, ∼25 643 adolescents and young adults living with HIV (AYA-PHIV) in Thailand were aged between 15 and 24 years.² Adolescence is a critical period marked by cognitive and emotional changes, increasing vulnerability to mental health challenges, particularly for those living with HIV. Mental health conditions such as depression and anxiety, compounded by social stigma, can significantly impact antiretroviral treatment (ART) adherence and overall well-being of AYA-PHIV.^3–5

AYA-PHIV are at higher risk for depression, anxiety, and other mental health problems due to some factors, including social support and stigma.⁶ HIV-related stigma could also have a negative impact on their mental health, social support, and ART adherence,⁷ which placed them at risk for developing HIV drug-resistance and unfavorable treatment outcomes.⁸

Several interventions, including education programs, counseling, and enhancing social support from peers, family members, and healthcare professionals, were reported as associated with improving adherence to ART in AYA-PHIV.^9,10 Effects of stigma reduction interventions,¹¹ mental health screening,¹² and self-esteem enhancing interventions¹³ on ART adherence were also reported.¹⁴ Addressing psychosocial stressors among Thai AYA-PHIV while providing HIV care is critical (Chantaratin et al, 2022).¹⁵

According to a systematic review and metaanalysis, music therapy was associated with a significant improvement in ART adherence and a reduction in depression and anxiety symptoms in people living with HIV/AIDS. Music therapy was proposed as a possible adjunct to standard ART care in promoting better treatment adherence and mental health outcomes for people living with HIV/AIDS.¹⁶ However, its effects in AYA-PHIV remain obscure. This study aimed to evaluate the efficacy of an integrated music and counseling intervention on mental health and ART adherence of AYA-PHIV with virologic nonsuppression.

Methods

Study Design and Participants

This randomized controlled trial was conducted in Chiang Mai, Thailand, from July 2023 to June 2024. AYA-PHIV were recruited from HIV clinics in Chiang Mai and neighboring provinces. Potential participants were informed about the study and invited to participate. Those who were interested in joining underwent an informed consent process. Inclusion criteria were (1) aged between 15 and 24 years, (2) documented perinatal HIV infection, (3) receiving dolutegravir-based ART regimens at enrollment, (4) had HIV-1 RNA > 50 copies/mL within 1-year prior to enrollment, and (5) having mild to moderate mental health scores (defined as Depression Anxiety Stress Scale (DASS-21) test results from 5 to 10 for depression, from 4 to 7 for Anxiety, from 8 to 12 for Stress, or self-esteem scores ≤20, or Thai Internalized HIV-related Stigma Scale (Thai-IHSS) scores ≥10). Exclusion criteria were (1) documented or suspected major psychological problems as determined by the investigator based on available medical records, (2) having severe or higher mental health scores at screening (the DASS-21 test score ≥ 11 for depression, ≥ 8 for anxiety, or ≥ 13 for stress), (3) having other clinically significant diseases that could interfere with the study participation, and (4) having any other condition or social circumstances that would make study participation unsafe, or complicate study outcome interpretation. Stratified random sampling was used in this study. At enrollment, participants were stratified into four study arms based on age, sex, HIV-1 RNA levels, and mental health scores (Figure 1).

Figure 1.

CONSORT 2025 flow diagram, flow diagram of the progress through the phases of a randomized trial of four groups (ie, enrolment, intervention allocation, follow-up, and data analysis).

Figure 2.

Q–Q plot of standardized residuals from the linear mixed-effects model for depression.

Figure 3.

Q–Q plot of standardized residuals from the linear mixed-effects model for anxiety.

Figure 4.

Q–Q plot of standardized residuals from the linear mixed-effects model for stress.

Sample Size Justification

The sample size for this study was determined based on feasibility considerations and the recommended size for group-based counseling interventions. Previous literature suggests that counseling groups consisting of approximately five to 10 participants are optimal for facilitating effective interaction, participation, and peer support within the group process.

Given the limited size of the target population and the feasibility of implementing multiple counseling groups within the study setting, a total of 24 participants were recruited and randomly assigned to four groups, with six participants per group. This group size falls within the recommended range for effective group counseling interventions.

Study Interventions

The study is composed of three interventions and a control group. Participants were assigned to the integrated music and counseling program (Arm 1), the music program (Arm 2), the counseling program (Arm 3), or the control group without intervention (Arm 4) by study nurses using the REDCap program on the central computer with limited access.

The intervention consisted of 10 weekly 90-minute sessions facilitated by two trained counseling psychologists. The Integrated Music and Counseling Program (Arm 1) combined Trotzer's concept-based group counseling¹⁷ with structured music interventions, including guided listening and chorus singing, guitar and keyboard playing, and improvisation activities. The music therapy sessions were facilitated by a music professional (a university music instructor). The Music Program (Arm 2) involved passive and active engagement with music, focusing on emotional expression, relaxation, and social connection. The counseling program (Arm 3) followed Trotzer's approach group counseling framework, emphasizing trust-building, self-reflection, and peer support (Niyonzima et al, 2019).¹⁷ Guided discussions and group-sharing activities to address psychological stressors and stigma were conducted.

Study Intervention

Arm 1: Integrated music and counseling program. To create the integrated music and counseling program as well as select the music pieces that were implemented in this study, the investigator evaluated the study participants’ characteristics, preferences, and interests, as well as the counseling group situations. The gathered participants’ characteristics, preferences, and interest data were analyzed. The program was created to be appropriate with the group of participants in Arm 1. The created program was reviewed and approved by the music therapy and Counseling experts.

The integrated music and counseling program was what the intervention participants were attending (90 min/visit, once a week, 10 visits), as shown in Table 1.

The music pieces that were implemented in this study were considered based on the study participants’ characteristics, preferences, and interests, as well as the counseling group situations with expert approval. Participants were asked to rank and select from groups of 3–5 musical pieces. The music pieces that were also considered for use in this study are categorized in Table 2.

Arm 2: Music program. To select the music therapy program as well as the music pieces that were implemented in this study, the investigator evaluated the study participants’ characteristics, preferences, and interests as well as the counseling group situations. The data was analyzed. The program was selected to be appropriate with the group of participants in intervention arm 2. The program created was reviewed and approved by the counseling expert and tested for construct validity in adolescents and young adults living with HIV before being used with research participants.

The music program is what the intervention participants attended (90 min/visit, once a week, 10 visits), as shown in Table 3.

Arm 3: Counseling program. To create the counseling program that was implemented in this study, the investigator reviewed the handbooks and relevant research on group counseling. The data was analyzed. The program was created to be appropriate with the group of participants in intervention arm 3.

Table 1.

The Integrated Music and Counseling Program.

Study Visit		The Integrated Music and Counseling Program
The program is an integrated music and Trotzer's group counseling program.¹⁸ Trotzer's concept is divided into five stages as follows: (1) Initial Stage or Security stage: - Acquaintance with Each Other - Interpersonal warm-up - Setting boundaries - Building trust (2) The acceptance stage: - Personal sharing - Giving Feedback - Building cohesiveness and closeness and Groupness - Learning to accept self. - Learning to accept others. (3) The responsibility stage: - Individual assessment - Recognizing ownership - Experiencing and building responsibility - Giving respect - Doing a fair share (4) Work stage: - Problem solving - Mobilizing group resources - Reality testing (5) The closing stage or terminating stage*: - Giving support - Affirming and confirming growth - Saying goodbye
The participant-centered music activities were integrated into each session as appropriate. The investigator selected the music pieces/activities. The pieces of music that were used in each session were selected based on the participants’ characteristics, preferences, and interest with expert approval. The ISO principle was implemented in the integrated music and counseling program. The ISO principle, also known as the “individual-social” or “I-Thou” principle, is a concept in group counseling that emphasizes the importance of balancing individual expression and interaction with others in a group setting. The idea behind the ISO principle is that group therapy sessions should provide a space for both individual and collective exploration and growth. In practice, the ISO principle can be achieved by allowing each member of the group to share their thoughts, feelings, and experiences, while also encouraging group members to listen actively and respond in a supportive and nonjudgmental way. This can help foster a sense of community and connection within the group and can provide individuals with the opportunity to gain new insights and perspectives from their peers. According to the ISO principle, music therapy intervention guidelines recommend that people with emotional disorders benefit from listening to music. In accordance with the ISO principle, clients should initially listen to music that reflects their present mood before gradually shifting to music that expresses the desired mental state. The prior study revealed that the movie clip made people melancholy. In comparison to the other groups, the group of participants who listened to the joyful music after listening to the sad music finally reported higher positive effects, higher emotional valence, and lower negative effects. Not all between-group differences, nevertheless, were statistically significant. The study's result was that the current emotional state is affected by the sequence in which songs with various emotional valences are delivered. ¹⁹
Session 1	Purpose: (1) To get people ready to form a group (2) To build relationships and provide members with the opportunity to learn more about each other (3) To create an environment of trust and ownership Activities: Building dyads with listening to popular music Procedures: (1) The group leader welcomed the group members, introduced himself, co-leader, and had each member introduce themselves. (2) The group leader explained the purpose and duration of the group. (3) The group leader introduced the exercise “Building Relationships” (Building Dyads). (4) The group leader asked the members to come together to share a summary of what they had learnt. (5) The group leader asked participants to bring a song that they felt best reflected who they were to be discussed during the next session. Results of the activities: Group members of the group felt less nervous, more at ease, and more acquainted. They asked questions and interacted happily because they were interested about their peers’ stories.
Session 2	Purpose: (1) To allow members to explore themselves in various aspects of their lives. (2) To allow members to reveal more about themselves. (3) To provide participants with the opportunity to explore their own emotions and experiences. To let participants recognize the worth that each individual possesses. (4) To make it possible for every member to recognize and value the contributions of others. (5) To make themselves and other people feel good. Activities: Symbol of Me with Song Analysis (Listening to the song that symbolized each member and analyzing the song's meaning and reasons) Procedures: (1) The group leader provided a summary of the previous period's activities. (2) The activity was introduced by the group leader. “How does the owner of the song feel?” (This is me.) (3) The group leader asked everyone to give a short summary of what they have learned. (4) To introduce the “Self-esteem” activity, the group leader asks participants to list what characteristics they find valuable in each of the following categories. (5) The leader asks participants to provide a summary of what they have learnt. Results of the activities: As a result of the activity, the members had the chance to become more self-aware of themselves. Members feel satisfied, pleased, and aware of their own value. They additionally recognize that everyone has value to both them and others. Furthermore, they support one another and participate in their peers’ enjoyment.
Sessions 3-6	Purpose: (1) Provide participants with the chance to explore their unresolved problems and significant life events to better understand and accept themselves. (2) Allow group members to share their experiences and learn how to handle different life challenges. (3) Allow members to assist one another in resolving issues, as this encourages a feeling of intimacy and acceptance. Activities: Road of life with listening to classical/popular music. Procedures: (1) To introduce the “Road of life” activity, the group leader asked participants to sketch a road that symbolizes their experiences from previous lives. (2) The participants shared their own life stories with the group (3) The group leader asked participants to provide a summary of their learning. Results of the activities: As a result of the activity, participants felt more comfortable sharing their personal stories and were able to see the problems with their own and their friends’ lives. It was encouraged that people overcome the challenges they encounter. Participants felt better about others when they knew they had friends who were always willing to listen and support them when they were struggling.
Session 7	Purpose: (1) Provide participants with the chance to recognize themselves, peers, family, and other supporters around them. (2) Allow group members to keep the mind that they always have support from their important people, especially themselves. Activities: Bedside Diary with listening to classical/popular music. Procedures: (1) Introduce the “Bedside Diary” activity. The group leader asked participants to create a four-page diary for themselves with the details for each page as follows. Page 1 of 4: Cover page filled with the Name they called themselves and the photo relevant to their name. Page 2 of 4: Draw the Picture of their favorite musical instrument. Write about the benefits of this musical instrument and how it is like themselves. Page 3 of 4: Draw a picture of an impressive event that every time they think about it, they will feel happy. Page 4 of 4: Write the following message on each of the 5 hearts: Heart 1: Write an externally visible favor that you find admirable from a friend seated on the left. Heart 2: List three positive characteristics of the friend seated on the right. Heart 3: Write an internal reflection of something you admire that is recognized from the left. Heart 4: Write encouraging words for the friend on the right-hand seat next to you. Heart 5: Write to express gratitude for being your lifelong best friend. (2) The group leader asked participants to provide a summary of their learning. Results of the activities: As a result of the activity, Participants reflected on themselves and their families, recalled their own important moments, recognized their own positive characteristics that their friends appreciated, and got support from both them and their group members. Purpose: (1) To get everyone warmed up (2) To guide participants to have trust in others (3) To allow participants to learn and care for each other. (4) To learn choral with the team to promote self-esteem Procedures: (1) The leader summarizes what has been done together in the past. (2) The group leader introduced the exercise “Trust Walk” by having the members line up in a single row and cover the eyes of the person in front of them, like a blind person. Then, the members chose a group leader to lead them down the corridor of the building and to the room where the group activities were held. (3) The group leader asked the members to summarize what they had learned. (4) The music teacher taught the participants to practice choral starting from listening to the song and analyzing the lyrics. Results of the activities: As a result of the activity, the members were excited about the activity because they had to walk outside the room and gently held hands to close their friends’ eyes and walk carefully. The participants who were blind received good care, felt safe, and trust each other more. The participants were excited about the choral learning.
Session 8	Activities: Trust Walk and Choral Learning
Session 9	Purpose: (1) To learn choral with the team to promote self-esteem Activities: Choral Learning Procedures: (1) The leader summarizes what has been done together in the past. (2) The leader asked the members to summarize what they had learned. (3) The music teacher taught the participants to practice choral starting from listening to the song and analyzing the lyrics. Results of the activities: As a result of the activity, the participants were proud of themselves for their choral learning.
Session 10	Purpose: (1) Members had the opportunity to explore experiences of good relationships between themselves and others. (2) Members received encouragement from the group to continue living. Activities: Choral Learning and Impression Go Round Procedures: (1) The group leader summarized what they had done together in the past period. (2) The group leader presented the exercise “Impression Go Round”. (3) The group members were asked to share their impressions of the good times they had together on the last opportunity to join the group. (4) The group leader asked the group members to summarize what they have learned from joining the group this time and say goodbye. Results of the activities: As a result of the activity, the participants were proud of themselves for their learning. As they had an opportunity to share their thoughts about their work, the group environment was cheerful. In addition to communicating their appreciation for the positive things their friends had done for them—such as listening to each other in the group, accepting each other, and giving them suggestions for problem-solving, the group also came to know and understand one another better.

Note: The group counseling stages in the above table are guidelines used to conduct group counseling in this study; however, group counseling in participants living with HIV is flexible according to the group's situation. Each stage may occur later or earlier than planned, or some stages may not occur as planned, such as in the responsibility and work stages.

Table 2.

Music Pieces Were Implemented in This Study.

No	Category	Title/Link	Artist
1	Baroque Period	Pachelbel – Canon in d https://youtu.be/Pppexz-KKig	Johann Pachelbe 1653-1706 Nippon New Philharmonic Orchestra/Conductor: Takashi Ueno. Arr/T.Ueno.
		https://youtu.be/JvNQLJ1_HQ0	Pachelbel Canon in D Major—the original and best version
		https://youtu.be/es_3F3TLJS0	Wedding String Quartet—Canon in D (Best Version) (Johann Pachelbel)
		https://youtu.be/wny1ojJli8k	Canon in D (Pachelbel) Piano Cover by Ryan Jones
2	Classical period	Piano Sonata in A Major, K. 331: III. Alla Turca https://youtu.be/quxTnEEETbo	Elisso Bolkvadze, Wolfgang Amadeus Mozart
3	Romantic period	CHOPIN – Nocturne No 2 in E Flat Major Op 9 No 2 https://youtu.be/TqyLnMa3DJw	Frank Levy
		Chopin Nocturne E Flat Major Op.9 No.2 https://youtu.be/tV5U8kVYS88	Valentina Lisitsa
		Für Elise https://www.youtube.com/watch?v = yGLpHxdsLtM https://www.youtube.com/watch?v = R3Y6X0dCk-o	Beethoven
		https://www.youtube.com/watch?v = Wur0Sii4T-U	Beethoven
4	Relaxing music	https://www.youtube.com/watch?v = hlWiI4xVXKY	Beautiful Relaxing Music • Peaceful Piano Music & Guitar Music \| Sunny Mornings by @PederBHellandMusic
		https://www.youtube.com/watch?v = y1dbbrfekAM&list = RDy1dbbrfekAM&index=1	Chopin – Relaxing Classical
		https://www.youtube.com/watch?v = Rb0UmrCXxVA&list = RDy1dbbrfekAM&index =	Mozart – Relaxing Classical
	Participants’ favorite songs: Selected songs to perform		Do you see that sky? Pick up the sun

Table 3.

The Music Program.

Study Visit	The Music Program
Session 1	Acquaintance with Each Other: Rapport building with the passive music activities. (Listening to Classic and/or popular music)
Session 2	Music Listening and Feelings Expression with passive music activities. (Listening to Classic and/or popular music)
Session 3	Body Movement and Muscle Relaxation with passive music activities. (Listening to Classic and/or popular music)
Session 4	Singing with active music activities. (Listening to Popular music and singing)
Session 5	Improvisation with active music activities. (Listening to Popular music and improvisational playing)
Session 6	Stories Creating with Music Listening with the active music activities. (Listening to Classical and/or Popular music and creating stories)
Session 7	Musical Imagery with active music activities. (Listening to Classical and/or Popular music and creating stories)
Session 8	Drawing in conjunction with Music Listening with the active music activities. (Listening to Classical and/or Popular music and creating stories)
Session 9	Psychological Drama with active music activities. (Listening to Classical and/or Popular music and creating stories)
Session 10	Farewell and Feelings Sharing

The counseling program was what the intervention participants attended (90 min/visit, once a week, 10 visits), as shown in Table 1, but not including the music.

Recruitment and Training of Study Psychological Counselor Staff

The investigator identified the target of two psychological counselor study staff, including the investigator, to be recruited and trained to lead the Intervention sessions. The psychological counselor staff had some experience in mental health and group counseling. Identified psychological counselors were assigned at random to lead the Intervention sessions and received the necessary training.

Training

Designated expert trainers provided training to all study psychological counselors staff to manage safety concerns among participants. The study psychological counselor study staff designated to lead the Music group sessions was trained separately from the psychological counselor staff. Separate training minimized the risk of contamination across treatment arms. The study psychological counselor staff, who were designated to lead the intervention group sessions practiced delivering the sessions during the training.

Study Psychological Counselor Staff Supervision

During psychological counselor supervision sessions, the study staff received ongoing supervision from the expert trainer. The expert trainer used the Zoom application to review any concerns raised during the sessions and problem-solve best practices.

Data Collection and Outcome Measures

At enrolment, baseline demographic data and ART adherence were collected using self-reported questionnaires. CD4 and HIV-1 RNA results were extracted from the National AIDS Program database. Mental health metrics included the Depression Anxiety Stress Scale (DASS-21), Rosenberg Self-Esteem Scale (RSES-10), and Thai Internalized HIV-related Stigma Scale (Thai-IHSS).

The Depression Anxiety Stress Scale (DASS-21) is a self-report scale for measuring depressive, anxious, and stressed emotional states.²⁰ The test demonstrated internal consistency and validity with the Cronbach's alpha coefficient of 0.82. The Thai version was adapted across cultures and validated. The Rosenberg Self-Esteem Scale (RSES-10) is a commonly utilized instrument that has been evaluated for reliability and validity in a variety of contexts.²¹ The revised Thai RSES showed equivalent reliability and improved construct validity compared to the original. The Thai Internalized HIV-related Stigma Scale (Thai-IHSS) is a tool developed by Thai researchers (Uthis et al, 2023).²² It is composed of 22 items which were considered as meaningful and suitable for the Thai context, including negative thoughts toward self, anticipated negative thoughts, effects of negative thoughts toward self, and effects of negative thoughts toward family and access to healthcare services. A 4-point rating scale (strongly disagree, disagree, agree, and strongly agree) was used; higher scores indicated higher HIV-related stigma.

Biospecimen Collection

We collected hair samples to determine DTG levels as markers of ART adherence at the study entry, week 10, and week 24. Hair sample was collected from the vortex posterior (back of the head) by study nurses in the clinic. At least 20 hair fibers, each 1 cm long, were obtained from the portion of the hair strands closest to the scalp. Participants were advised to avoid hair treatment, wash their hair, and have it dried before coming to the clinic. The hair samples were placed in a piece of foil, labeled, stored in a sealed plastic Ziplock bag with desiccant, and kept at room temperature with light protection prior to shipment.

This randomized controlled trial was registered at the Thai Clinical Trials Registry https://www.thaiclinicaltrials.org/. (Identifier: TCTR20241231002) on December 31, 2024.

The reporting of this randomized controlled trial follows the Consolidated Standards of Reporting Trials (CONSORT) guidelines.²³ The completed CONSORT checklist is provided as the Supplemental File.

Data Analysis

Demographic information was analyzed by using descriptive statistics (frequency distribution), mean (standard deviation, SD), or median (interquartile range, IQR) as appropriate. Intention-to-treat analysis was used for participant outcome measure analysis. The number of session attendance of all randomly selected participants was counted. Overall and within the Integrated Music and Counselling program Intervention and Control Arms, baseline variables (gender, age, severity of mental health ratings, CD4, and viral load suppression status) were compared between participants with and without VL suppression. The lost-to-follow-up rates were compared between arms. The primary outcome analysis focused on the mental health scores, ART adherence scores, and ART concentration in hair samples at week 10. As the data has a nonnormal distribution, nonparametric statistics, including Kruskal-Wallis, Friedman, and Mann-Whitney U-test, were utilized.

Statistical analyses were also performed using linear mixed-effects models (LMMs) to evaluate differences in depression, anxiety, stress, self-esteem, and internalized HIV-related stigma scores, CD4 count, HIV-1 RNA, saliva cortisol, and hair DTG concentration across treatment arms and changes over study visits. Treatment arm, study visit (entry, week 10, and week 24 for depression, anxiety, stress, self-esteem, and internalized HIV-related stigma scores, saliva cortisol, and hair DTG concentration; and entry and week 24 for CD4 count and HIV-1 RNA), and their interaction were specified as fixed effects. Participant-level random effects (random intercepts and random slopes for the arm-by-visit interaction) were included to account for within-subject correlation arising from repeated measurements.

Residuals were assumed to follow a normal distribution. Model parameters were estimated using maximum likelihood, and degrees of freedom for fixed effects were approximated using the Satterthwaite method. Overall effects of fixed factors were assessed using F-tests. Post hoc pairwise comparisons between treatment arms were conducted with Bonferroni adjustment for multiple testing.

Model assumptions were evaluated by inspecting residual diagnostics, including Kolmogorov–Smirnov and Shapiro–Wilk tests, Q–Q plots, and residual-versus-fitted value plots. Statistical significance was defined as a two-sided P-value < .05.

The SPSS Statistics program, version 29, and Jamovi Statistics program, version 2.7,^24–29 were used in data analysis of this study.

Results

Demographic Characteristics

Of the 24 enrolled participants, eight (33%) were male. Their median age was 21.5 years (range 15-24). Four (17%) were university graduates, and 20 (83% were either undergraduate students or had stopped schooling. Eight (33%) were working, 10 (42%) were studying, and six (25%) did neither. Their HIV-1 RNA levels ranged from 116 to 600, 500 copies/mL within the past 12 months before enrollment. There were 16 (42%) with virologic failure (HIV-1 RNA > 1000 copies/mL) (Table 4).

Table 4.

Demographic Characteristics of Study Participants (n = 24).

	Characteristic	Arm 1	Arm 2	Arm 3	Arm 4	Total N (%)
Biological Sex
	Male	2	1	4	1	8 (33)
Age (years)
Median						21 (4)
	15-19	3	3	0	0	6 (25)
	20-24	3	3	6	6	18 (75)
Education
	University graduated	1	1	1	1	4 (17)
	Undergraduate	5	5	5	5	20 (83)
Employment/educational status
	Working	2	2	2	2	8 (33)
	Studying	2	2	3	3	10 (42)
	Neither studying nor working	2	2	1	1	6 (25)
HIV-1 RNA
	</=3 Log₁₀	2	1	2	3	8 (33)
	3-<4 Log₁₀	1	2	1	1	5 (21)
	4-<6 Log₁₀	3	3	2	2	10 (42)
	>/=6 Log₁₀	0	0	1	0	1 (4)
ART Adherence
	Excellent	0	0	1	0	1 (4)
	Very good	1	1	1	1	4 (17)
	Good	3	2	2	2	9 (37)
	Fair	0	3	0	3	6 (25)
	Poor	1	0	2	0	3 (13)
	Very poor	1	0	0	0	1 (4)

Mental Health Outcomes

Depression, anxiety, and stress scores were analyzed using a linear mixed-effects model (LMM), with treatment arm, study visit (entry, week 10, and week 24), and their interaction specified as fixed effects, and participant-specific random effects to account for within-subject correlation arising from repeated measurements. A total of 72 observations from 24 participants were included in the analysis.

The omnibus tests of fixed effects indicated no significant main effect of treatment arm on depression scores (F = 0.168, P = .916), anxiety scores (F = 0.141, P = .933), and stress scores (F = 0.0902, P = .964) and no significant main effect of study visit (F = 0.390, P = .685 for depression, F = 5.353, P = .014 for anxiety, and F = 2.5652, P = .110 for stress). In addition, the treatment Arm × Study Visit interaction was not statistically significant (F = 1.811, P = .183 for depression, F = 1.373, P = .299 for anxiety, and F = 1.0106, P = .475 for stress), suggesting that the trajectories of depression anxiety, and stress scores over time did not differ significantly across treatment arms (Table 5).

Table 5.

Fixed Effects Omnibus Tests for Depression, Anxiety, and Stress.

	F	df	df (res)	P
Depression
Arm	0.168	3	10.4	.916
Study Visit	0.390	2	12.9	.685
Arm × Study Visit	1.811	6	11.5	.183
Anxiety
Arm	0.141	3	9.20	.933
Study Visit	5.353	2	18.91	.014
Arm × Study Visit	1.373	6	12.42	.299
Stress
Arm	0.0902	3	10.55	.964
Study Visit	2.5652	2	15.03	.110
Arm × Study Visit	1.0106	6	8.84	.475

Parameter estimates further indicated that mean depression, anxiety, and stress scores at baseline were ∼2.13, 2.78, and 3.06, respectively, with no statistically significant differences between treatment arms relative to the reference group. No significant changes in depression, anxiety, and stress scores were observed at week 10 or week 24 compared with baseline. Post hoc pairwise comparisons between treatment arms, adjusted using the Bonferroni method, revealed no statistically significant differences across any pairwise contrasts.

The random-effects component showed substantial between-subject variability, with a high intraclass correlation coefficient (ICC ≈ 0.97, 0.799, and 0.782), indicating strong within-subject consistency of depression, anxiety, and stress scores over time (Table 6).

Table 6.

Random Components for Depression, Anxiety, and Stress.

Groups	Name	Variance	SD	ICC
Depression
Participant ID	(Intercept)	2.1234	1.457	0.969
	Arm1:Study VisitEntry	0.2194	0.468
	Arm2:Study VisitEntry	7.4484	2.729
	Arm3:Study VisitEntry	3.2775	1.810
	Arm4:Study VisitEntry	0.1565	0.396
	Arm1:Study VisitWeek 10	3.8246	1.956
	Arm2:Study VisitWeek 10	9.7680	3.125
	Arm3:Study VisitWeek 10	3.5196	1.876
	Arm4:Study VisitWeek 10	0.4200	0.648
	Arm1:Study VisitWeek 24	9.5290	3.087
	Arm2:Study VisitWeek 24	6.3579	2.521
	Arm3:Study VisitWeek 24	2.9206	1.709
	Arm4:Study VisitWeek 24	0.1853	0.430
Residual		0.0687	0.262
Note: Number of Obs: 72, Number of groups: Participant ID 24
Anxiety
Participant ID	(Intercept)	2.289	1.513	0.799
	Arm1:Study VisitEntry	1.497	1.223
	Arm2:Study VisitEntry	0.823	0.907
	Arm3:Study VisitEntry	1.706	1.306
	Arm4:Study VisitEntry	4.446	2.109
	Arm1:Study VisitWeek 10	5.053	2.248
	Arm2:Study VisitWeek 10	1.386	1.177
	Arm3:Study VisitWeek 10	5.922	2.434
	Arm4:Study VisitWeek 10	2.233	1.494
	Arm1:Study VisitWeek 24	2.436	1.561
	Arm2:Study VisitWeek 24	3.322	1.823
	Arm3:Study VisitWeek 24	2.176	1.475
	Arm4:Study VisitWeek 24	0.430	0.656
Residual		0.576	0.759
Note: Number of Obs: 72, Number of groups: Participant ID 24
Stress
Participant ID	(Intercept)	3.361	1.833	0.782
	Arm1:Study VisitEntry	5.014	2.239
	Arm2:Study VisitEntry	2.154	1.468
	Arm3:Study VisitEntry	2.802	1.674
	Arm4:Study VisitEntry	2.026	1.423
	Arm1:Study VisitWeek 10	9.794	3.130
	Arm2:Study VisitWeek 10	8.970	2.995
	Arm3:Study VisitWeek 10	5.084	2.255
	Arm4:Study VisitWeek 10	1.553	1.246
	Arm1:Study VisitWeek 24	1.600	1.265
	Arm2:Study VisitWeek 24	8.913	2.986
	Arm3:Study VisitWeek 24	1.825	1.351
	Arm4:Study VisitWeek 24	3.486	1.867
Residual		0.935	0.967

Note: Number of Obs: 72, Number of groups: Participant ID 24.

The Q–Q plot of standardized residuals showed that most observations aligned closely with the reference line in the central quantiles, indicating approximate normality of residuals overall. Deviations from the reference line were observed in both tails, with a few points departing noticeably from linearity, suggesting the presence of mild tail departures and potential outliers (Figures 2–4).

Notably, the model exhibited convergence issues, with an overparameterized random-effects structure relative to the number of observations. This may limit the stability and identifiability of the variance component estimates. Therefore, the results should be interpreted with caution, and sensitivity analyses using a more parsimonious random-effects structure are warranted.

At enrollment, their self-esteem, IHSS, depression, anxiety, and stress scores were Similar across the four arms. AYA-PHIV in Arm 1 (integrated music and counseling program) had significantly increased self-esteem scores at week 10 and week 24. Increasing trends were observed for self-esteem scores of participants in Arm 2 (music program) and Arm 3 (counseling program), although the changes were not statistically significant. There was no change in the self-esteem scores among AYA-PHIV in Arm 4 (the control group).

Self-esteem scores (continuous outcome) were also analyzed using linear mixed-effects models (LMMs) to examine the effects of study arm, study visit (entry, week 10, and week 24), and their interaction. Fixed effects included arm, visit, and arm-by-visit interaction. Random effects were specified at the participant level (Participant ID) with a random intercept and arm-by-visit slopes to account for within-subject correlation of repeated measures. Models were fitted using maximum likelihood estimation with the nloptwrap optimizer, and degrees of freedom were computed using the Satterthwaite approximation. Wald 95% confidence intervals were reported (Tables 7 and 8).

Table 7.

Fixed Effects Omnibus Tests for Self-Esteem.

	F	df	df (res)	P
Arm	9.92	3	10.5	.002
Study Visit	14.33	2	10.9	<.001
Arm × Study Visit	9.46	6	11.5	<.001

Table 8.

Random Components for Self-Esteem.

Groups	Name	Variance	SD	ICC
Participant ID	(Intercept)	2.861	1.692	0.939
	Arm1:Study VisitEntry	1.900	1.379
	Arm2:Study VisitEntry	3.187	1.785
	Arm3:Study VisitEntry	34.245	5.852
	Arm4:Study VisitEntry	2.727	1.651
	Arm1:Study VisitWeek 10	3.835	1.958
	Arm2:Study VisitWeek 10	22.614	4.755
	Arm3:Study VisitWeek 10	1.019	1.009
	Arm4:Study VisitWeek 10	0.800	0.894
	Arm1:Study VisitWeek 24	7.510	2.741
	Arm2:Study VisitWeek 24	7.221	2.687
	Arm3:Study VisitWeek 24	10.728	3.275
	Arm4:Study VisitWeek 24	3.425	1.851
Residual		0.187	0.432

Note: Number of Obs: 72, Number of groups: Participant ID 24.

Residual normality was assessed using Kolmogorov–Smirnov and Shapiro–Wilk tests and Q–Q plots. Given the complexity of the random-effects structure relative to the sample size, model results were interpreted with caution, and sensitivity analyses using simplified random-effects structures were conducted to assess the robustness of findings (Figure 5).

Figure 5.

Q–Q plot of standardized residuals from the linear mixed-effects model for self-esteem.

A significantly decreased internalized HIV-related stigma was observed among participants in Arm 1 and Arm 3 at both week 10 and week 24 of study participation. No change was observed in Arm 2 AYA-PHIV and the control arm. No significant change in depression, anxiety, and stress among AYA-PHIV in any intervention arms. In addition, we found participants in Arm 4 (control group) had significantly increased depression scores at week 24 (Table 9).

Table 9.

Self-Esteem and Internalized HIV-Related Stigma in Study Participants.

Study Visit	Arm 1	Arm 2	Arm 3	Arm 4	P Value
Self-Esteem scores [Median (IQR)]
Entry(Week 0)	17.0(14.8-18.0)	17.5(16.8-19.5)	13.0(9.0-19.0)	16.5(15.8-19.0)	.27
Week 10	25.0(22.3-27.0)	25.5(20.8-27.0)	19.0(18.8-22.5)	16.5(14.5-18.5)	. 01
Week 24	21.0(19.0-24.8)	20.5(19.0-24.5)	18.0(16.0-22.8)	15.5(14.8-18.5)	.04
P value	0.01	0.13	0.21	0.69
Internalized HIV-related Stigma scores [Median (IQR)]
Entry(Week 0)	47.5 (44.0-52.8)	31.5 (24.5-41.0)	35.5 (30.5-43.8)	38.5 (23.5-45.8)	.09
Week 10	24.0 (21.7-27.5)	31.5 (23.5-36.8)	28.0 (22.8-37.5)	36.0 (22.8-48.8)	.27
Week 24	24.0 (22.0-36.0)	28.0 (25.0-35.0)	29.0 (23.5-35.5)	35.0 (19.8)	.64
P value	.01	.58	.01	1

Self-esteem scores (continuous outcome) were also analyzed using LMMs to examine the effects of study arm, study visit (entry, week 10, and week 24), and their interaction. Fixed effects included arm, visit, and arm-by-visit interaction. Random effects were specified at the participant level (Participant ID) with a random intercept and arm-by-visit slopes to account for within-subject correlation of repeated measures. Models were fitted using maximum likelihood estimation with the nloptwrap optimizer, and degrees of freedom were computed using the Satterthwaite approximation. Wald 95% confidence intervals were reported (Tables 10 and 11).

Table 10.

Fixed Effects Omnibus Tests for Self-Esteem.

	F	df	df (res)	P
Arm	9.92	3	10.5	.002
Study Visit	14.33	2	10.9	<.001
Arm × Study Visit	9.46	6	11.5	<.001

Table 11.

Random Components for Self-Esteem.

Groups	Name	Variance	SD	ICC
Participant ID	(Intercept)	2.861	1.692	0.939
	Arm1:Study VisitEntry	1.900	1.379
	Arm2:Study VisitEntry	3.187	1.785
	Arm3:Study VisitEntry	34.245	5.852
	Arm4:Study VisitEntry	2.727	1.651
	Arm1:Study VisitWeek 10	3.835	1.958
	Arm2:Study VisitWeek 10	22.614	4.755
	Arm3:Study VisitWeek 10	1.019	1.009
	Arm4:Study VisitWeek 10	0.800	0.894
	Arm1:Study VisitWeek 24	7.510	2.741
	Arm2:Study VisitWeek 24	7.221	2.687
	Arm3:Study VisitWeek 24	10.728	3.275
	Arm4:Study VisitWeek 24	3.425	1.851
Residual		0.187	0.432

Note: Number of Obs: 72, Number of groups: Participant ID 24.

Figure 6.

Q–Q plot of standardized residuals from the linear mixed-effects model for self-esteem.

The LMM showed no significant overall differences in HIV-1-related stigma between study arms, F(3, 9.96) = 0.117, P = .948. In contrast, stigma scores changed significantly over time, F(2, 16.35) = 66.623, P < .001, and a significant Arm × Study Visit interaction was observed, F(6, 11.83) = 24.855, P < .001, indicating differential stigma trajectories across study arms during follow-up (Table 12). Random-effects estimates demonstrated substantial between-participant variability in baseline stigma (intercept variance = 49.10; ICC = 0.965), with additional heterogeneity in individual change patterns across arms and study visits (Table 13).

Table 12.

Fixed Effects Omnibus Tests for HIV-1-Related Stigma.

	F	df	df (res)	P
Arm	0.117	3	9.96	.948
Study Visit	66.623	2	16.35	<.001
Arm × Study Visit	24.855	6	11.83	<.001

Table 13.

Random Components for HIV-1-Related Stigma.

Groups	Name	Variance	SD	ICC
Participant ID	(Intercept)	49.10	7.01	0.965
	Arm1:Study VisitEntry	47.48	6.89
	Arm2:Study VisitEntry	70.94	8.42
	Arm3:Study VisitEntry	12.29	3.51
	Arm4:Study VisitEntry	45.81	6.77
	Arm1:Study VisitWeek 10	14.93	3.86
	Arm2:Study VisitWeek 10	77.12	8.78
	Arm3:Study VisitWeek 10	7.22	2.69
	Arm4:Study VisitWeek 10	78.76	8.87
	Arm1:Study VisitWeek 24	74.98	8.66
	Arm2:Study VisitWeek 24	10.83	3.29
	Arm3:Study VisitWeek 24	7.39	2.72
	Arm4:Study VisitWeek 24	18.73	4.33
Residual		1.79	1.34

Note: Number of Obs: 72, Number of groups: Participant ID 24.

Model diagnostics supported the adequacy of distributional assumptions; inspection of the normal Q–Q plot of standardized residuals showed close alignment with the theoretical reference line, with only minor deviations at the upper tail, indicating approximate normality of residuals (Figure 7).

Figure 7.

Q–Q plot of standardized residuals from the linear mixed-effects model for HIV-1-related stigma.

ART Adherence Outcomes

Improved self-report ART adherence was seen during and after the intervention at week 10 and 24, respectively, among AYA-PHIV who participated in Arms 1, 2, and 3, although it did not reach a significant level. The DTG hair level of our study participants ranged from 200 to 4143 pg/mL. While there is no therapeutic range for DTG hair level, we found six participants (one, three, and two participants in Arms 1, 3, and 4, respectively) with increased DTG hair levels by 1.5 times or more when compared to baseline at the study enrollment. However, the decrease in median HIV-1 RNA level was observed in all study arms (Table 14). At week 24, 16 (67%) achieved virologic suppression.

Table 14.

HIV-1 RNA Level at the Study Entry and Week 24 of Participants in Each Arm.

Visit	HIV-1 RNA (copies/mL) [Median (IQR)]
Visit	Arm 1	Arm 2	Arm 3	Arm 4	P Value
Entry(Week 0)	16 050(397-14 4228)	13 250(5374-154 650)	23 125(361-600 500)	617(116-33 200)	0.86
Week 24	1205(20-37 436)	1205(20-9684)	65(20-4579)	20(20-46)	0.67
P value	<.05	<.05	<.05	<.05

LMMs showed no significant between-arm differences in HIV-1 RNA levels (F = 0.757, P = .612) and no significant main effect of study visit (F = 2.973, P = .142). The Arm × Study Visit interaction was also not statistically significant (F = 1.783, P = .359), indicating comparable viral load trajectories across study arms over follow-up (Table 15).

Table 15.

Fixed Effects Omnibus Tests for HIV-1 RNA.

	F	df	df (res)	P
Arm	0.757	3	2.00	.612
Study Visit	2.973	2	4.93	.142
Arm × Study Visit	1.783	2	2.00	.359

Fixed-effect estimates were not statistically different from zero, with wide confidence intervals, consistent with substantial between-participant heterogeneity (ICC ≈ 1.00), suggesting that interindividual variability dominated within-person change over time (Table 16).

Table 16.

Random Components for HIV-1 RNA.

Groups	Name	Variance	SD	ICC
Participant ID	(Intercept)	1.74 × 10⁹	41 680	1.000
	Arm1:Study VisitEntry	2.14 × 10⁹	46 280
	Arm2:Study VisitEntry	7.91 × 10⁹	88 913
	Arm3:Study VisitEntry	1.75 × 10¹¹	418 812
	Arm4:Study VisitEntry	2.06 × 10⁹	45 388
	Arm1:Study VisitWeek 10	1.98 × 10⁹	44548
	Arm2:Study VisitWeek 10	2.62 × 10⁹	51 185
	Arm3:Study VisitWeek 10	3.36 × 10⁹	57 970
	Arm4:Study VisitWeek 10	6.04 × 10⁹	77 716
	Arm1:Study VisitWeek 24	3.35 × 10⁹	57 836
	Arm2:Study VisitWeek 24	2.97 × 10⁸	17 248
	Arm3:Study VisitWeek 24	1.74 × 10⁹	41 676
	Arm4:Study VisitWeek 24	2.72 × 10⁹	52 157
Residual		395 610	629

Note: Number of Obs: 30, Number of groups: Participant ID 24.

Note: Random variances cannot be computed. Some variance components equal zero.

Using LMMs with study arm, visit (entry and week 24), and their interaction as fixed effects and participant-level random effects, we found no statistically significant differences in HIV-1 RNA (copies/mL) across study arms (F = 0.757, P = .612), no main effect of visit (F = 2.973, P = .142), and no Arm × Visit interaction (F = 1.783, P = .359).

None of the fixed-effect parameter estimates reached statistical significance (P < .05). Substantial between-participant variability was observed, as indicated by large random-effects variances. Visual inspection of the Q–Q plot suggested deviations from normality in the tails of the residual distribution, consistent with the Shapiro–Wilk test indicating potential departures from normality (Figure 8).

Figure 8.

Q–Q plot of standardized residuals from the linear mixed-effects model for HIV-1 RNA.

Discussion

This study demonstrated the benefits of the integrated music and counseling program intervention on self-esteem, internalized HIV-related stigma, and ART adherence among AYA-PHIV.

The effect on self-esteem was in line with previous research using music therapy intervention in participants with a diagnosis of “adjustment reaction to adolescence,” in which an increase in self-esteem was reported.³⁰ In this study we used different music activities, which could affect mental health of AYA-PHIV like another previous study documented increased self-expression and self-esteem from musical performance.³¹ Music allows participants to express their emotions and their personalities, which can lead to better mental health and self-esteem. This was documented in a study of a woman with severe bereavement who struggled to define her identity. She improved after music therapy.³² Engaging in a variety of musical activities can boost one's confidence and sense of self, which is noticeable in a therapeutic setting.³³

Low self-esteem among AYA-PHIV put them at more vulnerable to internalized stigma related to HIV/AIDS, which results from the self-deprecating effects of comparing oneself to stigmatizing social norms about HIV.³⁰ Our result demonstrated that the integrated music and counseling and the counseling programs were associated with significantly decreased internalized HIV-related stigma in AYA-PHIV. The results are in line with the previous research addressing the benefits of other elements of creativity and HIV, including applied sculpting, drawing, song, and dance,^34–38 while empowerment, health promotion, and HIV awareness were more powerful in tackling internalized stigma. More strategies for decreasing internalized stigma associated with HIV were mentioned in a recent systematic review.³⁹

In addition, the LMM suggested that self-esteem changed over time and that trajectories differed across study arms. However, model nonconvergence and identifiability issues arising from an overparameterized random-effects structure relative to the sample size warrant cautious interpretation of the quantitative estimates. Sensitivity analyses using simplified random-effects structures (reported in the Supplemental Material) yielded directionally consistent findings, supporting the qualitative conclusions regarding differential temporal patterns of self-esteem across arms.

No significant difference was observed in the depression, anxiety, and stress levels among our study participants after intervention. This finding was inconsistent with previous metaanalyses that documented improved anxiety and depressive symptoms.⁴⁰ This might be due to the low frequency of anxiety and depression at enrollment and the short follow-up period. Moreover, we could identify other psychosocial stressors during the study period, including breaking up with a boyfriend, failing the university entrance exam, being hospitalized for acute bronchitis, retiring from the university, and encountering physical and psychological abuse, which could explain worsening mental health scores.

Although the residual distribution appears approximately normal in the central range, the observed tail deviations may reflect mild nonnormality or the influence of outlying observations. In the context of a relatively small sample size and a complex mixed-effects structure, such departures could affect the precision of parameter estimates and confidence intervals. Nevertheless, the deviations were not severe, and the normality assumption was considered reasonably met for inferential purposes. These findings support the primary analyses while underscoring the importance of complementary sensitivity analyses and cautious interpretation.

LMMs accounting for within-participant repeated measures indicated temporal changes in HIV-1-related stigma, with model-estimated means showing improvement from baseline across follow-up visits. The magnitude and pattern of change varied by study arm, suggesting context-specific effects of the intervention or care pathways on stigma-related outcomes.

For young people, group intervention seems to be an effective method to promote mental health and general health. Like a study in South Africa using family-based psychosocial intervention in early adolescents living with HIV, improvement in mental health, behaviors, HIV treatment knowledge, stigma, communication, and ART adherence were noted.⁴¹

Reductions in HIV-1-related stigma over time carry important clinical and public health implications, given the well-established associations between stigma, care engagement, disclosure, and treatment adherence. The observed group-specific trajectories highlight the need for contextually tailored stigma-reduction strategies rather than uniform approaches. Such targeted interventions may enhance psychosocial well-being and support sustained engagement in HIV care.

Improved self-report ART adherence and HIV-1 RNA levels were observed during and after the intervention at week 10 and week 24, respectively, among AYA-PHIV in all three intervention arms. The result was consistent with a previous study in Thai AYA-PHIV of similar age using an empowerment intervention, which reported increased ART adherence.⁴² They did not have HIV-1 RNA outcome, while we did. Nevertheless, improved ART adherence following group therapy in another study among adolescents with HIV was reported.⁴⁴ Unlike the report from a South African study, which reported no significant effect of group counseling on the treatment outcomes.⁴¹

To our surprise, a significantly decreased HIV-1 RNA level was also observed in the control arm when compared to baseline. Study participation could be an explanation for improving adherence. AYA-PHIV might improve their ART adherence due to awareness of being observed, which was previously described as the Hawthorne effect.⁴³ In another word, this might be explained by the effect of social support, such as positive communication and adherence messages transferred by the study staff while being followed up. A positive provider-patient relationship was reported to be associated with higher ART adherence among people living with HIV.⁴⁴ This finding was in line with several previous reports, in which self-compassion was found to be associated with enhanced participation in behaviors that promote health⁴⁵ and improved physical well-being and health.^46,47

These contextual supports may have functioned as co-interventions, attenuating detectable differences between trial arms. Such “trial effects” or Hawthorne-like effects are well recognized in pragmatic HIV research, particularly in settings where participants receive intensified follow-up, counseling, and monitoring beyond routine care. Other unmeasured factors, including regression to the mean, concurrent programmatic improvements at study sites, or changes in clinical management during the study period, may also have contributed to parallel viral load trajectories across groups. Collectively, these findings underscore the importance of accounting for background psychosocial and service-level supports when interpreting intervention effects on virologic outcomes in real-world HIV care settings.

From an implementation perspective, heterogeneity in attendance and engagement across group sessions may have diluted intervention effects. Group-based formats, while scalable and resource-efficient, may be less responsive to individualized barriers commonly faced by people living with HIV, including stigma, disclosure concerns, family or partner dynamics, and co-occurring psychosocial stressors. These factors are particularly salient in real-world service delivery and may require more tailored or multicomponent approaches to achieve measurable changes in mental health and adherence outcomes.

Notably, the absence of statistically significant effects does not preclude clinical relevance. Small, nonsignificant improvements may still be meaningful in practice when interventions are low-cost, acceptable to clients, and feasible for integration into existing service platforms. Moreover, group counseling may confer additional benefits not fully captured by the primary outcome measures, such as enhanced peer support, normalization of experiences, and improved therapeutic alliance with care providers, which may contribute to longer-term engagement in HIV care.

Although the group counseling program integrated with music did not demonstrate statistically significant superiority over the comparison arms in improving ART adherence, the direction and pattern of change observed suggest a potentially meaningful relationship between participation in the music-enhanced group sessions and engagement in medication-taking behaviors. Music-based activities could facilitate emotional expression, stress reduction, and group cohesion, thereby creating a supportive and nonstigmatizing environment that fostered motivation and self-efficacy related to treatment adherence. Such mechanisms are consistent with psychosocial theories emphasizing the role of affecting regulation, social connectedness, and intrinsic motivation in sustaining long-term adherence behaviors among people living with HIV.

Importantly, the lack of statistically significant between-group differences should not be interpreted as evidence of no clinical or practical relevance. In the context of a small sample size and high levels of background psychosocial support provided to all participants, the incremental benefit of the music-integrated group counseling may have been diluted. Nevertheless, the observed trends underscore the potential value of incorporating creative, culturally sensitive modalities such as music into adherence support interventions, particularly for adolescents and young adults, for whom conventional counseling approaches may be perceived as less engaging. Future trials with an increased number of participants are warranted to further delineate the specific contribution of music-based group counseling to ART adherence and to identify subgroups most likely to benefit.

The study strength included study participants with similar characteristics in each intervention group using random assignment, hair DTG level measurement, and available HIV-RNA results, which reflected HIV treatment adherence and outcomes. Future research should explore larger samples, that is, more than one group assigned to the same intervention, extended follow-up periods, and integration with digital health interventions to sustain ART adherence improvements and strengthen self-esteem.

Conclusion

The study demonstrated that integrated music and counseling programs for AYA-PHIV were associated with improving self-esteem and reducing internalized HIV-related stigma. Study participants reported better antiretroviral treatment adherence and improved virologic outcomes.

We support including psychosocial interventions in comprehensive HIV care strategies, especially for high-risk individuals with mental health or adherence issues. A larger study with an increased sample size is warranted to confirm these findings, improve statistical power, and better delineate the clinical and programmatic implications of group-based counseling interventions within routine HIV care settings. Integrated music and counseling programs could be adopted as they offer promising avenues for enhancing self-esteem and reducing stigma among AYA-PHIV.

Supplemental Material

sj-pdf-1-jia-10.1177_23259582261439584 - Supplemental material for Effects of Music and Counseling Programs on Mental Health and Antiretroviral Treatment Adherence of Adolescents and Young Adults Living With HIV in Northern Thailand, a Randomized Trial

Supplemental material, sj-pdf-1-jia-10.1177_23259582261439584 for Effects of Music and Counseling Programs on Mental Health and Antiretroviral Treatment Adherence of Adolescents and Young Adults Living With HIV in Northern Thailand, a Randomized Trial by Chintana Khamrong, Rattika Thammalangka, Kalunyu Kotchawat, Nakhon Khamrong, Linda Aurpibul, Jakkrit Klaphajone, Teerawan Teerapong and Natthapol Kosashunhanan in Journal of the International Association of Providers of AIDS Care (JIAPAC)

Supplemental Material

sj-docx-3-jia-10.1177_23259582261439584 - Supplemental material for Effects of Music and Counseling Programs on Mental Health and Antiretroviral Treatment Adherence of Adolescents and Young Adults Living With HIV in Northern Thailand, a Randomized Trial

Supplemental material, sj-docx-3-jia-10.1177_23259582261439584 for Effects of Music and Counseling Programs on Mental Health and Antiretroviral Treatment Adherence of Adolescents and Young Adults Living With HIV in Northern Thailand, a Randomized Trial by Chintana Khamrong, Rattika Thammalangka, Kalunyu Kotchawat, Nakhon Khamrong, Linda Aurpibul, Jakkrit Klaphajone, Teerawan Teerapong and Natthapol Kosashunhanan in Journal of the International Association of Providers of AIDS Care (JIAPAC)

Supplemental Material

sj-docx-4-jia-10.1177_23259582261439584 - Supplemental material for Effects of Music and Counseling Programs on Mental Health and Antiretroviral Treatment Adherence of Adolescents and Young Adults Living With HIV in Northern Thailand, a Randomized Trial

Supplemental material, sj-docx-4-jia-10.1177_23259582261439584 for Effects of Music and Counseling Programs on Mental Health and Antiretroviral Treatment Adherence of Adolescents and Young Adults Living With HIV in Northern Thailand, a Randomized Trial by Chintana Khamrong, Rattika Thammalangka, Kalunyu Kotchawat, Nakhon Khamrong, Linda Aurpibul, Jakkrit Klaphajone, Teerawan Teerapong and Natthapol Kosashunhanan in Journal of the International Association of Providers of AIDS Care (JIAPAC)

Supplemental Material

sj-docx-5-jia-10.1177_23259582261439584 - Supplemental material for Effects of Music and Counseling Programs on Mental Health and Antiretroviral Treatment Adherence of Adolescents and Young Adults Living With HIV in Northern Thailand, a Randomized Trial

Supplemental material, sj-docx-5-jia-10.1177_23259582261439584 for Effects of Music and Counseling Programs on Mental Health and Antiretroviral Treatment Adherence of Adolescents and Young Adults Living With HIV in Northern Thailand, a Randomized Trial by Chintana Khamrong, Rattika Thammalangka, Kalunyu Kotchawat, Nakhon Khamrong, Linda Aurpibul, Jakkrit Klaphajone, Teerawan Teerapong and Natthapol Kosashunhanan in Journal of the International Association of Providers of AIDS Care (JIAPAC)

Footnotes

Acknowledgments

This study would not be possible without the dedication and diligent effort of 24 adolescents and young adults, their parents/guardians, and healthcare providers, Family clinic (IMPAACT) team, RIHES, CMU, my main advisor, Dr Natthapol Kosashunhanan, MD, my co-advisors, Dr Linda Aurpibul, MD, Assoc Prof Jakkrit Klaphajone, MD, and Assoc Prof Dr Teerawan Teerapong.

ORCID iDs

Chintana Khamrong

Rattika Thammalangka

Kalunyu Kotchawat

Nakhon Khamrong

Linda Aurpibul

Jakkrit Klaphajone

Teerawan Teerapong

Natthapol Kosashunhanan

Ethical Considerations

Ethical approval for this study was obtained from the Human Experimentation Committee (HEC), Research Institute for Health Sciences, Chiang Mai University (certificate approval number 30/2023). Written informed consent was obtained from all participants, consent and assent were also obtained from the caregiver and the participant, respectively, prior to enrollment. This randomized controlled trial was registered at the Thai Clinical Trials Registry. (Identifier: TCTR20241231002). The study was conducted in accordance with the ethical principles of the Declaration of Helsinki and reported following the CONSORT guidelines.

Consent to Participate

Before enrollment, participants provided written informed consent for study participation, inclusion, and use of data in publication; for those under 20 years of age, consent and assent were also obtained from the caregiver and the participant, respectively.

Funding

The authors disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: Chintana Khamrong received funding from the Research Institute for Health Sciences, Chiang Mai University, Chiang Mai, Thailand, with grant number (grant number 013/2566) that contributed to this study.

Declaration of Conflicting Interests

The authors declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Supplemental Material

Supplemental material for this article is available online.

References

UNICEF. Global and regional trends. Accessed August 18, 2023. https://data.unicef.org/topic/hivaids/global-regional-trends/.

State of HIV Epidemic in Thailand in 2022. Accessed August 18, 2023. https://hivhub.ddc.moph.go.th/Download/RRTTR/HIV_Factsheet2022_EN_Final_220966.pdf.

Piaget

. The psychology of the child. Basic Books; 1972.

Harter

. The construction of the self: Developmental and sociocultural foundations. Guilford Press; 2012.

Hankin

Abramson

Moffitt

, et al. Development of depression from preadolescence to young adulthood: Emerging gender differences in a 10-year longitudinal study. J Abnorm Psychol. 2010;119(5):811‐833.

Mallinson

Rajabiun

Coleman

. The provider role in client engagement in HIV care. AIDS Patient Care STDS. 2007;21(s1):S77‐S84.

Dowshen

D'Angelo

Cheng

, et al. Adolescents living with HIV: Challenges in transition to adult care. Mt Sinai J Med. 2010;77(2):161‐166.

Kang

Delzell

Chhabra

Oberdorfer

. Factors associated with high rates of antiretroviral medication adherence among youth living with perinatal HIV in Thailand. Int J STD AIDS. 2015;26(8):534‐541.

Boivin

Barlow-Mosha

Chernoff

, et al. Neuropsychological performance in African children with HIV enrolled in a multisite antiretroviral clinical trial. AIDS. 2018;32(2):189‐204.

10.

Van der Straten

Shapley-Quinn

Reddy

, et al. It's all in the timing: Acceptability of a smartphone-based medication adherence intervention among adults with HIV. J Assoc Nurses AIDS Care. 2019;30(4):420‐432.

11.

Machado

Galano

Padoin

. Shared decision-making in the treatment of adolescents living with HIV: A systematic review. J Pediatr Nurs. 2019;47:e3‐e12.

12.

Wilson

Lee

Michaud

, et al. Validation of a new three-item self-report measure for medication adherence. AIDS Behav. 2006;10(5):541‐548.

13.

Amico

Toro-Alfonso

Fisher

. An empirical test of the information-motivation-behavioral skills model of antiretroviral therapy adherence. AIDS Care. 2005;17(6):661‐673.

14.

Lolekha

Boon-Yasidhi

Leowsrisook

, et al. Mental health among perinatally HIV-infected adolescents in Thailand: A case-control study. J Int AIDS Soc. 2014;17(1):18914.

15.

Chantaratin

Trimetha

Werarak

, et al. Depression and anxiety in youth and young adults living with HIV: Frequency and associated factors in Thai setting. J Int Assoc Provid AIDS Care. 2022;21(1).

16.

Trotzers

. The counselor and the group: Integrating theory, training and practice. Brooks/Core Publishing Co; 1977.

17.

Niyonzima

Halim

Asante

. The effectiveness of music therapy on antiretroviral therapy (ART) adherence and depression among people living with HIV/AIDS: A systematic review and meta-analysis. J Assoc Nurses AIDS Care. 2019;30(5):e92‐e108.

18.

Trotzers

. The counselor and the group: Integrating theory, training and practice. Brooks/Core Publishing Co; 1977.

19.

Starcke

Mayr

von Georgi

. Emotion modulation through music after sadness induction – The ISO principle in a controlled experimental study. Int J Environ Res Public Health. 2021 Nov 26;18(23):12486.

20.

Lovibond

. Manual for the Depression Anxiety Stress Scales. 2nd. Ed. Psychology Foundation; 1995.

21.

Wongpakaran

. A comparison of reliability and construct validity between the original and revised versions of the Rosenberg Self-Esteem Scale. Psychiatry Investig. 2012;Mar;9(1):54‐58.

22.

Uthis

Thiamnak

Srimanee

, et al. The Thai internalized HIV-related Stigma Scale. Front Psychol. 2023;14:4471.

23.

Schulz

Altman

Moher

, CONSORT Group. CONSORT 2010 Statement: Updated guidelines for reporting parallel group randomized trials. Br Med J. 2010;340:c332.

24.

The jamovi project. jamovi. (Version 2.7) [Computer Software]. Retrieved from https://www.jamovi.org. 2025.

25.

R Core Team. R: A Language and environment for statistical computing. (Version 4.5) [Computer software]. Retrieved from https://cran.r-project.org. (R packages retrieved from CRAN snapshot 2025-05-25). 2025.

26.

Kim

. ppcor: Partial and Semi-Partial (Part) Correlation. [R package]. Retrieved from https://cran.r-project.org/package = ppcor. 2015.

27.

Gallucci

. GAMLj: General analyses for linear models (Version 3.6.5). [jamovi module]. Retrieved from https://gamlj.github.io/. 2019.

28.

Gallucci

. Model goodness of fit in GAMLj. link. 2020.

29.

Lüdecke

Ben-Shachar

, Patil & Makowski. Extracting, Computing and Exploring the Parameters of Statistical Models using R. CRAN. link. 2020

30.

Brown

Dutton

Cook

. From the top down: Self-esteem and self-evaluation. Cogn Emot. 2001;15(5):615-631.

31.

Heal

Wigram

. Music Therapy in Health and Education. Edited by Margaret. Jessica Kingsley. 1993. p. 292.

32.

Henderson

. Effects of a music therapy program upon awareness of mood in music, group cohesion, and self-esteem among hospitalized adolescent patients. J Music Ther. 1983;20(1):14-20.

33.

Gardstrom

. Positive peer culture: A working definition for the music therapist. Music Ther Perspect. 1987;4(1):19-23.

34.

Smeijsters

Hurk

. Music therapy helping to work through grief and finding a personal identity. J Music Ther. 1999;36(3):222‐252.

35.

Freed

. Song writing with the chemically dependent. Music Ther Perspect. 1987;4(1):13-18.

36.

Willis

Frewin

Miller

, et al. “My Story” – HIV positive adolescents tell their story through film. Child Youth Serv Rev. 2014;45(C):129-136.

37.

Willis

Mavhu

Wogrin

, et al. Understanding the experience and manifestation of depression in adolescents living with HIV in Harare, Zimbabwe. PLoS One. 2018; 13(1):e0190423.

38.

McAdam

Kambugu

. Thayer

McAdam

Twijukye, C

. Sidebar: The crerativity initiative: empowering people living with HIV to become change agents In: S.J.

Marlink RGaT

, editor. From the Ground Up: Building Comprehensive HIV/AIDS Care Programs in Resource-Limited Settings. Elizabeth Glazer Pediatric AIDS Foundation 2009; 264‐276.

39.

Ferris France

Lyons

Cioringa

, et al. Interventions to reduce self-stigma among people living with HIV: A systematic review. Stigma and Health. Advance online publication. 2025:5-13.

40.

Werner-Seidler

Perry

Calear

, et al. School-based depression and anxiety prevention programs for young people: A systematic review and meta-analysis. Clin Psychol Rev. 2017;51(1):30-47.

41.

Bhana

Mellins

Petersen

, et al. The VUKA family program: Piloting a family-based psychosocial intervention to promote health and mental health among HIV infected early adolescents in South Africa. AIDS Care. 2014;26(1):1‐11.

42.

Kaihin

Kasatpibal

Chitreechuer

, et al.

Effect of an empowerment intervention on antiretroviral drug adherence in Thai youth.

Behav Med (Washington, DC). 2015;41(4):186‐194.

43.

Van Onzenoort

Menger

Neef

, et al. Participation in a clinical trial enhances adherence and persistence to treatment: A retrospective cohort study. Hypertension. 2011;58(4):573‐578.

44.

Schneider

Kaplan

Greenfield

, et al. Better physician-patient relationships are associated with higher reported adherence to antiretroviral therapy in patients with HIV infection. J Gen Intern Med. 2004;19(11):1096‐1103.

45.

Dunne

Sheffield

Chilcot

. Brief report: Self-compassion, physical health and the mediating role of health-promoting behaviours. J Health Psychol. 2018;23(7):993‐999.

46.

Sirois

Molnar

Hirsch

JK.

Self-Compassion, stress, and coping in the context of chronic illness. Self Identity. 2015; 14(3):334‐347.

47.

Brown

Bryant

Brown

, et al. Self-compassion, attitudes to ageing and indicators of health and well-being among midlife women. Aging Ment Health. 2016;20(10):1035‐1043.