Trends in hepatitis C virus seroprevalence and associated risk factors among men who have sex with men in Montréal: results from three cross-sectional studies (2005, 2009, 2018) ================================================================================================================================================================================== * Charlotte Laniece Delaunay * Joseph Cox * Marina B. Klein * Gilles Lambert * Daniel Grace * Nathan Lachowsky * Mathieu Maheu-Giroux ## Abstract **Introduction** To eliminate the hepatitis C virus (HCV) by 2030, Canada must adopt a micro-elimination approach targeting priority populations, including gay, bisexual, and other men who have sex with men (MSM). HCV prevalence and risk factors among MSM populations are context-dependent, and accurately describing these indicators at the local level is essential if we want to design appropriate, targeted prevention and treatment interventions. We aimed first to estimate and investigate temporal trends in HCV seroprevalence between 2005-2018 among Montréal MSM, and then to identify the socio-economic, behavioural, and biological factors associated with HCV exposure among this population. **Methods** We used data from three bio-behavioural cross-sectional surveys conducted among Montréal MSM in 2005 (n=1,795), 2009 (n=1,258), and 2018 (n=1,086). To ensure comparability of seroprevalence estimates across time, we standardized the 2005 and 2009 time-location samples to the 2018 respondent-driven sample. Time trends overall and stratified by HIV status, history of injection drug use (IDU), and age were examined. Modified Poisson regression analyses with generalized estimating equations were used to identify factors associated with HCV seropositivity pooling all surveys. We used multiple imputation by chained equations for all missing values. **Results** Standardized HCV seroprevalence among all MSM remained stable from 7% (95% confidence interval (CI): 3-10%) in 2005, to 8% (95%CI: 1-9%) in 2009, and 8% (95%CI: 4-11%) in 2018. This apparent stability hides diverging temporal trends in seroprevalence between age groups, with a decrease among MSM <30 years old, and an increase among MSM aged ≥45 years. History of IDU was the strongest predictor for HCV seropositivity (adjusted prevalence ratio: 8.0; 95%CI: 5.5-11.5), and no association was found between HCV seroprevalence and the sexual risk factors studied (condomless anal sex with men of serodiscordant/unknown HIV status, number of sexual partners, and group sex), nor with biological markers of syphilis. **Conclusions** HCV seroprevalence remained stable among Montréal MSM between 2005-2018. Unlike other settings where HCV infection was strongly associated with sexual risk factors among MSM subgroups, IDU was the preeminent risk factor for HCV seropositivity. Understanding the intersection of IDU contexts, practices, and populations is essential to prevent HCV transmission among MSM. Keywords * hepatitis C elimination * gay, bisexual, and other men who have sex with men * bio-behavioural surveillance * prevalence study * key population * respondent-driven sampling ## Introduction Canada endorsed the *World Health Organization*’s (WHO) global prevention and treatment service coverage targets to eliminate the hepatitis C virus (HCV) as a public health threat by 2030 [1]. With approximately 250,000 individuals chronically infected with HCV [2], Canada is not on track for HCV elimination [3]. The micro-elimination approach aims to develop tailored interventions to achieve the WHO’s targets among priority populations disproportionately affected by HCV [4, 5]. Due to exposures to HCV-contaminated blood and blood products before the introduction of HCV screening in 1990, HCV prevalence is high among individuals born between 1945-1975 in Canada [6]. Immigrants originating from HCV endemic countries may also be at increased risk of being chronically infected with HCV [7]. In the absence of behavioural risk factors such as sharing drug injection equipment, the risk of transmitting HCV is low among these groups. Contrarily, people who inject or use drugs, Indigenous peoples, people with experience in the prison system, and gay, bisexual, and other men who have sex with men (MSM) may be at high risk of acquiring and/or transmitting HCV [8]. Due to shared routes of transmission, HCV-HIV co-infection is common among individuals with similar vulnerabilities or risk behaviours. Since the early 2000s, HCV epidemics have been reported among MSM living with HIV in industrialized countries [9]. Global HCV incidence has more than doubled in this population, from 2.6/1,000 person-years before 2000 to 8.1/1,000 person-years since 2010 [10], with high reinfection rates [9]. Because HCV can be cleared, anti-HCV antibody detection in the blood –or seropositivity– is a marker of lifetime exposure to the virus, while detection of HCV RNA indicates an active infection. Global HCV seroprevalence among MSM living with HIV was estimated at 6.4% from 2002-2015 [11]. Comparatively, the burden of HCV has remained stable among HIV-negative MSM, with estimates of HCV seroprevalence ranging from 0.0-3.4% [12-16], and a global incidence of 0.4/1,000 person-years [17]. Nevertheless, recent cases of acute HCV infection have been reported among HIV-negative MSM without a history of injection drug use (IDU) in Europe [18], the United States [19], and Vancouver, British Columbia (Canada) [20]. These data may mask wide variations across geographical settings, highlighting the importance of local contexts and related risks. In industrialized countries, most HCV infections have been attributed to the use of unsterile needles and drug paraphernalia, and receipt of blood or blood products before the introduction of HCV screening [21]. In parallel, evidence from Europe, the United States, and Vancouver shows that sexual transmission can account for a substantial proportion of incident HCV infections among MSM [9, 20, 22-24]. Socio-economic factors such as education and income can also influence the risk of acquiring HCV in this population [22, 25-27]. The impact of socio-economic factors and injecting and sexual behaviours on HCV transmission among MSM subgroups is context-dependent and, to date, there is little evidence of sexual transmission of HCV among MSM in Canada [15, 16, 20, 28, 29]. It is estimated that between 312,681-426,384 MSM live in Canada, representing approximately 2.6% of the adult male population [30]. The relatively large size of this key population, the variability in HCV incidence and prevalence among MSM populations, as well as the diversity of potential routes of transmission reflect the importance of using locally-valid epidemiological evidence to inform tailored interventions targeting micro-elimination. Scarce evidence exists on HCV prevalence and risk factors among MSM subgroups in Canadian cities. In Vancouver, British Columbia, HCV seroprevalence among all MSM was estimated at 4.9% in a cross-sectional study conducted in 2008-2009 [15]. In this city, 16.8% of MSM living with HIV tested HCV seropositive, against 10.4% in a similar study conducted in Toronto, Ontario in 2010-2012 [16]. Temporal trends in HCV seroprevalence among MSM have not been assessed in Canada, however, available cross-sectional data suggest potential within-country variability in seroprevalence. In this study, we aimed to 1) investigate temporal trends in HCV seroprevalence among MSM living in Montréal from 2005-2018, overall and stratified by HIV status, IDU status, and age group; and 2) identify the social, behavioural, and biological factors associated with HCV exposure among this heterogeneous population. ## Methods ### Data sources and harmonization Three cross-sectional surveys of HIV and other sexually transmitted and blood-borne infections were conducted among Montréal MSM in 2005, 2009, and 2018: Argus 1, Argus 2, and the first wave of the Engage cohort study, respectively. These surveys have been previously described [31, 32]. Briefly, they combined HIV, HCV, and syphilis surveillance and behavioural monitoring through self-administered questionnaires. Time-location sampling was used to recruit 1,957 MSM in 2005 and 1,873 in 2009. In 2018, a community-based sample of 1,179 MSM was recruited using respondent-driven sampling (RDS), a link-tracing, adaptive sampling method that uses information about the social networks of participants to estimate their probability of being recruited and to adjust for the biases associated with over/under-sampling of certain groups [33]. Eligibility criteria varied slightly across the three surveys (Additional file 1). To compare the three surveys, we first harmonized the eligibility criteria across datasets. Our analyses were restricted to cis-gender men aged ≥18 years who reported sexual activity with a man in the past six months (P6M) and resided on the island of Montréal, Québec (Canada). The sensitivity and specificity of the assays used for anti-HCV antibody detection in the three surveys were comparable (sensitivity: 100% in 2005, 2009, and 2018; specificity: 99.95% in 2005 and 2009, and 99.69% in 2018 [31, 32, 34]). ### Temporal trends in HCV seroprevalence We estimated HCV seroprevalence in 2005, 2009, and 2018, among all MSM, stratified by HIV status, lifetime history of IDU, and age group, and among HIV-negative MSM without a history of IDU. To ensure comparability of the samples across time, direct standardization was performed. The RDS-weighted 2018 survey was used as the standard and the other surveys were adjusted to yield the same distribution with regards to age, sexual orientation, annual income, and first language. This was achieved by first fitting multivariable logistic regression models with HCV seropositivity as the outcome (defined as a positive anti-HCV antibody test) to the 2005 and 2009 data, separately. Age was modeled using a restricted cubic spline [35], and we investigated the presence of multiplicative statistical interaction between each pair of covariates. For the stratified analyses, the same model was fitted separately among each subgroup of interest. Using these models, we then predicted the individual probability of testing HCV seropositive using the covariate patterns observed in the 2018 standard sample, among all MSM and by subgroup. Finally, we computed the weighted mean of these individual probabilities using Volz-Heckathorn weights [36] to estimate standardized HCV seroprevalence among the different groups at all time points. We obtained 95% confidence intervals (CIs) using cluster-level block bootstrap, with clusters corresponding to the recruitment site in 2005 (n=40 sites) and 2009 (n=39 sites) and to the initial seed of the recruitment chain in 2018 (n=27 seeds). We performed complete case analyses as only 5% of observations had any missing value for the variables used in these estimations [37]. ### Social, behavioural, and biological factors associated with HCV seropositivity Pooling data from the three surveys, we generated prevalence ratios as measures of the effect of selected potential determinants of HCV seropositivity. We used univariable and multivariable modified Poisson regression models with generalized estimating equations, specifying exchangeable correlation structures to allow robust variance estimation and take clustering into account [38-40]. The following predictors of HCV seropositivity were selected a priori based on their potential association with HCV: year of data collection, lifetime history of IDU, HIV seropositivity, history of syphilis (defined as a reactive serological test), age (categorical), birth outside of Canada, annual income ≥30,000 Canadian dollars (CAD), education level higher than high school, sexual orientation other than gay/homosexual, self-identification as Indigenous, first language other than French or English, and specific sexual practices reported in the past P6M: transactional sex (defined as having given/received money, drugs, or other goods or services in exchange for sex), condomless anal sex (CAS) with a man of serodiscordant/unknown HIV status, >5 male sexual partners, and group sex. We used multiple imputation by chained equations for missing values as 33% of observations had at least one missing value for the variables used in these analyses, principally information regarding (23%) the question related to CAS [41, 42]. We conducted additional analyses to investigate whether the association of sexual practices with HCV seropositivity varied by HIV status, and to quantify the joint association of lifetime IDU and sexual behaviours with the outcome. We assessed the presence of additive and multiplicative interactions between the following terms: i) HIV seropositivity and four sexual risk factors reported in the P6M (i.e., transactional sex, CAS with a man of serodiscordant/unknown HIV status, number of sexual partners, and group sex) and ii) history of IDU and the same four factors. To evaluate multiplicative interaction, we included the product term of each selected pair of factors in the main analysis model separately and examined the estimate and CI of each product term coefficient [43]. Additive interaction was assessed by transforming the parameter estimates to obtain the relative excess risk due to interaction (RERI) using the method developed by Hosmer and Lemeshow [44]. CIs for the RERIs were computed following Zou [45]. All analyses were conducted using R 3.5.3 [46], with the *RDS* [47], *geepack* [48], *mice* [49], and *dplyr* [50] packages. ### Ethical considerations For each survey, written informed consent was obtained from all participants. A $10 CAD compensation was offered to participants of the two Argus surveys. In Engage, both participation ($50 CAD) and recruitment of peers ($15 CAD per recruit) were incentivized. The research ethics board of the *Research Institute of the McGill University Health Center* approved all three surveys. ## Results ### Characteristics of participants (pre-standardization) A total of 4,139 MSM were included in our analyses (1,795 in 2005, 1,258 in 2009, and 1,086 in 2018) (Figure 1). Before standardization, the compositions of the three samples differed. The median age was higher in 2005 and 2009 than in 2018 (Table 1). The proportion of men born outside of Canada was twice as low in the 2005 and 2009 samples as in the 2018 sample. The 2009 sample further contrasted with the others regarding social factors, with participants reporting a higher income, a lower level of education, and being more likely to identify as gay/homosexual. The proportion of participants living with HIV increased from 13% (95%CI: 11-14%) in 2005, to 15% (95%CI: 13-17%) in 2009, and 18% (95%CI: 16-21%) in 2018. We observed no clear trend in the percentage of participants with a history of IDU: 7% in 2005 (95%CI: 6-8%), 12% (95%CI: 10-14%) in 2009, and 10% (95%CI: 8-12%) in 2018. However, the proportions of MSM who had a reactive serological test for syphilis, reported recent CAS with a man of serodiscordant/unknown HIV status, and reported >5 male sexual partners in the P6M increased over time. Conversely, recent engagement in transactional sex decreased over the study period. Across the three surveys, 61% of participants who reported transactional sex gave, 62% received, and 23% both gave and received money, drugs, or other goods and services in exchange for sex. The proportion of MSM who engaged in group sex in the P6M increased from 2005-2009, and then decreased in 2018. View this table: [Table 1.](http://medrxiv.org/content/early/2020/01/28/2020.01.27.20018994/T1) Table 1. Description of participants included in three cross-sectional bio-behavioral surveys of men who have sex with men conducted in Montréal, Québec (Canada; 2005-2018; before any imputation or standardization). ![Figure 1.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/01/28/2020.01.27.20018994/F1.medium.gif) [Figure 1.](http://medrxiv.org/content/early/2020/01/28/2020.01.27.20018994/F1) Figure 1. Study eligibility flowchart. Individuals missing values on any eligibility criterion were excluded from our analyses. ### Temporal trends in HCV seroprevalence (post-standardization) Standardized HCV seroprevalence among all MSM remained stable and ranged from 7% (95%CI: 3-10%) in 2005, to 8% (95%CI: 1-9%) in 2009, and 8% (95%CI: 4-11%) in 2018 (Figure 2, Additional file 2). ![Figure 2.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/01/28/2020.01.27.20018994/F2.medium.gif) [Figure 2.](http://medrxiv.org/content/early/2020/01/28/2020.01.27.20018994/F2) Figure 2. HCV seroprevalence estimates among gay, bisexual, and other men who have sex with men (MSM) overall, stratified by (a) HIV status, (b) lifetime history of injecting drug use (IDU), and (c) age group, and among (d) HIV-negative MSM without a history of IDU in 2005, 2009, and 2018 in Montréal, Québéc (Canada). Stratification by HIV and IDU status did not reveal consistent temporal trends in HCV seroprevalence among these subgroups. We observed a drop in estimates for MSM living with HIV and MSM with a history of IDU in 2009 and obtained relatively wide CIs among these two groups due to smaller numbers of observations. The highest HCV seroprevalence point estimates were obtained among participants with a history of IDU, with 74% (95%CI: 62-85%) in 2005, 37% (95%CI: 8-48%) in 2009, and 60% (95%CI: 42-77%) in 2019. In contrast, HCV seroprevalence among MSM without a history of IDU was estimated at 3% (95%CI: 1-4%) in 2005, 2% (95%CI: 0-3%) in 2009, and 2% (95%CI: 1-4%) in 2018. HCV seropositivity was also far more common among MSM living with HIV at 17% (95%CI: 6-25%) in 2005, 8% (95%CI: 2-16%) in 2009, and 17% (95%CI: 6-29%) in 2018, than among HIV-negative MSM at 6% (95%CI: 3-8%) in 2005, 8% (95%CI: 0-9%) in 2009, and 6% (95%CI: 3-9%) in 2018. Finally, HCV seroprevalence remained relatively stable among HIV-negative MSM without a history of IDU at 2% (95%CI: 1-3%) in 2005, 2% (95%CI: 0-4%) in 2009, and 1% (95%CI: 0-2%) in 2018. When stratifying by age group, we observed a decrease in HCV seroprevalence among MSM <30 years old over time, from 6% (95%CI: 2%-9%) in 2005, to 3% (95%CI: 0-3%) in 2009, and 0% (95%CI: 0-1%) in 2018. We noted no clear temporal trend in HCV seroprevalence among MSM aged 30-44 years (9%; 95%CI: 3-12% in 2005, 10%; 95%CI: 1-12% in 2009, and 8%; 95%CI: 2-15% in 2018). Among MSM aged ≥45 years, HCV seroprevalence increased from 6% (95%CI: 1-9%) in 2005 to 9% (95%CI: 0-12%) in 2009 and to 17% (95%CI: 9-24%) in 2018. ### Factors associated with HCV seropositivity Lifetime history of IDU was the strongest predictor of HCV seropositivity in our population (Table 2). Increased age, recent transactional sex, HIV seropositivity, and sexual orientation other than gay/homosexual were also associated with increased HCV seroprevalence. Our results suggest an association between socio-economic factors and the outcome: MSM with higher income and education level were less likely to test HCV-seropositive. Being born outside of Canada was also associated with lower HCV seroprevalence. Few participants identified as Indigenous (1-2% across surveys), leading to wide uncertainty around the aPR for this population group. We found no association between recent sexual behaviours studied and HCV seropositivity. Finally, first language other than French or English was not associated with a variation in the outcome. View this table: [Table 2.](http://medrxiv.org/content/early/2020/01/28/2020.01.27.20018994/T2) Table 2. HCV seroprevalence ratios obtained by pooling three cross-sectional surveys of men who have sex with men conducted in Montréal, Québec (Canada) (2005, 2009, 2018). ### Sensitivity analyses We detected negative additive and multiplicative interactions between HIV seropositivity and transactional sex in the P6M (Table 3). The presence of negative additive and multiplicative interactions means that the association with HCV seropositivity when both factors are present is smaller than, respectively, the sum and the product, of the individual associations of the two factors with the outcome. We also identified negative additive interactions between history of IDU and each one of the sexual risk factors studied, separately, and a positive multiplicative interaction between history of IDU and group sex in the P6M. A positive multiplicative interaction indicates that the association when both factors are present is stronger than the product of the individual associations. View this table: [Table 3.](http://medrxiv.org/content/early/2020/01/28/2020.01.27.20018994/T3) Table 3. Adjusted prevalence ratios of multiplicative interaction terms and relative excess risk due to interaction of selected potential risk factors for HCV seropositivity pooling three cross-sectional surveys of men who have sex with men in Montréal, Québéc (Canada; 2005, 2009, 2018). ## Discussion Standardizing and pooling three large cross-sectional surveys of MSM from 2005-2018, our results suggest that HCV seroprevalence has remained stable at approximately 8% in Montréal. Overall, the point estimates for HCV seroprevalence among all MSM are higher than those recently estimated in other North American settings: 5% in Vancouver [15], 4% in San Francisco, California [51], and 3% in New-York City [52]. Pre-standardization estimates of HCV seroprevalence among all MSM were also stable over the study period at approximately 5%, contrasting with the increase observed in the prevalence of HIV and biological markers of syphilis. Despite the overall stability in HCV seroprevalence through time, we observed diverging trends among age groups. Because Canadian public health authorities do not systematically collect HCV incidence data, monitoring changes in HCV seroprevalence among young adults has been proposed as an approach to gauge trends in HCV incidence in the general population [8]. The decrease observed in HCV seroprevalence among participants <30 years old may therefore indicate a decline in new HCV infections among MSM from 2005-2018. In contrast, the increase in HCV seroprevalence among MSM aged ≥45 years may reflect birth cohort effects (i.e., males born between 1945-1975 [2]), temporal changes in injecting and/or sexual behaviours among older MSM, or changes in inter-generational transmission of HCV over time. HCV seroprevalence was particularly high among MSM living with HIV, and even more so among MSM with a history of IDU. The drop in seroprevalence observed in 2009 among these groups could reflect temporal changes in injecting behaviours (e.g., increased number of new injectors that have never been exposed to HCV). However, we cannot rule out that these results could be due to remaining differences in the composition of the samples post-standardization. HCV seroprevalence among HIV-negative MSM in Montréal ranged from 6-8% between 2005-2018, on the high end of estimates from Canada [16], the United States [52], and Europe [13]. The seroprevalence estimated among MSM without a history of IDU (2-3%) was also twice as high as that estimated in Vancouver for the same group [15]. These numbers are greater than the 1% HCV seroprevalence estimated for the general Canadian population [2]. Yet, the seroprevalence observed among HIV-negative MSM without a history of IDU was stable in the range of 1-2%, comparable to that among the general population. The relatively high HCV seroprevalence observed among all MSM may therefore reflect frequent exposure to the virus among subgroups with particular vulnerabilities or risk behaviours such as MSM living with HIV, and MSM with a history of IDU. Reporting a history of IDU was the strongest predictor of HCV seropositivity among Montréal MSM over the study period. This is particularly concerning given that the proportion of MSM who reported lifetime IDU in Engage (6%, [32]) is twelve times higher than that reported among the general population in Montréal (0.5%, [53]). Aside from recent transactional sex, which encompasses having given/received drugs (or money, goods, or services) in exchange for sex and has been previously identified as an important risk factor for HCV among MSM in Canada [15], we found no evidence of an association between sexual behaviours and HCV seroprevalence. These results are consistent with those obtained in Toronto [16], and with a study conducted in 2001 that did not find evidence of sexual transmission of HCV among a cohort of Montréal MSM [28]. However, they contrast with emerging evidence from Vancouver [20] and with studies conducted in Europe [9] and the United States [17], where sexual risk factors were strongly associated with new HCV infections among MSM. One potentially important risk factor that we could not investigate is chemsex, a form of sexualized drug use that has been previously associated with a higher risk of HCV infection [54, 55]. Nevertheless, the positive multiplicative interaction observed between group sex and IDU could partly capture the association of chemsex practices with HCV seroprevalence, although we could not measure the concomitance of injecting and sexual behaviours. Altogether, higher socio-economic status (i.e., income and education) was associated with lower HCV seroprevalence. Our data did not allow us to specifically examine the relationship between birth in HCV endemic countries and HCV seropositivity. In this regard, interpreting the negative association between birth outside of Canada and HCV seroprevalence is challenging. Having a first language other than French or English, another factor used as a proxy for country of origin, was not associated with the outcome. Reporting a sexual orientation other than gay/homosexual (e.g. bisexual, heterosexual, queer) was positively associated with HCV seropositivity. In line with previous work, this finding suggests that sexual orientation and behaviours should be examined in conjunction in order to fully understand how sexually transmitted and blood-borne infections are distributed among MSM populations [56]. Our results should be interpreted considering the study’s main limitations. First, the use of cross-sectional data renders the interpretation of observed associations challenging due to a temporality bias. However, pooling the available cross-sectional data allowed us to produce robust, comparable estimates of HCV exposure across time, and to investigate the plausibility of different modes of transmission of the virus among these groups. Second, we could not investigate the relationship between lifelong sexual practices and HCV seropositivity due to data limitations. Yet, we were able to quantify the association of numerous recent sexual behaviours with HCV exposure among Montréal MSM, adjusting for important factors such as age. Strengths of this study include its large sample size, resulting from the pooling of three comparable survey instruments. No sampling framework is available for MSM populations, and the Engage survey used an RDS design to constitute a diverse sample of MSM that enabled us to adjust for potential sampling biases in the Argus time-location surveys. Additionally, recent epidemiological data on HCV exposure and determinants among MSM living in Montréal is lacking, and very few population-based studies have examined HCV prevalence trends among MSM subgroups globally. As such, our work provides a unique opportunity to fill these research gaps. ## Conclusions HCV seroprevalence remained high among Montréal MSM between 2005-2018, likely reflecting frequent exposure to the virus among MSM living with HIV and MSM with a history of IDU. Unlike other settings where HCV infection was strongly associated with sexual risk factors among MSM subgroups, history of IDU was the preeminent risk factor for HCV seropositivity among Montréal MSM. Understanding the intersection of IDU contexts, practices, and populations is essential to prevent HCV transmission among these overlapping populations. ## Data Availability The data referred to in the manuscript are not publicly available. ## Competing interests C.L.D., G.L., D.G., and N.L. have no conflicts of interest to declare. J.C. has received research grant funding from Gilead Sciences, Merck Canada, and ViiV Healthcare. He also received travel support and honoraria for advisory work from Gilead Sciences, Merck Canada and ViiV Healthcare. M.K. reports grants for investigator-initiated studies from ViiV Healthcare, Merck, and Gilead; research grants from Janssen; personal fees from ViiV Healthcare, Bristol-Myers Squibb, AbbVie, and Gilead, all outside the submitted work. M.M.-G. reports grant funding from Gilead Sciences. ## Authors’ contributions C.L.D., J.C., M.K., and M.M.-G. conceived and designed the study. J.C. and G.L. were involved in the study design and data collection of the Argus 1 and 2 surveys. J.C., G.L., D.G., and N.L. were involved in the study design and data collection of the Engage study. C.L.D. administered and processed the different databases. C.L.D. performed the statistical analyses with inputs from M.K. and M.M.-G. All authors contributed to results interpretation. C.L.D. drafted the manuscript and all authors critically reviewed it for important intellectual content. All authors approved the final version. ## Funding The Argus surveys were funded by the *Institut national de santé publique du Québec* and the *Public Health Agency of Canada*. Engage is funded by the *Canadian Institutes of Health Research*, the *Canadian Foundation for AIDS Research*, and the *Ontario HIV Treatment Network*. C.L.D. received a Ph.D. trainee fellowship from the *Canadian Network on Hepatitis C*. The *Canadian Network on Hepatitis C* is funded by a joint initiative of the *Canadian Institutes of Health Research* (NHC-142832) and the *Public Health Agency of Canada*. C.L.D. also received a doctoral training award from the *Fond de Recherche du Québec – Santé*. M.K. is supported by a *Tier I Canada Research Chair*. M.M.-G.’s research program is funded by a career award from the *Fonds de recherche du Québec – Santé*. The authors thank the *Réseau de recherche en santé des populations du Québec* (RRSPQ) for its contribution to the financing of this publication. ## Additional files Additional file 1: Eligibility criteria for participating in the Argus 1 (2005), Argus 2 (2009), and Engage (2018) studies conducted in Montréal, Québec (Canada). View this table: [Table4](http://medrxiv.org/content/early/2020/01/28/2020.01.27.20018994/T4) Additional file 2: Standardized HCV seroprevalence among all men who have sex with men (MSM) and stratified by HIV status, by injection drug use status, and by age group, in 2005, 2009, and 2018 in Montréal, Québec (Canada). View this table: [Table5](http://medrxiv.org/content/early/2020/01/28/2020.01.27.20018994/T5) ## Acknowledgements The principal investigators of Argus are Michel Alary, Chris Archibald, Joseph Cox, Louis-Robert Frigault, Marc-André Gadoury, Gilles Lambert, René Lavoie, Robert Remis, Paul Sandstrom, Cécile Tremblay, François Tremblay, and Jon Vincelette. We would like to thank all the men who agreed to participate in the Argus surveys, as well as the owners and managers of the establishments where the participants were recruited. We would also like to thank the gay community organizations in Montréal who supported this research. The principal investigators of Engage are Joseph Cox, Daniel Grace, Trevor Hart, Jody Jollimore, Nathan Lachowsky, Gilles Lambert, and David Moore. We would like to thank all the Engage participants, the community engagement committee members, and the affiliated community agencies. ## Footnotes * E-mail addresses of authors: CLD: charlotte.laniece{at}mail.mcgill.ca, JC: joseph.cox{at}mcgill.ca, MK: marina.klein{at}mcgill.ca, GL: gilles.lambert.ccsmtl{at}ssss.gouv.qc.ca, DG: daniel.grace{at}utoronto.ca, NL: nlachowsky{at}uvic.ca, MMG: mathieu.maheu-giroux{at}mcgill.ca ## List of abbreviations aPR : adjusted prevalence ratio CAD : Canadian dollars CAS : condomless anal sex CI : confidence interval HCV : hepatitis C virus HIV : human immunodeficiency virus MSM : gay, bisexual, and other men who have sex with men P6M : past six months RDS : respondent-driven sampling RERI : relative excess risk due to interaction RNA : ribonucleic acid WHO : World Health Organization * Received January 27, 2020. * Revision received January 27, 2020. * Accepted January 28, 2020. * © 2020, Posted by Cold Spring Harbor Laboratory The copyright holder for this pre-print is the author. All rights reserved. The material may not be redistributed, re-used or adapted without the author's permission. ## References 1. 1.World Health Organization. Combating hepatitis B and C to reach elimination by 2030. Advocay brief. World Health Organization; 2016. 2. 2.Trubnikov M, Yan P, Archibald C. Estimated prevalence of hepatitis C virus infection in Canada, 2011. Canada communicable disease report = Releve des maladies transmissibles au Canada. 2014;40(19):429–36. 3. 3.Waheed Y, Siddiq M, Jamil Z, Najmi MH. Hepatitis elimination by 2030: progress and challenges. World journal of gastroenterology. 2018;24(44):4959–61. 4. 4.Lazarus JV, Wiktor S, Colombo M, Thursz M. Micro-elimination - A path to global elimination of hepatitis C. Journal of hepatology. 2017;67(4):665–6. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jhep.2017.06.033&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F01%2F28%2F2020.01.27.20018994.atom) 5. 5.Lazarus JV, Safreed-Harmon K, Thursz MR, Dillon JF, El-Sayed MH, Elsharkawy AM, et al. The micro-elimination approach to eliminating hepatitis C: strategic and operational considerations. Seminars in liver disease. 2018;38(3):181–92. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1055/s-0038-1666841&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F01%2F28%2F2020.01.27.20018994.atom) 6. 6.Trubnikov M, Yan P, Njihia J, Archibald C. Identifying and describing a cohort effect in the national database of reported cases of hepatitis C virus infection in Canada (1991-2010): an age-period-cohort analysis. Canadian Medical Association journal open. 2014;2(4):E281–E7. 7. 7.Greenaway C, Thu Ma A, Kloda LA, Klein M, Cnossen S, Schwarzer G, et al. The seroprevalence of hepatitis C antibodies in immigrants and refugees from intermediate and high endemic countries: a systematic review and meta-analysis. PloS one. 2015;10(11):e0141715. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0141715&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26558905&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F01%2F28%2F2020.01.27.20018994.atom) 8. 8.The Canadian Network on Hepatitis C Blueprint Writing Committee and Working Groups. Blueprint to inform hepatitis C elimination efforts in Canada. Canadian Network on Hepatitis C; 2019. 9. 9.Hagan H, Jordan AE, Neurer J, Cleland CM. Incidence of sexually transmitted hepatitis C virus infection in HIV-positive men who have sex with men. AIDS. 2015;29(17):2335–45. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1097/QAD.0000000000000834&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26258525&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F01%2F28%2F2020.01.27.20018994.atom) 10. 10.Ghisla V, Scherrer AU, Nicca D, Braun DL, Fehr JS. Incidence of hepatitis C in HIV positive and negative men who have sex with men 2000–2016: a systematic review and meta-analysis. Infection. 2017;45(3):309–21. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F01%2F28%2F2020.01.27.20018994.atom) 11. 11.Platt L, Easterbrook P, Gower E, McDonald B, Sabin K, McGowan C, et al. Prevalence and burden of HCV co-infection in people living with HIV: a global systematic review and meta-analysis. The Lancet infectious diseases. 2016;16(7):797–808. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1473-3099(15)00485-5&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26922272&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F01%2F28%2F2020.01.27.20018994.atom) 12. 12.Han R, Zhou J, Francois C, Toumi M. Prevalence of hepatitis C infection among the general population and high-risk groups in the EU/EEA: a systematic review update. BioMed Central infectious diseases. 2019;19(1):655. 13. 13.Newsum AM, van Rooijen MS, Kroone M, Bruisten SM, Matser A, Hogewoning A, et al. Stable low hepatitis C virus antibody prevalence among HIV-negative men who have sex with men attending the sexually transmitted infection outpatient clinic in Amsterdam, 2007 to 2017. Sexually transmitted diseases. 2018;45(12):813–7. 14. 14.Falla AM, Hofstraat SHI, Duffell E, Hahné SJM, Tavoschi L, Veldhuijzen IK. Hepatitis B/C in the countries of the EU/EEA: a systematic review of the prevalence among at-risk groups. BioMed Central infectious diseases. 2018;18(1):79. 15. 15.Wong J, Moore D, Kanters S, Buxton J, Robert W, Gustafson R, et al. Seroprevalence of hepatitis C and correlates of seropositivity among men who have sex with men in Vancouver, Canada: a cross-sectional survey. Sexually transmitted infections. 2015;91(6):430–3. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoic2V4dHJhbnMiO3M6NToicmVzaWQiO3M6ODoiOTEvNi80MzAiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMC8wMS8yOC8yMDIwLjAxLjI3LjIwMDE4OTk0LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 16. 16.Remis RS, Liu J, Loutfy MR, Tharao W, Rebbapragada A, Huibner S, et al. Prevalence of sexually transmitted viral and bacterial infections in HIV-positive and HIV-negative men who have sex with men in Toronto. PloS one. 2016;11(7):e0158090. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0158090&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F01%2F28%2F2020.01.27.20018994.atom) 17. 17.Jin F, Matthews GV, Grulich AE. Sexual transmission of hepatitis C virus among gay and bisexual men: a systematic review. Sexual health. 2017;14(1):28–41. 18. 18.McFaul K, Maghlaoui A, Nzuruba M, Farnworth S, Foxton M, Anderson M, et al. Acute hepatitis C infection in HIV-negative men who have sex with men. Journal of viral hepatitis. 2015;22(6):535–8. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/jvh.12366&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25412826&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F01%2F28%2F2020.01.27.20018994.atom) 19. 19.Volk JE, Marcus JL, Phengrasamy T, Hare CB. Incident hepatitis C virus infections among users of HIV preexposure prophylaxis in a clinical practice setting. Clinical infectious diseases. 2015;60(11):1728–9. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/cid/civ129&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25694649&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F01%2F28%2F2020.01.27.20018994.atom) 20. 20.Lachowsky N, Stephenson K, Cui Z, Shurgold S, Grennan T, Wong J, et al. Prevalence and factors of HCV infection among HIV-negative and HIV-positive MSM. Conference on retroviruses and opportunistic infections. 2016. 21. 21.Word Health Organization. Global hepatitis report. Word Health Organization; 2017. 22. 22.Breskin A, Drobnik A, Pathela P, Chan C, Braunstein S, Bornschlegel K, et al. Factors associated with hepatitis C infection among HIV-infected men who have sex with men with no reported injection drug use in New York City, 2000-2010. Sexually transmitted diseases. 2015;42(7):382–6. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F01%2F28%2F2020.01.27.20018994.atom) 23. 23.Danta M, Brown D, Bhagani S, Pybus OG, Sabin CA, Nelson M, et al. Recent epidemic of acute hepatitis C virus in HIV-positive men who have sex with men linked to high-risk sexual behaviours. Aids. 2007;21(8):983–91. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1097/QAD.0b013e3281053a0c&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17457092&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F01%2F28%2F2020.01.27.20018994.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000246754100011&link_type=ISI) 24. 24.Witt MD, Seaberg EC, Darilay A, Young S, Badri S, Rinaldo CR, et al. Incident hepatitis C virus infection in men who have sex with men: a prospective cohort analysis, 1984-2011. Clinical infectious diseases. 2013;57(1):77–84. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/cid/cit197&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23532480&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F01%2F28%2F2020.01.27.20018994.atom) 25. 25.Schmidt AJ, Falcato L, Zahno B, Burri A, Regenass S, Mullhaupt B, et al. Prevalence of hepatitis C in a Swiss sample of men who have sex with men: whom to screen for HCV infection? BioMed Central public health. 2014;14:3. 26. 26.Saxton PJ, Hughes AJ, Robinson EM. Sexually transmitted diseases and hepatitis in a national sample of men who have sex with men in New Zealand. The New Zealand medical journal. 2002;115(1158):U106. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12362181&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F01%2F28%2F2020.01.27.20018994.atom) 27. 27.Chen YC, Wiberg KJ, Hsieh YH, Bansal A, Bolzan P, Guy JA, et al. Favorable socioeconomic status and recreational polydrug use are linked with sexual hepatitis C virus transmission among human immunodeficiency virus-infected men who have sex with men. Open forum infectious diseases. 2016;3(3):ofw137. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/ofid/ofw137&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.ncbi.nlm.nih.gov/pubmed/27703998&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F01%2F28%2F2020.01.27.20018994.atom) 28. 28.Alary M, Joly JR, Vincelette J, Lavoie R, Turmel B, Remis RS. Lack of evidence of sexual transmission of hepatitis C virus in a prospective cohort study of men who have sex with men. American journal of public health. 2005;95(3):502–5. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2105/AJPH.2003.020388&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15727984&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F01%2F28%2F2020.01.27.20018994.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000228511100034&link_type=ISI) 29. 29.Myers T, Allman D, Xu K, Remis RS, Aguinaldo J, Burchell A, et al. The prevalence and correlates of hepatitis C virus (HCV) infection and HCV-HIV co-infection in a community sample of gay and bisexual men. International journal of infectious diseases. 2009;13(6):730–9. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ijid.2008.11.015&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19186090&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F01%2F28%2F2020.01.27.20018994.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000271373400015&link_type=ISI) 30. 30.Joint United Nations Programme on HIV/AIDS. Men who have sex with men: Population size estimate. Joint United Nations Programme on HIV/AIDS; 2013. 31. 31.Public Health Agency of Canada. M-Track: Enhanced surveillance of HIV, sexually transmitted and blood-borne infections, and associated risk behaviours among men who have sex with men in Canada. Phase 1 report. Public Health Agency of Canada; 2011. 32. 32.Lambert G, Cox J, Messier-Peet M, Apelian H, Moodie E. Engage Montréal, Portrait de la santé sexuelle des hommes de la région métropolitaine de Montréal ayant des relations sexuelles avec des hommes, Cycle 2017-2018, Faits saillants. Direction Régionale de Santé Publique du CIUSSS du Centre-Sud-de-l’Île de Montréal; 2019. 33. 33.World Health Organization, Regional Office for the Eastern Mediterranean Region. Introduction to HIV/AIDS and sexually transmitted infection surveillance: module 4: introduction to respondent-driven sampling. World Health Organization; 2013. 34. 34.Alborino F, Burighel A, Tiller FW, van Helden J, Gabriel C, Raineri A, et al. Multicenter evaluation of a fully automated third-generation anti-HCV antibody screening test with excellent sensitivity and specificity. Medical microbiology and immunology. 2011;200(2):77–83. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00430-010-0171-0&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20865278&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F01%2F28%2F2020.01.27.20018994.atom) 35. 35.1. Harrell JFE Harrell FE. General aspects of fitting regression models. In: Harrell JFE, editor. Regression modeling strategies: with applications to linear models, logistic and ordinal regression, and survival analysis. Springer international publishing; 2015. p. 13–44. 36. 36.Volz E, Heckathorn DD. Probability based estimation theory for respondent driven sampling. Journal of official statistics. 2008;24(1):79. [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000268944900006&link_type=ISI) 37. 37.Little R, Rubin D. Complete-case and available-case analysis, incuding weighting methods. In Statistical analysis with missing data, Third Edition. 2019. 38. 38.Zou G. A modified poisson regression approach to prospective studies with binary data. American journal of epidemiology. 2004;159(7):702–6. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/aje/kwh090&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15033648&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F01%2F28%2F2020.01.27.20018994.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000220484900009&link_type=ISI) 39. 39.Barros AJ, Hirakata VN. Alternatives for logistic regression in cross-sectional studies: an empirical comparison of models that directly estimate the prevalence ratio. BioMed Central medical research methodology. 2003;3:21. 40. 40.Yelland LN, Salter AB, Ryan P. Performance of the modified Poisson regression approach for estimating relative risks from clustered prospective data. American journal of epidemiology. 2011;174(8):984–92. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/aje/kwr183&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21841157&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F01%2F28%2F2020.01.27.20018994.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000295679700014&link_type=ISI) 41. 41.White IR, Royston P, Wood AM. Multiple imputation using chained equations: Issues and guidance for practice. Statistics in medicine. 2011;30(4):377–99. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/sim.4067&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21225900&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F01%2F28%2F2020.01.27.20018994.atom) 42. 42.van Buuren S, Groothuis-Oudshoorn K. mice: Multivariate imputation by chained equations in R. 2011;45(3):67. 43. 43.VanderWeele Tyler J, Knol Mirjam J. A tutorial on interaction. Epidemiologic methods. 2014;3(1):33. 44. 44.Hosmer DW, Lemeshow S. Confidence interval estimation of interaction. Epidemiology. 1992;3(5):452–6. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1097/00001648-199209000-00012&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=1391139&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F01%2F28%2F2020.01.27.20018994.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1992JM78900012&link_type=ISI) 45. 45.Zou GY. On the estimation of additive interaction by use of the four-by-two table and beyond. American journal of epidemiology. 2008;168(2):212–24. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/aje/kwn104&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18511428&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F01%2F28%2F2020.01.27.20018994.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000257394800011&link_type=ISI) 46. 46.R Core Team. R: A language and environment for statistical computing. 3.5.3 ed: R Foundation for Statistical Computing; 2017. 47. 47.Handcock MS, Gile KJ, Fellows IE, Neely WW. RDS: Respondent-driven sampling. 0.9-0 ed 2019. 48. 48.Højsgaard S, Halekoh U, Yan J. geepack: Generalized estimating equation package. 1.2-1 ed 2016. 49. 49.van Buuren S, Groothuis-Oudshoorn K, Robitzsch A, Vink G, Doove L, Jolani S, et al. mice: Multivariable imputation by chained equations. 3.6.0 ed 2019. 50. 50.Wickham H, François R, Henry L, Müller K, RStudio. dplyr: A grammar of data manipulation. 0.8.3 ed2019. 51. 51.Facente SN, Grebe E, Burk K, Morris MD, Murphy EL, Mirzazadeh A, et al. Estimated hepatitis C prevalence and key population sizes in San Francisco: A foundation for elimination. PloS one. 2018;13(4):e0195575. 52. 52.Tieu H-V, Laeyendecker O, Nandi V, Rose R, Fernandez R, Lynch B, et al. Prevalence and mapping of hepatitis C infections among men who have sex with men in New York City. PloS one. 2018;13(7):e0200269. 53. 53. Institut de la Statistique du Québec. L’Enquête québécoise sur la santé de la population, 2008 : pour en savoir plus sur la santé des Québécois. Institut de la Statistique du Québec; 2010. 54. 54.Maxwell S, Shahmanesh M, Gafos M. Chemsex behaviours among men who have sex with men: A systematic review of the literature. The international journal on drug policy. 2019;63:74–89. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F01%2F28%2F2020.01.27.20018994.atom) 55. 55.Tomkins A, George R, Kliner M. Sexualised drug taking among men who have sex with men: a systematic review. Perspectives in public health. 2019;139(1):23–33. 56. 56.Everett BG. Sexual orientation disparities in sexually transmitted infections: examining the intersection between sexual identity and sexual behavior. Archives of sexual behavior. 2013;42(2):225–36. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s10508-012-9902-1&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22350122&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F01%2F28%2F2020.01.27.20018994.atom)