Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Molecular Epidemiology of HIV-1 Infection among Men who Have Sex with Men in Taiwan in 2012

  • Szu-Wei Huang,

    Affiliations Institute of Microbiology and Immunology, National Yang Ming University, Taipei 11221, Taiwan, Center for Infectious Disease and Cancer Research, Kaohsiung Medical University, Kaohsiung 80708, Taiwan

  • Sheng-Fan Wang,

    Affiliations Center for Infectious Disease and Cancer Research, Kaohsiung Medical University, Kaohsiung 80708, Taiwan, Department of Medical Laboratory Science and Biotechnology, Kaohsiung Medical University, Kaohsiung 80708, Taiwan

  • Ángel E. Cowó,

    Affiliation International Health Program, National Yang Ming University, Taipei 11221, Taiwan

  • Marcelo Chen,

    Affiliations Department of Urology, Mackay Memorial Hospital, Taipei 10449, Taiwan, Department of Cosmetic Applications and Management, Mackay Junior College of Medicine, Nursing and Management, Taipei 25245, Taiwan, School of Medicine, Mackay Medical College, New Taipei City 25245, Taiwan

  • Yu-Ting Lin,

    Affiliation Center for Infectious Disease and Cancer Research, Kaohsiung Medical University, Kaohsiung 80708, Taiwan

  • Chun-Po Hung,

    Affiliation Rainbow Queer Health and Culture Center, Living with Hope Organization, Taiwan Society of Preventive Medicine, Taipei 10084, Taiwan

  • Yi-Hsien Chen,

    Affiliation Center for Infectious Disease and Cancer Research, Kaohsiung Medical University, Kaohsiung 80708, Taiwan

  • Jyh-Yuan Yang,

    Affiliation Centers for Disease Control, Taipei 10050, Taiwan

  • Hung-Jen Tang,

    Affiliations Department of Medicine, Chi Mei Medical Center, Tainan 71067, Taiwan, Department of Health and Nutrition, Chia Nan University of Pharmacy and Science, Tainan 71710, Taiwan

  • Yi-Ming Arthur Chen

    arthur@kmu.edu.tw

    Affiliations Institute of Microbiology and Immunology, National Yang Ming University, Taipei 11221, Taiwan, Center for Infectious Disease and Cancer Research, Kaohsiung Medical University, Kaohsiung 80708, Taiwan, Department of Microbiology and Immunology, Institute of Medical Research and Institute of Clinical Medicine, College of Medicine, Kaohsiung Medical University, Kaohsiung 80708, Taiwan

Abstract

The number of men who have sex with men (MSM) infected with HIV-1 in Taiwan has increased rapidly in the past few years. The goal of this study was to conduct a molecular epidemiological study of HIV-1 infection among MSM in Taiwan to identify risk factors for intervention. Voluntary counseling program and anonymous testing were provided to patrons at 1 gay bar, 7 night clubs and 3 gay saunas in Taipei and New Taipei Cities in 2012. HIV-1 subtypes were determined using gag subtype-specific PCR and phylogenetic analysis by env sequences. Recent HIV-1 infection was determined using LAg-Avidity EIA. In-depth interviews and questionnaires were used to identify risk factors. The prevalence and incidence of HIV-1 among MSM in Taiwan were 4.38% (53/1,208) and 3.29 per 100 person-years, respectively. Of 49 cases genotyped, 48 (97.9%) were infected with subtype B and 1 with CRF01_AE (2%). Phylogenetic analysis of 46 HIV-1 strains showed that 25 (54.4%) subtype B strains formed 9 clusters with each other or with other local strains. The CRF01_AE case clustered with a reference strain from a Thai blood donor with bootstrap value of 99. Multivariate logistic regression analysis showed that risk factors associated with HIV-1 infection included use of oil-based solution as lubricant (vs. saliva or water-based lubricants, OR= 4.23; p <0.001); exclusively receptive role (vs. insertive role, OR= 9.69; p <0.001); versatile role (vs. insertive role, OR= 6.45; p= 0.003); oral sex (vs. insertive role, OR= 11.93; p= 0.044); times of sexual contact per week (2-3 vs. zero per week, OR= 3.41; p= 0.021); illegal drug use (OR= 4.12; p <0.001); and history of sexually transmitted diseases (OR= 3.65; p= 0.002). In conclusion, there was no new HIV-1 subtype or circulating recombinant form responsible for the increase of HIV-1 among MSM in Taiwan in 2012. Misuse of oil-based solution as lubricant is a new risk factor identified among MSM in Taiwan. The Taiwan’s Centers for Disease Control has created a video (www.youtube.com/watch?v=BinExvvOTMM&feature=iv&src_vid=BW81-PfmY3E&annotation_id=annotation_2436493705) to correct such misconception in its AIDS prevention campaign.

Introduction

According to the Joint United Nations Program on HIV/AIDS (UNAIDS) report, 35.3 million people worldwide were living with HIV in 2012 [1]. A recent review has shown that the global trend of HIV-1 infection among men who have sex with men (MSM) has continued to increase, especially in East Asia, Africa, and Russia [2]. By the end of 2012, there were 25,081 cases of HIV infection been reported to Taiwan’s Centers for Disease Control (CDC) and Taiwanese nationals accounted for 96.6% of those cases (Taiwan CDC report 2012; http://www.cdc.gov.tw/info.aspx?treeid=1f07e8862ba550cf&nowtreeid=6c5ea6d932836f74&tid=58114AFB86117F53). Among them, homosexual men accounted for 41.9% of cases, followed by injection drug users (IDUs, 27.6%), heterosexuals (20.3%) and bisexual men (8.3%). There was an explosive outbreak of circulating recombinant form (CRF) _07BC among Taiwanese IDUs from 2004 to 2006 [35]. However, the HIV epidemic took a turn in 2007 as HIV-1 infection among IDUs decreased rapidly due to the implementation of harm reduction programs [6] and infection among MSM kept increasing exponentially [7]. Therefore, special attention needs to be given to any new HIV-1 subtype or CRF transmitted to a Taiwanese MSM. However, studies of the HIV-1 subtypes in MSM in Taiwan after 2009 are lacking.

Risk factors for HIV infection in MSM include unprotected anal intercourse [811], illegal drug use [1113], a history of other sexually transmitted diseases (STDs) [12], condom use, multiple sex partners and role during anal intercourse [8, 11, 14, 15]. Previous studies have shown an association between high level of substance use, sexual risk behavior and risk of HIV infection in young MSM [1619]. A study by Hikada et al. reported that 45% of MSM ever used one type of drug and 19.6% ever used more than one type of drugs [20]. According to the San Francisco Young Men’s Study, 43% of the subjects had used multiple illegal drugs [21].

The objective of this report was to conduct a molecular epidemiological study of HIV-1 subtypes and CRFs among MSM in Taiwan in 2012. Viral subtypes and recombinants and information on participants’ behavioral characteristics were obtained in order to investigate possible associations between subtypes and potential risk factors for HIV infection among MSM in Taiwan.

Materials and Methods

Subjects

Before beginning the study, we contacted commercial gay venues (one gay bar, seven gay night clubs and three gay saunas) in Taipei City to inquire about their willingness to provide space for our study. Anonymous, self-completed questionnaires and blood samples were collected from patrons at these venues. The information asked in our questionnaires includes demographics, HIV testing history, illegal drug use and sexual risk behavior (questionnaire available online at http://cicar.kmu.edu.tw/index.php/zh-TW/). All study participants received pre-test counseling. All study participants also received post-test counseling when results were given. Those who were HIV antibody negative received information on HIV prevention and were encouraged to undergo follow-up every 3 months. Those found to be HIV antibody positive were provided with information on where to obtain clinical assessment and treatment. This study was conducted with the approval of the institutional review board of the Mackay Memorial Hospital, Taipei, Taiwan. Written informed consent was obtained from participants before the interview and testing.

Serologic tests for HIV-1, syphilis and recent HIV-1 infections

HIV diagnosis was determined with a recombinant HIV enzyme immunoassay (Murex Diagnostics Limited, Dartford, UK). Positive samples were confirmed by HIV Western blot 2.2 (Genelab Technologies, Inc., Singapore). Active syphilis infection was defined using the following criteria: seropositive in Treponema pallidum hemagglutination assay (FTI-SERODIA-TPPA, Fujirebio Taiwan Inc., Taoyuan, Taiwan) and antibody titer ≥ 160 in rapid plasma regain (Macro-Vue RPR Kits, Becton, Dickinson and Company, NJ, USA) [22]. Recent HIV-1 infection was determined using LAg-Avidity EIA (Sedia Biosciences Corporation Portland, OR, USA). This assay was performed according to the manufacturer’s instructions. Briefly, recent infections were determined by a calibrator-normalized OD value (ODn) below the cut off value of 1. The mean duration of recency (MDR) was 141 days according to the manufacturer’s description. The annual HIV incidence was calculated using a formula provided by the manufacturer [23].

HIV-1 subtyping and phylogenetic analysis

DNA was extracted from peripheral blood mononuclear cells (PBMC) using QIAmp blood extraction kits (QIAgen, Valencia, CA.). HIV-1 gag and env regions were used to determine the subtypes by nested multiplex polymerase chain reaction (PCR) and phylogenetic analysis, respectively. A gag subtype-specific PCR with specific primers modified from a method previously described by Wei, et al. was used to determine HIV-1 viral subtypes and/or recombinants [24]. Extracted DNA was amplified by nested PCR using the TaKaRa Ex Taq system (Takara Bio Inc, Japan) according to the manufacturer’s instructions. HIV-1 env (C2-V5) was obtained through direct sequencing of the PCR products using a DNA analyzer (ABI 3730, Applied Biosystems, Foster City, CA). The nucleotide sequences of the patients were aligned with different reference sequences from Los Alamos HIV Database (http://www.hiv.lanl.gov) by BioEdit [25]. Phylogenetic trees were constructed using the best fit nucleotide substitution model [26] and the Neighbor-joining (NJ) and maximum likelihood (ML) methods using MEGA 6.06 [27] Phylip 3.69 [28]. The substitution model TN93 + G model was used to calculate the evolutionary distance and was followed by bootstrap value analysis with 1000 replicates for NJ tree construction [29]. To determine whether or not the HIV-1 infection was associated with foreign countries, sequences representing isolates from different regions in the world, in particular Asian countries and local sequences from Taiwan were used as reference sequences in the analysis. The local sequences were from subjects who were MSM (15/34), heterosexual (16/34) or IDUs (3/34) collected in our lab in 2004–2012. The selection of a bootstrap value of 70% was based on other references as well as on our past experiences. According to a study reported by Hills et al., the selection of bootstrap value of 70 correlates to a 95% probability of a cluster being real [30]. Therefore, bootstrap values >70 indicated significance of the sequences clustered.

Statistical analysis

Pearson's χ2 test and Fisher’s exact test were performed in a univariate analysis to investigate differences in demographic data and consumer behavior patterns among MSM. A multivariate logistic regression analysis was performed to identify risk factors associated with HIV-1 infection. The multivariate analyses were performed using SAS version 9.4 (SAS Institute, Inc., Gary, NC).

Nucleotide sequence accession numbers

The HIV-1 env sequences of 34 local control cases and 46 MSM obtained in this study were deposited in GenBank with accession numbers KM668624—KM668657 and KM668658—KM668703, respectively. The 34 HIV-1 env sequences were identified in our previous studies of MSM, IDUs and heterosexuals since we used them as local controls in our phylogenetic tree analysis.

Results

Demographic characteristics

A total of 1,208 participants were recruited in this study. The mean age of participants was 28.6 years (range, 18–57 years). Most of the participants were single (94%), self-identified homosexuals (80.2%) and with educational levels ranging from college (66.1%) to graduate school (16.4%). Over half of the participants (54.5%) were office workers and 20.9% were students. Most of the participants were recruited from gay night clubs (53.6%) and gay bar (33.9%), and the remaining 12.6% from gay saunas (Table 1).

thumbnail
Table 1. Demographic data of patrons from different gay venues participated in this study.

https://doi.org/10.1371/journal.pone.0128266.t001

Prevalence and incidence rates of HIV and syphilis among MSM

Fifty-three participants were diagnosed with HIV‐1 infection in 2012. The overall HIV‐1 prevalence was 4.38%. The results of LAg-Avidity EIA showed that 16 of the 53 HIV-1 infected MSM were recent seroconverters and the resultant incidence was 3.29 per 100 person-years (Table 2). The prevalence and incidence were not statistically different between gay venues (p = 0.763 and 0.496 respectively). The prevalence of syphilis among MSM in our cohort was 2.15%.

thumbnail
Table 2. Prevalence and incidence rates of HIV-1 infection among MSM from different gay venues in Taiwan.

https://doi.org/10.1371/journal.pone.0128266.t002

Distribution of HIV-1 subtypes among MSM

HIV-1 subtypes were determined by gag subtype-specific PCR amplification, and confirmed by sequencing and phylogenetic analysis of the env gene. During the study period, all HIV‐1 positive samples were available for genotyping, and 92.5% (49/53) were successfully genotyped. In 4 of our patients the subtypes could not be determined due to the limited quantity of the specimens. Therefore, their subtypes were classified as indeterminate. The result of gag subtype-specific PCR showed that 97.9% (48/49) of patients were infected with subtype B and one patient from gay sauna was infected with CRF01_AE. Forty-six sequences were amplified successfully by env-specific PCR. The NJ distance method classified 97.8% (45/46) of patients as infected with HIV-1 subtype B. One HIV-1 patient was classified as CRF01_AE. In addition, three HIV-1 subtype B infected patients were genotyped successfully by gag subtype-specific PCR but not by env-specific PCR. However, the results of HIV-1 subtype genotyping showed a consensus between gag subtype-specific PCR and phylogenetic analysis of env region (Table 3).

Taken together, except for a patient from gay sauna who was infected with CRF01_AE, all other patients (97.9%, 48/49) were infected with HIV-1 subtype B.

The origin and dissemination pathways of HIV-1 subtypes among MSM

Phylogenetic analysis using env V3-V5 loop gene was conducted to elucidate the transmission relationship between Asian countries and Taiwanese reference sequences (local control sequences). Consistent tree topologies were observed in both the NJ and ML trees. Of these local controls, 22 patients were infected with subtype B and 12 patients were infected with CRF01_AE. The results showed that most subtype B strains from gay venues clustered among themselves or with local controls recruited from 2004 to 2012 from different regions of Taiwan (13 clusters with bootstrap value > 70). Only one of the 46 HIV-1 patients from gay saunas clustered in clade CRF01_AE (G249). This sequence significantly clustered with a reference strain from Thailand (accession number: JN860762) with a bootstrap value of 99%. G249 was a patron recruited from a gay sauna and reported himself as a homosexual (Fig 1). Therefore, we conclude that HIV-1 subtype B strains from gay bar, gay night club and gay sauna patrons were of local origin and no non-B subtypes circulated in the MSM population.

thumbnail
Fig 1. Phylogenetic analysis of HIV-1 strains identified from MSM in Taiwan in 2012.

Neighbor-Joining tree based on env (7077–7619 nucleotide residues of HXB2) nucleotide sequences. All branch lengths are drawn to scale. Red circle indicates MSM patients from a gay bar. Aqua circle indicates MSM patients from gay night clubs. Blue circle indicates MSM patients from gay saunas. Green triangle indicates local controls from Taiwan. Clusters with bootstrap value ≥ 70 were labeled from numbers 1 to number 9.

https://doi.org/10.1371/journal.pone.0128266.g001

The characteristics of clusters of MSM infected with HIV-1 are shown in Figs 1 and 2. Twenty-five subtype B strains formed 9 clusters with each other or with other local controls. When we analyzed the clusters by age and illegal drug use, the results showed that there was one cluster (No. 3) of two members, both of whom were young (between 25–34 years old) who used illegal drugs (Fig 2). In addition, there were two clusters (No. 1 and No.5) of young MSM, although not all reported illegal drug use. The rest of the clusters were composed of either young MSM with or without reported use of illegal drugs (No. 8) or included MSM of a wide age range with or without reported use of illegal drugs (No. 4, 7) (S1 Table).

thumbnail
Fig 2. Analyses of the characteristics of the clusters of the phylogenetic trees based on age or illegal drug use.

(A) Characteristics of the cluster based on age. Green circle indicates 15–24 years old. Blue circle indicates 25–34 years old. Brown circle indicates 35–44 years old. Red circle indicates data not available (NA). (B) Characteristics of the cluster based on illegal drug use. Yellow circle indicates illegal drug use. Purple circle indicates no illegal drug use. Red circle indicates data not available (NA). Clusters with bootstrap value ≥ 70 were labeled from numbers 1 to number 9 and number of individuals per cluster recruited in this study.

https://doi.org/10.1371/journal.pone.0128266.g002

Risk factors associated with HIV-1 infection among MSM

Univariate and multivariate risk factor analyses for HIV-1 infection are presented in Tables 4 and 5. According to univariate analysis, exclusively receptive role and versatile role (vs. exclusively insertive role, OR = 6.79 and 4.92; p = 0.003 and 0.009), times of sexual contact per week (2–3 vs. zero per week, OR = 3.34; p = 0.015), illegal drug use (vs. no use, OR = 4.28; p <0.001), use of 3 or more drugs (vs. single drug use, OR = 5.37; p = 0.002), use of oil-based lubricants during sexual intercourse (vs. saliva or water-based lubricants, OR = 3.18; p = 0.004), and history of STDs (vs. no history, OR = 4.24; p <0.001) were all associated with an increased risk of HIV-1 infection. Multivariate analysis was conducted of the risk factors mentioned above. The risk factors significantly associated with HIV-1 infection were exclusively receptive role, versatile role and oral sex (vs. exclusively insertive role, OR = 9.69, 6.45 and 11.93; p <0.001, 0.003 and 0.044, respectively), times of sexual contact per week (2–3 vs. zero per week, OR = 3.41; p = 0.021), illegal drug use (vs. no use, OR = 4.12; p <0.001), use of oil-based lubricants during sexual intercourse (vs. saliva or water-based lubricants, OR = 4.23; p <0.001) and history of STDs (vs. no history, OR = 3.65; p = 0.002) (Table 5).

thumbnail
Table 4. Univariate analysis of risk factors for HIV-1 infection among MSM in Taiwan.

https://doi.org/10.1371/journal.pone.0128266.t004

thumbnail
Table 5. Multivariate logistic regression of risk factors for HIV-1 infection among MSM in Taiwan.

https://doi.org/10.1371/journal.pone.0128266.t005

Illegal drug use among MSM

A total of 259 participants reported use of one or more types of illegal drugs. As shown in Table 6, the top three types of drugs abused by MSM were RUSH (Alkyl nitrites, 44.7%), MDMA (9.8%) and Ketamine (9.8%). HIV-1-infected MSM had significantly higher rate (17.9%, 5/28) of Ketamine, MDMA and Viagra use than HIV-1 negative MSM (5.2%, 12/231). The difference of the rates between these two groups was statistically significant (p = 0.011). It is worthy to note that there were two HIV-1-infected MSM reporting use of 6 drugs (Table 6).

thumbnail
Table 6. Distribution of different types or numbers of illegal drugs used among MSM in Taiwan.

https://doi.org/10.1371/journal.pone.0128266.t006

MSM using illegal drugs had a mean age of 30 years (range, 18–57 years), and most of them were single (93.1%), self-identified homosexuals (85.7%) and with educational levels ranging from college (62.9%) to graduate school (16.2%). Over half of the participants (58.7%) were office workers and 13.1% were students (Table 7). Univariate analysis showed that Ketamine, MDMA and Viagra use (vs. single drug use, OR = 6.46; p = 0.004) and use of oil-based lubricants during sexual intercourse (vs. saliva or water-based lubricants, OR = 4.18; p = 0.019) were significantly associated with HIV-1 infection (Table 8). Multivariate analysis demonstrated that the risk factors mentioned above remained significantly associated with HIV-1 infection in MSM who reported illegal drug use (Table 9). Furthermore, versatile role (vs. exclusively insertive role, OR = 6.38; p = 0.015) and Ketamine, MDMA and Nimetazepam use (vs. single drug use, OR = 60.12; p = 0.014) were also associated with HIV-1 infection in MSM using illegal drugs.

thumbnail
Table 7. Demographic data of patrons who reported illegal drug use from different gay venues participated in this study.

https://doi.org/10.1371/journal.pone.0128266.t007

thumbnail
Table 8. Univariate analysis of risk factors for HIV-1 infection among MSM who reported illegal drug use in Taiwan.

https://doi.org/10.1371/journal.pone.0128266.t008

thumbnail
Table 9. Multivariate logistic regression of risk factors for HIV-1 infection among MSM who reported illegal drug use in Taiwan.

https://doi.org/10.1371/journal.pone.0128266.t009

Discussion

In this study, we found that most HIV-1 seropositive MSM in Taiwan were infected with subtype B. In addition, phylogenetic analysis showed that they were clustered with local controls. Furthermore, we found that exclusively receptive, versatile roles and oral sex, 2–3 times of sexual contact per week and history of STDs, illegal drug use and use of oil-based solution as lubricants were important risk factors for HIV-1 infection among MSM in Taiwan.

The present study reported a 4.38% prevalence of HIV-1 in MSM in Taiwan in 2012. This figure differs from that of a similarly published report that indicated a variable but considerably higher prevalence of HIV-1 among MSM visiting gay saunas in Taiwan between 2001 and 2005 [14]. According to results of syphilis testing, 2.15% of participants in this study were infected with syphilis. Syphilis infections significantly increase the risk of HIV infection. Some studies have reported that after controlling for other risk factors, participants who had been infected with syphilis were more likely to be infected with HIV [31, 32].

The BED-CEIA is the most commonly used assay for estimating HIV-1 incidence [33]. However, several previous studies in different populations with different HIV-1 subtype infections showed that LAg-Avidity EIA had a lower false recent rate and higher accuracy than BED-CEIA [3436]. Therefore, we used LAg-Avidity EIA in our analysis. The results showed that 16 of 53 HIV-1 infected MSM were recent seroconverters and the resultant incidence was 3.29 per 100 person-years. False recent infection in HIV-1 seropositive patients can be due to antiretroviral therapy (ARV), elite controllers, terminal stage of AIDS, and non-B subtype infection (WHO/UNAIDS Technical Update on HIV incidence assays for surveillance and epidemic monitoring, 2013). However, none of the factors mentioned above are present in our study: 15 of 16 (93.8%) recent seroconverters detected in our study were treatment naïve patients with subtype B infection and the remaining one (6.3%) was a treatment naïve patient with CRF01_AE infection.

In this study, we found that 97.9% of all infections were caused by subtype B. This was not surprising given that previous reports had already established B as the most predominant subtype among Taiwanese MSM [14, 37, 38]. Of the non-subtype B infections, 2% was attributed to CRF01_AE. Although CRF01_AE is predominantly found among female heterosexuals, several studies have speculated on their possible transmission to other risk groups [5, 7, 39]. The results of HIV-1 subtyping were consistent between gag subtype-specific PCR (48 patients infected with subtype B and 1 patient infected with CRF01_AE) and phylogenetic analysis (45 patients infected with subtype B and 1 patient infected with CRF01_AE). The rationale for using both gag subtype-specific PCR and DNA sequencing of env for genotyping was because the former method can detect dual or even triple-subtype infection and the latter can be used for phylogenetic tree analysis [7, 39].

Phylogenetic analysis of viral gene sequences has successfully been used to construct direct or indirect epidemiological links in geographically defined populations with acute/primary or chronic HIV-1 infection [40, 41]. The epidemic of HIV-1 infection in MSM has expanded rapidly in Asian countries [42, 43]. Our previous study showed that CRF07_BC was transmitted into Taiwan from China and subsequently caused an epidemic in the Taiwanese IDU population [3, 5]. A previous study showed that CRF07_BC infection has spread to the heterosexual population [7]. In this study, 97.9% MSM were infected with subtype B and CRF07_BC did not transmit into this population. A recent study demonstrated that CRF07_BC had significantly lower viral loads than subtype B [44]. This may explain why CRF07_BC cannot transmit through sexual contact as efficiently as through needle sharing.

In this study, 13 clusters were identified in the phylogenetic analysis of 46 env sequences. In all subtype B clusters, MSM sequences aggregated among themselves or with the local control sequences. The CRF01_AE strain (G249) clustered with a reference sequence from Thailand with a bootstrap value of 99%. This sequence was from a Thai blood donor [45]. No further details about the Thai individual could be obtained to better assess the association between these two sequences. However, G249 claimed that he had joined a sex tour in Thailand. Since CRF01_AE is predominant among MSM in China and Thailand [4648], further surveillance of the prevalence of CRF01_AE among MSM in Taiwan is needed to monitor its transmission dynamics in this population.

Compared to heterosexual males, MSM are more susceptible to HIV-1 infection. Cultural, social, political and religious factors that lead to denial of MSM and their sexual practices increase their vulnerability [49]. These are important aspects that HIV prevention programs should take into account. MSM are also known to face additional challenges that may lead to sexual risk taking behavior and increase their chances for HIV seroconversion. This has been well-documented throughout the epidemic and confirmed in this study. Previous studies demonstrated that receptive and versatile roles during anal intercourse were associated with an increased risk of HIV seroconversion among MSM [32, 50]. Generally, the receptive partner was at greater risk of contracting HIV because the lining of the rectum is thin and may allow the virus to enter the body through semen exchange [32, 51]. In a previous study among MSM visiting gay saunas in Taiwan between 2001 and 2005, men assuming receptive and versatile roles during anal intercourse were also identified to be at higher risk of seroconversion [14].

Illegal drug use among MSM is one of the most cited risk factors associated with HIV-1 infection. In this study, the odds of HIV-1 seroconversion in MSM who used illegal drugs was 4.12 times higher than in those who did not use illegal drugs (OR = 4.12; P <0.001). The relationships between substance use, high-risk sexual behaviors and HIV infection are complex. It is believed that illegal drug use among MSM may be either a trigger or an excuse for engaging in unprotected sex [52, 53].

In this study, compared to illegal use of a single drug, alternate use of Ketamine and MDMA and either Viagra or Nimetazepam was associated with HIV-1 infection. Ketamine has powerful hallucinogenic qualities and makes the user feel disassociated from his or her body. The user may feel sleepy or sluggish, or confused and clumsy [54]. MDMA is a synthetic hallucinogenic stimulant. The immediate effect of MDMA is a feeling of euphoria. It makes the user feel more comfortable in social situations [55]. Viagra is an orally administered prescription medication to treat erectile dysfunction. Nimetazepam possesses hypnotic, anxiolytic, sedative, and skeletal muscle relaxant properties. MDMA is usually not taken with Ketamine and Nimetazepam due to the different effects of the drugs. Several studies has shown that use of two or more psychoactive drugs (polydrug use) places individuals at greater risk of engaging in unprotected sex and thus increases the risk of HIV seroconversion [5658]. However, participation in risky sexual behavior and risk of HIV seroconversion vary greatly by type of substances used. Polydrug use may increase the risk of HIV seroconversion since it can interfere with the adherence or the effectiveness of antiretroviral therapy [12, 59]. Therefore, having unprotected anal intercourse with an HIV-infected partner who frequently engages in polydrug use increases the probability of HIV acquisition as these persons may be highly infectious due to a reduction in the effectiveness of their treatment. Moreover, ingesting a combination of drugs may lead to transient immune suppression [60]. Several studies have also identified reasons for polydrug consumption among MSM, including the potentiating effect associated with mixing drugs [61].

Participants with a history of STDs were also more likely to be infected with HIV (OR = 3.65; p = 0.002). Research suggests that STDs can increase both an HIV-negative person’s risk of becoming infected with HIV and an HIV-positive person’s risk of transmitting HIV to someone else [50, 62]. Having an STD increases a person’s risk of becoming infected with HIV due to various biological mechanisms. STDs increase the concentration of dendritic cells and CD4 T-cells in the genito-anal region, which is the most usual site of STDs infection [63]. The presence of an STD in a co-infected person increases the chances of transmitting HIV through open genital lesions, which come into contact with the partner’s anal and oral mucous membranes during sexual activity [31].

The inappropriate use of lubricants among MSM is not a new phenomenon, as it was observed in high-income settings nearly two decades ago [64]. Surprisingly, we found a very high percentage of oil- or petroleum-based solution misuse as lubricants in this study (47.2% in HIV-1-infected MSM vs. 25.9% HIV-1 sero-negative MSM). Oil-based lubricants misused included Vaseline (14%), baby oil (6%), soap (4%) and other types of lotion (6%). Multivariate analysis showed that participants using such lubricants were more likely to be infected with HIV than those using saliva or water-based lubricants (OR = 4.23; p <0.001). In addition, data analysis combined with in-depth interview at post-test counseling showed that among 53 HIV-1 seroconverters, 12 (22.6%) MSM used condoms consistently but misused oil-based solution as lubricants. Since most of them played a receptive role (bottom) during anal intercourse, they were not aware of the breakage of condoms during or after sexual intercourse. This is a very important risk factor that has been neglected in AIDS education for a long time. Similar phenomena may exist among MSM in other countries and this should be taught in AIDS education campaign. The Taiwan’s Centers for Disease Control has produced a video to correct this misconception (www.youtube.com/watch?v=BinExvvOTMM&feature=iv&src_vid=BW81-PfmY3E&annotation_id=annotation_2436493705). Our previous study showed that uncircumcised, versatile role during anal intercourse, and having sex with more than one person during each sauna visit were main risk factors for HIV-1 infection in Taiwan from 2001 to 2005 [14]. A previous study by Ko et al showed that illegal drug use was the primary risk factor for HIV infection [12]. A study in China showed that the number of sex partners and frequency of anal sex were associated with HIV-1 infection among MSM [15]. A cohort study conducted between April 2006 and July 2012 in Thai MSM showed that drug use for sexual pleasure and receptive anal intercourse were associated with HIV-1 infection [65]. Risk factors that are consistently present in different studies were receptive role during anal intercourse and illegal drug use. Misuse of oil-based solutions as lubricant and illegal drug use are risk factors that have never being mentioned previously.

Overall, 4.38% of MSM participating in this study were infected with HIV. Among the infected participants, the majority were infected with subtype B. Multivariate analysis showed that risk factors of HIV-1 seroconversion among MSM included an exclusively receptive role during anal intercourse, multiple sexual contacts per week, use of illegal drugs, use of oil-based lubricants during sexual intercourse and history of STDs.

The burden of HIV among MSM in Taiwan is disproportionately high. The results of this study can be useful for the design and implementation of evidence-based interventions aimed at preventing and controlling HIV and STDs among MSM, supporting safe sexual behavior, promoting uptake of HIV testing, and promoting HIV and STD health care seeking behavior.

Supporting Information

S1 Fig. Maximum likelihood tree analysis based on Env nucleotide sequence of MSM in Taiwan in 2012.

https://doi.org/10.1371/journal.pone.0128266.s001

(TIF)

S1 Table. Characteristics of the clusters in phylogenetic tree analysis.

https://doi.org/10.1371/journal.pone.0128266.s002

(DOC)

Acknowledgments

We thank the MSM who participated in this study, peer educators and social workers from the Living with Hope Organization for their help in collecting the questionnaires, and staff of the AIDS Prevention and Research Center of National Yang-Ming University for their administrative support and technical assistance. The authors also acknowledge the High-throughput Genome Analysis Core Facility of National Core Facility Program for Biotechnology from National Yang-Ming University, Taiwan (MOST 103-2319-B-010-001), for sequencing.

Author Contributions

Conceived and designed the experiments: YMAC. Performed the experiments: SWH SFW AEC YTL. Analyzed the data: SWH YHC MC. Contributed reagents/materials/analysis tools: JYY HJT. Wrote the paper: SWH AEC. Recruited MSM: CPH.

References

  1. 1. Joint United Nations Programme on HIV/AIDS., United Nations. Global Report: UNAIDS Report on the Global AIDS epidemic: 2012. Geneva: UNAIDS; 2012. 103, A5 p. p.
  2. 2. van Griensven F, de Lind van Wijngaarden JW, Baral S, Grulich A. The global epidemic of HIV infection among men who have sex with men. Current opinion in HIV and AIDS. 2009;4(4):300–7. Epub 2009/06/18. pmid:19532068.
  3. 3. Chen YM, Kuo SH. HIV-1 in Taiwan. Lancet. 2007;369(9562):623–5. pmid:17321295.
  4. 4. Chen YM, Lan YC, Lai SF, Yang JY, Tsai SF, Kuo SH. HIV-1 CRF07_BC infections, injecting drug users, Taiwan. Emerging infectious diseases. 2006;12(4):703–5. pmid:16715583; PubMed Central PMCID: PMC3294689.
  5. 5. Lin YT, Lan YC, Chen YJ, Huang YH, Lee CM, Liu TT, et al. Molecular epidemiology of HIV-1 infection and full-length genomic analysis of circulating recombinant form 07_BC strains from injection drug users in Taiwan. The Journal of infectious diseases. 2007;195(9):1283–93. Epub 2007/04/03. pmid:17396997.
  6. 6. Lyu SY, Su LW, Chen YM. Effects of education on harm-reduction programmes. Lancet. 2012;379(9814):e28–30. pmid:21851974.
  7. 7. Chen YJ, Lee CM, Chen M, Chuang SY, Liu HF, Wong WW, et al. Molecular epidemiology of HIV-1 infection in Taiwan from 2005 to 2008: further spread of CRF07_BC and emergence of CRF07_BC/subtype B dual infection. J Acquir Immune Defic Syndr. 2012;59(5):438–46. Epub 2012/02/22. pmid:22343173.
  8. 8. Read TR, Hocking J, Sinnott V, Hellard M. Risk factors for incident HIV infection in men having sex with men: a case-control study. Sexual health. 2007;4(1):35–9. Epub 2007/03/27. pmid:17382036.
  9. 9. He Q, Xia Y, Raymond HF, Peng R, Yang F, Ling L. HIV trends and related risk factors among men having sex with men in mainland China: findings from a systematic literature review. The Southeast Asian journal of tropical medicine and public health. 2011;42(3):616–33. Epub 2011/06/29. pmid:21706940.
  10. 10. Piaseczna MA, Craib KJ, Li K, Chan K, Weber AE, Strathdee SA, et al. Longitudinal patterns of sexual behavior and condom use in a cohort of HIV-negative gay and bisexual men in Vancouver, British Columbia, Canada, 1995–2000. J Acquir Immune Defic Syndr. 2001;28(2):187–93. Epub 2001/10/06. pmid:11588514.
  11. 11. Weber AE, Chan K, George C, Hogg RS, Remis RS, Martindale S, et al. Risk factors associated with HIV infection among young gay and bisexual men in Canada. J Acquir Immune Defic Syndr. 2001;28(1):81–8. Epub 2001/10/02. pmid:11579281.
  12. 12. Ko NY, Lee HC, Hung CC, Tseng FC, Chang JL, Lee NY, et al. Trends of HIV and sexually transmitted infections, estimated HIV incidence, and risky sexual behaviors among gay bathhouse attendees in Taiwan: 2004–2008. AIDS and behavior. 2011;15(2):292–7. pmid:20614174.
  13. 13. Mayer KH, Skeer MR, O'Cleirigh C, Goshe BM, Safren SA. Factors Associated with Amplified HIV Transmission Behavior Among American Men Who Have Sex with Men Engaged in Care: Implications for Clinical Providers. Annals of behavioral medicine: a publication of the Society of Behavioral Medicine. 2013. Epub 2013/07/23. pmid:23873338.
  14. 14. Chen YJ, Lin YT, Chen M, Huang SW, Lai SF, Wong WW, et al. Risk factors for HIV-1 seroconversion among Taiwanese men visiting gay saunas who have sex with men. BMC infectious diseases. 2011;11:334. pmid:22136602; PubMed Central PMCID: PMC3295735.
  15. 15. Xu HL, Jia MH, Min XD, Zhang RZ, Yu CJ, Wang J, et al. Factors influencing HIV infection in men who have sex with men in China. Asian journal of andrology. 2013;15(4):545–9. pmid:23708455; PubMed Central PMCID: PMC3739247.
  16. 16. Crosby GM, Stall RD, Paul JP, Barrett DC. Alcohol and drug use patterns have declined between generations of younger gay-bisexual men in San Francisco. Drug and alcohol dependence. 1998;52(3):177–82. Epub 1998/12/05. pmid:9839143.
  17. 17. Drumright LN, Little SJ, Strathdee SA, Slymen DJ, Araneta MR, Malcarne VL, et al. Unprotected anal intercourse and substance use among men who have sex with men with recent HIV infection. J Acquir Immune Defic Syndr. 2006;43(3):344–50. Epub 2006/09/19. pmid:16980913.
  18. 18. Drumright LN, Gorbach PM, Little SJ, Strathdee SA. Associations between substance use, erectile dysfunction medication and recent HIV infection among men who have sex with men. AIDS and behavior. 2009;13(2):328–36. pmid:18064558.
  19. 19. Wilson PA, Cook S, McGaskey J, Rowe M, Dennis N. Situational predictors of sexual risk episodes among men with HIV who have sex with men. Sexually transmitted infections. 2008;84(6):506–8. Epub 2008/11/26. pmid:19028956; PubMed Central PMCID: PMC3645477.
  20. 20. Hidaka Y, Ichikawa S, Koyano J, Urao M, Yasuo T, Kimura H, et al. Substance use and sexual behaviours of Japanese men who have sex with men: a nationwide internet survey conducted in Japan. BMC public health. 2006;6:239. Epub 2006/09/28. pmid:17002800; PubMed Central PMCID: PMC1599727.
  21. 21. Greenwood GL, White EW, Page-Shafer K, Bein E, Osmond DH, Paul J, et al. Correlates of heavy substance use among young gay and bisexual men: The San Francisco Young Men's Health Study. Drug and alcohol dependence. 2001;61(2):105–12. Epub 2001/01/04. pmid:11137274.
  22. 22. Chen YM, Yu PS, Lin CC, Jen I. Surveys of HIV-1, HTLV-I, and other sexually transmitted diseases in female sex workers in Taipei City, Taiwan, from 1993 to 1996. Journal of acquired immune deficiency syndromes and human retrovirology: official publication of the International Retrovirology Association. 1998;18(3):299–303. pmid:9665510.
  23. 23. Huang YF, Yang JY, Nelson KE, Kuo HS, Lew-Ting CY, Yang CH, et al. Changes in HIV incidence among people who inject drugs in Taiwan following introduction of a harm reduction program: a study of two cohorts. PLoS medicine. 2014;11(4):e1001625. pmid:24714449; PubMed Central PMCID: PMC3979649.
  24. 24. Wei M, Guan Q, Liang H, Chen J, Chen Z, Hei F, et al. Simple subtyping assay for human immunodeficiency virus type 1 subtypes B, C, CRF01-AE, CRF07-BC, and CRF08-BC. Journal of clinical microbiology. 2004;42(9):4261–7. Epub 2004/09/15. pmid:15365021; PubMed Central PMCID: PMC516361.
  25. 25. Hall TA. BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symposium Series1999. p. 95–8. pmid:10780396
  26. 26. Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: Molecular Evolutionary Genetics Analysis version 6.0. Molecular biology and evolution. 2013;30(12):2725–9. pmid:24132122; PubMed Central PMCID: PMC3840312.
  27. 27. Saitou N, Nei M. The neighbor-joining method: a new method for reconstructing phylogenetic trees. Molecular biology and evolution. 1987;4(4):406–25. pmid:3447015.
  28. 28. Felsenstein J. Using the quantitative genetic threshold model for inferences between and within species. Philosophical transactions of the Royal Society of London Series B, Biological sciences. 2005;360(1459):1427–34. pmid:16048785; PubMed Central PMCID: PMC1569509.
  29. 29. David M. Hillis JJB. An Empirical Test of Bootstrapping as a Method for Assessing Confidence in Phylogenetic Analysis. Systematic Biology. 1993;42(2):11.
  30. 30. Hillis DM, Bull JJ. An Empirical Test of Bootstrapping as a Method for Assessing Confidence in Phylogenetic Analysis. Systematic Biology. 1993;42(2):182–92.
  31. 31. Fleming DT, Wasserheit JN. From epidemiological synergy to public health policy and practice: the contribution of other sexually transmitted diseases to sexual transmission of HIV infection. Sexually transmitted infections. 1999;75(1):3–17. Epub 1999/08/17. pmid:10448335; PubMed Central PMCID: PMC1758168.
  32. 32. Zhong F, Lin P, Xu H, Wang Y, Wang M, He Q, et al. Possible increase in HIV and syphilis prevalence among men who have sex with men in Guangzhou, China: results from a respondent-driven sampling survey. AIDS and behavior. 2011;15(5):1058–66. Epub 2009/10/15. pmid:19826942.
  33. 33. Parekh BS, Kennedy MS, Dobbs T, Pau CP, Byers R, Green T, et al. Quantitative detection of increasing HIV type 1 antibodies after seroconversion: a simple assay for detecting recent HIV infection and estimating incidence. AIDS research and human retroviruses. 2002;18(4):295–307. pmid:11860677.
  34. 34. Duong YT, Qiu M, De AK, Jackson K, Dobbs T, Kim AA, et al. Detection of recent HIV-1 infection using a new limiting-antigen avidity assay: potential for HIV-1 incidence estimates and avidity maturation studies. PloS one. 2012;7(3):e33328. pmid:22479384; PubMed Central PMCID: PMC3314002.
  35. 35. Hauser A, Santos-Hoevener C, Meixenberger K, Zimmermann R, Somogyi S, Fiedler S, et al. Improved testing of recent HIV-1 infections with the BioRad avidity assay compared to the limiting antigen avidity assay and BED Capture enzyme immunoassay: evaluation using reference sample panels from the German Seroconverter Cohort. PloS one. 2014;9(6):e98038. pmid:24892795; PubMed Central PMCID: PMC4043688.
  36. 36. Longosz AF, Serwadda D, Nalugoda F, Kigozi G, Franco V, Gray RH, et al. Impact of HIV subtype on performance of the limiting antigen-avidity enzyme immunoassay, the bio-rad avidity assay, and the BED capture immunoassay in Rakai, Uganda. AIDS research and human retroviruses. 2014;30(4):339–44. pmid:24083837; PubMed Central PMCID: PMC3976571.
  37. 37. Chen YM, Huang KL, Jen I, Chen SC, Liu YC, Chuang YC, et al. Temporal trends and molecular epidemiology of HIV-1 infection in Taiwan from 1988 to 1998. J Acquir Immune Defic Syndr. 2001;26(3):274–82. Epub 2001/03/10. pmid:11242201.
  38. 38. Lan YC, Elbeik T, Dileanis J, Ng V, Chen YJ, Leu HS, et al. Molecular epidemiology of HIV-1 subtypes and drug resistant strains in Taiwan. Journal of medical virology. 2008;80(2):183–91. pmid:18098141.
  39. 39. Chen YJ, Huang YH, Chuang SY, Kao DY, Lan YC, Yang JY, et al. Molecular epidemiology of HIV-1 subtype B, CRF01_AE, and CRF07_BC infection among injection drug users in Taiwan. J Acquir Immune Defic Syndr. 2010;53(4):425–39. pmid:20130472.
  40. 40. Brenner BG, Roger M, Routy JP, Moisi D, Ntemgwa M, Matte C, et al. High rates of forward transmission events after acute/early HIV-1 infection. The Journal of infectious diseases. 2007;195(7):951–9. pmid:17330784.
  41. 41. Hue S, Pillay D, Clewley JP, Pybus OG. Genetic analysis reveals the complex structure of HIV-1 transmission within defined risk groups. Proceedings of the National Academy of Sciences of the United States of America. 2005;102(12):4425–9. pmid:15767575; PubMed Central PMCID: PMC555492.
  42. 42. Beyrer C, Baral SD, van Griensven F, Goodreau SM, Chariyalertsak S, Wirtz AL, et al. Global epidemiology of HIV infection in men who have sex with men. Lancet. 2012;380(9839):367–77. pmid:22819660; PubMed Central PMCID: PMC3805037.
  43. 43. Baral S, Sifakis F, Cleghorn F, Beyrer C. Elevated risk for HIV infection among men who have sex with men in low- and middle-income countries 2000–2006: a systematic review. PLoS medicine. 2007;4(12):e339. pmid:18052602; PubMed Central PMCID: PMC2100144.
  44. 44. Huang SW, Wang SF, Lin YT, Yen CH, Lee CH, Wong WW, et al. Patients Infected with CRF07_BC Have Significantly Lower Viral Loads than Patients with HIV-1 Subtype B: Mechanism and Impact on Disease Progression. PloS one. 2014;9(12):e114441. pmid:25502811
  45. 45. Rutvisuttinunt W, Sirivichayakul S, Oota S, Assawadarachai V, Poltavee K, Savadsuk H, et al. Two unique recombinant forms identified in incident HIV type 1 infections in Thai blood donors. AIDS research and human retroviruses. 2012;28(12):1703–11. pmid:22587412.
  46. 46. Ye J, Xin R, Yu S, Bai L, Wang W, Wu T, et al. Phylogenetic and temporal dynamics of human immunodeficiency virus type 1 CRF01_AE in China. PloS one. 2013;8(1):e54238. pmid:23365653; PubMed Central PMCID: PMC3554705.
  47. 47. Zhao B, Han X, Dai D, Liu J, Ding H, Xu J, et al. New trends of primary drug resistance among HIV type 1-infected men who have sex with men in Liaoning Province, China. AIDS research and human retroviruses. 2011;27(10):1047–53. pmid:21417755.
  48. 48. C. Kittinunvorakoon PW, W. Leelawiwat, P. Poolviwat, O. Kongpechsatit, F. Mueanpai, A. Charoenwatanachokchai, A. Varangrat, S. Chaikummao, P. Akarasewi, F. van Griensven, J. McNicholl. Distribution of HIV-1 strains among men who have sex with men in Bangkok, Thailand. XVII International AIDS Conference: Abstract no THPE00242008.
  49. 49. Beyrer C, Sullivan PS, Sanchez J, Dowdy D, Altman D, Trapence G, et al. A call to action for comprehensive HIV services for men who have sex with men. Lancet. 2012;380(9839):424–38. pmid:22819663; PubMed Central PMCID: PMC3805059.
  50. 50. Zhang X, Wang C, Hengwei W, Li X, Li D, Ruan Y, et al. Risk factors of HIV infection and prevalence of co-infections among men who have sex with men in Beijing, China. AIDS. 2007;21 Suppl 8:S53–7. Epub 2008/01/10. pmid:18172392.
  51. 51. Ellis RJ, Childers ME, Cherner M, Lazzaretto D, Letendre S, Grant I. Increased human immunodeficiency virus loads in active methamphetamine users are explained by reduced effectiveness of antiretroviral therapy. The Journal of infectious diseases. 2003;188(12):1820–6. Epub 2003/12/16. pmid:14673760.
  52. 52. Rusch M, Lampinen TM, Schilder A, Hogg RS. Unprotected anal intercourse associated with recreational drug use among young men who have sex with men depends on partner type and intercourse role. Sexually transmitted diseases. 2004;31(8):492–8. Epub 2004/07/27. pmid:15273583.
  53. 53. Lee SJ, Galanter M, Dermatis H, McDowell D. Circuit parties and patterns of drug use in a subset of gay men. Journal of addictive diseases. 2003;22(4):47–60. Epub 2004/01/16. pmid:14723477.
  54. 54. Green SM, Roback MG, Kennedy RM, Krauss B. Clinical practice guideline for emergency department ketamine dissociative sedation: 2011 update. Annals of emergency medicine. 2011;57(5):449–61. pmid:21256625.
  55. 55. Greene SL, Kerr F, Braitberg G. Review article: amphetamines and related drugs of abuse. Emergency medicine Australasia: EMA. 2008;20(5):391–402. pmid:18973636.
  56. 56. Kelly BC. Conceptions of risk in the lives of club drug-using youth. Substance use & misuse. 2005;40(9–10):1443–59. Epub 2005/07/29. pmid:16048827.
  57. 57. Maxwell JC. Party drugs: properties, prevalence, patterns, and problems. Substance use & misuse. 2005;40(9–10):1203–40. Epub 2005/07/29. pmid:16048814.
  58. 58. Patterson TL, Semple SJ, Zians JK, Strathdee SA. Methamphetamine-using HIV-positive men who have sex with men: correlates of polydrug use. Journal of urban health: bulletin of the New York Academy of Medicine. 2005;82(1 Suppl 1):i120–6. Epub 2005/03/02. pmid:15738313; PubMed Central PMCID: PMC3456168.
  59. 59. Reback CJ, Larkins S, Shoptaw S. Methamphetamine abuse as a barrier to HIV medication adherence among gay and bisexual men. AIDS care. 2003;15(6):775–85. Epub 2003/11/18. pmid:14617499.
  60. 60. Romanelli F, Smith KM, Pomeroy C. Use of club drugs by HIV-seropositive and HIV-seronegative gay and bisexual men. Topics in HIV medicine: a publication of the International AIDS Society, USA. 2003;11(1):25–32. Epub 2003/04/30. pmid:12717048.
  61. 61. Leri F, Bruneau J, Stewart J. Understanding polydrug use: review of heroin and cocaine co-use. Addiction. 2003;98(1):7–22. Epub 2002/12/21. pmid:12492751.
  62. 62. Chesson HW, Pinkerton SD. Sexually transmitted diseases and the increased risk for HIV transmission: implications for cost-effectiveness analyses of sexually transmitted disease prevention interventions. J Acquir Immune Defic Syndr. 2000;24(1):48–56. Epub 2000/07/06. pmid:10877495.
  63. 63. Wasserheit JN. Epidemiological synergy. Interrelationships between human immunodeficiency virus infection and other sexually transmitted diseases. Sexually transmitted diseases. 1992;19(2):61–77. Epub 1992/03/01. pmid:1595015.
  64. 64. Martin DJ. Inappropriate lubricant use with condoms by homosexual men. Public Health Rep. 1992;107(4):468–73. Epub 1992/07/01. pmid:1641445; PubMed Central PMCID: PMC1403679.
  65. 65. van Griensven F, Thienkrua W, McNicholl J, Wimonsate W, Chaikummao S, Chonwattana W, et al. Evidence of an explosive epidemic of HIV infection in a cohort of men who have sex with men in Thailand. AIDS. 2013;27(5):825–32. pmid:23169330.