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

Background

Crohn’s disease is associated with gut dysbiosis. Independent studies have shown an increase in the abundance of certain bacterial species, particularly Escherichia coli with the adherent-invasive pathotype, in the gut. The role of these species in this disease needs to be elucidated.

Methods

We performed a metagenomic study investigating the gut microbiota of patients with Crohn’s disease. A metagenomic reconstruction of the consensus genome content of the species was used to assess the genetic variability.

Results

The abnormal shifts in the microbial community structures in Crohn’s disease were heterogeneous among the patients. The metagenomic data suggested the existence of multiple E. coli strains within individual patients. We discovered that the genetic diversity of the species was high and that only a few samples manifested similarity to the adherent-invasive varieties. The other species demonstrated genetic diversity comparable to that observed in the healthy subjects. Our results were supported by a comparison of the sequenced genomes of isolates from the same microbiota samples and a meta-analysis of published gut metagenomes.

Conclusions

The genomic diversity of Crohn’s disease-associated E. coli within and among the patients paves the way towards an understanding of the microbial mechanisms underlying the onset and progression of the Crohn’s disease and the development of new strategies for the prevention and treatment of this disease.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5306-5) contains supplementary material, which is available to authorized users.

[Overexpression of microRNAs miR-9, -98, and -199 Correlates with the Downregulation of HK2 Expression in Colorectal Cancer]

Glycolysis activation is one of the main features of energy metabolism in cancer cells that is associated with the increase in glycolytic enzyme synthesis, primarily, hexokinases (HKs), in many types of tumors. Conversely, in colorectal cancer (CRC) the decrease in the expression of HK2 gene, which encodes one of the key rate-limiting enzyme of glycolysis, was revealed, thus, the study of the mechanisms of its inhibition in CRC is of particular interest. To search for potential microRNAs, inhibiting the expression of HK2 in CRC, we have performed the analysis of data from "The Cancer Genome Atlas" (TCGA) and five microRNA-mRNA target interaction databases (TargetScan, DIANA microT, mirSVR (miRanda), PicTar, and miRTarBase) using original CrossHub software. Seven microRNAs containing binding site on mRNA HK2, which expression is negatively correlated with HK2 expression, were selected for further analysis. The expression levels of these microRNAs and mRNA HK2 were estimated by quantitative PCR on a set of CRC samples. It has been shown, that the expression of three microRNAs (miR-9-5p, -98-5p, and -199-5p) was increased and correlated negatively with mRNA level of HK2 gene. Thus, downregulation of HK2 gene may be caused by its negative regulation through microRNAs miR-9-5p, -98-5p, and -199-5p.

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.

Effect of lentivirus-mediated shRNA inactivation of HK1, HK2, and HK3 genes in colorectal cancer and melanoma cells

Background

The switch from oxidative phosphorylation to glycolysis in proliferating cancer cells, even under aerobic conditions, has been shown first in 1926 by Otto Warburg. Today this phenomenon is known as the “Warburg effect” and recognized as a hallmark of cancer. The metabolic shift to glycolysis is associated with the alterations in signaling pathways involved in energy metabolism, including glucose uptake and fermentation, and regulation of mitochondrial functions. Hexokinases (HKs), which catalyze the first step of glycolysis, have been identified to play a role in tumorigenesis of human colorectal cancer (CRC) and melanoma. However, the mechanism of action of HKs in the promotion of tumor growth remains unclear.

Results

The purpose of the present study was to investigate the effect of silencing of hexokinase genes (HK1, HK2, and HK3) in colorectal cancer (HT-29, SW 480, HCT-15, RKO, and HCT 116) and melanoma (MDA-MB-435S and SK-MEL-28) cell lines using short hairpin RNA (shRNA) lentiviral vectors. shRNA lentiviral plasmid vectors pLSLP-HK1, pLSLP-HK2, and pLSLP-HK3 were constructed and then transfected separately or co-transfected into the cells. HK2 inactivation was associated with increased expression of HK1 in colorectal cancer cell lines pointing to the compensation effect. Simultaneous attenuation of HK1 and HK2 levels led to decreased cell viability. Co-transfection with shRNA vectors against HK1, HK2, and HK3 mRNAs resulted in a rapid cell death via apoptosis.

Conclusions

We have demonstrated that simultaneous inactivation of HK1 and HK2 was sufficient to decrease proliferation and viability of melanoma and colorectal cancer cells. Our results suggest that HK1 and HK2 could be the key therapeutic targets for reducing aerobic glycolysis in examined cancers.

The mystery of the fourth clone: comparative genomic analysis of four non-typeable Streptococcus pneumoniae strains with different susceptibilities to optochin

Optochin-resistant pneumococci can be rarely caught in clinical microbiology laboratories because of the routine identification of all such strains as viridans group non-pneumococci. We were lucky to find four non-typeable Streptococcus pneumoniae clones demonstrating the different susceptibilities to optochin: one of them (Spn_13856) was resistant to optochin, while the other three (Spn_1719, Spn_27, and Spn_2298) were susceptible. Whole genome nucleotide sequences of these strains were compared to reveal the differences between the optochin-resistant and optochin-susceptible strains. Two adjacent genes coding maltose O-acetyltransferase and uridine phosphorylase which were presented in the genomes of all optochin-susceptible strains and missed in the optochin-resistant strain were revealed. Non-synonymous substitutions in 14 protein-coding genes were discovered, including the Ala49Ser mutation in the C-subunit of the F0 part of the ATP synthase rotor usually associated with pneumococcal optochin resistance. Modeling of a process of optochin interaction with the F0 part of the ATP synthase rotor indicates that the complex of optochin with "domain C" composed by wild-type C-subunits is more stable than the same complex composed of Ala49Ser mutant C-subunits.

Genome-wide Mycobacterium tuberculosis variation (GMTV) database: a new tool for integrating sequence variations and epidemiology

Background

Tuberculosis (TB) poses a worldwide threat due to advancing multidrug-resistant strains and deadly co-infections with Human immunodeficiency virus. Today large amounts of Mycobacterium tuberculosis whole genome sequencing data are being assessed broadly and yet there exists no comprehensive online resource that connects M. tuberculosis genome variants with geographic origin, with drug resistance or with clinical outcome.

Description

Here we describe a broadly inclusive unifying Genome-wide Mycobacterium tuberculosis Variation (GMTV) database, (http://mtb.dobzhanskycenter.org) that catalogues genome variations of M. tuberculosis strains collected across Russia. GMTV contains a broad spectrum of data derived from different sources and related to M. tuberculosis molecular biology, epidemiology, TB clinical outcome, year and place of isolation, drug resistance profiles and displays the variants across the genome using a dedicated genome browser. GMTV database, which includes 1084 genomes and over 69,000 SNP or Indel variants, can be queried about M. tuberculosis genome variation and putative associations with drug resistance, geographical origin, and clinical stages and outcomes.

Conclusions

Implementation of GMTV tracks the pattern of changes of M. tuberculosis strains in different geographical areas, facilitates disease gene discoveries associated with drug resistance or different clinical sequelae, and automates comparative genomic analyses among M. tuberculosis strains.

Human gut microbiota community structures in urban and rural populations in Russia

The microbial community of the human gut has a crucial role in sustaining host homeostasis. High-throughput DNA sequencing has delineated the structural and functional configurations of gut metagenomes in world populations. The microbiota of the Russian population is of particular interest to researchers, because Russia encompasses a uniquely wide range of environmental conditions and ethnogeographical cohorts. Here we conduct a shotgun metagenomic analysis of gut microbiota samples from 96 healthy Russian adult subjects, which reveals novel microbial community structures. The communities from several rural regions display similarities within each region and are dominated by the bacterial taxa associated with the healthy gut. Functional analysis shows that the metabolic pathways exhibiting differential abundance in the novel types are primarily associated with the trade-off between the Bacteroidetes and Firmicutes phyla. The specific signatures of the Russian gut microbiota are likely linked to the host diet, cultural habits and socioeconomic status.

Unusual large-scale chromosomal rearrangements in Mycobacterium tuberculosis Beijing B0/W148 cluster isolates

The Mycobacterium tuberculosis (MTB) Beijing family isolates are geographically widespread, and there are examples of Beijing isolates that are hypervirulent and associated with drug resistance. One-fourth of Beijing genotype isolates found in Russia belong to the B0/W148 group. The aim of the present study was to investigate features of these endemic strains on a genomic level. Four Russian clinical isolates of this group were sequenced, and the data obtained was compared with published sequences of various MTB strain genomes, including genome of strain W-148 of the same B0/W148 group. The comparison of the W-148 and H37Rv genomes revealed two independent inversions of large segments of the chromosome. The same inversions were found in one of the studied strains after deep sequencing using both the fragment and mate-paired libraries. Additionally, inversions were confirmed by RFLP hybridization analysis. The discovered rearrangements were verified by PCR in all four newly sequenced strains in the study and in four additional strains of the same Beijing B0/W148 group. The other 32 MTB strains from different phylogenetic lineages were tested and revealed no inversions. We suggest that the initial largest inversion changed the orientation of the three megabase (Mb) segment of the chromosome, and the second one occurred in the previously inverted region and partly restored the orientation of the 2.1 Mb inner segment of the region. This is another remarkable example of genomic rearrangements in the MTB in addition to the recently published of large-scale duplications. The described cases suggest that large-scale genomic rearrangements in the currently circulating MTB isolates may occur more frequently than previously considered, and we hope that further studies will help to determine the exact mechanism of such events.

MALINA: a web service for visual analytics of human gut microbiota whole-genome metagenomic reads

MALINA is a web service for bioinformatic analysis of whole-genome metagenomic data obtained from human gut microbiota sequencing. As input data, it accepts metagenomic reads of various sequencing technologies, including long reads (such as Sanger and 454 sequencing) and next-generation (including SOLiD and Illumina). It is the first metagenomic web service that is capable of processing SOLiD color-space reads, to authors’ knowledge. The web service allows phylogenetic and functional profiling of metagenomic samples using coverage depth resulting from the alignment of the reads to the catalogue of reference sequences which are built into the pipeline and contain prevalent microbial genomes and genes of human gut microbiota. The obtained metagenomic composition vectors are processed by the statistical analysis and visualization module containing methods for clustering, dimension reduction and group comparison. Additionally, the MALINA database includes vectors of bacterial and functional composition for human gut microbiota samples from a large number of existing studies allowing their comparative analysis together with user samples, namely datasets from Russian Metagenome project, MetaHIT and Human Microbiome Project (downloaded from http://hmpdacc.org). MALINA is made freely available on the web at http://malina.metagenome.ru. The website is implemented in JavaScript (using Ext JS), Microsoft .NET Framework, MS SQL, Python, with all major browsers supported.