Abstract
Accurate estimation of the sizes of most-at-risk populations (MARPs) for HIV infection is critical for effective intervention planning, resource allocation, and service coverage assessment. This study aimed to estimate the population sizes of female sex workers (FSWs), people who inject drugs (PWIDs), and men who have sex with men (MSM) in the 18–45-year-old population of Shiraz, Iran, using the network scale-up (NSU) method. A cross-sectional study was conducted involving 768 purposively selected individuals aged 18–45 years residing in Shiraz for at least 5 years. Data were collected through a self-administered questionnaire, utilizing 16 reference groups with known frequencies to estimate average social network size. The prevalence of MARPs was then estimated and adjusted using established correction factors. The average active social network size was 203 contacts (103 males and 100 females). Adjusted size estimations indicated the highest prevalence among FSW at 2,181 per 100,000 (95% confidence interval [CI]: 1,979–2,375), followed by PWID at 561 per 100,000 (95% CI: 488–633) and MSM at 418 per 100,000 (95% CI: 272–563). Age distribution analyses revealed the greatest PWID prevalence in the 35–39 age group, FSW in 30–34, and MSM in 25–29. The estimated sizes of these populations substantially exceeded current harm-reduction program coverage in the region. Application of the NSU method revealed significant gaps between the actual sizes of high-risk populations and those reached by existing interventions in Shiraz. Targeted expansion of prevention, testing, and treatment programs, alongside initiatives to combat stigma and improve outreach, is urgently needed to control HIV transmission among these hidden and vulnerable groups.
Introduction
In the Middle East and North Africa region, HIV infection is predominantly concentrated among subgroups exhibiting high-risk behaviors, including people who inject drugs (PWIDs), female sex workers (FSWs) and their clients, and men who have sex with men (MSM).1,2 These groups are classified as the most-at-risk populations (MARPs) for HIV. In countries such as Iran, where the HIV epidemic is concentrated, these subgroups play a pivotal role in the transmission of the virus to the general population.2–5
To effectively plan interventions, allocate resources, estimate the number of infected individuals, predict disease burden, measure service coverage, and evaluate health care services, reliable population estimates of these groups are essential at both local and national levels.6–9 Accurate measurement of the sizes of these populations is a critical component of HIV infection control, as highlighted in the Global CDC/UNAIDS and WHO HIV/AIDS Surveillance Meeting in Bangkok in March 2009. 9 However, due to the nature of their behaviors, along with stigmatization, discrimination, and fear of legal punishment, these groups remain hidden and are considered hard-to-reach populations. Consequently, accurate estimation of their sizes poses a significant challenge within national HIV/AIDS surveillance systems.2,5,8,10,11
Standard sampling methods and direct estimation techniques are not recommended for estimating the sizes of these hidden populations.12,13 Recent years have seen the development of various techniques and guidelines for measuring the sizes of these groups, including enumeration, census, capture-recapture, and network scale-up (NSU) methods. Each of these methods has its own advantages and disadvantages.6,9,14 Among them, NSU is one of the most significant and widely utilized methods, as it does not require direct contact with high-risk groups and allows for the estimation of multiple groups in a single study.11,14 However, this method is also subject to certain errors.15,16 The present study aims to estimate the sizes of MARPs for HIV among the 18–45-year-old population of Shiraz using the NSU method, while also refining existing errors.
Materials and Methods
Sampling and data collection
This cross-sectional study aimed to estimate the sizes of MARPs for HIV infection in Shiraz, the capital of Fars Province, utilizing the NSU method. The target population consisted of individuals aged 18–45 years who had resided in Shiraz for a minimum of 5 years. According to the latest National Population and Housing Census conducted by the Statistical Center of Iran, the total population of Shiraz is approximately 1,712,745, with 866,160 individuals falling within the 18–45 age range. A total of 768 eligible participants were interviewed using a self-administered questionnaire.
The requirement of at least 5 years of residence in Shiraz was applied to ensure comparability with previous NSU studies conducted in Iran and to increase the likelihood that respondents had developed stable and representative social networks within the study setting.8,10,14,17 Individuals with shorter residency may have smaller or less stable and geographically dispersed networks, which could introduce additional bias and increase transmission error in NSU-based estimations.
Participants were selected through purposive sampling from the 12 regional districts of Shiraz. Only those who provided written informed consent were included in the study. Given the limitations of random sampling in obtaining a representative sample of the target population, quota sampling was employed based on sex (50% male and 50% female) and socioeconomic status, with four regional districts allocated to each socioeconomic category.
To minimize sensitivity regarding certain questions, interviews were conducted by gender-matched trained interviewers in public areas such as parks and streets. Participants were informed about the study’s objectives and the voluntary nature of their participation, after which they signed informed consent forms approved by Shiraz University of Medical Sciences. If an individual declined to participate, a replacement with similar characteristics was recruited. The study questionnaire was divided into three sections: (1) demographic information, (2) questions designed to estimate the average social network size within the target population, and (3) questions aimed at estimating the frequency of individuals belonging to MARPs within the 18–45 age group.
NSU method
The NSU method is predicated on the assumption that the prevalence of target groups within a representative social network mirrors their prevalence in the overall population. In this method, respondents are queried about the high-risk behaviors of individuals within their social networks, which mitigates the risks of underreporting and underestimating sensitive and stigmatized behaviors. The primary formula utilized in the NSU method is as follows:
where
(m) represents the number of individuals known by the respondent in the target group (e.g., the number of PWID in the respondent’s network), (c) denotes the size of the respondent’s social network, (e) is the estimated size of the target group, and (t) is the total size of the target population. Accordingly, the average size of the individuals’ social network can be computed using the following formula:
where j and i represent the respondent’s group and the respondent, respectively.
In this study, the term “active social network” refers to the number of individuals with whom the respondent has direct contact, whether in person, via telephone, email, or any internet-based communication platforms (including social media such as Instagram, Snapchat, Discord, and similar services), and can connect with if necessary.
Calculation and refinement of C
To estimate the average social network size, we selected 16 first names (eight masculine and eight feminine) with known prevalence rates ranging from 0.1% to 4%. Names that were gender-neutral or combined were excluded. The frequency of the selected names was obtained from the National Organization for Civil Registration. We aggregated the data from all 16 names to calculate the value of (C), which was subsequently used to back-calculate the sizes of reference groups. The ratio of the actual size to the back-calculated size was computed for each name. Any ratio outside the range of 0.5–2 was deemed ineligible. During each iteration of the (C) calculation, the reference group with the least favorable ratio was excluded, and (C) was recalculated using the remaining reference groups. This process continued until all ratios fell within the acceptable range. Ultimately, the sizes of the MARPs were determined based on the final (C) value.
Crude estimation of size of MARPs
Using the final (C) value, the size of the target group was computed through the NSU method as follows:
We defined MARPs as individuals who engaged in the specified high-risk behaviors at least once in the past year. For example, an individual who injected drugs even once during that period was classified as a member of the PWID group.
Generalized NSU
Two primary sources of error affect the accuracy of the NSU method: transmission error and relative network size bias. The NSU method assumes that respondents possess complete awareness of the behaviors of individuals in their social networks. Transmission error occurs when respondents are unaware of the sensitive or stigmatized behaviors of individuals within their social network, leading to underreporting of the target population. To address this, a visibility factor is applied to correct for the probability that the behavior is not visible or known to the respondent.
The NSU method assumes that respondents possess complete awareness of the behaviors of individuals in their social networks. Relative network size bias arises because members of high-risk populations often have smaller social networks compared with the general population, which can distort size estimates. This bias is corrected using a popularity factor, which adjusts for differences in social network sizes between target groups and the general population.
By incorporating these correction factors, the generalized NSU method enhances the accuracy of population size estimations for hidden and hard-to-reach groups, ensuring more reliable data for public health planning and intervention.
We used the correction factors from the study conducted by Maghsoudi et al., in which adjustment coefficients for the two main sources of error in the NSU method—transmission error and relative network size bias—were calculated specifically for estimating the size of MARPs for HIV infection. 16
Finally, the size of the target group was computed using the following formula:
The 95% uncertainty intervals (UIs) of C and e were estimated by applying the bootstrap procedure using 1,000 iterations.
Results
From 768 participants, the sample was evenly split by sex (50% male and 50% female). The mean age was 30.11 years with a standard deviation of 8.57 years. Academic education was reported by 60% of males and 55% of females. Additional demographic details are provided in Table 1.
The Participants’ Demographic Characteristics
After computing the ratio of actual to back-calculated values, two reference groups were excluded from the final C estimates. The average active social network for individuals aged 18–45 years in Shiraz was found to be 203 contacts, comprising 103 males and 100 females. Accordingly, each individual aged 18–45 years in Shiraz knew an average of 203 people within the same age range (103 males and 100 females). The network sizes across different demographic strata are presented in Table 2.
Social Network Size Based on Participants’ Demographic Variables
SD, standard deviation.
Across the study, the total number of PWIDs, FSWs, and MSM known by participants was 976. The active social network consisted of 326 PWID, 584 FSW, and 66 MSM.
According to the most recent National Population and Housing Census of Iran, 866,160 individuals aged 18–45 years reside in Shiraz, including 429,555 males and 436,605 females. The crude and adjusted MARP sizes are shown in Tables 3 and 4, respectively. Tables 3 and 4 indicate that the highest MARP frequency was among FSW, while the lowest was among MSM.
Crude Size Estimation of MARPs
FSW, female sex worker; MARP, most-at-risk population; MSM, men who have sex with men; PWID, people who inject drug.
Adjusted Size Estimation of MARPs
CI, confidence interval; FSW, female sex worker; MARP, most-at-risk population; MSM, men who have sex with men; PWID, people who inject drug; UI, uncertainty interval.
Crude and adjusted MARP sizes by age groups are shown in Table 5.
Crude and Adjusted Size Estimation of MARPs Based on Age Groups
FSW, female sex worker; MARP, most-at-risk population; MSM, men who have sex with men; PWID, people who inject drug.
The results indicate that the highest frequencies of PWID, FSW, and MSM occur in the 35–39, 30–34, and 25–29 age groups, respectively. Conversely, the lowest frequency of PWID is in the 18–24 age group, while the lowest frequencies for FSW and MSM are in the 40–45 age group.
Discussion
Active social network
Based on comparisons across different methodologies, the present study employed the known population approach to estimate the size of active social networks. Sixteen first names were used as reference groups. Previous research has confirmed the high accuracy of this method for estimating social network size. 10
In our study, the average active social network size among individuals aged 18–45 years in Shiraz was 203 persons (103 males and 100 females). Reported estimates vary widely due to differences in target populations and methodological approaches. For example, Killworth reported a network size of 286 using the known population method. 18 In China, Gue et al. applied backward estimation across 19 reference groups and found an average of 310 for adults aged 18–60 years. 19 In Iran, the mean active social network size was estimated at 308 among individuals over 18 years of age. 20 Such discrepancies may arise from differences in sampling strategies, geographic coverage, and age groups under study. Whereas the cited studies were conducted nationally, our investigation focused on a single city with a narrower age range, which likely explains the smaller estimate. Shati et al. reported a C value of 259 in Tehran for individuals >18 years, 21 which is close to our estimate. This similarity may reflect the restriction of participants’ networks to a single city, while minor differences could be explained by age range, geographic, cultural, and demographic variations.
Size estimation of MARPs
People who inject drugs
Several studies have estimated the size of the population of PWID in Iran. In this study, the prevalence of PWID among 18–45-year-olds in Shiraz was 561 per 100,000 (0.5%). Shokoohi et al. reported a prevalence of 1.2% among the same age group in Kerman, 14 more than twice our estimate. This difference may be partly explained by gender composition, as most PWIDs in Iran are male, while our estimate included both sexes. Jafari et al. found a prevalence of 280 per 100,000 using the NSU method, 10 approximately half of our estimate. The higher estimate in our study may be related to age distribution, as injection drug use is less common in those <18 and >45 years. 22
Narouee et al. estimated a prevalence of 1,263 per 100,000 among adults in Iranshahr, Sistan, and Baluchistan province, 8 about twice the prevalence we observed. This region has socioeconomic vulnerabilities such as high unemployment and proximity to drug-trafficking routes, which contribute to higher rates of injection drug use. In Jahrom, Fars province, the NSU method yielded a prevalence of 0.26% among males aged 18–45 years. 17 The lower estimate may reflect the absence of correction for NSU-related biases, potentially leading to underestimation.
Female sex workers
Few studies have assessed the size of FSW in Iran. In the present study, the prevalence of FSW among women aged 18–45 years in Shiraz was 2,180 per 100,000 (about 2%). This indicates that nearly 2% of women in this age group reported at least one transactional sexual relationship in the past year. Sharifi et al. estimated the prevalence of FSW across 13 Iranian cities using the NSU method in 2014. Prevalence ranged from 0.14% in Tabriz to 2.44% in Isfahan, with 10 cities reporting rates between 1% and 2% and a national average of 1.54%. 5 In Shiraz, their study estimated a prevalence of 1.67% among women aged 15–49 years, 5 which is broadly consistent with our findings. Slight differences may reflect variations in age groups.
By contrast, Jafari et al. reported a prevalence of 709 per 100,000 among women over 18 years in Tabriz, 10 roughly one-third of our estimate. Similarly, Sharifi’s study showed the lowest national prevalence in Tabriz. 5 Such discrepancies likely reflect differences in age ranges, socioeconomic status, and sociocultural conditions between cities.
Men who have sex with men
Due to cultural sensitivities, stigma, and discrimination, few studies have assessed the size of MSM in Iran. Shokoohi estimated a prevalence of 0.5% among men aged 18–45 years in Kerman (95% UI: 0.2%–1%), 14 consistent with our estimate of 0.42%. This similarity may reflect comparable target populations and methodologies.
Jafari reported a prevalence of 247 per 100,000 among men over 18 years in Tabriz, 10 lower than our estimate of 420 per 100,000. Age distribution may account for this difference, as MSM are predominantly younger than 30. Regional sociocultural variations may also contribute. Hemayatkhah reported a prevalence of 0.34% (95% UI: 0.17%–0.85%) among men aged 18–45 years in Jahrom, 17 which is in line with our findings, again reflecting methodological and demographic similarities.
Implications for policy and practice
Our findings reveal that the estimated sizes of PWID, FSW, and MSM in Shiraz substantially exceed the current coverage of existing harm-reduction programs. In Shiraz, these services are primarily delivered through fixed and outreach facilities, including drop-in centers offering needle and syringe programs (NSPs), methadone maintenance therapy, and voluntary counseling and testing for HIV.23,24 While these programs exist, their coverage remains insufficient. Official reports suggest that current services reach only a fraction of the estimated MARPs.25,26 Previous research among PWID in Iran has documented that while awareness of NSPs reached 62.8%, actual utilization was reported by only 54.8% of this population. This awareness-to-uptake gap is even more pronounced for drug treatment services, with awareness at just 19.7% and utilization at a mere 9.1%. 26 Utilization is further hampered by profound stigma, fear of legal repercussions, and a lack of awareness about available services.1,23,24 This substantial gap between the estimated population size and the reach of current programs underscores the urgent need to expand services, implement mobile units, and integrate HIV services with other health care provisions to reduce stigma and improve access for these marginalized groups.
Limitations
This study has a limitation that should be considered when interpreting the findings. The potential transiency of MARPs may have influenced the accuracy of our population size estimates. Due to the stigmatized and often illegal nature of their behaviors, individuals in these groups frequently experience residential instability or migrate to larger cities to avoid legal consequences, seek employment opportunities, or access specialized services.27–30 Our estimates were based on a 5-year residency criterion, which may have captured a more stable subset of these populations while underestimating the total number of individuals who temporarily reside in or pass through Shiraz. Future studies should incorporate data on migration patterns and residential history to better understand the population dynamics of MARPs and to inform the design of targeted, flexible, and mobile health interventions.
Conclusion
This study applied the generalized NSU method to estimate the population sizes of the most-at-risk groups for HIV infection—namely FSW, PWID, and MSM—among adults aged 18–45 years in Shiraz, Iran. The findings revealed that the estimated sizes of these high-risk groups are substantially greater than the current coverage of harm-reduction programs. These results highlight the persistent gap in service coverage and underscore the urgent need to expand targeted HIV prevention and treatment initiatives. Implementing comprehensive strategies for increasing access to health care, reducing stigma, and strengthening surveillance measures is critical for controlling the HIV epidemic among hidden and marginalized populations. Policymakers and health authorities should utilize these robust estimates to inform resource allocation, develop more inclusive harm-reduction programs, and monitor progress toward national and regional HIV/AIDS control goals.
Ethical Approval
The study protocol was approved by the ethics committee of Shiraz University of Medical Sciences [Ethical code: (IR.SUMS.REC.1398.217)]. Before obtaining verbal informed consent, participants were informed about the anonymity of all collected data and the voluntary nature of their participation. Data were reported as aggregate data and did not identify subjects as individuals.
Data Availability
Data supporting the findings of this study are available from the corresponding author upon reasonable request. Access to restricted data may be provided under data-sharing agreements approved by the authors’ institution and relevant ethics boards.
Authors’ Contributions
All authors contributed significantly to the development of this research. J.H. conceptualized and designed the study and led the article writing. A.M. and M.R.B. assisted with data collection and data synthesis and contributed to writing and critical revisions of the article. All authors reviewed and approved the final version of the article and agreed to be accountable for all aspects of the work.
Footnotes
Acknowledgments
The authors would like to thank Ms. A. Keivanshekouh at the Research Improvement Center of Shiraz University of Medical Sciences for improving the use of English in the article.
Author Disclosure Statement
No competing financial interests exist.
Funding Information
This work was supported by the Research Vice-chancellor of Shiraz University of Medical Sciences under grant (97-01-04-17410).
