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.

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.

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.

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.

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.

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.

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 H2 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 H2 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 H2 was obtained, which is 3.64 times higher than the concentration in unoptimized conditions (13.31 mmol H2/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 H2 (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.

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.

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.

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.

FROGS: Find, Rapidly, OTUs with Galaxy Solution.

Motivation: Metagenomics leads to major advances in microbial ecology and biologists need user friendly tools to analyze their data on their own. Results: This Galaxy-supported pipeline, called FROGS, is designed to analyze large sets of amplicon sequences and produce abundance tables of Operational Taxonomic Units (OTUs) and their taxonomic affiliation. The clustering uses Swarm. The chimera removal uses VSEARCH, combined with original cross-sample validation. The taxonomic affiliation returns an innovative multi-affiliation output to highlight databases conflicts and uncertainties. Statistical results and numerous graphical illustrations are produced along the way to monitor the pipeline. FROGS was tested for the detection and quantification of OTUs on real and in silico datasets and proved to be rapid, robust and highly sensitive. It compares favorably with the widespread mothur, UPARSE and QIIME. Availability and implementation: Source code and instructions for installation: https://github.com/geraldinepascal/FROGS.git. A companion website: http://frogs.toulouse.inra.fr. Contact: geraldine.pascal@inra.fr. Supplementary information: Supplementary data are available at Bioinformatics online.

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.

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.

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.

Microbial ecosystems are dominated by specialist taxa.

Abundance and specificity are two key characteristics of species distribution and biodiversity. Theories of species assembly aim to reproduce the empirical joint patterns of specificity and abundance, with the goal to explain patterns of biodiversity across habitats. The specialist-generalist paradigm predicts that specialists should have a local advantage over generalists and thus be more abundant. We developed a specificity index to analyse abundance-specificity relationships in microbial ecosystems. By analysing microbiota spanning 23 habitats from three very different data sets covering a wide range of sequencing depths and environmental conditions, we find that habitats are consistently dominated by specialist taxa, resulting in a strong, positive correlation between abundance and specificity. This finding is consistent over several levels of taxonomic aggregation and robust to errors in abundance measures. The relationship explains why shallow sequencing captures similar ß-diversity as deep sequencing, and can be sufficient to capture the habitat-specific functions of microbial communities.

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.

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.

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.