SARS-CoV-2 transmission in different settings: Analysis of cases and close contacts from the Tablighi cluster in Brunei Darussalam

Surveillance and case identification

Brunei’s Ministry of Health has responsibility for communicable disease surveillance. Since January 23, case definitions have been in place for suspected COVID-19 cases. Initially, individuals with fever, respiratory symptoms, and a travel history to a high-risk area were tested for SARS-CoV-2. Over the next weeks, the program expanded to include: (i) contacts of a confirmed case (regardless of symptoms); (ii) individuals admitted to an inpatient facility with pneumonia; and (iii) individuals who present to a health facility with influenza-like illness (ILI) for the second time within 14 days. On March 21, Brunei started testing and isolating all travellers and returning residents. On March 25, SARS-CoV-2 sampling at selected sentinel health centers was introduced, followed by mandatory screening for selected groups of foreign workers on April 7.

A confirmed case is defined as a person who tested positive for SARS-CoV-2 through real-time reverse transcriptase polymerase chain reaction (RT-PCR) test on nasopharyngeal (NP) swab. The first positive case in Brunei was detected on March 9, having met the testing criteria as a person with fever and cough, and had a recent travel history to Kuala Lumpur (where he had attended the Tablighi event). This marked the beginning of Brunei’s first COVID-19 cluster.

Epidemiological investigation

Under the Infectious Disease Act, the Ministry of Health conducts epidemiological investigation and data collection and analysis of each case and close contact. Following the identification of the first case, we interviewed him for demographic characteristics, clinical symptoms, travel history, activity mapping, and contact history. We identified his attendance at the Tablighi event in Malaysia, and found that a number of other Bruneians had also participated in the same event. We subsequently obtained the contact details of all Bruneian participants.

NP swabs were collected from all identified event participants and tested with RT-PCR. Those who tested positive were admitted to the National Isolation Centre (NIC), while those tested negative were quarantined for 14 days from their return to the country at a designated community quarantine facility where symptom and temperature screening was conducted daily. Those who developed symptoms were re-tested. Activity mapping of the confirmed cases was conducted, and contact tracing initiated.

A close contact is defined as any person living in the same household, or someone within one meter of a confirmed case in an enclosed space for more than 15 minutes. All close contacts of confirmed cases were tested with RT-PCR. Those who tested positive were admitted to NIC, while those tested negative, were placed under home quarantine for 14 days from last exposure to the confirmed case. For individuals under home quarantine, their compliance and health status were monitored daily, through video calls or face-to-face assessments. Those who developed symptoms during home quarantine were re-tested.

Clinical management

All confirmed cases were treated and isolated at NIC. We obtained clinical information on their history, examination, laboratory and radiological results from digital inpatient records on the national health information system database. Cases were discharged following two consecutive negative specimens collected at ≥ 24-hour intervals.

Cases

We categorised cases into two groups: primary cases (those who were presumably infected at the Tablighi event in Malaysia), and non-primary cases (those who did not attend the Tablighi event but had an epidemiological link to the primary cases).

For each case, the symptom status was recorded and classified as: (i) symptomatic, if symptoms were reported at NP swab collection; (ii) presymptomatic, if symptoms were reported after NP sampling but during admission; or (iii) asymptomatic, if no symptoms were ever reported up to discharge.

Close contacts

We classified the type of close contact into five settings: household, relatives, workplace, social, and local religious gathering.

Household setting is defined as those living in the same household. Among the household members, their relationship with a case was also recorded and classified as spouse, child, and others (which includes grandchild, parent, grandparent, sibling, housekeeper, and relatives living in the same household). Relatives setting is defined as those who live outside the household but are related to the case. Workplace setting is defined as contacts encountered in the workplace or school. Social setting is defined as those encountered during travel or in social events. Lastly, local religious gathering is defined as those who attended a local religious event in Brunei on March 5, which ran throughout the night with participants staying overnight.

Data analysis

Group comparison was done between primary and non-primary cases using Chi-square, Fisher’s Exact, or Mann-Whitney’s tests as appropriate. Incubation period was calculated based on cases whose dates of exposure and symptom onset were clear. Serial interval (SI) was calculated by subtracting the date of symptom onset of the infectee from the infector; only symptomatic and pre-symptomatic infector-infectee pairs with clear epidemiological links were included.

For each setting, the attack rate (AR) was calculated by dividing the number of positive contacts by the total number of close contacts (that is, the proportion of contacts that tested positive). To identify risk factors of infection, log-binomial regression analysis was applied to estimate the risk ratio for gender, age, and setting. Household contacts were also divided into three subcategories: spouse, child, and others. Further stratification was done to assess differences in infector symptom status across settings. The 95% confidence interval (95% CI) was estimated using the normal-approximation method, or binomial method if the count was less than five.

The mean observed reproductive number, R, and the distribution of personal reproductive numbers in each setting were calculated from the number of infected close contacts caused by each primary case. The 95% CI was estimated based on the Poisson distribution.¹³

All analyses were conducted using Microsoft Excel and R (ver. 3·6·3).¹⁴ A p-value <0·05 was considered as statistically significant. Ethical approval was obtained from the University Research Ethics Committee, Universiti Brunei Darussalam (Ref: UBD/OAVCR/UREC/Apr2020-05).

Epidemiological characteristics

We identified 75 individuals in Brunei who attended the Tablighi event in Malaysia. Of these, 19 were positive for SARS-CoV-2 resulting in 52 additional cases transmitted locally, bringing the total cluster size to 71. The figure illustrates the epidemiological links in the cluster by generation in the transmission chain. There were 34 (47·9%), 13 (18·3%), and 5 (7·0%) cases in generations one, two, and three, respectively.

• Download figure
• Open in new tab

Figure – The Tablighi cluster in Brunei Darussalam, with epidemiological links illustrated

Table 1 shows the demographic and clinical characteristics of the cases in the cluster. The median age was 33·0 years [Interquartile range (IQR): 29·5 years], majority male (65·7%, n = 46), and 5 cases (7·1%) had pre-existing chronic conditions. Compared to non-primary cases, primary cases were significantly older and predominantly male. More than three-quarter of the cases (77·4%, n=55) were diagnosed and immediately admitted to the NIC within 5 days.

A substantial proportion of cases were presymptomatic (31·0%, n=22) or asymptomatic (12·7%, n = 9). 40 (56·3%) patients reported symptoms during contact tracing investigation. The most commonly reported symptoms were fever, cough, and sore throat. Only one (1·4%) and two (2·8%) cases were critical and severe, respectively.

The incubation period was calculated from eight cases that had confirmed epidemiological links and had attended the March 5 religious gathering in Brunei. Using their attendance date as the exposure date, the median incubation period was 4·5 days (IQR = 2·75 days, range: 1 to 11 days). Based on 35 symptomatic infector-infectee pairs, the mean (SD) of SI was 4·26 (4·27) days, ranging from −4 to 17 days. Four pairs (11·4%) had negative SI values. The SI distribution resembles a normal distribution (Supplementary Figure S1).

Transmission characteristics

Of the 1755 close contacts in the Tablighi cluster, 51 local transmissions were detected, giving an overall non-primary AR of 2·9% (95% CI: 2·2, 3·8). Case 121 (see Figure) was excluded from this analysis because he was not detected during contact tracing. The highest AR was observed among spouse (41·9% [95% CI: 24·1, 60·7]), followed by the local religious gathering (14·8% [95% CI: 7·1, 27·7]) and children (14·1% [95% CI: 7·8, 23·8]). The overall household AR is 10·6% (95% CI: 7·3, 15·1).

The multiple log-binomial regression model revealed that type of close contact was the only statistically significant variable (p<0·001) (Table 2). When compared to social contacts, spouses of positive cases had the highest adjusted risk ratio of getting the infection (45·20 [95% CI: 16·8, 156·1]), followed by the local religious gathering attendees (15·60 [95% CI: 4·8, 59·9]), and children of positive cases (14·09 [95% CI: 4·8, 51·5]).

ARs also differed by symptom status of the infector (Table 3). ARs in households where the infectors were symptomatic were higher than those who were asymptomatic or presymptomatic (14·4% vs. 5·4%). AR for the local religious gathering could not be calculated as three primary cases at the event had different symptom status; hence we could not ascertain how transmission occurred. No significant difference was observed in other settings.

Table 4 summarises the transmission characteristics in each setting. The mean observed R was highest in the household setting (0·67 [95% CI: 0·44, 0·96]). The distribution of the observed R in the household setting was skewed towards zero (Supplementary Figure S2). 71·4% of household infections (20 of 28 positive contacts) were from 16·7% of cases (7 of 42 possible links).

SSE

Within this cluster, 38% of all cases were participants at an SSE – 19 (26·7%) from the Tablighi event in Malaysia, and eight (11·3%) from the local religious gathering in Brunei. Notably, 19 of the 75 Brunei attendees at the Tablighi event tested positive. Assuming a representative sample, this would suggest an AR of 25%, implying that approximately 4,000 cases (of an estimated 16,000 participants) may have been infected at that single event alone. Moreover, we find that the highest overall non-primary AR and mean observed R was at the local religious gathering (AR = 14·8%, R = 2·67), which was higher than that observed in the household setting (AR = 10·6%, R = 0·67). These observations suggest a role for mass gatherings in facilitating SARS-CoV2 transmission.

Over the course of this investigation, we identified several common characteristics at both the local religious gathering and the Tablighi event in Malaysia.¹² First, significant numbers of people gathered in an enclosed area for a prolonged time period. Second, attendees had a recent travel history – the Tablighi event in Malaysia drew participants from across the world, while the local religious gathering had at least three individuals who had recently returned from Malaysia. Third, communal sharing of sleeping areas and toilets, and shared dining were observed. We propose that these three characteristics are hallmarks for the development of SSE for SARS-CoV-2 transmission and can be used as ‘red flags’ by health authorities in their risk assessment and mitigation strategies for preventing and detecting high risk activities including mass gatherings, and other institutional settings such as care homes, prisons and dormitories.

Variability of non-primary AR across different settings

To a lesser degree, our observations on the within-household transmission are similar to that observed for the two religious gatherings. Out of 16 household contacts who subsequently became first generation cases, 10 (62·5%) of these were from just three primary cases. As such, even within similar settings, we can expect wide variability in transmission patterns. This observation supports our finding of a moderately high household AR but an observed R of less than one, suggesting that transmission is driven by a relatively small number of cases.⁶ High ARs in spouses and children reflect intimate relationships with high degree of interaction, close proximity, and in the case of the spouse, sleeping in the same room. Concordant with our findings for the local religious gathering, we suggest that encounters among groups of people that involve close proximity in enclosed settings for prolonged time periods (at least for one night) is a main driver of SARS-CoV-2 transmission.

Our overall non-primary AR result of 10·6% in the household setting is comparable to other studies that used contract-tracing datasets.^15–18 A study near Wuhan, China¹⁹ reported a higher AR of 16·3%; however, this study detected 56·2% of their cases more than five days after symptom onset. By contrast, 77·4% of the cases in our study were detected and isolated within five days of symptom onset, suggesting that early case isolation can reduce AR.

We note the low non-primary AR (<1%) and mean observed R (<0·3) for workplace and social settings. While moderate physical distancing was implemented in Brunei following the identification of this cluster, there was no community quarantine or lockdown, public services and businesses remained open, and no internal movement restrictions were imposed.

Combined with our observations on the role of SSE in driving SARS-CoV-2 transmission, we suggest that full lockdown measures that adopt a blunt approach by restricting all movement and casual social encounters can be avoided, in favour of a more targeted approach that includes a combination of case isolation, contact tracing, and moderate levels of physical distancing that take into account the ‘red flags’ for mass gatherings that we identified earlier. In fact, given the relatively high non-primary AR observed in the household compared to workplace and social settings, we suggest that lockdowns, in the absence of case isolation, may even be counterproductive. Other studies support our observations.^20,21 If done effectively, both contact tracing and case isolation approaches are shown to control the COVID-19 outbreak during its early stage.²² Modelling studies using South Korean data showed that less extreme physical distancing measures can help to suppress the outbreak,²³ indicating less need for national lockdown measures.

Symptomatic vs. asymptomatic and presymptomatic transmission

We identified several environmental (settings) and behavioural factors that potentially account for higher ARs observed in mass gatherings and in the household. In order to assess the impact of host factors in driving transmission, we compared non-primary AR in symptomatic vs. asymptomatic and presymptomatic individuals, considering the high proportion of asymptomatic (12·7%) and presymptomatic (31·0%) of cases.

While there are case reports of presumptive asymptomatic and presymptomatic transmission,^24,25 observational studies quantifying such transmission are few. A study from Ningbo, China analysed the overall attack rates in symptomatic vs. asymptomatic cases and did not find significant difference between the two groups.²⁶ Another study, on the other hand, reinterpreted the same data and theorized that under certain conditions, symptomatic cases could be more transmissible than asymptomatic ones.²⁷ In fact, our overall risk ratio (RR) for symptomatic vs. asymptomatic and presymptomatic cases shows no significant difference between the two groups (1·28 [95% CI: 0·69, 2·37], Table 3). However, we suggest that the overall RR masks the true picture in transmissibility when different settings are taken into account.

In our study, we do not find a significant difference in AR in non-household settings. These settings usually practice some form of non-pharmaceutical interventions (NPI)—individuals with moderate and severe symptoms may be on medical leave, and it is reasonable to expect some physical distancing would be practiced by contacts of persons who display visible symptoms. This is less feasible within the household setting; hence we suggest that transmission occurs more frequently at the household level where control measures are less practical. We observed that the household AR for symptomatic cases (14·4%) is higher than that of asymptomatic or presymptomatic cases (5·4%), suggesting that the presence of symptoms is a host factor in driving transmission. To our knowledge, our study is the first to compare the transmissibility between these two groups across various settings.

The higher household AR observed among symptomatic cases suggest that testing for such cases should be prioritised, especially in low resource areas with limited testing capacity. Nonetheless, an AR of 5·4% in asymptomatic and presymptomatic cases is not negligible. As such, our findings have several implications for high resource areas with greater testing capacity. First, it strengthens the argument for testing household contacts in the absence of symptoms. Second, there is a need to allow for ‘slack’ in the surveillance system as the high proportion of asymptomatic cases pose challenges for rapid detection and isolation. We recommend that even in countries with highly developed testing and tracing capacities, moderate levels of social distancing should be implemented to account for this. Third, proactive testing of travellers and also attendees at ‘red flag’ events and institutional settings may be necessary to contain the spread of COVID-19.

This study has several limitations. First, we have not accounted for other potential environmental factors such as the relative size of the household, air ventilation, and indirect transmission through fomites. Second, we do not have information on NPIs practiced by the close contacts; presumably, individuals would take precautions during an outbreak. Third, symptom status of the cases were reported during their swab collection date. We assume this to be reflective of their actual condition when their close contacts were exposed, however, this may not be necessarily true for all cases.

The main strength of our study is the availability of a highly detailed and complete contact tracing dataset at the national level. Since all contacts were tested, it is reasonable to assume that this study more accurately detects SARS-CoV-2 transmission than those that only test symptomatic contacts. In conclusion, our analysis highlights the variability of SARS-CoV-2 transmission across different settings and in particular, the role of SSEs. We identify ‘red flags’ for the development of potential SSEs, and identify environmental, behavioural, and host factors that drive transmission. Overall, we provide evidence that a combination of case isolation, contact tracing, and moderate physical distancing measures is an effective approach to containment.