# SARS-CoV-2 variants and mutational patterns: relationship with risk of ventilator-associated pneumonia in critically ill COVID-19 patients in the era of dexamethasone

Keyvan Razazi, Anissa Martins Bexiga, Romain Arrestier, Bastien Peiffer, Guillaume Voiriot, Charles Édouard Luyt, Tomas Urbina, Julien Mayaux, Tài Pham, Damien Roux, et al.

► To cite this version:

Keyvan Razazi, Anissa Martins Bexiga, Romain Arrestier, Bastien Peiffer, Guillaume Voiriot, et al. SARS-CoV-2 variants and mutational patterns: relationship with risk of ventilator-associated pneumonia in critically ill COVID-19 patients in the era of dexamethasone. Scientific Reports, 2023, 13 (1), 10.1038/s41598-023-33639-5. [https://doi.org/10.1038/s41598-023-33639-5](https://doi.org/10.1038/s41598-023-33639-5)

HAL Id: hal-04106275

https://hal.science/hal-04106275

Submitted on 2 Jun 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

OPEN

# SARS-CoV-2 variants and mutational patterns: relationship with risk of ventilator-associated pneumonia in critically ill COVID-19 patients in the era of dexamethasone

Keyvan Razazi<fnref n="1" />,²,³,¹⁹✉, Anissa Martins Bexiga<fnref n="1" />,²,¹⁹, Romain Arrestier<fnref n="1" />,², Bastien Peiffer⁴, Guillaume Voiriot⁵, Charles-Edouard Luyt⁶, Tomas Urbina⁷, Julien Mayaux⁸, Tài Pham³,⁹,¹⁰, Damien Roux¹¹,¹², Raphael Bellaiche¹³, Zakaria Alt Hamou¹⁴, Stéphane Gaudry¹⁵, Elie Azoulay¹⁶, Armand Mekontso Dessap<fnref n="1" />,²,³, Christophe Rodriguez³,¹⁷,¹⁸, Jean-Michel Pawlotsky³,¹⁷,¹⁸, Slim Fourati³,¹⁷,¹⁸,²⁰ & Nicolas de Prost<fnref n="1" />,²,³,²⁰

We aimed to explore the relationships between specific viral mutations/mutational patterns and ventilator-associated pneumonia (VAP) occurrence in COVID-19 patients admitted in intensive care units between October 1, 2020, and May 30, 2021. Full-length SARS-CoV-2 genomes were sequenced by means of next-generation sequencing. In this prospective multicentre cohort study, 259 patients were included. 222 patients (47%) had been infected with pre-existing ancestral variants, 116 (45%) with variant α, and 21 (8%) with other variants. 153 patients (59%) developed at least one VAP. There was no significant relationship between VAP occurrence and a specific SARS CoV-2 lineage/sublineage or mutational pattern.

<footnote n="1">¹Médecine Intensive Réanimation, Hôpitaux Universitaires Henri Mondor, Assistance Publique-Hôpitaux de Paris (AP-HP), 51, Av de Lattre de Tassigny, 94000 Créteil Cedex, France. ²Groupe de Recherche Clinique CARMAS, Université Paris-Est-Créteil (UPEC), Créteil, France. ³Université Paris-Est-Créteil (UPEC), Créteil, France. ⁴DMU Medecine, Hôpitaux Universitaires Henri Mondor, Assistance Publique-Hôpitaux de Paris (AP-HP), Créteil, France. ⁵Médecine Intensive Réanimation, Hôpital Tenon, Assistance Publique-Hôpitaux de Paris, Paris, France. ⁶Sorbonne Université, Assistance Publique-Hôpitaux de Paris, Hôpital Pitié-Salpêtrière, Médecine Intensive Réanimation, and INSERM UMRS 1166-iCAN, Institute of Cardiometabolism and Nutrition, Paris, France. ⁷Sorbonne Université, Assistance Publique-Hôpitaux de Paris, Hôpital Saint-Antoine, Médecine Intensive Réanimation, 75571 Paris Cedex 12, France. ⁸Sorbonne Université, Assistance Publique-Hôpitaux de Paris, Hôpital Pitié-Salpêtrière, Médecine Intensive Réanimation, Paris, France. ⁹Service de Médecine Intensive-Réanimation, Assistance Publique-Hôpitaux de Paris, Hôpital de Bicêtre, DMU 4 CORREVE Maladies du Cœur et des Vaisseaux, FHU Sepsis, Le Kremlin-Bicêtre, France. ¹⁰Université Paris-Saclay, UVSQ, Univ. Paris-Sud, Inserm U1018, Equipe d'Epidémiologie Respiratoire Intégrative, CESP, 94807 Villejuif, France. ¹¹Université de Paris, APHP, Hôpital Louis Mourier, DMU ESPRIT, Service de Médecine Intensive Réanimation, Colombes, France. ¹²Department of Immunology, Infectiology and Hematology, INSERM U1151, CNRS UMR 8253, Institut Necker-Enfants Malades (INEM), Colombes, Paris, France. ¹³Département d'Anesthésie Réanimations Chirurgicales, Hôpitaux Universitaires Henri Mondor, Assistance Publique-Hôpitaux de Paris, Créteil, France. ¹⁴Médecine Intensive Réanimation, Hôpital Cochin, Assistance Publique-Hôpitaux de Paris, Paris, France. ¹⁵Service de Réanimation, Hôpital Avicenne, Assistance Publique-Hôpitaux de Paris, Bobigny, France. ¹⁶Médecine Intensive Réanimation, Hôpital Saint-Louis, Assistance Publique-Hôpitaux de Paris, Paris, France. ¹⁷Department of Virology, Hôpitaux Universitaires Henri Mondor, Assistance Publique-Hôpitaux de Paris, Créteil, France. ¹⁸INSERM U955, Team «Viruses, Hepatology, Cancer», Créteil, France. ¹⁹These authors contributed equally: Keyvan Razazi and Anissa Martins Bexiga. ²⁰These authors jointly supervised this work: Slim Fourati and Nicolas de Prost. ✉email: Keyvan.razazi@aphp.fr</footnote>
## Abbreviations

<table><tr><td>ARDS</td><td>Acute respiratory distress syndrome</td></tr><tr><td>ICU</td><td>Intensive care unit</td></tr><tr><td>ECMO</td><td>Extra-corporeal membrane oxygenation</td></tr><tr><td>VOC</td><td>Variant of concern</td></tr><tr><td>VOI</td><td>Variant of interest</td></tr><tr><td>VAP</td><td>Ventilator associated pneumonia</td></tr><tr><td>COVID 19</td><td>Coronavirus disease 2019</td></tr></table>

Current data on ventilator-associated pneumonia (VAP) complicating COVID-19-associated acute respiratory distress syndrome (ARDS) have been mainly obtained during the first wave of COVID-19¹. Yet, bacterial super-infection may be affected by the nature of the infecting viral variants and the widespread use of corticosteroids². There is scarce information in the literature on the rate of superinfections in COVID-19 patients, especially of VAP, according to variant types in the dexamethasone era. Potential explanations for the high incidence of VAP in COVID-19 patients include the long duration of invasive mechanical ventilation support, the high incidence of ARDS, and the administration of immune-suppressive treatment. Specific risk factors for VAP, including SARS-CoV-2-related pulmonary lesions and bacteria-virus interaction in lung microbiota, might also play a role in VAP pathogenesis. As the upper respiratory tract microbiome composition significantly changed as COVID-19 severity increased³, we hypothesize that specific SARS-CoV-2 mutations, either in the spike protein or in another viral protein, could have an impact on VAP prevalence.

In this study, we aimed to explore the relationships between SARS-CoV-2 variant lineage/sublineage and specific viral mutations/mutational patterns with VAP occurrence in COVID-19 patients in the dexamethasone era.

## Patients and methods

This is a prospective multicentre observational cohort study. Patients admitted between October 1, 2020 (week 40/20) and May 30, 2021 (week 21/21) in one of the 11 participating ICUs from Greater Paris area hospitals included in the ANTICOV study (NCT04733105) were eligible for study inclusion. Inclusion criteria were as follows: age ≥ 18 years, SARS-CoV-2 infection confirmed by a positive reverse transcriptase-polymerase chain reaction (RT-PCR), patient admitted in the ICU who required mechanical ventilation more than 48 h for acute respiratory failure. Patients with a low viral load (PCR cycle threshold (Ct) > 32) in nasopharyngeal samples were not included in the study.

The primary endpoint was the difference in prevalence of the first VAP according to SARS-CoV-2 variants.

**Ventilator associated pneumonia definition.** All VAP episodes prospectively diagnosed by attending physicians were evaluated and retained for analyses, provided that they fulfilled the following criteria (the definition of VAP was based on current guidelines⁴): (1) new or progressive persistent pulmonary infiltrates on chest X-ray; combined with (2) purulent tracheal secretions; (3) fever or hypothermia (body temperature ≥ 38.5 °C or ≤ 36.5 °C, respectively) and/or leukocytosis or leukopenia (white blood cells count ≥ 10.4/mL or ≤ 4 × 10.3/mL, respectively); and (4) a positive quantitative lower respiratory tract sample (endotracheal aspirate [ETA] ≥ 10.5 colony-forming unit [CFU]/mL, bronchoalveolar lavage [BAL] fluid ≥ 10.4 CFU/mL or plugged telescopic catheter [PTC] ≥ 10.3 CFU/mL). Respiratory tract secretions were cultured for VAP diagnosis and susceptibility profiles of recovered microorganisms were recorded. VAP onset was defined as the day on which the lung sample tested positive.

**SARS-CoV-2 genome sequence analysis.** Full-length SARS-CoV-2 genomes from all included patients were sequenced by means of next-generation sequencing. Viral RNA was extracted from nasopharyngeal swabs in viral transport medium using NucliSENS' easyMAG kit on EMAG device (bioMérieux, Marcy-l'Étoile, France). Sequencing was performed with the Illumina COVIDSeq Test (Illumina, San Diego, California), that uses 98-target multiplex amplifications along the full SARS-CoV-2 genome. The libraries were sequenced with NextSeq 500/550 High Output Kit v2.5 (75 Cycles) on a NextSeq 500 device (Illumina). The sequences were demultiplexed and assembled as full-length genomes by means of the DRAGEN COVIDSeq Test Pipeline on a local DRAGEN server (Illumina). Lineages and clades were interpreted using Pangolin and NextClade, before being submitted to the GISAID database (https://www.gisaid.org). Phylogeny was performed after full-length genome alignment with Muscle v3.8.31 (maximum-likelihood model GTR+I; 1000 bootstrap replicates), by means of IQ-Tree v1.3.11.1 and iTOL.

**Statistical analysis.** We used a competing risk model (cumulative incidence function of the Gray model)⁵,⁶ to properly estimate the effect of COVID-19 on VAP development, while considering death and ventilator weaning as competing events using cmprsk package developed by Gray in R software ([http://biowww.dfci.harvard.edu/~gray/cmprsk_2.1-4.tar.gz](http://biowww.dfci.harvard.edu/~gray/cmprsk_2.1-4.tar.gz)). Two-sided p values < 0.05 were considered significant.

**Ethical approval and consent to participate.** The study was approved by the Comité de Protection des Personnes Nord-Ouest IV (N° EudraCT/ID-RCB: 2020-A03009-30). Informed consent was obtained from all patients or their relatives. All study procedures and processes were conducted in accordance with the Declaration of Helsinki and with the ethical standards of the institutional and/or research committee.
## Results

Between October 1, 2020, and May 30, 2021, 845 patients were admitted in one of the 11 participating ICUs. Among them, 737 had at least one positive SARS-CoV-2 RNA RT-PCR performed in a nasopharyngeal swab sample in the hospital; 413 patients with a cycle threshold (Ct) ≤ 32 and a nasopharyngeal sample that could be used for full-length viral genome sequence analysis were included in this study. Among these 413 patients, 267 were mechanically ventilated more than 48 h. Patients' characteristics are shown in Table 1. Most patients had ARDS criteria, including 240 (93%) patients who needed prone positioning, and 52 (20%) requiring veno-venous ECMO support. Among these 259 patients, 122 (47%) had been infected with so-called ancestral "pre-existing" variants, i.e., variants circulating before the emergence of variant α, 116 (45%) with variant α (B.1.1.7), and 21 (8%) with other variants, including β (B.1.351) (n = 19, 4.6%), γ (P.1) (n = 2, 0.5%) and other variants of interest (n = 12, 2.9%). One hundred and fifty-three patients (59%) developed at least one VAP episode.

During the study period, variant α was the dominant VOC in the Greater Paris area. Thus, we first assessed whether patients infected with variant α had a higher risk of developing VAP than others (merging patients infected with preexisting variants and other variants) during the ICU stay. The Fine and Gray model showed that the probability of VAP was not significantly different in the two groups [sub-hazard ratio = 1.26 (0.62–2.54), p = 0.53, Fig. 1] after adjusting for weaning and death as competing events.

Pre-existing and emerging variants (VOC or VOI) are characterized by multiple lineage-specific deletions and amino acid substitutions in their viral genomes. The Spike glycoprotein undergoes evolutive convergence at several amino-acid signature sites found in VOC and VOI circulating at the time of the study period. These sites were selected a priori for in-depth analysis. We found no significant relationship between VAP occurrence and any specific mutation at these sites (Table 2). These results persisted after adjusting for age, sex, antibiotic at admission and SAPS II (Table 3). Although one mutation (A67V) was found to be significantly associated with

<table><caption>**Table 1.** Characteristics of patients with Coronavirus disease 19 (COVID-19) (n = 259). COPD chronic obstructive pulmonary disease, HIV human immunodeficiency virus, SAPS II Simplified Acute Physiology Score II, SOFA sequential organ failure assessment, ARDS acute respiratory distress syndrome, ICU intensive care unit; variables are shown as median [quartile 1-quartile 3] and number (percentage).</caption><thead><tr><th>Variables</th><th>Total (n = 259)</th></tr></thead><tbody><tr><td>Age—median [IQR]</td><td>65 [58–71]</td></tr><tr><td>Male gender</td><td>185 (71.4%)</td></tr><tr><td>Medical history</td><td></td></tr><tr><td>Diabetes mellitus</td><td>80 (30.9%)</td></tr><tr><td>Congestive heart failure</td><td>25 (9.7%)</td></tr><tr><td>Hypertension</td><td>138 (53.3%)</td></tr><tr><td>COPD</td><td>20 (7.7%)</td></tr><tr><td>Chronic renal failure</td><td>36 (13.9%)</td></tr><tr><td>Dialysis</td><td>13 (5.0%)</td></tr><tr><td>Liver cirrhosis</td><td>7 (2.7%)</td></tr><tr><td>Current smoking</td><td>31 (12.0%)</td></tr><tr><td>Solid cancer</td><td>15 (5.8%)</td></tr><tr><td>Blood cancer</td><td>13 (5.0%)</td></tr><tr><td>Organ transplant</td><td>15 (5.8%)</td></tr><tr><td>HIV infection</td><td>4 (1.5%)</td></tr><tr><td>Frailty scale—median [IQR]</td><td>3 [2–4]</td></tr><tr><td>Clinical characteristics upon ICU admission</td><td></td></tr><tr><td>SAPS II—median [IQR]</td><td>37 [30–48]</td></tr><tr><td>Baseline SOFA—median [IQR]</td><td>4 [3–8]</td></tr><tr><td>ARDS (Berlin definition)</td><td>230 (88.8%)</td></tr><tr><td>Norepinephrine</td><td>67 (25.9%)</td></tr><tr><td>Antibiotics</td><td>186 (71.8%)</td></tr><tr><td>Treatment during stay</td><td></td></tr><tr><td>Corticosteroids</td><td>242 (93.4%)</td></tr><tr><td>Tocilizumab</td><td>31 (12.0%)</td></tr><tr><td>Neuromuscular blockade</td><td>254 (98.1%)</td></tr><tr><td>Prone position</td><td>240 (92.7%)</td></tr><tr><td>Inhaled nitric oxide</td><td>26 (10.0%)</td></tr><tr><td>Extra-corporeal membrane oxygenation</td><td>52 (20.1%)</td></tr><tr><td>Death in ICU</td><td>115 (44%)</td></tr></tbody></table>
<figure><img src="image_3.png" /><figcaption>**Figure1.** Day-60 cumulative probability of ventilator-associated pneumonia (VAP) in patients infected with the variant of origin (red lines), variant α (black lines) or other variants (blue lines). Cumulative incidence estimated using the Kalbfleish and Prentice method considering time from intubation to VAP (dotted line), to death (continuous line) and to weaning (dashed line).</figcaption></figure>

<table><caption>**Table 2.** Association between relevant mutations [(substitutions, deletions (Del))] in spike selected a priori and risk of ventilator-acquired pneumonia by univariable Fine and Gray model analysis. SHR sub-hazard ratio, CI 95% 95% confidence interval; variables are shown as median [quartile 1-quartile 3] and number (percentage). Significant values are in [bold].</caption><thead><tr><th rowspan="3">Mutations</th><th rowspan="3">Missing data</th><th colspan="2">VAP</th><th rowspan="2">Total</th><th rowspan="3">SHR (CI 95%)</th></tr><tr><th>Yes</th><th>No</th></tr><tr><th>N=153</th><th>N=106</th><th>N=259</th></tr></thead><tbody><tr><td>del16970</td><td>0</td><td>61 (39.9%)</td><td>55 (51.9%)</td><td>116 (44.8%)</td><td>0.77 [0.56; 1.06]</td></tr><tr><td>del140145</td><td>0</td><td>61 (39.9%)</td><td>55 (51.9%)</td><td>116 (44.8%)</td><td>0.77 [0.56; 1.06]</td></tr><tr><td>del242244</td><td>0</td><td>8 (5.2%)</td><td>5 (4.7%)</td><td>13 (5.0%)</td><td>1.12 [0.54; 2.30]</td></tr><tr><td>L5F</td><td>0</td><td>4 (2.6%)</td><td>3 (2.8%)</td><td>7 (2.7%)</td><td>0.89 [0.36; 2.20]</td></tr><tr><td>L18F</td><td>0</td><td>4 (2.6%)</td><td>3 (2.8%)</td><td>7 (2.7%)</td><td>1.04 [0.37; 2.94]</td></tr><tr><td><strong>A67V</strong></td><td><strong>0</strong></td><td><strong>1 (0.7%)</strong></td><td><strong>0 (0.0%)</strong></td><td><strong>1 (0.4%)</strong></td><td><strong>5.70 [4.27; 7.60]</strong></td></tr><tr><td>D80AG</td><td>0</td><td>9 (5.9%)</td><td>5 (4.7%)</td><td>14 (5.4%)</td><td>1.26 [0.62; 2.57]</td></tr><tr><td>S98F</td><td>0</td><td>4 (2.6%)</td><td>3 (2.8%)</td><td>7 (2.7%)</td><td>1.04 [0.36; 3.05]</td></tr><tr><td>K417NT</td><td>0</td><td>10 (6.5%)</td><td>5 (4.7%)</td><td>15 (5.8%)</td><td>1.33 [0.69; 2.58]</td></tr><tr><td>N439K</td><td>0</td><td>0 (0.0%)</td><td>0 (0.0%)</td><td>0 (0.0%)</td><td>-</td></tr><tr><td>N450K</td><td>0</td><td>0 (0.0%)</td><td>0 (0.0%)</td><td>0 (0.0%)</td><td>-</td></tr><tr><td>L452R</td><td>0</td><td>3 (2.0%)</td><td>2 (1.9%)</td><td>5 (1.9%)</td><td>1.08 [0.34; 3.43]</td></tr><tr><td>S477N</td><td>0</td><td>38 (24.8%)</td><td>21 (19.8%)</td><td>59 (22.8%)</td><td>1.12 [0.79; 1.57]</td></tr><tr><td>T478K</td><td>0</td><td>0 (0.0%)</td><td>0 (0.0%)</td><td>0 (0.0%)</td><td>-</td></tr><tr><td>E484K</td><td>0</td><td>11 (7.2%)</td><td>5 (4.7%)</td><td>16 (6.2%)</td><td>1.49 [0.77; 2.88]</td></tr><tr><td>N501Y</td><td>0</td><td>71 (46.4%)</td><td>61 (57.5%)</td><td>132 (51.0%)</td><td>0.80 [0.59; 1.09]</td></tr><tr><td>A570D</td><td>0</td><td>60 (39.2%)</td><td>55 (51.9%)</td><td>115 (44.4%)</td><td>0.76 [0.55; 1.04]</td></tr><tr><td>D614G</td><td>0</td><td>149 (97.4%)</td><td>104 (98.1%)</td><td>253 (97.7%)</td><td>0.85 [0.35; 2.06]</td></tr><tr><td>H655Y</td><td>0</td><td>2 (1.3%)</td><td>1 (0.9%)</td><td>3 (1.2%)</td><td>1.32 [0.32; 5.37]</td></tr><tr><td>P681H</td><td>0</td><td>61 (39.9%)</td><td>53 (50.0%)</td><td>114 (44.0%)</td><td>0.81 [0.59; 1.11]</td></tr><tr><td>P681R</td><td>0</td><td>2 (1.3%)</td><td>3 (2.8%)</td><td>5 (1.9%)</td><td>0.57 [0.15; 2.24]</td></tr></tbody></table>

the risk of VAP (p < 0.05), the number of patients harboring this mutation was considered too small (n = 1) to draw clinically relevant conclusions.
<table><caption>**Table 3.** Multivariable Fine and Gray model analysis on risk of ventilator-acquired pneumonia. SHR sub-hazard ratio, CI 95% 95% confidence interval; <sup>a</sup>pre-existing variant or other variants.</caption><thead><tr><td>Mutations</td><td>Missing data</td><td>SHR</td><td>CI 95%</td></tr></thead><tbody><tr><td>Variant α</td><td>0</td><td>0.76</td><td>[ 0.54; 1.07 ]</td></tr><tr><td>Other SARS-CoV-2 variants <sup>a</sup></td><td>0</td><td>1.07</td><td>[ 0.62; 1.87 ]</td></tr><tr><td>Age</td><td>0</td><td>1.02</td><td>[ 1.01; 1.04 ]</td></tr><tr><td>Male sex</td><td>0</td><td>1.39</td><td>[ 0.98; 1.98 ]</td></tr><tr><td>Antibiotic at admission</td><td>0</td><td>0.84</td><td>[ 0.59; 1.18 ]</td></tr><tr><td>SAPS II</td><td>0</td><td>1.00</td><td>[ 0.99; 1.01 ]</td></tr></tbody></table>

## Discussion

We previously failed to find a significant relationship between variant lineages, including variant α, pre-existing and other variants, and mortality<sup>7</sup>. The main results of this study show the lack of statistically significant relationship between variant lineages or the presence of any of 17 relevant spike substitutions and/or deletions selected a priori and VAP occurrence. These findings suggest that the occurrence of VAP is not related to the virological nature of the infecting agent, but to non-virological causes.

The full-length SARS-CoV-2 genomes of 413 critically ill COVID-19 patients from 11 ICUs, who were predominantly infected with pre-existing and α (B.1.1.7) variants, were sequenced and the relationship between viral sequences and VAP occurrence was studied. The main results of this study show the lack of statistically significant relationship between variant lineages or the presence of any of 17 relevant spike substitutions and/or deletions selected a priori and VAP occurrence. Previous large-scale data have suggested that patients infected with variant α (B.1.1.7), identified using the SGTF proxy, had a higher risk of dying<sup>8-10</sup>. Yet, no specific comorbidity or clinical conditions was shown to interact with variant types in the previous studies. Interestingly, following SARS-CoV-2-infected patients with different disease severities through the pathway of disease management, Grint et al. reported a weakening association between variant α (B.1.1.7) infection and mortality, highest in the primary care population, lower in the hospitalized population, not significant in patients admitted in the ICU<sup>10</sup>. Moreover, in another cohort of 496 hospitalized patients in the United Kingdom, 341 samples positive for SARS-CoV-2 on PCR were sequenced and analyzed, including 198 with B.1.1.7 variant (α) and 143 with non-B.1.1.7 variant infection. No association with variant and death was found with or without adjustment for confounding variables<sup>11</sup>.

Together, our results indicate that, in patients with the most severe forms of COVID-19 who required ICU admission, there was no particular mutational pattern associated with VAP. Therefore, the occurrence of VAP is not related to the virological nature of the infecting agent, but to other causes.

**Limitations and strengths of the study.** We acknowledge that our study has some limitations. In this cohort, included patients were predominantly infected with pre-existing and α (B.1.1.7) variants, while the inclusion period ended before the emergence of the δ and o variants, precluding the generalizability of our findings to these more recent variants. The number of patients included was limited and thus the statistical power may have been too weak to show between-group differences. Yet, there was no clear trend regarding associations between variant groups/mutations and VAP, suggesting that increasing the number of patients in the cohort would not have changed the results.

Our study also has strengths, including the constitution of a prospective multicenter cohort of well-phenotyped critically ill patients, and the fact that we performed full-length SARS-CoV-2 genome sequencing in a large number of patients. The start of inclusion of the present study was October 1, 2020, which corresponded in France to the beginning of the second COVID-19 wave, a period when ICU admission policies, ICU patient load and demand, and COVID-19 management strategies were more homogeneous among centres than during the first COVID-19 wave. The high incidence of VAP described in our study is in line with that reported in previous studies<sup>1,12</sup>.

## Conclusions

We found no association between the variant status or any mutational pattern in SARS-CoV-2 viral genes and VAP, indicating that the occurrence of VAP is related to non-virological causes.

## Data availability

The datasets supporting the conclusions are included within the article. The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.

Received: 27 September 2022; Accepted: 16 April 2023

Published online: 24 April 2023
## References

1. Rouzé, A. et al. Relationship between SARS-CoV-2 infection and the incidence of ventilator-associated lower respiratory tract infections: A European multicenter cohort study. *Intensive Care Med* **47**, 188–198. https://doi.org/10.1007/s00134-020-06323-9 (2021).

2. RECOVERY Collaborative Group et al. Dexamethasone in hospitalized patients with Covid-19. *N. Engl. J. Med*. **384**, 693–704. https://doi.org/10.1056/NEJMoa2021436 (2021).

3. Shilts, M. H. et al. Severe COVID-19 is associated with an altered upper respiratory tract microbiome. *Front. Cell. Infect. Microbiol.* **11**, 781968. https://doi.org/10.3389/fcimb.2021.781968 (2021).

4. Leone, M. et al. Brief summary of French guidelines for the prevention, diagnosis and treatment of hospital-acquired pneumonia in ICU. *Ann. Intensive Care* **8**, 104. https://doi.org/10.1186/s13613-018-0444-0 (2018).

5. Fine, J. P. & Gray, R. J. A proportional hazards model for the subdistribution of a competing risk. *J. Am. Stat. Assoc.* **94**, 496–509. https://doi.org/10.1080/01621459.1999.10474144 (1999).

6. Gray, R. J. A class of K-sample tests for comparing the cumulative incidence of a competing risk. *Ann. Stat.* **16**, 1141–1154 (1988).

7. Fourati, S. et al. SARS-CoV-2 genomic characteristics and clinical impact of SARS-CoV-2 viral diversity in critically ill COVID-19 patients: A prospective multicenter cohort study. *Viruses* **14**, 1529. https://doi.org/10.3390/v14071529 (2022).

8. Davies, N. G. et al. Increased mortality in community-tested cases of SARS-CoV-2 lineage B.1.1.7. *Nature* **593**, 270–274. https://doi.org/10.1038/s41586-021-03426-1 (2021).

9. Patone, M. et al. Mortality and critical care unit admission associated with the SARS-CoV-2 lineage B.1.1.7 in England: An observational cohort study. *Lancet Infect. Dis.* **21**, 1518–1528. https://doi.org/10.1016/S1473-3099(21)00318-2 (2021).

10. Grint, D. J. et al. Severity of SARS-CoV-2 alpha variant (B.1.1.7) in England. *Clin. Infect. Dis. Off. Publ. Infect. Dis. Soc. Am.* **20**, 754. https://doi.org/10.1093/cid/ciab754 (2021).

11. Frampton, D. et al. Genomic characteristics and clinical effect of the emergent SARS-CoV-2 B.1.1.7 lineage in London, UK: A whole-genome sequencing and hospital-based cohort study. *Lancet Infect. Dis.* **21**, 1246–1256. https://doi.org/10.1016/S1473-3099(21)00170-5 (2021).

12. Razazi, K. et al. Risks of ventilator-associated pneumonia and invasive pulmonary aspergillosis in patients with viral acute respiratory distress syndrome related or not to Coronavirus 19 disease. *Crit. Care Lond. Engl.* **24**, 699. https://doi.org/10.1186/s13054-020-03417-0 (2020).

## Acknowledgements

The authors wish to thank the patients involved in the study and all the nurses and physicians involved in patient care in the intensive care units;

## Author contributions

K.R., A.M.B. and N.P. had full access to all of the study data and is deemed responsible for data integrity and accuracy of their analysis. K.R., A.M.B. and N.P. contributed to the study initial design, analysis, interpretation of data, drafting of the initial manuscript, critical revision of intellectual content, and approval of the article final version. R.A., B.P., G.V., C.E.L., T.U., J.M., T.P., D.R., R.B., Z.A.H., S.G., M.S., E.A., A.M.D., C.R., J.M.P., S.F. contributed to the study design and analysis, interpretation of data, drafting of initial manuscript, critical revision of intellectual content, and approval of the submitted version of the article.

## Funding

This research received funding from the Agence Nationale de la Recherche (ANR-21-COVR-0022).

## Competing interests

C.-E. L. has served as consultant for Bayer Healthcare, Carmat and Thermo Fisher Brahms, and received lecture fees from MSD, Aerogen, Advanzpharma and BioMérieux, outside the submitted work. J.-M.P. has served as an advisor or speaker for Gilead, Abbvie, Merck, Assembly Biosciences and Arbutus. Other authors have no conflict of interest to disclose.

## Additional information

Correspondence and requests for materials should be addressed to K.R.

Reprints and permissions information is available at www.nature.com/reprints.

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

<figure><img src="image_4.png" /></figure>

**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 Author(s) 2023