Fertility of Czech Gay and Straight Men, Women, and Their Relatives: Testing the Sexually Antagonistic Gene Hypothesis.

One proposal for the persistence of homosexuality in the human population is the sexually antagonistic gene hypothesis, which suggests that the lower fertility of homosexual individuals, especially men, may be compensated by higher fertility of their relatives of the opposite sex. To test this hypothesis, we have collected data from 7,312 heterosexual men, 459 gay men, 3,352 heterosexual women, and 79 lesbian women mainly from Czechia. In an online survey, participants answered questions regarding their own as well as their parents' and grandparents' fertility. For men, we obtained no significant results except for higher fertility of gay men's paternal grandmothers, but the magnitude of this effect was very small. For the female sample, we recorded lower fertility of lesbian women's mothers and fathers. In line with our expectations, both gay men and lesbian women had lower fertility rates than their heterosexual counterparts. Our results are consistent with recent studies which likewise do not support the sexually antagonistic gene hypothesis.

The Desirable Dad Hypothesis: Male Same-Sex Attraction as the Product of Selection for Paternal Care via Antagonistic Pleiotropy.

Same-sex attraction, a heritable trait with a reproductive cost, lacks a comprehensive evolutionary explanation. Here we build on a hypothesis invoking antagonistic pleiotropy, which suggests that genes linked to male same-sex attraction remain in the gene pool because they have conferred some fitness advantage to heterosexual men possessing them. We posit the "desirable dad hypothesis," which proposes that alleles linked to male non-heterosexual orientations increase traits conducive to childcare; heterosexual men possessing same-sex attracted alleles are more desirable mating partners as a function of possessing superior paternal qualities. We conducted three studies to test predictions from this hypothesis. Results were consistent with all three predictions. Study 1 (N = 1632) showed that heterosexual men with same-sex attracted relatives were more feminine than men without, as indicated by self-report measures of femininity (η2 = .007), warmth (η2 = .002), and nurturance (η2 = .004 - .006). In Study 2 (N = 152), women rated feminine male profiles as more romantically appealing than masculine ones (d = 0.83)-but less so than profiles possessing a combination of feminine and masculine traits. In Study 3 (N = 153), women perceived feminine male profiles as depicting the best fathers and masculine profiles the worst (d = 1.56): consistent with the idea that femininity is attractive for childcare reasons. Together, these findings are consistent with the idea that sexual selection for male parental care may be involved in the evolution of male same-sex attraction.

Association of Polymorphisms in NHEJ Pathway Genes with HIV-1 Infection and AIDS Progression in a Northern Chinese MSM Population.

Background and Aims: Men who have sex with men (MSM) are at high risk of HIV infection. The nonhomologous end joining (NHEJ) pathway is the main way of double-stranded DNA break (DSB) repair in the higher eukaryotes and can repair the DSB timely at any time in cell cycle. It is also indicated that the NHEJ pathway is associated with HIV-1 infection since the DSB in host genome DNA occurs in the process of HIV-1 integration. The aim of the present investigation was to evaluate associations of single-nucleotide polymorphisms (SNPs) in NHEJ pathway genes with susceptibility to HIV-1 infection and AIDS progression among MSM residing in northern China. Methods: A total of 481 HIV-1 seropositive men and 493 HIV-1 seronegative men were included in this case-control study. Genotyping of 22 SNPs in NHEJ pathway genes was performed using the SNPscan™ Kit. Results: Positive associations were observed between XRCC6 rs132770 and XRCC4 rs1056503 genotypes and the susceptibility to HIV-1 infection. In gene-gene interaction analysis, significant SNP-SNP interactions of XRCC6 and XRCC4 genetic variations were found to play a potential role in the risk of HIV-1 infection. In stratified analysis, XRCC5 rs16855458 was significantly associated with CD4+ T cell counts in AIDS patients, whereas LIG4 rs1805388 was linked to the clinical phases of AIDS patients. Conclusions: NHEJ gene polymorphisms can be considered to be risk factors of HIV-1 infection and AIDS progression in the northern Chinese MSM population.

Prevalence of high-risk human papillomavirus DNA and mRNA and its association with abnormal anal cytology in the Czech male anal cancer screening cohort.

BACKGROUND: Anal cancer (AC) screening is justified in high-risk populations, particularly HIV-positive men having sex with men (MSM). HR-HPV testing could improve the efficiency of cytologically based screening of AC, as in the screening of biologically analogical cervical cancer. The specificity of HR-HPV testing is influenced by the prevalence of HR-HPV infection in the screened population. Reported anal HR-HPV DNA prevalence in MSM is high, but HR-HPV mRNA reflects rather long-term infections and is more specific for high-grade lesions. However, no data were published about HR-HPV DNA and mRNA prevalence in the Czech AC screening population. METHOD: Results of liquid-based anal cytology of 203 predominantly HIV-positive MSM from the Czech AC screening cohort were correlated with results of DNA and E6/E7 mRNA testing of 14 HR-HPV types, and HPV16 genotyping. Eighty-one MSM underwent a standard anoscopy. RESULTS: A total of 109 (53.7%) samples had abnormal cytology, with 12 (5.9%) ASC-H/HSIL, 67 (33.0%) samples cytologically negative, and 27 (13.3%) unsatisfactory. HR-HPV DNA was detected in 134 (66.0%) and HR-HPV RNA in 72 (35.5%) anal smears. HR-HPV mRNA and HPV16 mRNA positivity were associated with abnormal cytology (p = .0037, p = .0021). No significant association was found between HR-HPV DNA or HPV16 DNA positivity and abnormal cytology. No high-grade lesions were revealed by anoscopy. CONCLUSION: Prevalence of anal HR-HPV DNA among Czech MSM is high, however, the prevalence of HR-HPV mRNA is half and associated with abnormal cytology. Our results indicate an increased efficiency of cytological screening when combined with HR-HPV mRNA testing.

Genome-Wide Linkage and Association Study of Childhood Gender Nonconformity in Males.

Male sexual orientation is influenced by environmental and complex genetic factors. Childhood gender nonconformity (CGN) is one of the strongest correlates of homosexuality with substantial familiality. We studied brothers in families with two or more homosexual brothers (409 concordant sibling pairs in 384 families, as well as their heterosexual brothers), who self-recalled their CGN. To map loci for CGN, we conducted a genome-wide linkage scan (GWLS) using SNP genotypes. The strongest linkage peaks, each with significant or suggestive two-point LOD scores and multipoint LOD score support, were on chromosomes 5q31 (maximum two-point LOD = 4.45), 6q12 (maximum two-point LOD = 3.64), 7q33 (maximum two-point LOD = 3.09), and 8q24 (maximum two-point LOD = 3.67), with the latter not overlapping with previously reported strongest linkage region for male sexual orientation on pericentromeric chromosome 8. Family-based association analyses were used to identify associated variants in the linkage regions, with a cluster of SNPs (minimum association p = 1.3 × 10-8) found at the 5q31 linkage peak. Genome-wide, clusters of multiple SNPs in the 10-6 to 10-8 p-value range were found at chromosomes 5p13, 5q31, 7q32, 8p22, and 10q23, highlighting glutamate-related genes. This is the first reported GWLS and genome-wide association study on CGN. Further increasing genetic knowledge about CGN and its relationships to male sexual orientation should help advance our understanding of the biology of these associated traits.

The Identification of A Novel HIV-1 Second-Generation Recombinant form (CRF01_AE/CRF07_BC) Among Men Who Have Sex with Men in Jiangsu, China.

BACKGROUND: CRF01_AE and CRF07_BC are the two major HIV-1 virus strains circulating in China. The proportion of dominant subtypes (CRF01_AE and CRF07_BC) among MSM in Jiangsu province was over 80%. A large number of URFs have been found in China in recent years. OBJECTIVE: This study aimed to report on novel HIV-1 recombinants. METHODS: We constructed Phylogenetic trees using the maximum likelihood (ML) method with 1000 bootstrap replicates in IQ-TREE 1.6.8 software and determined recombination breakpoints using SimPlot 3.5.1. RESULTS: We identified a novel, second-generation HIV-1 recombinant (JS020202) between CRF01_ AE and CRF07_BC. The analysis of near full-length genome (NFLG) showed there were at least 8 breakpoints in the virus, which differed from any previously identified CRF and URF around the world. CONCLUSION: Novel diverse CRF01_AE/07_BC suggested the complexity trends of HIV-1 genetics. The emergency situation of diverse recombinant strains should be monitored continuously.

Genomic analysis of the meningococcal ST-4821 complex-Western clade, potential sexual transmission and predicted antibiotic susceptibility and vaccine coverage.

INTRODUCTION: The ST-4821 complex (cc4821) is a leading cause of serogroup C and serogroup B invasive meningococcal disease in China where diverse strains in two phylogenetic groups (groups 1 and 2) have acquired fluoroquinolone resistance. cc4821 was recently prevalent among carriage isolates in men who have sex with men in New York City (USA). Genome-level population studies have thus far been limited to Chinese isolates. The aim of the present study was to build upon these with an extended panel of international cc4821 isolates. METHODS: Genomes of isolates from Asia (1972 to 2017), Europe (2011 to 2018), North America (2007), and South America (2014) were sequenced or obtained from the PubMLST Neisseria database. Core genome comparisons were performed in PubMLST. RESULTS: Four lineages were identified. Western isolates formed a distinct, mainly serogroup B sublineage with alleles associated with fluoroquinolone susceptibility (MIC <0.03 mg/L) and reduced penicillin susceptibility (MIC 0.094 to 1 mg/L). A third of these were from anogenital sites in men who have sex with men and had unique denitrification gene alleles. Generally 4CMenB vaccine strain coverage was reliant on strain-specific NHBA peptides. DISCUSSION: The previously identified cc4821 group 2 was resolved into three separate lineages. Clustering of western isolates was surprising given the overall diversity of cc4821. Possible association of this cluster with the anogenital niche is worthy of monitoring given concerns surrounding antibiotic resistance and potential subcapsular vaccine escape.

Non-androgenic testicular mediation of androphilia in male mice with global overexpression of androgen receptors.

Sexual attraction is robustly sexually differentiated among mammalian species. Gonadal androgens acting perinatally and in adulthood are required for male-typical preference for female sexual cues. Recent evidence suggests that at the high extent of AR signaling, male mice show an increased preference for same-sex odor cues. These findings were found only in mice that overexpress AR globally in all tissues (CMV-AR), whereas neural AR overexpression (Nestin-AR) did not affect sexual preference. The present studies investigated the endocrine basis of this phenotype and examined whether preference for male or female stimulus animals (partner preference) was also affected in these transgenic animals. We manipulated adult gonadal hormones in male mice that overexpress AR globally and males that overexpress AR only in neural tissue. We replicate the finding that androphilia is increased in gonadally intact CMV-AR males, and these males exhibited reduced neural activation in response to estrus female odors. Testosterone treatment of gonadectomized CMV-AR males was sufficient to induce a gynephilic olfactory preference, while a gynephilic partner preference was induced with gonadectomy alone. These findings suggest that altered sexual preference of CMV-AR male mice is mediated by inhibitory activational functions of the testes. Together, these results suggest that at the high extent of AR signaling, non-neural AR via the gonads, can promote androphilia.

Human Immunodeficiency Virus Type 1 Phylodynamics to Detect and Characterize Active Transmission Clusters in North Carolina.

BACKGROUND: Human immunodeficiency virus type 1 (HIV-1) phylodynamics can be used to monitor epidemic trends and help target prevention through identification and characterization of transmission clusters. METHODS: We analyzed HIV-1 pol sequences sampled in North Carolina from 1997 to 2014. Putative clusters were identified using maximum-likelihood trees and dated using Bayesian Markov Chain Monte Carlo inference. Active clusters were defined as clusters including internal nodes from 2009 to 2014. Effective reproductive numbers (Re) were estimated using birth-death models for large clusters that expanded ≥2-fold from 2009 to 2014. RESULTS: Of 14 921 persons, 7508 (50%) sequences were identified in 2264 clusters. Only 288 (13%) clusters were active from 2009 to 2014; 37 were large (10-36 members). Compared to smaller clusters, large clusters were increasingly populated by men and younger persons; however, nearly 60% included ≥1 women. Clusters with ≥3 members demonstrated assortative mixing by sex, age, and sample region. Of 15 large clusters with ≥2-fold expansion, nearly all had Re approximately 1 by 2014. CONCLUSIONS: Phylodynamics revealed transmission cluster expansion in this densely sampled region and allowed estimates of Re to monitor active clusters, showing the propensity for steady, onward propagation. Associations with clustering and cluster characteristics vary by cluster size. Harnessing sequence-derived epidemiologic parameters within routine surveillance could allow refined monitoring of local subepidemics.

Experienced homophobia and gene expression alterations in Black and Latino men who have sex with men in Los Angeles County.

Men who have sex with men (MSM) experience high rates of homophobic victimization, which is linked to myriad chronic physical and mental health disparities. Social adversity such as rejection, isolation, and racial discrimination can induce a conserved transcriptional response to adversity (CTRA) involving upregulation of proinflammatory genes and downregulation of type I interferon and antibody synthesis genes. This study specifically examines whether homophobic victimization is associated with expression of CTRA profiles in Black and Latino MSM living in Los Angeles. Analyses linked behavioral survey data with quantified RNA from leukocytes from blood samples of 70 participants over 12 months. CTRA gene expression was increased by 3.1-fold in MSM who experienced homophobic victimization while adjusting for major leukocyte subsets and sociodemographics. Accounting for all these factors, CTRA gene expression was significantly enhanced in MSM who identified as Black compared to Latino. Our findings identify experiences of homophobic victimization as drivers of inflammatory and type I interferon gene expression profiles, which can contribute to physical and mental health challenges in Black and Latino MSM.

On the Mathematics of the Fraternal Birth Order Effect and the Genetics of Homosexuality.

Mathematicians have always been attracted to the field of genetics. The mathematical aspects of research on homosexuality are especially interesting. Certain studies show that male homosexuality may have a genetic component that is correlated with female fertility. Other studies show the existence of the fraternal birth order effect, that is, the correlation of homosexuality with the number of older brothers. This article is devoted to the mathematical aspects of how these two phenomena are interconnected. In particular, we show that the fraternal birth order effect implies a correlation between homosexuality and maternal fecundity. Vice versa, we show that the correlation between homosexuality and female fecundity implies the increase in the probability of the younger brothers being homosexual.

Large-scale GWAS reveals insights into the genetic architecture of same-sex sexual behavior.

Twin and family studies have shown that same-sex sexual behavior is partly genetically influenced, but previous searches for specific genes involved have been underpowered. We performed a genome-wide association study (GWAS) on 477,522 individuals, revealing five loci significantly associated with same-sex sexual behavior. In aggregate, all tested genetic variants accounted for 8 to 25% of variation in same-sex sexual behavior, only partially overlapped between males and females, and do not allow meaningful prediction of an individual's sexual behavior. Comparing these GWAS results with those for the proportion of same-sex to total number of sexual partners among nonheterosexuals suggests that there is no single continuum from opposite-sex to same-sex sexual behavior. Overall, our findings provide insights into the genetics underlying same-sex sexual behavior and underscore the complexity of sexuality.

Biogenetic ties and parent-child relationships: The misplaced critique.

According to an almost axiomatic standard in bioethics, moral commitment should ground parents' relationship with their children, rather than biogenetic relatedness. This standard has been used lately to express skepticism about extending existing assisted reproductive treatments (ARTs) to same-sex couples and to research into novel fertility interventions for those couples, but this skepticism is misplaced on several grounds. As a matter of access and equity, same-sex couples seem presumptively entitled to genetic relatedness to their children as far as possible both in regard to existing ARTs and to novel ARTs under investigation. For those worried about the effects of trying to secure biogenetic relatedness for same-sex couples, it may be noted that same-sex couples will only ever be a fraction of the parents implicated in propping up "biologism," as the expectation of biogenetic relatedness it is sometimes called. The cultural force of biologism would survive almost intact even if no same-sex couples were ever to have genetically related children. It is therefore hard to see why same-sex couples should forfeit aspirations to biogenetic relationships with their children or enjoy less subsidy for ARTs than the subsidy given to different-sex couples. As matter of moral consistency, the full implications of the biologism critique have yet to be evaluated relative to different-sex couples.

Rapid and complicated HIV genotype expansion among high-risk groups in Guangdong Province, China.

BACKGROUND: Guangdong Province is one of the most developed and populous provinces in southern China, with frequent foreign exchanges and large transient population. The annual number of cases of HIV/AIDS reported in Guangdong has been higher than most of provinces in China for several successive years. HIV infection by heterosexual transmission occurs across the province, with transmission among men who have sex with men occurring mainly in larger urban centers. There is a lack of widespread and representative data on the distribution of HIV subtypes in Guangdong. This study aimed to thoroughly investigate and estimate the prevalence and distribution of HIV-1 subtypes using a city-based sampling strategy to better understand the characteristics of HIV transmission in Guangdong. METHODS: Archived plasma samples (n = 1205) from individuals diagnosed as HIV-1 infection in 2013 were selected randomly from all 21 cities in Guangdong Province. Genotypes were determined using env and/or gag sequences using phylogenetic analysis. The distributions of HIV genotypes in different risk groups and different cities were analyzed. RESULTS: A total of 15 genotypes, including six discordant genotypes, were identified. The four main HIV-1 subtypes in Guangdong were CRF01_AE (43.2%), CRF07_BC (26.3%), CRF55_01B (8.5%), and CRF08_BC (8.4%). CRF01_AE was the predominant subtype in all risk populations. The high mobility of people shaped the complexity of the HIV genotypes, while the switch of risk factors affected the distribution and future trend of HIV-1 genotypes in Guangdong. Another epicenter located in the western region in addition to the known epicenter cities in the Pearl River Delta region of Guangdong may exist. CONCLUSIONS: Our study provides a comprehensive molecular epidemiologic dataset to understand the diversity and distribution of HIV genotypes in Guangdong, as well as to clarify the unique region- and risk group-specific transmission dynamics. The results provide critical and insightful information for more effective intervention strategies to limit HIV transmission in the future.

Increasing proportions of HIV-1 non-B subtypes and of NNRTI resistance between 2013 and 2016 in Germany: Results from the national molecular surveillance of new HIV-diagnoses.

BACKGROUND: Molecular surveillance of newly diagnosed HIV-infections is important for tracking trends in circulating HIV-variants, including those with transmitted drug resistances (TDR) to sustain ART efficacy. METHODS: Dried serum spots (DSS) are received together with the statutory notification of a new diagnosis. 'Recent infections' (<155 days) classified by a 'recent infection test algorithm' (BED-CEIA and clinical data) are genotyped in HIV-protease (PR), reverse transcriptase (RT) and integrase (INT) to determine the HIV-1 subtype, to calculate prevalence and trends of TDR, to predict baseline susceptibility and to identify potential transmission clusters for resistant variants. RESULTS: Between January 2013 and December 2016, 1,885 recent infections were analysed regarding the PR/RT genomic region, with 43.5% of these also being subjected to the analysis of INT. The proportion of HIV-1 non-B viruses (31.3%; 591/1,885) increased from 21.6% to 36.0%, particularly the subtypes A (5.0% to 8.3%) and C (3.2% to 7.7%; all ptrends < 0.01). The subtype A increment is mainly due to transmissions within men who have sex with men (MSM) while subtype C transmissions are associated with heterosexuals and people who inject drugs. The prevalence of TDR was stable at 11.0% (208/1,885) over the study period. Resistances to nucleotide RT inhibitors (NRTI) and PR inhibitors (PI) were 4.5% and 3.2%, respectively, without identifiable trends. In contrast, resistances to non-NRTIs (NNRTI, 4.7%) doubled between 2014 and 2016 from 3.2% to 6.4% (ptrend = 0.02) mainly due to the K103N mutation (from 1.7% to 4.1%; ptrend = 0.03) predominantly detected in recently infected German MSM not linked to transmission clusters. Transmitted INSTI mutations were present in only one case (T66I) and resistance to dolutegravir was not identified at all. Reduced susceptibility to recommended first-line therapies was low with 1.0% for PIs, 1.3% for NRTIs and 0.7% for INSTIs, but high for the NNRTIs efavirence (4.9%) and rilpivirine (6.0%) due to the K103N mutation and the polymorphic mutation E138A. These trends in therapy-naïve individuals impact current first-line regimens and require awareness and vigilant surveillance.

Pre-exposure prophylaxis (PrEP) for MSM in low HIV incidence places: should high risk individuals be targeted?

Pre-exposure prophylaxis (PrEP) targeting high-risk men who have sex with men (MSM) has been shown to be a cost-effective HIV control measure. However, the approach could be a challenge in low HIV incidence places with a low proportion of high-risk MSM. To examine the impact of PrEP in such setting in Asia, we developed an epidemic model and conducted cost-effectiveness analysis using empirical multicentre clinical and HIV sequence data from HIV-infected MSM in Hong Kong, in conjunction with behavioural data of local MSM. Without PrEP, the HIV incidence (per 100 person-years) would increase from 1.1 to 1.6 between 2011 and 2021. PrEP could avert 3-63% of total new infections in a five-year period (2017-2021), the variability of which depends on the implementation strategies and combination with test-and-treat. However, under current market drug price in 2016, the incremental cost per quality-adjusted life-year gained (QALYG) of PrEP (USD1583136/QALYG) is almost 3 times higher than test-and-treat intervention alone (USD396874/QALYG). Assuming 93% fall of PrEP drug price and in combination with test-and-treat, putting 30% of MSM on non-targeting PrEP would be more feasible, cost-effective (USD268915/QALYG), and could avert more new infections (40%). PrEP could contribute to HIV epidemic control in a low incidence place.

First report of sexually transmitted multi-drug resistant Shigella sonnei infections in Switzerland, investigated by whole genome sequencing.

Shigella sonnei causes foodborne infections, but has recently also been described as a sexually transmitted infection (STI), with increased levels of antimicrobial resistance. We describe three cases of sexually acquired Shigella sonnei infection - the first report of this emerging infection in Switzerland. We used in-house whole genome sequencing (WGS) to investigate possible transmission routes and epidemiological correlations between the three cases. The genomic analysis demonstrated that two of three case isolates were very closely related, with only two single nucleotide polymorphism differences between them, despite being isolated from two unrelated patients at time-points six months apart, and the infections having been acquired at different geographic locations within Europe. All three isolates were found to fall within two of the clusters (1 and 7) defined within UK men who have sex with men (MSM) isolate populations, but with higher divergence, suggesting a more diverse pool circulating within Europe. Phenotypic testing confirmed the genotypic findings, with all three isolates azithromycin resistant, and two out of three resistant to quinolones. This report underlines the importance of the implementation of new sequencing technologies in the investigation of epidemiological aspects of this STI circulating in the population of MSM. In such cases, therapy should always be guided by antimicrobial susceptibility testing owing to increasing resistances. Greater awareness of this emerging sexually transmitted infection is needed.

What Is the most Important for Elite Control: Genetic Background of Patient, Genetic Background of Partner, both or neither? Description of Complete Natural History within a Couple of MSM.

BACKGROUND: We describe a homosexual man who strongly controlled HIV-1 for ten years despite lack of protective genetic background. METHODS: HIV-1 DNA was measured in blood and other tissues. Cell susceptibility was evaluated with various strains. HIV-1-specific (CD4 and CD8 activation markers and immune check points) and NK cells responses were assessed; KIRs haplotypes and HLA alleles were determined. FINDINGS: Two HIV-1 RNA copies/mL of plasma were detected in 2009, using an ultra-sensitive assay. HIV-DNA was detected at 1.1 and 2 copies/106 PBMCs in 2009 and 2015 respectively, at 1.2 copies/106 cells in rectal cells in 2011. WBs showed weak reactivity with antibodies to gp160, p55 and p25 from 2007 to 2014, remaining incomplete in 2017. CD4 T cells were susceptible to various strains including HIVKON, a primary isolate of his own CRF02_AG variant. CD8 T cells showed a strong poly-functional response against HIV-Gag, producing mainly IFN-γ; a robust capacity of antibody-dependant cell cytotoxicity (ADCC) was observed in NK cells. Case patient was group B KIR haplotype. Neutralizing antibodies were not detected. CD4 and CD8 blood T cells showed normal proportions without increased activation markers. Phylogenetic analyses identified the same CRF02_AG variant in his partner. The patient and his partner were heterozygous for the CCR5ΔD32 deletion and shared HLA-B*07, C*07 non-protective alleles. INTERPRETATION: This thorough description of the natural history of an individual controlling HIV-1 in various compartments for ten years despite lack of protective alleles, and of his partner, may have implications for strategies to cure HIV-1 infection.

The second molecular epidemiological study of HIV infection in Mongolia between 2010 and 2016.

OBJECTIVE: Our previous 2005-2009 molecular epidemiological study in Mongolia identified a hot spot of HIV-1 transmission in men who have sex with men (MSM). To control the infection, we collaborated with NGOs to promote safer sex and HIV testing since mid-2010. In this study, we carried out the second molecular epidemiological survey between 2010 and 2016 to determine the status of HIV-1 infection in Mongolia. METHODS: The study included 143 new cases of HIV-1 infection. Viral RNA was extracted from stocked plasma samples and sequenced for the pol and the env regions using the Sanger method. Near-full length sequencing using MiSeq was performed in 3 patients who were suspected to be infected with recombinant HIV-1. Phylogenetic analysis was performed using the neighbor-joining method and Bayesian Markov chain Monte Carlo method. RESULTS: MSM was the main transmission route in the previous and current studies. However, heterosexual route showed a significant increase in recent years. Phylogenetic analysis documented three taxa; Mongolian B, Korean B, and CRF51_01B, though the former two were also observed in the previous study. CRF51_01B, which originated from Singapore and Malaysia, was confirmed by near-full length sequencing. Although these strains were mainly detected in MSM, they were also found in increasing numbers of heterosexual males and females. Bayesian phylogenetic analysis estimated transmission of CRF51_01B into Mongolia around early 2000s. An extended Bayesian skyline plot showed a rapid increase in the effective population size of Mongolian B cluster around 2004 and that of CRF51_01B cluster around 2011. CONCLUSIONS: HIV-1 infection might expand to the general population in Mongolia. Our study documented a new cluster of HIV-1 transmission, enhancing our understanding of the epidemiological status of HIV-1 in Mongolia.