Microbiome spectra serve as critical clues to elucidate the evolutionary biology pathways, potential pathologies, and even behavioral patterns of the host organisms. Furthermore, exotic sources of microbiota represent an unexplored niche to discover microbial secondary metabolites. However, establishing the bacterial functionality is complicated by an intricate web of interactions inside the microbiome. Here we apply an ultrahigh-throughput (uHT) microfluidic droplet platform for activity profiling of the entire oral microbial community of the Siberian bear to isolate Bacillus strains demonstrating antimicrobial activity against Staphylococcus aureus . Genome mining allowed us to identify antibiotic amicoumacin A (Ami) as responsible for inhibiting the growth of S. aureus . Proteomics and metabolomics revealed a unique mechanism of Bacillus self-resistance to Ami, based on a subtle equilibrium of its deactivation and activation by kinase AmiN and phosphatase AmiO, respectively. We developed uHT quantitative single-cell analysis to estimate antibiotic efficacy toward different microbiomes and used it to determine the activity spectra of Ami toward human and Siberian bear microbiota. Thus, uHT microfluidic droplet platform activity profiling is a powerful tool for discovering antibiotics and quantifying external influences on a microbiome.

Klebsiella pneumoniae is one of the most important pathogens concerned with multidrug resistance in healthcare-associated infections. The treating of infections caused by this bacterium is complicated due to the emergence and rapid spreading of carbapenem-resistant strains, which are associated with high mortality rates. Recently, several hypervirulent and carbapenemase-producing isolates were reported that make the situation even more complicated. In order to better understand the resistance and virulence mechanisms, and, in turn, to develop effective treatment strategies for the infections caused by multidrug-resistant K. pneumoniae, more comprehensive genomic and phenotypic data are required. Here, we present the first detailed molecular epidemiology report based on second and third generation (long-read) sequencing for the clinical isolates of K. pneumoniae in the Russian Federation. The data include three schemes of molecular typing, phenotypic and genotypic antibiotic resistance determination, as well as the virulence and plasmid profiling for 36 K. pneumoniae isolates. We have revealed 2 new multilocus sequence typing (MLST)-based sequence types, 32 multidrug-resistant (MDR) isolates and 5 colistin-resistant isolates in our samples. Three MDR isolates belonged to a very rare ST377 type. The whole genome sequences and additional data obtained will greatly facilitate further investigations in the field of antimicrobial resistance studies.

The human gut microbiome plays an important role both in health and disease. Use of antibiotics can alter gut microbiota composition, which can lead to various deleterious events. Here we report a whole genome sequencing metagenomic/genomic study of the intestinal microbiota changes caused by Helicobacter pylori (HP) eradication therapy. Using approaches for metagenomic data analysis we revealed a statistically significant decrease in alpha-diversity and relative abundance of Bifidobacterium adolescentis due to HP eradication therapy, while the relative abundance of Enterococcus faecium increased. We have detected changes in general metagenome resistome profiles as well: after HP eradication therapy, the ermB, CFX group, and tetQ genes were overrepresented, while tetO and tetW genes were underrepresented. We have confirmed these results with genome-resolved metagenomic approaches. MAG (metagenome-assembled genomes) abundance profiles have changed dramatically after HP eradication therapy. Focusing on ermB gene conferring resistance to macrolides, which were included in the HP eradication therapy scheme, we have shown a connection between antibiotic resistance genes (ARGs) and some overrepresented MAGs. Moreover, some E. faecium strains isolated from stool samples obtained after HP eradication have manifested greater antibiotic resistance in vitro in comparison to other isolates, as well as the higher number of ARGs conferring resistance to macrolides and tetracyclines.

Background: The dual infection with SARS-CoV-2 is poorly described and is currently under discussion. We present a study of two strains of SARS-CoV-2 detected in the same patient during the same disease presentation. Case presentation: A patient in their 90 s was hospitalised with fever. Oropharyngeal swab obtained on the next day (sample 1) tested positive for SARS-CoV-2. Five days later, the patient was transferred to the ICU (intensive care unit) of the hospital specialising in the treatment of COVID-19 patients, where the patient's condition progressively worsened and continuous oxygen insufflation was required. Repeated oropharyngeal swab (sample 2), which was taken eight days after the first one, also tested positive for SARS-CoV-2. After 5 days of ICU treatment, the patient died. The cause of death was a coronavirus infection, which progressed unfavourably due to premorbid status. We have performed sequencing of full SARS-CoV-2 genomes from oropharyngeal swabs obtained eight days apart. Genomic analysis revealed the presence of two genetically distant SARS-CoV-2 strains in both swabs. Detected strains belong to different phylogenetic clades (GH and GR) and differ in seven nucleotide positions. The relative abundance of strains was 70% (GH) and 30% (GR) in the first swab, and 3% (GH) and 97% (GR) in the second swab. Conclusions: Our findings suggest that the patient was infected by two genetically distinct SARS-CoV-2 strains at the same time. One of the possible explanations is that the second infection was hospital-acquired. Change of the dominant strain ratio during disease manifestation could be explained by the advantage or higher virulence of the GR clade strain.

Being diverse and widely distributed globally, bats are a known reservoir of a series of emerging zoonotic viruses. We studied fecal viromes of twenty-six bats captured in 2015 in the Moscow Region and found 13 of 26 (50%) samples to be coronavirus positive. Of P. nathusii (the Nathusius’ pipistrelle), 3 of 6 samples were carriers of a novel MERS-related betacoronavirus. We sequenced and assembled the complete genome of this betacoronavirus and named it MOW-BatCoV strain 15-22. Whole genome phylogenetic analysis suggests that MOW-BatCoV/15-22 falls into a distinct subclade closely related to human and camel MERS-CoV. Unexpectedly, the phylogenetic analysis of the novel MOW-BatCoV/15-22 spike gene showed the closest similarity to CoVs from Erinaceus europaeus (European hedgehog). We suppose MOW-BatCoV could have arisen as a result of recombination between ancestral viruses of bats and hedgehogs. Molecular docking analysis of MOW-BatCoV/15-22 spike glycoprotein binding to DPP4 receptors of different mammals predicted the highest binding ability with DPP4 of the Myotis brandtii bat (docking score −320.15) and the E. europaeus (docking score –294.51). Hedgehogs are widely kept as pets and are commonly found in areas of human habitation. As this novel bat-CoV is likely capable of infecting hedgehogs, we suggest hedgehogs can act as intermediate hosts between bats and humans for other bat-CoVs.

Members of the jingmenviruses group have been found in arthropods and mammals on all continents except Australia and Antarctica. Two viruses of this group were isolated from patients with fever after a tick bite. Using a nested RT-PCR assay targeting a jingmenvirus polymerase gene fragment, we screened ticks collected in seven regions of Russia and found that the abundant jingmenvirus-positive were of Ixodes ricinus species, with the prevalence ranging from 19.8% to 34.3%. In all cases, DNase/RNase treatment suggested that the detected molecule was DNA and subsequent next generation sequencing (NGS) proved that the viral polymerase gene was integrated in the I. ricinus genome. The copy number of the integrated polymerase gene was quantified by qPCR relative to the ITS2 gene and estimated as 1.32 copies per cell. At least three different genetic variants of the integrated polymerase gene were found in the territory of Russia. Phylogenetic analysis of the integrated jingmenvirus polymerase gene showed the highest similarity with the sequence of the correspondent gene obtained in Serbia from I. ricinus.