key: cord-0862227-8qm6eqto authors: Semakula, Muhammed; Niragire, FranÇois; Umutoni, Angela; Nsanzimana, Sabin; Ndahindwa, Vedaste; Rwagasore, Edison; Nyatanyi, Thierry; Remera, Eric; Faes, Christel title: The secondary transmission pattern of COVID-19 based on contact tracing in Rwanda date: 2021-06-08 journal: BMJ Glob Health DOI: 10.1136/bmjgh-2020-004885 sha: adf123685c8e3c68d23460d8de356ff68fc33645 doc_id: 862227 cord_uid: 8qm6eqto INTRODUCTION: COVID-19 has shown an exceptionally high spread rate across and within countries worldwide. Understanding the dynamics of such an infectious disease transmission is critical for devising strategies to control its spread. In particular, Rwanda was one of the African countries that started COVID-19 preparedness early in January 2020, and a total lockdown was imposed when the country had only 18 COVID-19 confirmed cases known. Using intensive contact tracing, several infections were identified, with the majority of them being returning travellers and their close contacts. We used the contact tracing data in Rwanda for understanding the geographic patterns of COVID-19 to inform targeted interventions. METHODS: We estimated the attack rates and identified risk factors associated to COVID-19 spread. We used Bayesian disease mapping models to assess the spatial pattern of COVID-19 and to identify areas characterised by unusually high or low relative risk. In addition, we used multiple variable conditional logistic regression to assess the impact of the risk factors. RESULTS: The results showed that COVID-19 cases in Rwanda are localised mainly in the central regions and in the southwest of Rwanda and that some clusters occurred in the northeast of Rwanda. Relationship to the index case, being male and coworkers are the important risk factors for COVID-19 transmission in Rwanda. CONCLUSION: The analysis of contact tracing data using spatial modelling allowed us to identify high-risk areas at subnational level in Rwanda. Estimating risk factors for infection with SARS-CoV-2 is vital in identifying the clusters in low spread of SARS-CoV-2 subnational level. It is imperative to understand the interactions between the index case and contacts to identify superspreaders, risk factors and high-risk places. The findings recommend that self-isolation at home in Rwanda should be reviewed to limit secondary cases from the same households and spatiotemporal analysis should be introduced in routine monitoring of COVID-19 in Rwanda for policy making decision on real time. Rwanda was one of the first African country to prepare against the spread of COVID-19 in January 2020 by screening and recording everyone who travelled to China or in any other country that had a declared COVID-19 case. A total lockdown was imposed when the country had only 18 COVID-19 confirmed cases on 22 March 2020. Using intensive contact tracing, several COVID-19 cases were identified, with the majority of them among returning travellers and contacts of identified cases. Lockdown and contact tracing strategies interrupted the chain of transmission and provided time to the Rwanda COVID-19 Joint Task Force to establish treatment centres and allow readiness of response teams. The contact tracing has a long history in Rwanda and has played a major role in controlling various infectious diseases including Ebola virus disease, severe acute respiratory syndrome and Middle East upper respiratory syndrome. 1 Several of these more recent outbreaks have occurred on the African continent, allowing African scientists and public health officials to build their teams and capacity. 2 As part of Rwanda's COVID-19 response, the government applied a comprehensive surveillance approach that includes screening, testing and isolation (quarantine) on real time using innovative technology. It defines a suspect case, a confirmed case and a contact as follows. A suspect COVID-19 case is a patient with acute respiratory illness such as fever and at least one symptom of respiratory disease (cough, shortness of breath, fatigue, anorexia, myalgia, headache, rhinorrhoea, sore throat, vomiting, diarrhoea, smell and taste disorders, clinical signs of pneumonia and abnormal chest) or someone who has travelled in the last 14 days to an area that has reported confirmed case(s) of COVID-19. 3 A confirmed COVID-19 case is a person with laboratory confirmation of COVID-19 infection using reverse transcription PCR testing for SARS-CoV-2, irrespective of clinical signs and symptoms. A contact of a COVID-19 case was defined as any person who had contact with a COVID-19 case within a timeframe ranging from 72 hours before the onset of symptoms of the case to 14 days after the onset of symptoms. 4 Rwanda's contact tracing strategy was to test each suspect COVID-19 case and trace all contacts, the travel history and visited places in the past 14 days for each confirmed case by a protected epidemiologist. Each identified contact was tested at least twice within a period of 2 weeks. In this way, Rwanda has traced more than 10 000 contacts between 14 March 2020 and 20 July 2020, with the aim to identify potentially infected persons and prevent onward transmissions of secondary cases. The first 300 COVID-19 cases identified in Rwanda were predominately imported cases, and most of them were sent to quarantine sites immediately on arrival at airport. Despite early detection strategies in place, some imported cases were identified after being in contact with their families, relatives, friends and coworkers. For a better management and control of COVID-19, understanding the geographical risk regions is the primary and fundamental step to take. 5 Spatial analysis is of great help in understanding the spread of infectious diseases, spatial association and predicting local outbreaks. Spatial analysis is critical to develop public health policies during early stages of COVID-19. 6 Until now, studies evaluating the spatial spread of the COVID-19 pandemic in Rwanda are limited. However, this study contributes on identifying high-risk areas in Rwanda at subnational level. Furthermore, estimating the secondary attack rate (SAR) is important in epidemiology. Attack rates may contribute significantly in designing and developing infection control and prevention policies for COVID-19 that can contribute to limit the spread of the disease. The SAR is defined as the ratio of the number of cases occurring within the incubation period on exposure to a primary case and the total susceptible contacts. 7 In this study, we used contact tracing data in Rwanda to understand geographic patterns of COVID-19 to inform targeted interventions, estimate SARs and identify factors associated with increased transmission probability among contacts of COVID-19 index cases. Findings from our study will help identify high-risk areas and groups and provide evidence to develop targeted interventions to reduce the burden of COVID-19 in Rwanda and similar settings. We have used COVID-19 contact tracing data, recorded from 14 March 2020 to 20 July 2020, from the Rwanda health management information system (RHMIS). The data set contains 2216 index cases, among them 626 cases were identified and isolated in government-supervised isolation centre before meeting anyone at point of entries. Only 1590 have completed all details, and their contacts (11 809) were identified through contact tracing. Rwanda has a digitised real-time COVID-19 data system. Everyone tested COVID-19 positive, laboratory system automatically sends his or her details to contact tracing system that enables the tracers to start investigation immediately. 8 Investigation is performed by interviewing persons at the treatment centre or at home for those under home-based care programme. However, extra information about contacts are also generated through algorithms in place that are using combination of telephone numbers and network connectivity used by a positive case in the last 7 days prior to test positive. Once all contacts are identified, the tracers record their details in the system and are followed up for a period of 7-14 days from the last day they met a confirmed COVID-19 case. 4 All contacts are tested at least twice during the follow-up period. For both cases and contacts, detailed demographics, clinical diagnostics, laboratory details, travel history, employment details, residence and people living with him/her are recorded in the contact tracing data system hosted at Rwanda National Data Center under management of Rwanda Biomedical Centre. Figure 1 shows an example of the systematic contact tracing of one confirmed case that generated seven secondary cases: two secondary cases and five tertiary cases. The black cycles in figure 1 are positive cases and include the case numbers and the date that the cases tested positive. Red, yellow and green cycles are the contacts identified, investigated and tested, and correspond to high-risk, medium-risk and low-risk contacts. A similar exercise was done to each and every positive case in Rwanda. Rwanda focused much on contact tracing, but unfortunately, less attention has been paid to the risk factors related to COVID-19 transmission. The attack rate was estimated as the percentage of contacts who were later confirmed to be infected with SARS-CoV-2 among all secondary contacts of the index case. The estimation was presented by gender, age, the BMJ Global Health type of relationship to the index case, place of infection and symptoms. Moreover, we explored the age distribution of the index cases and their contacts, as well as of the secondary cases to see which age groups are most likely to interact. Standard non-parametric tests were used, mainly Fisher's exact test, the Kruskal-Wallis and the Mann-Whitney U test to compare characteristics between demographic groups. 9 We used Bayesian disease mapping models to assess the spatial pattern of COVID-19 and to identify areas characterised by unusually high or low relative risk. In addition, we used multiple variable conditional logistic regression models without interaction term (M1) and with interaction term (M2) to identify factors associated with COVID-19 infection. We calculated the standardised incidence ratio (SIR) in each area taking the age and sex structure in an area into consideration to investigate the COVID-19 risk in an area. [10] [11] [12] We extended our methods to take into account spatial dependence among neighbouring areas using Bayesian disease mapping approaches. A Bayesian disease mapping model consists of three components: the data model (ie, the distribution of the data given the parameters), the process model (ie, a description of underlying spatial trend) and the parameter model (ie, the prior distribution of the parameters to be estimated). 13 The data model is given by Yit ∼ Poisson (Eit θit); (1) where a Poisson distribution is appropriate since disease data are counts of number of cases and are non-negative. It is assumed that the mean is a product of the expected count Eit and the relative risk θit. The process model describes the underlying structure of the relative risks. We used the Besag-York-Mollie model, which is the CAR convolution model with two random effects, one spatially structured area-specific random effect and one unstructured area-specific random effect 14 15 log(θi) = α + ui + νi (2) where, ui is the spatially structured area-specific random effect that allows for smoothing among adjacent areas, namely 14 With δ_i and n_δi, respectively, the set of neighbours and number of neighbours for a specific area i. The unstructured component v_i is modelled using as a Gaussian process and allows for extra heterogeneity in the counts due to unobserved (and spatially unstructured) risk factors. We have used Integrated Nested Laplace approximation for estimation of this model. 16 We make use of excess probability to investigate the hotspots. 10 17 Furthermore, to identify individual risk factors, logistic regression is used with multiple explanatory variables, which took into consideration the interaction terms of the age and the type of relationship to the index case. 18 The response variable is whether the secondary case tested positive or negative to COVID-19, which is a binary outcome for each contact to the index case and which we denote as Z. The probability of testing COVID-19 positive is defined as π, and consequently, the probability of testing negative is 1-π. The multiple logistic regression equation can be written as Z~binomial(pi) with with p the number of independent variables (gender, age, relationship to the index cases and intersections between the independent variables). β0 is the intercept and quantifies the probability of infection in the reference group, that is, the group that takes on the value 0 for all variables; βi is the regression coefficient corresponding to risk factor xi , with exp(beta_i) the multiplicative increase or decrease in the odds to test positive for a one unit increase in xi , while all other covariates remain unchanged. The probability of testing positive can be written as The regression parameters are estimated using maximum likelihood. The regression coefficient βi are transformed into OR by taking exp ( βi ). The OR value reflects the strength of association between exposure and the outcome. 18 For the goodness of fit of multiple variable logistic models, we used likelihood ratio test. The likelihood ratio statistic G 2 (M0|M1) = −2(L 0 -L 1 )was defined to test whether model parameters are zero, given that M1 holds, by comparing the log likelihood L1 for fitted model M1 with L0 for a simpler model M0. Data used for this study were anonymised, deidentified and routinely collected during outbreak investigation and response data maintained by RHMIS, under Public Health Surveillance and Epidemic Preparedness and Response Division of Rwanda Biomedical Center. No participants were involved directly in data collection; therefore, their consent was waived by Rwanda National Ethics Committee, which also approved the use of outbreak and investigation team presented here. The Rwanda Ministry of Health also granted approval for data access and use to the principal investigator (MS) for the purpose of improving COVID-19 contact tracing in Rwanda. The results are presented in two parts. The first part includes COVID-19 attack rate by exposure routes. The second part provides a geographical analysis of the number of contacts and secondary cases. The mean number of contacts per index cases was 7, though the number of contacts per case ranges from 0 to 301 contacts (IQR 0-9, figure 2A ). We have observed five index cases with more than 200 contacts in less than 1 week. The mean number of secondary cases per index case was 0.25 (IQR 0-0 and range 0-6, figure 2B ). The individual SAR, defined as the ratio of the number of secondary cases to the number of contact, ranges from 0.0 to 0.5, with mean equal to 0.017 (IQR 0.0-0.0, figure 2C ). The index case that had 301 contacts was among one of the superspreaders of COVID-19, with five of his contacts tested positive. Another superspreader in this study was a school teacher with 96 contacts, resulting in six positive secondary cases. We have provided the demographic details of the index cases and their contacts in table 1. With 1590 index cases, we traced 11 809 contacts of which 209 are secondary cases and 11 600 uninfected contacts. The overall attack rate is 1.77% (95% CI 1.55% to 2.02%). The majority of both the index cases (76.13%) and the contacts (51.35%) were male and aged between 20 and 39 years. The SAR for males is 2.06% (95% CI 1.73% to 2.46%), slightly higher as compared with females with overall BMJ Global Health rate of 2.11% (95% CI 0.67% to 6.42%), and contact among coworkers (4.98%), with an attack rate of 2.55% (95% CI 1.54% to 4.19%). All the contacts were screened for symptoms and 278 of the 11 809 contacts (2.35%) had symptoms, with the most common symptoms being cough (85.25%) and chest pain (14.30%). However, only few of them tested positive, and among the 209 secondary cases, only 7 showed symptoms. Some care is needed with the above numbers, as these effects can be confounded by other covariates. Therefore, a multiple logistic regression is performed to investigate the joint effect of gender, age and the type of relationship. From this analysis (table 2), we can see that the OR for male index cases to have a secondary case is higher as compared with a female index case (OR 1.43, 95% CI 1.08 to 1.89). Contacts of index cases aged 50+ years had lower odds as compared with those below 20 years old in Rwanda (OR 0.45, 95% CI 0.26 to 0.79). The household members of index cases had nearly twice the risk of infection compared with contacts that met index cases in public places such open market, shops or pubs (OR 1.92, 95% CI 1.17 to 3.15). However, the estimates are different for different age groups within the same household, with the OR for index cases aged <20 years being lower (OR 0.87, 95% CI 0.21 to 3.64), and for index cases aged 20, 49 and 50+ years higher (OR 2.57, 95% 0.55 to 12.01 and OR2.9, 95% CI 0.36 to 23.24, respectively). While not statistically significant-possible due to limited number of secondary infections-this hints that there is a difference in contact behaviour depending on the type of contact and age. The overall OR of coworkers as compared with contacts at public places to be infected is 1.69 (95% CI 0.89 to 2.89), though (while not statistically significant) it is higher for index cases <20 years of age (OR 2.8, 95% CI 1.16 to 6.76) and smaller for index cases aged 20, 49 or >50 years (OR 0.49, 95% CI 0.15 to 1.53) and OR 0.39, 95% CI 0.04 to 3.67, respectively). Also contacts at social activities and on transport are of somewhat higher risk as contacts at public places, though not statistically different. Finally, we investigate the interaction of age of index case and age of contacts. Therefore, we compared the age groups (0-9, 10-19, 20-29, 30-39, 40-49, 50-59 and 60+ years) of contacts and index cases. Figure 3 shows the attack rates (ARs) for men and women, with the AR for males being highest among contacts aged 20-29 years and 30-39 years who interact more with index cases aged 50-59 years or above 60 years. BMJ Global Health Figure 4 shows that the age-age interactions on the risk for secondary infection varies in different contact settings, but was largely influenced by contacts in household settings and among coworkers. There are major geographical differences in the risk of making contact to an index case. Figure 5 in table 3 , together with the SIR. In some sectors, the risk of making contact to an index case was very high (SIR >10), mainly in the southwestern province. In Rusizi District, more than eight sectors had SIR of above 8. The map presented in figure 5 shows that contacts to index cases are clustered mainly in the City of Kigali, which is located in the centre of Rwanda and in the southwest of the country (Rusizi and Nyamasheke Districts). We observe that most COVID-19 cases are in the most densely populated regions. However, not all densely populated areas had COVID-19 confirmed cases. The relative risk (RR >1) of secondary cases in a sector as compared with the overall risk in the country (SIR) showed that the contacts who tested positive for COVID-19 were also clustered in the southwestern part of Rwanda, City of Kigali and northeastern part of Rwanda ( figure 5 (top right) ). The area-specific probability figure 5 (top right and bottom) are based on a spatial model and uses information only from within the area; therefore, it might give uncertain estimates for small areas. We estimated relative risks that take into account the uncertainty and spatial association. Using this approach and stratifying by sex, we identified new clusters in northeast of Rwanda for both men and women. Among women, we identified another cluster in the southeast of Rwanda in addition to the one identified using SIR. The administrative map of sectors of Rwanda is provided as online supplemental figure S1. Spread of COVID-19 in Rwanda largely occurred in Kigali (capital city). Age, male and sharing a home or office with an index case were factors driving the risk for secondary cases. In Rwanda, men accounted for 60% of all cases. At population level, markets/pubs/ shops were the places where the largest number of secondary cases occurred due to the large number of contacts being made at this place. Spread seemed to be driven by different age groups mixing, and younger individuals (20-29 and 30-39 years) seemed to drive the spread of the epidemic. We found that in Rwanda, index cases and their contacts are clustered mainly in the City of Kigali, which is located in the central of Rwanda, in the southwest of country, Rusizi and Nyamasheke Districts. Among the first hundred COVID-19 cases in Rwanda, most were imported and identified as having recently travelled from Dubai, United Arab Emirates, which triggered contact tracing in City of Kigali. Regarding the Western province, the areas bordering Republic Democratic of Congo and Rwanda, Rusizi and Nyamashake had a SARS-CoV-2 community transmission due to legal and illegal cross-border movement during lockdown with inconsistency of using masks, keeping distance and washing hands regularly. Rusizi and Nyamasheke were the highest risk areas for COVID-19 in Rwanda. The contact tracing data showed that in the highest risk areas, the number of contacts per index case and per areas were higher as compared with other areas. There was an increased transmission of infection among areas that were identified late as COVID-19 high-risk areas. The sectors of Rusizi and Nyamasheke Districts, such as Kamembe, Gihundwe and Mururu, recorded the highest number of identified contacts, as figure 5 (top left) showed in the map. Those sectors remained in lockdown for more than 3 months, while temporary localised lockdown duration in Rwanda was about 3 weeks maximum. We estimated relative risks and findings provided two new extra clusters in northeast for male and southeast for females. The spatial models can be useful in identifying clusters in low-risk COVID-19 areas. Spatial analysis is of great help in understanding the spread of infectious diseases. The spatial association was identified to be key during the early stages of the COVID-19 pandemic in mainland China. 6 Five index cases with more than 200 contacts were identified. The first superspreader identified was a school teacher with 96 contacts, with among them 6 tested positive. Another one was taxi driver with 301 contacts of which five tested positive. The overall ratio of index case per secondary cases in this study is 8:1, while for the first superspreader is 1:6. Some other studies 19 reported rapid increase of new cases linked to superspreaders and showed that with progress of diagnostic technology, potential superspreaders may be discovered in the future. 20 We have observed that only 2.87% of the secondary cases had symptoms. In order to prevent the spread of SARS-CoV-2, asymptomatic spreaders need to be studied to detect them early. Though the overall AR is less than 2% among secondary cases, this varied by a number of factors. We found that AR was higher in men and people under 30 years old, household members living with index case and people attending social activities or workmates of index case. The AR in Rwanda was smaller as compared with the AR in cities in China and USA (above 10%). 21 Unlike findings from other studies in Asia, 22 the attack rate among the elderly was lower (0.13%) in Rwanda for people aged above 60 years old. This might be partially explained by the fact that contacts aged 60 years and above could have less physical activities and movement as compared with people age 20-29 years old in Rwandan setting. In addition, Rwanda's population as a whole is quite young. Over 50% of Rwandans are under 20 years old, and the median age is 22.7 years old. 23 We found that men experienced nearly twice the odds of infection compared with women. This finding differs from studies conducted elsewhere that observed female contacts that were more likely to be infected by SARS-CoV-2 than male contacts. 9 22 24 Despite the enforcement of prevention measures of COVID-19 in Rwanda, especially in public places, we observed that household members of index are two times at risk as compared with contacts that met index cases in public places such as market, shops or pubs (OR 1.92, 95% CI 1.17 to 3.15). This could be explained by the fact that most of COVID-19 cases in Rwanda were asymptomatic, and several studies reported that both symptomatic and asymptomatic person are infectious. 25 26 Therefore, most infections might have been occurred unknowingly in families and relatives living in the same households. Though Rwandans comply with several public health measures in place such as hand washing, universal face mask and physical distance, we have realised that more transmission occurred at markets. However, the individual risk appeared greater within the household (in spite of lockdowns). In summary, it is crucial to understand the transmission channels and situation of SARS-CoV-2 as a new virus to prevent and control its spreading in the community. The analysis of contact tracing data using spatial modelling allowed us to identify high-risk areas at subnational level in Rwanda. Estimating risk factors for infection with SARS-CoV-2 is vital in identifying the clusters in low spread of SARS-CoV-2 subnational level. In a pandemic period, it is imperative to understand the interactions between the index case and contacts to identify superspreaders, risk factors and high-risk places. Understanding spread of the virus in a specific country leads to contextually appropriate interventions that slow the spread of SARS-CoV-2 and saves lives. The findings recommend that self-isolation at home in Rwanda should be reviewed to limit secondary cases from the same households. The contact tracing strategy should include routine surveillance and testing of highrisk groups to limit community transmission. Spatial and spatiotemporal analysis should be introduced in routine monitoring of COVID-19 in Rwanda for policy making decision on real time. Contributors MS and CF conceived the presented idea. MS developed the theory and performed the analysis. CF verified methods and validated the analysis. NF, AU, SN, VD, ER, TN and ER reviewed the manuscript. CF supervised the findings of this work. All authors discussed the results and contributed to the final manuscript. Funding This study was funded by European Union's innovation programmeproject EpiPose (101003688). Map disclaimer The depiction of boundaries on the map(s) in this article does not imply the expression of any opinion whatsoever on the part of BMJ (or any member of its group) concerning the legal status of any country, territory, jurisdiction or area or of its authorities. The map(s) are provided without any warranty of any kind, either express or implied. Competing interests None declared. Patient consent for publication Not required. Provenance and peer review Not commissioned; externally peer reviewed. Data availability statement Data are available in a public, open access repository. All data relevant to the study are included in the article or uploaded as supplementary information. Data sets used are available for public after manuscript is being accepted for publication. 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. Open access 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/. Eric Remera http:// orcid. org/ 0000-0002-6084-8877 How Africa has tackled covid-19 COVID-19: shining the light on Africa Ministry of Health. Coronavirus disease 2019, National preparedness and response plan RBC. Rwanda COVID-19 clinical managment guidelines Spatial modeling, risk mapping, change detection, and outbreak trend analysis of coronavirus (COVID-19) in Iran (days between Spatial epidemic dynamics of the COVID-19 outbreak in China Secondary attack rate of COVID-19 in household contacts: a systematic review QJM Lessons learned from Rwanda: innovative strategies for prevention and containment of COVID-19 Household secondary attack rate of COVID-19 and associated determinants in Guangzhou Bayesian spatio-temporal modeling of malaria risk in Rwanda Spatial and Spatio temporal Bayesian models with R-INLA Bayesian disease mapping for public health Bayesian image restoration, with two applications in spatial statistics Spatial statistics and Bayesian computation Comparing INLA and OpenBUGS for hierarchical Poisson modeling in disease mapping Using non-exceedance probabilities of policy-relevant malaria prevalence thresholds to identify areas of low transmission in Somalia Categorical data analysis UK Understanding of COVID-19 based on current evidence Transmission of 2019 nCov infection from an asymptomatic in Germany Contact settings and risk for transmission in 3410 close contacts of patients with COVID-19 in Guangzhou, China, a prospective cohort study Risk factors associated with COVID-19 infection: a retrospective cohort study based on contacts tracing Population of Rwanda:Results from 4th census Epidemiology and transmission of COVID-19 in 391 cases and 1286 of their close contacts in Shenzhen, China: a retrospective cohort study Transmission of 2019-nCoV infection from an asymptomatic contact in Germany A familial cluster of pneumonia associated with the 2019 novel coronavirus indicating personto-person transmission: a study of a family cluster