Article Text

Original research
Effect of screening young men for Chlamydia trachomatis on the rates among women: a network modelling study for high-prevalence communities
  1. Zhuolin Qu1,2,
  2. Asma Azizi3,
  3. Norine Schmidt4,
  4. Megan Clare Craig-Kuhn4,
  5. Charles Stoecker5,
  6. James Mac Hyman1,
  7. Patricia J Kissinger4
  1. 1Department of Mathematics, Tulane University, New Orleans, Louisiana, USA
  2. 2Department of Mathematics, The University of Texas at San Antonio, San Antonio, Texas, USA
  3. 3Simon A Levin Mathematical Computational Modeling Science Center, Arizona State University, Tempe, Arizona, USA
  4. 4Department of Epidemiology, Tulane University, New Orleans, Louisiana, USA
  5. 5Department of Health Policy and Management, Tulane University, New Orleans, Louisiana, USA
  1. Correspondence to Professor Patricia J Kissinger; kissing{at}tulane.edu

Abstract

Objective Chlamydia trachomatis (Ct) is the most commonly reported sexually transmitted infection in the USA and causes important reproductive morbidity in women. The Centers for Disease Control and Prevention recommend routine screening of sexually active women under age 25 but not among men. Despite three decades of screening women, chlamydia prevalence in women remains high. Untested and untreated men can serve as a reservoir of infection in women, and male-screening based intervention can be an effective strategy to reduce infection in women. We assessed the impact of screening men on the Ct prevalence in women.

Design We created an individual-based network model to simulate a realistic chlamydia epidemic on sexual contact networks for a synthetic population (n=5000). The model is calibrated to the ongoing routine screening among African American (AA) women in the USA and detailed a male-screening programme, Check It, that bundles best practices for Ct control. We used sensitivity analysis to quantify the relative importance of each intervention component.

Setting Community-based venues in New Orleans, Louisiana, USA.

Participants Heterosexual AA men, aged 15 to 24, who had sex with women in the past 2 months.

Intervention Venue-based screening, expedited index treatment, expedited partner treatment and rescreening.

Results We estimate that by annually screening 7.5% of the AA male population in the age-range, the chlamydia prevalence would be reduced relatively by 8.1% (95% CI 5.9% to 10.4%) in AA women and 8.8% (95% CI 6.9% to 10.8%) in AA men. Each man screened could prevent 0.062 (95% CI 0.030 to 0.094) cases in men and 0.204 (95% CI 0.143 to 0.267) cases in women. The model suggested the importance of intervention components ranked from high to low as venue-based screening, expedited index treatment, expedited partner treatment and rescreening.

Conclusion The findings indicated that male-screening has the potential to substantially reduce the prevalence among women in high-prevalence communities.

  • statistics & research methods
  • public health
  • infection control
http://creativecommons.org/licenses/by-nc/4.0/

This is an open access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited, appropriate credit is given, any changes made indicated, and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/.

Statistics from Altmetric.com

Request Permissions

If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.

Strengths and limitations of this study

  • The network-based model captures the complex assortative mixing in sexual partnerships among an African American (AA) population in New Orleans using the survey data on sexual behaviours.

  • We present a novel modelling study to assess the impact of a new male-screening programme that bundles the best practices for chlamydia control on the current chlamydia prevalence in AA women, given the ongoing routine screening among AA women in the USA.

  • The proposed stochastic, heterosexual and individual-based model provided a flexible framework for a public health team to answer ‘what if’ questions that are hard to address in the field.

  • The model is parametrised using data from two surveys on sexual behaviours, which reflects our best understanding of the AA population in the New Orleans area.

  • The model assumed a closed and stable population in an area, and the quantitative results are only valid for a short-term prediction that does not consider impacts by the external factors, such as behaviour changes, natural disasters and changes in mixing patterns.

Introduction

Chlamydia trachomatis (Ct) is the most commonly reported infectious disease in the USA, with over 1.7 million cases each year.1 It is a major cause of infertility, pelvic inflammatory disease and ectopic pregnancy among women2 and has been associated with increased HIV acquisition.3 Because women experience the most severe sequelae, the focus of Ct prevention in the USA has been on screening sexually active women <25 years old, providing her and her partner(s) with treatment and rescreening. Despite three decades of screening women, chlamydia prevalence remains high in the sexually active young women in the USA.1

There is no recommendation for Ct screening among men in the USA. In 2007, an expert panel at the Centers for Disease Control and Prevention (CDC) concluded that the evidence is not sufficient to recommend routine screening for Ct in sexually active young men.4 The conclusion was based on the Ct prevalence in 2007 and the feasibility, efficacy and cost-effectiveness of screening men. However, since then, evidence of the benefit of screening young men for Ct in high prevalence areas has been mounting. More recent modelling studies indicate that screening men in high prevalence populations can be cost-effective due to averted cases among women.5 6

The community-based programme, Check It,7 Ct screens African American (AA) men aged 15 to 24 in New Orleans. The core of this intervention is Ct screening for men. We hypothesise that men are an important reservoir of infection for women and therefore need to be targeted for intervention.7 The Check It programme bundles several key Ct control strategies (table 1):

  • Venue-based screening (VBS) of participants at non-clinical community venues, such as barbershops, colleges and universities, in high-prevalence neighbourhoods characterised by similar demographic and geographical factors. This venue-based enrolment is enhanced with marketing strategies, such as the distribution of flyers, web education, social media and informational cards.

  • Expedited treatment by providing medication for the Ct-positive men (or index) (expedited index treatment or EIT) and his sexual partner(s) (expedited partner treatment or EPT) via partnering community pharmacies without a medical examination to speed up the treatment of his sexual partners and reduce reinfection rates in the index.

  • Rescreening of Ct-positive men; retesting for infection 3 months after treatment.

  • Social network peer referral (SNPR) that encourages men to refer young AA men in their social network to Check It via flyers, social media or text messages to promote the programme and increase the total enrolment.

Table 1

Summary of interventions involved in the Check It programme

Mathematical models create frameworks for understanding the underlying epidemiology of disease and help test the potential effectiveness of different approaches to bring the epidemic under control. Our modelling effort created a detailed simulation to model the specific practices implemented by the Check It programme. Most chlamydia simulations studies use differential equation-based compartmental models. These models are parametrised at the population level and assume homogeneous mixing of the individuals in the same compartment. We formulated an agent-based model at a more granular level that captures the complex assortative mixing pattern in sexual partnership networks. The structure of the sexual network affects the spread of the infection and the effectiveness of the mitigation efforts. We calibrated the agent-based model based on the sexual behaviour surveys7 8 and assessed the impact of these male-screening-based strategies on mitigating the Ct epidemic among AA women aged 15 to 24 years old in New Orleans.

Methods

We modelled the Ct transmission among a synthetic population connected through heterosexual partnership networks. These sexual networks were generated based on the data from two survey studies7 8 that investigated the sexual behaviour of young AA men and women in New Orleans. The model parameterised the transmission pathways on the individual level, and it simulated the sexual behaviour and kept track of the infection status for each individual in the synthetic population over time. We then modelled the intervention strategies, including both the standard preventive healthcare for women (routine Ct screening) and the Check It intervention, for each individual to study the impact of male-screening on the Ct prevalence in women.

Generation of a synthetic population over dynamic sexual networks

We constructed a closed 5000-member population, where the heterosexual partnerships were represented by bipartite sexual networks. The heterosexual networks captured the assortative mixing pattern among our targeted AA population by matching the population-level quantities (the degree distribution and joint-degree distribution, online supplemental appendix A.1) from two surveys: the ongoing Check It study7 and the ‘You Geaux Girl!’ (YGG) study,8 which enrolled 1318 AA men (May 2017 to April 2019, ongoing, age range 15 to 24) and 473 AA women (September 2012 to December 2015, completed, age range 18 to 19), respectively, in New Orleans. All participants gave written informed consent before taking part in the studies.

We assumed that the population is closed and men and women in our model only have sexual partnerships within this population cohort. This assumption closely agrees with our data sets, where about 92% of the male participants (Check It study) and 94% of the female participants (YGG study) have partners in the same age cohort (15 to 24 years old). The surveys also found that about 90% of the partners for male participants (Check It study) and 95% of the partners for female participants (YGG study) were also AA. We included the technical details of network descriptions in online supplemental appendix A.1.

To further characterise sexual behaviours, we categorised one’s partner as either primary or casual partner, depending on the survey responses to questions such as level of commitment and the duration of the relationship (detailed criteria in online supplemental appendix A.1.2 and the distributions in online supplemental appendix A.2). These categorisations can be asymmetric: A is B’s primary partner, but B may be A’s casual partner, and there are three types of partnerships in the networks: primary-primary, primary-casual and casual-casual partnership.

When modelling the epidemics over a long time (several months or years), we updated the sexual partner(s) every 2 months (the time-frame covered by the Check It survey) and created a series of dynamic sexual networks that are evolving in time. To simulate the partner updating behaviour, we employed a detailed realistic social contact network that modelled the daily activities for 150 000 people in New Orleans.9 Then the heterosexual networks were updated through people’s contacts in the social network. Specifically, we assumed that the primary-primary partnerships were preserved throughout the simulation period and half of the primary-casual and all of the casual-casual partnerships are replaced by one’s social contacts every 2 months (details in online supplemental appendix A.1.3). There were, on average, 80% of the sexual partnerships coming from one’s heterosexual social network, which is in agreement with the Check It data (76% of partners were reported as social contacts from school, work, the neighbourhood and so on).

Chlamydia epidemic on dynamic sexual networks

We modelled the infection status of each individual using the Susceptible-Infectious-Susceptible (SIS) framework. All uninfected individuals are susceptible to being infected, and all infected people recover to this susceptible state after either spontaneous recovery or treatment. A susceptible person can be infected by his/her infectious partner. The force of infection is the probability (per day) that a susceptible person will be infected. This probability was estimated by considering risk factors of how many partners the person has, the type of partnership with each partner, the probability of having sexual contact per partner per day (contact rate) and the probability of using a condom.

We estimated the sexual contact rates for the primary and casual partnerships from the data sets (see table A.3). We observed that the contact rate was higher for primary partnerships than for casual partnerships, and there was a decreasing trend in the per partner contact rate when the number of partners increases. Moreover, the data sets gave a higher probability of using a condom with a casual partner than with a primary one (table 2). We introduced a condom failure rate to include cases when the condom is not used properly.

Table 2

Model parameters for the AA population aged 15–24

The details on the configuration of the epidemic model over networks are fully described in (online supplemental appendix A.2).

Modelling the intervention strategies

Our goal was to investigate the net impact of screening men through the Check It programme given the existing screening policy for women and the ongoing endemic Ct epidemic. To this end, we outlined the intervention strategies for both women and men in separate flow charts in figure 1. The baseline scenario accounts for Ct screenings completed at women’s annual exams as part of regular preventive healthcare. We also included the Ct screenings prompted by symptomatic infections (clinical visits) in the baseline scenario for both men and women.

Figure 1

Flowchart for Ct intervention strategies in men and women. The solid lines are the new practices incorporated in the male-screening programme, Check It, and the dashed lines are the existing interventions implemented in the healthcare system. These current interventions include women’s annual screening and screenings prompted by symptomatic infections in both men and women. Our modelling effort assessed the net impact of the male-screening programme to help control the Ct epidemic. The Check It programme targets the male population and uses venue-based enrolment, expedited index treatment, expedited partner treatment, rescreening and social network peer referral (see table 1). The intervention parameters are marked along the routes indicating the rates of compliance and delays, obtained from either literature or Check It data (see table 2). Ct,Chlamydia trachomatis.

Existing intervention strategies

The current (baseline) Ct mitigation efforts for women (the right side of figure 1) include Ct screening during routine annual exams and clinical visits for symptomatic infections. The model assumes that the same fraction (Embedded Image) of women return for a physical exam each year (more details in online supplemental appendix A.3). Symptoms can appear in a small fraction of infected women (Embedded Image), and we assumed that a fraction Embedded Image of these women get medical care within an average of Embedded Image days, including incubation period and appointment scheduling process, after infection.10 After the diagnostic test, we assumed that all the positive cases get index treatment in an average of Embedded Image days. Moreover, the CDC recommends EPT for infected women by providing treatment to the patient to bring to her partner(s) without first examining the partner(s).11 We assumed that a fraction Embedded Image of the partners are treated with an average delay of Embedded Image days. The EPT fraction Embedded Image is the product of (1) the fraction of the physicians practising EPT as recommended and (2) the fraction of compliance from the notified partners. Last, women diagnosed with Ct infection should be retested after the initial treatment.12 Thus, we assumed a fraction Embedded Image of treated women are retested for infection Embedded Image days after the initial treatment.

Since routine male-screening is not recommended and is rarely practised,13 we only account for screenings prompted by symptomatic infections. Similar to the process in women, we assumed that there is a small fraction (Embedded Image) of infected men that develop symptomatic Ct infections, and, with a delay of Embedded Image days, a fraction Embedded Image of them get screening and treatment. Moreover, follow-up interventions, such as partner treatment or rescreening after the index treatment, are not commonly implemented in clinics for men.13 Therefore, follow-up interventions for men were not included in the model.

Check It intervention strategies

The Check It programme recruited participants in community venues, including community colleges, historically black colleges and universities, barbershops and other community-based organisations. This venue-based enrolment is enhanced with marketing strategies (distribution of flyers, web education, social media and informational cards) and enrols a fraction Embedded Image of the target male population for Ct screenings.

Some participants learnt about the programme through their social networks, such as text messages and information cards sent by friends or word-of-mouth. We accounted for this peer impact as SNPR and included it as a source of enrolment (see table A.4). On average, the proportion between the non-peer enrolled men and peer-recruited men is 1:ρ, and we assumed that peer-referred men are enrolled in the programme with an average delay of Embedded Image days from the time the referring man enrolled.

We modelled the non-peer enrolment process as a random sampling from the entire male population. Meanwhile, we modelled the peer-referred enrolment (SNPR) by searching the background social network of each non-peer enrolled man and randomly sampling among the eligible candidates.

The rest of the intervention practice was modelled similarly to women’s cases: a fraction Embedded Image of the screened and infected men receive EIT after a delay of Embedded Image days. A fraction Embedded Image of these men’s partners receive EPT with a delay of Embedded Image days. Finally, a fraction Embedded Image of the index men return for rescreening Embedded Image days after the initial infection. If the rescreened men are infected with Ct, then they are treated as index cases, and the intervention process is repeated. Details on the estimates of compliance rates and delays were summarised in table A.5 and A.6.

Model parameters and calibration

The model parameters used in the simulations (table 2) represent our best knowledge of the current Ct epidemic and mitigation efforts among the AA population in the age range 15 to 24 in the USA. The values for the parameter Ct transmission probabilities per sexual contact between men and women (Embedded Image and Embedded Image) are not clear in the current biological literature. Quinn et al14 gives the estimate of infection frequency per partnership. There are modelling studies15 16 that provide rough estimates with large variations on the probability of transmission per contact. Because of the large uncertainty with these two parameters, they may vary in a wide range and are treated as tuning parameters for model calibration.

We calibrated the model to fit the current Ct prevalence in New Orleans among the AA population aged 15 to 24 (10.2% in men and 13.5% in women).17 This prevalence reflects both the high infection rate in the region and the ongoing mitigation efforts among women. To have a realistic initial infection within the sexual network, we set the initial infection population to be consistent with the distribution of an emerging epidemic. We obtained such a quasi-steady-state balanced initial condition by starting a small epidemic in the past and letting it grow to the current (pre-Check It) Ct endemic state (online supplemental appendix A.4). This initialisation process considered the existing (baseline) combination of Ct interventions (dashed routes in figure 1) to give a comprehensive approximation of the current Ct control before the launch of Check It.

Approximately half of the Ct infections in women are cleared naturally by the first year after being infected and 80% are cleared after 2 years.18 We fitted an exponential distribution for the average time of natural recovery Embedded Image days to the data from,18 and we assumed the same distribution for men and women. Our numerical simulations suggest that the model results are not sensitive to this assumption (see Discussion and online supplemental appendix A.5). We also assumed the recovery time with treatment also follows an exponential distribution Embedded Image in days.

Sensitivity analysis

The model parameters in table 2 represent the best-guess estimates for practical scenarios, and we used local and extended sensitivity analysis to quantify the most significant model parameters.19 To check the impact of each component of the Check It intervention bundle one-at-a-time, we conducted the extended sensitivity analysis by varying each intervention parameter (the parameter of interest), while fixing the other parameters. We then checked its corresponding impact on the Ct prevalence for women, men and the overall population (the quantities of interest).

For the local sensitivity analysis, we defined the relative local sensitivity index of a quantity of interest, q, with respect to the parameter of interest, p, as Embedded Image. This normalised sensitivity index, Embedded Image, measures the percentage change in an output quantity given the percentage change in an input parameter. If an input parameter, p, changes by Embedded Image, then the output quantity, q, changes by Embedded Image.

To further investigate the synergistic effect of the intervention components beyond the current levels, which could be limited by the protocols and available resources of Check It, we conducted the global sensitivity analysis by varying two intervention parameters together while fixing all other parameters. We then predicted the impact under different combinations of intervention parameters.

Patient and public involvement

There was no patients or public involvement in designing, conducting, reporting on or dissemination of information related to our modelling study.

Results

Impact of male-screening programme at existing intervention level

We quantified the impact of the male-screening programme implemented at the current level of intervention intensity (as shown in table 2). Figure 2 shows the change in Ct prevalences after the launch of the male-screening intervention for men at year 0 with a balanced initial condition (online supplemental appendix A.4). Around year 5 of the programme, the Ct prevalences are controlled at much lower levels (quasi-steady states): 12.4% (95% CI 12.1% to 12.7%) in women and 9.3% (95% CI 9.1% to 9.5%) in men, which are reduced relatively by 8.1% (95% CI 5.9% to 10.4%) in women and 8.8% (95% CI 6.9% to 10.8%) in men. When stopping the male-screening programme later around year 14, the Ct prevalences return back to the baseline scenario in about 5.5 years.

Figure 2

Impact of male-screening programme implemented at the existing intervention level. The curves are the mean of 50 stochastic simulations, and the bands around the curves indicate the one SD. The baseline Ct prevalences (before year 0) are 13.5% and 10.2% in women and men, respectively. At year 0, the male-screening intervention is turned on. Around year 5, the Ct prevalences reach quasi-steady states, which are 12.4% in women and 9.3% in men. The prevalences are reduced by 8.1% in women and 8.8% in men relatively. When the male-screening programme is stopped around year 14, the Ct prevalences return to the baseline levels in about 5.5 years. Ct,Chlamydia trachomatis.

In fact, near the lower quasi-steady states, our model predicts the following annual statistics from the programme based on a population of 5000:

  • Each year, the programme conducts in total 174 screenings, including 42 from peer-recruited participants (SNPR), and achieves 13 treatments for index men (EIT) and 8 treatments for partners of those men (EPT).

  • Among all the screened men found to be Ct-positive, the average number of partners within the past 2 months is 2.30 (95% CI 2.25 to 2.35).

  • Compared with the scenario without Check It, the programme prevents 10.8 (95% CI 5.3 to 16.3) cases in men and 35.6 (95% CI 24.8 to 46.4) cases in women per year.

  • Roughly, for each man screened, it could prevent 0.062 (95% CI 0.030 to 0.094) cases in men and 0.204 (95% CI 0.143 to 0.267) cases in women.

Significance of the components of the intervention

The sensitivity analysis quantified the relative significance of the intervention components in the programme. In the results presented below, we have considered the quantity of interest to be Ct prevalence and omit the upper index in the sensitivity index for the simplicity of the presentation.

From figures 3 and 4 (left), the Ct prevalences have an almost linear response to the intervention parameters, and the sensitivity at the current level of Check It intensity is ranked from high to low as VBS, EIT, EPT, rescreening, where the sensitivity of VBS is close to EIT. Moreover, when increasing the coverage of VBS from 7.5% to 40% (figure 4), the magnitude of the corresponding sensitivity index for EPT becomes seven times larger (-0.117 vs −0.017), which suggests that EPT would be much more effective in reducing the prevalence with high male-screening coverage.

Figure 3

Local and extended sensitivity analysis on Check It intervention parameters (on the x-axis) against the Ct prevalence (y-axis). For each plot, the parameter of interest is varied while the other model parameters are fixed as in table 2. The Ct prevalences for men, women and the entire population are plotted, which are averaged over the time-frame year 4.5~5.5 of 30 simulations. The error bars give the one SD above and below the average. The local sensitivity indices (q=quantity of interest is Ct prevalence) are given in the titles. Ct,Chlamydia trachomatis; EIT,expedited index treatment; VBS, venue-based screening.

Figure 4

Local and extended sensitivity analysis on expedited partner treatment (on the x-axis) against the Ct prevalence (y-axis). Left: the analysis at the current level of Check It intervention intensity, screening 7.5% of the target male population. Together with the results in figure 3 and at the current Check It level, the significance of intervention components is ranked as venue-based screening ≈expedited index treatment >expedited partner treatment >rescreening. Right: the analysis at a much higher 40% male-screening rate while fixing other intervention parameters. The magnitude of the local sensitivity index is almost seven times larger (-0.117 vs −0.017), which suggests that the partner treatment becomes more important in reducing prevalence when increasing the screening coverage in men. Ct,Chlamydia trachomatis; EIT,expedited index treatment.

We then conducted a global sensitivity analysis using the two most significant parameters from the local and extended sensitivity analysis: VBS and EPT. The response plot for women’s Ct prevalence (figure 5) shows that the male-screening strategy has the potential to reduce the Ct prevalence in women substantially, and it predicts the effectiveness under different combinations of intervention intensities. For example, the model estimates the combination of VBS=30% of the target men and EPT=40% of their partners will give a 30% reduction in Ct prevalence among women after 5 years of intervention.

Figure 5

Global sensitivity analysis of Ct prevalence in women (marked in contour lines) against two intervention parameters: venue-based screening (VBS, x-axis) and expedited partner treatment (EPT, y-axis) on a uniform 5×5 grid. At each grid point, the Ct prevalence is averaged over the time-frame year 4.5~5.5 of 10 simulations, and the contour surface is smoothed by a least-square fit of a two-dimensional quadratic polynomial to the grid values. The baseline scenario (VBS=0% and EPT=0%) and current intervention level (VBS=7.5% and EPT=27%) are marked in crosses, which shows an 8.1% relative reduction in women’s prevalence. To achieve a 30% relative reduction in women’s Ct prevalence, the combined intervention levels required are marked by the dashed contour line. For example, with a coverage of 30% male-screening and 40% partner treatment, the model predicts that the Ct prevalence in AA women will be reduced to 9.45%. AA,African American; Ct, Chlamydia trachomatis.

Discussion

Our model provides a framework for public health workers to ask ‘what if’ questions that are hard to evaluate in the field. We simulated the Ct epidemic over sexual networks based on a young AA population in New Orleans, where the ongoing female-screening interventions are not sufficient to bring down the high Ct prevalence in women. From the simulation results at the current intervention level (figure 2), by annually screening men in the target population (AA, sexually active, aged 15 to 24), together with other best practices in Ct control, there is the potential to further mitigate Ct prevalence from the baseline (female-screening intervention only) among this high-prevalence cohort. Moreover, once the male-screening programme is suspended, the prediction shows that the prevalence will return to the baseline level given no significant changes in the underlying force of infection in the population. This suggests that, similar to the existing recommendation of annual screening policy for high-risk women, the male-screening intervention also requires long-term efforts to maintain its effectiveness over time. With the joint force of the interventions for both men and women, Ct could be controlled at much lower prevalences.

From the local and extended sensitivity analysis (figure 3), we quantified the impact of each intervention component. For all the intervention components involved, the Ct prevalences have linear responses to the variation of each intervention parameter. The sensitivity indices quantify the impact and show that the two most significant intervention parameters are (1) the coverage of venue-based screening for men, which identifies and treats the Ct-positive men, and (2) the coverage of expedited partner treatment, which prevents the potential reinfection between the couple. Figure 4 shows the impact of partner treatment on the Ct prevalences but with much higher coverage of male-screening among the target population (increase from 7.5% to 40%). The larger magnitude in the sensitivity index implies that the expedited partner treatment will be much more effective in reducing Ct prevalence when increasing the screening coverage in men.

By using the global sensitivity analysis on these two most significant intervention parameters, the model further explores the potential impact on a larger parameters space. Figure 5 gives the prediction of the Ct prevalences under different combinations of intervention parameters. At the current Check It intervention level, the impact could be limited by the resources and the capacity of the programme and the policy imposed on the protocol. In the adapted Check It protocol (as of February 2020), the expedited partner treatment incorporated the practices of patient delivered partner therapy and mail delivery, which improved the fraction of partner treatment from 27% to nearly 40%. Moreover, if male-screening could be part of the standard recommendation as it is for women and if a moderate compliance rate of 30% could be achieved, then the model predicts that this joint screening programme will result in 30% relative reduction in the Ct prevalence for the AA young women. After three decades of routinely screening women with no improvement in Ct rates, alternative interventions must be considered. This model provides support for changing the recommendations to include routine screening for men in addition to the current recommendations and best Ct control practices like EPT.

Our model considered a comprehensive picture of the current Ct epidemic in a young AA community and the control strategies implemented in men and women, separately. The model is parametrised using the survey data from the target population in New Orleans. This allows a detailed reconstruction of the assortative sexual mixing among the population and gives realistic predictions on the spread of the Ct epidemic and assessments on the mitigation efforts.

Although the model is parametrised based on the data sets from New Orleans, they represent similar cities that have high Ct rates (urban, southern and largely impoverished AA community) in many ways. Our study presented a robust and flexible model framework that could be adapted for the Ct epidemics in other similar populations. Nevertheless, the quantitative results presented in this study can only be interpreted for other cities after carefully examining the differences between sexual behaviour and assortative mixing patterns in two populations.

Many of our model limitations are associated with the scope of the available data sets. The uncertainty and bias in the model parameters and model assumptions can affect the reliability of the quantitative predictions. For example, due to the ethical challenges, rigorous estimates from longitudinal studies on the natural history of chlamydia are not available.20 We have assumed the same natural recovery period for men and women, and we studied the impact of this uncertainty in the (online supplemental Appendix A.5). Our assumptions on the partner updating frequency were made due to a lack of longitudinal data on the subject. Still, our numerical simulations (in online supplemental appendix A.1.3) suggested that, on the model recalibration, the relative ranking of the responses in the sensitivity analysis and the trend in prediction pattern is robust and insensitive. Also, Check It is still an ongoing study, and some of our model parameters may be further refined using the incoming data.

The dynamic sexual networks approximated the sexual behaviour between the high-risk age cohorts from Check It and YGG studies (men 15 to 24 and women 18 to 19). Meanwhile, we noticed a minority of women (6%) in the YGG study reported partners from different age cohort (>24 years old). These partnerships, which may serve as a residual reservoir, could not be detected by Check It programme that targets age range 15 to 24. However, due to the lack of further data on quantifying these sexual behaviours, we did not include the intergenerational partnership with large age differences.

In the current study, we did not incorporate the demographic characteristics of the individual. This may lead to some bias in modelling the mixing (eg, age-mixing) structure for the sexual partnership. We reduced this bias by embedding the sexual networks into a grand social network, which considers an individual’s age, ethnicity, social groups, economic status and geographical location. When updating the sexual networks, people will be more likely to establish sexual partnerships with their frequent social contacts. Another model assumption due to the lack of age is we didn’t include the ageing effect: people will not be removed as they get older than eligible age range, and no new susceptible young people will come into the population. Thus, the caveat is the model will only be a good approximation for a limited time, and it is not suitable for simulations over a long period. Our future modelling study will be focussed on improving the model, such as considering the ageing process in the population and accounting for sexual partnerships outside the current range. These improvements will help us better quantify the impact of a male-screening programme.

We recognise that mathematical models are a simplistic representation of the real world. The quantitative predictions from our model may need further validation using incoming data stream on Ct prevalences, however, the qualitative results of our analysis adds to the evidence that male-screening interventions, together with other best practices in Ct control, has the potential to further mitigate the Ct prevalence among women in a high-prevalence community. Our future work will include further model calibration and validation using the incoming data from the Check It programme and evaluating the cost-effectiveness of the male-screening programme.

Acknowledgments

We thank William Geisler (University of Alabama at Birmingham) for his helpful insight into the natural history of Chlamydia in men and women. We are also grateful to M Jacques Nsuami (Louisiana State University Health Sciences Center) whose comments and suggestions helped improve and clarify this manuscript.

References

Supplementary materials

  • Supplementary Data

    This web only file has been produced by the BMJ Publishing Group from an electronic file supplied by the author(s) and has not been edited for content.

Footnotes

  • Contributors All authors meet the ICMJE criteria for authorship. ZQ, AA and JH contributed to the design of the mathematical model. ZQ conducted data analysis, undertook numerical simulations and visualisation and interpreted results. NS and MCC-K contributed to data collection and analysis and reviewed model design and interpretation. CS reviewed the model design and interpretation. JH and PK contributed to the study, reviewed the model, data analysis and results interpretations and oversaw and coordinated the investigation. ZQ wrote the first draft of the article. All authors contributed to the writing and review of the draft and approved final manuscripts.

  • Funding This work was supported by the grants from the National Institutes of Health National Institute of Child Health and Human Development (R01HD086794) and the endowment for the Evelyn and John G Phillips Distinguished Chair in Mathematics at Tulane University. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.

  • Competing interests None declared.

  • Patient consent for publication Not required.

  • Ethics approval Tulane University Institutional Review Board.

  • Provenance and peer review Not commissioned; externally peer reviewed.

  • Data availability statement No data are available.

  • Supplemental material This content has been supplied by the author(s). It has not been vetted by BMJ Publishing Group Limited (BMJ) and may not have been peer-reviewed. Any opinions or recommendations discussed are solely those of the author(s) and are not endorsed by BMJ. BMJ disclaims all liability and responsibility arising from any reliance placed on the content. Where the content includes any translated material, BMJ does not warrant the accuracy and reliability of the translations (including but not limited to local regulations, clinical guidelines, terminology, drug names and drug dosages), and is not responsible for any error and/or omissions arising from translation and adaptation or otherwise.