Abstract
The Omicron variant of the severe acute respiratory syndrome coronavirus 2 (SARS‑CoV‑2) infected a substantial proportion of Chinese population, and understanding the factors underlying the severity of the disease and fatality is valuable for future prevention and clinical treatment. We recruited 64 patients with invasive ventilation for COVID-19 and performed metatranscriptomic sequencing to profile host transcriptomic profiles, plus viral, bacterial, and fungal content, as well as virulence factors and examined their relationships to 28-day mortality were examined. In addition, the bronchoalveolar lavage fluid (BALF) samples from invasive ventilated hospital/community-acquired pneumonia patients (HAP/CAP) sampled in 2019 were included for comparison. Genomic analysis revealed that all Omicron strains belong to BA.5 and BF.7 sub-lineages, with no difference in 28-day mortality between them. Compared to HAP/CAP cohort, invasive ventilated COVID-19 patients have distinct host transcriptomic and microbial signatures in the lower respiratory tract; and in the COVID-19 non-survivors, we found significantly lower gene expressions in pathways related viral processes and positive regulation of protein localization to plasma membrane, higher abundance of opportunistic pathogens including bacterial Alloprevotella, Caulobacter, Escherichia-Shigella, Ralstonia and fungal Aspergillus sydowii and Penicillium rubens. Correlational analysis further revealed significant associations between host immune responses and microbial compositions, besides synergy within viral, bacterial, and fungal pathogens. Our study presents the relationships of lower respiratory tract microbiome and transcriptome in invasive ventilated COVID-19 patients, providing the basis for future clinical treatment and reduction of fatality.