Skip to main content

Nasopharyngeal angiotensin converting enzyme 2 (ACE2) expression as a risk-factor for SARS-CoV-2 transmission in concurrent hospital associated outbreaks

Abstract

Background

Widespread human-to-human transmission of the severe acute respiratory syndrome coronavirus two (SARS-CoV-2) stems from a strong affinity for the cellular receptor angiotensin converting enzyme two (ACE2). We investigate the relationship between a patient’s nasopharyngeal ACE2 transcription and secondary transmission within a series of concurrent hospital associated SARS-CoV-2 outbreaks in British Columbia, Canada.

Methods

Epidemiological case data from the outbreak investigations was merged with public health laboratory records and viral lineage calls, from whole genome sequencing, to reconstruct the concurrent outbreaks using infection tracing transmission network analysis. ACE2 transcription and RNA viral load were measured by quantitative real-time polymerase chain reaction. The transmission network was resolved to calculate the number of potential secondary cases. Bivariate and multivariable analyses using Poisson and Negative Binomial regression models was performed to estimate the association between ACE2 transcription the number of SARS-CoV-2 secondary cases.

Results

The infection tracing transmission network provided n = 76 potential transmission events across n = 103 cases. Bivariate comparisons found that on average ACE2 transcription did not differ between patients and healthcare workers (P = 0.86). High ACE2 transcription was observed in 98.6% of transmission events, either the primary or secondary case had above average ACE2. Multivariable analysis found that the association between ACE2 transcription (log2 fold-change) and the number of secondary transmission events differs between patients and healthcare workers. In health care workers Negative Binomial regression estimated that a one-unit change in ACE2 transcription decreases the number of secondary cases (β = -0.132 (95%CI: -0.255 to -0.0181) adjusting for RNA viral load. Conversely, in patients a one-unit change in ACE2 transcription increases the number of secondary cases (β = 0.187 (95% CI: 0.0101 to 0.370) adjusting for RNA viral load. Sensitivity analysis found no significant relationship between ACE2 and secondary transmission in health care workers and confirmed the positive association among patients.

Conclusion

Our study suggests that ACE2 transcription has a positive association with SARS-CoV-2 secondary transmission in admitted inpatients, but not health care workers in concurrent hospital associated outbreaks, and it should be further investigated as a risk-factor for viral transmission.

Peer Review reports

Introduction

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) causes potentially life threatening lower respiratory and systemic inflammatory disease defined as COVID-19 [1,2,3]. SARS-CoV-2 has spread widely since late 2019 causing a global pandemic. Two previous public health emergencies have provided the opportunity to study human-to-human transmission of pathogenic coronaviruses. In 2003, the severe acute respiratory syndrome emerged infecting eight-thousand four-hundred and twenty-two people [4]. In 2012 and 2015 the Middle East respiratory syndrome coronavirus (MERS-CoV) caused epidemics in the Kingdom of Saudi Arabia and South Korea with two-thousand five hundred and sixty-two laboratory confirmed cases reported to the World Health Organization [5]. SARS- and MERS-CoV spread predominately in health care settings and community spread was controlled by public health interventions or self-limited [4, 5]. The transmission pattern of SARS-CoV-2 contrasts that of the other pathogenic coronaviruses (SARS-CoV and MERS-CoV). SARS-CoV-2 has predominately spread within the community, the virus has evolved to become more infectious over time and vaccination does not protect against re-infection or transmission of newer strains [6]. Transmission of respiratory viruses like SARS-CoV-2 involves a complex interplay of social, environmental, and biological variables [7]. Numerous observational studies have described risk factors of SARS-CoV-2 transmission across a variety of settings including hospitals [8], households [9], and schools [10, 11]. Cumulative evidence suggests that the increased human-to-human transmissibility of SARS-CoV-2 in comparison to SARS-CoV-1 stems from stronger affinity for the primary host receptor angiotensin converting enzyme 2 (ACE2) [12,13,14]. SARS-CoV-2 binds to ACE2 through its spike glycoprotein, the spike protein attaches to ACE2, and undergoes proteolytic cleavage via transmembrane serine protease 2 (TMPRSS2) or a similar host protease prompting endocytosis and membrane fusion [13]. The Omicron variant of SARS-CoV-2 has a stronger affinity for ACE2 and less reliance on host proteases biasing cell entry towards endocytosis, consequently less cell fusion and syncytia formation occur [15]. In addition to ACE2, neuropilin-1 (NRP1) has also been identified as a possible host receptor for SARS-CoV-2. Enhancement of SARS-CoV-2 infectivity by NRP1 stems from increased viral entry and may occur in both the presence of ACE2 and TMPRSS2 [16, 17]. ACE2 is highly expressed in the upper and lower respiratory tract; upregulation may occur in response to stimulation with interferon [18]. SARS-CoV-2 exploits upregulation of ACE2 in an interferon rich environment (like the early stages of an innate immune response) to spread from cell-to-cell producing a high viral titre [19]. In the later stages of replication, SARS-CoV-2 downregulates ACE2, preventing any benefits to the host mediated by interferon stimulated upregulation [20, 21]. SARS-CoV-2 has a shifting relationship with ACE2 and like viral load, ACE2 expression fluctuates through the course of viral infection [18, 22, 23]. Regardless, SARS-CoV-2 replicates most efficiently in tissues rich with ACE2, the strong positive association between viral load and ACE2 expression at the time of diagnostic testing increase the degree of viral shedding and thereby, the risk of ongoing transmission [24].

To understand the role of ACE2 expression in SARS-CoV-2 transmission, we conducted an infection tracing network analysis of positive testers in concurrent, single-site, hospital associated COVID-19 outbreaks in British Columbia during late 2020 and early 2021. The study aims to: (i) describe the outbreaks in context of person, place, and time, (ii) visualize an infection tracing transmission network to infer the number of secondary cases per primary case and (iii) quantify the relationship between the number of potential secondary cases (outcome) and nasopharyngeal transcription of transmembrane ACE2 (exposure) while adjusting for viral RNA load, and interaction by case designation (patient or healthcare worker). We build from previous work which demonstrated no relationship between transmembrane nasopharyngeal ACE2 transcription and age, biological sex or TMPRSS2 transcription in a sample of COVID-19 negative persons [24]. This study contributes novel evidence of how nasopharyngeal ACE2 transcription may drive SARS-CoV-2 transmission and highlights the potential role of respiratory masks and infection prevention and control measures in limiting the viruses’ nosocomial spread.

Methods

Study design and participants

We conducted surveillance for COVID-19 cases in a British Columbian hospital over a series of outbreaks, from declaration of the first facility associated outbreak on 07-11-2020 to the date the last outbreak was declared over on 04-01-2021 (dd-mm-yyyy). Participants included patients and health care workers who tested positive for COVID-19 upon self-reporting or showing indicative symptoms (e.g. fever, fatigue, cough and loss of taste or smell) or through asymptomatic point prevalence testing during the outbreak investigations. Inclusion criteria were applied to select study participants whose diagnostic specimens were tested centrally at the British Columbia Center for Disease Control Public Health Laboratory, underwent SARS-CoV-2 whole genome sequencing, and had adequate remaining volume of nucleic acid extract to assay ACE2 transcription. In the case that a participant was tested for COVID-19 more than once, their first diagnostic specimen from the study period was sampled. Participants who met the inclusion criteria (n = 202) were excluded from the study if their specimen collection container was not identifiable (n = 22), whole genome sequencing failed and was not able to classify the viral lineage (n = 63) or ACE2 transcription was unmeasurable (n = 14) (Figure S1). Some collection containers were not identifiable because of damage to their labels during storage or there was discordance between the line list and labels, due to clerical error. An analytic dataset of n = 103 participants was used for analysis. Demographic variables of age, biological sex, case description (patient or health care worker) and site of infection/ site of transmission (hospital unit) were drawn from public health laboratory data or the outbreak report of the investigating epidemiologist. Outbreaks were defined by the presence of one or more confirmed cases of COVID-19 which were epidemiologically linked within hospital units (Table S3). The laboratory procedures were performed in a laboratory accredited by the College of American Pathologists and BC’s Diagnostic Accreditation Program using validated RT-qPCR and whole genome sequencing protocols [25,26,27].

Ethics statement

Ethical approval for the study was sought from the University of British Columbia human ethics board (#H20-01110) which was harmonized with the Fraser Health Authority. Written informed consent was not required. All data was de-identified prior to analysis and the results were not linked back to any identifying records. This study was deemed as minimal risk to the participants involved. To ensure privacy the site of the outbreak series will not be disclosed.

Procedures

Participant’s diagnostic specimens were collected by nasopharyngeal swab and stored in Universal Transport Medium™ (UTM®) (COPAN®) and stored at 4°C before RNA extraction. RNA extraction was performed on the MagMAX-96™ platform with the Viral RNA isolation kit (ThermoFisher Scientific). Host and viral gene transcription was assayed by quantitative real-time polymerase chain reaction (RT-qPCR) on the Applied Biosystems 7500 Real-Time PCR platform using TaqMan™ FastVirus 1-step polymerase (ThermoFisher Scientific). Reaction volumes were 20 μl, with 5 μl of RNA template, 1 μl of 20X primer/probe, 5 μl FastVirus and 9 μl of nuclease free water. Cycling conditions were set to: 50°C for 5 min, 95°C for 20s followed by 40 cycles of 95°C for 15s and 60°C for 1 min. As previously described, multiplex RT-qPCR reactions were used to detect host (ACE2, GAPDH, RNaseP) and viral (E gene) transcription [28]. Participants were diagnosed COVID-19 positive with an E gene cycle threshold value < 38 and a RnaseP gene cycle threshold value < 40. The E gene Ct values were transformed to genome equivalents per millilitre using a standard curve of SARS-CoV-2 synthetic RNA (MN908947.3) (Twist Bioscience) [24]. Relative gene transcription for ACE2 was calculated in proportion to the control gene GAPDH using the 2−ΔΔCt method, reported in log2 fold-change [29]. SARS-CoV-2 whole genome sequencing was performed on all diagnostic specimens which tested positive for SARS-CoV-2 by RT-qPCR. Viral genomes were amplified using a 1200 bp amplicon scheme and sequenced on an Illumina NextSeq instrument [25]. Genome assembly was performed via a modified ARTIC Nextflow pipeline [30]. Quality control did not pass any genomes with < 85% completeness or < 10X depth of coverage. Viral lineages were assigned using the PANGOLIN tool (Version 1.15.1) [31]. All molecular and genomic testing was performed at the British Columbia Centre for Disease Control Public Health Laboratory (BCCDC-PHL).

Variable definition

Demographic variables of interest were drawn from public health laboratory data, or the outbreak reports provided by the investigating epidemiologist. The measures included in the study are age, biological sex, case description, diagnostic specimen collection date (dd-mm-yyyy), SARS-CoV-2 E gene cycle-threshold (Ct) value, ACE2 gene transcription, SARS-CoV-2 lineage (PANGO lineage), site of transmission (hospital unit) and site of infection (hospital unit). Age, E gene Ct value (transformed to log10 GE/ml [24]) and ACE2 gene transcription are continuous numeric variables. Biological sex, case description, viral lineage, collection date and site of transmission/infections are categorical variables. We assumed that collection date equals date of symptom onset for this study because, once the hospital declared an outbreak, patients and health care workers were required to be tested at the onset of symptoms. ACE2 gene transcription was transformed to a categorical value using the mean transcription value (X̄ = 0.00, SD = 1.08) of COVID-19 negative testers, from a previous study which collected a random sample of n = 212 nasopharyngeal specimens in British Columbia during 2020 [25]. The number of potential secondary cases per primary case was determined by infection tracing transmission network analysis. ACE2 gene transcription and number of potential secondary cases are the exposure and outcome of interest. Viral load and case description were included in multivariable analysis because, they meet the definition of a confounder and are a common cause of the exposure (ACE2 gene transcription) and outcome (number of secondary cases) [33]. Viral load and ACE2 gene transcription share a strong positive association, people with high ACE2 gene transcription were found to have high viral loads [34]. Viral load relates to transmission in that people with high viral loads may shed more infectious virus [35, 36]. Case description (patient or healthcare worker) relates to ACE2 gene transcription because patients may have comorbidities or receive medications which affect ACE2 transcription [37]. Factors that affect differential upper respiratory tract ACE2 expression are not well understood causing unmeasured confounding [38]. Health care workers have additional social connections relative to in-patients, as they can leave work and be exposed to SARS-CoV-2 outside of clinical areas or in the community [8]. Viral lineage was not considered a confounder or effect modifier in the transmission model as none of the SARS-CoV-2 variants circulating at the time of the study were deemed variants of concern, implying a possible phenotypic difference [39].

Statistical analysis

Descriptive statistics

The analytic data (n = 103) was used for bivariate analysis. The variables of interest were stratified by case description and categorical ACE2 transcription. An epidemiological curve was visualized for the series of concurrent hospital outbreaks from 07-11-2020 to 04-01-2021 (dd-mm-yyyy). Parametric statistical tests were used given the large sample size of our study [39].

Infection tracing transmission network analysis

A matrix of all possible transmission pairs was permuted using SARS-CoV-2 lineage classifications, for two cases to form a pair they had to share the same lineage. Transmission networks are degenerate meaning that a single case may form one, multiple or no pairs with secondary cases. As such the total number of cases and total number of transmission pairs are not additive but factorial (e.g., 4 cases could produce 24 transmission pairs) [40]. SARS-CoV-2 lineages circulating in the hospital were diverse enough that the minimal difference of single nucleotide polymorphisms (SNPs) between them was unlikely to have occurred by mutation during a single transmission event, suggestion multiple independent introductions (Table S4). Four assumptions (i-iv) were used to select potential transmission pairs from all possible permutations. The first assumption (i) stipulated that specimen collection date (symptom onset date) of the primary case was before that of the secondary case. The second assumption (ii) affirms that the primary case’s unit of transmission equals the secondary case’s unit of infection. The third assumption (ii) holds that the collection date of the secondary case is at least one-serial interval (5 days) from that of the primary case [41]. The fourth assumption (iv) dictates that the collection date of the secondary case is not more than three-serial intervals (15 days) from the primary case [42]. Potential transmission pairs, which met the four-assumptions, were stratified by categorical ACE2 transcription (High/Low) or viral lineage and plotted overtime.

Primary analysis

Multivariable analysis was performed to estimate the relationship between ACE2 transcription and number of secondary cases using a Poisson generalized linear regression model. Variable importance was assessed conceptually using the common cause criterion and statistically by the partial F-test. Collinearity was measured using the variable inflation factor, variables with a value greater than five were excluded [43]. Effect modification terms were included if they were statistically significant and supported conceptually. The two assumptions of Poisson models; overdispersion and zero inflation were checked using a ratio of residual deviance to degrees of freedom (theta) [44] and the score test [45], respectively. If the specified Poisson model failed to meet either of these assumptions than an alternative quasi-Poisson or Negative Binomial (NB) model was tested, and the assumptions re-examined. Model fit was measured by the Akaike information criterion [46].

Sensitivity analysis

Sensitivity analysis was performed to measure the impact of assumptions (ii, iii and iv) used to construct the transmission network on the relationship between ACE2 transcription and number of secondary cases. The assumptions were excluded independently, and the transmission network was iteratively reconstructed with (n-1) assumptions. Additionally, the fourth assumption (iv) was challenged by changing the infectious period from 15 days (three serial intervals) to 10 days (two serial intervals). The number of secondary cases was computed for each transmission network and used in multivariable analysis.

Statistical analysis was performed in R version 4.04 using the packages: readxl, tidyverse, dataexplorer, ggpubr, car, ggsci, stringr, tableone, rio, remotes, lubridate, dplyr, epicontacts, AER, devtools, rlang, DHARMa, MASS, pscl, epiR, EpiCurve [47].

Results

Descriptive statistics

The analytic dataset contains (n = 103) cases of COVID-19 associated with a single-site concurrent series of hospital outbreaks from 07-11-2020 to 04-01-2021 (dd-mm-yyyy) in British Columbia, Canada. Participant characteristics are displayed in Table 1, where the variables of interest are stratified by case description (Patients or Health Care Workers) (Table 1). The bivariate relationship between variables and categorical ACE2 transcription was also investigated (Table S1). Most cases occurred in patients n = 57 (55%), while n = 48 infections were observed in health care workers (45%). The mean age of patients associated with the hospital outbreaks was 78.61 years, health care workers were significantly younger with a mean age of 41.11 years (P < 0.001). Most of the healthcare workers (n = 41, 89%) and patients were biologic females (n = 33, 58%). ACE2 transcription and viral load did not differ between case descriptions (P = 0.86, P = 0.30). Six SARS-CoV-2 lineages among cases were characterized by whole genome sequencing, two of which did not have more than one case (B.1.128 and B.1.36.38). The highest proportion of observed SARS-CoV-2 infections were either SARS-CoV-2 lineage AL.1 (n = 56) or B.1.2 (n = 40) (Table 1), these two lineages have at least a five SNP difference between them (Table S4).

Table 1 Characteristics of analytic data stratified by case designation (n = 103)

An epidemiological curve was made to show the incidence of positive testers per week over the course of the hospital associate outbreaks. The outbreaks spanned ten epidemiological weeks (Week-#45-2020 to Week-#1-2021) and possesses a non-normal distribution (Fig. 1A & B). We are cautious to directly interpret the incidence curve as infection prevention and control practices implemented during the outbreaks (like point prevalence testing and stoppage of admissions) likely biases the observed incidence of SARS-CoV-2 cases overtime. A first peak in cases occurred on Week #46-2020 and a second on Week #52-2020. In the first phase of the outbreaks, many cases occurred on the second floor of the hospital in units 2A, 2B and 2C. Later transmission was predominately observed on floors four or five in units 4B, 4C, 4D and 5A (Fig. 1A). Multiple introductions of SARS-CoV-2 into the hospital environment occurred over the surveillance period; however, timely and effective infection prevention and control measures limited transmission of several variants like B.1.36.38, B.1.128 and B.1.279 (Table S3, Figure S2).

Fig. 1
figure 1

Incidence curve of several hospital associated SARS-CoV-2 outbreaks in British Columbia by epidemiological week from 07-11-2020 to 04-01-2021 (n = 103 laboratory confirmed cases). SARS-CoV-2 incidence follows a non-normal distribution where an initial peak caseload was observed in W46-2020 and a second lesser case surge occurred in W52-2020. A) SARS-CoV-2 cases are stratified by hospital floor, early in the surveillance period transmission occurred mostly on floor number two and later moved to floors four and five. If the floor where transmission occurred was not determined by the investigating epidemiologist, then it was coded as ‘unknown’ for our analysis. B) SARS-CoV-2 cases are stratified by viral variant, n = 7 unique SARS-CoV-2 variants were identified by whole genome sequencing participants diagnostic specimens. SARS-CoV-2 variants AL.1 and B.1.2 caused the highest percentage of cases: AL.1, 55/103, 53% and B.1.2, 39/103, 38%

Infection tracing transmission network analysis

Permuting all possible transmission pairs in which the primary and secondary case share the same SARS-CoV-2 lineage yielded n = 9389 combinations. Applying the first assumption (i) filtered the possible transmission pairs to n = 2167. The second assumption (ii) eliminated another n = 1480 possible pairs for a total of n = 687. The third and fourth assumptions restricted the possible pairs to n = 382 and n = 76, respectively. A directed transmission network was constructed from n = 103 cases with n = 76 possible transmission pairs. The n = 103 cases include n = 55 isolates of

viral lineage AL.1 and n = 39 for B.1.2; case counts of the other viral lineages was either too low (only one case reported) or did not meet the assumption criteria (n = 9) (Fig. 2B). The nodes of the transmission network are shown stratified by categorical ACE2 transcription level (Low or High) (Fig. 2A). The distribution of ACE2 transcription among hospital associated cases was tested with a single sample proportion test, 82% of cases had above average ACE2 transcription (84/103, P < 0.001). Transmission pairs (n = 76) were further stratified into primary (n = 34) and secondary cases (n = 28). Primary cases possessed predominately high ACE2 transcription (29/34, P < 0.001). ACE2 transcription was similarly enriched in secondary cases (23/28, P < 0.01). Overall, 98.6% of transmission events involved at least one case with high ACE2 transcription (Table S2).

Fig. 2
figure 2

Infection tracing transmission network analysis of a series of hospital associated SARS-CoV-2 outbreaks in British Columbia by epidemiological week from 07-11-2020 to 04-01-2021 (dd-mm-yyyy) (n = 103 laboratory confirmed cases). Nodes represent confirmed cases while the edges show the direction of outgoing transmission, link a primary case with at least one secondary case. (A) Nodes in the transmission network are stratified by ACE2 transcription (High or Low), 82% of cases (84/103) had above average, high ACE2 transcription (P < 0.001). (B) Nodes in the transmission network are classified by SARS-CoV-2 viral variant. The viral variants AL.1 and B.1.2 caused multiple secondary infections thorough the outbreaks, in the first phase (before Dec.15th, 2020) of the study period most cases were AL.1 and later transitioned to B.1.2 in phase two (after Dec.15th, 2020)

Primary analysis

Multivariable NB regression estimated that a one-unit change in ACE2 transcription decreases the number of secondary cases in health care workers by β = -0.13 (95%CI: -0.255 to -0.0181) adjusting for RNA viral load and case description (Table 2). Effect modification was observed between ACE2 transcription and case description; in patients, a one-unit change in ACE2 increased the number of secondary cases by β = 0.187 (95% CI: 0.0101 to 0.370). Therefore, a 5 unit increase in ACE2 transcription (log2 fold change) may lead to approximately more than one (R0 = 1.21) SARS-CoV-2 secondary cases per primary patient case. Poisson and quasi-Poisson regression models were constructed using the same variables, neither provided a better fit to the NB option (Table 2).

Table 2 Multivariable analysis of the relationship between ACE2 transcription and the number of SARS-CoV-2 secondary cases

Sensitivity analysis

A sensitivity analysis was performed to measure the impact of assumptions ii, iii and iv on construction of the transmission network and number of secondary cases. The first assumption- that the symptom onset date of the primary case occurs before the secondary case-was not included in the sensitivity analysis, as breaking this assumption would make the secondary case the primary case. Omission of the second assumption (ii), so that the primary case did not share a hospital unit with the secondary case, increased the number of possible transmission pairs to n = 538 (Figure S3). The network was not clearly resolved and no significant relationship between ACE2 transcription and the number of secondary cases was estimated using a Poisson (β = -0.17 [95%CI: -0.903 to 0.558], P = 0.64) or NB regression model ( β = -0.01 [95%CI: -0.082 to 0.0613], P = 0.8) (Fig. 3). Both models had a worse fit then when used in the primary analysis (PoissonAIC = 705 and NBAIC = 691). Leaving out the third assumption (iii), to appreciate that the secondary case could have been infected in less time than the average serial interval, increased the number of possible transmission pairs to n = 246 (Figure S4). A significant relationship between ACE2 transcription and the number of secondary cases was estimated using a Poisson ( β = -0.60 [95%CI: -1.08 to -0.125], P = 0.015) and NB regression model ( β = -0.11 [95%CI: -0.215 to -0.0189], P = 0.03) (Fig. 3). However, neither model fit better then when used in the primary analysis (PoissonAIC = 618 and NBAIC = 551). Re-analysis without the fourth assumption that the infectious period of SARS-CoV-2 does not eclipse 15 days increased the possible transmission pairs to n = 158 (Figure S5), no significant relationship between ACE2 transcription and the number of secondary cases was estimated using a Poisson ( β = 0.0143 [95%CI: -0.284 to 0.313], P = 0.64) or NB regression model ( β = 3.90 × 10− 3 [95%CI: -0.0849 to 0.0911], P = 0.8) (Fig. 3). Additionally, assumption four (iv) was challenged by decreasing the infectious period to 9.5 days or ~ 2 serial intervals (Figure S6, transmission pairs n = 48). Applying a shorter infectious period provided the following estimates using a Poisson ( β = -0.270 [95%CI: -0.454 to -0.0852], P = 0.00) and NB regression model (B = -0.160 [95%CI: -0.324 to -0.0118], P = 0.09) (Fig. 3). Importantly, the effect modification term was significant in the NB model, and it provided a better fit than when used in the primary analysis (AIC = 269). In patients a one-unit change in ACE2 transcription increases the number of secondary cases by 1.307 (95% CI: 1.03 to 1.67, P = 0.04), when adjusting for RNA viral load (Fig. 3).

Fig. 3
figure 3

Forest plot of point estimates from the association between nasopharyngeal ACE2 transcription of the primary case and the SARS-CoV-2 secondary cases. Poisson and Negative Binomial (NB) regression models were used for the primary analysis and sensitivity analysis of the assumptions used to construct the infection tracing networks. (A) βeta-coefficients (β) and 95% confidence intervals for the association between nasopharyngeal ACE2 transcription and secondary cases within health care workers, adjusting for viral load. The best fit NB model shows no association in the primary analysis or sensitivity analysis (ii, iv or v). A significant association between nasopharyngeal ACE2 transcription and the number of SARS-CoV-2 secondary cases was found for healthcare workers in sensitivity analysis iii, where the serial interval of 5.6 days was not used to build the transmission network. (B) βeta-coefficientsand 95% confidence intervals for the association between nasopharyngeal ACE2 transcription and secondary cases within patients, adjusting for viral load. The best fit NB model finds a positive association between in the primary analysis, and sensitivity analysis iii, and v. No association between nasopharyngeal ACE2 transcription and the number of SARS-CoV-2 secondary cases was found for patients in sensitivity analysis ii and iv, where assumptions of transmission occurring within hospital units and a 15.6-day infectious period were excluded from constructing the transmission network

Discussion

This study describes case-to-case transmission within a series of hospital associated SARS-CoV-2 outbreaks in British Columbia during late 2020 and early 2021. We used a combination of epidemiological data from outbreak reports, laboratory measurements of ACE2 transcription, and SARS-CoV-2 lineages derived from whole genome sequencing to recreate the outbreaks using infection tracing transmission network analysis [48]. The transmission network was resolved to provide the possible number of transmission events (secondary cases) per primary case. Bivariate and multivariable analysis was employed to estimate the relationship between nasopharyngeal ACE2 transcription and the number of potential secondary cases.

Bivariate analysis found that nasopharyngeal ACE2 transcription did not differ between patients and health care workers. However, ACE2 transcription was enriched in cases, more than 50% of the cases had ‘High’ ACE2. The transmission network constructed with all four assumptions provided n = 76 transmission events for analysis. The simplest recreation of the outbreaks was provided by a sensitivity analysis where the fourth assumption (iv) was changed to stipulate an infectious period of 9.5 days (~ 2 serial intervals), which yielded n = 48 transmission events. Multivariable analysis with Poisson regression models was inferior to NB regression as the number of secondary cases was over dispersed. In the primary analysis, NB regression estimated that a one-unit change in nasopharyngeal ACE2 transcription decreases the number of secondary cases in health care workers by β = -0.132 (95%CI: -0.255 to -0.0181, P = 0.04) and increases transmission in patients by β = 0.187 (95% CI: 0.0101 to 0.370, P = 0.04), adjusting for SARS-CoV-2 viral load. An association between primary case ACE2 expression and secondary transmission after adjusting for viral load suggests that RNA viral load may not serve as a proxy for infectiousness or that “High” ACE2 expression extends the duration of shedding overtime [49]. Sensitivity analysis using a shorter infectious period of 9.5 days, provided comparable results, no association was found between nasopharyngeal ACE2 transcription and SARS-CoV-2 transmission from health care workers (P = 0.09). In patients, nasopharyngeal ACE2 expression was associated with secondary transmission β = 0.27 (95% CI: 0.0315 to 0.518, P = 0.04) adjusting for RNA viral load (Fig. 3A & B).

Viral transmission exemplifies a complex system where understanding the mechanism of transmission requires deeper analysis than observing the sum of its components. This observational study describes nasopharyngeal ACE2 transcription as a component of SARS-CoV-2 transmission. In a series of hospital associated outbreaks in British Columbia nasopharyngeal ACE2 transcription was positively associated with the number of secondary cases (ongoing transmission) in patients but not health care workers. The reason for this difference may stem from health care workers wearing personal protective equipment while providing care. The Public Health Agency of Canada recommends that health care workers providing care to COVID-19 positive patients, as a minimum, adhere to droplet and contact precautions [50]. If care requires aerosol generating procedures, then the health care worker dons a N95 respirator with a gown, gloves, and eye protection [50]. Surgical masks and N95 respirators provide a physical barrier between the mucosal membrane of the nasopharyngeal passage and virus carrying droplets or particles suspended in the air. The design and intended use of respirators and surgical masks differ. Respirators are designed to prevent inhalation of airborne particles and must fit tightly to the user’s face. Surgical masks are designed to protect others from aerosolized droplet production from the wearer’s upper respiratory system [51]. In either case, the physical barrier provided by a mask can prevent SARS-CoV-2 viral particles from binding to ACE2 proteins on the surface of host cells and initiating the viral replication cycle. Masking may not completely prevent exposure to SARS-CoV-2 yet still decreases the infectious dose below the level necessary for infection [52]. Various studies have demonstrated the utility of surgical masks and N95 respirators at preventing COVID-19 or infection with other respiratory pathogens [52,53,54]. Use of surgical masks and social distancing was associated with a reduction of 44.9 COVID-19 cases per 1000 students and staff in Boston area school districts. This reduction represents a 4.49% decrease in the incidence of SARS-CoV-2 [55]. A clinical trial investigating the use of surgical masks or N95 respirators to prevent SARS-CoV-2 infection in health care workers found no difference in the masks’ efficacy to protect from infection [56]. This suggests that when worn by health care workers, surgical masks and N95 respirators have comparable efficacy to reduce respiratory infections transmitted by small, aerosolized particles < 5 μm. Taken together, this evidence supports the hypothesis that health care workers in our study were sufficiently protected by their personal protective equipment to negate the association between nasopharyngeal ACE2 transcription and secondary transmission. In contrast, hospitalized patients did not wear PPE while receiving care and nasopharyngeal ACE2 transcription possessed a positive association with secondary transmission. The principle of a mask protecting the wearer from causing transmission instead of protecting contacts from acquiring infection parallels how high nasopharyngeal ACE2 transcription in an infected person may produce more secondary cases. High nasopharyngeal ACE2 transcription may promote viral shedding in the upper respiratory tract and thus the production of finer aerosols from talking, singing, or sneezing. Fine aerosols can spread further than the large ones, produced lower in the respiratory tract, and pose an increased risk of infection to contacts [57]. As PPE use was not measured in our study competing hypotheses could explain the apparent difference in the relationship between ACE2 transcription and transmission among health care workers and patients. An unmeasured confounding variable present for health care workers but not patients could bias the relationship between ACE2 transcription and secondary transmission. For example, health care workers have variable shift lengths and leave the hospital daily after work; therefore, their contact time with an infected patient differs from that of an inpatient who shares the same hospital unit. One assumption used to construct the transmission network was that primary cases’ unit of transmission equaled the secondary cases’ unit of acquisition which is more relevant for patients than healthcare workers as some health care workers may routinely move between units. This example also serves to demonstrate the complexity of measuring SARS-CoV-2 and other respiratory pathogen’s transmission. We acknowledge that SARS-CoV-2 transmission occurs due to an intricate balance of social, biological, and environmental risk-factors of which nasopharyngeal ACE2 transcription may contribute a singular role [7].

Regulators of ACE2 transcription in the lower respiratory tract are well documented a lesser evidence base exists for the upper respiratory tract. Importantly, upper, and lower respiratory tract ACE2 transcription are not strongly corelated indicating that they may occur independently of each other [58]. In previous work we found that nasopharyngeal transcription of ACE2 did not differ by age (within adults over the age of eighteen) or biological sex [28], a finding which was further supported by this study. High nasopharyngeal ACE2 expression has been associated with long term inhaled corticosteroid use and exposure to fine particulate matter PM2.5. Below average, low nasopharyngeal ACE2 transcription is associated with type 2 inflammation mediated by the cytokines IL-4, IL-5, and IL-13 [59, 60]. Therefore, individuals with asthma or allergic rhinitis may have lower upper respiratory tract ACE2 expression than those that do not. Further studies are required to understand additional predictors of upper respiratory tract ACE2 expression and if it relates to that of the lower respiratory tract.

Strengths and limitations

The described study has several limitations in design, data collection and analysis. Sampling cases of SARS-CoV-2 but not their contacts prevented us from constructing a more robust transmission network [48]. Contact tracing data would benefit our analysis by increasing power of the study and allow estimation of the relationship between nasopharyngeal ACE2 transcription and lack of transmission. Estimating the relationship between nasopharyngeal ACE2 transcription and secondary transmission in only confirmed cases may have resulted in selection bias, where we selected for ACE2 measurements in persons that had been exposed and infected by SARS-CoV-2 [61]. Sampling from a population of hospitalized patients may also have introduced Berkson’s bias (admission rate bias) into the study and restricts our ability to generalize the results to other sub-groups or the community [61]. The hospitalized patients could share an unmeasured exposure, comorbidity or drug treatment which unknowingly affected their nasopharyngeal ACE2 transcription. Though we would expect a greater propensity for comorbidities in patients as they are hospitalized for a reason, healthcare workers may also have comorbidities that affect their upper respiratory tract ACE2 expression. Performing a sensitivity analysis helped us understand the importance of our assumptions in building the infection tracing transmission network. The assumption that the primary and secondary case shared the same hospital unit, proved the most important. Without knowing the place of exposure, the transmission network would not have resolved transmission pairs. Whole genome sequencing provided limited specificity of transmission events when the SARS-CoV-2 outbreak was short and genomes were classified by lineage and not at the resolution of single nucleotide polymorphisms (SNPs). Incorporating SNPs into the analysis may increase resolution of the transmission network by more accurately determining transmission pairs. Using RT-qPCR to measure ACE2 transcription and RNA viral load overapproximates available ACE2 protein and viral particles [24]. To categorize ACE2 transcription we used a reference group of SARS-CoV-2 negative specimens that were collected using a different collection medium (Hologic® Aptima®) [28], than that used to collect specimens from positive testers in the hospital outbreaks (COPAN UTM®). Universal transport medium (UTM) is a sub-optimal collection medium for molecular testing as cellular process and viral replication may still take place in the clinical specimen until cold storage at -20 °C [62]. In contrast, molecular grade medium like Hologic® Aptima® contains detergents, chelators and chaotropic agents which arrest cellular processes, denature enzymes, and inactivate infectious viral particles, making gene expression measurements more accurate [63]. During the COVID-19 pandemic many common testing supplies were limited or exhausted, comparison of different swabs and transport mediums found marginal differences for measuring viral gene expression [64]. Future work should aim to re-estimate the association between nasopharyngeal ACE2 transcription and secondary cases using data from a SARS-CoV-2 household transmission study. The design of a household transmission study has the benefits of selecting participants from a non-hospitalized population, across a wide range of ages while observing secondary infections and contacts within family clusters in a context with limited use of personal protective equipment.

Conclusion

We estimate the association between nasopharyngeal ACE2 transcription and secondary SARS-CoV-2 transmission in a series of hospital associated outbreaks in British Columbia from late 2020 to early 2021. Analysis shows that 98.6% of transmission pairs and 85% of primary cases had high nasopharyngeal ACE2 transcription. Multivariable analysis adjusting for RNA viral load and interaction by case description found that nasopharyngeal ACE2 transcription was positively associated with SARS-CoV-2 transmission in hospital patients but not health care workers. We postulate that use of masks among health care workers explains this difference, transmission from health care workers with high ACE2 was interrupted by barrier protection. SARS-CoV-2 transmission remains a complex and unresolved phenomenon driven by biological, environmental, and social risk factors. Differential nasopharyngeal ACE2 transcription has been observed in COVID-19 negative persons and having high- above average- expression of ACE2 may serve as a risk factor for transmission of SARS-CoV-2.

Data availability

The datasets generated and/or analysed during the current study are available in the “figshare” repository, https://doi.org/10.6084/m9.figshare.24158637.v1.

References

  1. Long B, Brady WJ, Koyfman A, Gottlieb M. Cardiovascular complications in COVID-19. Am J Emerg Med. 2020;38:1504–7. https://doi.org/10.1016/j.ajem.2020.04.048

    Article  PubMed  PubMed Central  Google Scholar 

  2. Dhar D, Dey T, Samim MM, Padmanabha H, Chatterjee A, Naznin P, Chandra SR, Mallesh K, Shah R, Siddiqui S, et al. Systemic inflammatory syndrome in COVID-19–SISCoV study: systematic review and meta-analysis. Pediatr Res. 2022;91:1334–49. https://doi.org/10.1038/s41390-021-01545-z

    Article  CAS  PubMed  Google Scholar 

  3. Guo T, Fan Y, Chen M, Wu X, Zhang L, He T, Wang H, Wan J, Wang X, Lu Z. Cardiovascular implications of fatal outcomes of patients with coronavirus disease 2019 (COVID-19). JAMA Cardiol. 2020;5:811. https://doi.org/10.1001/jamacardio.2020.1017

    Article  PubMed  Google Scholar 

  4. Chan-Yeung M, Xu R-H, SARS. Epidemiol Respirol. 2003;8:9–S14. https://doi.org/10.1046/j.1440-1843.2003.00518.x

    Article  Google Scholar 

  5. Zhang A-R, Shi W-Q, Liu K, Li X-L, Liu M-J, Zhang W-H, Zhao G-P, Chen J-J, Zhang X-A, Miao D, et al. Epidemiology and evolution of Middle East respiratory syndrome coronavirus, 2012–2020. Infect Dis Poverty. 2021;10. https://doi.org/10.1186/s40249-021-00853-0

  6. Singanayagam A, Hakki S, Dunning J, Madon KJ, Crone MA, Koycheva A, Derqui-Fernandez N, Barnett JL, Whitfield MG, Varro R, et al. Community transmission and viral load kinetics of the SARS-CoV-2 Delta (B.1.617.2) variant in vaccinated and unvaccinated individuals in the UK: a prospective, longitudinal, cohort study. Lancet Infect Dis. 2022;22:183–95. https://doi.org/10.1016/S1473-3099(21)00648-4

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Leung NHL. Transmissibility and transmission of respiratory viruses. Nat Rev Microbiol. 2021;19:528–45. https://doi.org/10.1038/s41579-021-00535-6

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Lindsey BB, Villabona-Arenas CJ, Campbell F, Keeley AJ, Parker MD, Shah DR, Parsons H, Zhang P, Kakkar N, Gallis M, et al. Characterising within-hospital SARS-CoV-2 transmission events using epidemiological and viral genomic data across two pandemic waves. Nat Commun. 2022;13. https://doi.org/10.1038/s41467-022-28291-y

  9. Madewell ZJ, Yang Y, Longini IM, Halloran ME, Dean NE. Household transmission of SARS-CoV-2. JAMA Netw Open. 2020;3:e2031756. https://doi.org/10.1001/jamanetworkopen.2020.31756

    Article  PubMed  PubMed Central  Google Scholar 

  10. Cowger TL, Murray EJ, Clarke J, Bassett MT, Ojikutu BO, Sánchez SM, Linos N, Hall KT. Lifting universal masking in schools — Covid-19 incidence among students and staff. N Engl J Med. 2022. https://doi.org/10.1056/NEJMoa2211029

    Article  PubMed  PubMed Central  Google Scholar 

  11. Cordery R, Reeves L, Zhou J, Rowan A, Watber P, Rosadas C, Crone M, Storch M, Freemont P, Mosscrop L, et al. Transmission of SARS-CoV-2 by children to contacts in schools and households: a prospective cohort and environmental sampling study in London. Lancet Microbe. 2022;3:e814–23. https://doi.org/10.1016/S2666-5247(22)00124-0

    Article  PubMed  PubMed Central  Google Scholar 

  12. Li F, Li W, Farzan M, Harrison SC. Structure of SARS coronavirus spike receptor-binding domain complexed with receptor. Sci (1979). 2005;309:1864–8. https://doi.org/10.1126/science.1116480

    Article  CAS  Google Scholar 

  13. Mannar D, Saville JW, Zhu X, Srivastava SS, Berezuk AM, Tuttle KS, Marquez AC, Sekirov I, Subramaniam S. SARS-CoV-2 Omicron variant: antibody evasion and Cryo-EM structure of spike protein–ACE2 complex. Science (1979). 2022;375:760–764, https://doi.org/10.1126/science.abn7760

  14. Shang J, Ye G, Shi K, Wan Y, Luo C, Aihara H, Geng Q, Auerbach A, Li F. Structural basis of receptor recognition by SARS-CoV-2. Nature. 2020;581:221–4. https://doi.org/10.1038/s41586-020-2179-y

    Article  CAS  PubMed  PubMed Central  ADS  Google Scholar 

  15. Meng B, Abdullahi A, Ferreira IATM, Goonawardane N, Saito A, Kimura I, Yamasoba D, Gerber PP, Fatihi S, Rathore S, et al. Altered TMPRSS2 usage by SARS-CoV-2 Omicron impacts infectivity and fusogenicity. Nature. 2022;603:706–14. https://doi.org/10.1038/s41586-022-04474-x

    Article  CAS  PubMed  PubMed Central  ADS  Google Scholar 

  16. Cantuti-Castelvetri L, Ojha R, Pedro LD, Djannatian M, Franz J, Kuivanen S, van der Meer F, Kallio K, Kaya T, Anastasina M, et al. Neuropilin-1 facilitates SARS-CoV-2 cell entry and infectivity. Sci (1979). 2020;370:856–60. https://doi.org/10.1126/science.abd2985

    Article  CAS  Google Scholar 

  17. Daly JL, Simonetti B, Klein K, Chen K-E, Williamson MK, Antón-Plágaro C, Shoemark DK, Simón-Gracia L, Bauer M, Hollandi R, et al. Neuropilin-1 is a host factor for SARS-CoV-2 infection. Sci (1979). 2020;370:861–5. https://doi.org/10.1126/science.abd3072

    Article  CAS  Google Scholar 

  18. Ziegler CGK, Allon SJ, Nyquist SK, Mbano IM, Miao VN, Tzouanas CN, Cao Y, Yousif AS, Bals J, Hauser BM, et al. SARS-CoV-2 receptor ACE2 is an interferon-stimulated gene in human airway epithelial cells and is detected in specific cell subsets across tissues. Cell. 2020;181:1016–1035e19. https://doi.org/10.1016/j.cell.2020.04.035

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Sui Y, Li J, Venzon DJ, Berzofsky JA. SARS-CoV-2 spike protein suppresses ACE2 and type I interferon expression in primary cells from macaque lung bronchoalveolar lavage. Front Immunol. 2021;12. https://doi.org/10.3389/fimmu.2021.658428

  20. Verdecchia P, Cavallini C, Spanevello A, Angeli F. The pivotal link between ACE2 deficiency and SARS-CoV-2 infection. Eur J Intern Med. 2020;76:14–20. https://doi.org/10.1016/j.ejim.2020.04.037

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Groß S, Jahn C, Cushman S, Bär C, Thum T. SARS-CoV-2 receptor ACE2-dependent implications on the cardiovascular system: from basic science to clinical implications. J Mol Cell Cardiol. 2020;144:47–53. https://doi.org/10.1016/j.yjmcc.2020.04.031

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Jia HP, Look DC, Shi L, Hickey M, Pewe L, Netland J, Farzan M, Wohlford-Lenane C, Perlman S, McCray PB. ACE2 receptor expression and severe acute respiratory syndrome coronavirus infection depend on differentiation of human airway epithelia. J Virol. 2005;79:14614–21. https://doi.org/10.1128/JVI.79.23.14614-14621.2005

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Ni W, Yang X, Yang D, Bao J, Li R, Xiao Y, Hou C, Wang H, Liu J, Yang D, et al. Role of angiotensin-converting enzyme 2 (ACE2) in COVID-19. Crit Care. 2020;24:422. https://doi.org/10.1186/s13054-020-03120-0

    Article  PubMed  PubMed Central  Google Scholar 

  24. Nikiforuk AM, Kuchinski KS, Twa DDW, Lukac CD, Sbihi H, Basham CA, Steidl C, Prystajecky NA, Jassem AN, Krajden M, et al. The contrasting role of nasopharyngeal angiotensin converting enzyme 2 (ACE2) transcription in SARS-CoV-2 infection: a cross-sectional study of people tested for COVID-19 in British Columbia, Canada. EBioMedicine. 2021;66:103316. https://doi.org/10.1016/j.ebiom.2021.103316

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Freed NE, Vlková M, Faisal MB, Silander OK. Rapid and inexpensive whole-genome sequencing of SARS-CoV-2 using 1200 Bp tiled amplicons and Oxford nanopore rapid barcoding. Biol Methods Protoc. 2020;5. https://doi.org/10.1093/biomethods/bpaa014

  26. Messmann R, Dietl A, Wagner S, Domenig O, Jungbauer C, Luchner A, Maier LS, Schopka S, Hirt S, Schmid C, et al. Alterations of the renin angiotensin system in human end-stage heart failure before and after mechanical cardiac unloading by LVAD support. Mol Cell Biochem. 2020;472:79–94. https://doi.org/10.1007/s11010-020-03787-7

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. LeBlanc JJ, Gubbay JB, Li Y, Needle R, Arneson SR, Marcino D, Charest H, Desnoyers G, Dust K, Fattouh R, et al. Real-time PCR-based SARS-CoV-2 detection in Canadian laboratories. J Clin Virol. 2020;128:104433. https://doi.org/10.1016/j.jcv.2020.104433

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Nikiforuk AM, Kuchinski KS, Twa DDW, Lukac CD, Sbihi H, Basham CA, Steidl C, Prystajecky NA, Jassem AN, Krajden M et al. The contrasting role of nasopharyngeal angiotensin converting enzyme 2 (ACE2) transcription in SARS-CoV-2 infection: a cross-sectional study of people tested for COVID-19 in British Columbia, Canada. EBioMedicine. 2021;66:103316, https://doi.org/10.1016/j.ebiom.2021.103316

  29. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2– ∆∆CT method. Methods. 2001;25:402–8. https://doi.org/10.1006/meth.2001.1262

    Article  CAS  PubMed  Google Scholar 

  30. British columbia centre for disease control Ncov2019-Artic-Nf (accessed on 22 Nov 2022). https://github.com/BCCDC-PHL/ncov-tools-nf

  31. O’Toole Á, Scher E, Underwood A, Jackson B, Hill V, McCrone JT, Colquhoun R, Ruis C, Abu-Dahab K, Taylor B, et al. Assignment of epidemiological lineages in an emerging pandemic using the Pangolin tool. Virus Evol. 2021. https://doi.org/10.1093/ve/veab064

    Article  PubMed  PubMed Central  Google Scholar 

  32. VanderWeele TJ. Principles of Confounder Selection. Eur J Epidemiol. 2019;34:211–9. https://doi.org/10.1007/s10654-019-00494-6

    Article  PubMed  PubMed Central  Google Scholar 

  33. Gutiérrez-Chamorro L, Riveira-Muñoz E, Barrios C, Palau V, Nevot M, Pedreño-López S, Senserrich J, Massanella M, Clotet B, Cabrera C, et al. SARS-CoV-2 infection modulates ACE2 function and subsequent inflammatory responses in swabs and plasma of COVID-19 patients. Viruses. 2021;13:1715. https://doi.org/10.3390/v13091715

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Marc A, Kerioui M, Blanquart F, Bertrand J, Mitjà O, Corbacho-Monné M, Marks M, Guedj J. Quantifying the relationship between SARS-CoV-2 viral load and infectiousness. Elife. 2021;10. https://doi.org/10.7554/eLife.69302

  35. Carraturo F, del Giudice C, Morelli M, Cerullo V, Libralato G, Galdiero E, Guida M. Persistence of SARS-CoV-2 in the environment and COVID-19 transmission risk from environmental matrices and surfaces. Environ Pollut. 2020;265:115010. https://doi.org/10.1016/j.envpol.2020.115010

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Onweni CL, Zhang YS, Caulfield T, Hopkins CE, Fairweather DL, Freeman WD. ACEI/ARB therapy in COVID-19: the double-edged sword of ACE2 and SARS-CoV-2 viral docking. Crit Care. 2020;24:475. https://doi.org/10.1186/s13054-020-03195-9

    Article  PubMed  PubMed Central  Google Scholar 

  37. Ortiz ME, Thurman A, Pezzulo AA, Leidinger MR, Klesney-Tait JA, Karp PH, Tan P, Wohlford-Lenane C, McCray PB, Meyerholz DK. Heterogeneous expression of the SARS-Coronavirus-2 receptor ACE2 in the human respiratory tract. EBioMedicine. 2020;60:102976, https://doi.org/10.1016/j.ebiom.2020.102976

  38. World health organisation. Tracking SARS-CoV-2 variants. https://www.who.int/activities/tracking-SARS-CoV-2-variants

  39. Kwak SG, Kim JH. Central limit theorem: the cornerstone of modern statistics. Korean J Anesthesiol. 2017;70. https://doi.org/10.4097/kjae.2017.70.2.144

  40. Nagarajan K, Muniyandi M, Palani B, Sellappan S. Social network analysis methods for exploring SARS-CoV-2 contact tracing data. BMC Med Res Methodol. 2020;20:233. https://doi.org/10.1186/s12874-020-01119-3

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Alene M, Yismaw L, Assemie MA, Ketema DB, Gietaneh W, Birhan TY. Serial interval and incubation period of COVID-19: a systematic review and meta-analysis. BMC Infect Dis. 2021;21:257. https://doi.org/10.1186/s12879-021-05950-x

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Byrne AW, McEvoy D, Collins AB, Hunt K, Casey M, Barber A, Butler F, Griffin J, Lane EA, McAloon C, et al. Inferred duration of infectious period of SARS-CoV-2: rapid scoping review and analysis of available evidence for asymptomatic and symptomatic COVID-19 cases. BMJ Open. 2020;10:e039856. https://doi.org/10.1136/bmjopen-2020-039856

    Article  PubMed  PubMed Central  Google Scholar 

  43. Craney TA, Surles JG. Model-dependent variance inflation factor cutoff values. Qual Eng. 2002;14:391–403. https://doi.org/10.1081/QEN-120001878

    Article  Google Scholar 

  44. Cameron AC, Trivedi PK. Regression-based tests for overdispersion in the Poisson model. J Econom. 1990;46:347–64. https://doi.org/10.1016/0304-4076(90)90014-K

    Article  MathSciNet  Google Scholar 

  45. van den Broek J. A score test for zero inflation in a poisson distribution. Biom. 1995;51:738–43. https://doi.org/10.2307/2532959

    Article  MathSciNet  PubMed  Google Scholar 

  46. Portet S. A primer on model selection using the Akaike information criterion. Infect Dis Model. 2020;5:111–28. https://doi.org/10.1016/j.idm.2019.12.010

    Article  PubMed  PubMed Central  Google Scholar 

  47. Posit team. RStudio: integrated development environment for R. Posit software, PBC, boston, MA. 2024. http://www.posit.co/.

  48. Keeling MJ, Eames KTD. Networks and epidemic models. J R Soc Interface. 2005;2:295–307. https://doi.org/10.1098/rsif.2005.0051

    Article  PubMed  PubMed Central  Google Scholar 

  49. Hakki S, Zhou J, Jonnerby J, Singanayagam A, Barnett JL, Madon KJ, Koycheva A, Kelly C, Houston H, Nevin S, et al. Onset and window of SARS-CoV-2 infectiousness and temporal correlation with symptom onset: a prospective, longitudinal, community cohort study. Lancet Respir Med. 2022;10:1061–73. https://doi.org/10.1016/S2213-2600(22)00226-0

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Public health agency of canada infection prevention and control for COVID-19: interim guidance for acute healthcare settings. https://www.canada.ca/en/public-health/services/diseases/2019-novel-coronavirus-infection/health-professionals/infection-prevention-control.html

  51. Vincent M, Edwards P. Disposable surgical face masks for preventing surgical wound infection in clean surgery. Cochrane Database of Systematic Reviews. 2016;2016. https://doi.org/10.1002/14651858.CD002929.pub3

  52. Leung NHL, Chu DKW, Shiu EYC, Chan K-H, McDevitt JJ, Hau BJP, Yen H-L, Li Y, Ip DKM, Peiris JSM, et al. Respiratory virus shedding in exhaled breath and efficacy of face masks. Nat Med. 2020;26:676–80. https://doi.org/10.1038/s41591-020-0843-2

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Andrejko KL, Pry J, Myers JF, Openshaw J, Watt J, Birkett N, DeGuzman JL, Barbaduomo CM, Dong ZN, Fang AT, et al. Predictors of severe acute respiratory syndrome coronavirus 2 infection following high-risk exposure. Clin Infect Dis. 2022;75:e276–88. https://doi.org/10.1093/cid/ciab1040

    Article  PubMed  Google Scholar 

  54. Fennelly KP. Particle sizes of infectious aerosols: implications for infection control. Lancet Respir Med. 2020;8:914–24. https://doi.org/10.1016/S2213-2600(20)30323-4

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. Cowger TL, Murray EJ, Clarke J, Bassett MT, Ojikutu BO, Sánchez SM, Linos N, Hall KT. Lifting universal masking in schools — Covid-19 incidence among students and staff. N Engl J Med. 2022;387:1935–46. https://doi.org/10.1056/NEJMoa2211029

    Article  PubMed  Google Scholar 

  56. Loeb M, Bartholomew A, Hashmi M, Tarhuni W, Hassany M, Youngster I, Somayaji R, Larios O, Kim J, Missaghi B, et al. Medical masks versus N95 respirators for preventing COVID-19 among health care workers. Ann Intern Med. 2022;175:1629–38. https://doi.org/10.7326/M22-1966

    Article  PubMed  Google Scholar 

  57. Wang CC, Prather KA, Sznitman J, Jimenez JL, Lakdawala SS, Tufekci Z, Marr LC. Airborne transmission of respiratory viruses. Sci (1979). 2021;373. https://doi.org/10.1126/science.abd9149

  58. Hou YJ, Okuda K, Edwards CE, Martinez DR, Asakura T, Dinnon KH, Kato T, Lee RE, Yount BL, Mascenik TM, et al. SARS-CoV-2 reverse genetics reveals a variable infection gradient in the respiratory tract. Cell. 2020;182:429–446e14. https://doi.org/10.1016/j.cell.2020.05.042

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  59. Kimura H, Francisco D, Conway M, Martinez FD, Vercelli D, Polverino F, Billheimer D, Kraft M. Type 2 inflammation modulates ACE2 and TMPRSS2 in airway epithelial cells. J Allergy Clin Immunol. 2020;146:80–88e8. https://doi.org/10.1016/j.jaci.2020.05.004

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  60. Takabayashi T, Yoshida K, Imoto Y, Schleimer RP, Fujieda S. Regulation of the expression of SARS-CoV-2 receptor angiotensin-converting enzyme 2 in nasal mucosa. Am J Rhinol Allergy. 2022;36:115–22. https://doi.org/10.1177/19458924211027798

    Article  PubMed  Google Scholar 

  61. Westreich D. Berkson’s bias, selection bias, and missing data. Epidemiology. 2012;23:159–64. https://doi.org/10.1097/EDE.0b013e31823b6296

    Article  PubMed  PubMed Central  Google Scholar 

  62. COPAN universal transport Medium. https://www.copanusa.com/sample-collection-transport-processing/utm-viral-transport/

  63. Aptima HOLOGIC. SARS-CoV-2 assay. https://www.hologic.com/sites/default/files/2018-06/AW-14334-001_003_01.pdf

  64. Garnett L, Bello A, Tran KN, Audet J, Leung A, Schiffman Z, Griffin BD, Tailor N, Kobasa D, Strong JE. Comparison analysis of different swabs and transport mediums suitable for SARS-CoV-2 testing following shortages. J Virol Methods. 2020;285:113947. https://doi.org/10.1016/j.jviromet.2020.113947

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Funding

This work was funded by the Canadian Institutes of Health Research (#434951) and Genome British Columbia (COV-55).

Author information

Authors and Affiliations

Authors

Contributions

A.M.N., K.S.K., K.S. and I.S. conceived and designed the study. Epidemiological reporting, data recording and stewardship was performed by K.S. and S.R. Laboratory data reporting and stewardship was performed by N.P. and A.N.J. Laboratory measurements were taken by A.M.N and K.S.K. Data analysis was performed by A.M.N and K.S.K., D.M.P and M.A.I assisted in interpretation of the results. The manuscript was written by A.M.N, all authors helped with editing and provided their approval for publication.

Corresponding author

Correspondence to Inna Sekirov.

Ethics declarations

Ethics approval and consent to participate

Ethical approval for the study was sought from the University of British Columbia Human Ethics Board (#H20-01110) which was harmonized with the Fraser Health Authority. Written informed consent was not required and waived by the University of British Columbia Human Ethics Board. All data was de-identified prior to analysis and the results were not linked back to any identifying records. This study was deemed as minimal risk to the participants involved. To ensure privacy the site of the outbreak series will not be disclosed.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary Material 1

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Nikiforuk, A.M., Kuchinski, K.S., Short, K. et al. Nasopharyngeal angiotensin converting enzyme 2 (ACE2) expression as a risk-factor for SARS-CoV-2 transmission in concurrent hospital associated outbreaks. BMC Infect Dis 24, 262 (2024). https://doi.org/10.1186/s12879-024-09067-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12879-024-09067-9

Keywords