Oak genome reveals facets of long lifespan

Oaks are an important part of our natural and cultural heritage. Not only are they ubiquitous in our most common landscapes¹ but they have also supplied human societies with invaluable services, including food and shelter, since prehistoric times². With 450 species spread throughout Asia, Europe and America³, oaks constitute a critical global renewable resource. The longevity of oaks (several hundred years) probably underlies their emblematic cultural and historical importance. Such long-lived sessile organisms must persist in the face of a wide range of abiotic and biotic threats over their lifespans. We investigated the genomic features associated with such a long lifespan by sequencing, assembling and annotating the oak genome. We then used the growing number of whole-genome sequences for plants (including tree and herbaceous species) to investigate the parallel evolution of genomic characteristics potentially underpinning tree longevity. A further consequence of the long lifespan of trees is their accumulation of somatic mutations during mitotic divisions of stem cells present in the shoot apical meristems. Empirical⁴ and modelling⁵ approaches have shown that intra-organismal genetic heterogeneity can be selected for⁶ and provides direct fitness benefits in the arms race with short-lived pests and pathogens through a patchwork of intra-organismal phenotypes⁷. However, there is no clear proof that large-statured trees consist of a genetic mosaic of clonally distinct cell lineages within and between branches. Through this case study of oak, we demonstrate the accumulation and transmission of somatic mutations and the expansion of disease-resistance gene families in trees.

Gut bacteria are critical for optimal muscle function: a potential link with glucose homeostasis

Gut microbiota is involved in the development of several chronic diseases, including diabetes, obesity, and cancer, through its interactions with the host organs. It has been suggested that the cross talk between gut microbiota and skeletal muscle plays a role in different pathological conditions, such as intestinal chronic inflammation and cachexia. However, it remains unclear whether gut microbiota directly influences skeletal muscle function. In this work, we studied the impact of gut microbiota modulation on mice skeletal muscle function and investigated the underlying mechanisms. We determined the consequences of gut microbiota depletion after treatment with a mixture of a broad spectrum of antibiotics for 21 days and after 10 days of natural reseeding. We found that, in gut microbiota-depleted mice, running endurance was decreased, as well as the extensor digitorum longus muscle fatigue index in an ex vivo contractile test. Importantly, the muscle endurance capacity was efficiently normalized by natural reseeding. These endurance changes were not related to variation in muscle mass, fiber typology, or mitochondrial function. However, several pertinent glucose metabolism markers, such as ileum gene expression of short fatty acid chain and glucose transporters G protein-coupled receptor 41 and sodium-glucose cotransporter 1 and muscle glycogen level, paralleled the muscle endurance changes observed after treatment with antibiotics for 21 days and reseeding. Because glycogen is a key energetic substrate for prolonged exercise, modulating its muscle availability via gut microbiota represents one potent mechanism that can contribute to the gut microbiota-skeletal muscle axis. Taken together, our results strongly support the hypothesis that gut bacteria are required for host optimal skeletal muscle function.

Characterization of an extensive rainbow trout miRNA transcriptome by next generation sequencing

Background

MicroRNAs (miRNAs) have emerged as important post-transcriptional regulators of gene expression in a wide variety of physiological processes. They can control both temporal and spatial gene expression and are believed to regulate 30 to 70 % of the genes. Data are however limited for fish species, with only 9 out of the 30,000 fish species present in miRBase. The aim of the current study was to discover and characterize rainbow trout (Oncorhynchus mykiss) miRNAs in a large number of tissues using next-generation sequencing in order to provide an extensive repertoire of rainbow trout miRNAs.

Results

A total of 38 different samples corresponding to 16 different tissues or organs were individually sequenced and analyzed independently in order to identify a large number of miRNAs with high confidence. This led to the identification of 2946 miRNA loci in the rainbow trout genome, including 445 already known miRNAs. Differential expression analysis was performed in order to identify miRNAs exhibiting specific or preferential expression among the 16 analyzed tissues. In most cases, miRNAs exhibit a specific pattern of expression in only a few tissues. The expression data from sRNA sequencing were confirmed by RT-qPCR. In addition, novel miRNAs are described in rainbow trout that had not been previously reported in other species.

Conclusion

This study represents the first characterization of rainbow trout miRNA transcriptome from a wide variety of tissue and sets an extensive repertoire of rainbow trout miRNAs. It provides a starting point for future studies aimed at understanding the roles of miRNAs in major physiological process such as growth, reproduction or adaptation to stress. These rainbow trout miRNAs repertoire provide a novel resource to advance genomic research in salmonid species.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2505-9) contains supplementary material, which is available to authorized users.

Deciphering intra-species bacterial diversity of meat and seafood spoilage microbiota using gyrB amplicon sequencing: A comparative analysis with 16S rDNA V3-V4 amplicon sequencing

Meat and seafood spoilage ecosystems harbor extensive bacterial genomic diversity that is mainly found within a small number of species but within a large number of strains with different spoilage metabolic potential. To decipher the intraspecies diversity of such microbiota, traditional metagenetic analysis using the 16S rRNA gene is inadequate. We therefore assessed the potential benefit of an alternative genetic marker, gyrB, which encodes the subunit B of DNA gyrase, a type II DNA topoisomerase. A comparison between 16S rDNA-based (V3-V4) amplicon sequencing and gyrB-based amplicon sequencing was carried out in five types of meat and seafood products, with five mock communities serving as quality controls. Our results revealed that bacterial richness in these mock communities and food samples was estimated with higher accuracy using gyrB than using16S rDNA. However, for Firmicutes species, 35% of putative gyrB reads were actually identified as sequences of a gyrB paralog, parE, which encodes subunit B of topoisomerase IV; we therefore constructed a reference database of published sequences of both gyrB and pare for use in all subsequent analyses. Despite this co-amplification, the deviation between relative sequencing quantification and absolute qPCR quantification was comparable to that observed for 16S rDNA for all the tested species. This confirms that gyrB can be used successfully alongside 16S rDNA to determine the species composition (richness and evenness) of food microbiota. The major benefit of gyrB sequencing is its potential for improving taxonomic assignment and for further investigating OTU richness at the subspecies level, thus allowing more accurate discrimination of samples. Indeed, 80% of the reads of the 16S rDNA dataset were represented by thirteen 16S rDNA-based OTUs that could not be assigned at the species-level. Instead, these same clades corresponded to 44 gyrB-based OTUs, which differentiated various lineages down to the subspecies level. The increased ability of gyrB-based analyses to track and trace phylogenetically different groups of strains will generate improved resolution and more reliable results for studies of the strains implicated in food processes.

Comparative metatranscriptomic analysis of anaerobic digesters treating anionic surfactant contaminated wastewater

Three distinct biological reactors fed with synthetic medium (UASB_Control), synthetic medium and linear alkylbenzene sulfonate (LAS; UASB_SL), and real laundry wastewater (UASB_LW) were compared using a metatranscriptomic approach to determine putative bioindicator genes and taxonomies associated to all steps of anaerobic LAS biodegradation pathway. A homemade bioinformatics pipeline combined with an R workflow was developed to perform the RNAseq data analysis. UASB_SL and UASB_LW showed similar values of LAS biological degradation (~47%) and removal (53-55%). Rarefaction analysis revealed that 1-2 million reads were sufficient to access the whole functional capacity. In the first step of LAS biodegradation pathway, fumarate reductase subunit C was detected and taxonomically assigned to the genus Syntrophobacter (0.002% - UASB_SL; 0.0015% - UASB_LW; not detected - UASB_Control). In the second step, many enzymes related to beta-oxidation were observed and most of them with low relative abundance in UASB Control and taxonomically related with Smithella, Acinetobacter and Syntrophorhabdus. For the ring cleavage step, the abundance of 6 OCH CoA hydrolase putative gene was ten times higher in UASB_SL and UASB_LW when compared to UASB_Control, and assigned to Desulfomonile and Syntrophorhabdus. Finally, the adenylylsulfate reductase, taxonomically related with Desulfovibrio and Desulfomonile, was observed in the desulfonation step with the highest relative abundance in UASB_LW.

Taxon Appearance From Extraction and Amplification Steps Demonstrates the Value of Multiple Controls in Tick Microbiota Analysis

Background

The development of high-throughput sequencing technologies has substantially improved analysis of bacterial community diversity, composition, and functions. Over the last decade, high-throughput sequencing has been used extensively to identify the diversity and composition of tick microbial communities. However, a growing number of studies are warning about the impact of contamination brought along the different steps of the analytical process, from DNA extraction to amplification. In low biomass samples, e.g., individual tick samples, these contaminants may represent a large part of the obtained sequences, and thus generate considerable errors in downstream analyses and in the interpretation of results. Most studies of tick microbiota either do not mention the inclusion of controls during the DNA extraction or amplification steps, or consider the lack of an electrophoresis signal as an absence of contamination. In this context, we aimed to assess the proportion of contaminant sequences resulting from these steps. We analyzed the microbiota of individual Ixodes ricinus ticks by including several categories of controls throughout the analytical process: homogenization, DNA extraction, and DNA amplification.

Results

Controls yielded a significant number of sequences (1,126-13,198 mean sequences, depending on the control category). Some operational taxonomic units (OTUs) detected in these controls belong to genera reported in previous tick microbiota studies. In this study, these OTUs accounted for 50.9% of the total number of sequences in our samples, and were considered contaminants. Contamination levels (i.e., the percentage of sequences belonging to OTUs identified as contaminants) varied with tick instar and sex: 76.3% of nymphs and 75% of males demonstrated contamination over 50%, while most females (65.7%) had rates lower than 20%. Contamination mainly corresponded to OTUs detected in homogenization and extraction reagent controls, highlighting the importance of carefully controlling these steps.

Conclusion

Here, we showed that contaminant OTUs from sample laboratory processing steps can represent more than half the total sequence yield in sequencing runs, and lead to unreliable results when characterizing tick microbial communities. We thus strongly advise the routine use of negative controls in tick microbiota studies, and more generally in studies involving low biomass samples.

Mining the equine gut metagenome: poorly-characterized taxa associated with cardiovascular fitness in endurance athletes

Emerging evidence indicates that the gut microbiome contributes to endurance exercise performance. Still, the extent of its functional and metabolic potential remains unknown. Using elite endurance horses as a model system for exercise responsiveness, we built an integrated horse gut gene catalog comprising ~25 million unique genes and 372 metagenome-assembled genomes. This catalog represents 4179 genera spanning 95 phyla and functional capacities primed to exploit energy from dietary, microbial, and host resources. The holo-omics approach shows that gut microbiomes enriched in Lachnospiraceae taxa are negatively associated with cardiovascular capacity. Conversely, more complex and functionally diverse microbiomes are associated with higher glucose concentrations and reduced accumulation of long-chain acylcarnitines and non-esterified fatty acids in plasma, suggesting increased ß-oxidation capacity in the mitochondria. In line with this hypothesis, more fit athletes show upregulation of mitochondrial-related genes involved in energy metabolism, biogenesis, and Ca²⁺ cytosolic transport, all of which are necessary to improve aerobic work power, spare glycogen usage, and enhance cardiovascular capacity. The results identify an associative link between endurance performance and gut microbiome composition and gene function, laying the basis for nutritional interventions that could benefit horse athletes.

An integrated gene catalog of the gut microbiome in elite endurance horses is build. The holo-omics analyses identify an associative link between endurance performance and gut microbiome composition and gene function.

Insights into the cyanosphere: capturing the respective metabolisms of cyanobacteria and chemotrophic bacteria in natural conditions?

Specific interactions have been highlighted between cyanobacteria and chemotrophic bacteria within the cyanosphere, suggesting that nutrients recycling could be optimized by cyanobacteria/bacteria exchanges. In order to determine the respective metabolic roles of the cyanobacterial and bacterial consortia (microbiome), a day-night metatranscriptomic analysis was performed on Dolichospermum sp. (N₂ -fixer) and Microcystis sp. (non N₂ -fixer) natural blooms occurring successively within a French peri-urban lake. The taxonomical and functional analysis of the metatranscriptoms have highlighted specific association of bacteria within the cyanosphere, driven by the cyanobacteria identity, without strongly modifying the functional composition of the microbiomes, suggesting functional redundancy within the cyanosphere. Moreover, the functional composition of these active communities was driven by the living mode. During the two successive bloom events, it appeared that NH₄ ⁺ (newly fixed and/or allochthonous) was preferentially transformed into amino acids for the both the microbiome and the cyanobacteria, while phosphate metabolism was enhanced, suggesting that due to a high cellular growth, P limitation might take place within the cyanosphere consortium.

FROGS: a powerful tool to analyse the diversity of fungi with special management of internal transcribed spacers

Fungi are present in all environments. They fulfil important ecological functions and play a crucial role in the food industry. Their accurate characterization is thus indispensable, particularly through metabarcoding. The most frequently used markers to monitor fungi are ITSs. These markers are the best documented in public databases but have one main weakness: polymerase chain reaction amplification may produce non-overlapping reads in a significant fraction of the fungi. When these reads are filtered out, traditional metabarcoding pipelines lose part of the information and consequently produce biased pictures of the composition and structure of the environment under study. We developed a solution that enables processing of the entire set of reads including both overlapping and non-overlapping, thus providing a more accurate picture of fungal communities. Our comparative tests using simulated and real data demonstrated the effectiveness of our solution, which can be used by both experts and non-specialists on a command line or through the Galaxy-based web interface.

The ileal fungal microbiota is altered in Crohn's disease and is associated with the disease course

Introduction

Fungal microbiota's involvement in the pathogenesis of Crohn's disease (CD) is incompletely understood. The terminal ileum is a predilection site both for primary involvement and recurrences of CD. We, therefore, assessed the mucosa-associated mycobiota in the inflamed and non-inflamed ileum in patients with CD.

Methods

The mucosa-associated mycobiota was assessed by ITS2 sequencing in a total of 168 biopsies sampled 5 and 15 cm proximal of the ileocecal valve or ileocolic anastomosis in 44 CD patients and 40 healthy controls (HC). CD patients with terminal ileitis, with endoscopic inflammation at 5 cm and normal mucosa at 15 cm and no history of upper CD involvement, were analyzed separately. The need for additional CD treatment the year following biopsy collection was recorded.

Results

CD patients had reduced mycobiota evenness, increased Basidiomycota/Ascomycota ratio, and reduced abundance of Chytridiomycota compared to HC. The mycobiota of CD patients were characterized by an expansion of Malassezia and a depletion of Saccharomyces, along with increased abundances of Candida albicans and Malassezia restricta. Malassezia was associated with the need for treatment escalation during follow-up. Current anti-TNF treatment was associated with lower abundances of Basidiomycota. The alpha diversity of the inflamed and proximal non-inflamed mucosa within the same patients was similar. However, the inflamed mucosa had a more dysbiotic composition with increased abundances of Candida sake and reduced abundances of Exophiala equina and Debaryomyces hansenii.

Conclusions

The ileal mucosa-associated mycobiota in CD patients is altered compared to HC. The mycobiota in the inflamed and proximal non-inflamed ileum within the same patients harbor structural differences which may play a role in the CD pathogenesis. Increased abundance of Malassezia was associated with an unfavorable disease course.

The Rhinobiome of Exacerbated Wheezers and Asthmatics: Insights From a German Pediatric Exacerbation Network

Although the nose, as a gateway for organism–environment interactions, may have a key role in asthmatic exacerbation, the rhinobiome of exacerbated children with asthma was widely neglected to date. The aim of this study is to understand the microbiome, the microbial immunology, and the proteome of exacerbated children and adolescents with wheeze and asthma. Considering that a certain proportion of wheezers may show a progression to asthma, the comparison of both groups provides important information regarding clinical and phenotype stratification. Thus, deep nasopharyngeal swab specimens, nasal epithelial spheroid (NAEsp) cultures, and blood samples of acute exacerbated wheezers (WH), asthmatics (AB), and healthy controls (HC) were used for culture (n = 146), 16 S-rRNA gene amplicon sequencing (n = 64), and proteomic and cytokine analyses. Interestingly, Proteobacteria were over-represented in WH, whereas Firmicutes and Bacteroidetes were associated with AB. In contrast, Actinobacteria commonly colonized HCs. Moreover, Staphylococcaceae, Enterobacteriaceae, Burkholderiaceae, Xanthobacteraceae, and Sphingomonadaceae were significantly more abundant in AB compared to WH and HC. The α-diversity analyses demonstrated an increase of bacterial abundance levels in atopic AB and a decrease in WH samples. Microbiome profiles of atopic WH differed significantly from atopic AB, whereby atopic samples of WH were more homogeneous than those of non-atopic subjects. The NAEsp bacterial exposure experiments provided a disrupted epithelial cell integrity, a cytokine release, and cohort-specific proteomic differences especially for Moraxella catarrhalis cultures. This comprehensive dataset contributes to a deeper insight into the poorly understood plasticity of the nasal microbiota, and, in particular, may enforce our understanding in the pathogenesis of asthma exacerbation in childhood.

A Probiotic Mixture Induces Anxiolytic- and Antidepressive-Like Effects in Fischer and Maternally Deprived Long Evans Rats

A role of the gut microbiota in psychiatric disorders is supported by a growing body of literature. The effects of a probiotic mixture of four bacterial strains were studied in two models of anxiety and depression, naturally stress-sensitive Fischer rats and Long Evans rats subjected to maternal deprivation. Rats chronically received either the probiotic mixture (1.10⁹ CFU/day) or the vehicle. Anxiety- and depressive-like behaviors were evaluated in several tests. Brain monoamine levels and gut RNA expression of tight junction proteins (Tjp) and inflammatory markers were quantified. The gut microbiota was analyzed in feces by 16S rRNA gene sequencing. Untargeted metabolite analysis reflecting primary metabolism was performed in the cecal content and in serum. Fischer rats treated with the probiotic mixture manifested a decrease in anxiety-like behaviors, in the immobility time in the forced swimming test, as well as in levels of dopamine and its major metabolites, and those of serotonin metabolites in the hippocampus and striatum. In maternally deprived Long Evans rats treated with the probiotic mixture, the number of entries into the central area in the open-field test was increased, reflecting an anxiolytic effect. The probiotic mixture increased Tjp1 and decreased Ifnγ mRNA levels in the ileum of maternally deprived rats. In both models, probiotic supplementation changed the proportions of several Operational Taxonomic Units (OTU) in the gut microbiota, and the levels of certain cecal and serum metabolites were correlated with behavioral changes. Chronic administration of the tested probiotic mixture can therefore beneficially affect anxiety- and depressive-like behaviors in rats, possibly owing to changes in the levels of certain metabolites, such as 21-deoxycortisol, and changes in brain monoamines.

Large-scale multivariate dataset on the characterization of microbiota diversity, microbial growth dynamics, metabolic spoilage volatilome and sensorial profiles of two industrially produced meat products subjected to changes in lactate concentration and packaging atmosphere

Data in this article provide detailed information on the diversity of bacterial communities present on 576 samples of raw pork or poultry sausages produced industrially in 2017. Bacterial growth dynamics and diversity were monitored throughout the refrigerated storage period to estimate the impact of packaging atmosphere and the use of potassium lactate as chemical preservative. The data include several types of analysis aiming at providing a comprehensive microbial ecology of spoilage during storage and how the process parameters do influence this phenomenon. The analysis includes: the gas content in packaging, pH, chromametric measurements, plate counts (total mesophilic aerobic flora and lactic acid bacteria), sensorial properties of the products, meta-metabolomic quantification of volatile organic compounds and bacterial community metagenetic analysis. Bacterial diversity was monitored using two types of amplicon sequencing (16S rRNA and GyrB encoding genes) at different time points for the different conditions (576 samples for gyrB and 436 samples for 16S rDNA). Sequencing data were generated by using Illumina MiSeq. The sequencing data have been deposited in the bioproject PRJNA522361. Samples accession numbers vary from SAMN10964863 to SAMN10965438 for gyrB amplicon and from SAMN10970131 to SAMN10970566 for 16S.

Detection of an amplification bias associated to Leuconostocaceae family with a universal primer routinely used for monitoring microbial community structures within food products

Objectives

Sequencing of 16S rDNA V3–V4 region is widely applied for food community profiling. However, two different universal forward primers (named here MUYZER-primer1 and KLINDWORTH-primer2) targeting an identical conservative sequence upstream of the V3 region of 16S rRNA gene, and only distinguished by a single mismatch are both used. This study was carried out to compare whether the accuracy of food microbiota analysis would depend on the choice of one of these two primers.

Results

Alignment of both primers with common food-borne bacteria 16S sequences revealed that the mismatch between both primers might specifically affect the amplification of Leuconostoc, Oenococcus and Fructobacillus species but not Weissella species. Food products containing either Leuconostoc and/or Weissella were selected for a detection test. As expected from our in silico analysis, our study showed that this mismatch induced a strong biased amplification specifically associated to the OTUs belonging to the genus Leuconostoc but not to the genus Weissella. In presence of Muyzer-primer1, none of the sequences expected for Leuconostoc genus was detected whereas those sequences were correctly amplified with Klindworth-primer2. Since Leuconostoc is an important genus in food, agro-environments and in digestive tract of animals, we recommend that Muyzer-primer1 should thus be abandoned for the bacterial characterization of their associated microbiota.

Microbial community dispersal from wheat grains to sourdoughs : a contribution of participatory research

Understanding microbial dispersal is critical to understand the dynamics and evolution of microbial communities. However, microbial dispersal is difficult to study because of uncertainty about their vectors of migration. This applies to both microbial communities in natural and human-associated environments. Here, we studied microbial dispersal along the sourdoughs bread making chain using a participatory research approach. Sourdough is a naturally fermented mixture of flour and water. It hosts a community of bacteria and yeasts whose origins are only partially known. We analysed the potential of wheat grains and flour to serve as an inoculum for sourdough microbial communities using 16S rDNA and ITS1 metabarcoding. First, in an experiment involving farmers, a miller and bakers, we followed the microbiota from grains to newly initiated and propagated sourdoughs. Second, we compared the microbiota of 46 sourdough samples collected everywhere in France, and of the flour used for their backslopping. The core microbiota detected on the seeds, in the flour and in the sourdough was composed mainly of microbes known to be associated with plants and not living in sourdoughs. No sourdough yeast species were detected on grains and flours. Sourdough lactic acid bacteria were rarely found in flour. When they were, they did not have the same amplicon sequence variant (ASV) as found in the corresponding sourdough. However, the low sequencing depth for bacteria in flour did not allow us to draw definitive conclusion. Thus, our results showed that sourdough yeasts did not come from flour, and suggest that neither do sourdough LAB.

First step of odorant detection in the olfactory epithelium and olfactory preferences differ according to the microbiota profile in mice

We have previously provided the first evidence that the microbiota modulates the physiology of the olfactory epithelium using germfree mice. The extent to which changes to the olfactory system depend on the microbiota is still unknown. In the present work, we explored if different microbiota would differentially impact olfaction. We therefore studied the olfactory function of three groups of mice of the same genetic background, whose parents had been conventionalized before mating with microbiota from three different mouse strains. Caecal short chain fatty acids profiles and 16S rRNA gene sequencing ascertained that gut microbiota differed between the three groups. We then used a behavioural test to measure the attractiveness of various odorants and observed that the three groups of mice differed in their attraction towards odorants. Their olfactory epithelium properties, including electrophysiological responses recorded by electro-olfactograms and expression of genes related to the olfactory transduction pathway, also showed several differences. Overall, our data demonstrate that differences in gut microbiota profiles are associated with differences in olfactory preferences and in olfactory epithelium functioning.

Virulent Phages Isolated from a Smear-Ripened Cheese Are Also Detected in Reservoirs of the Cheese Factory

Smear-ripened cheeses host complex microbial communities that play a crucial role in the ripening process. Although bacteriophages have been frequently isolated from dairy products, their diversity and ecological role in such this type of cheese remain underexplored. In order to fill this gap, the main objective of this study was to isolate and characterize bacteriophages from the rind of a smear-ripened cheese. Thus, viral particles extracted from the cheese rind were tested through a spot assay against a collection of bacteria isolated from the same cheese and identified by sequencing the full-length small subunit ribosomal RNA gene. In total, five virulent bacteriophages infecting Brevibacterium aurantiacum, Glutamicibacter arilaitensis, Leuconostoc falkenbergense and Psychrobacter aquimaris species were obtained. All exhibit a narrow host range, being only able to infect a few cheese-rind isolates within the same species. The complete genome of each phage was sequenced using both Nanopore and Illumina technologies, assembled and annotated. A sequence comparison with known phages revealed that four of them may represent at least new genera. The distribution of the five virulent phages into the dairy-plant environment was also investigated by PCR, and three potential reservoirs were identified. This work provides new knowledge on the cheese rind viral community and an overview of the distribution of phages within a cheese factory.

Investigation of the Phageome and Prophages in French Cider, a Fermented Beverage

Phageomes are known to play a key role in the functioning of their associated microbial communities. The phageomes of fermented foods have not been studied thoroughly in fermented foods yet, and even less in fermented beverages. Two approaches were employed to investigate the presence of phages in cider, a fermented beverage made from apple, during a fermentation process of two cider tanks, one from an industrial producer and one from a hand-crafted producer. The phageome (free lytic phages) was explored in cider samples with several methodological developments for total phage DNA extraction, along with single phage isolation. Concentration methods, such as tangential flow filtration, flocculation and classical phage concentration methods, were employed and tested to extract free phage particles from cider. This part of the work revealed a very low occurrence of free lytic phage particles in cider. In parallel, a prophage investigation during the fermentation process was also performed using a metagenomic approach on the total bacterial genomic DNA. Prophages in bacterial metagenomes in the two cider tanks seemed also to occur in low abundance, as a total of 1174 putative prophages were identified in the two tanks overtime, and only two complete prophages were revealed. Prophage occurrence was greater at the industrial producer than at the hand-crafted producer, and different dynamics of prophage trends were also observed during fermentation. This is the first report dealing with the investigation of the phageome and of prophages throughout a fermentation process of a fermented beverage.

Comparison of Fecal Sample Collection Methods for Microbial Analysis Embedded within Colorectal Cancer Screening Programs

BACKGROUND

Colorectal cancer screening programs with fecal sample collection may provide a platform for population-based gut microbiome disease research. We investigated sample collection and storage method impact on the accuracy and stability of the V3-V4 region of the 16S rRNA genes and bacterial quantity across seven different collection methods [i.e., no solution, two specimen collection cards, and four types of fecal immunochemical test (FIT) used in four countries] among 19 healthy volunteers.

METHODS

Intraclass correlation coefficients (ICC) were calculated for the relative abundance of the top three phyla, the most abundant genera, alpha diversity metrics, and the first principal coordinates of the beta diversity matrices to estimate the stability of microbial profiles after storage for 7 days at room temperature, 4°C or 30°C, and after screening for the presence of occult blood in the stool. In addition, accuracy was estimated for samples frozen immediately compared to samples with no solution (i.e., the putative gold standard).

RESULTS

When compared with the putative gold standard, we observed significant variation for all collection methods. However, interindividual variability was much higher than the variability introduced by the collection method. Stability ICCs were high (≥0.75) for FIT tubes that underwent colorectal cancer screening procedures. The relative abundance of Actinobacteria (0.65) was an exception and was lower for different FIT tubes stored at 30°C (range, 0.41-0.90) and room temperature (range, 0.06-0.94).

CONCLUSIONS

Paper-based collection cards and different types of FIT are acceptable tools for microbiome measurements.

IMPACT

Our findings inform on the utility of commonly used fecal sample collection methods for developing microbiome-focused cohorts nested within screening programs.

Integration of Multiomic Data to Characterize the Influence of Milk Fat Composition on Cantal-Type Cheese Microbiota

A previous study identified differences in rind aspects between Cantal-type cheeses manufactured from the same skimmed milk, supplemented with cream derived either from pasture-raised cows (P) or from cows fed with maize silage (M). Using an integrated analysis of multiomic data, the present study aimed at investigating potential correlations between cream origin and metagenomic, lipidomic and volatolomic profiles of these Cantal cheeses. Fungal and bacterial communities of cheese cores and rinds were characterized using DNA metabarcoding at different ripening times. Lipidome and volatolome were obtained from the previous study at the end of ripening. Rind microbial communities, especially fungal communities, were influenced by cream origin. Among bacteria, Brachybacterium were more abundant in P-derived cheeses than in M-derived cheeses after 90 and 150 days of ripening. Sporendonema casei, a yeast added as a ripening starter during Cantal manufacture, which contributes to rind typical aspect, had a lower relative abundance in P-derived cheeses after 150 days of ripening. Relative abundance of this fungus was highly negatively correlated with concentrations of C18 polyunsaturated fatty acids and to concentrations of particular volatile organic compounds, including 1-pentanol and 3-methyl-2-pentanol. Overall, these results evidenced original interactions between milk fat composition and the development of fungal communities in cheeses.

In-depth characterization of a selection of gut commensal bacteria reveals their functional capacities to metabolize dietary carbohydrates with prebiotic potential

The microbial utilization of dietary carbohydrates is closely linked to the pivotal role of the gut microbiome in human health. Inherent to the modulation of complex microbial communities, a prebiotic implies the selective utilization of a specific substrate, relying on the metabolic capacities of targeted microbes. In this study, we investigated the metabolic capacities of 17 commensal bacteria of the human gut microbiome toward dietary carbohydrates with prebiotic potential. First, in vitro experiments allowed the classification of bacterial growth and fermentation profiles in response to various carbon sources, including agave inulin, corn fiber, polydextrose, and citrus pectin. The influence of phylogenetic affiliation appeared to statistically outweigh carbon sources in determining the degree of carbohydrate utilization. Second, we narrowed our focus on six commensal bacteria representative of the Bacteroidetes and Firmicutes phyla to perform an untargeted high-resolution liquid chromatography-mass spectrometry metabolomic analysis: Bacteroides xylanisolvens, Bacteroides thetaiotaomicron, Bacteroides intestinalis, Subdoligranulum variabile, Roseburia intestinalis, and Eubacterium rectale exhibited distinct metabolomic profiles in response to different carbon sources. The relative abundance of bacterial metabolites was significantly influenced by dietary carbohydrates, with these effects being strain-specific and/or carbohydrate-specific. Particularly, the findings indicated an elevation in short-chain fatty acids and other metabolites, including succinate, gamma-aminobutyric acid, and nicotinic acid. These metabolites were associated with putative health benefits. Finally, an RNA-Seq transcriptomic approach provided deeper insights into the underlying mechanisms of carbohydrate metabolization. Restricting our focus on four commensal bacteria, including B. xylanisolvens, B. thetaiotaomicron, S. variabile, and R. intestinalis, carbon sources did significantly modulate the level of bacterial genes related to the enzymatic machinery involved in the metabolization of dietary carbohydrates. This study provides a holistic view of the molecular strategies induced during the dynamic interplay between dietary carbohydrates with prebiotic potential and gut commensal bacteria.

IMPORTANCE

This study explores at a molecular level the interactions between commensal health-relevant bacteria and dietary carbohydrates holding prebiotic potential. We showed that prebiotic breakdown involves the specific activation of gene expression related to carbohydrate metabolism. We also identified metabolites produced by each bacteria that are potentially related to our digestive health. The characterization of the functional activities of health-relevant bacteria toward prebiotic substances can yield a better application of prebiotics in clinical interventions and personalized nutrition. Overall, this study highlights the importance of identifying the impact of prebiotics at a low resolution of the gut microbiota to characterize the activities of targeted bacteria that can play a crucial role in our health.

Holistic integration of omics data reveals the drivers that shape the ecology of microbial meat spoilage scenarios

Background

The use of omics data for monitoring the microbial flow of fresh meat products along a production line and the development of spoilage prediction tools from these data is a promising but challenging task. In this context, we produced a large multivariate dataset (over 600 samples) obtained on the production lines of two similar types of fresh meat products (poultry and raw pork sausages). We describe a full analysis of this dataset in order to decipher how the spoilage microbial ecology of these two similar products may be shaped differently depending on production parameter characteristics.

Methods

Our strategy involved a holistic approach to integrate unsupervised and supervised statistical methods on multivariate data (OTU-based microbial diversity; metabolomic data of volatile organic compounds; sensory measurements; growth parameters), and a specific selection of potential uncontrolled (initial microbiota composition) or controlled (packaging type; lactate concentration) drivers.

Results

Our results demonstrate that the initial microbiota, which is shown to be very different between poultry and pork sausages, has a major impact on the spoilage scenarios and on the effect that a downstream parameter such as packaging type has on the overall evolution of the microbial community. Depending on the process, we also show that specific actions on the pork meat (such as deboning and defatting) elicit specific food spoilers such as Dellaglioa algida, which becomes dominant during storage. Finally, ecological network reconstruction allowed us to map six different metabolic pathways involved in the production of volatile organic compounds involved in spoilage. We were able connect them to the different bacterial actors and to the influence of packaging type in an overall view. For instance, our results demonstrate a new role of Vibrionaceae in isopropanol production, and of Latilactobacillus fuchuensis and Lactococcus piscium in methanethiol/disylphide production. We also highlight a possible commensal behavior between Leuconostoc carnosum and Latilactobacillus curvatus around 2,3-butanediol metabolism.

Conclusion

We conclude that our holistic approach combined with large-scale multi-omic data was a powerful strategy to prioritize the role of production parameters, already known in the literature, that shape the evolution and/or the implementation of different meat spoilage scenarios.