Snapper: high-sensitive detection of methylation motifs based on Oxford Nanopore reads.

MOTIVATION: The Oxford Nanopore technology has a great potential for the analysis of methylated motifs in genomes, including whole-genome methylome profiling. However, we found that there are no methylation motifs detection algorithms, which would be sensitive enough and return deterministic results. Thus, the MEME suit does not extract all Helicobacter pylori methylation sites de novo even using the iterative approach implemented in the most up-to-date methylation analysis tool Nanodisco. RESULTS: We present Snapper, a new highly sensitive approach, to extract methylation motif sequences based on a greedy motif selection algorithm. Snapper does not require manual control during the enrichment process and has enrichment sensitivity higher than MEME coupled with Tombo or Nanodisco instruments that was demonstrated on H.pylori strain J99 studied earlier by the PacBio technology and on four external datasets representing different bacterial species. We used Snapper to characterize the total methylome of a new H.pylori strain A45. At least four methylation sites that have not been described for H.pylori earlier were revealed. We experimentally confirmed the presence of a new CCAG-specific methyltransferase and inferred a gene encoding a new CCAAK-specific methyltransferase. AVAILABILITY AND IMPLEMENTATION: Snapper is implemented using Python and is freely available as a pip package named "snapper-ont." Also, Snapper and the demo dataset are available in Zenodo (10.5281/zenodo.10117651).

16S rRNA gene sequencing data of the upper respiratory tract microbiome in the SARS-CoV-2 infected patients.

The SARS-CoV-2 pandemic is a big challenge for humanity. The COVID-19 severity differs significantly from patient to patient, and it is important to study the factors protecting from severe forms of the disease. Respiratory microbiota may influence the patient's susceptibility to infection and disease severity due to its ability to modulate the immune system response of the host organism. This data article describes the microbiome dataset from the upper respiratory tract of SARS-CoV-2 positive patients from Russia. This dataset reports the microbial community profile of 335 human nasopharyngeal swabs collected between 2020-05 and 2021-03 during the first and the second epidemic waves. Samples were collected from both inpatients and outpatients in 4 cities of the Russian Federation (Moscow, Kazan, Irkutsk, Nizhny Novgorod) and sequenced using the 16S rRNA gene amplicon sequencing of V3-V4 region. Data contains information about the patient such as age, sex, hospitalization status, percent of damaged lung tissue, oxygen saturation (SpO2), respiratory rate, need for supplemental oxygen, chest computer tomography severity score, SARS-CoV-2 lineage, and also information about smoking and comorbidities. The amplicon sequencing data were deposited at NCBI SRA as BioProject PRJNA751478.

Novel Klebsiella pneumoniae K23-Specific Bacteriophages From Different Families: Similarity of Depolymerases and Their Therapeutic Potential.

Antibiotic resistance is a major public health concern in many countries worldwide. The rapid spread of multidrug-resistant (MDR) bacteria is the main driving force for the development of novel non-antibiotic antimicrobials as a therapeutic alternative. Here, we isolated and characterized three virulent bacteriophages that specifically infect and lyse MDR Klebsiella pneumoniae with K23 capsule type. The phages belonged to the Autographiviridae (vB_KpnP_Dlv622) and Myoviridae (vB_KpnM_Seu621, KpS8) families and contained highly similar receptor-binding proteins (RBPs) with polysaccharide depolymerase enzymatic activity. Based on phylogenetic analysis, a similar pattern was also noted for five other groups of depolymerases, specific against capsule types K1, K30/K69, K57, K63, and KN2. The resulting recombinant depolymerases Dep622 (phage vB_KpnP_Dlv622) and DepS8 (phage KpS8) demonstrated narrow specificity against K. pneumoniae with capsule type K23 and were able to protect Galleria mellonella larvae in a model infection with a K. pneumoniae multidrug-resistant strain. These findings expand our knowledge of the diversity of phage depolymerases and provide further evidence that bacteriophages and phage polysaccharide depolymerases represent a promising tool for antimicrobial therapy.

Genome Complexity Browser: Visualization and quantification of genome variability.

Comparative genomics studies may be used to acquire new knowledge regarding genome architecture, which defines the rules for combining sets of genes in the genome of living organisms. Hundreds of thousands of prokaryotic genomes have been sequenced and assembled. However, computational tools capable of simultaneously comparing large numbers of genomes are lacking. We developed the Genome Complexity Browser, a tool that allows the visualization of gene contexts, in a graph-based format, and the quantification of variability for different segments of a genome. The graph-based visualization allows the inspection of changes in gene contents and neighborhoods across hundreds of genomes, simultaneously, which may facilitate the identification of conserved and variable segments of operons or the estimation of the overall variability associated with a particular genome locus. We introduced a measure called complexity, to quantify genome variability. Intraspecies and interspecies comparisons revealed that regions with high complexity values tended to be located in areas that are conserved across different strains and species.

Correction to: Draft genome sequences of Hirudo medicinalis and salivary transcriptome of three closely related medicinal leeches.

An amendment to this paper has been published and can be accessed via the original article.

Draft genome sequences of Hirudo medicinalis and salivary transcriptome of three closely related medicinal leeches.

BACKGROUND: Salivary cell secretion (SCS) plays a critical role in blood feeding by medicinal leeches, making them of use for certain medical purposes even today. RESULTS: We annotated the Hirudo medicinalis genome and performed RNA-seq on salivary cells isolated from three closely related leech species, H. medicinalis, Hirudo orientalis, and Hirudo verbana. Differential expression analysis verified by proteomics identified salivary cell-specific gene expression, many of which encode previously unknown salivary components. However, the genes encoding known anticoagulants have been found to be expressed not only in salivary cells. The function-related analysis of the unique salivary cell genes enabled an update of the concept of interactions between salivary proteins and components of haemostasis. CONCLUSIONS: Here we report a genome draft of Hirudo medicinalis and describe identification of novel salivary proteins and new homologs of genes encoding known anticoagulants in transcriptomes of three medicinal leech species. Our data provide new insights in genetics of blood-feeding lifestyle in leeches.

Shifts in the Human Gut Microbiota Structure Caused by Quadruple Helicobacter pylori Eradication Therapy.

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.

Proteomic dataset: Profiling of cultivated Echerichia coli isolates from Crohn's disease patients and healthy individuals.

One of the dysbioses often observed in Crohn's disease (CD) patients is an increased abundance of Escherichia coli (10-100 fold compared to healthy individuals) (Gevers et al., 2014). The data reported is a large-scale proteome profile for E. coli isolates collected from CD patients and healthy individuals. 43 isolates were achieved from 30 CD patients (17 male, 12 female, median age 30) and 19 isolates from 7 healthy individuals (7 male, median age 19). Isolates were cultivated on LB medium at aerobic conditions up to medium log phase. Protein extraction was performed with sodium deoxycholate (DCNa) and urea, alcylation with tris(2-carboxyethyl)phosphine and iodacetamide. Protein trypsinolysis was performed as described in (Matyushkina et al., 2016). Total cell proteomes were analysed by shotgun proteomics with HPLC-MS/MS on a maXis qTOF mass-spectrometer. The data including HPLC-MS/MS raw files and exported Mascot search results was deposited to the PRIDE repository project accession: PXD010920, project https://doi.org/10.6019/PXD010920.

VERA: agent-based modeling transmission of antibiotic resistance between human pathogens and gut microbiota.

MOTIVATION: The resistance of bacterial pathogens to antibiotics is one of the most important issues of modern health care. The human microbiota can accumulate resistance determinants and transfer them to pathogenic microbiota by means of horizontal gene transfer. Thus, it is important to develop methods of prediction and monitoring of antibiotics resistance in human populations. RESULTS: We present the agent-based VERA model, which allows simulation of the spread of pathogens, including the possible horizontal transfer of resistance determinants from a commensal microbiota community. The model considers the opportunity of residents to stay in the town or in a medical institution, have incorrect self-treatment, treatment with several antibiotics types and transfer and accumulation of resistance determinants from commensal microorganism to a pathogen. In this model, we have also created an assessment of optimum observation frequency of infection spread among the population. Investigating model behavior, we show a number of non-linear dependencies, including the exponential nature of the dependence of the total number of those infected on the average resistance of a pathogen. As the model infection, we chose infection with Shigella spp., though it could be applied to a wide range of other pathogens. AVAILABILITY AND IMPLEMENTATION: Source code and binaries VERA and VERA.viewer are freely available for download at github.com/lpenguin/microbiota-resistome. The code is written in Java, JavaScript and R for Linux platform. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Long-term impact of fecal transplantation in healthy volunteers.

BACKGROUND: Fecal microbiota transplantation (FMT) has been recently approved by FDA for the treatment of refractory recurrent clostridial colitis (rCDI). Success of FTM in treatment of rCDI led to a number of studies investigating the effectiveness of its application in the other gastrointestinal diseases. However, in the majority of studies the effects of FMT were evaluated on the patients with initially altered microbiota. The aim of our study was to estimate effects of FMT on the gut microbiota composition in healthy volunteers and to monitor its long-term outcomes. RESULTS: We have performed a combined analysis of three healthy volunteers before and after capsule FMT by evaluating their general condition, adverse clinical effects, changes of basic laboratory parameters, and several immune markers. Intestinal microbiota samples were evaluated by 16S rRNA gene and shotgun sequencing. The data analysis demonstrated profound shift towards the donor microbiota taxonomic composition in all volunteers. Following FMT, all the volunteers exhibited gut colonization with donor gut bacteria and persistence of this effect for almost ∼1 year of observation. Transient changes of immune parameters were consistent with suppression of T-cell cytotoxicity. FMT was well tolerated with mild gastrointestinal adverse events, however, one volunteer developed a systemic inflammatory response syndrome. CONCLUSIONS: The FMT leads to significant long-term changes of the gut microbiota in healthy volunteers with the shift towards donor microbiota composition and represents a relatively safe procedure to the recipients without long-term adverse events.

Genetic diversity of Escherichia coli in gut microbiota of patients with Crohn's disease discovered using metagenomic and genomic analyses.

BACKGROUND: Crohn's disease is associated with gut dysbiosis. Independent studies have shown an increase in the abundance of certain bacterial species, particularly Escherichia coli with the adherent-invasive pathotype, in the gut. The role of these species in this disease needs to be elucidated. METHODS: We performed a metagenomic study investigating the gut microbiota of patients with Crohn's disease. A metagenomic reconstruction of the consensus genome content of the species was used to assess the genetic variability. RESULTS: The abnormal shifts in the microbial community structures in Crohn's disease were heterogeneous among the patients. The metagenomic data suggested the existence of multiple E. coli strains within individual patients. We discovered that the genetic diversity of the species was high and that only a few samples manifested similarity to the adherent-invasive varieties. The other species demonstrated genetic diversity comparable to that observed in the healthy subjects. Our results were supported by a comparison of the sequenced genomes of isolates from the same microbiota samples and a meta-analysis of published gut metagenomes. CONCLUSIONS: The genomic diversity of Crohn's disease-associated E. coli within and among the patients paves the way towards an understanding of the microbial mechanisms underlying the onset and progression of the Crohn's disease and the development of new strategies for the prevention and treatment of this disease.

Ultrahigh-throughput functional profiling of microbiota communities.

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.

Draft genomes of Enterococcus faecium strains isolated from human feces before and after eradication therapy against Helicobacter pylori.

The abundance of Enterococci in the human intestinal microbiota environment is usually < 0.1% of the total bacterial fraction. The multiple resistance to antibiotics of the opportunistic Enterococcus spp. is alarming for the world medical community because of their high prevalence among clinically significant strains of microorganisms. Enterococci are able to collect different mobile genetic elements and transmit resistance to antibiotics to wide range of Gram-positive and Gram-negative species of microorganisms, including the transmission of vancomycin resistance to methicillin-resistant strains of Staphylococcus aureus. The number of infections caused by antibiotics resistant strains of Enterococcus spp. is increasing. Here we present a draft genomes of Enterococcus faecium strains. These strains were isolated from human feces before and after (1 month) Helicobacter pylori eradication therapy. The samples were subject to whole-genome sequencing using Illumina HiSeq. 2500 platform. The data is available at NCBI https://www.ncbi.nlm.nih.gov/bioproject/PRJNA412824.

Identification of unusual peptides with new Cys frameworks in the venom of the cold-water sea anemone Cnidopus japonicus.

Sea anemones (Actiniaria) are intensely popular objects of study in venomics. Order Actiniaria includes more than 1,000 species, thus presenting almost unlimited opportunities for the discovery of novel biologically active molecules. The venoms of cold-water sea anemones are studied far less than the venoms of tropical sea anemones. In this work, we analysed the molecular venom composition of the cold-water sea anemone Cnidopus japonicus. Two sets of NGS data from two species revealed molecules belonging to a variety of structural classes, including neurotoxins, toxin-like molecules, linear polypeptides (Cys-free), enzymes, and cytolytics. High-throughput proteomic analyses identified 27 compounds that were present in the venoms. Some of the toxin-like polypeptides exhibited novel Cys frameworks. To characterise their function in the venom, we heterologously expressed 3 polypeptides with unusual Cys frameworks (designated CjTL7, CjTL8, and AnmTx Cj 1c-1) in E. coli. Toxicity tests revealed that the CjTL8 polypeptide displays strong crustacean-specific toxicity, while AnmTx Cj 1c-1 is toxic to both crustaceans and insects. Thus, an improved NGS data analysis algorithm assisted in the identification of toxins with unusual Cys frameworks showing no homology according to BLAST. Our study shows the advantage of combining omics analysis with functional tests for active polypeptide discovery.

Genome analysis of E. coli isolated from Crohn's disease patients.

BACKGROUND: Escherichia coli (E. coli) has been increasingly implicated in the pathogenesis of Crohn's disease (CD). The phylogeny of E. coli isolated from Crohn's disease patients (CDEC) was controversial, and while genotyping results suggested heterogeneity, the sequenced strains of E. coli from CD patients were closely related. RESULTS: We performed the shotgun genome sequencing of 28 E. coli isolates from ten CD patients and compared genomes from these isolates with already published genomes of CD strains and other pathogenic and non-pathogenic strains. CDEC was shown to belong to A, B1, B2 and D phylogenetic groups. The plasmid and several operons from the reference CD-associated E. coli strain LF82 were demonstrated to be more often present in CDEC genomes belonging to different phylogenetic groups than in genomes of commensal strains. The operons include carbon-source induced invasion GimA island, prophage I, iron uptake operons I and II, capsular assembly pathogenetic island IV and propanediol and galactitol utilization operons. CONCLUSIONS: Our findings suggest that CDEC are phylogenetically diverse. However, some strains isolated from independent sources possess highly similar chromosome or plasmids. Though no CD-specific genes or functional domains were present in all CD-associated strains, some genes and operons are more often found in the genomes of CDEC than in commensal E. coli. They are principally linked to gut colonization and utilization of propanediol and other sugar alcohols.

Microfluidic droplet platform for ultrahigh-throughput single-cell screening of biodiversity.

Ultrahigh-throughput screening (uHTS) techniques can identify unique functionality from millions of variants. To mimic the natural selection mechanisms that occur by compartmentalization in vivo, we developed a technique based on single-cell encapsulation in droplets of a monodisperse microfluidic double water-in-oil-in-water emulsion (MDE). Biocompatible MDE enables in-droplet cultivation of different living species. The combination of droplet-generating machinery with FACS followed by next-generation sequencing and liquid chromatography-mass spectrometry analysis of the secretomes of encapsulated organisms yielded detailed genotype/phenotype descriptions. This platform was probed with uHTS for biocatalysts anchored to yeast with enrichment close to the theoretically calculated limit and cell-to-cell interactions. MDE-FACS allowed the identification of human butyrylcholinesterase mutants that undergo self-reactivation after inhibition by the organophosphorus agent paraoxon. The versatility of the platform allowed the identification of bacteria, including slow-growing oral microbiota species that suppress the growth of a common pathogen, Staphylococcus aureus, and predicted which genera were associated with inhibitory activity.

Specific pools of endogenous peptides are present in gametophore, protonema, and protoplast cells of the moss Physcomitrella patens.

BACKGROUND: Protein degradation is a basic cell process that operates in general protein turnover or to produce bioactive peptides. However, very little is known about the qualitative and quantitative composition of a plant cell peptidome, the actual result of this degradation. In this study we comprehensively analyzed a plant cell peptidome and systematically analyzed the peptide generation process. RESULTS: We thoroughly analyzed native peptide pools of Physcomitrella patens moss in two developmental stages as well as in protoplasts. Peptidomic analysis was supplemented by transcriptional profiling and quantitative analysis of precursor proteins. In total, over 20,000 unique endogenous peptides, ranging in size from 5 to 78 amino acid residues, were identified. We showed that in both the protonema and protoplast states, plastid proteins served as the main source of peptides and that their major fraction formed outside of chloroplasts. However, in general, the composition of peptide pools was very different between these cell types. In gametophores, stress-related proteins, e.g., late embryogenesis abundant proteins, were among the most productive precursors. The Driselase-mediated protonema conversion to protoplasts led to a peptide generation "burst", with a several-fold increase in the number of components in the latter. Degradation of plastid proteins in protoplasts was accompanied by suppression of photosynthetic activity. CONCLUSION: We suggest that peptide pools in plant cells are not merely a product of waste protein degradation, but may serve as important functional components for plant metabolism. We assume that the peptide "burst" is a form of biotic stress response that might produce peptides with antimicrobial activity from originally functional proteins. Potential functions of peptides in different developmental stages are discussed.