Demographic Characteristics
At our tertiary care health center, only 264 HNI individuals were included in the analysis. There were 84 women (31.81%) and 180 men (68.18%). The patient's ages ranged from 9-72, with a mean of 40.62 (standard deviation: 17.58) years. Female patients ranged in age from 9-65 (mean 48.45, SD 18.566) years, whereas male patients ranged in age from 16-72 (mean 49.7, SD 17.1). HNIs were most common in those aged 25–50, followed by those above 50. The submandibular area, parotid area, buccal/oral area, sino-nasal area, otitis media, neck region, and lateral pharyngeal region were targeted for sample collection.
Microbiological Analysis and Bacterial Isolations
In total, 218 (82.58%) positive bacterial isolates were obtained from 264 head and neck infection samples, and the rest 46 (17.42%), did not meet our inclusion criteria. Therefore, 218 cultures were processed for biochemical identification, and 46 samples were excluded from the study. The biochemical analysis identified 8 different types of bacteria, of which the abundance of P. aeruginosa (26.15%) was the highest followed by K. pneumoniae (21.56%), S. aureus (17.89%), A. baumanii (16.51%), E. faecalis (8.26%), S. pyogenes (5.05%), C. freundii (3.21%) and K. oxytoca (1.38%) (Figure1).
Antibiotic susceptibility
The antibiotic susceptibility pattern was analysed, where the resistance pattern was determined according to the zone of inhibition produced by different antibiotics. The susceptibility pattern of phenotypically characterized bacteria was noted differently for gram-negative and gram-positive bacteria. Sixteen different antibiotics were used for both groups of bacteria. Ticarcillin was highly resistant (70.21%) to K. pneumoniae in the case of gram-negative bacteria and the beta-lactam group of antibiotics (Ticarcillin and Piperacillin), followed by monobactam (aztreonam) 65.96% and carbapenem (meropenem) 53.19%. P. aeruginosa revealed the highest resistance to tetracycline at 85.45% and tigecycline at 70.91%. A. baumannii was resistant to monobactam (aztreonam) (76.47%) and 73.53%, whereas C. freundii and K. oxytoca were highly resistant to tetracycline 71.43% and beta-lactam group (ticarcillin) 100%, respectively (Table 1).
According to the antibiotic susceptibility pattern of gram-positive bacteria, the cephalosporin group of antibiotics was highly resistant (97.44%) against S. aureus, amoxicillin had a 93.75% resistant pattern against E. faecalis, and the cyclic-lipopeptide group of antibiotics (daptomycin) had 100% resistance against S. pyogenes (Table 2). According to the statistical analysis, P. aeruginosa (26.15%) is the most resistant isolate, followed by K. oxytocawith the least resistance (1.38%), with a mean and standard deviation 28.94±8.13 and 0.63±0.81 respectively (Table 3).
Bacterial identification results obtained from Vitek 2 procedure remained as usual for the biochemical identification. But the susceptibility pattern revealed higher resistance activity with the same antibiotics in Vitek 2 MIC analysis, where a particular concentration of antibiotics showed their mode of action towards the bacteria (Table 4).
Molecular identification and phylogenetic analysis
The 25 clinical HNI samples were amplified with primers specific to the 16S rRNA gene, to produce 604 bp amplicons (Figure 2). GenBank accession numbers OP881606- OP889581 were assigned to the corresponding sequences (Table 5). The similarities among the existing isolates were found to be between 99-100% in terms of the nucleic acid sequence.
For all four bacteria, the Hd and π ranged from 0.558-0.828 and 0.03236-0.28428, respectively (Supplementary file 1). The S. aureus population had the highest Hd (0.828), followed by K. pneumoniae (0.776), P. aeruginosa (0.598), and A. baumannii (0.558), respectively. In this case, K. pneumoniae (0.28428) has the highest nucleotide diversity followed by A. baumannii (0.14254), S. aureus (0.06180), and P. aeruginosa (0.03236) genes. As the value of selection pressure analysis (θ = dN/dS) < 1, all bacterium of the present study has undergone purifying selection. The negative value of Tajima's D test of neutrality revealed these four bacteria laid on a period of evolutionary expansion. Four species-specific phylogenetic trees were constructed (Figure 3, a-d) utilizing the sequences submitted during this study along with fifty global isolates retrieved from GenBank. Interestingly all the present isolates derived during this study related to Chinese isolates indicating the transboundary potential of these bacteria.