Hostname: page-component-78c5997874-j824f Total loading time: 0 Render date: 2024-11-10T08:56:18.379Z Has data issue: false hasContentIssue false

A network meta-analysis of secondary attack rates of COVID-19 in different contact environments

Published online by Cambridge University Press:  05 October 2021

Xunying Zhao
Affiliation:
Department of Epidemiology and Health Statistics, West China School of Public Health and West China Fourth Hospital, Sichuan University, Chengdu 610041, China
Ziqiong Shen
Affiliation:
Department of Epidemiology and Health Statistics, West China School of Public Health and West China Fourth Hospital, Sichuan University, Chengdu 610041, China
Litao Sun
Affiliation:
Centre for Medical Education, Cardiff University, Cardiff, UK
Long Cheng
Affiliation:
Department of Epidemiology and Health Statistics, West China School of Public Health and West China Fourth Hospital, Sichuan University, Chengdu 610041, China
Mengyuan Wang
Affiliation:
Department of Epidemiology and Health Statistics, West China School of Public Health and West China Fourth Hospital, Sichuan University, Chengdu 610041, China
Xiaofan Zhang
Affiliation:
Department of Epidemiology and Health Statistics, West China School of Public Health and West China Fourth Hospital, Sichuan University, Chengdu 610041, China
Bin Xu
Affiliation:
Department of Epidemiology and Health Statistics, West China School of Public Health and West China Fourth Hospital, Sichuan University, Chengdu 610041, China
Lulu Tian
Affiliation:
Department of Epidemiology and Health Statistics, West China School of Public Health and West China Fourth Hospital, Sichuan University, Chengdu 610041, China
Yunqi Miao
Affiliation:
Department of Epidemiology and Health Statistics, West China School of Public Health and West China Fourth Hospital, Sichuan University, Chengdu 610041, China
Xueyao Wu
Affiliation:
Department of Epidemiology and Health Statistics, West China School of Public Health and West China Fourth Hospital, Sichuan University, Chengdu 610041, China
Kun Zou*
Affiliation:
Department of Child, Adolescent and Maternal Health, West China School of Public Health and West China Fourth Hospital, Sichuan University, Chengdu, China West China Research Centre for Rural Health Development, Sichuan University, Chengdu, China HEOA Group, Institute for Healthy Cities, Sichuan University, Chengdu, China
Jiayuan Li*
Affiliation:
Department of Epidemiology and Health Statistics, West China School of Public Health and West China Fourth Hospital, Sichuan University, Chengdu 610041, China
*
Authors for correspondence: Kun Zou, E-mail: zoukun.jy@outlook.com; Jiayuan Li, E-mail: lijiayuan@scu.edu.cn
Authors for correspondence: Kun Zou, E-mail: zoukun.jy@outlook.com; Jiayuan Li, E-mail: lijiayuan@scu.edu.cn
Rights & Permissions [Opens in a new window]

Abstract

As the corona virus disease 2019 (COVID-19) pandemic continues around the world, understanding the transmission characteristics of COVID-19 is vital for prevention and control. We conducted the first study aiming to estimate and compare the relative risk of secondary attack rates (SARs) of COVID-19 in different contact environments. Until 26 July 2021, epidemiological studies and cluster epidemic reports of COVID-19 were retrieved from SCI, Embase, PubMed, CNKI, Wanfang and CBM in English and Chinese, respectively. Relative risks (RRs) were estimated in pairwise comparisons of SARs between different contact environments using the frequentist NMA framework, and the ranking of risks in these environments was calculated using the surface under the cumulative ranking curve (SUCRA). Subgroup analysis was performed by regions. Thirty-two studies with 68 260 participants were identified. Compared with meal or gathering, transportation (RR 10.55, 95% confidence interval (CI) 1.43–77.85), medical care (RR 11.68, 95% CI 1.58–86.61) and work or study places (RR 10.15, 95% CI 1.40–73.38) had lower risk ratios for SARs. Overall, the SUCRA rankings from the highest to the lowest were household (95.3%), meal or gathering (81.4%), public places (58.9%), daily conversation (50.1%), transportation (30.8%), medical care (18.2%) and work or study places (15.3%). Household SARs were significantly higher than other environments in the subgroup of mainland China and sensitive analysis without small sample studies (<100). In light of the risks, stratified personal protection and public health measures need to be in place accordingly, so as close contacts categorising and management.

Type
Review
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution and reproduction, provided the original article is properly cited.
Copyright
Copyright © The Author(s), 2021. Published by Cambridge University Press

Introduction

The world is still in the pandemic of corona virus disease 2019 (COVID-19) since its outbreak in December 2019 [1]. Although vaccination is available and implemented in many countries, public health measures are still essential for disease control and prevention. The transmission of COVID-19 is primarily through air and saliva during close contact [Reference Wiersinga2]. Long-term contact with infected persons and short-term contact with symptomatic patients are associated with high risks of transmission [Reference Chu3].

Understanding the secondary attack rates (SARs) among diverse contact environments is meaningful for making corresponding close contacts' management plans. Epidemiological studies have been conducted on the SARs in different contact environments during the COVID-19 epidemic. However, the reported SARs varied largely across original studies. For example, the reported household SAR ranged from 1.87% to 55.56% [Reference Hui-hui4, Reference Hui-Jun5]. Several systematic reviews and traditional meta-analyses have also been conducted. Some reviews reported the estimate of pooled SARs of COVID-19 in households, which varied from 16.6% to 21.1% [Reference Madewell6Reference Fung8]. Besides, Tian and Huo [Reference Tian and Huo9] and Thompson et al. [Reference Thompson10] also estimated SARs in other contact environments, such as social gathering settings, schools, workplaces, healthcare facilities and so on, with pooled SARs varying from 1% to 6%. These systematic reviews estimated the pooled SARs in different contact environments by subgroup analysis. However, pairwise comparisons and relative risks (RRs) of SARs among multiple contact environments have not been systematically estimated yet. Therefore, the understanding and evidence are lacking on whether there is a statistical difference between the SARs of two or more contact environments.

Network meta-analysis (NMA) can be an effective method to solve this problem. Compared with traditional meta-analysis, it can directly and indirectly compare multiple treatments and interventions, and more accurately evaluate the relative effects of different comparison groups to draw conclusions [Reference Andrea11]. NMA can not only compare the effects of treatments and interventions in pairs, but also get the relative ranking of these effects [Reference Shim12]. For example, Hutton et al. [Reference Hutton13] conducted an NMA on the risk of lymphoid tissue disorders and other cancers in kidney transplant patients by using the results of randomised controlled trials and observational studies. Most epidemiological studies of infectious diseases are observational studies. When a certain research problem in observational studies has multiple groups, the comparison of multi-arm groups can also be analysed using the principle of NMA.

Therefore, this study aimed to comprehensively compare the risk of SARs between different contact environments during the period of COVID-19 epidemic using the NMA.

Materials and methods

Data sources and searches

We searched English databases including PubMed, Embase and Web of Science, and Chinese databases including China National Knowledge Infrastructure (CNKI), Wanfang and China Biology Medicine disc (CBM) until 26 July 2021. Search terms included ‘novel corona*’ OR ‘new corona*’ OR ‘CoV’ OR ‘nCoV’ OR ‘COVID*’ OR ‘SARS-cov-2’ OR ‘SARS2’ combined with ‘Close contact*’ OR ‘family cluster’ OR ‘secondary attack rate’ OR ‘SAR’ OR ‘Transmission’ OR ‘contact trac*’. Subject words and free words (in title and abstract) were used in the search. The search strategy was first formulated in PubMed and then modified for other databases. The search strategies were presented in Supplementary Table S1. A hand search of the reference lists of included studies was made for more potentially eligible studies.

Study selection

To be included, the primary studies must be retrospective analyses of epidemiological characteristics and clustered epidemic reports. The contact follow-up data of COVID-19 patients and their close contacts must exceed 7–14 days, including the total number of COVID-19 patients, the number of follow-up cases and contact environments (at least two contact environments). The SARs of the close contacts must be reported in the articles or could be calculated by manual statistics (the number of positive cases of the close contacts/the total number of the close contacts). The following types of studies were excluded, including duplicate reports (in which case the study with the largest sample size was included), meeting abstract, expert consensus and suggestions.

The inclusion or exclusion of each study was conducted by three reviewers (X.-Y. Zhao, Z.-Q. Shen or L.-T. Sun). First, abstracts and titles were read to exclude the obvious irrelevant studies. Then, full texts for potentially eligible studies were obtained and read against the inclusion and exclusion criteria. Finally, cross-checking was made among the three reviewers. Disagreements were determined by group discussions.

Data extraction and quality assessment

A pre-designed data extraction table was used to extract data from studies including the first author, time of publication, study site, contact environments, the number of close contact management, the number of follow-up cases and the number of cases infected. The contact environments in each study were classified into seven well-defined categories.

The Joanna Briggs Institute (JBI) Critical Appraisal Tool for Case Series Study was used to evaluate the quality of studies [14]. Due to potential attrition bias introduced by the loss of follow-up, studies that did not conduct further tracking and polymerase chain reaction (PCR) testing of all identified close contacts were considered as low quality.

Data extraction and quality assessment were completed independently and cross-checked by two reviewers (X.-Y. Zhao, Z.-Q. Shen or L.-T. Sun). The classification of contact environments was defined according to the Guideline of COVID-19 Control and Prevention-Close Contact Management Procedure (7th ed.) and was adjusted by included studies [15, Reference Lin16]. Finally, seven categories of contact environments were included in our study: household settings (Hou), public places (Pub), meal or gathering settings (Mea), transportation (Tra), daily conversation (Dai), work or study places (Wor) and medical care (Med) (Supplementary Table S3). The contact environments were classified by two researchers (X.-Y. Zhao and Z.-Q. Shen) independently. If contact environments reported in the study cannot be classified into one of the above seven types, the data were discarded. Any disagreement was resolved by discussion to reach a consensus.

Definitions of measures

SAR is defined as the probability that close contacts of index cases become confirmed cases of a disease [Reference Liu17]. Relative risk (RR) is defined as the estimated relative risk between SARs of two different contact environments in pairs, calculated by the ratio of SARs between the two groups in the NMA [Reference Dias and Caldwell18] (the formula was as follows):

$$RR = \displaystyle{{SAR\;of\;contact\;setting\;A\;} \over {SAR\;of\;contact\;setting\;B}}$$

Data synthesis and analysis

Pair-wise RRs of the SARs between two contact environments were estimated using frequentist NMA [Reference Shim12]. First, the consistency test was used to test the global inconsistency across the network [Reference Watt19, Reference Higgins20]. If P < 0.05, the inconsistency model was selected, otherwise the consistency model was applied to calculate the RRs. The side-splitting (node-splitting) method was further used to compare whether there was a significant difference between the direct comparison and the indirect comparison of each RR in two contact environments [Reference Yu-Kang21]. The surface under the cumulative ranking curve (SUCRA) was calculated and the probability of being the highest SAR among SARs of all contact environments was sorted in the descending order [Reference Andrea11]. The Q test and I-squared (I 2) metrics were used to evaluate the heterogeneity of the studies of each pairwise comparison. And publication bias was examined by funnel plots and Egger's test [Reference Shim12]. Subgroup analysis had been planned to assess whether there were differences in the RRs and ranking of the SARs among different regions. For sensitivity analysis, small sample studies (<100) and studies with low quality were excluded.

Traditional meta-analysis was also conducted as an additional analysis to estimate the pooled SARs for each contact environment. Q test was used to test the heterogeneity among studies, and I 2 was applied to quantify the heterogeneity. If I 2 ≤ 50% or P > 0.1, the fixed-effect model was used, otherwise the random-effect model was selected [Reference Bo22]. Publication bias was examined by funnel plots and Egger's test [Reference Dawson23].

All analyses were performed using Stata 15.1 SE, employing network package [Reference Shim12] and metan package [Reference Harris24]. The significance level α was 0.05. The PROSPERO number of this study is CRD42020206576.

Results

Characteristics of included studies

A total of 2606 citations were retrieved from the literature search, including 857 citations from Chinese databases and 1749 from English databases. Among them, 1738 citations remained after removing duplicates. After screening the titles and abstracts, 1261 studies were excluded. The full texts of the remaining 328 studies were obtained and read. Finally, 296 studies were excluded by reading full text, and 32 studies were included (Fig. 1).

Fig. 1. Study selection process.

The characteristics of the included studies are presented in Supplementary Table S2. After removing types of contact environments that could not be classified, 32 studies involved 68 260 close contacts were included in the analysis, in which 2406 were infected. The research period of included studies was between January 2020 and March 2021. Among them, 21 studies were from mainland China, 2 were from the United States, 2 were from South Korea and the remaining 9 studies were from Australia, Singapore, Japan, Yemen, Rwanda, India and Taiwan, China, respectively.

Risk of bias in included studies

Among them, three studies reported the SARs of seven types of contact environments, 13 studies reported the SARs of three types of contact environments. Four studies have not conducted further tracking and PCR testing of all identified close contacts and were considered as low quality. The results of the study quality evaluation are presented in Supplementary File S4.

In 21 pairs of pairwise comparisons, 13 pairs of I 2 values exceed 50%. Heterogeneity of the pairwise comparisons between households and other contact environments was the highest (the values of I 2 were around 85%). The details of heterogeneity are presented in Supplementary Table S5. For all analyses, the Egger's test showed that there was no publication bias (P > 0.1). The results of the Egger's test and the funnel plots are presented in Supplementary File S9.

Outcomes of NMA for all included studies

Figure 2 shows the network geometry of eligible comparisons for SARs of the seven contact environments, with each contact environment represented by a node. The line between the nodes indicated evidence of a direct comparison between the two contact environments. All contact environments were included in the direct network comparison. Due to the result of the consistency test (P = 0.0022), the inconsistency model was applied for fitting the data. The results of side-splitting showed that the results of direct and indirect comparisons were consistent in each pairwise comparison (Supplementary File S6).

Fig. 2. NMA of SAR in seven contact environments.

Note: (A) All studies; (B) subgroup in mainland China, (C) studies without small sample sizes and (D) studies without low quality. Width of the line is directly proportional to the number of included studies. Size of the node is proportional to the sample size.

Table 2 shows the RRs of the NMA between the seven contact environments. The SAR of meal or gathering settings was nearly 10 times higher than that of transportation (RR 10.55, 95% confidence interval (CI) 1.43–77.85), work or study places (RR 11.68, 95% CI 1.58–86.61) and medical care (RR 10.15, 95% CI 1.40–73.38). As shown in Table 1, the SARs of the SUCRA ranked from the highest to lowest were: household settings (95.3%), meal or gathering settings (81.4%), public places (58.9%), daily conversation (50.1%), transportation (30.8%), medical care (18.2%) and work or study places (15.3%). Household SAR of SUCRA ranked the highest, which means household SAR had the greatest probability of having the highest SAR among all contact environments. But, due to the wide ranges of CIs, there were no statistically significant differences between the SAR of the household and other contact environments.

Table 1. Meta-analysis of SARs in different contact environments

SUCRA, the surface under the cumulative ranking curve; SAR, secondary attack rates.

a P for Q test.

For the results of the traditional meta-analysis, the pooled SAR of each contact environment is listed in Table 1. The pooled SARs was the highest for household SAR as 11.7% (95% CI 9.7%–13.8%), and the lowest one was medical care SAR as 0.7% (95% CI 0.4%–1.1%).

Outcomes of NMA for studies from mainland China

The network comparisons for SARs in the subgroup of mainland China are also shown in Figure 2. After excluding the 11 studies from regions other than mainland China, the data met the consistency assumption (P = 0.1401), and the consistency model was used for fitting.

For the 21 studies conducted in mainland China, the number of RRs of pairwise comparisons with statistical significance increased from three to eight. The SAR of households compared with other environments were all significantly higher (Table 2), in which RRs varied from 2.48 (95% CI 1.35–4.57) for meal or gathering to 9.71 (95% CI 4.70–20.05) for transportation. Besides, the SAR of meal or gathering settings was nearly three times higher than that of transportation and medical care (RR 3.91, 95% CI 1.85–8.27; RR 3.13, 95% CI 1.34–7.29). Compared with the SUCRA for all 32 studies, the relative rankings of SUCRA only exchanged between work or study places and transportation (Table 1).

Table 2. Pairwise comparisons of SARs for seven contact environments

Hou, household settings; Pub, public places; Mea, meal or gathering settings; Tra, transportation; Dai, daily conversation; Wor, work or study places; Med, medical care.

(A) The lower left corner shows the pairwise comparison results of all 32 studies, and the upper right corner shows the results of 21 mainland China. (B) The lower left corner shows the pairwise comparison results of studies without small sample sizes, and the upper right corner shows the results of studies without low quality.

The effect size represents the relative risk (RR, 95% CI) of the SAR of the contact environments on the left relative to the right. Bold and * indicate that the result is statistically significant.

Sensitivity analysis

Sensitivity analysis was conducted by excluding small sample studies (<100) and studies with low quality, respectively.

After removing the small sample studies, the P-value of the consistency test changed to 0.0759 and the consistency model was used for fitting. The results were quite different compared with the results of 32 studies, and there were 11 RRs of pairwise comparison showing significant differences. The RRs of household SAR compared with others were all significantly beyond 1, which was similar to the results of subgroup (Table 2). Among them, household SAR was 2.01 times that of meal or gathering (95% CI 1.12–3.60) and 11.04 times that of transportation (RR 11.04, 95% CI 5.53–22.04). Moreover, meal or gathering SAR was 2.56 times that of public places (RR 2.56, 95% CI 1.24–5.29), 3.37 times that of medical care (RR 3.37, 95% CI 1.55–7.33) and 5.50 times that of transportation (RR 5.50, 95% CI 2.58–11.72). When taking transportation as the reference group, the SAR of daily conversation and work or study places were close to three times that of transportation. In addition, the SUCRA rankings also have changed (Table 1). The daily conversation and the work or study places rose to the 3rd (59.4%) and the 4th (52%), and the public places and the transportation dropped to the 5th (35.1%) and the 7th (2.9%).

After excluding a low-quality study, the 4th to 7th SUCRA rankings changed, and the ranking of the work or study places rose to the 4th. However, the point estimations of the RRs were similar with the results of overall studies (Table 2).

Discussion

Main findings

This systematic review included currently the largest number of eligible observational studies involving 68 260 participants of SAR of COVID-19 in varied contact environments. The SARs of seven types of contact environments were compared and the pairwise RRs of them were estimated using the NMA approach for the first time. For all studies, the SAR of meal or gathering was significantly higher than the other three contact environments: Mea vs. Tra (RR 10.55, 95% CI 1.43–77.85), Mea vs. Wor (RR 11.68, 95% CI 1.58–86.61) and Mea vs. Med (RR 10.15, 95% CI 1.40–73.38). The SUCRA rankings of SARs from the highest to the lowest were household, meal or gathering, public places, daily conversation, transportation, medical care and work or study places. Although the SAR of the household settings was higher than other environments, their RRs were not statistically significant.

Our study suggested that the estimates of RRs from mainland China and studies without small sample sizes were more stable, supported by the significant RRs with narrower 95% CIs in NMA and larger P-value (>0.05) in consistency test. After excluding 11 non-mainland China studies and five studies with small sample sizes, the statistically significant RRs increased to 8 and 11, respectively. And the RRs between the household SARs and all other environments SARs became statistically significant in both analyses. Besides the household and meal or gathering settings, the SUCRA ranking of other contact environments also changed. These changes were possibly caused by the changes in SARs in diverse contact environments and the removing of small-study effect. The results of the subgroup analysis and sensitivity analysis showed that the SARs are significantly higher among multiple contact environments, especially for household and meal or gathering settings. Therefore, stratified personal protection and public health measures need to be in place accordingly, so as close contacts categorising and management.

Findings in our study could be explained by the characteristics of transmission, and were in line with previous studies. SARS-CoV-2 is mainly spread through respiratory droplets in face-to-face contact [Reference Wiersinga2]. Wearing a mask and keeping social distance can effectively reduce the possibility of transmission [Reference Wiersinga2]. People living together spend more time together and have a closer social distance, which may lead to higher SAR [Reference Tian and Huo9, Reference Koh25]. In healthcare environments, people use more protective measures, such as masks and disinfection, which reduces the probability of infection [Reference Chu3].

There were some potential confounders in our study, which may lead to inaccurate results, especially for the analysis of the overall studies. Multiple factors might influence the results, such as the population of the country, criteria for determining close contacts and lockdown status [Reference Shah7]. These regional factors could increase the heterogeneity of our study. Also, there were some factors considered as the confounder factors of SARs, such as asymptomatic detection rates and the status of herd immunity [Reference Madewell6, Reference Thompson10, Reference Koh25]. Some reviews suggested that the SAR of symptomatic cases was higher than asymptomatic cases because symptomatic cases release more viruses during the onset of symptoms, making close contacts more susceptible to infection [Reference Madewell6, Reference Thompson10, Reference Koh25]. Since the asymptomatic detection rates were almost equal in diverse contact environments in each study, the effects of the factor were partly offset and attenuated when calculating the RRs of SARs. For herd immunity, although whether the antibodies were good to prevent re-infection still needed to be confirmed, the status of herd immunity could be a potential confounder for SARs [Reference Poonia and Kottilil26]. But, when the status of herd immunity was similarly in the compared contact environments, the confounder effects may also be offset and attenuated.

The differences in regions

In our study, the results in mainland China were different from the overall studies. There may be several reasons. Most of the studies included in this study were from mainland China, which was possibly because mainland China has more detailed and unified management and classification of close contacts than many other countries. After excluding 11 studies from non-mainland China, results became consistent, estimates became significant and more accurate as the 95% CIs narrowed, indicating higher stability and reliability, and SUCRA ranking and the RRs of pairwise comparisons changed. These differences may be caused by different standards for determining close contacts and management programmes in China and other regions.

There were some studies in other countries that reported close contacts, but their management classification of exposure to the environments varied and often simplified, which mainly focused on household SAR [Reference Ng27Reference Wilkinson29]. China has conducted strict and detailed epidemiological information surveys of each close contact, and classified the contact environments into seven types according to national guidelines [15]. The World Health Organization's guidance for contact tracing classification was similar to the Chinese guideline, which has six categories including household contacts, closed settings, healthcare settings, professional contacts, public or shared transport and other well-defined settings [30]. The US Centre for Disease Control and Prevention provides risk assessment information and guidelines for close contacts, including school work environment, travel, community and medical care environment [31]. In Taiwan China, the Ministry of Health and Welfare suggested that contact environments be classified into four types, including hospital contacts, airplane contacts, school contacts and others [32]. The health departments of Singapore, South Korea and Australia did not have a clear classification of the settings of close contacts [33Reference Zastrow35]. Besides, the detection rate by contact environments between different regions may also influence the results. In China, almost all the close contacts were conducted further tracking and PCR testing, but the situation in other countries was different, such as Australia, India, Rwanda, Japan and United States, which might cause potential attrition bias introduced by the loss of follow-up [Reference Draper36Reference Sundar and Bhaskar40].

Meanwhile, in the early stage of the epidemic, the Chinese government conducted strict social isolation measures, which might lead to more contact frequency and opportunities in households than in other environments [Reference Wang and Wang41]. This might also be a potential cause of significant results when calculating the RRs between household SAR and other environments in mainland China.

Sensibility analysis of the study

In our study, the results excluding small sample studies were more stable and credible. More RRs showed significance, and their 95% CIs became much narrower. The risk of household SAR was significantly higher than that of other contact environments. The order of 3rd to 7th in the SUCRA ranking has also changed. Small sample studies may lead to bias because of potential publication bias [Reference Ioannidis42] and instability of effect value of SARs [Reference Pereira43]. Compared with large sample studies, the method of correcting zero events in the STATA software network package might cause a larger biased estimate for small sample studies, which may also cause NMA results including small sample studies to be more unstable [Reference Shim12]. Also, the inconsistency model might lead to a wider 95% CI than the consistency model, which may also contribute to the insignificant results [Reference Higgins20].

Besides the small-study effect, other biases of studies included in our studies were acceptable. Except for the higher risk of bias in clinical features (no clinical features of follow-up cases reported), which was not the focus of this study, the included studies have a lower risk of bias. Most studies might have a bias in one respect. After excluding low-quality studies, the 4th–7th SUCRA rankings changed, and the ranking of the work or study places group rose to 4th. But, the RRs were similar. Due to the decrease in the number of studies and participants included, the uncertainty of the results has increased and the SUCRA ranking has changed [Reference Wang and Carter44, Reference Jansen45]. It indicated the instability of the SUCRA ranking once again, and suggested that a conclusion should not be made on the results of SUCRA only, but with the effect size of the pairwise comparison [Reference Wang and Carter44].

Limitations of the study

This study has some limitations. First, due to the lack of specific quality evaluation tools for retrospective epidemiological studies on infectious diseases, after examining several risks of bias assessment tools, the JBI tool was selected to evaluate the biases, and an adaptive explanation was made for each question according to general principles of epidemiological studies [14]. But, this may still lead to certain errors or inaccuracies of quality evaluation. Therefore, to obtain more rigorous evidence through the NMA, methodological improvement is still needed.

Second, we planned to do subgroup analysis by region, age and other confounder factors. However, there were few studies reporting available demographic characteristics. Some only reported the demographic characteristics of the follow-up cases instead of all the close contacts. Some reported the demographic characteristics of overall close contact rather than that of each contact environment group. Due to the lack of information, we only carried out subgroup analysis by region eventually.

Finally, due to the lack of information on the classification of contact environments, we excluded some studies that could not be included in the seven classifications. Meanwhile, due to the method limitations of the NMA, we only included studies that reported at least two groups of exposure routes; in other words, studies that reported only one exposure route and total SAR were not included in our study, which might miss some information. And, NMA could not estimate the pooled SAR in each contact environment, so traditional meta-analysis was conducted to calculate the pooled SARs as Supplementary results.

Conclusion

Overall, household SAR has the probability of being the highest SAR among SARs of all contact environments, followed by SARs of meal or gathering, public places, daily conversation, transportation, medical care and work or study places. Household SARs were significantly higher than other environments supported by more accurate results using studies from China and sensitive analysis without small sample studies (<100). This supported the necessity of classifying and managing close contacts using stratified strategies in the COVID-19 epidemic.

Supplementary material

The supplementary material for this article can be found at https://doi.org/10.1017/S0950268821002223

Acknowledgements

This study was funded by the Sichuan Provincial Department of Science and Technology with grant number 2020YFS0559. The funding source did not involve in the research process. Any opinions in this study represent only the views of the authors and were not necessarily related to the Department of Science and Technology.

Conflict of interest

No conflict of interest exists in the submission of this manuscript.

Data availability statement

All the data were from the published studies. The dataset is available from the corresponding author by reasonable request (Jiayuan Li; E-mail: ).

Footnotes

*

The first two authors contributed equally to this work.

References

Listings of WHO's response to COVID-19. Available at https://www.who.int/news/item/29-06-2020-covidtimeline (Accessed 29 June 2020).Google Scholar
Wiersinga, WJ et al. (2020) Pathophysiology, transmission, diagnosis, and treatment of coronavirus disease 2019 (COVID-19): a review. JAMA 324, 782793. Published online: 11 July 2020.CrossRefGoogle ScholarPubMed
Chu, DK et al. (2020) Physical distancing, face masks, and eye protection to prevent person-to-person transmission of SARS-CoV-2 and COVID-19: a systematic review and meta-analysis. Lancet (London, England) 395, 19731987. Published online: 5 June 2020.CrossRefGoogle ScholarPubMed
Hui-hui, Z et al. (2020) Epidemic characteristics of close contacts of coronavirus disease 2019 in Xi'an. Journal of Xi'an Jiaotong University (Medical Sciences) 41, 502505.Google Scholar
Hui-Jun, L et al. (2020) Treatment and analysis of a clustering pneumonia caused by 2019-nCoV in Wuzhi county. Henan Journal of Preventive Medicine 31, 573576.Google Scholar
Madewell, ZJ et al. (2020) Household transmission of SARS-CoV-2: a systematic review and meta-analysis. JAMA Network Open 3, e2031756. Published online: 15 December 2020.CrossRefGoogle ScholarPubMed
Shah, K et al. (2020) Secondary attack rate of COVID-19 in household contacts: a systematic review. QJM: Monthly Journal of the Association of Physicians 113, 841850. Published online: 30 July 2020.CrossRefGoogle ScholarPubMed
Fung, HF et al. (2020) The household secondary attack rate of SARS-CoV-2: a rapid review. Clinical Infectious Diseases 73, S138S145.CrossRefGoogle Scholar
Tian, T and Huo, X (2020) Secondary attack rates of COVID-19 in diverse contact settings, a meta-analysis. The Journal of Infection in Developing Countries 14, 13611367.CrossRefGoogle ScholarPubMed
Thompson, HA et al. (2021) SARS-CoV-2 setting-specific transmission rates: a systematic review and meta-analysis. Clinical Infectious Diseases 73, e754e764.CrossRefGoogle ScholarPubMed
Andrea, C et al. (2013) Conceptual and technical challenges in network meta-analysis. Annals of Internal Medicine 159, 130137.Google Scholar
Shim, S et al. (2017) Network meta-analysis: application and practice using Stata. Epidemiology and Health 39, e2017047. Published online: 3 November 2017.CrossRefGoogle ScholarPubMed
Hutton, B et al. (2014) Checking whether there is an increased risk of post-transplant lymphoproliferative disorder and other cancers with specific modern immunosuppression regimens in renal transplantation: protocol for a network meta-analysis of randomized and observational studies. Systematic Reviews 3, 16. Published online: 25 February 2014.CrossRefGoogle ScholarPubMed
The Joanna Briggs Institute. Joanna Briggs institute reviewers' manual. In: Institute TJB, editor. Australia 2016.Google Scholar
The Guideline of COVID-19 (2020) Control and Prevention-Close Contact Management Procedure (7th ed).Google Scholar
Lin, Z et al. (2020) Incidence analysis of 1403 close contacts of corona virus disease 2019 patients in different contact modes. Journal of Shandong University (Health Sciences) 58, 5861, Published online: 4 August 2020.Google Scholar
Liu, Y et al. (2020) Secondary attack rate and superspreading events for SARS-CoV-2. Lancet (London, England) 395, e47. Published online: 3 March 2020.CrossRefGoogle Scholar
Dias, S and Caldwell, DM (2019) Network meta-analysis explained. Archives of Disease in Childhood-Fetal and Neonatal Edition 104, F8F12. Published online: 15 November 2018.CrossRefGoogle ScholarPubMed
Watt, J et al. (2019) Research techniques made simple: network meta-analysis. Journal of Investigative Dermatology 139, 412.e11. Published online: 24 December 2018.CrossRefGoogle ScholarPubMed
Higgins, JP et al. (2012) Consistency and inconsistency in network meta-analysis: concepts and models for multi-arm studies. Research Synthesis Methods 3, 98110. Published online: 1 June 2012.CrossRefGoogle ScholarPubMed
Yu-Kang, T (2016) Node-splitting generalized linear mixed models for evaluation of inconsistency in network meta-analysis. Value in Health 19, 957963. Published online: 19 December 2016.CrossRefGoogle ScholarPubMed
Bo, Z et al. Meta-Analysis of Rates and Software Implementation. doi: 10.7507/1672-2531.20140166CrossRefGoogle Scholar
Dawson, DV et al. (2016) Understanding and evaluating meta-analysis. Journal of the American Dental Association 147, 264270. Published online: 27 December 2015.CrossRefGoogle ScholarPubMed
Harris, RJ et al. (2008) Metan: fixed- and random-effects meta-analysis. The Stata Journal 8, 328.CrossRefGoogle Scholar
Koh, WC et al. (2020) What do we know about SARS-CoV-2 transmission? A systematic review and meta-analysis of the secondary attack rate and associated risk factors. PLoS One 15, e0240205.CrossRefGoogle ScholarPubMed
Poonia, B and Kottilil, S. (2020) Immune correlates of COVID-19 control. Frontiers in Immunology 11, 569611. Published online: 3 November 2020.CrossRefGoogle ScholarPubMed
Ng, OT et al. (2020) SARS-CoV-2 seroprevalence and transmission risk factors among high-risk close contacts: a retrospective cohort study. Lancet Infectious Diseases 21, 333343.CrossRefGoogle ScholarPubMed
Li, W et al. (2020) Characteristics of household transmission of COVID-19. Clinical Infectious Diseases 71, 19431946. Published online: 18 April 2020.CrossRefGoogle ScholarPubMed
Wilkinson, K et al. (2021) Secondary attack rate of COVID-19 in household contacts in the Winnipeg Health Region, Canada. Canadian Journal of Public Health Revue Canadienne de Santé Publique 112, 1216. Published online: 19 November 2020.CrossRefGoogle Scholar
Organization WH (2021) Contact tracing in the context of COVID-19 interim guidance.Google Scholar
Centers for Disease Control and Prevention. Available at https://www.cdc.gov/coronavirus/2019-ncov/index.html (Accessed 17 August 2021).Google Scholar
Ministry of Health and Welfare. Guidelines for investigation and contact tracing of severe special infectious pneumonia. Available at https://www.cdc.gov.tw/File/Get/B0h3BBa-bXGwwBwkxoUb1g (Accessed 22 January 2020).Google Scholar
What is home quarantine order?. Available at https://www.moh.gov.sg/covid-19/general/resources (Accessed 17 August 2021).Google Scholar
Australian Health Sector Emergency Response Plan for Novel Coronavirus (COVID-19). Available at https://www.health.gov.au/sites/default/files/documents/2020/02/australian-health-sector-emergency-response-plan-for-novel-coronavirus-covid-19_2.pdf (Accessed 7 February 2020).Google Scholar
Zastrow, M. How South Korea prevented a coronavirus disaster – and why the battle isn't over. Available at https://www.nationalgeographic.com/science/article/how-south-korea-prevented-coronavirus-disaster-why-battle-is-not-over (Accessed 12 May 2020).Google Scholar
Draper, AD et al. (2020) The first 2 months of COVID-19 contact tracing in the northern territory of Australia, March–April 2020. Communicable Diseases Intelligence (2018) 44. doi: http://dx.doi.org/10.33321/cdi.2020.44.53.Google Scholar
Burke, RM et al. (2020) Enhanced contact investigations for nine early travel-related cases of SARS-CoV-2 in the United States. PLoS One 15, e0238342. Published online: 2020/09/03.CrossRefGoogle ScholarPubMed
Akaishi, T et al. (2021) COVID-19 transmission in group living environments and households. Scientific Reports 11, 11616. Published online: 4 June 2021.CrossRefGoogle ScholarPubMed
Semakula, M et al. (2021) The secondary transmission pattern of COVID-19 based on contact tracing in Rwanda. BMJ Global Health 6, e004885.CrossRefGoogle ScholarPubMed
Sundar, V and Bhaskar, E (2021) Low secondary transmission rates of SARS-CoV-2 infection among contacts of construction laborers at open air environment. Germs 11, 128131. Published online: 27 April 2021.CrossRefGoogle ScholarPubMed
Wang, J and Wang, Z (2020) Strengths, Weaknesses, Opportunities and Threats (SWOT) analysis of China's prevention and control strategy for the COVID-19 epidemic. International Journal of Environmental Research and Public Health 17, 2235.CrossRefGoogle ScholarPubMed
Ioannidis, JP et al. (1998) Issues in comparisons between meta-analyses and large trials. JAMA 279, 10891093. Published online: 18 April 1998.CrossRefGoogle ScholarPubMed
Pereira, TV et al. (2012) Empirical evaluation of very large treatment effects of medical interventions. JAMA 308, 16761684. Published online: 25 October 2012.CrossRefGoogle ScholarPubMed
Wang, Z and Carter, RE (2018) Ranking of the most effective treatments for cardiovascular disease using SUCRA: is it as sweet as it appears? European Journal of Preventive Cardiology 25, 842843. Published online: 24 March 2018.CrossRefGoogle Scholar
Jansen, JP et al. (2014) Indirect treatment comparison/network meta-analysis study questionnaire to assess relevance and credibility to inform health care decision making: an ISPOR-AMCP-NPC good practice task force report. Value in Health 17, 157173.CrossRefGoogle ScholarPubMed
Figure 0

Fig. 1. Study selection process.

Figure 1

Fig. 2. NMA of SAR in seven contact environments.Note: (A) All studies; (B) subgroup in mainland China, (C) studies without small sample sizes and (D) studies without low quality. Width of the line is directly proportional to the number of included studies. Size of the node is proportional to the sample size.

Figure 2

Table 1. Meta-analysis of SARs in different contact environments

Figure 3

Table 2. Pairwise comparisons of SARs for seven contact environments

Supplementary material: File

Zhao et al. supplementary material

Zhao et al. supplementary material

Download Zhao et al. supplementary material(File)
File 1.8 MB