Projected HIV and Bacterial STI Incidence Following COVID-Related Sexual Distancing and Clinical Service Interruption
=====================================================================================================================

* Samuel M. Jenness
* Adrien Le Guillou
* Christina Chandra
* Laura M. Mann
* Travis Sanchez
* Daniel Westreich
* Julia L. Marcus

## ABSTRACT

**Background** The global COVID-19 pandemic has the potential to indirectly impact the transmission dynamics and prevention of HIV and other sexually transmitted infections (STI). Studies have already documented reductions in sexual activity (“sexual distancing”) and interruptions in HIV/STI services, but it is unknown what combined impact these two forces will have on HIV/STI epidemic trajectories.

**Methods** We adapted a network-based model of co-circulating HIV, gonorrhea, and chlamydia for a population of approximately 103,000 men who have sex with men (MSM) in the Atlanta area. Model scenarios varied the timing, overlap, and relative extent of COVID-related sexual distancing in casual and one-time partnership networks and service interruption within four service categories (HIV screening, HIV PrEP, HIV ART, and STI treatment).

**Results** A 50% relative decrease in sexual partnerships and interruption of all clinical services, both lasting 18 months, would generally offset each other for HIV (total 5-year population impact for Atlanta MSM: −227 cases), but have net protective effect for STIs (−23,800 cases). Greater relative reductions and longer durations of service interruption would increase HIV and STI incidence, while greater relative reductions and longer durations of sexual distancing would decrease incidence of both. If distancing lasted only 3 months but service interruption lasted 18 months, the total 5-year population impact would be an additional 890 HIV cases and 57,500 STI cases.

**Conclusions** The counterbalancing impact of sexual distancing and clinical service interruption depends on the infection and the extent and durability of these COVID-related changes. If sexual behavior rebounds while service interruption persists, we project an excess of hundreds of HIV cases and thousands of STI cases just among Atlanta MSM over the next 5 years. Immediate action to limit the impact of service interruptions is needed to address the indirect effects of the global COVID pandemic on the HIV/STI epidemic.

KEYWORDS
*   Mathematical model
*   Sexual networks
*   HIV
*   STI
*   COVID-19
*   Men who have sex with men

## INTRODUCTION

The 2019 novel coronavirus (COVID-19) global pandemic has directly resulted in substantial morbidity and mortality, but has also indirectly impacted the transmission of other infectious diseases.1 For HIV and other sexually transmitted infections (STIs), behavioral responses to COVID-19 have included major reductions in social contacts that also entailed reductions in sexual activity (“sexual distancing”).2–4 The pandemic also has interrupted the provision of clinical services for HIV/STIs.5 One critical question is how these two phenomena — distancing that could decrease HIV/STI transmission, and service interruption that could increase transmission — will impact the overall incidence of HIV and STIs immediately and in the post-COVID era.

In the United States, men who have sex with men (MSM) are a key population for HIV/STI prevention.6 COVID-19 has already resulted in sexual distancing for MSM.4,7–9 Behavioral changes that began during COVID-related restrictions (March 2020) have included reductions in the number of sexual partners (40–60% of MSM) and lower frequency of sexual activity within partnerships (20%). While some social activities have already rebounded in some areas,10 resumption back to “normal” levels may be delayed until the widespread distribution of a SARS-CoV-2 vaccine. Use of clinical HIV/STI services has also declined among MSM in the U.S. Categories of reduced services include HIV/STI diagnostic screening (25–85%), use of HIV preexposure prophylaxis (PrEP) (20–72%), and retention in HIV care (25%–45%).7,9 Providers have partially addressed these interruptions by replacing in-person clinical visits with telehealth services,11 but these tools may be less available to people with the greatest need.12 Local health departments have also redeployed STI services towards COVID-19 contact tracing efforts.13 It is unknown what the immediate and longer-term impact of these service disruptions will be for uniquely for HIV versus STI incidence.

Parallel reductions in sexual behavior and HIV/STI screening rates make it challenging to understand whether any declines in reported case data reflect true changes in disease incidence or gaps in disease surveillance. Transmission modeling can disentangle these two explanations. In this study, we used a stochastic network-based transmission model to project the impact of sexual distancing versus HIV/STI service interruptions driven by COVID-19. We evaluated how changes to the sexual partnership networks among MSM in Atlanta may reduce the disease incidence in a model that represents the cocirculation of HIV, gonorrhea, and chlamydia. We also explored how COVID-related service disruptions in four categories (HIV screening, HIV PrEP, HIV treatment, and STI treatment) could increase the incidence of these infections.

## METHODS

### Study Design

This model of HIV/STI transmission dynamics for U.S. MSM was built with the *EpiModel* platform,14 which simulates epidemics over dynamic contact networks using temporal exponential random graph models (TERGMs).15 This builds on our previous applied HIV/STI modeling of co-circulating HIV, *Neisseria gonorrhoeae* (NG), and *Chlamydia trachomatis* (CT) among MSM.16 For this study, we implemented time points for the start and end of COVID-related sexual distancing and clinical service interruptions. Full methodological details are provided in the **Supplemental Appendix**.

Our model represented main, casual, and one-time sexual partnerships for Black, Hispanic, and white/other MSM, aged 15 to 65, in Atlanta. The starting network size in the model simulations was 10,000 MSM, which could stochastically increase or decrease over time based on arrival (sexual debut) and departure (mortality or sexual cessation).

### HIV/STI Epidemic Model

The epidemic model consisted of 5 main components: 1) statistical network models (TERGMs) to generate dynamic sexual partnerships; 2) statistical models to predict behavior within partnerships; 3) simulation of pathogen transmission across active partnerships; 4) simulation of disease progression and other natural history features; and 5) simulation of prevention and treatment service engagement.

To fit the network models, we used data from ARTnet, a web-based egocentric network study conducted in 2017–2019 of MSM in the US.17 Parameters were weighted by census-based race/ethnicity and age distributions to account for ARTnet sampling biases. Multivariate predictors of partnership formation included partnership type, degree distributions by partnership types, heterogeneity in network degree, assortative mixing by race/ethnicity and age, and mixing by sexual position. Generalized linear m models were then fit to ARTnet to predict the frequency of acts and the probability of condom use as a multivariate function of race/ethnicity, age, diagnosed HIV status, and partnership type and duration.

MSM could be screened for HIV and initiate anti-retroviral therapy (ART), which would lower their HIV viral load (VL) and increase their longevity. MSM progressed through HIV disease with VLs represented continuously. Lower VL with sustained ART use was associated with a reduced probability of HIV transmission per act. Other factors modifying the HIV transmission probability per act included PrEP use, condom use, sexual position, circumcision, and a prevalent STI.

For HIV services, we represented an integrated HIV continuum of antiretroviral-based prevention and care, with HIV screening as the gateway to both.18 MSM engaged in HIV screening at regular intervals, calibrated to local surveillance data on the proportion of MSM with HIV who were diagnosed.19 MSM screening HIV-positive could then enter the HIV care continuum (linkage and retention in ART) while MSM who screened negative could enter the HIV prevention continuum (PrEP initiation, adherence, and persistence). MSM were linked to ART, and could cycle off and back on ART based on rates calibrated to local surveillance of care entry and VL suppression.19 MSM on ART achieved suppression after 3 months, with a rebound to set-point VL after halting ART.

The HIV prevention continuum consisted of initiation, adherence, and persistence in PrEP care for daily oral tenofovir/emtricitabine.20 MSM who tested HIV-negative and met indications for PrEP based on CDC guidelines were eligible to start.21 Eligible MSM then started PrEP based on an initiation probability generating a coverage level of 15%, consistent with Atlanta estimates.22 Heterogeneous PrEP adherence was modeled, with 78.4% meeting a high-adherence level that resulted in a 99% relative reduction in HIV acquisition risk.23 PrEP discontinuation was based on secondary estimates of the proportion of MSM who were retained in PrEP care at 6 months (57%).24 PrEP care consisted of routine HIV and STI screening.

Gonorrhea and chlamydia transmission were simulated along the same partnership network as HIV, but with disease recovery through either natural clearance or antibiotic treatment.25 STI transmission was directional and site-specific during anal intercourse at the rectal and urogenital sites. Men could be infected at both anatomical sites and with both gonorrhea and chlamydia. The symptomatic status of the newly acquired infections depended on site of infection, with most rectal infections asymptomatic and most urethral infections symptomatic.26 STI symptoms influenced the probability of testing and treatment, which reduced mean time to clearance. Successful treatment for an STI at one site resulted in clearance at the other site for dual-site infections.

### COVID-Related Impact on Behavior and Services

Experimental scenarios implemented reductions to the partnership network structure and HIV/STI service utilization individually and jointly. All scenarios were simulated for a period of 5 years in weekly time steps. A base scenario kept all parameters constant over that simulation window. Experimental scenarios first simulated one year of no change, followed by the initiation of either sexual distancing only, service interruption only, or both combined. In primary scenarios, service interruption lasted for 18 months and sexual distancing either for 18 months or 3 months. The 18-month window was selected based on predictions of the timeline for COVID-19 vaccine deployment; the 3-month distancing window was selected based on empirical data suggesting a more rapid behavioral rebound of sexual partnerships as early as June 2020.9,10 Sensitivity analyses varied the duration of both.

Distancing was modeled to varying levels by reducing network degree for casual partnerships and incidence for one-time partners; degree for main partnerships remained unchanged given minimal expected impact on cohabitating partners. Service interruption was reflected in four types of HIV/STI interventions: HIV screening, HIV PrEP (reduction in new users and increase in discontinuation for current users), HIV ART (through retention in care), and linked bacterial STI screening and treatment. These were reduced on a relative scale individually and jointly. We selected the joint 50% relative reduction in sexual behavior and clinical services as a key scenario to highlight in figures based on empirical data suggesting upwards of this level of change.7,9

### Calibration, Simulation, and Analysis

Given uncertainty in some model parameters, we calibrated the model with a Bayesian approach that defined prior distributions for these parameters and fit the model to empirical surveillance-based estimates of diagnosed HIV, NG, and CT incidence for the target population. This involved projecting incidence estimates to 2019 based on historical data and rising trends in cases over the past decade. After calibration, for each study scenario, we simulated the model 500 times and summarized the distribution of results with medians and 95% simulation intervals.

Three primary outcomes were: 1) standardized HIV and STI incidence per 100 person-years at risk (PYAR) at 2.5 years (or 18 months after start of COVID-related response); 2) standardized cumulative incidence over 5 years per 1000 disease-susceptible MSM; and 3) the total 5-year population impact, which was calculated in two steps. We first multiplied the standardized cumulative incidence by estimates of the total susceptible population size of MSM in the Atlanta metropolitan area (102,642 sexually active MSM for STI outcomes, and 87,723 sexually active HIV-negative MSM for HIV outcomes27) to quantify the total population 5-year incidence. We then subtracted this total population incidence for each scenario from the value in the base scenario to obtain an absolute difference.

## RESULTS

**Figure 1** visualizes the primary scenarios of 18 months of service interruption and 18 months of sexual distancing. The panels show standardized incidence rates of HIV and STIs for scenarios in which sexual behavior and services were jointly reduced by 50% during the COVID period. **Panel A** shows that distancing in the absence of service reduction (green line) was associated with a decrease in HIV incidence, whereas service reduction only (red line) was associated with an increase in HIV incidence. HIV incidence changes in both scenarios persisted after the resumption of baseline behavior and services at year 2.5 (week 130). In the combined scenario (blue line), these relative changes in behavior and services effectively largely each other, resulting in minimal difference in HIV incidence compared with the base (no change) scenario. **Panel B** shows the impact of the same scenarios on combined gonorrhea and chlamydia incidence. Sexual distancing had a strong and sustained reduction in STI incidence. With only service interruption but no distancing, STI incidence increased substantially. In the combined scenario, the projected incidence was slightly lower initially, rebound after sexual distancing ends, and then continued to decline again through year 5 (week 260).

![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/10/21/2020.09.30.20204529/F1.medium.gif)

[Figure 1.](http://medrxiv.org/content/early/2020/10/21/2020.09.30.20204529/F1)

Figure 1. 
HIV and combined STI incidence before, during, and after an 18-month period of clinical service interruption and 18-month period of sexual distancing. Model scenarios within each panel compare a 50% relative reduction in behavior only, service interruption only, or both jointly against the base (no change) model. Thick lines show median values and bands show inter-quartile range of values across 500 simulations per scenario.

**Table 1** quantifies the estimated impact of the sexual-distancing-only scenarios with an 18-month duration. A relative 50% reduction in casual and one-time networks jointly (corresponding to the green line scenario in Figure 1) resulted in a decrease in HIV incidence from 1.23 per 100 PYAR to 0.79, a 36% relative reduction. The 5-year cumulative incidence per 1,000 susceptible MSM was 51.4 cases in that scenario, compared to 62.3 cases in the base scenario, corresponding with a 17% reduction in HIV incidence. The impact on the point incidence at year 2.5 was more extreme than the impact on the cumulative incidence across 5 years because the point incidence rebounded to the baseline level after distancing and service interruptions ended. The total population impact on HIV incidence for Atlanta MSM was projected in that scenario to be −966 cases, or 193 fewer cases per year, compared to the base scenario. The relative size of reductions in casual partnerships correlated with stronger relative declines in HIV incidence than comparable relative reductions in one-time partnerships. This was driven by the underlying behavioral patterns, in which coital frequency was higher and condom use lower for casual partnerships. The mechanistic impact on network degree statistics is provided in **Supplemental Table 14**.

View this table:
[Table 1.](http://medrxiv.org/content/early/2020/10/21/2020.09.30.20204529/T1)

Table 1. 
Changes in HIV and Combined Gonorrhea & Chlamydia Incidence Following an 18-Month Period of Sexual Distancing, with no Associated Changes in HIV/STI Prevention or Treatment Service Utilization

The patterns were generally similar for STI incidence across the same **Table 1** scenarios of sexual distancing only with an 18-month duration. The baseline STI incidence rate was 19.39 (per 100 PYAR), but 50% joint network distancing resulted in an estimated incidence of 5.85 at year 2.5 (a 70% reduction) and a cumulative incidence of 446 cases compared to 960 in the base scenario (a 54% reduction). The 5-year population impact of this scenario was a reduction in total STI cases by over 50,000. Also compared to HIV, 50% reductions within each partnership sub-network were associated with similar declines in STI incidence, whereas HIV was more strongly impacted by casual network reductions. Individual STI outcomes are provided in **Supplemental Table 15**.

**Table 2** quantifies the impact of HIV/STI service reductions in the absence of sexual distancing. The scenario in which all services were reduced by 50% corresponds to the red lines in the **Figure 1**. Reductions in PrEP had a more moderate impact on HIV incidence than reductions in ART. This was because baseline PrEP coverage was relatively low (15%) in Atlanta, and COVID-related changes only applied to the PrEP-indicated population. In comparison, changes in ART retention would impact the care of the entire HIV-diagnosed population, subsequently affecting transmission through lower levels of VL suppression. For STIs, reductions in STI treatment had a dramatic impact on STI incidence, which was held in equilibrium in the base scenario by high levels of routine screening. Gaps in STI screening, which were projected to dramatically increase STI incidence, also had downstream effects on HIV incidence through the biological relationship between prevalent STI infection and HIV acquisition risk. Process outcomes associated for these scenarios are provided in **Supplemental Table 16** and individual STI outcomes in **Supplemental Table 17**.

View this table:
[Table 2.](http://medrxiv.org/content/early/2020/10/21/2020.09.30.20204529/T2)

Table 2. 
Changes in HIV and Combined Gonorrhea & Chlamydia Incidence Following an 18-Month Period of Reduced HIV/STI Prevention or Treatment Service Utilization, with no Associated Changes in Behavior

The top half of **Table 3** shows the impact of combined sexual distancing and service reduction for 18 months. The 50% network reduction and 50% service reduction scenario corresponds to the blue lines in **Figure 1**. For HIV, commensurate relative reductions in behavior and services by the same amounts generally kept incidence similar to the base scenario. The population impact of paired scaled-down scenarios had a minor protective effect at lower reductions (25%/25%: −69; 50%/50%: −227) but resulted in slightly higher incidence for the most extreme reductions (90%/90%: 215). For STIs, in contrast, declines in behavior strongly overwhelmed the declines in services, with a net reduction of 23,800 cases in the 50%/50% paired scaled-down scenario.

View this table:
[Table 3.](http://medrxiv.org/content/early/2020/10/21/2020.09.30.20204529/T3)

Table 3. 
Changes in HIV and Combined Gonorrhea & Chlamydia Incidence Following Joint Behavioral Change and Reduced HIV/STI Prevention and Treatment Service Utilization

**Figure 2** demonstrates how the interaction of the duration of sexual distancing and service interruptions impacts the standardized cumulative incidence outcomes. Similar to **Figure 1**, the scenarios here reflected a 50% relative reduction in both behavior and services during the eligible change period. HIV and STI incidence were lowest (53 and 533 cases per 1,000 susceptible, respectively) when service interruption lasted for 3 months and distancing lasted for 18 months. Both HIV and combined STI incidence were highest (73 and 1,529 cases, respectively) when services were interrupted for 18 months but sexual distancing occurred for 3 months.

![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/10/21/2020.09.30.20204529/F2.medium.gif)

[Figure 2.](http://medrxiv.org/content/early/2020/10/21/2020.09.30.20204529/F2)

Figure 2. 
Relationship between the duration of sexual distancing and service interruption on cumulative (5-year) incidence of HIV and STIs per 1000 susceptible MSM. Individual box and whiskers display the distribution of cumulative incidence across 500 simulations within each scenario.

Finally, we explored that last scenario further in **Figure 3** and the bottom of **Table 3**. Here the resumption time for sexual behavior was varied independently from the resumption time for services. For HIV, 3-month distancing in the absence of service change had no substantive impact on the trajectory of HIV incidence (green line), and therefore was unable to counterbalance the effects of service interruption in the combined scenario (blue line). Higher point incidence lasted through the end of 5 years. Over 5 years, this resulted in 890 excess HIV cases. For STIs, incidence in the combined scenario followed a similar pattern as for HIV, but to a more extreme level. STI incidence in the combined scenario (blue line) more closely tracked the scenario with service interruption only (red line), and STI incidence did not return to baseline values through 5 years. This resulted in an excess of 57,500 cases over 5 years.

![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/10/21/2020.09.30.20204529/F3.medium.gif)

[Figure 3.](http://medrxiv.org/content/early/2020/10/21/2020.09.30.20204529/F3)

Figure 3. 
HIV and combined STI incidence before, during, and after an 18-month period of clinical service interruption and 3-month period of sexual distancing. Model scenarios within each panel compare a 50% relative reduction in behavior only, service interruption only, or both jointly against the base (no change) model. Thick lines show median values and bands show inter-quartile range of values across 500 simulations per scenario.

## DISCUSSION

This study projected the 5-year impact of COVID-related sexual behavior changes (“sexual distancing”) and interruption of clinical services on the incidence of HIV, gonorrhea, and chlamydia among MSM in the Atlanta area. We found that the magnitude and timing of epidemiological impact depended on the infection and on the relative extent and durability of the COVID-related changes. Durable sexual distancing could offset (for HIV) or overcome (for STIs) the excess incidence attributable to an equal period of clinical service interruption. However, sexual activity rebounding while service interruption persists would lead to higher incidence for both HIV and STIs. Based on current estimates of behavioral and clinical change, and future predictions of an 18-month clinical service interruption duration, we project an excess of nearly 900 HIV cases and over 57,000 STI cases just among Atlanta MSM over the next 5 years. Our findings suggest that immediate action is needed to address the indirect effects of the COVID-19 pandemic on the HIV/STI epidemic.

The protective effects of sexual distancing on both HIV and STI incidence highlight the need to understand and address the extent and durability of behavior change during the COVID-19 pandemic. In our model, sexual partner reduction had a protective effect on HIV and STI incidence. Empirical data suggests a wide heterogeneity of changes in sexual activity during initial COVID restrictions in March 2020.7,8 However, fewer studies have characterized the timing of behavioral rebounds; some have suggested a return to pre-COVID levels starting as early as June 2020.9,10 Our model suggests that such a transient change will have no substantive impact for HIV and only a minor impact for STIs. The epidemiological impact of behavioral responses to COVID-19 highlight the need for effective sexual health messaging that supports sustainable behavior changes throughout the pandemic, which could have a protective effect on both HIV/STIs and SARS-CoV-2 transmission. Local health departments, such as the New York City Department of Health, have developed such guidelines.28 Sexual health messaging should adopt a harm reduction framework.29

The projected detrimental effects on HIV and STI incidence of clinical service interruption in the model demonstrate the critical importance of maintaining sexual health services amidst the COVID-19 pandemic response. Our model evaluated four clinical services for which there was already evidence of interruption to the relative degree in our scenarios.7 In some jurisdictions, health department staff assigned to HIV/STI partner services have been redeployed for COVID contact tracing.13 Interruption of ART care for persons living with HIV had the largest impact on projected excess HIV incidence in our model. This finding is consistent with two other models of the impact of COVID-19 on HIV outcomes, focused on low/middle-income countries.1,30 Minimizing service interruption will require innovative approaches to ensure access to clinical services and overcome common barriers to care during the COVID pandemic, including travel limitations and gaps in health insurance. These approaches will remain important even as sexual health services return to pre-COVID capacity and long-lasting impacts on health care access affect re-engagement in services.

### Limitations

The primary limitation of our model projections is the future uncertainty about the types and durability of both sexual distancing and service reduction. While our sensitivity analysis (**Figure 2**) explored durability, both components may change or rebound in different ways from those modeled. Modeling changes in the COVID-19 era required mapping data on broader aggregate reductions onto individual rate-based model parameters. Future empirical research on sexual distancing and service interruption should measure these changes with more individual-level specificity within and across persons. Second, we assumed that there was no correlation between changes in individual behavior and changes in service engagement. This decision was based on a re-analysis of one published study,7 provided in **Supplemental Table 4**.**1**, showing limited evidence of such individual-level correlation between change in sexual partner numbers and change in service access (small correlations were observed for HIV and STI screening only). Additionally, some group-level age-related changes in both behavior and services have been noted, but we did not include that heterogeneity in the model because it did not emerge in individual-level correlations. Third, our model did not explicitly represent the transmission dynamics of COVID-19 that may result in changes to sexual networks based on real or perceived COVID risk (e.g., partner selection by COVID status or risk factors, or other disruptions to travel to meet sex partners that could result in concentrated transmission clusters). This is an important topic for future data collection and modeling studies. Finally, our target population was MSM in the Atlanta area, a population with lower baseline access to HIV/STI services and a higher baseline of infection than other populations and areas. Therefore, our standardized results may not be scalable to populations with a different epidemiological context.

### Conclusions

The global COVID-19 pandemic could present substantial challenges to the prevention and control of other infectious diseases, including HIV and STIs, as a result of clinical service interruptions. While sexual behavior change may offset service interruptions, this will depend on overlap in the timing of these changes. With transient behavior change but persistent service interruptions, we project major increases in HIV and STI incidence that may take years to return to pre-COVID levels. This calls for improved sexual health messaging and innovative approaches to addressing service gaps during the ongoing COVID-19 pandemic.

## Supporting information

Supplemental Appendix [[supplements/204529_file02.pdf]](pending:yes)

## Data Availability

Model code and data are available at the Github repositories linked in the Supplemental Appendix document.

## Footnotes

*   **FUNDING** This work was supported by grants from National Institutes of Health grant (R01 AI138783 and P30 AI050409) and the MAC AIDS Fund. JLM is supported in part by the National Institute of Allergy and Infectious Diseases (grant K01 AI122853).

*   **DISCLOSURES** JLM has previously consulted for Kaiser Permanente Northern California on a research grant from Gilead Sciences. DW has consulted for Sanofi-Pasteur on unrelated topics.

*   Edited text throughout for clarity. Fixed plotting of temporal outcomes in Figures 1 and 3.

*   Received September 30, 2020.
*   Revision received October 21, 2020.
*   Accepted October 21, 2020.


*   © 2020, Posted by Cold Spring Harbor Laboratory

This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/)

## REFERENCES

1.  Hogan AB, Jewell BL, Sherrard-Smith E, et al. Potential impact of the COVID-19 pandemic on HIV, tuberculosis, and malaria in low-income and middle-income countries: a modelling study. Lancet Glob Health 2020; 8: e1132– 41.
    
    

2.  Andersen M. Early Evidence on Social Distancing in Response to COVID-19 in the United States. 2020; published online April 5. DOI:10.2139/ssrn.3569368.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2139/ssrn.3569368&link_type=DOI) 

3.  Coombe J, Kong F, Bittleston H, et al. Love during lockdown: findings from an online survey examining the impact of COVID-19 on the sexual practices of people living in Australia. medRxiv 2020; : 2020.08.10.20171348.
    
    

4.  Hensel DJ, Rosenberg M, Luetke M, Fu T, Herbenick D. Changes in Solo and Partnered Sexual Behaviors during the COVID-19 Pandemic: Findings from a U.S. Probability Survey. medRxiv 2020; published online June 12. DOI:10.1101/2020.06.09.20125609.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMC4wNi4wOS4yMDEyNTYwOXYyIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMTAvMjEvMjAyMC4wOS4zMC4yMDIwNDUyOS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

5.  Spinner T. Report from the Field: The Impact of COVID-19 on Local Health Department HIV, STI, and Hepatitis Programs - NACCHO. 2020; published online April 15. [https://www.naccho.org/blog/articles/report-from-the-field-the-impact-of-covid-19-on-local-health-department-hiv-sti-and-hepatitis-programs](https://www.naccho.org/blog/articles/report-from-the-field-the-impact-of-covid-19-on-local-health-department-hiv-sti-and-hepatitis-programs) (accessed Sept 29, 2020).
    
    

6.  Centers for Disease Control and Prevention. STDs in Men Who Have Sex with Men - 2018 Sexually Transmitted Diseases Surveillance. 2019; published online Oct 8. [https://www.cdc.gov/std/stats18/msm.htm](https://www.cdc.gov/std/stats18/msm.htm) (accessed Sept 22, 2020).
    
    

7.  Sanchez TH, Zlotorzynska M, Rai M, Baral SD. Characterizing the Impact of COVID-19 on Men Who Have Sex with Men Across the United States in April, 2020. AIDS Behav 2020; published online April 29. DOI:10.1007/s10461-020-02894-2.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s10461-020-02894-2&link_type=DOI) 

8.  McKay T, Henne J, Gonzales G, Quarles R, Gavulic KA, Gallegos SG. The COVID-19 Pandemic and Sexual Behavior among Gay and Bisexual Men in the United States. SocArXiv 2020; published online May 30. DOI:10.31235/osf.io/8fzay.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.31235/osf.io/8fzay&link_type=DOI) 

9.  Pampati S, Emrick K, Siegler A, Jones, Jeb. Changes in sexual behavior and access to sexual health services due to the COVID-19 pandemic among PrEP-using MSM in the South. In: STD Prevention Conference 2020. Atlanta, GA, 2020.
    
    

10. Stephenson R, Chavanduka TMD, Rosso MT, et al. Sex in the Time of COVID-19: Results of an Online Survey of Gay, Bisexual and Other Men Who Have Sex with Men’s Experience of Sex and HIV Prevention During the US COVID-19 Epidemic. AIDS Behav 2020; published online Sept 2. DOI:10.1007/s10461-020-03024-8.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s10461-020-03024-8&link_type=DOI) 

11. Rogers BG, Coats CS, Adams E, et al. Development of Telemedicine Infrastructure at an LGBTQ+Clinic to Support HIV Prevention and Care in Response to COVID-19, Providence, RI. AIDS Behav 2020; 24: 2743–7.
    
    

12. Jaffe DH, Lee L, Huynh S, Haskell TP. Health Inequalities in the Use of Telehealth in the United States in the Lens of COVID-19. Population Health Management 2020; published online Aug 18. DOI:10.1089/pop.2020.0186.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1089/pop.2020.0186&link_type=DOI) 

13. Centers for Disease Control and Prevention. Guidance and Resources During Disruption of STD Clinical Services. 2020; published online Sept 8. [https://www.cdc.gov/std/prevention/disruptionGuidance.htm](https://www.cdc.gov/std/prevention/disruptionGuidance.htm) (accessed Sept 21, 2020).
    
    

14. Jenness SM, Goodreau SM, Morris M. EpiModel: An R Package for Mathematical Modeling of Infectious Disease over Networks. J Stat Softw 2018; 84: 1–47.
    
    

15. Krivitsky PN, Handcock MS. A Separable Model for Dynamic Networks. J R Stat Soc Series B Stat Methodol 2014; 76: 29–46.
    
    

16. Jenness SM, Weiss KM, Goodreau SM, et al. Incidence of Gonorrhea and Chlamydia Following Human Immunodeficiency Virus Preexposure Prophylaxis Among Men Who Have Sex With Men: A Modeling Study. Clin Infect Dis 2017; 65: 712–8.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/cid/cix439&link_type=DOI) 

17. Weiss KM, Goodreau SM, Morris M, et al. Egocentric sexual networks of men who have sex with men in the United States: Results from the ARTnet study. Epidemics 2020; 30: 100386.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.epidem.2020.100386&link_type=DOI) 

18. Horn T, Sherwood J, Remien RH, Nash D, Auerbach JD. Towards an integrated primary and secondary HIV prevention continuum for the United States: A cyclical process model. J Int AIDS Soc 2016; 19: 21263.
    
    

19. Georgia Department of Public Health, HIV/AIDS Epidemiology Section. HIV Surveillance Summary, Georgia 2017. 2019 [https://dph.georgia.gov/data-fact-sheet-summaries](https://dph.georgia.gov/data-fact-sheet-summaries) (accessed April 20, 2020).
    
    

20. Jenness SM, Maloney KM, Smith DK, et al. Addressing Gaps in HIV Preexposure Prophylaxis Care to Reduce Racial Disparities in HIV Incidence in the United States. Am J Epidemiol 2019; 188: 743–752.
    
    

21. Centers for Disease Control and Prevention. Preexposure prophylaxis for the prevention of HIV infection in the United States–2017 Update: A Clinical Practice Guideline. 2017 [https://www.cdc.gov/hiv/pdf/risk/prep/cdc-hiv-prep-guidelines-2017.pdf](https://www.cdc.gov/hiv/pdf/risk/prep/cdc-hiv-prep-guidelines-2017.pdf).
    
    

22. Kelley CF, Kahle E, Siegler A, et al. Applying a PrEP Continuum of Care for Men who Have Sex with Men in Atlanta, GA. Clin Infect Dis 2015; 61: 1590–7.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/cid/civ664&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26270691&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F21%2F2020.09.30.20204529.atom) 

23. Grant RM, Anderson PL, McMahan V, et al. Uptake of pre-exposure prophylaxis, sexual practices, and HIV incidence in men and transgender women who have sex with men: a cohort study. Lancet Infect Dis 2014; 14: 820–9.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1473-3099(14)70847-3&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25065857&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F21%2F2020.09.30.20204529.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000341342700025&link_type=ISI) 

24. Chan PA, Mena L, Patel R, et al. Retention in care outcomes for HIV pre-exposure prophylaxis implementation programmes among men who have sex with men in three US cities. J Int AIDS Soc 2016; 19: 20903.
    
    

25. 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: 48–56.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1097/00126334-200005010-00009&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10877495&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F21%2F2020.09.30.20204529.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000087734400007&link_type=ISI) 

26. Kent CK, Chaw JK, Wong W, et al. Prevalence of rectal, urethral, and pharyngeal chlamydia and gonorrhea detected in 2 clinical settings among men who have sex with men: San Francisco, California, 2003. Clin Infect Dis 2005; 41: 67–74.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1086/430704&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15937765&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F21%2F2020.09.30.20204529.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000229530400011&link_type=ISI) 

27. Rosenberg ES, Grey JA, Sanchez TH, Sullivan PS. Rates of Prevalent HIV Infection, Prevalent Diagnoses, and New Diagnoses Among Men Who Have Sex With Men in US States, Metropolitan Statistical Areas, and Counties, 2012-2013. JMIR Public Health Surveill 2016; 2: e22.
    
    

28. New York City Department of Health and Mental Hygiene. Safer Sex and COVID-19. 2020 [https://www1.nyc.gov/assets/doh/downloads/pdf/imm/covid-sex-guidance.pdf](https://www1.nyc.gov/assets/doh/downloads/pdf/imm/covid-sex-guidance.pdf).
    
    

29. Marcus J. Americans Aren’t Getting the Advice They Need. The Atlantic. 2020; published online May 28. [https://www.theatlantic.com/ideas/archive/2020/05/no-one-telling-americans-how-reopen-their-lives/612172/](https://www.theatlantic.com/ideas/archive/2020/05/no-one-telling-americans-how-reopen-their-lives/612172/) (accessed Sept 29, 2020).
    
    

30. Jewell BL, Mudimu E, Stover J, et al. Potential effects of disruption to HIV programmes in sub-Saharan Africa caused by COVID-19: results from multiple mathematical models. Lancet HIV 2020; 7: e629–40.