Insight into the Role of Gut Microbiota in Duchenne Muscular Dystrophy: An Age-Related Study in Mdx Mice

Dystrophin deficiency alters the sarcolemma structure, leading to muscle dystrophy, muscle disuse, and ultimately death. Beyond limb muscle deficits, patients with Duchenne muscular dystrophy have numerous transit disorders. Many studies have highlighted the strong relationship between gut microbiota and skeletal muscle. The aims of this study were: i) to characterize the gut microbiota composition over time up to 1 year in dystrophin-deficient mdx mice, and ii) to analyze the intestine structure and function and expression of genes linked to bacterial-derived metabolites in ileum, blood, and skeletal muscles to study interorgan interactions. Mdx mice displayed a significant reduction in the overall number of different operational taxonomic units and their abundance (α-diversity). Mdx genotype predicted 20% of β-diversity divergence, with a large taxonomic modification of Actinobacteria, Proteobacteria, Tenericutes, and Deferribacteres phyla and the included genera. Interestingly, mdx intestinal motility and gene expressions of tight junction and Ffar2 receptor were down-regulated in the ileum. Concomitantly, circulating inflammatory markers related to gut microbiota (tumor necrosis factor, IL-6, monocyte chemoattractant protein-1) and muscle inflammation Tlr4/Myd88 pathway (Toll-like receptor 4, which recognizes pathogen-associated molecular patterns) were up-regulated. Finally, in mdx mice, adiponectin was reduced in blood and its receptor modulated in muscles. This study highlights a specific gut microbiota composition and highlights interorgan interactions in mdx physiopathology with gut microbiota as the potential central metabolic organ.

Microbial communities of a variety of 75 homemade fermented vegetables

Fermentation is an ancient practice of food preservation. Fermented vegetables are popular in Eastern European and Asian countries. They have received a growing interest in Western countries, where they are mainly manufactured at domestic and artisanal scales and poorly characterized. Our aim was to investigate the microbial communities and the safety of French homemade fermented vegetables, in the frame of a citizen science project. Fermented vegetables and the data associated with their manufacture were collected from citizens and characterized for pH, NaCl concentration, and microbiology by culturomics and 16S DNA metabarcoding analysis. Lactic acid bacteria (LAB) and yeast isolates were identified by 16S rRNA gene sequencing and D1/D2 domains of the large subunit of the rRNA gene, respectively. The 75 collected samples contained 23 types of vegetables, mainly cabbage, followed by carrots and beets, and many mixtures of vegetables. They were 2 weeks to 4 years old, and their median pH was 3.56, except for two samples with a pH over 4.5. LAB represented the dominant viable bacteria. LAB concentrations ranged from non-detectable values to 8.7 log colony-forming units (CFU)/g and only depended on the age of the samples, with the highest most frequently observed in the youngest samples (<100 days). The 93 LAB isolates identified belonged to 23 species, the two mains being Lactiplantibacillus pentosus/plantarum and Levilactobacillus brevis. The other microbial groups enumerated (total aerobic bacteria, halotolerant bacteria, Gram-negative bacteria, and acetic acid bacteria) generally showed lower concentrations compared to LAB concentrations. No pathogenic bacteria were detected. Viable yeasts were observed in nearly half the samples, at concentrations reaching up to 8.0 log CFU/g. The 33 yeast clones identified belonged to 16 species. Bacterial metabarcoding showed two main orders, namely, Lactobacillales (i.e., LAB, 79% of abundance, 177 of the 398 total ASVs) and Enterobacterales (19% of abundance, 191 ASVs). Fifteen LAB genera were identified, with Lactiplantibacillus and Levilactobacillus as the most abundant, with 41 and 12% of total reads, respectively. Enterobacterales members were mainly represented by Enterobacteriaceae and Yersiniaceae. This study is the first wide description of the microbiota of a large variety of homemade fermented vegetables and documents their safety.

Microbiota members from body sites of dairy cows are largely shared within individual hosts throughout lactation but sharing is limited in the herd

BACKGROUND

Host-associated microbes are major determinants of the host phenotypes. In the present study, we used dairy cows with different scores of susceptibility to mastitis with the aim to explore the relationships between microbiota composition and different factors in various body sites throughout lactation as well as the intra- and inter-animal microbial sharing.

RESULTS

Microbiotas from the mouth, nose, vagina and milk of 45 lactating dairy cows were characterized by metataxonomics at four time points during the first lactation, from 1-week pre-partum to 7 months post-partum. Each site harbored a specific community that changed with time, likely reflecting physiological changes in the transition period and changes in diet and housing. Importantly, we found a significant number of microbes shared among different anatomical sites within each animal. This was between nearby anatomic sites, with up to 32% of the total number of Amplicon Sequence Variants (ASVs) of the oral microbiota shared with the nasal microbiota but also between distant ones (e.g. milk with nasal and vaginal microbiotas). In contrast, the share of microbes between animals was limited (< 7% of ASVs shared by more than 50% of the herd for a given site and time point). The latter widely shared ASVs were mainly found in the oral and nasal microbiotas. These results thus indicate that despite a common environment and diet, each animal hosted a specific set of bacteria, supporting a tight interplay between each animal and its microbiota. The score of susceptibility to mastitis was slightly but significantly related to the microbiota associated to milk suggesting a link between host genetics and microbiota.

CONCLUSIONS

This work highlights an important sharing of microbes between relevant microbiotas involved in health and production at the animal level, whereas the presence of common microbes was limited between animals of the herd. This suggests a host regulation of body-associated microbiotas that seems to be differently expressed depending on the body site, as suggested by changes in the milk microbiota that were associated to genotypes of susceptibility to mastitis.

Diversity of novel archaeal viruses infecting methanogens discovered through coupling of stable isotope probing and metagenomics

Diversity of viruses infecting non-extremophilic archaea has been grossly understudied. This is particularly the case for viruses infecting methanogenic archaea, key players in the global carbon biogeochemical cycle. Only a dozen of methanogenic archaeal viruses have been isolated so far. In the present study, we implemented an original coupling between stable isotope probing and complementary shotgun metagenomic analyses to identify viruses of methanogens involved in the bioconversion of formate, which was used as the sole carbon source in batch anaerobic digestion microcosms. Under our experimental conditions, the microcosms were dominated by methanogens belonging to the order Methanobacteriales (Methanobacterium and Methanobrevibacter genera). Metagenomic analyses yielded several previously uncharacterized viral genomes, including a complete genome of a head-tailed virus (class Caudoviricetes, proposed family Speroviridae, Methanobacterium host) and several near-complete genomes of spindle-shaped viruses. The two groups of viruses are predicted to infect methanogens of the Methanobacterium and Methanosarcina genera and represent two new virus families. The metagenomics results are in good agreement with the electron microscopy observations, which revealed the dominance of head-tailed virus-like particles and the presence of spindle-shaped particles. The present study significantly expands the knowledge on the viral diversity of viruses of methanogens.

Genomic evidence for global ocean plankton biogeography shaped by large-scale current systems

Biogeographical studies have traditionally focused on readily visible organisms, but recent technological advances are enabling analyses of the large-scale distribution of microscopic organisms, whose biogeographical patterns have long been debated. Here we assessed the global structure of plankton geography and its relation to the biological, chemical, and physical context of the ocean (the ‘seascape’) by analyzing metagenomes of plankton communities sampled across oceans during the Tara Oceans expedition, in light of environmental data and ocean current transport. Using a consistent approach across organismal sizes that provides unprecedented resolution to measure changes in genomic composition between communities, we report a pan-ocean, size-dependent plankton biogeography overlying regional heterogeneity. We found robust evidence for a basin-scale impact of transport by ocean currents on plankton biogeography, and on a characteristic timescale of community dynamics going beyond simple seasonality or life history transitions of plankton.

Oceans are brimming with life invisible to our eyes, a myriad of species of bacteria, viruses and other microscopic organisms essential for the health of the planet. These ‘marine plankton’ are unable to swim against currents and should therefore be constantly on the move, yet previous studies have suggested that distinct species of plankton may in fact inhabit different oceanic regions. However, proving this theory has been challenging; collecting plankton is logistically difficult, and it is often impossible to distinguish between species simply by examining them under a microscope. However, within the last decade, a research schooner called Tara has travelled the globe to gather thousands of plankton samples. At the same time, advances in genomics have made it possible to identify species based only on fragments of their DNA sequence.

To understand the hidden geography of plankton communities in Earth’s oceans, Richter et al. pored over DNA from the Tara Oceans expedition. This revealed that, despite being unable to resist the flow of water, various planktonic species which live close to the surface manage to occupy distinct, stable provinces shaped by currents. Different sizes of plankton are distributed in different sized provinces, with the smallest organisms tending to inhabit the smallest areas. Comparing DNA similarities and speeds of currents at the ocean surface revealed how these might stretch and mix plankton communities.

Plankton play a critical role in the health of the ocean and the chemical cycles of planet Earth. These results could allow deeper investigation by marine modellers, ecologists, and evolutionary biologists. Meanwhile, work is already underway to investigate how climate change might impact this hidden geography.

Prebiotic Supplementation During Gestation Induces a Tolerogenic Environment and a Protective Microbiota in Offspring Mitigating Food Allergy

Food allergy is associated with alterations in the gut microbiota, epithelial barrier, and immune tolerance. These dysfunctions are observed within the first months of life, indicating that early intervention is crucial for disease prevention. Preventive nutritional strategies with prebiotics are an attractive option, as prebiotics such as galacto-oligosaccharides and inulin can promote tolerance, epithelial barrier reinforcement, and gut microbiota modulation. Nonetheless, the ideal period for intervention remains unknown. Here, we investigated whether galacto-oligosaccharide/inulin supplementation during gestation could protect offspring from wheat allergy development in BALB/cJRj mice. We demonstrated that gestational prebiotic supplementation promoted the presence of beneficial strains in the fecal microbiota of dams during gestation and partially during mid-lactation. This specific microbiota was transferred to their offspring and maintained to adulthood. The presence of B and T regulatory immune cell subsets was also increased in the lymph nodes of offspring born from supplemented mothers, suggestive of a more tolerogenic immune environment. Indeed, antenatal prebiotic supplementation reduced the development of wheat allergy symptoms in offspring. Our study thus demonstrates that prebiotic supplementation during pregnancy induces, in the offspring, a tolerogenic environment and a microbial imprint that mitigates food allergy development.

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.

Glucose but Not Fructose Alters the Intestinal Paracellular Permeability in Association With Gut Inflammation and Dysbiosis in Mice

A causal correlation between the metabolic disorders associated with sugar intake and disruption of the gastrointestinal (GI) homeostasis has been suggested, but the underlying mechanisms remain unclear. To unravel these mechanisms, we investigated the effect of physiological amounts of fructose and glucose on barrier functions and inflammatory status in various regions of the GI tract and on the cecal microbiota composition. C57BL/6 mice were fed chow diet and given 15% glucose or 15% fructose in drinking water for 9 weeks. We monitored caloric intake, body weight, glucose intolerance, and adiposity. The intestinal paracellular permeability, cytokine, and tight junction protein expression were assessed in the jejunum, cecum, and colon. In the cecum, the microbiota composition was determined. Glucose-fed mice developed a marked increase in total adiposity, glucose intolerance, and paracellular permeability in the jejunum and cecum while fructose absorption did not affect any of these parameters. Fructose-fed mice displayed increased circulation levels of IL6. In the cecum, both glucose and fructose intake were associated with an increase in Il13, Ifnγ, and Tnfα mRNA and MLCK protein levels. To clarify the relationships between monosaccharides and barrier function, we measured the permeability of Caco-2 cell monolayers in response to IFNγ+TNFα in the presence of glucose or fructose. In vitro, IFNγ+TNFα-induced intestinal permeability increase was less pronounced in response to fructose than glucose. Mice treated with glucose showed an enrichment of Lachnospiracae and Desulfovibrionaceae while the fructose increased relative abundance of Lactobacillaceae. Correlations between pro-inflammatory cytokine gene expression and bacterial abundance highlighted the potential role of members of Desulfovibrio and Lachnospiraceae NK4A136 group genera in the inflammation observed in response to glucose intake. The increase in intestinal inflammation and circulating levels of IL6 in response to fructose was observed in the absence of intestinal permeability modification, suggesting that the intestinal permeability alteration does not precede the onset of metabolic outcome (low-grade inflammation, hyperglycemia) associated with chronic fructose consumption. The data also highlight the deleterious effects of glucose on gut barrier function along the GI tract and suggest that Desulfovibrionaceae and Lachnospiraceae play a key role in the onset of GI inflammation in response to glucose.

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.

Microbial and functional characterization of an allochthonous consortium applied to hydrogen production from Citrus Peel Waste in batch reactor in optimized conditions

Energy recovery from lignocellulosic waste has been studied as an alternative to the problem of inappropriate waste disposal. The present study aimed at characterizing the microbial community and the functional activity of reactors applied to H₂ production through lignocellulosic waste fermentation in optimized conditions. The latter were identified by means of Rotational Central Composite Design (RCCD), applied to optimize allochthonous inoculum concentration (2.32-5.68 gTVS/L of granular anaerobic sludge), pH (4.32-7.68) and Citrus Peel Waste (CPW) concentration (1.55-28.45 g/L). After validation, the conditions identified for optimal H₂ production were 4 gSTV/L of allochthonous inoculum, 29.8 g/L of CPW (substrate) and initial pH of 8.98. In these conditions, 48.47 mmol/L of H₂ was obtained, which is 3.64 times higher than the concentration in unoptimized conditions (13.31 mmol H₂/L using 15 g/L of CPW, 2 gTVS/L of allochthonous inoculum, pH 7.0). Acetogenesis was the predominant pathway, and maximal concentrations of 3,731 mg/L of butyric acid and 3,516 mg/L of acetic acid were observed. Regarding the metataxonomic profile, Clostridium genus was dramatically favored in the optimized condition (79.78%) when compared to the allochthonous inoculum (0.43%). It was possible to identify several genes related to H₂ (i.e dehydrogenases) and volatile fatty acids (VFA) production and with cellulose degradation, especially some CAZymes from the classes Auxiliary Activities, Glycoside Hydrolases and Glycosyl Transferase. By means of differential gene expression it was observed that cellulose degradation and acetic acid production pathways were overabundant in samples from the optimized reactors, highlighting endo-β-1,4-glucanase/cellulose, endo-β-1,4-xylanase, β-glucosidase, β-mannosidase, cellulose β-1,4-cellobiosidase, cellobiohydrolase, and others, as main the functions.

Temporal patterns in Ixodes ricinus microbial communities: an insight into tick-borne microbe interactions

BACKGROUND

Ticks transmit pathogens of medical and veterinary importance and are an increasing threat to human and animal health. Assessing disease risk and developing new control strategies requires identifying members of the tick-borne microbiota as well as their temporal dynamics and interactions.

METHODS

Using high-throughput sequencing, we studied the Ixodes ricinus microbiota and its temporal dynamics. 371 nymphs were monthly collected during three consecutive years in a peri-urban forest. After a Poisson lognormal model was adjusted to our data set, a principal component analysis, sparse network reconstruction, and differential analysis allowed us to assess seasonal and monthly variability of I. ricinus microbiota and interactions within this community.

RESULTS

Around 75% of the detected sequences belonged to five genera known to be maternally inherited bacteria in arthropods and to potentially circulate in ticks: Candidatus Midichloria, Rickettsia, Spiroplasma, Arsenophonus and Wolbachia. The structure of the I. ricinus microbiota varied over time with interannual recurrence and seemed to be mainly driven by OTUs commonly found in the environment. Total network analysis revealed a majority of positive partial correlations. We identified strong relationships between OTUs belonging to Wolbachia and Arsenophonus, evidence for the presence of the parasitoid wasp Ixodiphagus hookeri in ticks. Other associations were observed between the tick symbiont Candidatus Midichloria and pathogens belonging to Rickettsia. Finally, more specific network analyses were performed on TBP-infected samples and suggested that the presence of pathogens belonging to the genera Borrelia, Anaplasma and Rickettsia may disrupt microbial interactions in I. ricinus.

CONCLUSIONS

We identified the I. ricinus microbiota and documented marked shifts in tick microbiota dynamics over time. Statistically, we showed strong relationships between the presence of specific pathogens and the structure of the I. ricinus microbiota. We detected close links between some tick symbionts and the potential presence of either pathogenic Rickettsia or a parasitoid in ticks. These new findings pave the way for the development of new strategies for the control of ticks and tick-borne diseases. Video abstract.

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.

The Poisson-Lognormal Model as a Versatile Framework for the Joint Analysis of Species Abundances

Joint Species Distribution Models (JSDM) provide a general multivariate framework to study the joint abundances of all species from a community. JSDM account for both structuring factors (environmental characteristics or gradients, such as habitat type or nutrient availability) and potential interactions between the species (competition, mutualism, parasitism, etc.), which is instrumental in disentangling meaningful ecological interactions from mere statistical associations. Modeling the dependency between the species is challenging because of the count-valued nature of abundance data and most JSDM rely on Gaussian latent layer to encode the dependencies between species in a covariance matrix. The multivariate Poisson-lognormal (PLN) model is one such model, which can be viewed as a multivariate mixed Poisson regression model. Inferring such models raises both statistical and computational issues, many of which were solved in recent contributions using variational techniques and convex optimization tools. The PLN model turns out to be a versatile framework, within which a variety of analyses can be performed, including multivariate sample comparison, clustering of sites or samples, dimension reduction (ordination) for visualization purposes, or inferring interaction networks. This paper presents the general PLN framework and illustrates its use on a series a typical experimental datasets. All the models and methods are implemented in the R package PLNmodels, available from cran.r-project.org.

Exploring short k-mer profiles in cells and mobile elements from Archaea highlights the major influence of both the ecological niche and evolutionary history

BACKGROUND

K-mer-based methods have greatly advanced in recent years, largely driven by the realization of their biological significance and by the advent of next-generation sequencing. Their speed and their independence from the annotation process are major advantages. Their utility in the study of the mobilome has recently emerged and they seem a priori adapted to the patchy gene distribution and the lack of universal marker genes of viruses and plasmids. To provide a framework for the interpretation of results from k-mer based methods applied to archaea or their mobilome, we analyzed the 5-mer DNA profiles of close to 600 archaeal cells, viruses and plasmids. Archaea is one of the three domains of life. Archaea seem enriched in extremophiles and are associated with a high diversity of viral and plasmid families, many of which are specific to this domain. We explored the dataset structure by multivariate and statistical analyses, seeking to identify the underlying factors.

RESULTS

For cells, the 5-mer profiles were inconsistent with the phylogeny of archaea. At a finer taxonomic level, the influence of the taxonomy and the environmental constraints on 5-mer profiles was very strong. These two factors were interdependent to a significant extent, and the respective weights of their contributions varied according to the clade. A convergent adaptation was observed for the class Halobacteria, for which a strong 5-mer signature was identified. For mobile elements, coevolution with the host had a clear influence on their 5-mer profile. This enabled us to identify one previously known and one new case of recent host transfer based on the atypical composition of the mobile elements involved. Beyond the effect of coevolution, extrachromosomal elements strikingly retain the specific imprint of their own viral or plasmid taxonomic family in their 5-mer profile.

CONCLUSION

This specific imprint confirms that the evolution of extrachromosomal elements is driven by multiple parameters and is not restricted to host adaptation. In addition, we detected only recent host transfer events, suggesting the fast evolution of short k-mer profiles. This calls for caution when using k-mers for host prediction, metagenomic binning or phylogenetic reconstruction.

Unraveling the history of the genus Gallus through whole genome sequencing

The genus Gallus is distributed across a large part of Southeast Asia and has received special interest because the domestic chicken, Gallus gallus domesticus, has spread all over the world and is a major protein source for humans. There are four species: the red junglefowl (G. gallus), the green junglefowl (G. varius), the Lafayette's junglefowl (G. lafayettii) and the grey junglefowl (G. sonneratii). The aim of this study is to reconstruct the history of these species by a whole genome sequencing approach and resolve inconsistencies between well supported topologies inferred using different data and methods. Using deep sequencing, we identified over 35 million SNPs and reconstructed the phylogeny of the Gallus genus using both distance (BioNJ) and maximum likelihood (ML) methods. We observed discrepancies according to reconstruction methods and genomic components. The two most supported topologies were previously reported and were discriminated by using phylogenetic and gene flow analyses, based on ABBA statistics. Terminology fix requested by the deputy editor led to support a scenario with G. gallus as the earliest branching lineage of the Gallus genus, instead of G. varius. We discuss the probable causes for the discrepancy. A likely one is that G. sonneratii samples from parks or private collections are all recent hybrids, with roughly 10% of their autosomal genome originating from G. gallus. The removal of those regions is needed to provide reliable data, which was not done in previous studies. We took care of this and additionally included two wild G. sonneratii samples from India, showing no trace of introgression. This reinforces the importance of carefully selecting and validating samples and genomic components in phylogenomics.

Function-Driven Design of Lactic Acid Bacteria Co-cultures to Produce New Fermented Food Associating Milk and Lupin

Designing bacterial co-cultures adapted to ferment mixes of vegetal and animal resources for food diversification and sustainability is becoming a challenge. Among bacteria used in food fermentation, lactic acid bacteria (LAB) are good candidates, as they are used as starter or adjunct in numerous fermented foods, where they allow preservation, enhanced digestibility, and improved flavor. We developed here a strategy to design LAB co-cultures able to ferment a new food made of bovine milk and lupin flour, consisting in: (i) in silico preselection of LAB species for targeted carbohydrate degradation; (ii) in vitro screening of 97 strains of the selected species for their ability to ferment carbohydrates and hydrolyze proteins from milk and lupin and clustering strains that displayed similar phenotypes; and (iii) assembling strains randomly sampled from clusters that showed complementary phenotypes. The designed co-cultures successfully expressed the targeted traits i.e., hydrolyzed proteins and degraded raffinose family oligosaccharides of lupin and lactose of milk in a large range of concentrations. They also reduced an off-flavor-generating volatile, hexanal, and produced various desirable flavor compounds. Most of the strains in co-cultures achieved higher cell counts than in monoculture, suggesting positive interactions. This work opens new avenues for the development of innovative fermented food products based on functionally complementary strains in the world-wide context of diet diversification.

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.

SimkaMin: fast and resource frugal de novo comparative metagenomics

MOTIVATION

De novo comparative metagenomics is one of the most straightforward ways to analyze large sets of metagenomic data. Latest methods use the fraction of shared k-mers to estimate genomic similarity between read sets. However, those methods, while extremely efficient, are still limited by computational needs for practical usage outside of large computing facilities.

RESULTS

We present SimkaMin, a quick comparative metagenomics tool with low disk and memory footprints, thanks to an efficient data subsampling scheme used to estimate Bray-Curtis and Jaccard dissimilarities. One billion metagenomic reads can be analyzed in <3 min, with tiny memory (1.09 GB) and disk (≈0.3 GB) requirements and without altering the quality of the downstream comparative analyses, making of SimkaMin a tool perfectly tailored for very large-scale metagenomic projects.

AVAILABILITY AND IMPLEMENTATION

https://github.com/GATB/simka.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Detection of selection signatures in Limousin cattle using whole-genome resequencing

Limousin, a renowned beef breed originating from central France, has been selectively bred over the last 100 years to improve economically important traits. We used whole-genome sequencing data from 10 unrelated Limousin bull calves to detect polymorphisms and identify regions under selection. A total of 13 943 766 variants were identified. Moreover, 311 852 bi-allelic SNPs and 92 229 indels located on autosomes were fixed for the alternative allele in all sequenced animals, including the previously reported missense deleterious F94L mutation in MSTN. We performed a whole-genome screen to discover genomic regions with excess homozygosity, using the pooled heterozygosity score and identified 171 different candidate selective sweeps. In total, 68 candidate genes were found in only 57 of these regions, indicating that a large fraction of the genome under selection might lie in non-coding regions and suggesting that a majority of adaptive mutations might be regulatory in nature. Many QTL were found within candidate selective sweep regions, including QTL associated with shear force or carcass weight. Among the putative selective sweeps, we located genes (MSTN, NCKAP5, RUNX2) that potentially contribute to important phenotypes in Limousin. Several candidate regions and genes under selection were also found in previous genome-wide selection scans performed in Limousin. In addition, we were able to pinpoint candidate causative regulatory polymorphisms in GRIK3 and RUNX2 that might have been under selection. Our results will contribute to improved understanding of the mechanisms and targets of artificial selection and will facilitate the interpretation of GWASs performed in Limousin.

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.

Incorporating Phylogenetic Information in Microbiome Differential Abundance Studies Has No Effect on Detection Power and FDR Control

We consider the problem of incorporating evolutionary information (e.g., taxonomic or phylogenic trees) in the context of metagenomics differential analysis. Recent results published in the literature propose different ways to leverage the tree structure to increase the detection rate of differentially abundant taxa. Here, we propose instead to use a different hierarchical structure, in the form of a correlation-based tree, as it may capture the structure of the data better than the phylogeny. We first show that the correlation tree and the phylogeny are significantly different before turning to the impact of tree choice on detection rates. Using synthetic data, we show that the tree does have an impact: smoothing p-values according to the phylogeny leads to equal or inferior rates as smoothing according to the correlation tree. However, both trees are outperformed by the classical, non-hierarchical, Benjamini–Hochberg (BH) procedure in terms of detection rates. Other procedures may use the hierarchical structure with profit but do not control the False Discovery Rate (FDR) a priori and remain inferior to a classical Benjamini–Hochberg procedure with the same nominal FDR. On real datasets, no hierarchical procedure had significantly higher detection rate that BH. Intuition advocates that the use of hierarchical structures should increase the detection rate of differentially abundant taxa in microbiome studies. However, our results suggest that current hierarchical procedures are still inferior to standard methods and more effective procedures remain to be invented.

Shifts in the temperature-sensitive periods for spring phenology in European beech and pedunculate oak clones across latitudes and over recent decades

Spring phenology of temperate trees has advanced worldwide in response to global warming. However, increasing temperatures may not necessarily lead to further phenological advance, especially in the warmer latitudes because of insufficient chilling and/or shorter day length. Determining the start of the forcing phase, that is, when buds are able to respond to warmer temperatures in spring, is therefore crucial to predict how phenology will change in the future. In this study, we used 4,056 leaf-out date observations during the period 1969-2017 for clones of European beech (Fagus sylvatica L.) and pedunculate oak (Quercus robur L.) planted in 63 sites covering a large latitudinal gradient (from Portugal ~41°N to Norway ~63°N) at the International Phenological Gardens in order to (a) evaluate how the sensitivity periods to forcing and chilling have changed with climate warming, and (b) test whether consistent patterns occur along biogeographical gradients, that is, from colder to warmer environments. Partial least squares regressions suggest that the length of the forcing period has been extended over the recent decades with climate warming in the colder latitudes but has been shortened in the warmer latitudes for both species, with a more pronounced shift for beech. We attribute the lengthening of the forcing period in the colder latitudes to earlier opportunities with temperatures that can promote bud development. In contrast, at warmer or oceanic climates, the beginning of the forcing period has been delayed, possibly due to insufficient chilling. However, in spite of a later beginning of the forcing period, spring phenology has continued to advance at these areas due to a faster satisfaction of heat requirements induced by climate warming. Overall, our results support that ongoing climate warming will have different effects on the spring phenology of forest trees across latitudes due to the interactions between chilling, forcing and photoperiod.

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.

Genetic and metabolic signatures of Salmonella enterica subsp. enterica associated with animal sources at the pangenomic scale

BACKGROUND

Salmonella enterica subsp. enterica is a public health issue related to food safety, and its adaptation to animal sources remains poorly described at the pangenome scale. Firstly, serovars presenting potential mono- and multi-animal sources were selected from a curated and synthetized subset of Enterobase. The corresponding sequencing reads were downloaded from the European Nucleotide Archive (ENA) providing a balanced dataset of 440 Salmonella genomes in terms of serovars and sources (i). Secondly, the coregenome variants and accessory genes were detected (ii). Thirdly, single nucleotide polymorphisms and small insertions/deletions from the coregenome, as well as the accessory genes were associated to animal sources based on a microbial Genome Wide Association Study (GWAS) integrating an advanced correction of the population structure (iii). Lastly, a Gene Ontology Enrichment Analysis (GOEA) was applied to emphasize metabolic pathways mainly impacted by the pangenomic mutations associated to animal sources (iv).

RESULTS

Based on a genome dataset including Salmonella serovars from mono- and multi-animal sources (i), 19,130 accessory genes and 178,351 coregenome variants were identified (ii). Among these pangenomic mutations, 52 genomic signatures (iii) and 9 over-enriched metabolic signatures (iv) were associated to avian, bovine, swine and fish sources by GWAS and GOEA, respectively.

CONCLUSIONS

Our results suggest that the genetic and metabolic determinants of Salmonella adaptation to animal sources may have been driven by the natural feeding environment of the animal, distinct livestock diets modified by human, environmental stimuli, physiological properties of the animal itself, and work habits for health protection of livestock.

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.

Fructose malabsorption induces cholecystokinin expression in the ileum and cecum by changing microbiota composition and metabolism

Current fructose consumption levels often overwhelm the intestinal capacity to absorb fructose. We investigated the impact of fructose malabsorption on intestinal endocrine function and addressed the role of the microbiota in this process. To answer this question, a mouse model of moderate fructose malabsorption [ketohexokinase mutant (KHK)^-/-] and wild-type (WT) littermate mice were used and received a 20%-fructose (KHK-F and WT-F) or 20%-glucose diet. Cholecystokinin (Cck) mRNA and protein expression in the ileum and cecum, as well as preproglucagon (Gcg) and neurotensin (Nts) mRNA expression in the cecum, increased in KHK-F mice. In KHK-F mice, triple-label immunohistochemistry showed major up-regulation of CCK in enteroendocrine cells (EECs) that were glucagon-like peptide-1 (GLP-1)⁺/Peptide YY (PYY^-) in the ileum and colon and GLP-1^-/PYY^- in the cecum. The cecal microbiota composition was drastically modified in the KHK-F in association with an increase in glucose, propionate, succinate, and lactate concentrations. Antibiotic treatment abolished fructose malabsorption-dependent induction of cecal Cck mRNA expression and, in mouse GLUTag and human NCI-H716 cells, Cck mRNA expression levels increased in response to propionate, both suggesting a microbiota-dependent process. Fructose reaching the lower intestine can modify the composition and metabolism of the microbiota, thereby stimulating the production of CCK from the EECs possibly in response to propionate.-Zhang, X., Grosfeld, A., Williams, E., Vasiliauskas, D., Barretto, S., Smith, L., Mariadassou, M., Philippe, C., Devime, F., Melchior, C., Gourcerol, G., Dourmap, N., Lapaque, N., Larraufie, P., Blottière, H. M., Herberden, C., Gerard, P., Rehfeld, J. F., Ferraris, R. P., Fritton, J. C., Ellero-Simatos, S., Douard, V. Fructose malabsorption induces cholecystokinin expression in the ileum and cecum by changing microbiota composition and metabolism.

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.

Inference of Adaptive Shifts for Multivariate Correlated Traits

To study the evolution of several quantitative traits, the classical phylogenetic comparative framework consists of a multivariate random process running along the branches of a phylogenetic tree. The Ornstein-Uhlenbeck (OU) process is sometimes preferred to the simple Brownian motion (BM) as it models stabilizing selection toward an optimum. The optimum for each trait is likely to be changing over the long periods of time spanned by large modern phylogenies. Our goal is to automatically detect the position of these shifts on a phylogenetic tree, while accounting for correlations between traits, which might exist because of structural or evolutionary constraints. We show that, in the presence of shifts, phylogenetic Principal Component Analysis fails to decorrelate traits efficiently, so that any method aiming at finding shifts needs to deal with correlation simultaneously. We introduce here a simplification of the full multivariate OU model, named scalar OU, which allows for noncausal correlations and is still computationally tractable. We extend the equivalence between the OU and a BM on a rescaled tree to our multivariate framework. We describe an Expectation-Maximization (EM) algorithm that allows for a maximum likelihood estimation of the shift positions, associated with a new model selection criterion, accounting for the identifiability issues for the shift localization on the tree. The method, freely available as an R-package (PhylogeneticEM) is fast, and can deal with missing values. We demonstrate its efficiency and accuracy compared to another state-of-the-art method ($\ell$1ou) on a wide range of simulated scenarios and use this new framework to reanalyze recently gathered data sets on New World Monkeys and Anolis lizards.

Phages infecting Faecalibacterium prausnitzii belong to novel viral genera that help to decipher intestinal viromes

BACKGROUND

Viral metagenomic studies have suggested a role for bacteriophages in intestinal dysbiosis associated with several human diseases. However, interpretation of viral metagenomic studies is limited by the lack of knowledge of phages infecting major human gut commensal bacteria, such as Faecalibacterium prausnitzii, a bacterial symbiont repeatedly found depleted in inflammatory bowel disease (IBD) patients. In particular, no complete genomes of phages infecting F. prausnitzii are present in viral databases.

METHODS

We identified 18 prophages in 15 genomes of F. prausnitzii, used comparative genomics to define eight phage clades, and annotated the genome of the type phage of each clade. For two of the phages, we studied prophage induction in vitro and in vivo in mice. Finally, we aligned reads from already published viral metagenomic data onto the newly identified phages.

RESULTS

We show that each phage clade represents a novel viral genus and that a surprisingly large fraction of them (10 of the 18 phages) codes for a diversity-generating retroelement, which could contribute to their adaptation to the digestive tract environment. We obtained either experimental or in silico evidence of activity for at least one member of each genus. In addition, four of these phages are either significantly more prevalent or more abundant in stools of IBD patients than in those of healthy controls.

CONCLUSION

Since IBD patients generally have less F. prausnitzii in their microbiota than healthy controls, the higher prevalence or abundance of some of its phages may indicate that they are activated during disease. This in turn suggests that phages could trigger or aggravate F. prausnitzii depletion in patients. Our results show that prophage detection in sequenced strains of the microbiota can usefully complement viral metagenomic studies.

OH1 from Orf Virus: A New Tyrosine Phosphatase that Displays Distinct Structural Features and Triple Substrate Specificity

Viral tyrosine phosphatases such as VH1 from Vaccinia and Variola virus are recognized as important effectors of host-pathogen interactions. While proteins sharing sequence to VH1 have been identified in other viruses, their structural and functional characterization is not known. In this work, we determined the crystal structure of the VH1 homolog in the Orf virus, herein named OH1. Similarly to Variola and Vaccinia VH1, the structure of OH1 shows a dimer with the typical dual-specificity phosphatase fold. In contrast to VH1, the OH1 dimer is covalently stabilized by a disulfide bond involving residue Cys15 in the N-terminal helix alpha-1 of both monomers, and Cys15 is a conserved residue within the Parapoxvirus genus. The in vitro functional characterization confirms that OH1 is a dual-specificity phosphatase and reveals its ability to dephosphorylate phosphatidylinositol 3,5-bisphosphate, a new activity potentially relevant in phosphoinositide recycling during virion maturation.

Analysis of large 16S rRNA Illumina data sets: Impact of singleton read filtering on microbial community description

Next-generation sequencing technologies give access to large sets of data, which are extremely useful in the study of microbial diversity based on 16S rRNA gene. However, the production of such large data sets is not only marred by technical biases and sequencing noise but also increases computation time and disc space use. To improve the accuracy of OTU predictions and overcome both computations, storage and noise issues, recent studies and tools suggested removing all single reads and low abundant OTUs, considering them as noise. Although the effect of applying an OTU abundance threshold on α- and β-diversity has been well documented, the consequences of removing single reads have been poorly studied. Here, we test the effect of singleton read filtering (SRF) on microbial community composition using in silico simulated data sets as well as sequencing data from synthetic and real communities displaying different levels of diversity and abundance profiles. Scalability to large data sets is also assessed using a complete MiSeq run. We show that SRF drastically reduces the chimera content and computational time, enabling the analysis of a complete MiSeq run in just a few minutes. Moreover, SRF accurately determines the actual community diversity: the differences in α- and β-community diversity obtained with SRF and standard procedures are much smaller than the intrinsic variability of technical and biological replicates.

Identification of proteins involved in the anti-inflammatory properties of Propionibacterium freudenreichii by means of a multi-strain study

Propionibacterium freudenreichii, a dairy starter, can reach a population of almost 10⁹ propionibacteria per gram in Swiss-type cheese at the time of consumption. Also consumed as a probiotic, it displays strain-dependent anti-inflammatory properties mediated by surface proteins that induce IL-10 in leukocytes. We selected 23 strains with varied anti-inflammatory potentials in order to identify the protein(s) involved. After comparative genomic analysis, 12 of these strains were further analysed by surface proteomics, eight of them being further submitted to transcriptomics. The omics data were then correlated to the anti-inflammatory potential evaluated by IL-10 induction. This comparative omics strategy highlighted candidate genes that were further subjected to gene-inactivation validation. This validation confirmed the contribution of surface proteins, including SlpB and SlpE, two proteins with SLH domains known to mediate non-covalent anchorage to the cell-wall. Interestingly, HsdM3, predicted as cytoplasmic and involved in DNA modification, was shown to contribute to anti-inflammatory activity. Finally, we demonstrated that a single protein cannot explain the anti-inflammatory properties of a strain. These properties therefore result from different combinations of surface and cytoplasmic proteins, depending on the strain. Our enhanced understanding of the molecular bases for immunomodulation will enable the relevant screening for bacterial resources with anti-inflammatory properties.

Unprecedented large inverted repeats at the replication terminus of circular bacterial chromosomes suggest a novel mode of chromosome rescue

The first Lactobacillus delbrueckii ssp. bulgaricus genome sequence revealed the presence of a very large inverted repeat (IR), a DNA sequence arrangement which thus far seemed inconceivable in a non-manipulated circular bacterial chromosome, at the replication terminus. This intriguing observation prompted us to investigate if similar IRs could be found in other bacteria. IRs with sizes varying from 38 to 76 kbp were found at the replication terminus of all 5 L. delbrueckii ssp. bulgaricus chromosomes analysed, but in none of 1373 other chromosomes. They represent the first naturally occurring very large IRs detected in circular bacterial genomes. A comparison of the L. bulgaricus replication terminus regions and the corresponding regions without IR in 5 L. delbrueckii ssp. lactis genomes leads us to propose a model for the formation and evolution of the IRs. The DNA sequence data are consistent with a novel model of chromosome rescue after premature replication termination or irreversible chromosome damage near the replication terminus, involving mechanisms analogous to those proposed in the formation of very large IRs in human cancer cells. We postulate that the L. delbrueckii ssp. bulgaricus-specific IRs in different strains derive from a single ancestral IR of at least 93 kbp.

Olfactory epithelium changes in germfree mice

Intestinal epithelium development is dramatically impaired in germfree rodents, but the consequences of the absence of microbiota have been overlooked in other epithelia. In the present study, we present the first description of the bacterial communities associated with the olfactory epithelium and explored differences in olfactory epithelium characteristics between germfree and conventional, specific pathogen-free, mice. While the anatomy of the olfactory epithelium was not significantly different, we observed a thinner olfactory cilia layer along with a decreased cellular turn-over in germfree mice. Using electro-olfactogram, we recorded the responses of olfactory sensitive neuronal populations to various odorant stimulations. We observed a global increase in the amplitude of responses to odorants in germfree mice as well as altered responses kinetics. These changes were associated with a decreased transcription of most olfactory transduction actors and of olfactory xenobiotic metabolising enzymes. Overall, we present here the first evidence that the microbiota modulates the physiology of olfactory epithelium. As olfaction is a major sensory modality for most animal species, the microbiota may have an important impact on animal physiology and behaviour through olfaction alteration.

Beneficial metabolic effects of selected probiotics on diet-induced obesity and insulin resistance in mice are associated with improvement of dysbiotic gut microbiota

Alterations in gut microbiota composition and diversity were suggested to play a role in the development of obesity, a chronic subclinical inflammatory condition. We here evaluated the impact of oral consumption of a monostrain or multi-strain probiotic preparation in high-fat diet-induced obese mice. We observed a strain-specific effect and reported dissociation between the capacity of probiotics to dampen adipose tissue inflammation and to limit body weight gain. A multi-strain mixture was able to improve adiposity, insulin resistance and dyslipidemia through adipose tissue immune cell-remodelling, mainly affecting macrophages. At the gut level, the mixture modified the uptake of fatty acids and restored the expression level of the short-chain fatty acid receptor GPR43. These beneficial effects were associated with changes in the microbiota composition, such as the restoration of the abundance of Akkermansia muciniphila and Rikenellaceae and the decrease of other taxa like Lactobacillaceae. Using an in vitro gut model, we further showed that the probiotic mixture favours the production of butyrate and propionate. Our findings provide crucial clues for the design and use of more efficient probiotic preparations in obesity management and may bring new insights into the mechanisms by which host-microbe interactions govern such protective effects.

Quantitative prediction of genome-wide resource allocation in bacteria

Predicting resource allocation between cell processes is the primary step towards decoding the evolutionary constraints governing bacterial growth under various conditions. Quantitative prediction at genome-scale remains a computational challenge as current methods are limited by the tractability of the problem or by simplifying hypotheses. Here, we show that the constraint-based modeling method Resource Balance Analysis (RBA), calibrated using genome-wide absolute protein quantification data, accurately predicts resource allocation in the model bacterium Bacillus subtilis for a wide range of growth conditions. The regulation of most cellular processes is consistent with the objective of growth rate maximization except for a few suboptimal processes which likely integrate more complex objectives such as coping with stressful conditions and survival. As a proof of principle by using simulations, we illustrated how calibrated RBA could aid rational design of strains for maximizing protein production, offering new opportunities to investigate design principles in prokaryotes and to exploit them for biotechnological applications.

Mutations and genomic islands can explain the strain dependency of sugar utilization in 21 strains of Propionibacterium freudenreichii

Background

Propionibacterium freudenreichii (PF) is an actinobacterium used in cheese technology and for its probiotic properties. PF is also extremely adaptable to several ecological niches and can grow on a variety of carbon and nitrogen sources. The aim of this work was to discover the genetic basis for strain-dependent traits related to its ability to use specific carbon sources. High-throughput sequencing technologies were ideal for this purpose as they have the potential to decipher genomic diversity at a moderate cost.

Results

21 strains of PF were sequenced and the genomes were assembled de novo. Scaffolds were ordered by comparison with the complete reference genome CIRM-BIA1, obtained previously using traditional Sanger sequencing. Automatic functional annotation and manual curation were performed. Each gene was attributed to either the core genome or an accessory genome. The ability of the 21 strains to degrade 50 different sugars was evaluated. Thirty-three sugars were degraded by none of the sequenced strains whereas eight sugars were degraded by all of them. The corresponding genes were present in the core genome. Lactose, melibiose and xylitol were only used by some strains. In this case, the presence/absence of genes responsible for carbon uptake and degradation correlated well with the phenotypes, with the exception of xylitol. Furthermore, the simultaneous presence of these genes was in line the metabolic pathways described previously in other species. We also considered the genetic origin (transduction, rearrangement) of the corresponding genomic islands. Ribose and gluconate were degraded to a greater or lesser extent (quantitative phenotype) by some strains. For these sugars, the phenotypes could not be explained by the presence/absence of a gene but correlated with the premature appearance of a stop codon interrupting protein synthesis and preventing the catabolism of corresponding carbon sources.

Conclusion

These results illustrate (i) the power of correlation studies to discover the genetic basis of binary strain-dependent traits, and (ii) the plasticity of PF chromosomes, probably resulting from horizontal transfers, duplications, transpositions and an accumulation of mutations. Knowledge of the genetic basis of nitrogen and sugar degradation opens up new strategies for the screening of PF strain collections to enable optimum cheese starter, probiotic and white biotechnology applications.

Electronic supplementary material

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

Identification of amino acids in mitochondrially encoded proteins that correlate with lifespan

Animals show a huge diversity in their lifespan that can vary from a few weeks to over a hundred years in vertebrates. Size is a key element in this variation and the positive correlation between size and maximum lifespan can be observed in each class of vertebrate. Some groups and species clearly stand out in this size-lifespan relationship and the ones with exceptionally long lifespan have been studied to understand the biological causes of their low aging rate. Among the potential explanations of animals' lifespan variations, mitochondria and mitochondrially encoded genes have drawn attention because of their importance in the aging process. To understand both the extent of lifespan variations and their dependence to genes and amino acid variations in mitochondrial genes and DNA (mtDNA), we analyze in a systematic way all 13 proteins encoded by mitochondria in all vertebrates for which we had information on weight, maximum lifespan and mtDNA sequence. This comparison allows us to visualize positions, and even specific amino acids, in these sequences that correlate with lifespan. With this approach, we draw a map of 356 amino acid residues, at 296 positions within the sequence, that correlate with longer or shorter lifespan. We also compared this map with the human mitochondrial polymorphism to determine its potential as a predictive tool.

Local and systemic immune mechanisms underlying the anti-colitis effects of the dairy bacterium Lactobacillus delbrueckii

Several probiotic bacteria have been proposed for treatment or prevention of inflammatory bowel diseases (IBD), showing a protective effect in animal models of experimental colitis and for some of them also in human clinical trials. While most of these probiotic bacteria are isolated from the digestive tract, we recently reported that a Lactobacillus strain isolated from cheese, L. delbrueckii subsp. lactis CNRZ327 (Lb CNRZ327), also possesses anti-inflammatory effects in vitro and in vivo, demonstrating that common dairy bacteria may be useful in the treatment or prevention of IBD. Here, we studied the mechanisms underlying the protective effects of Lb CNRZ327 in vivo, in a mouse dextran sodium sulfate (DSS) colitis model. During colitis, Lb CNRZ327 modulated the production of TGF-β, IL-6, and IL-12 in colonic tissue and of TGF-β and IL-6 in the spleen, and caused an expansion of CD4+Foxp3+ regulatory T cells in the cecal lymph nodes. Moreover, a strong tendency to CD4+Foxp3+ expansion was also observed in the spleen. The results of this study for the first time show that orally administered dairy lactobacilli can not only modulate mucosal but also systemic immune responses and constitute an effective treatment of IBD.

Analyses of the radiation of birnaviruses from diverse host phyla and of their evolutionary affinities with other double-stranded RNA and positive strand RNA viruses using robust structure-based multiple sequence alignments and advanced phylogenetic methods

Background

Birnaviruses form a distinct family of double-stranded RNA viruses infecting animals as different as vertebrates, mollusks, insects and rotifers. With such a wide host range, they constitute a good model for studying the adaptation to the host. Additionally, several lines of evidence link birnaviruses to positive strand RNA viruses and suggest that phylogenetic analyses may provide clues about transition.

Results

We characterized the genome of a birnavirus from the rotifer Branchionus plicalitis. We used X-ray structures of RNA-dependent RNA polymerases and capsid proteins to obtain multiple structure alignments that allowed us to obtain reliable multiple sequence alignments and we employed “advanced” phylogenetic methods to study the evolutionary relationships between some positive strand and double-stranded RNA viruses. We showed that the rotifer birnavirus genome exhibited an organization remarkably similar to other birnaviruses. As this host was phylogenetically very distant from the other known species targeted by birnaviruses, we revisited the evolutionary pathways within the Birnaviridae family using phylogenetic reconstruction methods. We also applied a number of phylogenetic approaches based on structurally conserved domains/regions of the capsid and RNA-dependent RNA polymerase proteins to study the evolutionary relationships between birnaviruses, other double-stranded RNA viruses and positive strand RNA viruses.

Conclusions

We show that there is a good correlation between the phylogeny of the birnaviruses and that of their hosts at the phylum level using the RNA-dependent RNA polymerase (genomic segment B) on the one hand and a concatenation of the capsid protein, protease and ribonucleoprotein (genomic segment A) on the other hand. This correlation tends to vanish within phyla. The use of advanced phylogenetic methods and robust structure-based multiple sequence alignments allowed us to obtain a more accurate picture (in terms of probability of the tree topologies) of the evolutionary affinities between double-stranded RNA and positive strand RNA viruses. In particular, we were able to show that there exists a good statistical support for the claims that dsRNA viruses are not monophyletic and that viruses with permuted RdRps belong to a common evolution lineage as previously proposed by other groups. We also propose a tree topology with a good statistical support describing the evolutionary relationships between the Picornaviridae, Caliciviridae, Flaviviridae families and a group including the Alphatetraviridae, Nodaviridae, Permutotretraviridae, Birnaviridae, and Cystoviridae families.

Taxon influence index: assessing taxon-induced incongruities in phylogenetic inference

Understanding the evolutionary history of species is at the core of molecular evolution and is done using several inference methods. The critical issue is to quantify the uncertainty of the inference. The posterior probabilities in Bayesian phylogenetic inference and the bootstrap values in frequentist approaches measure the variability of the estimates due to the sampling of sites from genes and the sampling of genes from genomes. However, they do not measure the uncertainty due to taxon sampling. Taxa that experienced molecular homoplasy, recent selection, a spur of evolution, and so forth may disrupt the inference and cause incongruences in the estimated phylogeny. We define a taxon influence index to assess the influence of each taxon on the phylogeny. We found that although most taxa have a weak influence on the phylogeny, a small fraction of influential taxa strongly alter it even in clades only loosely related to them. We conclude that highly influential taxa should be given special attention and sampling them more thoroughly can lead to more dependable phylogenies.