MINTIA: a metagenomic INserT integrated assembly and annotation tool

The earth harbors trillions of bacterial species adapted to very diverse ecosystems thanks to specific metabolic function acquisition. Most of the genes responsible for these functions belong to uncultured bacteria and are still to be discovered. Functional metagenomics based on activity screening is a classical way to retrieve these genes from microbiomes. This approach is based on the insertion of large metagenomic DNA fragments into a vector and transformation of a host to express heterologous genes. Metagenomic libraries are then screened for activities of interest, and the metagenomic DNA inserts of active clones are extracted to be sequenced and analysed to identify genes that are responsible for the detected activity. Hundreds of metagenomics sequences found using this strategy have already been published in public databases. Here we present the MINTIA software package enabling biologists to easily generate and analyze large metagenomic sequence sets, retrieved after activity-based screening. It filters reads, performs assembly, removes cloning vector, annotates open reading frames and generates user friendly reports as well as files ready for submission to international sequence repositories. The software package can be downloaded from https://github.com/Bios4Biol/MINTIA.

Investigating host-microbiome interactions by droplet based microfluidics

BACKGROUND

Despite the importance of the mucosal interface between microbiota and the host in gut homeostasis, little is known about the mechanisms of bacterial gut colonization, involving foraging for glycans produced by epithelial cells. The slow pace of progress toward understanding the underlying molecular mechanisms is largely due to the lack of efficient discovery tools, especially those targeting the uncultured fraction of the microbiota.

RESULTS

Here, we introduce an ultra-high-throughput metagenomic approach based on droplet microfluidics, to screen fosmid libraries. Thousands of bacterial genomes can be covered in 1 h of work, with less than ten micrograms of substrate. Applied to the screening of the mucosal microbiota for β-N-acetylgalactosaminidase activity, this approach allowed the identification of pathways involved in the degradation of human gangliosides and milk oligosaccharides, the structural homologs of intestinal mucin glycans. These pathways, whose prevalence is associated with inflammatory bowel diseases, could be the result of horizontal gene transfers with Bacteroides species. Such pathways represent novel targets to study the microbiota-host interactions in the context of inflammatory bowel diseases, in which the integrity of the mucosal barrier is impaired.

CONCLUSION

By compartmentalizing experiments inside microfluidic droplets, this method speeds up and miniaturizes by several orders of magnitude the screening process compared to conventional approaches, to capture entire metabolic pathways from metagenomic libraries. The method is compatible with all types of (meta)genomic libraries, and employs a commercially available flow cytometer instead of a custom-made sorting system to detect intracellular or extracellular enzyme activities. This versatile and generic workflow will accelerate experimental exploration campaigns in functional metagenomics and holobiomics studies, to further decipher host-microbiota relationships. Video Abstract.

PNPase is involved in the coordination of mRNA degradation and expression in stationary phase cells of Escherichia coli

Background

Exoribonucleases are crucial for RNA degradation in Escherichia coli but the roles of RNase R and PNPase and their potential overlap in stationary phase are not well characterized. Here, we used a genome-wide approach to determine how RNase R and PNPase affect the mRNA half-lives in the stationary phase. The genome-wide mRNA half-lives were determined by a dynamic analysis of transcriptomes after transcription arrest. We have combined the analysis of mRNA half-lives with the steady-state concentrations (transcriptome) to provide an integrated overview of the in vivo activity of these exoribonucleases at the genome-scale.

Results

The values of mRNA half-lives demonstrated that the mRNAs are very stable in the stationary phase and that the deletion of RNase R or PNPase caused only a limited mRNA stabilization. Intriguingly the absence of PNPase provoked also the destabilization of many mRNAs. These changes in mRNA half-lives in the PNPase deletion strain were associated with a massive reorganization of mRNA levels and also variation in several ncRNA concentrations. Finally, the in vivo activity of the degradation machinery was found frequently saturated by mRNAs in the PNPase mutant unlike in the RNase R mutant, suggesting that the degradation activity is limited by the deletion of PNPase but not by the deletion of RNase R.

Conclusions

This work had identified PNPase as a central player associated with mRNA degradation in stationary phase.

Electronic supplementary material

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

Large-Scale Measurement of mRNA Degradation in Escherichia coli: To Delay or Not to Delay

In this study, we compared different computational methods used for genome-wide determination of mRNA half-lives in Escherichia coli with a special focus on the impact on considering a delay before the onset of mRNA decay after transcription arrest. A wide variety of datasets were analyzed coming from different technical methods for mRNA quantification (microarrays, RNA-seq, and RT-qPCR) and different bacterial growth conditions. The exponential decay of mRNA levels was fitted using both linear and exponential models and with or without a delay. We showed that for all the models, independently of mRNA quantification methods and growth conditions, ignoring the delay resulted in only a modest overestimation of the half-life. For approximately 80% of the mRNAs, differences in mRNA half-life values were less than 34s. The correlation between half-lives estimated with and without a delay was extremely high. However, the slope of the linear regression between the half-lives with and without a delay tended to decrease with the delay. For the few mRNAs for which taking into account the delay influenced the estimated half-life, the impact was dependent on the model and the growth condition. The smallest impact was obtained for the linear model.

Dual role of transcription and transcript stability in the regulation of gene expression in Escherichia coli cells cultured on glucose at different growth rates

Microorganisms extensively reorganize gene expression to adjust growth rate to changes in growth conditions. At the genomic scale, we measured the contribution of both transcription and transcript stability to regulating messenger RNA (mRNA) concentration in Escherichia coli. Transcriptional control was the dominant regulatory process. Between growth rates of 0.10 and 0.63 h(-1), there was a generic increase in the bulk mRNA transcription. However, many transcripts became less stable and the median mRNA half-life decreased from 4.2 to 2.8 min. This is the first evidence that mRNA turnover is slower at extremely low-growth rates. The destabilization of many, but not all, transcripts at high-growth rate correlated with transcriptional upregulation of genes encoding the mRNA degradation machinery. We identified five classes of growth-rate regulation ranging from mainly transcriptional to mainly degradational. In general, differential stability within polycistronic messages encoded by operons does not appear to be affected by growth rate. We show here that the substantial reorganization of gene expression involving downregulation of tricarboxylic acid cycle genes and acetyl-CoA synthetase at high-growth rates is controlled mainly by transcript stability. Overall, our results demonstrate that the control of transcript stability has an important role in fine-tuning mRNA concentration during changes in growth rate.

Functional metagenomics reveals novel pathways of prebiotic breakdown by human gut bacteria

The human intestine hosts a complex bacterial community that plays a major role in nutrition and in maintaining human health. A functional metagenomic approach was used to explore the prebiotic breakdown potential of human gut bacteria, including non-cultivated ones. Two metagenomic libraries, constructed from ileum mucosa and fecal microbiota, were screened for hydrolytic activities on the prebiotic carbohydrates inulin, fructo-oligosaccharides, xylo-oligosaccharides, galacto-oligosaccharides and lactulose. The DNA inserts of 17 clones, selected from the 167 hits that were identified, were pyrosequenced in-depth, yielding in total 407, 420 bp of metagenomic DNA. From these sequences, we discovered novel prebiotic degradation pathways containing carbohydrate transporters and hydrolysing enzymes, for which we provided the first experimental proof of function. Twenty of these proteins are encoded by genes that are also present in the gut metagenome of at least 100 subjects, whatever are their ages or their geographical origin. The sequence taxonomic assignment indicated that still unknown bacteria, for which neither culture conditions nor genome sequence are available, possess the enzymatic machinery to hydrolyse the prebiotic carbohydrates tested. The results expand the vision on how prebiotics are metabolized along the intestine, and open new perspectives for the design of functional foods.

Mining for hemicellulases in the fungus-growing termite Pseudacanthotermes militaris using functional metagenomics

BACKGROUND

The metagenomic analysis of gut microbiomes has emerged as a powerful strategy for the identification of biomass-degrading enzymes, which will be no doubt useful for the development of advanced biorefining processes. In the present study, we have performed a functional metagenomic analysis on comb and gut microbiomes associated with the fungus-growing termite, Pseudacanthotermes militaris.

RESULTS

Using whole termite abdomens and fungal-comb material respectively, two fosmid-based metagenomic libraries were created and screened for the presence of xylan-degrading enzymes. This revealed 101 positive clones, corresponding to an extremely high global hit rate of 0.49%. Many clones displayed either β-d-xylosidase (EC 3.2.1.37) or α-l-arabinofuranosidase (EC 3.2.1.55) activity, while others displayed the ability to degrade AZCL-xylan or AZCL-β-(1,3)-β-(1,4)-glucan. Using secondary screening it was possible to pinpoint clones of interest that were used to prepare fosmid DNA. Sequencing of fosmid DNA generated 1.46 Mbp of sequence data, and bioinformatics analysis revealed 63 sequences encoding putative carbohydrate-active enzymes, with many of these forming parts of sequence clusters, probably having carbohydrate degradation and metabolic functions. Taxonomic assignment of the different sequences revealed that Firmicutes and Bacteroidetes were predominant phyla in the gut sample, while microbial diversity in the comb sample resembled that of typical soil samples. Cloning and expression in E. coli of six enzyme candidates identified in the libraries provided access to individual enzyme activities, which all proved to be coherent with the primary and secondary functional screens.

CONCLUSIONS

This study shows that the gut microbiome of P. militaris possesses the potential to degrade biomass components, such as arabinoxylans and arabinans. Moreover, the data presented suggests that prokaryotic microorganisms present in the comb could also play a part in the degradation of biomass within the termite mound, although further investigation will be needed to clarify the complex synergies that might exist between the different microbiomes that constitute the termitosphere of fungus-growing termites. This study exemplifies the power of functional metagenomics for the discovery of biomass-active enzymes and has provided a collection of potentially interesting biocatalysts for further study.

A high-throughput screening system for the evaluation of biomass-hydrolyzing glycoside hydrolases

To implement a protein engineering strategy for the improvement of enzyme performance on biomass, a straightforward, robust high-throughput method was devised and tested with recombinant GH11 xylanase as acting on wheat straw. The method requires automated liquid handling equipment, but avoids the need for specialized milling and powder weighing devices and the use of labour intensive steps such as manual cutting of pipette tips. After expression in Escherichia coli cells grown in microtiter plates, recombinant xylanase was released into the culture medium and used directly for biomass hydrolysis. Reactions were monitored using a micro-3,5-dinitrosalicylic acid assay. The cumulative error of the method was less than 15%. To validate the method, randomly generated xylanase mutants were analyzed. This allowed the detection of one mutant, which produced a 74% increase in hydrolysis compared to the parental enzyme. Closer analysis revealed that this increase in activity was correlated with a twofold increase in xylanase expression.

Capecitabine plus oxaliplatin (XELOX) versus 5-fluorouracil/folinic acid plus oxaliplatin (FOLFOX-4) as second-line therapy in metastatic colorectal cancer: a randomized phase III noninferiority study

BACKGROUND

To demonstrate the noninferiority of capecitabine plus oxaliplatin (XELOX) versus 5-fluorouracil/folinic acid and oxaliplatin (FOLFOX-4) as second-line therapy in patients with metastatic colorectal cancer after prior irinotecan-based chemotherapy.

PATIENTS AND METHODS

A total of 627 patients were randomly assigned to receive XELOX (n = 313) or FOLFOX-4 (n = 314) following disease progression/recurrence or intolerance to irinotecan-based chemotherapy. The primary end point was progression-free survival (PFS).

RESULTS

PFS for XELOX was noninferior to FOLFOX-4 [hazard ratio (HR) = 0.97; 95% confidence interval (CI) 0.83-1.14] in the intention-to-treat (ITT) population. Median PFS was 4.7 months with XELOX versus 4.8 months with FOLFOX-4. The robustness of the primary analysis was supported by multivariate and subgroup analyses. Median overall survival in the ITT population was 11.9 months with XELOX versus 12.5 months with FOLFOX-4 (HR = 1.02; 95% CI 0.86-1.21). Treatment-related grade 3/4 adverse events occurred in 50% of XELOX- and 65% of FOLFOX-4-treated patients. Whereas grade 3/4 neutropenia (35% versus 5% with XELOX) and febrile neutropenia (4% versus < 1%) were more common with FOLFOX-4, grade 3/4 diarrhea (19% versus 5% with FOLFOX-4) and grade 3 hand-foot syndrome (4% versus < 1%) were more common with XELOX.

CONCLUSION

XELOX is noninferior to FOLFOX-4 when administered as second-line treatment in patients with metastatic colorectal cancer.