Quantification of microbial species in solid state fermentation samples using signature genomic sequences

Solid state fermentation processes are mediated by the collective metabolism of specialized microbial communities. Monitoring the relative abundance of dominating species is a critical task in quality control, which is traditionally done by wet lab techniques, such as quantitative PCR (qPCR). In this study, we developed a computational method to quantify microbial species in metagenomes based on their signature genomic sequences, i.e., unique k-mers. Bacterial species found in fermentation starters of a Chinese liquor producer were used as examples to demonstrate the development and application of the method. A database was constructed, comprising 562 complete genome sequences of 93 bacterial species that had been found in relevant fermentation samples. K-mers in length of 12 were extracted from each species and compared against each other to identify the ones that were unique to each species. The quantity of a species was determined by the average frequencies of unique k-mers encountered in the metagenome. Six dominating bacterial species were chosen as reporter species to test the quantification method. Four metagenome datasets were simulated, which contained various portions of sequence reads generated from the genomes of the reporter species. The amount of reads sampled from each reporter species followed a pre-determined ratio, i.e., a known relationship in relative abundance. For each simulated dataset, the cell number of each reporter species was computed based on the unique k-mers found in the metagenome. In all datasets, the computed quantities of the reporter species reflected the expected relative abundance by displaying a linear relationship with the pre-determined ratio. This demonstrates that quantification based on a set of unique k-mers is a reliable way to detect relative abundance among species. Besides industrial fermentation, this method may also be applied to areas such as wastewater treatment, microbiota analysis, etc.

[1]  R Zhang,et al.  Aroma characteristics of Moutai‐flavour liquor produced with Bacillus licheniformis by solid‐state fermentation , 2013, Letters in applied microbiology.

[2]  Patrick D. Schloss,et al.  Reducing the Effects of PCR Amplification and Sequencing Artifacts on 16S rRNA-Based Studies , 2011, PloS one.

[3]  B. Han,et al.  Baijiu (白酒), Chinese liquor: History, classification and manufacture , 2016 .

[4]  Guangyuan Jin,et al.  Mystery behind Chinese liquor fermentation , 2017 .

[5]  Patrick D. Schloss,et al.  The Effects of Alignment Quality, Distance Calculation Method, Sequence Filtering, and Region on the Analysis of 16S rRNA Gene-Based Studies , 2010, PLoS Comput. Biol..

[6]  Zhaohui Xu,et al.  Employment of Near Full-Length Ribosome Gene TA-Cloning and Primer-Blast to Detect Multiple Species in a Natural Complex Microbial Community Using Species-Specific Primers Designed with Their Genome Sequences , 2016, Molecular Biotechnology.

[7]  Bernhard Y. Renard,et al.  Metagenomic abundance estimation and diagnostic testing on species level , 2012, Nucleic acids research.

[8]  Guo-li Gong,et al.  Microorganisms in Daqu: a starter culture of Chinese Maotai-flavor liquor , 2008 .

[9]  Liu Xiu,et al.  Determination of microbial diversity in Daqu, a fermentation starter culture of Maotai liquor, using nested PCR-denaturing gradient gel electrophoresis , 2012, World Journal of Microbiology and Biotechnology.

[10]  Jian-Gang Yang,et al.  Comparison of microbial communities in the fermentation starter used to brew Xiaoqu liquor , 2017 .

[11]  Xiaofei Ding,et al.  Characterisation of microbial communities in Chinese liquor fermentation starters Daqu using nested PCR-DGGE , 2014, World Journal of Microbiology and Biotechnology.

[12]  Jin-Ho Seo,et al.  A competitive quantitative polymerase chain reaction method for characterizing the population dynamics during kimchi fermentation , 2014, Journal of Industrial Microbiology & Biotechnology.

[13]  Marcel H Zwietering,et al.  Complex microbiota of a Chinese "Fen" liquor fermentation starter (Fen-Daqu), revealed by culture-dependent and culture-independent methods. , 2012, Food microbiology.

[14]  Jennifer M. Fettweis,et al.  The truth about metagenomics: quantifying and counteracting bias in 16S rRNA studies , 2015, BMC Microbiology.

[15]  Derrick E. Wood,et al.  Kraken: ultrafast metagenomic sequence classification using exact alignments , 2014, Genome Biology.

[16]  Yue-Qin Tang,et al.  Characterization of the microbial community in three types of fermentation starters used for Chinese liquor production , 2015 .

[17]  C. Huttenhower,et al.  Metagenomic microbial community profiling using unique clade-specific marker genes , 2012, Nature Methods.

[18]  Liangqiang Chen,et al.  Improving flavor metabolism of Saccharomyces cerevisiae by mixed culture with Bacillus licheniformis for Chinese Maotai-flavor liquor making , 2015, Journal of Industrial Microbiology & Biotechnology.

[19]  Lauren M. Bragg,et al.  Metagenomics using next-generation sequencing. , 2014, Methods in molecular biology.

[20]  Hongqiang Ren,et al.  Metagenomic analysis of bacterial community composition and antibiotic resistance genes in a wastewater treatment plant and its receiving surface water. , 2016, Ecotoxicology and environmental safety.

[21]  K. Pollard,et al.  Toward Accurate and Quantitative Comparative Metagenomics , 2016, Cell.

[22]  Chittibabu Guda,et al.  MetaID: A novel method for identification and quantification of metagenomic samples , 2013, BMC Genomics.