Non-human peptides revealed in blood reflect the composition of intestinal microbiota.

BACKGROUND: The previously underestimated effects of commensal gut microbiota on the human body are increasingly being investigated using omics. The discovery of active molecules of interaction between the microbiota and the host may be an important step towards elucidating the mechanisms of symbiosis. RESULTS: Here, we show that in the bloodstream of healthy people, there are over 900 peptides that are fragments of proteins from microorganisms which naturally inhabit human biotopes, including the intestinal microbiota. Absolute quantitation by multiple reaction monitoring has confirmed the presence of bacterial peptides in the blood plasma and serum in the range of approximately 0.1 nM to 1 µM. The abundance of microbiota peptides reaches its maximum about 5 h after a meal. Most of the peptides correlate with the bacterial composition of the small intestine and are likely obtained by hydrolysis of membrane proteins with trypsin, chymotrypsin and pepsin - the main proteases of the gastrointestinal tract. The peptides have physicochemical properties that likely allow them to selectively pass the intestinal mucosal barrier and resist fibrinolysis. CONCLUSIONS: The proposed approach to the identification of microbiota peptides in the blood, after additional validation, may be useful for determining the microbiota composition of hard-to-reach intestinal areas and monitoring the permeability of the intestinal mucosal barrier.

Identification and Genetic Characterization of MERS-Related Coronavirus Isolated from Nathusius' Pipistrelle (Pipistrellus nathusii) near Zvenigorod (Moscow Region, Russia).

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.

Consistent Stool Metagenomic Biomarkers Associated with the Response To Melanoma Immunotherapy.

The human gut microbiome plays an important role in both health and disease. Recent studies have demonstrated a strong influence of the gut microbiome composition on the efficacy of cancer immunotherapy. However, available studies have not yet succeeded in finding reliable and consistent metagenomic markers that are associated with the response to immunotherapy. Therefore, the reanalysis of the published data may improve our understanding of the association between the composition of the gut microbiome and the treatment response. In this study, we focused on melanoma-related metagenomic data, which are more abundant than are data from other tumor types. We analyzed the metagenomes of 680 stool samples from 7 studies that were published earlier. The taxonomic and functional biomarkers were selected after comparing the metagenomes of patients showing different treatment responses. The list of selected biomarkers was also validated on additional metagenomic data sets that were dedicated to the influence of fecal microbiota transplantation on the response to melanoma immunotherapy. According to our analysis, the resulting cross-study taxonomic biomarkers included three bacterial species: Faecalibacterium prausnitzii, Bifidobacterium adolescentis, and Eubacterium rectale. 101 groups of genes were identified to be functional biomarkers, including those potentially involved in the production of immune-stimulating molecules and metabolites. Moreover, we ranked the microbial species by the number of genes encoding functionally relevant biomarkers that they contained. Thus, we put together a list of potentially the most beneficial bacteria for immunotherapy success. F. prausnitzii, E. rectale, and three species of bifidobacteria stood out as the most beneficial species, even though some useful functions were also present in other bacterial species. IMPORTANCE In this study, we put together a list of potentially the most beneficial bacteria that were associated with a responsiveness to melanoma immunotherapy. Another important result of this study is the list of functional biomarkers of responsiveness to immunotherapy, which are dispersed among different bacterial species. This result possibly explains the existing irregularities between studies regarding the bacterial species that are beneficial to melanoma immunotherapy. Overall, these findings can be utilized to issue recommendations for gut microbiome correction in cancer immunotherapy, and the resulting list of biomarkers might serve as a good stepping stone for the development of a diagnostic test that is aimed at predicting patients' responses to melanoma immunotherapy.

UniqPy: A tool for estimation of short-chain fatty acids composition by gas-chromatography/mass-spectrometry with headspace extraction.

Short-chain fatty acids are metabolites widely presented in many natural sources, including human feces and blood. Estimation of their composition is a common procedure, usually performed using nuclear magnetic resonance or gas chromatography with a flame ionization detector. However, the commonly used methods often depend on specific sample preparation, such as filtration and homogenization. The gas-chromatography/mass-spectrometry (GC/MS) method with headspace extraction allows sample preparation to be kept to a minimum regardless of the physical state of the sample, which can be potentially useful in metabolomics research of complex natural samples such as blood or feces. In this work, we have demonstrated the applicability of Headspace GC-MS for estimating short chain fatty acid (SCFA) composition. The main problem here is the complex, non-linear dependence between the composition of the compounds in the source phase and the relative pressures in the vapor phase, which are directly measured by this method. We have implemented a thermodynamic model that performs the reverse transformation of relative abundances in the vapor phase to relative concentrations in the liquid phase, and have tested it on some synthetic SCFA mixtures. The developed method is available as a pip package called UniqPy and can be used to describe liquid-vapor equilibrium for any multicomponent system if a sufficient amount of training data is provided. The gas chromatography method with headspace extraction in conjunction with the UniqPy data transformation showed satisfactory quantification accuracy for propionic acid, butyric acid, isobutyric acid, and valeric acid (R-squared > 0.96). The applicability of the method was additionally demonstrated on a series of fecal samples.

The effects of Levilactobacillus brevis on the physiological parameters and gut microbiota composition of rats subjected to desynchronosis.

BACKGROUND: All living organisms have developed during evolution complex time-keeping biological clocks that allowed them to stay attuned to their environments. Circadian rhythms cycle on a near 24 h clock. These encompass a variety of changes in the body ranging from blood hormone levels to metabolism, to the gut microbiota composition and others. The gut microbiota, in return, influences the host stress response and the physiological changes associated with it, which makes it an important determinant of health. Lactobacilli are traditionally consumed for their prophylactic and therapeutic benefits against various diseases, namely, the inflammatory bowel syndrome, and even emerged recently as promising psychobiotics. However, the potential role of lactobacilli in the normalization of circadian rhythms has not been addressed. RESULTS: Two-month-old male rats were randomly divided into three groups and housed under three different light/dark cycles for three months: natural light, constant light and constant darkness. The strain Levilactobacillus brevis 47f was administered to rats at a dose of 0.5 ml per rat for one month and The rats were observed for the following two months. As a result, we identified the biomarkers associated with intake of L. brevis 47f. Changing the light regime for three months depleted the reserves of the main buffer in the cell-reduced glutathione. Intake of L. brevis 47f for 30 days restored cellular reserves of reduced glutathione and promoted redox balance. Our results indicate that the levels of urinary catecholamines correlated with light/dark cycles and were influenced by intake of L. brevis 47f. The gut microbiota of rats was also influenced by these factors. L. brevis 47f intake was associated with an increase in the relative abundance of Faecalibacterium and Roseburia and a decrease in the relative abundance of Prevotella and Bacteroides. CONCLUSIONS: The results of this study show that oral administration of L. brevis 47f, for one month, to rats housed under abnormal lightning conditions (constant light or constant darkness) normalized their physiological parameters and promoted the gut microbiome's balance.

Gene Networks Underlying the Resistance of Bifidobacterium longum to Inflammatory Factors.

As permanent residents of the normal gut microbiota, bifidobacteria have evolved to adapt to the host's immune response whose priority is to eliminate pathogenic agents. The mechanisms that ensure the survival of commensals during inflammation and maintain the stability of the core component of the normal gut microbiota in such conditions remain poorly understood. We propose a new in vitro approach to study the mechanisms of resistance to immune response factors based on high-throughput sequencing followed by transcriptome analysis. This approach allowed us to detect differentially expressed genes associated with inflammation. In this study, we demonstrated that the presence of the pro-inflammatory cytokines IL-6 and TNFα to the growth medium of the B. longum subsp. longum GT15 strain changes the latter's growth rate insignificantly while affecting the expression of certain genes. We identified these genes and performed a COG and a KEGG pathway enrichment analysis. Using phylogenetic profiling we predicted the operons of genes whose expression was triggered by the cytokines TNFα and IL-6 in vitro. By mapping the transcription start points, we experimentally validated the predicted operons. Thus, in this study, we predicted the genes involved in a putative signaling pathway underlying the mechanisms of resistance to inflammatory factors in bifidobacteria. Since bifidobacteria are a major component of the human intestinal microbiota exhibiting pronounced anti-inflammatory properties, this study is of great practical and scientific relevance.

Identification of Clinically Significant Prostate Cancer by Combined PCA3 and AMACR mRNA Detection in Urine Samples.

PURPOSE: Preclinical evaluation of PCA3 and AMACR transcript simultaneous detection in urine to diagnose clinical significant prostate cancer (prostate cancer with Gleason score ≥7) in a Russian cohort. PATIENTS AND METHODS: We analyzed urine samples of patients with a total serum PSA ≥2 ng/mL: 31 men with prostate cancer scheduled for radical prostatectomy, 128 men scheduled for first diagnostic biopsy (prebiopsy cohort). PCA3, AMACR, PSA and GPI transcripts were detected by multiplex reverse transcription quantitative polymerase chain reaction, and the results were used for scores for calculation and statistical analysis. RESULTS: There was no significant difference between clinically significant and nonsignificant prostate cancer PCA3 scores. However, there was a significant difference in the AMACR score (patients scheduled for radical prostatectomy p=0.0088, prebiopsy cohort p=0.029). We estimated AUCs, optimal cutoffs, sensitivities and specificities for PCa and csPCa detection in the prebiopsy cohort by tPSA, PCA3 score, PCPT Risk Calculator and classification models based on tPSA, PCA3 score and AMACR score. In the clinically significant prostate cancer ROC analysis, the PCA3 score AUC was 0.632 (95%CI: 0.511-0.752), the AMACR score AUC was 0.711 (95%CI: 0.617-0.806) and AUC of classification model based on the PCA3 score, the AMACR score and total PSA was 0.72 (95%CI: 0.58-0.83). In addition, the correlation of the AMACR score with the ratio of total RNA and RNA of prostate cells in urine was shown (tau=0.347, p=6.542e-09). Significant amounts of nonprostate RNA in urine may be a limitation for the AMACR score use. CONCLUSION: The AMACR score is a good predictor of clinically significant prostate cancer. Significant amounts of nonprostate RNA in urine may be a limitation for the AMACR score use. Evaluation of the AMACR score and classification models based on it for clinically significant prostate cancer detection with larger samples and a follow-up analysis is promising.

Deep Functional Profiling Facilitates the Evaluation of the Antibacterial Potential of the Antibiotic Amicoumacin.

The global spread of antibiotic resistance is forcing the scientific community to find new molecular strategies to counteract it. Deep functional profiling of microbiomes provides an alternative source for the discovery of novel antibiotic producers and probiotics. Recently, we implemented this ultrahigh-throughput screening approach for the isolation of Bacillus pumilus strains efficiently producing the ribosome-targeting antibiotic amicoumacin A (Ami). Proteomics and metabolomics revealed essential insight into the activation of Ami biosynthesis. Here, we applied omics to boost Ami biosynthesis, providing the optimized cultivation conditions for high-scale production of Ami. Ami displayed a pronounced activity against Lactobacillales and Staphylococcaceae, including methicillin-resistant Staphylococcus aureus (MRSA) strains, which was determined using both classical and massive single-cell microfluidic assays. However, the practical application of Ami is limited by its high cytotoxicity and particularly low stability. The former is associated with its self-lactonization, serving as an improvised intermediate state of Ami hydrolysis. This intramolecular reaction decreases Ami half-life at physiological conditions to less than 2 h, which is unprecedented for a terminal amide. While we speculate that the instability of Ami is essential for Bacillus ecology, we believe that its stable analogs represent attractive lead compounds both for antibiotic discovery and for anticancer drug development.

Links of gut microbiota composition with alcohol dependence syndrome and alcoholic liver disease.

BACKGROUND: Alcohol abuse has deleterious effects on human health by disrupting the functions of many organs and systems. Gut microbiota has been implicated in the pathogenesis of alcohol-related liver diseases, with its composition manifesting expressed dysbiosis in patients suffering from alcoholic dependence. Due to its inherent plasticity, gut microbiota is an important target for prevention and treatment of these diseases. Identification of the impact of alcohol abuse with associated psychiatric symptoms on the gut community structure is confounded by the liver dysfunction. In order to differentiate the effects of these two factors, we conducted a comparative "shotgun" metagenomic survey of 99 patients with the alcohol dependence syndrome represented by two cohorts-with and without liver cirrhosis. The taxonomic and functional composition of the gut microbiota was subjected to a multifactor analysis including comparison with the external control group. RESULTS: Alcoholic dependence and liver cirrhosis were associated with profound shifts in gut community structures and metabolic potential across the patients. The specific effects on species-level community composition were remarkably different between cohorts with and without liver cirrhosis. In both cases, the commensal microbiota was found to be depleted. Alcoholic dependence was inversely associated with the levels of butyrate-producing species from the Clostridiales order, while the cirrhosis-with multiple members of the Bacteroidales order. The opportunist pathogens linked to alcoholic dependence included pro-inflammatory Enterobacteriaceae, while the hallmarks of cirrhosis included an increase of oral microbes in the gut and more frequent occurrence of abnormal community structures. Interestingly, each of the two factors was associated with the expressed enrichment in many Bifidobacterium and Lactobacillus-but the exact set of the species was different between alcoholic dependence and liver cirrhosis. At the level of functional potential, the patients showed different patterns of increase in functions related to alcohol metabolism and virulence factors, as well as pathways related to inflammation. CONCLUSIONS: Multiple shifts in the community structure and metabolic potential suggest strong negative influence of alcohol dependence and associated liver dysfunction on gut microbiota. The identified differences in patterns of impact between these two factors are important for planning of personalized treatment and prevention of these pathologies via microbiota modulation. Particularly, the expansion of Bifidobacterium and Lactobacillus suggests that probiotic interventions for patients with alcohol-related disorders using representatives of the same taxa should be considered with caution. Taxonomic and functional analysis shows an increased propensity of the gut microbiota to synthesis of the toxic acetaldehyde, suggesting higher risk of colorectal cancer and other pathologies in alcoholics.

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.

State of the Art of Chromosome 18-Centric HPP in 2016: Transcriptome and Proteome Profiling of Liver Tissue and HepG2 Cells.

A gene-centric approach was applied for a large-scale study of expression products of a single chromosome. Transcriptome profiling of liver tissue and HepG2 cell line was independently performed using two RNA-Seq platforms (SOLiD and Illumina) and also by Droplet Digital PCR (ddPCR) and quantitative RT-PCR. Proteome profiling was performed using shotgun LC-MS/MS as well as selected reaction monitoring with stable isotope-labeled standards (SRM/SIS) for liver tissue and HepG2 cells. On the basis of SRM/SIS measurements, protein copy numbers were estimated for the Chromosome 18 (Chr 18) encoded proteins in the selected types of biological material. These values were compared with expression levels of corresponding mRNA. As a result, we obtained information about 158 and 142 transcripts for HepG2 cell line and liver tissue, respectively. SRM/SIS measurements and shotgun LC-MS/MS allowed us to detect 91 Chr 18-encoded proteins in total, while an intersection between the HepG2 cell line and liver tissue proteomes was ∼66%. In total, there were 16 proteins specifically observed in HepG2 cell line, while 15 proteins were found solely in the liver tissue. Comparison between proteome and transcriptome revealed a poor correlation (R2 ≈ 0.1) between corresponding mRNA and protein expression levels. The SRM and shotgun data sets (obtained during 2015-2016) are available in PASSEL (PASS00697) and ProteomeExchange/PRIDE (PXD004407). All measurements were also uploaded into the in-house Chr 18 Knowledgebase at http://kb18.ru/protein/matrix/416126 .

RNA-Seq gene expression profiling of HepG2 cells: the influence of experimental factors and comparison with liver tissue.

BACKGROUND: Human hepatoma HepG2 cells are used as an in vitro model of the human liver. High-throughput transcriptomic sequencing is an advanced approach for assessing the functional state of a tissue or cell type. However, the influence of experimental factors, such as the sample preparation method and inter-laboratory variation, on the transcriptomic profile has not been evaluated. RESULTS: The whole-transcriptome sequencing of HepG2 cells was performed using the SOLiD platform and validated using droplet digital PCR. The gene expression profile was compared to the results obtained with the same sequencing method in another laboratory and using another sample preparation method. We also compared the transcriptomic profile HepG2 cells with that of liver tissue. Comparison of the gene expression profiles between the HepG2 cell line and liver tissue revealed the highest variation, followed by HepG2 cells submitted to two different sample preparation protocols. The lowest variation was observed between HepG2 cells prepared by two different laboratories using the same protocol. The enrichment analysis of the genes that were differentially expressed between HepG2 cells and liver tissue mainly revealed the cancer-associated gene signature of HepG2 cells and the activation of the response to chemical stimuli in the liver tissue. The HepG2 transcriptome obtained with the SOLiD platform was highly correlated with the published transcriptome obtained with the Illumina and Helicos platforms, with moderate correspondence to microarrays. CONCLUSIONS: In the present study, we assessed the influence of experimental factors on the HepG2 transcriptome and identified differences in gene expression between the HepG2 cell line and liver cells. These findings will facilitate robust experimental design in the fields of pharmacology and toxicology. Our results were supported by a comparative analysis with previous HepG2 gene expression studies.

Chromosome 18 transcriptoproteome of liver tissue and HepG2 cells and targeted proteome mapping in depleted plasma: update 2013.

We report the results obtained in 2012-2013 by the Russian Consortium for the Chromosome-centric Human Proteome Project (C-HPP). The main scope of this work was the transcriptome profiling of genes on human chromosome 18 (Chr 18), as well as their encoded proteome, from three types of biomaterials: liver tissue, the hepatocellular carcinoma-derived cell line HepG2, and blood plasma. The transcriptome profiling for liver tissue was independently performed using two RNaseq platforms (SOLiD and Illumina) and also by droplet digital PCR (ddPCR) and quantitative RT-PCR. The proteome profiling of Chr 18 was accomplished by quantitatively measuring protein copy numbers in the three types of biomaterial (the lowest protein concentration measured was 10(-13) M) using selected reaction monitoring (SRM). In total, protein copy numbers were estimated for 228 master proteins, including quantitative data on 164 proteins in plasma, 171 in the HepG2 cell line, and 186 in liver tissue. Most proteins were present in plasma at 10(8) copies/µL, while the median abundance was 10(4) and 10(5) protein copies per cell in HepG2 cells and liver tissue, respectively. In summary, for liver tissue and HepG2 cells a "transcriptoproteome" was produced that reflects the relationship between transcript and protein copy numbers of the genes on Chr 18. The quantitative data acquired by RNaseq, PCR, and SRM were uploaded into the "Update_2013" data set of our knowledgebase (www.kb18.ru) and investigated for linear correlations.

Chromosome 18 transcriptome profiling and targeted proteome mapping in depleted plasma, liver tissue and HepG2 cells.

The final goal of the Russian part of the Chromosome-centric Human Proteome Project (C-HPP) was established as the analysis of the chromosome 18 (Chr 18) protein complement in plasma, liver tissue and HepG2 cells with the sensitivity of 10(-18) M. Using SRM, we have recently targeted 277 Chr 18 proteins in plasma, liver, and HepG2 cells. On the basis of the results of the survey, the SRM assays were drafted for 250 proteins: 41 proteins were found only in the liver tissue, 82 proteins were specifically detected in depleted plasma, and 127 proteins were mapped in both samples. The targeted analysis of HepG2 cells was carried out for 49 proteins; 41 of them were successfully registered using ordinary SRM and 5 additional proteins were registered using a combination of irreversible binding of proteins on CN-Br Sepharose 4B with SRM. Transcriptome profiling of HepG2 cells performed by RNAseq and RT-PCR has shown a significant correlation (r = 0.78) for 42 gene transcripts. A pilot affinity-based interactome analysis was performed for cytochrome b5 using analytical and preparative optical biosensor fishing followed by MS analysis of the fished proteins. All of the data on the proteome complement of the Chr 18 have been integrated into our gene-centric knowledgebase ( www.kb18.ru ).

Application of matrix-assisted laser desorption/ionization time-of-flight mass spectrometry for the study of Helicobacter pylori.

The characteristics of matrix-assisted laser desorption/ionization time-of-flight (MALDI-TOF) mass spectrometry based investigation of extremely variable bacteria such as Helicobacter pylori were studied. H. pylori possesses a very high natural variability. Accurate tools for species identification and epidemiological characterization could help the scientific community to better understand the transmission pathways and virulence mechanisms of these bacteria. Seventeen clinical as well as two laboratory strains of H. pylori were analyzed by the MALDI Biotyper method for rapid species identification. Mass spectra collected were found containing 7-13 significant peaks per sample, and only six protein signals were identical for more than half of the strains. Four of them could be assigned to ribosomal proteins RL32, RL33, RL34, and RL36. The reproducible peak with m/z 6948 was identified as a histidine-rich metal-binding polypeptide by tandem mass spectrometry (MS/MS). In spite of the evident protein heterogeneity of H. pylori the mass spectra collected for a particular strain under several cultivations were highly reproducible. Moreover, all clinical strains were perfectly identified as H. pylori species through comparative analysis using the MALDI Biotyper software (Bruker Daltonics, Germany) by pattern matching against a database containing mass spectra from different microbial strains (n = 3287) including H. pylori 26695 and J99. The results of this study allow the conclusion that the MALDI-TOF direct bacterial profiling is suited for H. pylori identification and could be supported by mass spectra fragmentation of the observed polypeptide if necessary.