Broadening symptom criteria improves early case identification in SARS-CoV-2 contacts
- Hamish Houston1,9,
- Seran Hakki1,9,
- Timesh D Pillay1,
- Kieran Madon1,
- Nieves Derqui-Fernandez1,
- Aleksandra Koycheva1,
- Anika Singanayagam2,
- Joe Fenn1,
- Rhia Kundu1,
- Emily Conibear1,
- Robert Varro1,
- Jessica Cutajar1,
- Valerie Quinn1,
- Lulu Wang1,
- Janakan S Narean1,
- Mica R Tolosa-Wright1,
- Jack Barnett1,
- Onn Min Kon1,3,
- Richard Tedder4,
- Graham Taylor5,
- Maria Zambon2,
- Neil Ferguson6,
- Jake Dunning2,7,
- Jonathan J Deeks8,9 and
- Ajit Lalvani1,10⇑
- 1NIHR Health Protection Research Unit in Respiratory Infections, Imperial College London, London, UK
- 2National Infection Service, Public Health England, London, UK
- 3Tuberculosis Service, Imperial College Healthcare NHS Trust, London, UK
- 4Molecular Diagnostic's Unit, Imperial College London, London, UK
- 5Section of Virology, Department of Infectious Disease, Imperial College London, London, UK
- 6Department of Infectious Disease Epidemiology, Faculty of Medicine, Imperial College London, St Mary's Campus, Norfolk Place, London, UK
- 7NIHR Health Protection Research Unit in Emerging and Zoonotic Infections, University of Oxford, Oxford, UK
- 8Test Evaluation Research Group, Institute of Applied Health Research, University of Birmingham, Birmingham, UK
- 9HH & SH contributed equally
- 10JDe & AL contributed equally
- Corresponding author: Ajit Lalvani (a.lalvani{at}imperial.ac.uk)
Abstract
Introduction The success of case isolation and contact tracing for the control of SARS-CoV-2 transmission depends on the accuracy and speed of case identification. We assessed whether inclusion of additional symptoms alongside three canonical symptoms (CS) - fever; cough; loss or change in smell or taste – could improve case definitions and accelerate case identification in SARS-CoV-2 contacts.
Methods Two prospective longitudinal London-based cohorts of community SARS-CoV-2 contacts, recruited within 5 days of exposure, provided independent training and test datasets. Infected and uninfected contacts completed daily symptom diaries from the earliest possible time-points. Diagnostic information gained by adding symptoms to the CS was quantified using likelihood ratios and AUC-ROC. Improvements in sensitivity and time-to-detection were compared to penalties in terms of specificity and number-needed-to-test.
Results Of 529 contacts within two cohorts, 164 (31%) developed PCR-confirmed infection and 365 (69%) remained uninfected. In the training dataset (n=168), 29% of infected contacts did not report the CS. Four symptoms (sore throat, muscle aches, headache and appetite loss) were identified as early-predictors (EP) which added diagnostic value to the CS. The broadened symptom criterion “≥1 of the CS, or ≥2 of the EP” identified PCR-positive contacts in the test dataset on average 2 days earlier after exposure (p=0.07) than “≥1 of the CS”, with only modest reduction in specificity (5.7%).
Conclusions Broadening symptom criteria to include individuals with at least 2 of muscle aches, headache, appetite loss and sore throat identifies more infections and reduces time-to-detection, providing greater opportunities to prevent SARS-CoV-2 transmission.
INTRODUCTION
SARS-CoV-2 transmission is not entirely prevented by current vaccines. [1] As vaccination coverage increases, blanket isolation rules for SARS-CoV-2 contacts (i.e. test-trace-isolate) become less acceptable to society. Rapid identification and isolation of contacts who become infected is an increasingly important alternative strategy for prevention and containment. [2] Effectiveness depends crucially on how quickly such cases are detected and initiate self-isolation, [3] because individuals are most infectious early in the course of infection. [4, 5] However, case definitions must also be sufficiently specific to avoid overwhelming testing capacity. [6, 7]
There is considerable international heterogeneity in policy for COVID-19 community testing within the general population (S1). Most criteria include fever, cough, and loss or change in smell or taste (hereafter referred to as the canonical symptoms [CS]) alongside a range of other symptoms. Some countries are currently considering altering their case definitions. [13] There is thus an urgent need for empirical data to identify whether additional symptoms (we call early-predictors [EPs]) can augment the CS within community case definitions.
Surprisingly, empirical longitudinal data from recently exposed SARS-CoV-2 contacts are scarce. Recent large-scale cross-sectional studies of community testing data support adding more symptoms to the CS. [7, 14] Descriptive longitudinal retrospective studies of SARS-CoV-2 infections also exist. [15, 16] However, a high-resolution longitudinal evaluation of symptom combinations for differentiating infected SARS-CoV-2 contacts from exposed but uninfected controls has not, to our knowledge, been performed before.
Using data from two prospective longitudinal cohorts of SARS-CoV-2 contacts, we aimed to establish definitively whether broadening symptom criteria beyond the CS can accelerate and improve case detection without weakening specificity. Rapid recruitment following clearly defined exposure enabled optimal symptom criteria to be identified. Daily contemporaneously recorded symptom diaries ensured symptom onset times were recorded and time-savings measured with maximum precision. Through direct study of relevant community-based cohorts we provide generalisable evidence-based criteria for effective case definitions to rapidly identify and isolate infectious cases.
METHODS
Recruitment and study procedures
INSTINCT (Integrated Network for Surveillance, Trials and Investigations into COVID-19 Transmission) and ATACCC (Assessment of Transmission And Contagiousness of COVID-19 in Contacts) were two community-based cohort studies in which contacts of COVID-19 cases in Greater London were identified and recruited from 10th May 2020 through 31st March 2021.
Index cases, or contacts identified by the UK contact tracing system (NHS Test and Trace [NTAT]), were referred from Public Health England (PHE). Initially, referrals were also received from the Royal College of General Practitioners Research and Surveillance Centre (RCGP-RSC) network. Contacts referred within five days of their index case symptom onset (ISO), who provided valid informed consent were enrolled within our recruitment capacity until the end of the 2nd pandemic wave in the UK. Ethics approval was granted by the Health Research Authority (Research Ethics Committee reference: 20/NW/0231).
In INSTINCT, household-contacts living with their index cases were enrolled at home by research nurses (day 0) and visited again on days 7, 14 and 27. Date of ISO was recorded at enrolment and served as a proxy for exposure. Combined nose and throat swabs (CNTS) for RT-PCR testing and blood samples for serology were taken by research nurses at each visit and an additional CNTS by participants on day 4. Samples were processed at the Molecular Diagnostics Unit, Imperial College London. Antibody (IgM and IgG) to SARS-CoV-2 receptor binding domain (anti-RBD) was measured using a two-step double antigen binding assay (DABA) with recombinant S1 antigen on the solid-phase and labelled recombinant RBD as detector in the fluid-phase. [16] In ATACCC, household and non-household-contacts (i.e. not residing with their index) were enrolled. Dates of ISO (household-contacts) or exposure event (non-household-contacts) were provided by NTAT. After nurse-delivered training, participants self-sampled CNTS daily for 14 consecutive days. SARS-CoV-2 RT-PCR testing was performed at the Virus Reference Department, PHE Colindale.
At enrolment, demographic information was collected and participants recorded the onset date of prior symptoms. After enrolment, participants completed a daily symptom diary which assessed 20 symptoms (S2). Loss or change in smell or taste was recorded as one item (hereafter referred to as anosmia).
Definitions and reference standards
INSTINCT data were used as the training dataset. “Current infection” was set as the target condition and a rigorous composite reference standard was constructed to establish its presence or absence with maximum accuracy. [17] Contacts were assigned to the “infected” group if they were PCR-positive at day 0, 4 or 7. Contacts were assigned to the “uninfected” group if they were PCR-negative and had undetectable SARS-CoV-2 antibodies at all time-points. Participants were excluded if they had no serology results or were PCR-negative at all time-points but had detectable SARS-CoV-2 antibodies at study day 0, 7 or 27.
ATACCC data were used as the test dataset. In this cohort, daily PCR results were available but serological testing was not performed routinely. Contacts were assigned to the “PCR-positive” group if they had a positive PCR result by 7 days after enrolment and PCR-negative if all results were negative. Participants who became PCR-positive after study day 7 or had no PCR results were excluded from the analysis. Participants with only one positive PCR result with a high Ct value (>28) were excluded to minimise false-positives caused by recent rather than current infection.
In both cohorts, participants were made aware of their PCR results as they became available. Participants with missing ISO or exposure dates were excluded from analyses requiring these data. The study flowchart (figure 1) depicts participant numbers included in each analysis.
Statistical analyses
Time-to-event analysis
We used time-to-event analysis to describe the onset of COVID-19-related symptoms relative to ISO in INSTINCT (figure 1 cohort A). Briefly, we used symptoms reported by “uninfected” contacts to define baseline time-dependent hazards, and the difference between “infected” and “uninfected” contacts to define COVID-19-related hazards for each symptom (S3 for detailed methodology). Symptoms with a probability of occurring due to COVID-19 of >15% by 10 days post-ISO were selected as candidate symptoms for further evaluation.
Spiegelhalter knill-Jones models
We aimed to quantify any additional diagnostic value gained by adding each of the candidate symptoms to the CS using likelihood ratios (LRs) estimated for individual symptoms within combinations of symptoms. The Spiegelhalter Knill-Jones (SKJ) method was used rather than the independence Bayes approach in order to adjust for dependency caused by symptom co-occurrence. [18] This method is summarised in S4, having been described in detail previously. [18–22] Symptoms were considered as a series of binary tests based on their occurrence by each study day (e.g. fever by day 3 would be regarded as positive if fever had been reported on study day 2). Persistent cough and productive cough were combined into a single cough variable. We compared models using the CS to those with an additional symptom. AUC-ROC allows evaluation of model discrimination in training and test datasets. Candidate symptoms with useful LRs after adjustment for dependency with the CS and whose addition improved AUC-ROC across multiple early time-points were considered “early-predictors” (EP).
Evaluating simple case definitions
To assess real-world impact through readily applicable case definitions, each of the EP were added to the CS individually and together as a list requiring more than 1 to be positive by using the words “at least”. Diagnostic performance was assessed against the serial PCR reference standard in the test dataset (cohort D, figure 1) at each day post-exposure. We used time-to-event analysis to measure how quickly broadened case definitions would identify PCR-positive individuals and log-rank tests to make comparisons with the CS. Finally, we quantified the prevalence-dependent trade-off between true-positives and false-positives by calculating the NNT: the number of false-positives for every true-positive plus 1.
Software
Statistical analyses were performed in Stata version 17.0 (StataCorp, LLC, College Station TX) and R (R Core Team, 2020).
RESULTS
Patient cohorts
53011 referrals were received via three recruitment pathways (S5). Of 529 contacts within two cohorts, 164 (31%) developed PCR-confirmed infection and 365 (69%) remained uninfected. S6 shows demographic details for INSTINCT and ATACCC. Sex, ethnicity and BMI were similar between cohorts. Participants were slightly older in ATACCC than in INSTINCT (median 38 versus 34 years, p<0.001).
Participants were enrolled a median 3 days (IQR 2–4) post-ISO in INSTINCT, 4 days (IQR 4–5) post-ISO in household-contacts in ATACCC and 5 days (IQR 4–6) post-exposure-event in non-household-contacts. In INSTINCT 22/168 (13.1%) contacts were linked to index cases identified through the RCGP-RSC network with suspected, but not confirmed COVID-19 and were included to avoid selection bias. Symptom diaries were completed in INSTINCT for median 26 days (IQR 7–27) by “infected” and 7 days (IQR 0–27) by “uninfected” contacts, and in ATACCC for 13 days (IQR 0–20) by “PCR-positives” and 7 days (IQR 0–14) by “PCR-negatives”.
Four participants with confirmed infection required hospitalisation. Whilst more than 90% of “infected” contacts in INSTINCT (68/73, 93.1%) reported at least one of the 20 symptoms by day 7, over a quarter (21/73, 28.8%) did not report fever, cough or anosmia by day 7.
Sequence of onset of COVID-19 related symptoms
Time-to-event analysis of symptom onset following exposure (S7–8) showed that fever preceded anosmia and persistent cough preceded productive cough. Sore throat and rhinitis occurred early, and breathlessness later. Fatigue was commonly reported by “uninfected” contacts.
Thirteen symptoms had a probability of occurring due to COVID-19 of >15% by 10 days post-ISO (fever, persistent cough, productive cough, anosmia, headache, muscle aches, sore throat, rhinitis, appetite loss, breathlessness, diarrhoea, nausea and abdominal pain). 9 of these 13 symptoms are not included in the CS and were denoted candidate symptoms in further analyses. Other than the CS: rhinitis, sore throat, headache, muscle aches, and appetite loss had the largest cumulative COVID-19-related hazards.
Additional diagnostic value of candidate symptoms
Raw counts of participants who had reported each symptom by each study day in the training cohort are presented in S9. Used alone, cough, rhinitis, headache and muscle aches were the most sensitive symptoms whilst nausea and abdominal pain were insensitive (S10). Anosmia, fever, and appetite loss were highly specific symptoms.
The crude LRs (S11) show that any of the symptoms will affect post-test odds when they are used alone. However, when used in combination with other symptoms, their LRs after adjustment using the SKJ approach (figure 2, table 1, Table S12) were all less extreme than their crude LRs, indicating considerable dependency between symptoms.
When cough was used in combination with fever and anosmia it's adjusted LRs were closer to 1 than those of anosmia or fever (figure 2A, Figure S12). This was most likely due to the higher specificity of anosmia and fever (S10). When combined with the CS the presence or absence of nausea did not independently affect post-test odds, its adjusted LRs lying close to 1 (figure 2). Breathlessness was more common in the “infected” group. However, whilst breathlessness was reported without fever, anosmia or cough by “uninfected” contacts, this was rare in the “infected”, explaining why its adjusted positive LRs are below 1 and negative LRs above 1.
In training and test datasets, AUC-ROCs increased with study day reflecting improved discrimination afforded by greater accumulation of symptoms by later study days in “infected” contacts (figure 3). AUC-ROC was often greater in test data than training data, likely reflecting the longer median time to recruitment in ATACCC.
Between study days 0 to 3, the addition of headache, sore throat, muscle aches and appetite loss to the CS yielded the greatest improvements in AUC-ROC in the test dataset. When combined with the CS, appetite loss, headache, sore throat and muscle aches all consistently had positive adjusted LRs above 1 and negative adjusted LRs below 1 showing that both their presence and their absence added to the CS's ability to discriminate between the infected and uninfected. These symptoms were therefore considered “early-predictors” (EP).
Evaluating simple case definitions
Each of the four EP were combined individually with the CS using an “or” operator, as well as together using “or” and “at least” operators (Box 1).
Box 1 Construction of easily comprehensible case definitions using Boolean operators
Case definitions can contain Boolean logical operators such as: AND, OR and at least (see supplementary box 1 for examples in current use) | |
Symptom groups: | Referred to in text as: |
fever, cough, anosmia | Canonical Symptoms (CS) |
headache, sore throat, muscle aches, appetite loss | Early-predictors (EP) |
The “early-predictors” were each combined individually with the CS using an “OR” operator. | |
fever OR cough OR anosmia OR headache | CS or headache |
fever OR cough OR anosmia OR sore throat | CS or sore throat |
fever OR cough OR anosmia OR muscle aches | CS or muscle aches |
fever OR cough OR anosmia OR appetite loss | CS or appetite loss |
The “early-predictors were all combined together with the CS using “at least” and “OR” operators. | |
fever OR cough OR anosmia OR AT LEAST 1 OF headache, sore throat, muscle aches and appetite loss | ≥1 of the CS, or ≥1 of the EP |
fever OR cough OR anosmia OR AT LEAST 1 OF headache, sore throat, muscle aches and appetite loss | ≥1 of the CS, or ≥2 of the EP |
fever OR cough OR anosmia OR AT LEAST 1 OF headache, sore throat, muscle aches and appetite loss | ≥1 of the CS, or ≥3 of the EP |
fever OR cough OR anosmia OR AT LEAST 1 OF headache, sore throat, muscle aches and appetite loss | ≥1 of the CS, or ≥4 of the EP |
The addition of any symptom to the CS using an “or” operator increased sensitivity (figure 4A) whilst reducing specificity (figure 4B). The addition of appetite loss produced the smallest changes compared to the CS.
The CS identified 50% of PCR-positives by 6 days post exposure (figure 5, table 2, Table S14). Adding headache yields the greatest increase in sensitivity (figure 4A, Figure S13A) and would identify PCR-positives on average 2 days earlier (p=0.02), but causes the largest reduction in specificity (15.2% at 5 days post-exposure; figure 4B, Figure S13A). In contrast, CS or sore throat only reduced specificity by 5.7% at 5 days and identified PCR-positive cases earlier than the CS, by 1 day on average. This change was not statistically significant given the small number of PCR-positive participants in the cohort (p=0.1, n=91).
When all 4 EP are added to the CS, if all 4 are required, there is very little difference to the CS. In contrast, the case definition ≥1 of the CS, or ≥1 of the EP would increase sensitivity and identify PCR-positive cases a median 2 days earlier than the CS (p=0.002). However, the corresponding reduction in specificity by 5 days post-exposure (19.7%, figure 4B, Figure S13B) would lead 25% of PCR-negative individuals to be inappropriately identified (table 2, figure 5B). ≥1 of the CS, or ≥2 of the EP identified PCR-positive cases a median 2 days earlier than the CS (p=0.07) with a reduction in specificity of only 5.7% at 5 days post-exposure. This reduction is smaller than that caused by moving from the CS to various other international case definitions (S16). None of the EP were dispensable from this proposed criterion (S17).
The number of individuals identified in order to yield a single PCR-positive case, the NNT, increases rapidly immediately after exposure, reflecting an initial accumulation of false-positives because no one has yet developed symptoms actually caused by infection (figure 4C). NNT plateaus around 4–5 days following exposure, reflecting the incubation period. At 25.6% prevalence, ≥1 of the CS, or ≥2 of the EP had a NNT at 5 days post-exposure of 1.78 compared to 1.61 for the CS, indicating 17 additional individuals identified for every 100 infected individuals identified.
DISCUSSION
This, to our knowledge, is the first study to use daily symptom data prospectively collected from recently exposed infected and uninfected SARS-CoV-2 contacts to evaluate the diagnostic performance of symptom combinations for detecting infection.
Using this definitive study design, we found that 29% of individuals with PCR-confirmed COVID-19 did not report any of the CS, but 93% reported at least one symptom from a broader list of 20. We identified 4 EP symptoms (sore throat, headache, muscle aches, and appetite loss) providing additional early predictive power for identifying SARS-CoV-2-infected contacts. The case definition ≥1 of the CS, or ≥2 of the EP identified PCR-positive contacts 2 days earlier after exposure than the CS alone (p=0.07). This time-saving is critical given that shortening the delay from infectiousness to self-isolation from 2.6 to 1.2 days has been estimated to reduce transmission by 47%. [3] Moreover, the proportion of “symptomatic” infections and time-to-symptom onset are critical parameters in studies modelling effectiveness of testing and isolation strategies for contacts. [24]
Consistent with previous studies, headache and sore throat were sensitive symptoms, [7] which occurred early in the course of infection, [15] and were prevalent in our relatively young participants. [24] Importance of these symptoms will increase as vaccination of older age groups increases the proportion of infections occurring in the young. In agreement with the Real-time Assessment of Community Transmission-1 (REACT-1) study, we found that headache, muscle aches and appetite loss improved discrimination within statistical prediction models. [14] We add a crucial evaluation of readily applicable case definitions. We observed that both the structure of symptom criteria (e.g. use of the Boolean operator “at least”) and time-from-exposure had a considerable effect on diagnostic performance.
The SKJ approach enabled another important new observation. Although an important indicator of disease severity, [26] breathlessness was not a useful additional symptom for identifying early and mild infections because a hierarchy of symptoms exists. Breathlessness is unlikely to occur due to COVID-19 without prior fever, cough or anosmia and its inclusion reduces specificity.
Further strengths include day-by-day measurement of diagnostic performance following exposure and prospective data collection which mitigates recall bias. The rigorous reference standard employed in our training cohort maximised accuracy for the target condition and ensured only the most useful symptoms were taken forward to the test data. Neither serology nor PCR have 100% sensitivity for SARS-CoV-2 infection. [27] Using both serology and PCR at multiple time-points to define the absence of infection we minimised false-negatives. False-positive PCR results caused by recent rather than current infection were likely less common in our longitudinal study of recently exposed contacts than in studies involving random community-sampling. [14]
Limitations include modest sample size, largely White British population, minor differences between training and test cohorts and the potential for tick-box and behavioural biases. Study participants were usually highly motivated and attentiveness to mild symptoms (e.g. rhinitis) may have been increased by awareness of exposure, frequent study visits and co-residence with other participants. Contacts could not be blinded to their PCR results or those of their index.
Since we studied community-based COVID-19 contacts identified through NTAT, our findings are very likely generalizable. As large-scale cross-sectional data replicate our findings in smaller-scale daily-resolution longitudinal data, the combined evidence-base is now sufficient to influence policy. Broadening symptom criteria for use in the general population would likely identify more infections and reduce time-to-detection, reducing transmission. We propose that symptom criteria within case definitions to prompt symptomatic isolation and testing of SARS-CoV-2 contacts should include headache, sore throat, muscle aches and appetite loss as well as the canonical symptoms to optimise sensitivity. Two of these additional symptoms should be required to maximise specificity.
As highly vaccinated regions transition to lower COVID-19 incidence, investment in RT-PCR testing capacity will make such broader case definitions feasible. As societies develop alternatives to test-trace-isolate, application of evidence-based symptom criteria alongside judicious testing will be critical for early discrimination of infected and uninfected contacts. Accordingly, our findings should inform development of evidence-based national testing policies in many parts of the world now and in subsequent phases of the pandemic.
Acknowledgements
We would like to thank all the participants who were involved in the study and the support of our study research nurses Kristel Timcang, Sarah Hammett, Eimear McDermott, Constanta Luca, Sam Bremang Jada Samuel. We thank Myra McClure and Ellie Parker on behalf of the Molecular Diagnostic's Unit, Imperial College London. We thank Joanna Ellis and the Virus Reference Department Staff, and Lucy Mosscrop, Carolina Rosadas and the Molecular Diagnostic's Unit staff for providing SARS-CoV-2 PCR data. We thank PHE England staff for facilitating recruitment into the study. We thank Michael Whitfield, Mohammed Essoussi, Chitra Tejpal, Guilia Miserocchi, Harriet Catchpole and Anjeli Ketkar for conducting data entry and quality control, Holly Grey and Megan Davies for providing further logistical support and Charlotte Williams and Michael Whitfield for performing a literature search of the SARS-CoV-2 symptom case definitions used worldwide.
Footnotes
This article has supplementary material available from erj.ersjournals.com
Reporting standards: Results are presented in accordance with STARD (Standards for Reporting of Diagnostic Accuracy Studies) and STROBE (Strengthening Reporting of Observational Studies in Epidemiology) guidelines.
Funding Statement: This work is supported by a Department of Health and Social Care COVID-19 Fighting Fund award (ATACCC), by the NIHR Health Protection Research Unit in Respiratory Infections at Imperial College London in partnership with Public Health England (award NIHR200927) and by MRC/UKRI; nCoV: Serological detection of past SARS-CoV-2 infection by non-invasive sampling for field epidemiology and quantitative antibody detection. UK Research and Innovation; DOI: http://dx.doi.org/10.13039/100014013; Grant: CV 220-111; Department of Health and Social Care; DOI: http://dx.doi.org/10.13039/501100000276; Grant: COVID-19 Fighting Fund award (ATACCC); Public Health England; DOI: http://dx.doi.org/10.13039/501100002141; Grant: NIHR200927; National Institute for Health Research; DOI: http://dx.doi.org/10.13039/501100000272; Grant: NIHR200927; Medical Research Council; DOI: http://dx.doi.org/10.13039/501100000265; Grant: CV 220-111.
Author contributions: HH, SH & AL made substantial contribution to the conception of the work. TP, RK, RV, JF, HH, SH and AL contributed to study design and data collection methodology. HH, SH, AL, NF & JDe made substantial contribution to the design of the analyses. RV, MT-W, JB, VQ, JC, JDu, EC, LW, JN, JF, RK, SH, ND-F, AK, GT & MZ delivered and supported the administrative and technical study logistics. SH, ND-F, AK, & KM carried out the initial data acquisition. RT provided the serology data. GT & MZ provided the RT-PCR data. HH, KM & JDe analysed the data. Statistical analyses were checked by JDe. HH, SH, KM, JDe, AS, JDu, OMK & AL contributed to data interpretation. HH & SH drafted the manuscript. AL was the Principal Investigator for INSTINCT. AL and JDu were joint Principal Investigators for ATACCC. All authors contributed to revising the manuscript critically for important intellectual content, approved the final manuscript and are accountable for all aspects of the work.
Conflict of interest: No competing interests were declared by any of the study authors.
- Received August 22, 2021.
- Accepted November 11, 2021.
- Copyright ©The authors 2021
This version is distributed under the terms of the Creative Commons Attribution Licence 4.0.