Influence of Internalised Homonegativity on Sexual Risk Behaviour of Men Who Have Sex with Men in Spain

In a sample of men who have sex with men (MSM) (N = 3436) in Spain who bear intrinsic HIV risk, we investigated how internalised homonegativity (IH) is associated with the number of non-steady male partners with condomless intercourse (as a proxy of sexual risk behaviour). Using structural equation modelling (SEM), we examined the relationship between IH and sexual risk behaviour, and mediating effects of HIV/PrEP knowledge and substance use during sex on this relationship. We found no direct association between IH and sexual risk behaviour, nor did IH influence substance use during sex. In line with our hypothesis, association between IH and sexual risk behaviour was significant when mediated by HIV/PrEP knowledge. We found that as IH increased, sexual risk behaviour decreased, because higher IH was associated with lower HIV/PrEP knowledge while higher HIV/PrEP knowledge was associated with increased non-condom use with non-steady partners. Substance use during sex was significantly associated with sexual risk behaviour. Our results emphasize the continuing importance of prevention strategies focused on behavioural changes and community level interventions, especially targeting substance use.


Introduction
Forty years into the epidemic, sex between men remains the predominant mode of HIV transmission in the countries of the European Economic Area, accounting for 38.7% of all new HIV diagnoses in 2019 (European Centre for Disease Prevention and Control & World Health Organization, 2020). In Spain, 3381 new HIV 1 3 Influence of Internalised Homonegativity on Sexual Risk… diagnoses were reported in 2017, with more than half (54.3%) of these in men who have sex with men (MSM) (Centro Nacional de Epidemiología, 2018). A recent biobehavioural study conducted in 13 European cities showed that Barcelona (Spain) has one of the highest rates of HIV among MSM (Mirandola et al., 2018). In addition to the spread of HIV, increasing numbers of sexually transmitted infections (STIs) have been reported in Spain after 2005(Centro Nacional de Epidemiología, 2020. Sexual transmission risk among MSM is mediated by multiple social and structural factors that influence individuals' sexual practices (Baral et al., 2013). These factors not only influence risk behaviours, but also directly limit options for accessing prevention services for HIV and other STIs and jeopardise prevention efforts (Andrinopoulos et al., 2015;Velter et al., 2015). Meyer (2003) developed the minority stress model, which refers to the "excess stress to which individuals from stigmatized social categories are exposed as a result of their social, often a minority, position", and which may partly explain behaviours that increase HIV transmission risk, such as substance use and condomless sex (Meyer, 2003). Internalised Homonegativity (IH), defined as negative feelings about one's homosexuality (Herek, 2004), is the product of social and political stigma and bias instead of a response which stems from within individuals. IH is one of the minority stressors that has been expanded upon in Meyer's (2003) minority stress model, and has a documented damaging effect on the mental health and well-being of sexual minorities . However, inconsistent research results point to uncertainties about the extent to which IH influences engagement in sexual risk behaviours Puckett et al., 2017). In Catalonia, an autonomous community of Spain, previous research among gay, bisexual and other MSM found that IH was an independent predictor of sexual risk behaviours (SRB). This has also been shown for European (Berg et al., 2015) and non-European countries . However, other studies have not found significant associations between IH and SRB (Dudley et al., 2004;Kashubeck-West & Szymanski, 2008a). These inconsistencies on the literature may reflect, in part, the existence of potential mediating variables that affect the relation between IH and sexual risk behaviours (Kashubeck-West & Szymanski, 2008a;Newcomb & Mustanski, 2010b).
Similar inconsistent research results exist regarding the relation between IH and drug use. Some researchers found that IH was associated with higher levels of drug use, whereas others found either no significant relation with drug use or negative associations (Puckett et al., 2017;Ross et al., 2001). The effect of sexualized drug use on SRB has also been studied extensively. Previous studies in Spain found evidence that a higher prevalence of drug use consumption was associated with a higher prevalence of condomless anal sex or sex with multiple partners (Fernández-Dávila & Zaragoza Lorca, 2009;Folch et al., 2006Folch et al., , 2010González-Baeza et al., 2018).
Men with high level of IH are less likely to be involved in the gay community and likewise more isolated from getting informed about HIV prevention and risk reduction programs. Among 569 gay and bisexual men Huebner et al. (Huebner et al., 2002) found that IH was negatively associated with the number of HIV-related services that these men were aware of. Further, of the 443 MSM who had heard about at least one service, IH was not a significant predictor for participation in these services, when controlled for education levels. Therefore, high IH was a barrier for men to be aware of these services in the first place. On the other hand, the results of the 2010 wave of EMIS among more than 144,000 MSM across 38 countries in Europe provide evidence that IH was positively associated with less knowledge about HIV and HIV testing (Berg et al., 2013). Likewise, in a sample of substance using HIV-negative and unknown status gay and bisexual men in New York City, it was found that community connection was protective against sexual risk and drug use, especially among younger men (Lelutiu-Weinberger et al., 2013). Previous research also documented that IH can hinder gay men's connection to the gay community (Goldbach et al., 2015;Moody et al., 2018), which can, in part, explain a possible link between IH and SRB, through a lack of gay community acculturation and where the targeted information is available for gay and bisexual men (Williamson, 2000).
An improved understanding of the impact of critical factors that mediate the relation between IH and SRB would be important in general, and in Spain specifically, in order to tailor community services to those MSM with higher levels of IH and in turn at higher risk of HIV infection. Therefore, using structural equation modelling (SEM), our study aimed to disentangle the possible influence of drug use and knowledge regarding HIV and PrEP on the relation between IH and SRB in a national sample of MSM living in Spain and recruited online. We had three hypotheses. First, we tested the 'IH will be positively associated with sexual risk behaviour (SRB)' hypothesis. As abovementioned, we argue that inconsistencies in previous research regarding the relationship between IH and SRB may reflect particular roles of mediator variables, such as frequency of sex under the influence of substances and HIV/PrEP knowledge. Thus, we tested our second and third hypotheses; HIV/ PrEP knowledge will strongly mediate the relationship between IH and SRB' and 'sex under the influence of substances will strongly mediate the relationship between IH and SRB.'

Sample and Data
We used data from the 2017 wave of the European MSM Internet Survey (EMIS-2017). The detailed methods have been reported elsewhere (Weatherburn et al., 2020). In summary, EMIS-2017 was a 33-language, internet-based, self-completion survey for men living in Europe who have sex with men and/or are sexually attracted to other men. No financial incentives were given to participants; no personal identifying information (including IP addresses) were collected. More background information is available at www. emis2 017. eu. The sub-sample of MSM living in Spain consisted of 10,652 respondents, including men living in the autonomous provinces of Canarias, Ceuta, Melilla (geographically outside Europe) as well as men living in the Principality of Andorra (but not in the British Overseas Territory of Gibraltar), with 92.1% using the Spanish (Castilian) version of the survey (no other co-official languages of Spain, such as Catalan/Valencian, Galician, or Basque were offered).
The Short Internalised Homonegativity Scale (SIHS) was randomly distributed to half of respondents (N = 5310) of which 4632 answered all SIHS items. This random distribution has been done to avoid losing participants because of asking too many questions. MSM who did not provide answers to all seven items were excluded. We also excluded 583 HIV-diagnosed MSM who reported having undetectable viral load, and 78 PrEP users, because condomless anal intercourse among men with undetectable viral load or using PrEP does not bear any intrinsic HIV risk. The analytic sample thus consists of 3436 MSM living in Spain.

Measures
Internalised homonegativity-To assess IH, we used the 7-item SIHS (Berg et al., 2013;Tran et al., 2018). This term was defined by Ross and colleagues (Berg et al., 2013), and the construct of IH commonly refers to internalization of homophobic attitudes within lesbian, gay, and bisexual individuals (Newcomb & Mustanski, 2011), thus an attachment of external homonegativity (often incorrectly referred to as 'homophobia') to the sense of self (Malyon, 1982;Stein & Cohen, 1986). EMIS-2017 participants answered 7 items on a 7-point disagree-agree (with does-notapply) scale. These items were "social situations with gay men make me feel uncomfortable"; "homosexuality is morally acceptable to me"; "even if I could change my sexual orientation, I wouldn't"; "I feel comfortable in gay bars"; "I feel comfortable being seen in public with an obviously gay person"; "I feel comfortable being a homosexual man"; "even if I could change my sexual orientation, I wouldn't". The validity and reliability of SIHS were also confirmed across 38 European countries, with multigroup validation for 7-item scale fit indices showing good fit to data from 38 country groups (CFI = 0.982, TLI = 0.983, and RMSEA = 0.032) (see, Tran et al. (2018) for further statistics).
Sexual Risk Behaviour (SRB)-SRB of the respondents was assessed with a single question: "how many non-steady male partners have you had intercourse without a condom with in the last 12 months?" Here, participants were informed that non-steady partners mean "men you have had sex with once only, and men you have sex with more than once but who you don't think of as a steady partner (including one-night stands, anonymous and casual partners, regular sex buddies)". The possible answer options for this question in the survey ranged from 0 to 15; with numbers 0 to 10 equivalent to their values, and numbers 11 to 15 indicating 11-20, 21-30, 31-40, 41-50, and more than 50 partners respondents had condomless intercourse with. We recoded this variable into seven categories: 0; 1; 2-5; 6-10; 11-20; 20-50; and more than 50. We would like to highlight that our definition of the risk behaviour is related to HIV risk and does not relate directly to other STIs.
HIV/PrEP Knowledge -To construct the HIV/PrEP knowledge latent variable, two measures were used: HIV knowledge and PrEP knowledge. HIV knowledge was constructed from seven items, assessed with a 5-point knowledge response set, with possible answers including "I do not believe this", "I wasn't sure about this", and "I knew this already". These items were "AIDS is caused by a virus called HIV"; "if someone becomes infected with HIV it may take several weeks before it can be detected in a test"; "you cannot be confident about whether someone has HIV or not from their appearance"; "there is a medical test that can show whether or not you have HIV"; "There is currently no cure for HIV infection"; "HIV infection can be controlled with medicines so that its impact on health is much less"; "a person with HIV who is on effective treatment (called 'undetectable viral load') cannot pass their virus to someone else during sex." PrEP knowledge included three items assessed with the same response set: "Pre-Exposure Prophylaxis (PrEP) involves someone who does not have HIV taking pills before as well as after sex to prevent them getting HIV"; "PrEP can be taken as a single daily pill if someone does not know in advance when they will have sex"; "If someone knows in advance when they will have sex, PrEP needs to be taken as a double dose approximately 24 h before sex and then at both 24 and 48 h after the double dose." Each of these 10 items were recoded into a dummy variable, with value 1 indicating "I knew this already," and value 0 indicating all the other answers. Then, we created an additive scale with these 10 items, ranging from 0 to 10. With each factual knowledge (I knew this already) of each question, respondents scored one point in the additive scale.
Substance Use-We used six observed variables for the substance use latent variable, based on how long ago respondents used substances in any context. The six substances (see Table 1) were assessed with an 8-point recency scale, ranging from (1) "never" to (8) "in the past 24 h" (after inverting the original scale).
Sex Under the Influence of Substances (SUIS) -For this variable, the respondents were asked, "in the last 12 months, how much of the sex you've had with men has been under the influence of alcohol or any other drug?" The possible answers for this question ranged from (1) "none of it" to (7) "all of it".

Statistical Analysis
We use RStudio and the 'lavaan' package (Rosseel, 2012) to analyse the hypothesized structural equation model. Prior to the analysis, the data was checked for multicollinearity, missing data, departures from normality and distributions. Multicollinearity was not present. Missing data were handled with pairwise deletion.
First, we examined descriptive statistics and correlation among variables used. Second, we estimated the confirmatory factor analysis (CFA) model. Finally, we estimated the hypothesized SEM presented in Fig. 1. Since we will be testing mediation effects, we follow the study of Shrout and Bolger, (2002) and use bias-corrected bootstrap method for estimating our model. This estimation method allows interval estimated without relying on a distribution assumption. Bias-corrected bootstrap estimation adjusts for possible bias and problematic skewness, if any, in the bootstrap samples' distribution (Beaujean, 2014). Therefore, we estimated our SEM using a bootstrapped MLM estimator. The SEM had 16 observed variables, and 136 known parameters ((16*17) * ½ = 136). The total number of unknowns were 36 Influence of Internalised Homonegativity on Sexual Risk… and constitutive of; 4 covariances, 16 factor loadings, 16 error variances. Thus, our model was over-identified with 100 (136-36) degrees of freedom, and over-identification is a necessary factor in structural model's fit to data (Weston & Gore, 2006).
Proceeding to the structural model, we examine its fit to the data. Evaluation of the structural model's fit to the data is not a simple procedure (Raykov et al., 1991), and there are no universally accepted fit indices (Raykov et al., 1991;Saris et al., 2009;Thoemmes et al., 2018;Weston & Gore, 2006). To assess our model's fit to our data, we examine both global and local fit indices. We do not pay attention to the significance of the χ 2 statistic, as in large samples (N > 500) the χ 2 is affected by sample size (Martens, 2005;Raykov et al., 1991;Weston & Gore, 2006). Furthermore, χ 2 test statistics and global fit indices are asymmetrically sensitive to different misspecifications (Saris et al., 2009), thus we also look at the local fit indices; expected parameters change (EPC) and modification index (MI) test (Saris et al., 1987). These tests are done with examining the power of the EPC's and significance of MI's; whereas the combination of high-power EPC's and nonsignificant MI's indicate no misspecification, and the combination of low-power EPC's and significant MI's indicate the misspecification of the parameter(s).
For our SEM, we provide both unstandardized and standardized estimates of coefficients and errors. Unstandardized estimates do not depend on equal variances from our specific sample, therefore they serve as more generalizable estimates of the relationships (Grace & Bollen, 2005) and interpret the unstandardized coefficients of the estimates and errors of the hypothesized measurement model.  Table 1 for detailed summary of the variables 1 3

Results
In this section we provide results from our descriptive statistics, correlations matrix, CFA, and SEM. From Table 1, we see that possible scores for this scale ranged from 0 to 6, with higher scores indicating greater IH, and MSM in our sample had a mean score of 1.34 (SD = 1.22). For SRB, out of 3,694 MSM 59.7% (N = 2,205) reported no no-steady partners that they had condomless sex with, while 0.5% (N = 17) reported more than 50 partners in the last twelve months. In the HIV/PrEP additive scale, out of 3,838 MSM in Spain in our sample 432 respondents has acquired a score of 10, and 3 has acquired a score of 0; the mean score was 7.15. Almost half of respondents (46.2%, N = 1,717) reported no SUIS in the last twelve months, while a small portion reported almost all of it and all of it (4.74%, N = 176).
In Table 2 correlational analysis among main study variables is presented in order to examine the relationship between variables. Due to the nature of the variables used, some of them are highly correlated. For this reason, we use a p-value of 0.001 in order to reduce the possible Type I Error threat, following the study of Kashubeck-West and Szymanski (2008). From Table 2 we note that we did not find a relationship between IH and SRB, nor with each of the items of IH (except for the first item, IH1). On the other hand, we found that HIV/PrEP knowledge (r = 0.11) and SUIS (r = 0.18) were positively related to SRB, meaning that when these variables increased, SRB also increased. Similarly, more recent use of each substance was positively related to SRB, expect for alcohol. We will further test these relationships with using CFA and SEM.

Confirmatory Factor Analysis
As the first step in SEM, we test the model's fit to the data using a CFA. For both CFA and SEM, there are copious measures to test model fit to the data (Meuleners et al., 2003). As commonly used global fit indices and based on Raykov et al. (1991), Hu and Bentler (1999) and Martens (2005), we used the following; (a) the Comparative Fit Index (CFI); (b) the Tucker-Lewis Index (TLI); (c) the root mean square error of approximation (RMSEA), and (d) Standardized root mean squared residual (SRMR). The majority of studies suggest that values higher than 0.95 for CFI and TLI indicate good fit (Hu & Bentler, 1999;Weston & Gore, 2006); while others suggest that CFI and TLI > 0.90 indicate a good fit to data (Moonie et al., 2009), as these fit indicates are susceptible to factors such as estimators and complexity (Xia & Yang, 2019). We use the cut-off value set by Hu and Bentler (1999). Further, Hu and Bentler (1999) suggested that values of RMSEA and SRMR < 0.06 are acceptable, and these values are widely accepted (Weston & Gore, 2006).
The results of the CFA proved that the model is a good fit to the data; CFI = 0.95, TLI = 0.95, RMSEA = 0.044 (90% confidence interval [CI] for the RMSEA lower bound = 0.041 and upper bound = 0.047), and SRMR = 0.03. A summary of these (and also for SEM, see below) fit indices can be found in Table 3.

Structural Equation Modelling
The SEM output for these global fit indices suggested that the measurement model (Fig. 1) was a good fit to the data; CFI = 0.96, TLI = 0.95, RMSEA = 0.044 (90% confidence interval [CI] for the RMSEA lower bound = 0.041 and upper bound = 0.047), and SRMR = 0.03. All of the parameter estimates were significant. Further, expected parameters change (EPC) and modification index (MI) test's fit indices recommendations and cut-off values can be found in Table 3. All EPC's and MI's of the measurement model's parameters are not mis-specified and each value meets the criteria suggested by Saris et al., (1987Saris et al., ( , 2009. The unstandardized estimation results of SEM are presented in Fig. 2 and Table 4. In Fig. 2, we present the estimated SEM results, with upper numbers presenting the unstandardized path coefficients, numbers in brackets presenting standardized coefficients, and numbers in parenthesis presenting robust standard errors. In Table 4, we present the estimated results for the defined parameters (i.e. mediation effects) from the structural model.
The direct path from IH to SRB (dashed line, Fig. 2) was statistically insignificant. Similarly, the defined parameter for the direct effect of IH to SRB (Table 4) was also insignificant. We did not find evidence for our first hypothesis, IH will be positively associated with SRB, at least not within a single country, in the specified SEM model.
As can be seen in Fig. 2, IH was negatively associated with HIV/PrEP knowledge, each unit increase in the IH latent variable was significantly associated with 0.39 decrease in the HIV/PrEP knowledge units. On the other hand, HIV/PrEP knowledge increases SRB, each unit increase in the HIV/PrEP knowledge was significantly associated with 0.04 increase in the SRB units. The indirect relationship between IH and SRB through HIV/PrEP knowledge was statistically significant  (Hu & Bentler, 1999;Weston & Gore, 2006) b (Saris et al., 1987(Saris et al., , 2009 (Table 4). An increase of 0.38 units in the IH latent variable was associated with 0.01 decrease SRB units, when mediated by the HIV/PrEP knowledge. Therefore, we found that as IH increased, the number of condomless intercourse with nonsteady partners reported (SRB) decreased; because higher IH decreases HIV/PrEP knowledge while higher HIV knowledge predicted increased SRB. Thus, we found  With respect to the relationship among IH, SUIS and SRB, we see that the relationship between IH and SUIS is insignificant, while the association between SUIS latent variable and SRB was positive and significant. Expectedly, Fig. 2 shows that a unit increase in SUIS was significantly associated with a 0.35 increase in the SRB units. The indirect path between IH and SRB mediated by SUIS was statistically insignificant (Table 4). Therefore, we did not find evidence for our third hypothesis that, sex under the influence of substances will fully mediate the relationship between IH and sexual risk behaviour.
In Fig. 2, we found that as IH decreased the recency of substance use. On the other hand, recency of substance use predicted SUIS significantly and positively. As can be seen in Table 4 (row 3), the indirect path from IH to SUIS when mediated by substance use recency, however, was statistically significant. As IH increased, SUIS decreased; because higher IH predicted less recent use of substances and more recent use of substances predicted higher numbers of sex under the influence of substances (SUIS). When looking at the variances, as expected, the role of alcohol in explaining the variances in SRB in comparison to rest of the substances is relatively small (0.032). The analysis has shown that particular substances in the SEM had greater influence than the rest, with 'ecstasy', regardless of its form, explaining the biggest share of the variances in the model (pill = 0.728; crystal = 0.701).
In another model (not shown here, available upon request), we added total number of partners (regardless of condom use or intercourse) in the last 12 months for controlling the HIV/PrEP knowledge latent variable. We found that the total number of partners is associated with increased HIV/PrEP knowledge significantly, and the rest of the results reported above stayed same or changed slightly. However, this model's fit indices were very different than the accepted criteria (CFI = 0.89, TLI = 0.87, RMSEA = 0.073, SRMR = 0.064). Therefore, we did not include this variable in our SEM model. Implications of this variable on our current results will be discussed.

Discussion
In this study, we investigated the relation between two variables, IH and SRB, and the possible mediating effects of HIV/PrEP knowledge and SUIS in a sample of non-PrEP using MSM living in Spain, who are HIV-negative or have been diagnosed with HIV but have detectable or unknown viral load. Contrary to our hypothesis that IH will be positively associated with SRB, we did not find any direct effect of IH on SRB. Some studies found that higher levels of IH leads to SRB (Folch et al., 2009;Puckett et al., 2017), while some studies did not find evidence for this direct relationship (Dawson et al., 2019;Puckett et al., 2017), including ours. Therefore, our study further contributes to the literature which suggest that there might be potential mediating variables that affect the relation between IH and SRB (Kashubeck-West & Szymanski, 2008a;Newcomb & Mustanski, 2010b). On one hand, in line with this and our expectations, results of our SEM analysis showed that HIV/PrEP knowledge mediated the relationship between IH and SRB. On the other hand, the results indicated that while SUIS was significantly associated with SRB, it did not mediate the relationship between IH and SRB.
We found that the more knowledgeable men are about how HIV is transmitted and PrEP works, the higher the numbers of condomless sexual intercourse with non-steady partners reported. One explanation for this surprising result might be that these men, in particular, take their sexual health seriously and are confident in knowing when to use condom and in their condom negotiation skills with nonsteady partners. For example, Klein (2013) found that condom use self-efficacy, which refers to an individual's self-confidence in their ability to use condoms, was significantly and positively associated with HIV knowledge. Therefore, these men may be more comfortable initiating conversation about how long ago they have been tested, whether their non-steady partner is using PrEP, or negotiate safeness that go into sex (i.e. knowing that no HIV transmissions from the HIV-positive partner to the HIV-negative partner would occur if their viral load is undetectable ("U = U; Undetectable = Untransmittable")). Our additional analysis (not shown in the paper, available upon request) has shown that the total number of partners may, in return, influence HIV/PrEP knowledge. Similarly, men who have a higher number of partners might have a better connection to the gay community and an increased chance of acquiring more knowledge about protecting themselves. Empirical theories about the link between behaviour and knowledge suggest that self-perceived vulnerability to HIV is probably the main factor underlying SRB, more so than knowledge (McKusik et al., 1985). That is, when people perceive that they are less vulnerable to HIV, they would be more likely to engage in potential sexual risk behaviour, independent of their knowledge about HIV. Moreover, HIV knowledge showed to be necessary, but not sufficient, to motivate individuals to avoid HIV-related risks (Pando et al., 2013). Therefore, HIV prevention programs should consider focusing on communicating what actually makes one less vulnerable to HIV, which is the knowledge that goes into the ability to negotiate safe sex.
We found that as IH increased, the number of condomless intercourse with nonsteady partners reported (SRB) decreased; because higher IH decreases HIV/PrEP knowledge while higher HIV knowledge predicted increased SRB. A high level of IH may serve as a barrier to participation in HIV testing and other health-promoting behaviours, and less contact with prevention and educational interventions. From previous studies we know that men with higher IH showed a reduced perception of their self-efficacy for condom use, even after intervention (Huebner et al., 2002). This finding may be indicative of certain men with high IH who are not confident and knowledgeable enough to communicate safer sex practices with non-steady partners. Alternatively, it is possible that the men with higher IH in our sample are less likely to be involved in sexual intercourse with non-steady partners in the first place. However, our conclusions are in consistency with the notion that IH is most likely harmful; higher levels of IH may lead to reduced number of non-steady partners, but it also reduces relevant knowledge to protect oneself , or is associated with condom use frequency in general, and is not directly associated with the number of non-steady partners. As discussed, with the development of PrEP the meaning of 'risky' sex may also change for MSM; with the number of partners not relating to risk if one knows how to protect oneself. Therefore, future research and prevention programs should consider different indicators for 'risky sex' in an aim to target those who not only have condomless sex with random partners, but those who do not know how to protect themselves.
MSM in our sample who were more prone to use alcohol or any other drug during sexual intercourse reported higher frequency of anal sex with non-steady partners without using condom. Provided they are, or their partner is seropositive (in the case of not knowing one's own or partner's HIV status), use of substances during sex can be of immediate relevance for risk of HIV exposure, for each individual. Further, higher IH was not associated with increased SUIS, nor was SUIS a variable through which IH influenced sexual risk behaviour. Connection to the gay community may promote unhealthy behaviours through submersion into a subculture that promotes drug use and provides easier access to drugs (Halkitis et al., 2005). High levels of IH may serve as a barrier to engagement with the gay community, thereby also serving as a barrier from community level factors that lead to greater substance use (Moody et al., 2018). These results emphasize the continuing importance of community level interventions that address substance use in the gay community.
Our results show that SUIS predicted the higher numbers of condomless intercourse with non-steady partners reported, without IH influencing SUIS. We found that substance use recency mediated the relationship between IH and SUIS. As IH increased, SUIS decreased; because higher IH predicted less recent use of substances and more recent use of substances predicted higher numbers of SUIS. When looking at the variances, as expected, the role of alcohol in explaining the variances in SRB in comparison to rest of the substances is relatively small (0.032). The analysis has shown that particular substances in the SEM had greater influence than the rest, with ecstasy explaining the largest share of the variances in the model (pill = 0.728; crystal = 0.701)-which is in line with how the question was asked (i.e. the question was asked as the most recent use of substances, and not the recent use of substances in amounts that would affect one's thinking or decision-making).

Strength and Limitations
Our study has several methodological strengths, including being the largest dataset of MSM living in Spain and use of SEM. We also note its limitations. We use recency time formats (when did you last…) for substance use variables, which further reduces the chance of recall bias-unlike frequency formats, which is not natural for most people. Similarly, there is no recall bias in knowledge, or the proportion of sex under the influence of substances. Accurately reporting partner numbers is generally a problem; but that does not affect our conclusions. Another strength of this survey lies in its anonymous character, through which the risk of social desirability bias is reduced (as opposed to the interview setting). While we used a large, diverse sample of MSM, the data come from a non-probability sample, potentially limiting generalizability, especially to those who are older, have lower education, fewer LGBT community attachments, or are more likely to conceal their sexual orientation (Prah et al., 2016). Yet, our analyses assume that the distribution of variables in the EMIS-2017 sample matches the distribution of these variables in the population. Further, non-probability sampling can also lead to higher estimates of sexual risk, drug use, or knowledge among MSM. These concerns are somewhat attenuated given that the present study was not focused on establishing population estimates or risk behaviours, HIV/PrEP knowledge, or IH, but instead sought to examine associations among variables, for which non-probability sampling is more appropriate (Meyer & Wilson, 2009). Further, probability-based studies typically include relatively small numbers of sexual minorities in one country only, and thus would not have provided an adequate sample size across numerous countries with which to evaluate our research aims.

Conclusion
The impacts of IH on the sexual risk behaviours of gay, bisexual and other MSM have been extensively studied, and we extend knowledge of the nature of the relation by attempting to disentangle the potential influence of HIV-related knowledge and substance use. Our SEM results suggest that IH is not directly implicated in the path to SRB and that HIV/PrEP knowledge, but not SUIS, mediate the relationship between IH and SRB. Future prevention strategies should also target specific counselling for MSM with low IH, and who are relatively knowledgeable about HIV risks and how PrEP works, in order to ensure that they are included within the prevention messages. Similarly, future interventions should consider addressing particular problems at the community level, such as substance use in general, sex under the influence of substances, and social homophobia which is exercised structurally against individuals.