Bioinformatic tools for marker gene sequencing data analysis are continuously and rapidly evolving, thus integrating most recent techniques and tools is challenging. We present an R package for data analysis of 16S and ITS amplicons based sequencing. This workflow is based on several R functions and performs automatic treatments from fastq sequence files to diversity and differential analysis with statistical validation. The main purpose of this package is to automate bioinformatic analysis, ensure reproducibility between projects, and to be flexible enough to quickly integrate new bioinformatic tools or statistical methods. rANOMALY is an easy to install and customizable R package, that uses amplicon sequence variants (ASV) level for microbial community characterization. It integrates all assets of the latest bioinformatics methods, such as better sequence tracking, decontamination from control samples, use of multiple reference databases for taxonomic annotation, all main ecological analysis for which we propose advanced statistical tests, and a cross-validated differential analysis by four different methods. Our package produces ready to publish figures, and all of its outputs are made to be integrated in Rmarkdown code to produce automated reports.

Streptococcus thermophilus is of major importance for cheese manufacturing to ensure rapid acidification; however, studies indicate that intensive use of commercial strains leads to the loss of typical characteristics of the products. To strengthen the link between the product and its geographical area and improve the sensory qualities of cheeses, cheese-producing protected designations of origin (PDO) are increasingly interested in the development of specific autochthonous starter cultures. The present study is therefore investigating the genetic and functional diversity of S. thermophilus strains isolated from a local cheese-producing PDO area. Putative S. thermophilus isolates were isolated and identified from milk collected in the Saint-Nectaire cheese-producing PDO area and from commercial starters. Whole genomes of isolates were sequenced, and a comparative analysis based on their pan-genome was carried out. Important functional properties were studied, including acidifying and proteolytic activities. Twenty-two isolates representative of the diversity of the geographical area and four commercial strains were selected for comparison. The resulting phylogenetic trees do not correspond to the geographical distribution of isolates. The clustering based on the pan-genome analysis indicates that isolates are divided into five distinct groups. A Kyoto Encyclopedia of Genes and Genomes (KEGG) functional annotation of the accessory genes indicates that the accessory gene contents of isolates are involved in different functional categories. High variability in acidifying activities and less diversity in proteolytic activities were also observed. These results indicate that high genetic and functional variabilities of the species S. thermophilus may arise from a small (1,800 km²) geographical area and may be exploited to meet demand for use as autochthonous starters.

In dairy cattle, teat disinfection at the end of milking is commonly applied to limit colonization of the milk by pathogenic microorganisms via the teat canal. The post-milking products used can irritate the teat skin and unbalance its microbial population. Our study aimed to assess the impact of different milking products on the balance of the microbial communities on the teat skin of cows and in their milk. For 12 weeks at the end of each milking operation, three groups of seven Holstein dairy cows on pasture received either a chlorhexidine gluconate-based product (G) or a hydrocolloidal water-in-oil emulsion (A), or no post-milking product (C). The composition of the bacterial and fungal communities on the teat skin and in the milk were characterized using a culture-dependent method and by high-throughput sequencing of marker genes to obtain amplicon sequence variants (ASVs). The individual microbiota on the cows’ teat skin was compared for the first time to that of a cow pool. In contrast to the milk, the post-milking treatment influenced the microbiota of the teat skin, which revealed a high microbial diversity. The water-in-oil emulsion appeared to slightly favour lactic acid bacteria and yeasts and to limit the development of undesirable bacteria such as Pseudomonas and Staphylococcus.

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