Bacterial microbiota of Kazakhstan cheese revealed by single molecule real time (SMRT) sequencing and its comparison with Belgian, Kalmykian and Italian artisanal cheeses

In Kazakhstan, traditional artisanal cheeses have a long history and are widely consumed. The unique characteristics of local artisanal cheeses are almost completely preserved. However, their microbial communities have rarely been reported. The current study firstly generated the Single Molecule, Real-Time (SMRT) sequencing bacterial diversity profiles of 6 traditional artisanal cheese samples of Kazakhstan origin, followed by comparatively analyzed the microbiota composition between the current dataset and those from cheeses originated from Belgium, Russian Republic of Kalmykia (Kalmykia) and Italy.

Our results together suggest that the diversity of bacterial communities in different groups is stratified by geographic region. This study does not only provide novel information on the bacterial microbiota of traditional artisanal cheese of Kazakhstan at species level, but also interesting insights into the bacterial diversity of artisanal cheeses of various geographical origins.

Across the Kazakhstan cheese samples, a total of 238 bacterial species belonging to 14 phyla and 140 genera were identified. Lactococcus lactis (28.93%), Lactobacillus helveticus (26.43%), Streptococcus thermophilus (12.18%) and Lactobacillus delbrueckii (12.15%) were the dominant bacterial species for these samples. To further evaluate the cheese bacterial diversity of Kazakhstan cheeses in comparison with those from other geographic origins, 16S rRNA datasets of 36 artisanal cheeses from Belgium, Russian Republic of Kalmykia (Kalmykia) and Italy were retrieved from public databases. The cheese bacterial microbiota communities were largely different across sample origins. By principal coordinate analysis (PCoA) and multivariate analysis of variance (MANOVA), the structure of the Kazakhstan artisanal cheese samples was found to be different from those of the other geographic origins. Furthermore, the redundancy analysis (RDA) identified 16 bacterial OTUs as the key variables responsible for such microbiota structural difference.