If you don't remember your password, you can reset it by entering your email address and clicking the Reset Password button. You will then receive an email that contains a secure link for resetting your password
If the address matches a valid account an email will be sent to __email__ with instructions for resetting your password
COVID-19 epidemic monitoring after non-pharmaceutical interventions: The use of time-varying reproduction number in a country with a large migrant population
Rt can be used to closely monitor the non-pharmaceutical interventions (NPI) of the COVID-19 epidemic.
•
Activity of the epidemic in Oman is examined.
•
Other factors, like the return of overseas students, have increased the epidemic activity.
•
Responses to NPI are different between migrants and natives.
Abstract
Background
COVID-19’s emergence carries with it many uncertainties and challenges, including strategies to manage the epidemic. Oman has implemented non-pharmaceutical interventions (NPIs) to mitigate the impact of COVID-19. However, responses to NPIs may be different across different populations within a country with a large number of migrants, such as Oman. This study investigated the different responses to NPIs, and assessed the use of the time-varying reproduction number (Rt) to monitor them.
Methods
Polymerase chain reaction (PCR) laboratory-confirmed COVID-19 data for Oman, from February 24 to June 3, 2020, were used alongside demographic and epidemiological information. Data were arranged into pairs of infector–infectee, and two main libraries of R software were used to estimate reproductive number (Rt). Rt was calculated for both Omanis and non-Omanis.
Findings
A total of 13,538 cases were included, 44.9% of which were Omanis. Among all these cases we identified 2769 infector–infectee pairs for calculating Rt. There was a sharp drop in Rt from 3.7 (95% confidence interval [CI] 2.8–4.6) in mid-March to 1.4 (95% CI 1.2–1.7) in late March in response to NPIs. Rt then decreased further to 1.2 (95% CI 1.1–1.3) in late April after which it rose, corresponding to the easing of NPIs. Comparing the two groups, the response to major public health controls was more evident in Omanis in reducing Rt to 1.09 (95% CI 0.84–1.3) by the end of March.
Interpretation
Use of real-time estimation of Rt allowed us to follow the effects of NPIs. The migrant population responded differently than the Omani population.
The emergence of a new infectious disease carries with it many challenges and uncertainties regarding its natural history, clinical course, transmissibility, and methods of control. With these in mind, the leading resort for control when there is no treatment or vaccine is to apply well-known and strict public-health measures to mitigate and prevent its spread (
). These measures were implemented by the Chinese government after the reporting of 27 cases of corona virus disease 2019 (COVID-19) in Wuhan, in China’s Hubei province, on 30 December 2019. Measures included the isolation of confirmed and suspected cases, and restriction of movement inside the country (
). During the two months following the appearance of COVID-19, and despite the measures that were taken, international travel spread cases to 26 countries around the world (
With the increasing number of cases and the introduction of multiple NPIs by most countries, there is a need to assess and monitor the transmissibility of the disease as a measurement of the effectiveness of control measures. The effective reproduction number, Rt, defined as the average number of secondary cases from a partially susceptible population per infectious case (
). The effects of the control and the dynamics of transmission of COVID-19 were expected to be different between these two populations—Omanis and non-Omanis. Similar situations that have led to increasing numbers of cases have been documented in many neighbouring countries of the Gulf Cooperation Council (GCC) (
The government of Oman responded to the COVID-19 pandemic much like most countries by implementing multiple NPIs in phases to control the disease, beginning in mid-March 2020. Examples ranged from restricting flights from infected countries, to closing schools and commercial activities (Table 1).
Table 1Timeline for preparation phases and non-pharmaceutical intervention taken in the sultanate to prevent and contain the COVID-19 outbreak.
Time frame
Non-pharmaceutical interventions
Phase 1: Preparatory phase
Mid-February to mid-March
Recommendation to avoid travelling to China unless absolutely necessary
Suspension of flights to and from China
Quarantine on arrivals from China
Quarantine application on arrivals from Singapore, South Korea, Japan, and Iran (February 2)
Suspension of flights from Iran
Avoidance of travel to the above countries
Suspension of flights to and from Italy
Quarantine on arrivals from Italy
Quarantine on arrivals from Egypt
Formation of the supreme committee (March 10)
Phase 2: Major public health control phase
Starting from mid-March to mid-April
Closure of schools and higher education institutions
Closure of port of entries (airports, seaports, and ground crossings)
Suspension of all public transport
Closure of mosques and all non-Muslim places of worship
Closure of all movie theatres
Prohibition of public gatherings and conferences
Closure of all stores in shopping malls, with the exception of food and consumer catering shops, clinics, pharmacies, and optical stores
Closure of traditional markets, such as Mutrah, Nizwa, Rustaq, and Sinaw, and closure of popular markets, such as Wednesday markets, Thursday markets, and Friday markets
Prohibition of the provision of food in restaurants and cafes, including those in hotels, with the exception of external requests
Suspension of sports activities of all kinds, and closure of sports and cultural clubs
Closure of gyms and health clubs, and men’s and women’s barbers/hairdressers and beauty shops
Reduction of the number of employees in workplaces in government agencies to no more than 30% of the total number of employees
Closure of all exchange shops, provided that the banks provide exchange services,
Return of all overseas students (March 23)
Activation of the role of the military and security services in reducing the movement of residents and citizens
Isolation of the Governorate of Muscat until April 1
Closure of Muscat and Mutrah (April 10)
Closure of Jaalan Bani Bu Ali (April 16)
Phase 3: Gradual easing up
From mid-April to end of May
Expanded opening of some activities in grades (April 22 and May 18)
Allowing 50% of governmental employees to return to work (May 31)
In this study, we analysed the dynamics of COVID-19 infection transmissibility in Oman in the different populations (Omani and non-Omani), and the effects of the introduction of non-pharmaceutical measures on disease transmissibility.
Methods
Data source: national COVID-19 data
Data for laboratory-confirmed COVID-19 cases in Oman, from February 24 to June 3, 2020, were included in this study. Suspected cases of COVID-19 (as per the Ministry of Health case definition) (
) were confirmed by PCR testing. Subsequently, all confirmed cases underwent epidemiological investigation, generating the daily final line list, which included information on each patient’s demography (age, gender, residency, nationality), epidemiology (source of infection if known, date of onset, primary case [infector] or secondary case [infectee], designation [cluster or sporadic]). The type of surveillance, whether active (proactive case finding) or passive (regular reporting by health-care institutions), was also included in the data, with only the passive surveillance cases used in our study.
Time-varying reproduction number
The basic reproductive number before mitigation starts is called R0. The reproduction number after mitigation starts, Rt, is a measure of the transmissibility of the infection, and is defined as the average number of new infections one case can produce (
). An Rt of more than 1 means that the infection is spreading, with more cases generated, whereas an Rt of less than 1 means that the spread of infection is decreasing. Theoretically, we need information about the generation time—defined as the time period between the infection of the index and the next case. However, this information is usually difficult to ascertain, and therefore information regarding the serial interval (defined as the interval between disease onset in the index and the next case) distribution in the data is used instead.
Using R software, we utilised two main libraries to estimate Rt—Epicontact and EpiEstim (
) is to arrange the data and help in estimating the distribution parameters of the serial interval for our data. The estimate-R function in the EpiEstim library (
) assumes a gamma distribution of the SI and models the transmission of the infection using a Poisson likelihood to calculate the instantaneous reproduction number.
We arranged our line list into two parts, as required by Epicontact library (
)—the main line list data and the contact data. The daily list of all new PCR-positive cases contains all the demographic and exposure data for each case, whereas the contact data contain information about the transmission of infection between each identifier number. As case-by-case transmission data are available for many cases on the line list, we used the Epicontact R library to find out the serial interval for our data and its distribution for the entire population, and classified these data by nationality.
We used the serial interval distribution calculated by the Epicontact library to estimate the Rt for the entire population, classified by nationality, for the period between February 24 and June 3, 2020. We used the estimate_R function in the EpiEstim library to calculate Rt, given the distribution of the serial interval (obtained by the Epicontacts library) and incidence time series. This was done through a sliding window of 7 days using the Poisson transmission model (
The comparison of Rt values between Omanis and non-Omanis was performed after extracting Rt values from each estimate_R object and plotting the two time series against each other.
To investigate the differences in behaviour of transmission for the two groups, the epidemic curve trend was investigated according to cluster/sporadic types for Omanis versus non-Omanis using the geom_smooth function in a ggplot2 library (
). The comparison of Rt values for the different groups was performed using a simple line plot.
Results
As of June 3, 2020, a total of 13,538 PCR laboratory-confirmed COVID-19 cases had been included in this study. Only passive cases presented to the health institutions were included, with the active surveillance cases (totalling 1974) removed from the dataset. Of the included cases, 44.9% were Omanis, with the majority of these (71.1%) in the Muscat governorate, where Oman’s capital city, Muscat, is located. Among all included cases, 2769 infector–infectee pairs data were identified and included in the contact dataset. The median serial interval was estimated to be 6 (IQR 3–14).
The daily epidemic curve, by date of onset, for passive laboratory-confirmed COVID-19 cases classified by nationality—Omani and non-Omani—is shown in (Figure 1). The first cases registered in Oman were in two women returning from a visit to Iran on February 24. There was a doubling of cases from the middle of March to the end of May, noticed in both nationality groups.
Figure 1Daily COVID-19 epidemic curve by date of collection for Omanis compared with non-Omanis, February 24 to June 3, 2020.
Figure 2 shows the epidemic curve with the corresponding reproduction numbers between February 24 and June 3, 2020 for the total population. There was a sharp drop in Rt from 3.7 (95% CI 2.8–4.6) in mid-March to 1.4 (95% CI 1.2–1.7) towards the end of March. Subsequently, Rt started to rise again and fluctuated between 2.1 (95% CI 1.8–2.1) at the start of April and 1.3 (95% CI 1.2–1.5) in mid-April. This effect corresponded with the implementation of phase 2 of NPIs (major public health control—Table 1). Figure 2 also shows that Rt decreased to approximately 1.2 (95% CI 1.1–1.3) towards the end of April, when it started to rise again, corresponding with the start of phase 3 and the gradual easing of restrictions.
Figure 2Daily COVID-19 epidemic curve with corresponding reproduction number by date of onset for Oman, February 24 to June 3, 2020. The phases of public health control measures are also shown.
Analysis of the trend in Rt for the different nationality groups is presented in Figure 3, with the events corresponding with those in Table 1. Generally, there was a difference in Rt trends between the two groups. The response to major public health control (phase 2) was more evident in Omanis in reducing Rt to 1.09 (95% CI 0.84–1.3) by the end of March. Nevertheless, closure of Muscat and Mutrah (the old market area) led a marked reduction in Rt for the non-Omani group, reaching 0.9 (95% CI 0.8–1.1) by mid-April.
Figure 3Rt trend variation between Omanis and non-Omanis, February 24 to June 3, 2020. The phases of public health control measures and events are also shown. The arrows indicate the start of public health control measures. The mean value corresponds to the start date of the 7-day moving window.
Figure 4 illustrates an explanation for the differences in Rt between the two groups, showing a trend towards a cluster in the non-Omani population peaking around May 9. In addition, the number of cases belonging to a cluster tended to remain stable in the Omani group and to decline in the non-Omani group, while the number of sporadic cases not linked to a known cluster increased in both groups (Figure 4).
Figure 4The trend of incidence per 100,000 population (with 95% CI) for daily cases in Omanis and non-Omanis. After June 1, the number of cases that could not be linked to known clusters increased rapidly. A cluster of cases is where more than two cases can be linked to a common index case. A sporadic case is one that cannot be linked to an existing cluster.
With the increasing numbers of COVID-19 cases in Oman, our study showed the feasibility of using the time-varying Rt to assess and explain transmissibility dynamics and epidemic progression. We showed that there was a marked reduction in the reproduction number for COVID-19 infections in Oman in response to the major public health control introduced by the government. However, this reduction was less evident in the non-Omani group compared with the Omanis. Moreover, the closure of Muscat (specifically the old market area) drastically decreased Rt, presumably because of of high proportion of immigrants in Muscat. The estimated median serial interval in our study was 6 (IQR 3–14). Previous studies have estimated this parameter to be 4.6 days (95% CrI: 3.5–5.9) and 4.9 days (95% CI: 3.6–6.2) (
Estimating the serial interval of the novel coronavirus disease (COVID-19): a statistical analysis using the public data in Hong Kong from January 16 to February 15, 2020.
The daily number of new cases is known to be influenced by testing capacity. Moreover, a study from China demonstrated that changes in reporting rates substantially affect estimates of Rt (
Preliminary estimation of the basic reproduction number of novel coronavirus (2019-nCoV) in China, from 2019 to 2020: a data-driven analysis in the early phase of the outbreak.
). Thus, the initial Rt is influenced by testing frequency, contact tracing, and reporting of mild cases outside hospital. For these reasons, Rt will fluctuate from country to country. However, the use of Rt in assessing the transmissibility dynamics and epidemic progression in this Oman study was a crucial tool for establishing the influences of mitigation measures.
The development of the pandemic in Oman has followed the trend seen in other countries that could not curb the outbreak before the spread of COVID-19 in the community. Oman introduced most of the same mitigation measures implemented in other countries—stopping international travel by closing the points of entry (airports, seaports, and ground crossings), closing schools and shops, sending public- and private-sector workers home, introducing the use of masks in public spaces, and restricting movement in districts with particularly high densities of people, such as the old markets in Muscat. A modelling study from China showed how the comprehensive package of mitigation interventions China implemented in Hubei province caused Rt to decline rapidly to less than 1 in areas where NPIs were implemented (
First-wave COVID-19 transmissibility and severity in China outside Hubei after control measures, and second-wave scenario planning: a modelling impact assessment.
). Another study based on data from Wuhan, China found that mitigation measures reduced the median number of infections by more than 92% (IQR 66–97) and by 24% (13–90) at the peak and at the post peak of the outbreak, respectively (
On the basis of exponential curve prediction, and the assumption that the duration of infection ranges from 15 to 20 days, a study from Italy estimated the R0 to be between 2.8 and 3.3 (
First-wave COVID-19 transmissibility and severity in China outside Hubei after control measures, and second-wave scenario planning: a modelling impact assessment.
). These numbers are similar to those found in our study.
Throughout the outbreak, the number of cases has been higher in the migrant population compared with Omani nationals. This reflects differences in living conditions, where migrant workers often share dormitories, sleeping in close proximity—conditions that favour transmission of a respiratory virus.
Many of the disease clusters in the Omani population are present because Omanis tend to live in extended families (5–12 per household), thereby increasing the possible number of contacts for each primary case. A study from the UK found that a 74% reduction in the average daily number of contacts observed per participant (from 10.8 to 2.8) would be sufficient to reduce Rt from 2.6 prior to lockdown to 0.62 (95% CI 0.37–0.89) after the lockdown, based on all types of contact (
The rise in Rt by the end of March in our study was due to a considerable number of Omani students (roughly 8000) returning from abroad, even though they were placed in institutional or home quarantine. One study looked at the effects of travel restrictions in China and found that the travel quarantine imposed in Wuhan delayed the overall epidemic progression by only 3–5 days in mainland China, but had a more marked effect on an international scale, where case importation was reduced by nearly 80% until mid-February (
). The travel restrictions in Oman were effective, except for food supply trucks, which continued over the border with the United Arab Emirates, and between the locked-down governorates.
The situation in Oman is complicated by having two major and socio-economically different population groups—Omani nationals and migrants. So far, countries with this population structure have been unable to control the outbreak as efficiently as countries with homogeneous populations, with Taiwan, South Korea, and Hong Kong as prime examples (
), with 731 cases per 100,000 population confirmed by June 16. Qatar, another country with a large migrant population, has reported 3080 positive cases per 100,000 population since the start of the outbreak (
Easing mitigation efforts in Oman increased the number of cases, and it can be argued that the opening of small businesses and shops happened too early. However, the overall situation was complex because many small shop owners and self-employed had had no income for months. Nevertheless, with this transmissibility and the way in which Oman handles communicable diseases, the crude mortality rate from COVID-19 in Oman remains comparably low, at 0.4%.
Isolation and contact tracing reduce the time during which cases are infectious in the community, thereby reducing Rt. The overall impact of isolation and contact tracing, however, is uncertain and highly dependent on the number of asymptomatic cases (
Effectiveness of isolation, testing, contact tracing, and physical distancing on reducing transmission of SARS-CoV-2 in different settings: a mathematical modelling study.
). As case numbers rise, the burden of contact tracing and quarantining increases. A modeling study from the US indicated that, in such a situation, emphasis should be on physical distancing and contact tracing, whereas isolation should be prioritised for persons with a high risk of transmission to others (
An analysis of the four major clusters in South Korea estimated Rt to be 1.5 (95% CI 1.4–1.6). The intrinsic growth rate was estimated to be 0.6 (95% CI 0.6–0.7), and the scaling of growth parameter was estimated to be 0.8 (95% CI 0.7–0.8), indicating sub-exponential growth dynamics for COVID-19. The results indicate an early sustained transmission of COVID-19 in South Korea, and support the implementation of social distancing measures to rapidly control the outbreak (
One model found that the effects of physical distancing strategies varied across age categories, with the reduction in incidence being highest among school children and older individuals, and the lowest among working-age adults. Children are at a similar risk of infection to the general population, but less likely to have severe symptoms; hence they should be considered in analyses of transmission and control (
Another limitation of this study is that the testing capacity increased as the pandemic progressed. Using daily time series of COVID-19 incidence, epidemic curves of reported cases may not always reflect the true epidemic growth rate due to changes in testing rates, with limited diagnostic testing capacity during the early epidemic phase (
A third limitation of the study is the increasing number of sporadic cases by the end of the study period, indicating a lag in the identification and hence classification of the source of infection. This is likely due to the overwhelmed and fatigued public health workforce in the country. There is, therefore, a call for increasing and training this workforce to help it cope with current and future epidemics. In the short term, the introduction of advanced technologies (such as artificial intelligence and location tracking systems) will help public health professionals in this battle against COVID-19.
Conclusion
The use of real-time estimation of Rt allowed us to follow the effects of the mitigation strategies adopted by the government. Our analysis shows that the migrant population behaves differently from the native population, and that the COVID-19 infection is spreading more rapidly in non-Omanis mainly because of their special living conditions.
Authors’ contributions
Adil Al Wahaibi conducted the statistical analysis and wrote the draft manuscript. Abdullah Al Manji wrote the introduction and Table 1. Amal Al-Maani, Fatma Alyaquobi, Khalid Said Al Harthy, Bader Al Rawahi, and Amina Al Jardani reviewed and edited the manuscript. Eskild Petersen wrote the discussion and contributed to the overall manuscript. Seif Al-Abri supervised the study and participated in all stages of the manuscript.
Conflicts of interest
The authors declare that they have no conflicts of interest.
Funding
This study was funded by the Ministry of Health’s Directorate General for Disease Surveillance and Control, as part of its operational research.
Ethical approval
This study was approved by the Directorate General for Disease Surveillance and Control. There was no need for patients’ consent because the study was anonymous and used data produced for public health purposes.
Acknowledgment
The authors would like to thank Lesley Carson for her editorial assistance in finalizing the manuscript.
References
Bi Q.
Wu Y.
Mei S.
et al.
Epidemiology and transmission of COVID-19 in 391 cases and 1286 of their close contacts in Shenzhen, China: a retrospective cohort study.
Effectiveness of isolation, testing, contact tracing, and physical distancing on reducing transmission of SARS-CoV-2 in different settings: a mathematical modelling study.
First-wave COVID-19 transmissibility and severity in China outside Hubei after control measures, and second-wave scenario planning: a modelling impact assessment.
Estimating the serial interval of the novel coronavirus disease (COVID-19): a statistical analysis using the public data in Hong Kong from January 16 to February 15, 2020.
Preliminary estimation of the basic reproduction number of novel coronavirus (2019-nCoV) in China, from 2019 to 2020: a data-driven analysis in the early phase of the outbreak.
30668-8&id=gr1.jpg){kind=link}
30668-8&id=gr2.jpg){kind=link}
30668-8&id=gr3.jpg){kind=link}
30668-8&id=gr4.jpg){kind=link}