Search

Highlights • Central components of subway samples are soil and human skin microbes. • Microbial diversity is correlated with passenger traffic. • Co-occurring groups of microbes resonate with other studies of subway microbiome. • No substantial evidence of major pathogens presence was detected., Microbial ecosystems of the built environments have become key mediators of health as people worldwide tend to spend large amount of time indoors. Underexposure to microbes at an early age is linked to increased risks of allergic and autoimmune diseases. Transportation systems are of particular interest, as they are globally the largest space for interactions between city-dwellers. Here we performed the first pilot study of the Moscow subway microbiome by analyzing swabs collected from 5 types of surfaces at 4 stations using high-throughput 16S rRNA gene sequencing. The study was conducted as a part of The Metagenomics and Metadesign of the Subways and Urban Biomes (MetaSUB) project. The most abundant microbial taxa comprising the subway microbiome originated from soil and human skin. Microbiome diversity was positively correlated with passenger traffic. No substantial evidence of major human pathogens presence was found. Co-occurrence analysis revealed clusters of microbial genera including combinations of microbes likely originating from different niches. The clusters as well as the most abundant microbes were similar to ones obtained for the published data on New-York City subway microbiome. Our results suggest that people are the main source and driving force of diversity in subway-associated microbiome. The data form a basis for a wider survey of Moscow subway microbiome to explore its longitudinal dynamics by analyzing an extended set of sample types and stations. Complementation of methods with viability testing, “shotgun” metagenomics, sequencing of bacterial isolates and culturomics will provide insights for public health, biosafety, microbial ecology and urban design.

Probiotics and vitamins can impact immune responses and modulate gut microbiome. We evaluated the effects of consuming a yogurt fortified with vitamins and probiotic Lacticaseibacillus casei and rhamnosus on upper respiratory tract infections frequency and gut microbiome during a 3-month intervention. The study included 2 case groups (consuming different flavours) and a placebo group of healthy adults (n = 158–160 in each group). The effects on URTI-related parameters in both case groups were gender-specific. The female subjects had lower URTI incidence and frequency, while for the males no significant differences were found. The URTI duration was shorter in one of the probiotic groups for females and in both such groups - for males. The observed changes in microbiome composition, blood and stool parameters were not different from those observed in the placebo group. Consumption of fortified fermented dairy foods is promising for improving immunity status within the general population.

Hyperglycemia may contribute to the progression of carcinomas by triggering epithelial-to-mesenchymal transition (EMT). Some proteostasis systems are involved in metastasis; in this paper, we sought to explore the mechanism of Hsp70 chaperone in EMT. We showed that knockdown of Hsp70 reduced cell migration capacity concomitantly with levels of mRNA of the Slug, Snail, and Twist markers of EMT, in colon cancer cells incubated in high glucose medium. Conversely, treatment of cells with Hsp70 inducer U-133 were found to elevate cell motility, along with the other EMT markers. To prove that inhibiting Hsp70 may reduce EMT efficiency, we treated cells with a CL-43 inhibitor of the HSF1 transcription factor, which lowered Hsp70 and HSF1 content in the control and induced EMT in carcinoma cells. Importantly, CL-43 reduced migration capacity, EMT-linked transcription factors, and increased content of epithelial marker E-cadherin in colon cancer cells of three lines, including one derived from a clinical sample. To prove that Hsp70 chaperone should be targeted when inhibiting the EMT pathway, we treated cancer cells with 2-phenylethynesulfonamide (PES) and demonstrated that the compound inhibited substrate-binding capacity of Hsp70. Furthermore, PES suppressed EMT features, cell motility, and expression of specific transcription factors. In conclusion, the Hsp70 chaperone machine efficiently protects mechanisms of the EMT, and the safe inhibitors of the chaperone are needed to hamper metastasis at its initial stage.

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. Electronic supplementary material The online version of this article (10.1186/s12864-018-5306-5) contains supplementary material, which is available to authorized users.

Background Viliuisk encephalomyelitis (VE) is a rare endemic neurodegenerative disease occurring in the Yakut population of Northeastern Siberia. The main clinical features of VE are spasticity, dysarthria, dementia, central paresis and paralysis, and cortical atrophy observed via MRI. Many hypotheses have been proposed regarding its etiology, including infectious agents, genetics, environmental factors, and immunopathology. Each of these hypotheses has been supported to some extent by epidemiological and experimental data. Nevertheless, none of them has been decisively proven. Gut microbiome is one of the factors that might be involved in VE pathogenesis. Results Here we performed a pilot survey of the stool microbiomes of Yakut subjects with VE (n = 6) and without VE (n = 11). 16S rRNA sequencing showed that in comparison with the control group, the Yakuts with VE had increased proportions of Methanobrevibacter and Christensenella, which are reported to be linked to body mass index, metabolism, dietary habits and potentially to neurodegenerative disorders. The identified associations suggest that the microbiome may be involved in VE. Overall, the Yakut microbiome was quite specific in comparison with other populations, such as metropolitan Russians and native inhabitants of the Canadian Arctic. Conclusions Describing the gut microbiome of indigenous human populations will help to elucidate the impact of dietary and environmental factors on microbial community structure and identify risks linked to the lifestyles of such groups as well as endemic diseases.

Genome-wide association studies have led to a significant progress in identification of genomic loci affecting coronary artery disease (CAD) risk. However, revealing the causal genes responsible for the observed associations is challenging. In the present study, we aimed to prioritize CAD-relevant genes based on cumulative evidence from the published studies and our own study of colocalization between eQTLs and loci associated with CAD using SMR/HEIDI approach. Prior knowledge of candidate genes was extracted from both experimental and in silico studies, employing different prioritization algorithms. Our review systematized information for a total of 51 CAD-associated loci. We pinpointed 37 genes in 36 loci. For 27 genes we infer they are causal for CAD, and for 10 further genes we judge them most likely causal. Colocalization analysis showed that for 18 out of these loci, association with CAD can be explained by changes in gene expression in one or more CAD-relevant tissues. Furthermore, for 8 out of 36 loci, existing evidence suggested additional CAD-associated genes. For the remaining 15 loci, we concluded that evidence for gene prioritization remains inconsistent, insufficient, or absent. Our results provide deeper insights into the genetic etiology of CAD and demonstrate knowledge gaps where further research is warranted.

We aimed to replicate, in a specific athletic event cohort (only track and field) and in two different ethnicities (Japanese and East European, i.e. Russian and Polish), original findings showing the association of the angiotensin-II receptor type-2 gene (AGTR2) rs11091046 A>C polymorphism with athlete status. We compared genotypic frequencies of the AGTR2 rs11091046 polymorphism among 282 track and field sprint/power athletes (200 men and 82 women), including several national record holders and Olympic medallists (214 Japanese, 68 Russian and Polish), and 2024 control subjects (842 men and 1182 women) (804 Japanese, 1220 Russian and Polish). In men, a meta-analysis from the two combined cohorts showed a significantly higher frequency of the C allele in athletes than in controls (odds ratio: 1.62, P=0.008, heterogeneity index I 2 =0%). With regard to respective cohorts, C allele frequency was higher in Japanese male athletes than in controls (67.7% vs. 55.9%, P=0.022), but not in Russian/Polish male athletes (61.9% vs. 51.0%, P=0.172). In women, no significant results were obtained by meta-analysis for the two cohorts combination (P=0.850). The AC genotype frequency was significantly higher in Russian/Polish women athletes than in controls (69.2% vs. 42.1%, P=0.022), but not in Japanese women athletes (P=0.226). Our results, in contrast to previous findings, suggested by meta-analysis that the C allele of the AGTR2 rs11091046 polymorphism is associated with sprint/power track and field athlete status in men, but not in women. Sin financiación 2.202 JCR (2018) Q2, 33/83 Sport Sciences 0.941 SJR (2018) Q1, 64/301 Orthopedics and Sports Medicine No data IDR 2018 UEM

© 2017Alcoholism is associated with significant changes in gut microbiota composition. Metagenomic sequencing allows to assess the altered abundance levels of bacterial taxa and genes in a culture-independent way. We collected 99 stool samples from the patients with alcoholic dependence syndrome (n=72) and alcoholic liver cirrhosis (n=27). Each of the samples was surveyed using “shotgun” (whole-genome) sequencing on SOLiD platform. The reads are deposited in the ENA (project ID: PRJEB18041).

Surveys of environmental microbial communities using metagenomic approach produce vast volumes of multidimensional data regarding the phylogenetic and functional composition of the microbiota. Faced with such complex data, a metagenomic researcher needs to select the means for data analysis properly. Data visualization became an indispensable part of the exploratory data analysis and serves a key to the discoveries. While the molecular-genetic analysis of even a single bacterium presents multiple layers of data to be properly displayed and perceived, the studies of microbiota are significantly more challenging. Here we present a review of the state-of-art methods for the visualization of metagenomic data in a multi-level manner: from the methods applicable to an in-depth analysis of a single metagenome to the techniques appropriate for large-scale studies containing hundreds of environmental samples.

Accumulated data suggests that the gut microbiome can rapidly respond to changes in diet. Consumption of fermented dairy products (FDP) fortified with probiotic microbes may be associated with positive impact on human health. However, the extent and details of the possible impact of FDP consumption on gut community structure tends to vary across individuals. We used microbiome analysis to characterize changes in gut microbiota composition after 30 days of oral intake of a yoghurt fortified with Bifidobacterium animalis subsp. lactis BB-12. 16S rRNA gene sequencing was used to assess the gut microbial composition before and after FDP consumption in healthy adults (n = 150). Paired comparison of gut microbial content demonstrated an increase in presence of potentially beneficial bacteria, particularly, Bifidobacterium genus, as well as Adlercreutzia equolifaciens and Slackia isoflavoniconvertens. At a functional level, an increased capacity to metabolize lactose and synthesize amino acids was observed accompanied by a lowered potential for synthesis of lipopolysaccharides. Cluster analysis revealed that study volunteers segregated into two groups with post-intervention microbiota response that was dependent on the baseline microbial community structure.

Gut microbiota is an up-to date theme that is widely studied. Now it is obvious that microbiota could be considered as one of the factors affecting aging process. Some reports suggested the potential mechanisms by which gut microbes may influence human health span. Aging-related conditions are accompanied by decline in immune system functioning and low-grade inflammation, and our microbes play a certain role in these processes. Gut microbiota composition changes during the aging process may indicate a more pronounced aging-related health loss.

Purpose of the study. To study the relationship between the gut microbiota composition, the state of the vascular wall and nutrient intake in patients without clinical manifestation of cardiovascular disease.Materials and methods. The study included 92 residents from Moscow and the Moscow region, men and women aged 25 to 76 y.o., without clinical manifestations of chronic diseases, with no medical treatment but with possible presence of cardiovascular risk factors. All the participants were examined with duplex scanning of the carotid arteries with the definition of the thickness of the intima-media complex (IMC), the presence of atherosclerotic plaques, measuring carotid-femoral pulse wave velocity (PWV), the measurement of interleukin-6 (IL-6) and C-reactive protein (C-RP) and sequencing of the V3-V4 variable regions of 16S rRNA genes of the gut microbiota.Diet was studied by quantitative assessment of consumer products using a standardized computer program «Analysis of Human Nutrition states» (version 1.2.4 Research Institute of Nutrition 2003-2006).Results. The study revealed that in the investigated cohort IMC thickness was greater in those with a large number of bacterial genera Serratia and Blautia. Bifidobacterium representation was increased and Blautia representation was decreased in those who consumed a lot of starch.Conclusion. The composition of the gut microbiota is associated with both vascular wall stiffness and atherosclerosis. In addition, the bacteria associated with subclinical atherosclerosis were more represented among those who consumed smaller amounts of complex carbohydrates.

Urothelial carcinoma (UC) is the most common type of bladder cancer and is the second most frequently diagnosed genitourinary tumor. The identification of fusion genes in bladder cancer might provide new perspectives for its classification and significance. In this study, we present a thorough search on three UC samples for novel fusion transcripts in bladder cancer using high-throughput RNA sequencing. We used stringent requirements for 819 fusion candidates and nominated 10 candidate fusion transcripts. Among them four novel fusion genes SEPT9/CYHR, IGF1R/TTC23, SYT8/TNNI2 and CASZ1/DFFA were validated and characterized in 48 formalin-fixed paraffin-embedded (FFPE) specimens of bladder cancer. Chromosomal rearrangements of regions 17q25, 15q26.3 and 1p36.22 resulting in the fusion transcripts SEPT9/CYHR, IGF1R/TTC23 and CASZ1/DFFA, appeared to be rare or unique events because they were not detected in the 48 UC samples. In contrast, the SYT8/TNNI2 fusion transcript resulting from transcription-induced chimerism by read-through mechanisms was a rather common and tumor-specific event occurring in 37.5% (18/48) of the UC specimens. Further investigation of functional and clinical relevance of novel fusion genes remains to be elucidated to reveal their role in bladder carcinogenesis.

Nowadays the advent of innovative high-throughput sequencing allows obtaining high-quality microbiome profiling. However, PCR-based tests are still considered the "golden standard" for many clinical applications. Here, we designed a qPCR-based platform with fluorescent-labeled oligonucleotide probes for assessing human gut microbiome composition. The system allows conducting qualitative and semiquantitative analysis for 12 prokaryotic taxa that are prevalent in the human gut and associated with diseases, diet, age and other factors. The platform was validated by comparing microbiome profile data obtained with two different methods - the platform and high-throughput 16S rRNA sequencing - across 42 stool samples. The test can form the basis for precise and cost-efficient microbiome assay for large-scale surveys including clinical trials with interventions related to diet and disease risks.

The aim of this paper was to study gut microbiota composition in patients with different metabolic statuses. Methods: 92 participants aged 25&ndash, 76 years (26 of whom were men), with confirmed absence of cardiovascular and other chronic diseases (but with the possible presence of cardiovascular risk factors) were included. Carotid ultrasound examinations, 16S rRNA sequencing of stool samples and diet assessments were performed. Statistical analysis was performed using R programming language, 3.1.0. Results: Enterotyping yielded two clusters differentiated by alpha-diversity. Intima-media thickness was higher in the cluster with lower diversity (adj. p &lt, 0.001). Obesity was associated with higher Serratia (adj. p = 0.003) and Prevotella (adj. p &lt, 0.0003) in relative abundance. Abdominal obesity was associated with higher abundance of Serratia (adj. p = 0.004) and Prevotella (adj. p = 0.0008) and lower levels of Oscillospira (adj. p = 0.0005). Glucose metabolism disturbances were associated with higher Blautia (adj. p = 0.0007) and Serratia (adj. p = 0.003) prevalence. Arterial hypertension was associated with high Blautia levels (adj. p = 0.002). The Blautia genus strongly correlated with low resistant starch consumption (adj. p = 0.007). A combination of high-fat diet and elevated Blautia levels was very common for diabetes mellitus type 2 patients (adj. p = 0.0001). Conclusion: The results show that there is a relationship between metabolic changes and higher representation of opportunistic pathogens and low diversity of gut microbiota even in apparently healthy participants.

Personalized nutrition is of increasing interest to individuals actively monitoring their health. The relations between the duration of diet intervention and the effects on gut microbiota have yet to be elucidated. Here we examined the associations of short-term dietary changes, long-term dietary habits and lifestyle with gut microbiota. Stool samples from 248 citizen-science volunteers were collected before and after a self-reported 2-week personalized diet intervention, then analyzed using 16S rRNA sequencing. Considerable correlations between long-term dietary habits and gut community structure were detected. A higher intake of vegetables and fruits was associated with increased levels of butyrate-producing Clostridiales and higher community richness. A paired comparison of the metagenomes before and after the 2-week intervention showed that even a brief, uncontrolled intervention produced profound changes in community structure: resulting in decreased levels of Bacteroidaceae, Porphyromonadaceae and Rikenellaceae families and decreased alpha-diversity coupled with an increase of Methanobrevibacter, Bifidobacterium, Clostridium and butyrate-producing Lachnospiraceae- as well as the prevalence of a permatype (a bootstrapping-based variation of enterotype) associated with a higher diversity of diet. The response of microbiota to the intervention was dependent on the initial microbiota state. These findings pave the way for the development of an individualized diet.

SummaryMetagenomic surveys of human microbiota are becoming increasingly widespread in academic research as well as in food and pharmaceutical industries and clinical context. Intuitive tools for exploration of experimental data are of high interest to researchers. Knomics-Biota is a Web-based resource for exploratory analysis of human gut metagenomes. Users can generate analytical reports that correspond to common experimental schemes (like case-control study or paired comparison). Statistical analysis and visualizations of microbiota composition are provided in association with the external factors and in the context of thousands of publicly available datasets.Availability and ImplementationThe Web-service is available at https://biota.knomics.ru.Contactanna.popenko@knomics.ru or a.tyakht@gmail.com.Supplementary informationSupplementary figures are available at Bioinformatics online.

Long Interspersed Nuclear Element-1 (LINE-1, L1) is a mobile genetic element active in human genomes. L1-encoded ORF1 and ORF2 proteins bind L1 RNAs, forming ribonucleoproteins (RNPs). These RNPs interact with diverse host proteins, some repressive and others required for the L1 lifecycle. Using differential affinity purifications, quantitative mass spectrometry, and next generation RNA sequencing, we have characterized the proteins and nucleic acids associated with distinctive, enzymatically active L1 macromolecular complexes. Among them, we describe a cytoplasmic intermediate that we hypothesize to be the canonical ORF1p/ORF2p/L1-RNA-containing RNP, and we describe a nuclear population containing ORF2p, but lacking ORF1p, which likely contains host factors participating in target-primed reverse transcription., eLife digest Our genome consists of about two percent genes, while around 60 to 70 percent are made up of hundreds of thousands of copies of very similar DNA sequences. These repeats have accumulated over time due to specific genetic elements called transposons. Transposons are often referred to as ‘jumping genes’, as they can move within the genome and thereby create mutations that may lead to cancer or other genetic diseases. LINE-1 is the only remaining active transposon in humans, and it expands by copying and pasting itself to new locations. To do so, it is first transcribed into RNA – the molecules that help to make proteins – and then converted back into identical DNA sequences. In a never-ending battle, our cells have been fighting to keep LINE-1 and its ancestors from replicating, and so evolved various defense mechanisms. Yet, LINE-1 has learned to circumvent these barriers, and continues to replicate and cause disease. Our understanding of these defenses and of how LINE-1 evades them is limited. Previous research has shown that the LINE-1 RNA and its two encoded proteins, called ORF1p and ORF2p, interact with a series of other proteins, with which they can form different types of complexes. Now, Taylor, Altukhov, Molloy et al. used human embryonic kidney cells grown in the laboratory with different LINE-1 mutations to identify how they affect the bound proteins and RNAs. The results showed that LINE-1 can form at least two different sets of complexes with other proteins. The complex containing ORF1p and ORF2p and several other proteins was located in the cytoplasm, the fluid that fills the cells. However, the experiments also revealed a new complex in the cell nucleus, which contained ORF2p and proteins involved in DNA replication and repair, but not ORF1p. The results suggest ORF1p delivers RNPs to the nucleus around the time the cell divides. Another group of researchers has looked more closely at what happens during cell division. A next step will be to study how exactly LINE-1 contributes to cancer. In the future, overactive LINE-1 proteins could be targeted to kill cancer cells, to identify cancer early, or to see if the cancer has come back. LINE-1 may also provide clues on how the genome has evolved.

Here we present the first metagenomic study of gut microbiota in patients with alcohol dependence syndrome (ADS) performed in the whole-genome ("shotgun") format. Taxonomic analysis highlighted changes in community "drivers" abundance previously associated with inflammatory processes (including increase in Ruminococcus gnavus and torques, as well as decrease in Faecalibacterium and Akkermansia). Microbiota of alcoholics manifested presence of specific opportunistic pathogens rarely detected in healthy control subjects of the world. Differential analysis of metabolic potential basing on changes in KEGG Orthology groups abundance revealed increase in pathways associated with response to oxidative stress. Analysis of two specific gene groups--alcohol metabolism and virulence factors--also showed increase in comparison with the control groups. We suggest that gut microbiota distinct in alcoholics by both taxonomic and functional composition plays role in modulating the effect of alcohol on host organism.Provedeno pervoe metagenomnoe issledovanie mikrobioty kishechnika u patsientov s sindromom alkogol'noĭ zavisimosti (SAZ) v polnogenomnom (angl. “shotgun”) formate. Taksonomicheskiĭ analiz vyiavil izmeneniia v otnositel'noĭ predstavlennosti preobladaiushchikh bakteriĭ, assotsiirovannye s vospaleniem (v tom chisle povyshenie urovnia vidov Ruminococcus gnavus i R. torques, a takzhe ponizhenie urovnia rodov Faecalibacterium i Akkermansia). V mikrobiote patsientov s SAZ obnaruzheno prisutstvie opportunisticheskikh patogenov, redko detektiruemykh v metagenomakh zdorovykh individov iz raznykh stran mira. Soglasno izmeneniiam otnositel'noĭ predstavlennosti grupp genov KEGG Orthology, sravnitel'nyĭ analiz metabolicheskogo potentsiala vyiavil povyshenie urovnia puteĭ, assotsiirovannykh s otvetom na okislitel'nyĭ stress. U patsientov s SAZ obnaruzheno povyshenie dvukh spetsifichnykh grupp genov, kodiruiushchikh fermenty, uchastvuiushchie v metabolizme alkogolia, a takzhe faktory virulentnosti. Predpolozhitel'no, mikrobiota kishechnika bol'nykh alkogolizmom, demonstriruiushchaia izmeneniia kak taksonomicheskogo, tak i funktsional'nogo sostava, igraet rol' v moduliatsii vozdeĭstviia alkogolia na organizm khoziaina.

Ovarian cancer ascites is a native medium for cancer cells that allows investigation of their secretome in a natural environment. This medium is of interest as a promising source of potential biomarkers, and also as a medium for cell-cell communication. The aim of this study was to elucidate specific features of the malignant ascites metabolome and proteome. In order to omit components of the systemic response to ascites formation, we compared malignant ascites with cirrhosis ascites. Metabolome analysis revealed 41 components that differed significantly between malignant and cirrhosis ascites. Most of the identified cancer-specific metabolites are known to be important signaling molecules. Proteomic analysis identified 2096 and 1855 proteins in the ovarian cancer and cirrhosis ascites, respectively; 424 proteins were specific for the malignant ascites. Functional analysis of the proteome demonstrated that the major differences between cirrhosis and malignant ascites were observed for the cluster of spliceosomal proteins. Additionally, we demonstrate that several splicing RNAs were exclusively detected in malignant ascites, where they probably existed within protein complexes. This result was confirmed in vitro using an ovarian cancer cell line. Identification of spliceosomal proteins and RNAs in an extracellular medium is of particular interest; the finding suggests that they might play a role in the communication between cancer cells. In addition, malignant ascites contains a high number of exosomes that are known to play an important role in signal transduction. Thus our study reveals the specific features of malignant ascites that are associated with its function as a medium of intercellular communication.

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. Electronic supplementary material The online version of this article (10.1186/s40168-017-0359-2) contains supplementary material, which is available to authorized users.

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. Electronic supplementary material The online version of this article (doi:10.1186/s12864-017-3917-x) contains supplementary material, which is available to authorized users.

1.SummaryLong Interspersed Nuclear Element-1 (LINE-1, L1) is a mobile genetic element active in human genomes. L1-encoded ORF1 and ORF2 proteins bind L1 RNAs, forming ribonucleoproteins (RNPs). These RNPs interact with diverse host proteins, some repressive and others required for the L1 lifecycle. Using differential affinity purifications and quantitative mass spectrometry, we have characterized the proteins associated with distinctive L1 macromolecular complexes. Our findings support the presence of multiple L1-derived retrotransposition intermediates in vivo. Among them, we describe a cytoplasmic intermediate that we hypothesize to be the canonical ORF1p/ORF2p/L1-RNA-containing RNP, and we describe a nuclear population containing ORF2p, but lacking ORF1p, which likely contains host factors participating in template-primed reverse transcription.

The aim of this research was to study the association between the gut microbiota composition and arterial wall properties. The study included 92 participants, men and women aged 25–76 years old without clinical manifestation of chronic diseases but with the possible presence of cardiovascular risk factors. Carbohydrate metabolism examination, duplex scanning of the carotid arteries with the measurement of the intima-media thickness (IMT), the carotid-femoral pulse wave velocity (PWV) measurement, and 16S rRNA (V3 V4 regions) sequencing of the gut microbiota were performed in all participants. Higher Serratia abundance was associated with the increased IMT and CRP levels. Blautia representation was associated with IMT. Higher Bacteroides representation was associated with higher pulse wave velocity in non-diabetic subjects. Although this study had some limitations, we have demonstrated that the composition of the gut microbiota was associated both with atherosclerotic and arterial stiffness markers.

The gut microbiota is essentially a multifunctional bioreactor within a human being. The exploration of its enormous metabolic potential provides insights into the mechanisms underlying microbial ecology and interactions with the host. The data obtained using "shotgun" metagenomics capture information about the whole spectrum of microbial functions. However, each new study presenting new sequencing data tends to extract only a little of the information concerning the metabolic potential and often omits specific functions. A meta-analysis of the available data with an emphasis on biomedically relevant gene groups can unveil new global trends in the gut microbiota. As a step toward the reuse of metagenomic data, we developed a method for the quantitative profiling of user-defined groups of genes in human gut metagenomes. This method is based on the quick analysis of a gene coverage matrix obtained by pre-mapping the metagenomic reads to a global gut microbial catalogue. The method was applied to profile the abundance of several gene groups related to antibiotic resistance, phages, biosynthesis clusters and carbohydrate degradation in 784 metagenomes from healthy populations worldwide and patients with inflammatory bowel diseases and obesity. We discovered country-wise functional specifics in gut resistome and virome compositions. The most distinct features of the disease microbiota were found for Crohn's disease, followed by ulcerative colitis and obesity. Profiling of the genes belonging to crAssphage showed that its abundance varied across the world populations and was not associated with clinical status. We demonstrated temporal resilience of crAssphage and the influence of the sample preparation protocol on its detected abundance. Our approach offers a convenient method to add value to accumulated "shotgun" metagenomic data by helping researchers state and assess novel biological hypotheses.

Metagenomics, the application of high-throughput DNA sequencing for surveys of environmental samples, has revolutionized our view on the taxonomic and genetic composition of complex microbial communities. An enormous richness of microbiota keeps unfolding in the context of various fields ranging from biomedicine and food industry to geology. Primary analysis of metagenomic reads allows to infer semi-quantitative data describing the community structure. However, such compositional data possess statistical specific properties that are important to be considered during preprocessing, hypothesis testing and interpreting the results of statistical tests. Failure to account for these specifics may lead to essentially wrong conclusions as a result of the survey. Here we present a researcher introduced to the field of metagenomics with the basic properties of microbial compositional data including statistical power and proposed distribution models, perform a review of the publicly available software tools developed specifically for such data and outline the recommendations for the application of the methods.

To investigate the association between multiple single-nucleotide polymorphisms (SNPs), aerobic performance and elite endurance athlete status in Russians. By using GWAS approach, we examined the association between 1,140,419 SNPs and relative maximal oxygen consumption rate ([Formula: see text]O2max) in 80 international-level Russian endurance athletes (46 males and 34 females). To validate obtained results, we further performed case-control studies by comparing the frequencies of the most significant SNPs (with P < 10(-5)-10(-8)) between 218 endurance athletes and opposite cohorts (192 Russian controls, 1367 European controls, and 230 Russian power athletes). Initially, six 'endurance alleles' were identified showing discrete associations with [Formula: see text]O2max both in males and females. Next, case-control studies resulted in remaining three SNPs (NFIA-AS2 rs1572312, TSHR rs7144481, RBFOX1 rs7191721) associated with endurance athlete status. The C allele of the most significant SNP, rs1572312, was associated with high values of [Formula: see text]O2max (males: P = 0.0051; females: P = 0.0005). Furthermore, the frequency of the rs1572312 C allele was significantly higher in elite endurance athletes (95.5%) in comparison with non-elite endurance athletes (89.8%, P = 0.0257), Russian (88.8%, P = 0.007) and European (90.6%, P = 0.0197) controls and power athletes (86.2%, P = 0.0005). The rs1572312 SNP is located on the nuclear factor I A antisense RNA 2 (NFIA-AS2) gene which is supposed to regulate the expression of the NFIA gene (encodes transcription factor involved in activation of erythropoiesis and repression of the granulopoiesis). Our data show that the NFIA-AS2 rs1572312, TSHR rs7144481 and RBFOX1 rs7191721 polymorphisms are associated with aerobic performance and elite endurance athlete status.

The avian bacterial pathogen Mycoplasma gallisepticum is a good model for systems studies due to small genome and simplicity of regulatory pathways. In this study, we used RNA-Seq and MS-based proteomics to accurately map coding sequences, transcription start sites (TSSs) and transcript 3′-ends (T3Es). We used obtained data to investigate roles of TSSs and T3Es in stress-induced transcriptional responses. We identified 1061 TSSs at a false discovery rate of 10% and showed that almost all transcription in M. gallisepticum is initiated from classic TATAAT promoters surrounded by A/T-rich sequences. Our analysis revealed the pronounced operon structure complexity: on average, each coding operon has one internal TSS and T3Es in addition to the primary ones. Our transcriptomic approach based on the intervals between the two nearest transcript ends allowed us to identify two classes of T3Es: strong, unregulated, hairpin-containing T3Es and weak, heat shock-regulated, hairpinless T3Es. Comparing gene expression levels under different conditions revealed widespread and divergent transcription regulation in M. gallisepticum. Modeling suggested that the core promoter structure plays an important role in gene expression regulation. We have shown that the heat stress activation of cryptic promoters combined with the hairpinless T3Es suppression leads to widespread, seemingly non-functional transcription.

Muscle fibre type is a heritable trait and can partly predict athletic success. It has been proposed that polymorphisms of genes involved in the regulation of muscle fibre characteristics may predispose the muscle precursor cells of a given individual to be predominantly fast or slow. In the present study, we examined the association between 15 candidate gene polymorphisms and muscle fibre type composition of the vastus lateralis muscle in 55 physically active, healthy men. We found that rs11091046 C allele carriers of the angiotensin II type 2 receptor gene (AGTR2; involved in skeletal muscle development, metabolism and circulatory homeostasis) had a significantly higher percentage of slow-twitch fibres than A allele carriers [54.2 (11.1) versus 45.2 (10.2)%; P = 0.003]. These data indicate that 15.2% of the variation in muscle fibre composition of the vastus lateralis muscle can be explained by the AGTR2 genotype. Next, we investigated the frequencies of the AGTR2 alleles in 2178 Caucasian athletes and 1220 control subjects. The frequency of the AGTR2 C allele was significantly higher in male and female endurance athletes compared with power athletes (males, 62.7 versus 51.7%, P = 0.0038; females, 56.6 versus 48.1%, P = 0.0169) and control subjects (males, 62.7 versus 51.0%, P = 0.0006; elite female athletes, 65.1 versus 55.2%, P = 0.0488). Furthermore, the frequency of the AGTR2 A allele was significantly over-represented in female power athletes (51.9%) in comparison to control subjects (44.8%, P = 0.0069). We also found that relative maximal oxygen consumption was significantly greater in male endurance athletes with the AGTR2 C allele compared with AGTR2 A allele carriers [n = 28; 62.3 (4.4) versus 57.4 (6.0) ml min(-1) kg(-1); P = 0.0197]. Taken together, these results demonstrate that the AGTR2 gene C allele is associated with an increased proportion of slow-twitch muscle fibres, endurance athlete status and aerobic performance, while the A allele is associated with a higher percentage of fast-twitch fibres and power-oriented disciplines.

Exercise-induced oxidative stress is a state that primarily occurs in athletes involved in high-intensity sports when pro-oxidants overwhelm the antioxidant defense system to oxidize proteins, lipids, and nucleic acids. During exercise, oxidative stress is linked to muscle metabolism and muscle damage, because exercise increases free radical production. The T allele of the Ala16Val (rs4880 C/T) polymorphism in the mitochondrial superoxide dismutase 2 (SOD2) gene has been reported to reduce SOD2 efficiency against oxidative stress. In the present study we tested the hypothesis that the SOD2 TT genotype would be underrepresented in elite athletes involved in high-intensity sports and associated with increased values of muscle and liver damage biomarkers. The study involved 2664 Caucasian (2262 Russian and 402 Polish) athletes. SOD2 genotype and allele frequencies were compared to 917 controls. Muscle and liver damage markers [creatine kinase (CK), creatinine, alanine transaminase (ALT), aspartate transaminase (AST), alkaline phosphatase (ALP)] were examined in serum from 1444 Russian athletes. The frequency of the SOD2 TT genotype (18.6%) was significantly lower in power/strength athletes (n = 524) compared to controls (25.0%, p = 0.0076) or athletes involved in low-intensity sports (n = 180; 33.9%, p < 0.0001). Furthermore, the SOD2 T allele was significantly associated with increased activity of CK (females: p = 0.0144) and creatinine level (females: p = 0.0276; males: p = 0.0135) in athletes. Our data show that the SOD2 TT genotype might be unfavorable for high-intensity athletic events. © 2014 Informa UK, Ltd.

A multitude of metagenomic studies has brought to light an enormous richness of human gut microbiota compositions. In this space of possible configurations, clinical specialists are trying to mine the markers of healthy microbiota via case-control and longitudinal studies. We have discovered potentially beneficial communities while examining the microbial diversity in rural Russians in comparison with the urban dwellers. In this addendum, we further examine the data by elaborating on some of the less common types and suggesting the possible co-metabolism of their drivers. In the light of the first validated clinically effective bacterial transplantation, we discuss the concept of a reference healthy microbiota, outline the problems encountered on the way to its restoration in the developed world, and speculate if rural communities can serve as a source for its prototype.

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.

Fragilysin (BFT) is metalloprotease that is secreted by enterotoxigenic Bacteroides fragilis. Studying the mechanism of BFT interaction with intestinal epithelial cells requires a pure protein sample. In this study, we cloned DNA-fragments coding for the catalytic domain of fragilysin-2 and profragilysin-2 into an E. coli expression vector. Purification methods for the recombinant fragilysin-2 catalytic domain and profragilysin-2 were developed. In addition, we obtained mature active fragilysin-2 from recombinant proprotein by limited tryptic digestion. We tested the biological activity of the recombinant protein samples and revealed that E-cadherin was cleaved when HT-29 cells were treated with mature fragilysin-2 but not with profragilysin-2. Azocoll, azocasein and gelatin were not proteolytically cleaved by mature fragilysin-2. Proteins released in culture medium after HT-29 cells treatment with mature active BFT-2 were identified.

SummaryWe created ResistoMap — a Web-based interactive visualization of the presence of genetic determinants conferring resistance to antibiotics, biocides and heavy metals in human gut microbiota. ResistoMap displays the data about more than 1600 published gut metagenomes of the world populations including both healthy subjects and patients. Multiparameter display filters allow visual assessment of the associations between the meta-data and proportions of resistome. The geographic map navigation layer allows to state hypotheses regarding the global trends of antibiotic resistance and correlate the gut resistome variations with the national clinical guidelines on antibiotics application.Availability and implementationResistoMap is publicly available athttp://resistomap.datalaboratory.ru

Beijing B0/W148, a “successful” clone of Mycobacterium tuberculosis, is widespread in the Russian Federation and some countries of the former Soviet Union. Here, we used label-free gel-LC-MS/MS shotgun proteomics to discover features of Beijing B0/W148 strains that could explain their success. Qualitative and quantitative proteome analyses of Beijing B0/W148 strains allowed us to identify 1,868 proteins, including 266 that were differentially abundant compared with the control strain H37Rv. To predict the biological effects of the observed differences in protein abundances, we performed Gene Ontology analysis together with analysis of protein-DNA interactions using a gene regulatory network. Our results demonstrate that Beijing B0/W148 strains have increased levels of enzymes responsible for long-chain fatty acid biosynthesis, along with a coincident decrease in the abundance of proteins responsible for their degradation. Together with high levels of HsaA (Rv3570c) protein, involved in steroid degradation, these findings provide a possible explanation for the increased transmissibility of Beijing B0/W148 strains and their survival in host macrophages. Among other, we confirmed a very low level of the SseA (Rv3283) protein in Beijing B0/W148 characteristic for all «modern» Beijing strains, which could lead to increased DNA oxidative damage, accumulation of mutations and potentially facilitate the development of drug resistance.

Background A rapidly increasing flow of genomic data requires the development of efficient methods for obtaining its compact representation. Feature extraction facilitates classification, clustering and model analysis for testing and refining biological hypotheses. “Shotgun” metagenome is an analytically challenging type of genomic data - containing sequences of all genes from the totality of a complex microbial community. Recently, researchers started to analyze metagenomes using reference-free methods based on the analysis of oligonucleotides (k-mers) frequency spectrum previously applied to isolated genomes. However, little is known about their correlation with the existing approaches for metagenomic feature extraction, as well as the limits of applicability. Here we evaluated a metagenomic pairwise dissimilarity measure based on short k-mer spectrum using the example of human gut microbiota, a biomedically significant object of study. Results We developed a method for calculating pairwise dissimilarity (beta-diversity) of “shotgun” metagenomes based on short k-mer spectra (5≤k≤11). The method was validated on simulated metagenomes and further applied to a large collection of human gut metagenomes from the populations of the world (n=281). The k-mer spectrum-based measure was found to behave similarly to one based on mapping to a reference gene catalog, but different from one using a genome catalog. This difference turned out to be associated with a significant presence of viral reads in a number of metagenomes. Simulations showed limited impact of bacterial genetic variability as well as sequencing errors on k-mer spectra. Specific differences between the datasets from individual populations were identified. Conclusions Our approach allows rapid estimation of pairwise dissimilarity between metagenomes. Though we applied this technique to gut microbiota, it should be useful for arbitrary metagenomes, even metagenomes with novel microbiota. Dissimilarity measure based on k-mer spectrum provides a wider perspective in comparison with the ones based on the alignment against reference sequence sets. It helps not to miss possible outstanding features of metagenomic composition, particularly related to the presence of an unknown bacteria, virus or eukaryote, as well as to technical artifacts (sample contamination, reads of non-biological origin, etc.) at the early stages of bioinformatic analysis. Our method is complementary to reference-based approaches and can be easily integrated into metagenomic analysis pipelines. Electronic supplementary material The online version of this article (doi:10.1186/s12859-015-0875-7) contains supplementary material, which is available to authorized users.

There are strong genetic components to cardiorespiratory fitness and its response to exercise training. It would be useful to understand the differences in the genomic profile of highly trained endurance athletes of world class caliber and sedentary controls. An international consortium (GAMES) was established in order to compare elite endurance athletes and ethnicity-matched controls in a case-control study design. Genome-wide association studies were undertaken on two cohorts of elite endurance athletes and controls (GENATHLETE and Japanese endurance runners), from which a panel of 45 promising markers was identified. These markers were tested for replication in seven additional cohorts of endurance athletes and controls: from Australia, Ethiopia, Japan, Kenya, Poland, Russia and Spain. The study is based on a total of 1520 endurance athletes (835 who took part in endurance events in World Championships and/or Olympic Games) and 2760 controls. We hypothesized that world-class athletes are likely to be characterized by an even higher concentration of endurance performance alleles and we performed separate analyses on this subsample. The meta-analysis of all available studies revealed one statistically significant marker (rs558129 at GALNTL6 locus, p = 0.0002), even after correcting for multiple testing. As shown by the low heterogeneity index (I2 = 0), all eight cohorts showed the same direction of association with rs558129, even though p-values varied across the individual studies. In summary, this study did not identify a panel of genomic variants common to these elite endurance athlete groups. Since GAMES was underpowered to identify alleles with small effect sizes, some of the suggestive leads identified should be explored in expanded comparisons of world-class endurance athletes and sedentary controls and in tightly controlled exercise training studies. Such studies have the potential to illuminate the biology not only of world class endurance performance but also of compromised cardiac functions and cardiometabolic diseases. Fondo de Investigaciones Sanitarias (grant # PI12/00914). Research in sports genetics to A. Lucia in Spain was funded by grants from the Consejo Superior de Deportes (CSD) 2.806 JCR (2016) Q1, 15/64 Multidisciplinary Sciences No data SJR 2016 No data IDR 2016 UEM

We created ResistoMap—a Web-based interactive visualization of the presence of genetic determinants conferring resistance to antibiotics, biocides and heavy metals in human gut microbiota. ResistoMap displays the data on more than 1500 published gut metagenomes of world populations including both healthy subjects and patients. Multiparameter display filters allow visual assessment of the associations between the meta-data and proportions of resistome. The geographic map navigation layer allows to state hypotheses regarding the global trends of antibiotic resistance and correlates the gut resistome variations with the national clinical guidelines on antibiotics application. Availability and Implementation ResistoMap was implemented using AngularJS, CoffeeScript, D3.js and TopoJSON. The tool is publicly available at http://resistomap.rcpcm.org. Supplementary information Supplementary data are available at Bioinformatics online.

This report describes a new method for desorption of low-molecular weight (LMW) peptides from abundant blood proteins for use in subsequent mass spectrometry analyses. Heating of diluted blood serum to 98 °C for 15 min resulted in dissociation of LMW peptides from the most abundant blood proteins. Application of blood plasma/serum fractionation using magnetic beads with a functionalized surface followed by heating of the resultant fractions significantly increases the number of LMW peptides detected by MALDI–TOF MS, enhances the general reproducibility of mass spectrometry profiles and considerably increases the number of identified blood serum peptides by LC–MS/MS using an Agilent 6520 Accurate-Mass Q-TOF.

Using ClinProt magnetic beads with reverse-phase (MB-HIC 8 and HB-HIC 18), weak cation exchange (MB-WCX) and metal affinity (MB-IMAC Cu) surfaces fractions of peptides and proteins were isolated from human sera for their profiling by MALDI-TOF mass spectrometry. Proteome profiling of sera from basically healthy women (47 subjects, average age 49) and from women with verified ovarian cancer (stages 1-IV, 47 patients, average age 51) by means of MB-WCX beads allowed to generate the best diagnostic models based on Genetic Algorithm and Supervised Neural Network classifiers; these models demonstrated 100% sensitivity and specificity during analysis of the test set. Introduction of additional sera from patients with colorectal cancer (19) and ulcerous colitis (5) to the statistical model confirmed 100% ovarian cancer recognition. Statistical analysis of mass-spectrometry peak areas included to the diagnostic classifiers showed 3 peaks characteristic for ovarian cancer and 4 peak areas exhibiting changes associated with both ovarian and colorectal cancer.

Motivation: High-throughput metagenomic sequencing has revolutionized our view on the structure and metabolic potential of microbial communities. However, analysis of metagenomic composition is often complicated by the high complexity of the community and the lack of related reference genomic sequences. As a start point for comparative metagenomic analysis, the researchers require efficient means for assessing pairwise similarity of the metagenomes (beta-diversity). A number of approaches were used to address this task, however, most of them have inherent disadvantages that limit their scope of applicability. For instance, the reference-based methods poorly perform on metagenomes from previously unstudied niches, while composition-based methods appear to be too abstract for straightforward interpretation and do not allow to identify the differentially abundant features. Results: We developed MetaFast, an approach that allows to represent a shotgun metagenome from an arbitrary environment as a modified de Bruijn graph consisting of simplified components. For multiple metagenomes, the resulting representation is used to obtain a pairwise similarity matrix. The dimensional structure of the metagenomic components preserved in our algorithm reflects the inherent subspecies-level diversity of microbiota. The method is computationally efficient and especially promising for an analysis of metagenomes from novel environmental niches. Availability and Implementation: Source code and binaries are freely available for download at https://github.com/ctlab/metafast. The code is written in Java and is platform independent (tested on Linux and Windows x86_64). Contact: ulyantsev@rain.ifmo.ru Supplementary information: Supplementary data are available at Bioinformatics online.

Type 2 diabetes (T2D) is a serious disease. The gut microbiota (GM) has recently been identified as a new potential risk factor in addition to well-known diabetes risk factors. To investigate the GM composition in association with the dietary patterns in patients with different glucose tolerance, we analyzed 92 patients: with normal glucose tolerance (n=48), prediabetes (preD, n=24), and T2D (n=20). Metagenomic analysis was performed using 16S rRNA sequencing. The diet has been studied by a frequency method with a quantitative evaluation of food intake using a computer program. Microbiota in the samples was predominantly represented by Firmicutes, in a less degree by Bacteroidetes. Blautia was a dominant genus in all samples. The representation of Blautia, Serratia was lower in preD than in T2D patients, and even lower in those with normal glucose tolerance. After the clustering of the samples into groups according to the percentage of protein, fat, carbohydrates in the diet, the representation of the Bacteroides turned to be lower and Prevotella abundance turned to be higher in carbohydrate cluster. There were more patients with insulin resistance, T2D in the fat–protein cluster. Using the Calinski–Harabasz index identified the samples with more similar diets. It was discovered that half of the patients with a high-fat diet had normal tolerance, the others had T2D. The regression analysis showed that these T2D patients also had a higher representation of Blautia. Our study provides the further evidence concerning the structural modulation of the GM in the T2DM pathogenesis depending on the dietary patterns.

Background One of the most challenging tasks in genomic analysis nowadays is metagenomics. Biomedical applications of metagenomics give rise to datasets containing hundreds and thousands of samples from various body sites for hundreds of patients. Inherently metagenome is by far more complex than a single genome as it varies in time by the amount of bacteria comprising it. Other levels of data complexity include geography of the samples and phylogenetic distance between the genomes of the same operational taxonomic unit (OTU). We have developed the visualization concept for the representation of multilayer metagenomics data – the bacterial rose garden. The approach allows to display the taxonomic distance between the representatives of the same OTU in different samples and use variety of the metadata for display. Results We have developed the principle of visualization allowing for multilayer information representation. We have incorporated data on OTU diversity across metagenomes and origin of the samples. The visual representation we have called “rose” is focused on the phylogenetic distance between the representatives of the same OTU. The visual representation is realized as interactive data chart which allows user to interact with data and explore variables. It is known that classical representation of the taxonomic tree is a reduction of information from original pairwise distance matrix. The visualization presented is a way to save all the information available through projection of distance matrix into single dimensional space of one sample. It could serve as a basis for further more complex information representation. We have used the principle proposed for visualization of 101 bacterial OTUs phylogenetic distances, finally we provide open code for the web page generation. Conclusions Bacterial rose garden is a versatile visualization principle coping with the major difficulties of metagenomic big-data visualization without loss of data. The method proposed is showing the interconnectedness of variables and is realized as user-friendly web page allowing for dynamic data exploration. The concept provided serves as one of the original approaches for metagenomic data representation and sharing. Full functional prototype could be found at http://rosegarden.datalaboratory.ru Electronic supplementary material The online version of this article (doi:10.1186/s13040-015-0045-5) contains supplementary material, which is available to authorized users.