ART use was dichotomized as use of any ART since the last study visit versus no therapy used

The classification of child-appealing was informed by laws in other states and constructed with authors’ own understanding, which may not reflect California lawmakers’ intention or completely align with recently released new regulations. Further, there might be inevitable measurement errors even after two field workers discussed and resolved discrepancies between them. Lastly, this study only gathered data on RMDs in closest proximity to public schools. Results may not be generalizable to RMDs around private schools or children’s homes. Marijuana use is common among persons living with HIV as studies have reported prevalence rates of current marijuana use between 24 and 56 % as compared to approximately 7 % in the general United States population. Men who have sex with men report higher rates of current and past-year marijuana use than their heterosexual counterparts. Several studies report that persons living with HIV use marijuana to alleviate stress, anxiety, depression, HIV-related symptoms and side effects of antiretroviral therapy. In one recent study, among HIV-seropositive persons who inject drugs and who recently seroconverted, heavy cannabis use was associated with lower plasma viral load levels. The therapeutic effects of marijuana are proposed to be mediated via the actions of active cannabinoid chemicals in cannabis grow equipment—cannabidiol—at specific receptor sites: cannabinoid receptors located mainly on cells and tissues of the immune system. In contrast the primary psychoactive cannabinoid in marijuana: tetrahydrocannabidiol binds to and activates another receptor site: cannabinoid receptor located mainly in areas of the brain to produce the euphoric and cognitive impairing effects of marijuana.

Accordingly, there are concerns that marijuana use may be associated with poorer HIV treatment outcomes. Previous studies have found marijuana use to be associated with decreased cognitive function as well as reduced ART adherence, which is crucial for persons living with HIV as optimal adherence to ART medications is required for long-term viral suppression. Effective prevention strategies to reduce unhealthy or harmful marijuana use require an indepth understanding of subgroups with different patterns of use. Despite the published evidence that marijuana use is common among HIV+ individuals and MSM and the potential adverse health outcomes associated with its use in these populations, very little is known about the patterns of marijuana use or how patterns of marijuana use may change over time in these populations. Developmental research suggests different rather than similar pathways via which individuals initiate and progress to unhealthy or problem substance use over the life course. For instance, individuals who start using substances at an early age have increased risk of progressing to problem use and developing use disorders. Among HIV+ women, depressive symptoms and the presence of hepatitis C infection was associated with a pattern of persistent heavy drinking over time. Another study found that low income and concurrent substance use were factors that predicted consistent hazardous drinking among HIV + MSM. Therefore, understanding the natural history of marijuana use and the identification of different trajectories of use over time is important in order for intervention programs to be most effective. For instance, the identification of different patterns of marijuana use over time can help characterize subgroups of individuals with the greatest risk of progressing to heavy patterns of marijuana use and reveal unique predictors of such patterns of use which can be used to inform targeted intervention programs.

To the best of our knowledge, there has not been a published report that has followed HIV+ individuals and MSM longitudinally over an extended period to characterize the natural history of their marijuana use. Past studies on substance use patterns in these populations have often focused on alcohol or heavy episodic drinking, cigarette smoking or stimulant use. Therefore, the objectives of the current study are to characterize the longitudinal trajectories of marijuana use in a sample of HIV-seropositive and HIVseronegative MSM over a period of 29 years, and to identify factors associated with unique trajectories of marijuana use, as well as those that can change over time that may modify the course of the trajectory.The MACS study protocols were approved by the institutional review boards at the respective recruitment centers and their community affiliates and informed consent was obtained from all participants. MACS participants return every 6 months for a physical examination, collection of blood specimens and completion of a detailed interview and questionnaires. The interview and questionnaires include demographic, psychosocial, behavioral and medical history data. The questions about their recreational drug use, including marijuana, alcohol, poppers, cocaine, crack, heroin, methamphetamine, ecstasy, injection drug use as well as smoking history since their last visit were collected using Audio Computer Assisted Self-Interviewing , an approach previously demonstrated to provide more accurate assessments of ‘sensitive behaviors’ than interview-administered questionnaires among MSM. This analysis included data collected from standardized marijuana use questions from semiannual study visits 1 through visit 59 . The study sample included 3658 participants who had data about marijuana use for at least 25 % of their possible study visits during the follow-up period.

Specifically, the men enrolled in 1984–1985 and 1987–1991 had 15 and 13 visits or more respectively, whereas, the men enrolled in 2001–2003 had 6 or more visits. The median years of follow-up was 11.5 years . Participants were classified as never, former and current smokers of cigarettes at each study visit. Participants were asked two questions including: whether they ever smoked cigarettes and whether they smoke cigarettes now. Participants were considered to be current smokers if they responded ‘Yes’ to both questions. Participants were categorized as former smokers if they answered ‘Yes’ to ever smoking cigarettes and ‘No’ to the smoking cigarettes now. Never smokers included participants who answered ‘No’ to both questions. In addition, among current smokers, pack-years of smoking at initial visit and at each subsequent visit was calculated using participants’ responses to questions about the number of packs of cigarettes smoked per day. HIV serostatus was assessed using an enzyme-linked immunosorbent assay with confirmatory Western blot tests on all MACS participants at each participant’s initial visit and at each study visit for participants who were initially HIV−. Standardized flow cytometry was used to quantify CD4 + T-lymphocyte subset levels by each MACS site and categorized as ≤200/mm3 , 201–500/mm3 , and >500/mm3 . Levels of plasma HIV RNA were measured using either the standard reverse transcription-polymerase chain reaction assay or with the Roche ultrasensitive assay were used to create a dichotomous variable to denote detectable versus undetectable viral load. Hepatitis C virus infection status was categorized as HCV negative if HCV antibody testing was negative. Participants were classified at each MACS study visit as HCV positive if they were found to be in the process of seroconversion, acute infection, chronic infection, clearing , or previously HCV positive, but now clear of HCV RNA. In addition to the aforementioned covariates,horticulture rack we considered that the trajectories of marijuana use over time among HIV+ participants may be influenced by factors specific to HIV-infection. We used participant’s self-reported frequency of marijuana use across the follow-up period to identify trajectories using a semi-parametric group-based mixture model: PROC TRAJ SAS procedure. This approach sorts each participant’s frequency of marijuana use over their follow-up period into ‘clusters’ and estimates a single model—consisting of distinct trajectories. The procedure calculates the probability of each participant belonging to each trajectory group and assigns individuals into trajectories based on their highest probability of trajectory membership. Participants were followed from the time of enrollment until either the time of death, lost to follow up or until the end of the study period . We began by fitting a series of models with two to five trajectories by assuming linear, quadratic and cubic shape of the trajectory group curves. Several factors were considered in determining model fit and the optimal number of trajectory groups that best represented the heterogeneity of groups within the data: including, a priori knowledge from previous research on trajectories of marijuana use, model fit statistics including Bayesian information criterion, Akaike Information Criterion, average posterior probability of group membership, significance of the shape of the trajectory group curves , and size of the group membership .

Model fitting was an iterative process, starting with a quadratic specification for the shape of the trajectory group curves and assessing whether an additional group resulted in a better model fit based on the aforementioned criteria. We then estimated higher order shapes of the trajectory group curves and subsequently dropped non-significant terms. Models used a zeroinflated Poisson distribution to account for the large number of participants who reported not using marijuana. After the optimal number of trajectory groups and shape of trajectory change were selected, we included covariates of interest to the trajectory models. For this analysis, two types of covariates were considered: time-fixed/risk factors of trajectory group membership and time varying covariates. These time-fixed/risk factors comprise characteristics established before or at the time of the initial period of trajectories that may serve to predict membership in a given trajectory. Time-varying covariates measured during the course of the trajectory provide trajectory group-specific estimates of whether these covariates alter the course of the trajectory. One advantage of the PROC TRAJ software is that it allows for joint estimation of the parameters that describe the shape of the trajectory group curves, adjusted odds ratio and the coefficient estimates . We estimated models for all participants as well as by HIV serostatus. The analysis of all participants was adjusted for sociodemographic characteristics, depressive symptoms, substance use variables, hepatitis C infection status, attrition variables, and HIV serostatus. To account for potential differences in marijuana use by geographic region/site and MACS enrollment cohort, all models included variables for MACS center and enrollment cohort. In the analysis restricted to HIV+ participants, we included other clinical factors relevant to HIV+ status such as ART use, CD4 counts, and viral load detectability. All analysis was performed in SAS 9.4 . Using data for the entire sample, participants’ self-reported frequency of marijuana use across the follow-up period identified four groups with distinct trajectories of marijuana use. We chose a four-group solution based on model parsimony, interpretability of trajectories, BIC and AIC values, significance of the polynomial growth terms, average posterior probabilities and trajectory group size membership . Model fit information and average posterior probabilities of all models are displayed in supplemental Tables 2 to 5. Figure 1 displays the trajectories of marijuana use of these four groups, which we labelled as: “Abstainer/Infrequent”, “Decreasers”, “Increasers” and “Chronic high” trajectory groups. The abstainer or infrequent use group was characterized by a group of men who abstained from or infrequently used marijuana during the follow-up period. The decreaser group consisted of a group of men who reduced their marijuana use from nearly weekly use to infrequent use over the follow-up period. The increaser group comprised a group of men who initially decreased their marijuana use during the first 10 years of follow-up, after which they began to increase their use over time. The chronic high group represents a group of men who persistently used marijuana nearly daily over the follow-up period. Figure 2 displays trajectories of marijuana use among HIV+ participants: 61 % were in the abstainer/ infrequent use group, 14 % were in the decreaser group, 14 % in the increaser group, and 11 % in the chronic high group. Table 2 displays the baseline characteristics of the entire sample by the four identified trajectory groups. The median number of visits was lower among those in the increaser trajectory group. Participants in the abstainer/infrequent use group were older at baseline compared to the other groups. Frequency of marijuana use at baseline varied across the marijuana trajectory groups: as the proportion of daily users were < 1 % in the abstainer/ infrequent, 3 % in the decreasers, 10 % in the increasers, and 54 % in the chronic high groups. Racial status , detectable HIV viral load and CD4 counts were similar across the marijuana trajectory groups. This study utilized data from the MACS cohort to assess different patterns of marijuana use and to examine both risk factors and time-varying correlates associated with the different trajectories of marijuana use. Our analysis revealed that MSM in the MACS exhibited four distinct trajectories of marijuana use over time, including: abstainer/infrequent, decreasers, increasers and chronic high groups.