SARS-CoV-2 genetic variation and bacterial communities of naso-oropharyngeal samples in middle-aged and elderly COVID-19 patients in West Java, Indonesia

OBJECTIVE: The number of COVID-19 cases in Indonesia reflects the disease severity and rapid dissemination. In response to the mounting threat, SARS-CoV-2 genomic surveillance and the investigation of naso-oropharyngeal bacterial communities in West Java were conducted, as dysbiosis of the upper respiratory tract microbiota might adversely affect the clinical condition of patients. METHODS: We utilized the Oxford Nanopore sequencing platform to analyze genetic variation of 43 samples of SARS-CoV-2 and 11 selected samples for 16S rRNA gene sequencing, using samples collected from May to August 2021. RESULTS: The prevalence of AY.23 (>82%) predominated among five virus lineages in the populations (AY.23, AY.24, AY.26, AY.42, B.1.1.7). The region in the SARS-CoV-2 genome found to have the highest number of mutations was the spike (S) protein (>20%). There was no association between SARS-CoV-2 lineages, mutation frequency, patient profile, and COVID-19 rapid spread-categorized cases. There was no association of bacterial relative abundance, alpha-beta diversity, and linear discriminant analysis effect size analysis with patient profile and rapid spread cases. MetagenomeSeq analysis showed eight differential abundance species in individual patient profiles, including Pseudomonas aeruginosa and Haemophilus parainfluenzae. CONCLUSIONS: The data demonstrated relevant AY.23 dominance (the Delta variant) in West Java during that period supporting the importance of surveillance program in monitoring disease progression. The inconsistent results of the bacterial communities suggest that a complex multifactor process may contribute to the progression of bacterial-induced disease in each patient.