Skip to main content
Advertisement
  • Loading metrics

Seroprevalence of SARS-CoV-2 and risk factors in Bantul Regency in March-April 2021, Yogyakarta, Indonesia

  • Riris Andono Ahmad ,

    Roles Conceptualization, Data curation, Funding acquisition, Investigation, Methodology, Supervision, Validation, Writing – original draft, Writing – review & editing

    risandono_ahmad@ugm.ac.id

    Affiliations Faculty of Medicine, Public Health, Center for Tropical Medicine, Nursing Universitas Gadjah Mada, Yogyakarta, Indonesia, Faculty of Medicine, Public Health, Department of Biostatistics, Epidemiology and Population Health, Nursing Universitas Gadjah Mada, Yogyakarta, Indonesia

  • Citra Indriani,

    Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Project administration, Software, Validation, Writing – review & editing

    Affiliations Faculty of Medicine, Public Health, Center for Tropical Medicine, Nursing Universitas Gadjah Mada, Yogyakarta, Indonesia, Faculty of Medicine, Public Health, Department of Biostatistics, Epidemiology and Population Health, Nursing Universitas Gadjah Mada, Yogyakarta, Indonesia

  • Risalia Reni Arisanti,

    Roles Conceptualization, Data curation, Investigation, Methodology, Project administration, Validation, Writing – review & editing

    Affiliations Faculty of Medicine, Public Health, Center for Tropical Medicine, Nursing Universitas Gadjah Mada, Yogyakarta, Indonesia, Faculty of Medicine, Public Health, Department of Biostatistics, Epidemiology and Population Health, Nursing Universitas Gadjah Mada, Yogyakarta, Indonesia

  • Ratih Oktri Nanda,

    Roles Data curation, Formal analysis, Investigation, Project administration, Software, Writing – review & editing

    Affiliation Faculty of Medicine, Public Health, Center for Tropical Medicine, Nursing Universitas Gadjah Mada, Yogyakarta, Indonesia

  • Yodi Mahendradhata,

    Roles Conceptualization, Funding acquisition, Investigation, Methodology, Supervision, Validation, Writing – review & editing

    Affiliations Faculty of Medicine, Public Health, Center for Tropical Medicine, Nursing Universitas Gadjah Mada, Yogyakarta, Indonesia, Faculty of Medicine, Public Health, Department of Health Policy and Management, Nursing Universitas Gadjah Mada, Yogyakarta, Indonesia

  • Tri Wibawa

    Roles Conceptualization, Funding acquisition, Investigation, Methodology, Supervision, Validation, Writing – review & editing

    Affiliation Faculty of Medicine, Public Health, Department of Microbiology, Nursing Universitas Gadjah Mada, Yogyakarta, Indonesia

Abstract

COVID-19 case counts in Indonesia inevitably underestimate the true cumulative incidence of infection due to limited diagnostic test availability, barriers to testing accessibility and asymptomatic infections. Therefore, community-based serological data is essential for understanding the true prevalence of infections. This study aims to estimate the seroprevalence of SARS-CoV-2 infection and factors related to the seropositivity in Bantul Regency, Yogyakarta, Indonesia. A cross-sectional study involving 425 individuals in 40 clusters was conducted between March and April 2021. Participants were interviewed using an e-questionnaire developed in the Kobo toolbox to collect information on socio-demographic, COVID-19 suggestive symptoms, history of COVID-19 diagnosis and COVID-19 vaccination status. A venous blood sample was collected from each participant and tested for immunoglobulin G (Ig-G) SARS-CoV-2 antibody titers using the enzyme-linked immunosorbent assay (ELISA). Seroprevalence was 31.1% in the Bantul Regency: 34.2% in semi-urban and 29.9% in urban villages. Participants in the 55–64 age group demonstrated the highest seroprevalence (43.7%; p = 0.00), with a higher risk compared to the other age group (aOR = 3.79; 95% CI, 1.46–9.85, p<0.05). Seroprevalence in the unvaccinated participants was 29.9%. Family clusters accounted for 10.6% of the total seropositive cases. No significant difference was observed between seropositivity status, preventive actions, and mobility. Higher seroprevalence in semi-urban rather than urban areas indicates a gap in health services access. Surveillance improvement through testing, tracing, and treatment, particularly in areas with lower access to health services, and more robust implementation of health protocols are necessary.

Introduction

Coronavirus disease 2019 (COVID-19) is a respiratory illness caused by the newly discovered severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), leading to a global pandemic, including Indonesia. Since the first confirmed SARS-CoV-2 on March 2 2020 [1], Indonesia has experienced exponential growth of COVID-19 cases in 34 provinces. As of May 2 2021, the government of Indonesia reported 1.677.274 confirmed COVID-19 cases and 45.796 deaths [2] and became one of the countries with high cumulative and incidence cases of COVID-19. The number of COVID-19 cases in Indonesia has been increasing rapidly since January, leading to the highest peak thus far, with provinces in Java making up 64% of the national cases. Two notable waves were observed from March 2020 to December 2021 (Fig 1). Based on the daily trend of COVID-19 cases, the first wave of COVID-19 cases was observed in January, with the highest number of reported cases reaching 14518 cases on January 16. After that, the cases fluctuated and reached the highest during the third wave in September 2022, with 56.757 recorded cases.

thumbnail
Fig 1. COVID-19 cases daily trend in Indonesia March 2020 –December 2021 (data source: https://ourworldindata.org/coronavirus/country/indonesia, accessed on October 4 2022).

https://doi.org/10.1371/journal.pgph.0000698.g001

However, these case counts inevitably underestimate the true cumulative incidence of infection because of limited diagnostic test availability, barriers to testing accessibility, and asymptomatic infections, not to mention underreporting [2]. As a consequence, the national prevalence of SARS-CoV-2 remains unknown.

The lack of clarity about the number of SARS-CoV-2 infections across Indonesia limits the Indonesian government’s ability to plan appropriately, prepare and respond to this epidemic. Monitoring the incidence of newly diagnosed cases of severe COVID-19 and the case fatality rate is critical to address the demands on the healthcare system.

One of the epidemiological investigations used to determine the level of disease spread is to conduct a seroprevalence survey. According to the Centers for Disease Control and Prevention (CDC), this survey uses serological tests to detect antibodies in the blood, indicating an infection [3]. This test uses an enzyme-linked immunosorbent (ELISA) where the antigen used is purified SARS-CoV-2 S protein (without live virus). Population-based serological testing provides better estimates of the cumulative incidence of infection by complementing diagnostic testing of acute illness and helping to inform the public health response to COVID-19. A seroprevalence study could also be a powerful tool to detect subclinical infections and improve policy-making in the country [4,5]. Furthermore, as the world moves through the vaccine and variant era, synthesizing seroepidemiology findings is increasingly important to track the spread of infection, identify disproportionately affected groups, and measure progress towards herd immunity [2].

Seroprevalence varies geographically; the denser urban areas have higher seropositivity rates than rural areas [6]. A study in East Java, Indonesia, in the second semester of 2020 showed a higher prevalence in Surabaya (13.1%), which is an urban area, than in Jombang (9,9%), a rural area [7]. The epidemiological trend also implicates SARS-CoV-2 spread among rural communities only later in the epidemic wave [8,9], which would require sound anticipatory interventions. A seroprevalence survey involving more diverse groups of people among urban and rural communities is necessary to grasp the overall picture of SARS-CoV-2 infection. Bantul Regency has become one area that contributes to many cases and leads to the high transmission of COVID-19 in the Yogyakarta Provinces. Therefore, this study aims to estimate the SARS-CoV-2 seroprevalence, seropositive risk factors, and COVID-19 vaccine acceptance in Bantul Regency, Yogyakarta.

Materials and methods

Study setting

Bantul Regency is located in the southern region of Yogyakarta Province, covering 506.85 km2 (Fig 2). The regency consists of 17 sub-districts and 75 villages, 30% semi-urban. Most adults are engaged in the non-formal sector, such as farming, trade and the service industry [10]. The daily mobility of residents between districts to and from Bantul is high, which may increase the risk of disease transmission [11].

thumbnail
Fig 2. Map of Bantul Regency, Special Region of Yogyakarta, Indonesia.

https://doi.org/10.1371/journal.pgph.0000698.g002

Study design and sampling

A cross-sectional study was conducted from March to April 2021. The study population consisted of individuals who lived in the Bantul regency for at least six months, with 1,018,402 inhabitants [10]. EpiInfo was used to calculate the sample size.

Sampling was determined using a multistage cluster random sample adopted from the WHO/EPI rapid survey [12]. The calculated sample size was 414, which accounted for the estimated prevalence of 20%, 6% sampling error, a significance level of 0.05 with a design effect of 2 and a non-response rate of 15%. There were 40 clusters, with 11 participants per cluster. The cluster was based on the smallest administrative area called RT (rukun tetangga or household group). The RT consists of approximately 50–70 households.

We used systematic random sampling to select households at the selected cluster. In the selected households, we recruited all eligible household members, i.e. anyone who lived under the same building and resided in the study area for at least six months, with a minimum age of 5 years old.

Participants must meet the following inclusion criteria: residing in the study location at least six months before the survey commences, age of 5 years old or older, able to communicate verbally, give written consent to participate in the research or consent from parents/ guardians for respondents under 18 years old. Participants aged 13–17 were interviewed by their parents/guardians. However, for participants under 13 years old, interviews were conducted with the parents. Meanwhile, the exclusion criteria include people with existing chronic illnesses (including immunocompromised individuals with a history of blood disorders and people with mental disorders).

Data collection

A door-to-door visit was conducted to collect the primary data. Written informed consent was obtained from each study participant before data collection. Fifteen enumerators, including a phlebotomist, were involved in the data collection. Two supervisors were assigned to ensure the methodology and conduct the spot check. All study teams were involved in four training days regarding methods, data entry, phlebotomy and ethics in research.

Using an aseptic procedure, two millilitres of blood were drawn from a cubital vein with a disposable sterile syringe. The specimen was then kept in an EDTA tube. Before transportation, blood samples were placed in a cool box with an ice pack and transported within the same day to the laboratory. Blood samples were further stored in the refrigerator. Each participant’s data included a unique identifier (barcode label) linked to their blood sample and data for tracking and confidentiality. The blood sample was transported, examined and stored in the Laboratory of Microbiology, Faculty of Medicine, Public Health, and Nursing UGM. Plasma was then tested for Ig-G anti-SARS-Cov2 procedures using (Human Anti-2019 n-CoV(N) IgG ELISA Kit V1.5 FineTest [13].

Risk factors information at household and individual levels were obtained using an electronic structured questionnaire developed in the KoBo toolbox [14,15]. Variables included in the household questionnaire were socio-demographic information (age, gender, relationship with the head of household, and household income). Individual questionnaires were used to gather socio-demographic information (age, gender, highest education, occupation), COVID-19 vaccination status, previous diagnosis of COVID-19 and symptoms related to COVID-19 within the last six months, preventive actions taken, and mobility in the previous two weeks. The data manager validated data daily. A day after the data entry, the data manager sent the supervisor feedback and confirmation regarding completeness and data consistency.

Updated cumulative data on notified cases were obtained through the COVID-19 surveillance system conducted by Bantul District Health Office Bantul from March to April 2021. COVID-19 was defined as a symptomatic or asymptomatic person with a positive PCR result tested for SARS-CoV-2. Population numbers per sub-district were obtained through Bantul Statistical Bureau for calculating incidence per 1000 population [10]. In addition, data on urban-rural classification was obtained from the National Statistical Bureau [16].

Statistical analysis

Proportion and percentage describe the seroprevalence and socio-demographic characteristics of study participants. A map developed in ArcGIS spatially compared seroprevalence and incidence of notified cases by district and urban-rural status [17,18]. Bivariate analysis was conducted to identify the association between presumed risk factors (age, gender, occupation, comorbidities, prevention taken and mobility over the last two weeks) and anti-SARS CoV-2 seropositivity. After considering collinearity, a multiple logistic regression was used to evaluate risks between seropositive and seronegative groups with adjustment for sex. The P-value would be considered statistically significant at p < 0.05. All statistical analysis was performed using STATA 14.0. Raw data used for analysis is provided in (S1 Data).

Ethics

The Medical and Health Research Ethics Committee of the Faculty of Medicine, Public Health and Nursing Universitas Gadjah Mada approved the study (Ref: KE/1242/12/2020). Written informed consent was obtained from adult respondents and parents of enrolled children. Confidentiality of information from the respondents was upheld with utmost care throughout data collection, processing and analysis for all data collected. Therefore, their names were included in the notes only for traceability and referral during the data analysis.

Results

Characteristics of study participants

A total of 425 people participated and were tested during the survey. See 1 File. The number of females (59.1%) is higher than that of males (40.9%); 47.2% of the participants were unemployed/students/housewives, and 72.5% lived in urban areas. The majority of the participants had no comorbidities and were not vaccinated.

The prevalence of SARS-CoV-2 seropositivity among the participants in this study was 31.1% (n = 132/425). A significant difference was observed in the seroprevalence among age groups (p = 0.000), with the highest proportion reported in the 55–64 age group (43.7%; n = 31/71). Meanwhile, the under-15 years age group showed no seropositivity.

A significant difference was also observed among occupation groups (p = 0.009). The highest seroprevalence was demonstrated by participants working as daily workers/farmers (37.2%), followed by professional/health workers (34.6%), and unemployed/students/housewives (26%). Seroprevalence did not differ between semi-urban and urban areas, even though we observed that semi-urban areas had higher seroprevalence (34.2%; n = 40/117) than urban areas (29.9%; n = 92/308). Females demonstrated higher seroprevalence; nevertheless, no association was found between genders and seropositivity. Seropositivity also does not differ in occupation category (Table 1).

thumbnail
Table 1. Distribution of study participants by background characteristics.

https://doi.org/10.1371/journal.pgph.0000698.t001

This study also explored seroprevalence based on the presence of selected chronic diseases and the history of comorbidities. Participants with a history of chronic conditions accounted for 22.6% of the total participants. Among those with comorbidities, seroprevalence for participants with obesity, diabetes mellitus, and hypertension was 50%, 40% and 35.4%, respectively.

Participants with a prior diagnosis of COVID-19 in the last six months accounted for 1 out of 132 seropositive cases. Seroprevalence among participants with no previous experience of COVID-19-related symptoms was found to be 31.1%. Seropositivity was observed in 14 participants from seven households, making family clusters account for 10.6% of the total positive cases.

Most of the participants in this study reported decreased social interactions and mobility during this period. Study participants stated that they attended fewer traditional and religious activities (90.8%), had reduced visits to relatives and friends (90.1%, decreased visits to markets (60.9%), and spent more time at home, except for essential tasks ((73.7%). If participants did leave home, they reported wearing masks (89.4%), regular washing of hands for 20 seconds with running water (80.9%) and maintaining physical distancing in public areas (64.5%). However, no significant differences were observed between individual preventive actions and the level of mobility with seropositivity status.

Seroprevalence of anti-SARS-CoV-2 based on the geographical distribution

The highest seroprevalence was observed in 3 semi-urban areas (Pajangan, Dlingo, Sanden). Most semi-urban areas have relatively higher seroprevalence than urban areas. However, we found no statistical significance (Fig 3). The highest distribution of confirmed cases acquired from routine regional data reported higher cumulative cases in urban areas, such as Banguntapan, Bantul, Sewon, and Jetis.

thumbnail
Fig 3. Seroprevalence and reported cumulative incidence of SARS-Cov-2 in April 2021 Based on District.

https://doi.org/10.1371/journal.pgph.0000698.g003

Risk factors associated with SARS-CoV-2 seropositivity

This study further explored the risk factors associated with the seropositivity of SARS-CoV-2, adjusting for age, sex, occupation, comorbidities, and vaccination status. A significant association was observed within specific age groups. The odds of SARS-CoV-2 seropositivity are higher in the age 55–64 (adjusted odds ratio [aOR] = 3.79; 95% CI 1.46–9.85, p = 0.006).

Females demonstrated higher seroprevalence. Nevertheless, no association was found between genders and seropositivity. Seropositivity also does not differ in occupation category (aOR = 0.81; 95% CI 0.50–1.31, p>0.005) (Table 2).

thumbnail
Table 2. Factors associated with SARS-Cov-2 seropositivity in Bantul Regency (N = 425).

https://doi.org/10.1371/journal.pgph.0000698.t002

Comorbidities were presumed as one of the risk factors of SARS-CoV-2 seropositivity. The odds of SARS-CoV-2 seropositivity were higher in participants with obesity (aOR = 2.46; 95% CI 0.57–10.65) and diabetes mellitus (aOR = 1.12; 95% CI 0.38–3.31) than participants without each respective comorbidity. However, this finding is not statistically significant (p-value>0.05).

Although vaccination had not been made available to the general public during the data collection period, most participants (91.5%) were aware that the government would provide vaccinations in the future and stated their willingness (79.1%) to be vaccinated. However, a small proportion (19.1%) said they did not want to be vaccinated (Table 3). Reasons for this included fear of adverse effects, concerns about the safety and effectiveness of the vaccines, and religious beliefs.

thumbnail
Table 3. Vaccine acceptance of the participants (N = 425).

https://doi.org/10.1371/journal.pgph.0000698.t003

Discussion

This SARS-CoV-2 seroprevalence study found that the prevalence of IgG antibodies against SARS-CoV-2 was 31.1% in Bantul, significantly higher than the 11.4% seropositive which was found in a study in East Java in late 2020 [7]. On the other hand, surveillance data of Bantul Regency DHO reported a considerably lower cumulative incidence of 1.1% during the study. This finding is in line with previous research suggesting that the estimates obtained from seroprevalence were 18.1 times higher than the corresponding cumulative incidence of COVID-19 infections, implying that confirmed cases are a poor indicator of the extent of the disease spread [19]. While antigen testing and Polymerase Chain Reaction (PCR) to monitor COVID-19 detect the presence of specific antigens or genetic material indicating current viral infection, serological tests for COVID-19 are used to detect antibodies that the body produces in response to the virus. Serological tests can detect antibodies in people who have been infected with COVID-19 in the past, even if they were asymptomatic or had mild symptoms. Those can help identify people who may have been infected but were never diagnosed. In population-level studies, serological testing is used to estimate the proportion of people who have been infected with COVID-19 in a given area.

Seroprevalence varies geographically, and previous research indicated that urban areas mainly reported a higher number of seropositivity [6]. The reported COVID-19 cases in the Bantul regency also showed a similar trend, in which more cases were reported in urban areas (Bantul, Banguntapan, Sewon, and Jetis). However, our study showed that the three highest seropositivities were observed in semi-urban areas (Pajangan, Dlingo and Sanden). Meanwhile, areas with the highest reported COVID-19 cases showed lower seropositivity. Thus, this study suggest that the higher transmission observed in semi-urban areas was due to lower access to health services, including the lack individual testing, resulting in unreported positive cases. These discrepancies can be explained by the limited and centralized diagnosis facilities when the study was conducted.

Urban areas with dense populations were impacted earlier at the pandemic’s beginning. A previous study on epidemiological modelling predicted that COVID-19 spread faster in urban than rural areas [20]. Despite the rapid transmission, urban areas are usually supported by better surveillance management, such as more accessible testing, tracing, and infection control. Consequently, cases in urban areas were reported faster and earlier.

Meanwhile, rural areas usually have shortages in resources and workforces, which affects their ability to detect, respond, prevent, and control infectious disease outbreaks [21,22]. In this study, a large proportion of seropositivity was reported in participants with no prior diagnosis of COVID-19, indicating that the infection might be undiagnosed due to a lack of individual testing. Thus, this study suggests that the higher transmission observed in semi-urban areas was lower access to health services, including lack of individual testing, resulting in unreported positive cases.

Higher seroprevalence in this study was observed in females. However, no statistically significant gender difference was reported. A previous study also suggested that SARS-CoV-2 IgG/IgM dynamic is mainly affected by age and disease severity, not sex [23]. Higher seropositivity was observed in the age group 55–64, and this finding supports many previous studies that indicate older age as one of the risk factors for COVID-19.

This research found a portion of seropositivity among people who primarily stayed at home. Additionally, occupations requiring less mobility, such as housewives, students, and the unemployed, also accounted for seropositivity. The research also discovered a tenth of family clusters, suggesting that the transmission may have occurred at the household level. Previous research contends that the primary mode of COVID-19 infection is through the household spread [24,25].

Our study found no difference between high and low mobility, suggesting they were already widespread community transmission, even in rural areas. Preventive measures did appear to prevent acquiring an infection, which might be because it is challenging to do preventive measures in a household setting or to practice preventive measures consistently outside [26].

Surveillance and containment measures such as large-scale social restrictions and other anticipatory prevention tools become a priority to curb transmission, primarily focusing on areas with lower access to health services. Earlier in the pandemic, outbreak containment was concentrated in the dense urban areas, which may have reduced the number of cases. However, due to high mobility, the transmission slowly moved to semi-urban and rural communities. People living in rural communities might develop a false sense of security and take fewer precautions at the beginning than the urban communities [8,9].

One of the most effective ways to prevent fatalities caused by COVID-19 is through vaccination. Due to the lack of a vaccination program for general populations during the data collection period, it is reasonable that higher seroprevalence was found in unvaccinated groups. As this study explored the seroprevalence when most people were still unvaccinated, a better picture of the infection spread in Bantul was obtained since there was no implication from vaccine-induced antibodies.

Despite the findings, this study has the following limitations. Following a natural infection, antibody titers peak and begin to wane in various manners, with some in shorter duration [27]. However, this study did not consider the assay performance concerning the waning immunity. Thus it may underestimate the true prevalence. Furthermore, the analysis results only showed the association, not causality. Finally, as a cross-sectional analysis, this study only analyzed the variables at one point of time and did not explore the seroprevalence changes over a longer period. Therefore, future research is needed to conduct a periodic or longitudinal survey to determine the prevalence in the longer term.

Conclusion

This serosurvey demonstrated a higher seroprevalence than reported data in the same period. Based on the findings, it is strongly recommended that the local government strengthen the surveillance and 3T (testing, tracing, and treatment) efforts by involving the task force at the neighbourhood community and village levels throughout Bantul Regency, particularly in areas with lower access to health services. Besides that, it needs to increase awareness and implementation of health protocols for high-mobility individuals to prevent transmission within the household. As vaccination program for the general population are being rolled out, it is crucial to provide adequate implementation information, including health resources and logistics support. This study can be implemented in other areas, both at the district/city and provincial levels, to understand the seroprevalence of SARS-CoV-2 in Indonesia better. This study provides a district-level view of the extent of COVID-19 spread and a different approach to conducting serosurvey among diverse populations in various regions to fit the gaps in understanding COVID-19’s spread globally.

Supporting information

S1 Table. The response rate of each district.

https://doi.org/10.1371/journal.pgph.0000698.s001

(PDF)

S1 Data. Data of Seroprevalence of SARS-CoV-2 in Yogyakarta.

https://doi.org/10.1371/journal.pgph.0000698.s002

(XLSX)

Acknowledgments

We thank all study participants involved in the study, Bantul DHO, for supporting the research and all of the enumerators and the supervisor for the great work during the data collection.

References

  1. 1. World Health Organization. Coronavirus disease 2019 (COVID-19): situation report, 1 [Internet]. World Health Organization; 2020 Mar [cited 2023 April 10]. Available from: https://cdn.who.int/media/docs/default-source/searo/indonesia/covid19/who-indonesia-situation-report-1.pdf?sfvrsn=6be5b359_0.
  2. 2. World Health Organization. Coronavirus disease 2019 (COVID-19): situation report, 51 [Internet]. World Health Organization; 2020 Mar [cited 2021 December 23]. Available from: https://apps.who.int/iris/handle/10665/331475.
  3. 3. CDC. Cases, Data, and Surveillance [Internet]. Centers for Disease Control and Prevention. 2020 [cited 2021 December 23]. Available from: https://www.cdc.gov/coronavirus/2019-ncov/covid-data/serology-surveillance/index.html.
  4. 4. West R, Kobokovich A, Connell N, Gronvall GK. COVID-19 Antibody Tests: A Valuable Public Health Tool with Limited Relevance to Individuals. Trends Microbiol. 2021 Mar;29(3):214–23. pmid:33234439.
  5. 5. Ojeda DS, Ledesma MMGL, Pallarés HM, Navarro GSC, Sanchez L, Perazzi B, et al. Emergency response for evaluating SARS-CoV-2 immune status, seroprevalence and convalescent plasma in Argentina. PLOS Pathogens. 2021 Jan 14;17(1):e1009161. pmid:33444413.
  6. 6. Mohanan M, Malani A, Krishnan K, Acharya A. Prevalence of SARS-CoV-2 in Karnataka, India. JAMA. 2021 Mar 9;325(10):1001–3. pmid:33538774.
  7. 7. Megasari NLA, Utsumi T, Yamani LN, Juniastuti , Gunawan E, Furukawa K, et al. Seroepidemiological study of SARS-CoV-2 infection in East Java, Indonesia. PLOS ONE. 2021 May 6;16(5):e0251234. pmid:33956869.
  8. 8. Peters DJ. Community Susceptibility and Resiliency to COVID-19 Across the Rural-Urban Continuum in the United States. The Journal of Rural Health. 2020;36(3):446–56. pmid:32543751.
  9. 9. Murhekar MV, Bhatnagar T, Selvaraju S, Saravanakumar V, Thangaraj JWV, Shah N, et al. SARS-CoV-2 antibody seroprevalence in India, August–September 2020: findings from the second nationwide household serosurvey. The Lancet Global Health. 2021 Mar 1;9(3):e257–66. pmid:33515512.
  10. 10. BPS Bantul. Bantul Regency in Figure 2021. Badan Pusat Statistik Kabupaten Bantul; 2021.
  11. 11. Sadali MI, Sadali MI. Mobilitas Pekerja Pada Kawasan Industri Piyungan di Kabupaten Bantul, Daerah Istimewa Yogyakarta [Internet]. Universitas Gadjah Mada; 2014 [cited 2021 Dec 23]. Available from: http://etd.repository.ugm.ac.id/penelitian/detail/75081.
  12. 12. WHO. The EPI Coverage Survey: Module: Training for mid-level managers (MLM). World Health Organization; 2008.
  13. 13. Van Honacker E, Coorevits L, Boelens J, Verhasselt B, Van Braeckel E, Bauters F, et al. Sensitivity and specificity of 14 SARS-CoV-2 serological assays and their diagnostic potential in RT-PCR negative COVID-19 infections. Acta Clin Belg.: 1–6. pmid:33350362.
  14. 14. Dizon TJR, Saito N, Reñosa MDC, Bravo TAP, Silvestre CJ, Endoma VF, et al.: Experiences in Using KoBo Collect and KoBo Toolbox in a Cross-Sectional Dog Population and Rabies Knowledge and Practices Household Survey in the Philippines. Stud Health Technol Inform 2022, 290:1082–1083. pmid:35673216
  15. 15. The KoboToolbox Software [https://www.kobotoolbox.org/about-us/software].
  16. 16. BPS. Klasifikasi Perkotaan dan Perdesaan di Indonesia: Peraturan Kepala Badan Pusat Statistik Nomor 37 Tahun 2010. Badan Pusat Statistik; 2010.
  17. 17. Zamani M, Poustchi H, Mohammadi Z, Dalvand S, Sharafkhah M, Motevalian SA, et al.: Seroprevalence of SARS-CoV-2 antibody among urban Iranian population: findings from the second large population-based cross-sectional study. BMC Public Health 2022, 22(1):1031. pmid:35606743
  18. 18. ArcGIS Overview [https://www.esri.com/en-us/arcgis/about-arcgis/overview].
  19. 19. Bobrovitz N, Arora RK, Cao C, Boucher E, Liu M, Donnici C, et al. Global seroprevalence of SARS-CoV-2 antibodies: A systematic review and meta-analysis. PLOS ONE. 2021 Jun 23;16(6):e0252617. pmid:34161316.
  20. 20. Stier AJ, Berman MG, Bettencourt LMA. COVID-19 attack rate increases with city size [Internet]. 2020 Apr [cited 2021 December 23] p. 2020.03.22.20041004. Available from: https://www.medrxiv.org/content/10.1101/2020.03.22.20041004v2.
  21. 21. Turner-Musa J, Ajayi O, Kemp L. Examining Social Determinants of Health, Stigma, and COVID-19 Disparities. Healthcare (Basel). 2020 Jun 12;8(2):168. pmid:32545647.
  22. 22. Patel L, Elliott A, Storlie E, Kethireddy R, Goodman K, Dickey W. Ethical and Legal Challenges During the COVID-19 Pandemic: Are We Thinking About Rural Hospitals? The Journal of Rural Health. 2021;37(1):175–8. pmid:32282953.
  23. 23. Luo C, Liu M, Li Q, Zheng X, Ai W, Gong F, et al. Dynamic changes and prevalence of SARS-CoV-2 IgG/IgM antibodies: Analysis of multiple factors. International Journal of Infectious Diseases. 2021 July 1;108:57–62. pmid:33932603.
  24. 24. Shen M, Peng Z, Guo Y, Rong L, Li Y, Xiao Y, et al. Assessing the effects of metropolitan-wide quarantine on the spread of COVID-19 in public space and households. International Journal of Infectious Diseases. 2020 Jul 1;96:503–5. pmid:32416146
  25. 25. Madewell ZJ, Yang Y, Longini IM Jr, Halloran ME, Dean NE. Household Transmission of SARS-CoV-2: A Systematic Review and Meta-analysis. JAMA Network Open. 2020 Dec 14;3(12):e2031756 PMCID: PMC7737089. pmid:33315116
  26. 26. Peak CM, Kahn R, Grad YH, Childs LM, Li R, Lipsitch M, et al. Individual quarantine versus active monitoring of contacts for the mitigation of COVID-19: a modelling study. The Lancet Infectious Diseases. 2020 Sep 1;20(9):1025–33 pmid:32445710.
  27. 27. Seow J, Graham C, Merrick B, Acors S, Pickering S, Steel KJA, et al. Longitudinal observation and decline of neutralizing antibody responses in the three months following SARS-CoV-2 infection in humans. Nat Microbiol. 2020 Dec;5(12):1598–1607. pmid:33106674.