Genomic evolution of Saccharomyces cerevisiae under Chinese rice wine fermentation

Analysis of Key Genes Responsible for Low Urea Production in Saccharomyces cerevisiae JH301

There is a potential safety risk with ethyl carbamate (EC) in Hongqu Huangjiu production; 90% of the EC in rice wine is produced by the reaction of the urea with the alcohol of Saccharomyces cerevisiae. In our previous experiments, we screened and obtained a S. cerevisiae strain JH301 that offered low urea production. However, the key genes responsible for low urea production of strain JH301 remain unclear. Here, the whole genome sequencing of S. cerevisiae strain JH301 was accomplished via a next-generation high-throughput sequencing and long-read sequencing technology. There are six main pathways related to the urea metabolism of strain JH301 based on KEGG pathway mapping. Three species-specific genes are related to the urea metabolism pathways and were found in comparative genome analysis between strains JH301 and S288c during Hongqu Huangjiu production for the first time. Finally, the ARG80 gene was found to be likely a key gene responsible for low urea production of S. cerevisiae strain JH301, as determined by PCR and qRT-PCR check analyses from DNA and RNA levers. In conclusion, the results are useful for a scientific understanding of the mechanism of low urea production by Saccharomyces cerevisiae during Hongqu Huangjiu fermentation. It also is important to control the urea and EC contents in Hongqu Huangjiu production.

Bio-Organic Fertilizer: A Green Technology to Reduce Synthetic N and P Fertilizer for Rice Production

Decomposed organic materials, in combination with plant growth-promoting bacteria (PGPB), are environmentally friendly and reduce synthetic fertilizer use in rice production. A bio-organic fertilizer (BoF) was prepared using kitchen waste (79%), chita-dhan (unfilled rice grain) biochar (15%), rock phosphate (5%), and a consortium of 10 PGPB (1%) to supplement 30% nitrogen and to replace triple superphosphate (TSP) fertilizer in rice production with an improvement of soil health. PGPB were local isolates and identified using 16S ribosomal RNA partial gene sequences as Bacillus mycoides, Proteus sp., Bacillus cereus, Bacillus subtilis, Bacillus pumilus, Paenibacillus polymyxa, and Paenibacillus spp. Isolates could fix N₂ by 0.7-1.4 g kg^-1, solubilize 0.1-1.2 g kg^-1 phosphate, and produce 0.1-40 g kg^-1 indoleacetic acid. The performance of BoF was evaluated by 16 field experiments and 18 farmers' field demonstration trials during the year 2017-2020 in different parts of Bangladesh. Performances of BoF were evaluated based on control (T₁), full synthetic fertilizer dose of N, P, and K (T₂), BoF (2 t ha^-1) + 70% N as urea + 100% K as muriate of potash (T₃), 70% N as urea + 100% P as TSP + 100% K as muriate of potash (T₄), and 2 t ha^-1 BoF (T₅) treatments. At the research station, average grain yield improved by 10-13% in T₃ compared with T₂ treatment. Depending on seasons, higher agronomic N use efficiency (19-30%), physiological N use efficiency (8-18%), partial factor productivity (PFP)_N (114-150%), recovery efficiency (RE)_N (3-31%), N harvest index (HI_N) (14-24%), agronomic P use efficiency (22-25%), partial factor productivity of P (9-12%), ARE_P (15-23%), and HI_P (3-6%) were obtained in T₃ compared with T₂ treatment. Research results were reflected in farmers' field, and significant (P < 0.05) higher plant height, tiller, panicle, grain yield, partial factor productivity of N and P were obtained in the same treatment. Application of BoF improved soil organic carbon by 6-13%, along with an increased number of PGPB as compared with full synthetic fertilizer dose. In conclusion, tested BoF can be considered as a green technology to reduce 30% synthetic N and 100% TSP requirements in rice production with improved soil health.

Chinese Yellow Rice Wine Processing with Reduced Ethyl Carbamate Formation by Deleting Transcriptional Regulator Dal80p in Saccharomyces cerevisiae

Ethyl carbamate (EC) is a potential carcinogen that forms spontaneously during Chinese rice wine fermentation. The primary precursor for EC formation is urea, which originates from both external sources and arginine degradation. Urea degradation is suppressed by nitrogen catabolite repression (NCR) in Saccharomyces cerevisiae. The regulation of NCR is mediated by two positive regulators (Gln3p, Gat1p/Nil1p) and two negative regulators (Dal80p/Uga43p, Deh1p/Nil2p/GZF3p). DAL80 revealed higher transcriptional level when yeast cells were cultivated under nitrogen-limited conditions. In this study, when DAL80-deleted yeast cells were compared to wild-type BY4741 cells, less urea was accumulated, and genes involved in urea utilization were up-regulated. Furthermore, Chinese rice wine fermentation was conducted using dal80Δ cells; the concentrations of urea and EC were both reduced when compared to the BY4741 and traditional fermentation starter. The findings of this work indicated Dal80p is involved in EC formation possibly through regulating urea metabolism and may be used as the potential target for EC reduction.

Isolation of biogenic amine-negative lactic acid bacteria for Chinese rice wine fermentation based on molecular marker reverse screening

BACKGROUND

Chinese rice wine (CRW; a traditional alcoholic beverage in China with unique flavor and high nutritional value) containing high level of biogenic amines (BAs) may be deleterious to human health. The processes of rice soaking, primary fermentation and secondary fermentation were found to be the major factors for accumulation of BAs during industrial CRW production.

RESULTS

To reduce the risk of the formation of BAs in CRW production, Enterococcus durans AR315, a BA-negative lactic acid bacterium, was isolated from CRW samples by PCR-based molecular marker reverse screening in this work. With addition of AR315 during steeping rice phase, the level of total BAs was significantly decreased by 45.1% in comparison with the control. Moreover, the final BA concentration with the addition of AR315 was 27.6% lower than that of the control during fermentation phase.

CONCLUSIONS

To our knowledge, this is the first report of decreased accumulation of BAs in CRW production using a BA-negative lactic acid bacterium. Hence, using a BA-negative lactic acid bacterium as a starter culture could be an efficient strategy for significantly reducing the formation of BAs, which has the potential for industrial application in CRW production. © 2020 Society of Chemical Industry.

Adaptation to Industrial Stressors Through Genomic and Transcriptional Plasticity in a Bioethanol Producing Fission Yeast Isolate

Schizosaccharomyces pombe is a model unicellular eukaryote with ties to the basic research, oenology and industrial biotechnology sectors. While most investigations into S. pombe cell biology utilize Leupold’s 972h^- laboratory strain background, recent studies have described a wealth of genetic and phenotypic diversity within wild populations of S. pombe including stress resistance phenotypes which may be of interest to industry. Here we describe the genomic and transcriptomic characterization of Wilmar-P, an S. pombe isolate used for bioethanol production from sugarcane molasses at industrial scale. Novel sequences present in Wilmar-P but not in the laboratory S. pombe genome included multiple coding sequences with near-perfect nucleotide identity to Schizosaccharomyces octosporus sequences. Wilmar-P also contained a ∼100kb duplication in the right arm of chromosome III, a region harboring ght5⁺, the predominant hexose transporter encoding gene. Transcriptomic analysis of Wilmar-P grown in molasses revealed strong downregulation of core environmental stress response genes and upregulation of hexose transporters and drug efflux pumps compared to laboratory S. pombe. Finally, examination of the regulatory network of Scr1, which is involved in the regulation of several genes differentially expressed on molasses, revealed expanded binding of this transcription factor in Wilmar-P compared to laboratory S. pombe in the molasses condition. Together our results point to both genomic plasticity and transcriptomic adaptation as mechanisms driving phenotypic adaptation of Wilmar-P to the molasses environment and therefore adds to our understanding of genetic diversity within industrial fission yeast strains and the capacity of this strain for commercial scale bioethanol production.

Effect of milk types on the attributes of a glutinous rice wine-fermented yogurt-like product

This study was conducted to analyze the effect of milk types on the attributes of the glutinous rice wine-fermented yogurt-like product named Kouwan Lao (KWL). Four types of raw milks were used in this study, including high temperature, long time (HTLT: H milk), HTLT milk supplemented with 3% skim milk powder (S milk), pasteurized milk (P milk), and ultra-high temperature milk (U milk). Microbiological compositions of the fermented glutinous rice and KWL at different stages were analyzed using PCR-denaturing gradient gel electrophoresis and gene sequencing based on 16S rRNA and 26S rRNA. The physicochemical properties of KWL samples were determined, and textural properties of those were analyzed using a texture analyzer (Jiawei Innovation and Technology Co. Ltd., Zhejiang Province, China). The microstructure of KWL samples was observed using scanning electron microscopy. The results showed that the milk types had significant influences on the bacterial composition of KWL. In the curdling process, the predominant bacteria of H, S, P, and U KWL samples were Lactobacillus brevis, Janthinobacterium sp., Lactobacillus casei, and Streptococcus agalactiae, respectively. In the ripening process, the main strains in H KWL were Enterococcus faecium and Pediococcus pentosaceus. Lactobacillus casei and Lactobacillus paracasei were the dominant bacteria of U KWL. Lactobacillus casei was the main strain of P KWL, and no bacteria were detected in S KWL. Saccharomyces cerevisiae was the dominant fungus of KWL, and no significant effect of milk types on fungal composition of KWL was found. The results of physicochemical properties showed that significant differences in protein contents were found in the KWL samples, and P KWL had the highest protein content. The fat content of U KWL was significantly higher than that of samples from the other 3 groups. The U KWL and P KWL showed lower moisture than that of the other 2 KWL samples. In addition, no significant difference in pH value was found in all samples. The results of texture analysis and microstructure showed that compared with other 3 types of KWL samples, the related mass parameters of U KWL were more advantageous and improved significantly with the increase of the heat treatment temperature of raw milk and the addition of skim milk powder. Our findings revealed the effects of milk types on microbial composition, physicochemical properties, textural properties, and microstructure of KWL, and provided a basic theory for the optimization and industrial production of KWL.

Comprehensive evaluation of structural variation detection algorithms for whole genome sequencing

BACKGROUND

Structural variations (SVs) or copy number variations (CNVs) greatly impact the functions of the genes encoded in the genome and are responsible for diverse human diseases. Although a number of existing SV detection algorithms can detect many types of SVs using whole genome sequencing (WGS) data, no single algorithm can call every type of SVs with high precision and high recall.

RESULTS

We comprehensively evaluate the performance of 69 existing SV detection algorithms using multiple simulated and real WGS datasets. The results highlight a subset of algorithms that accurately call SVs depending on specific types and size ranges of the SVs and that accurately determine breakpoints, sizes, and genotypes of the SVs. We enumerate potential good algorithms for each SV category, among which GRIDSS, Lumpy, SVseq2, SoftSV, Manta, and Wham are better algorithms in deletion or duplication categories. To improve the accuracy of SV calling, we systematically evaluate the accuracy of overlapping calls between possible combinations of algorithms for every type and size range of SVs. The results demonstrate that both the precision and recall for overlapping calls vary depending on the combinations of specific algorithms rather than the combinations of methods used in the algorithms.

CONCLUSION

These results suggest that careful selection of the algorithms for each type and size range of SVs is required for accurate calling of SVs. The selection of specific pairs of algorithms for overlapping calls promises to effectively improve the SV detection accuracy.

Chromosomal genome assembly of the ethanol production strain CBS 11270 indicates a highly dynamic genome structure in the yeast species Brettanomyces bruxellensis

Here, we present the genome of the industrial ethanol production strain Brettanomyces bruxellensis CBS 11270. The nuclear genome was found to be diploid, containing four chromosomes with sizes of ranging from 2.2 to 4.0 Mbp. A 75 Kbp mitochondrial genome was also identified. Comparing the homologous chromosomes, we detected that 0.32% of nucleotides were polymorphic, i.e. formed single nucleotide polymorphisms (SNPs), 40.6% of them were found in coding regions (i.e. 0.13% of all nucleotides formed SNPs and were in coding regions). In addition, 8,538 indels were found. The total number of protein coding genes was 4897, of them, 4,284 were annotated on chromosomes; and the mitochondrial genome contained 18 protein coding genes. Additionally, 595 genes, which were annotated, were on contigs not associated with chromosomes. A number of genes was duplicated, most of them as tandem repeats, including a six-gene cluster located on chromosome 3. There were also examples of interchromosomal gene duplications, including a duplication of a six-gene cluster, which was found on both chromosomes 1 and 4. Gene copy number analysis suggested loss of heterozygosity for 372 genes. This may reflect adaptation to relatively harsh but constant conditions of continuous fermentation. Analysis of gene topology showed that most of these losses occurred in clusters of more than one gene, the largest cluster comprising 33 genes. Comparative analysis against the wine isolate CBS 2499 revealed 88,534 SNPs and 8,133 indels. Moreover, when the scaffolds of the CBS 2499 genome assembly were aligned against the chromosomes of CBS 11270, many of them aligned completely, some have chunks aligned to different chromosomes, and some were in fact rearranged. Our findings indicate a highly dynamic genome within the species B. bruxellensis and a tendency towards reduction of gene number in long-term continuous cultivation.

Complete genome sequence and analysis of the industrial Saccharomyces cerevisiae strain N85 used in Chinese rice wine production

Chinese rice wine is a popular traditional alcoholic beverage in China, while its brewing processes have rarely been explored. We herein report the first gapless, near-finished genome sequence of the yeast strain Saccharomyces cerevisiae N85 for Chinese rice wine production. Several assembly methods were used to integrate Pacific Bioscience (PacBio) and Illumina sequencing data to achieve high-quality genome sequencing of the strain. The genome encodes more than 6,000 predicted proteins, and 238 long non-coding RNAs, which are validated by RNA-sequencing data. Moreover, our annotation predicts 171 novel genes that are not present in the reference S288c genome. We also identified 65,902 single nucleotide polymorphisms and small indels, many of which are located within genic regions. Dozens of larger copy-number variations and translocations were detected, mainly enriched in the subtelomeres, suggesting these regions may be related to genomic evolution. This study will serve as a milestone in studying of Chinese rice wine and related beverages in China and in other countries. It will help to develop more scientific and modern fermentation processes of Chinese rice wine, and explore metabolism pathways of desired and harmful components in Chinese rice wine to improve its taste and nutritional value.

Whole-Genome Sequence and Variant Analysis of W303, a Widely-Used Strain of Saccharomyces cerevisiae

The yeast Saccharomyces cerevisiae has emerged as a superior model organism. Selection of distinct laboratory strains of S. cerevisiae with unique phenotypic properties, such as superior mating or sporulation efficiencies, has facilitated advancements in research. W303 is one such laboratory strain that is closely related to the first completely sequenced yeast strain, S288C. In this work, we provide a high-quality, annotated genome sequence for W303 for utilization in comparative analyses and genome-wide studies. Approximately 9500 variations exist between S288C and W303, affecting the protein sequences of ∼700 genes. A listing of the polymorphisms and divergent genes is provided for researchers interested in identifying the genetic basis for phenotypic differences between W303 and S288C. Several divergent functional gene families were identified, including flocculation and sporulation genes, likely representing selection for desirable laboratory phenotypes. Interestingly, remnants of ancestor wine strains were found on several chromosomes. Finally, as a test of the utility of the high-quality reference genome, variant mapping revealed more accurate identification of accumulated mutations in passaged mismatch repair-defective strains.

Elucidating and Regulating the Acetoin Production Role of Microbial Functional Groups in Multispecies Acetic Acid Fermentation

Acetoin (3-hydroxy-2-butanone) formation in vinegar microbiota is crucial for the flavor quality of Zhenjiang aromatic vinegar, a traditional vinegar produced from cereals. However, the specific microorganisms responsible for acetoin formation in this centuries-long repeated batch fermentation have not yet been clearly identified. Here, the microbial distribution discrepancy in the diacetyl/acetoin metabolic pathway of vinegar microbiota was revealed at the species level by a combination of metagenomic sequencing and clone library analysis. The results showed that Acetobacter pasteurianus and 4 Lactobacillus species (Lactobacillus buchneri, Lactobacillus reuteri, Lactobacillus fermentum, and Lactobacillus brevis) might be functional producers of acetoin from 2-acetolactate in vinegar microbiota. Furthermore, A. pasteurianus G3-2, L. brevis 4-22, L. fermentum M10-3, and L. buchneri F2-5 were isolated from vinegar microbiota by a culture-dependent method. The acetoin concentrations in two cocultures (L. brevis 4-22 plus A. pasteurianus G3-2 and L. fermentum M10-3 plus A. pasteurianus G3-2) were obviously higher than those in monocultures of lactic acid bacteria (LAB), while L. buchneri F2-5 did not produce more acetoin when coinoculated with A. pasteurianus G3-2. Last, the acetoin-producing function of vinegar microbiota was regulated in situ via augmentation with functional species in vinegar Pei After 72 h of fermentation, augmentation with A. pasteurianus G3-2 plus L. brevis 4-22, L. fermentum M10-3, or L. buchneri F2-5 significantly increased the acetoin content in vinegar Pei compared with the control group. This study provides a perspective on elucidating and manipulating different metabolic roles of microbes during flavor formation in vinegar microbiota.

IMPORTANCE

Acetoin (3-hydroxy-2-butanone) formation in vinegar microbiota is crucial for the flavor quality of Zhenjiang aromatic vinegar, a traditional vinegar produced from cereals. Thus, it is of interest to understand which microbes are driving the formation of acetoin to elucidate the microbial distribution discrepancy in the acetoin metabolic pathway and to regulate the metabolic function of functional microbial groups in vinegar microbiota. Our study provides a perspective on elucidating and manipulating different metabolic roles of microbes during flavor formation in vinegar microbiota.

Integration of new alternative reference strain genome sequences into the Saccharomyces genome database

The Saccharomyces Genome Database (SGD; http://www.yeastgenome.org/) is the authoritative community resource for the Saccharomyces cerevisiae reference genome sequence and its annotation. To provide a wider scope of genetic and phenotypic variation in yeast, the genome sequences and their corresponding annotations from 11 alternative S. cerevisiae reference strains have been integrated into SGD. Genomic and protein sequence information for genes from these strains are now available on the Sequence and Protein tab of the corresponding Locus Summary pages. We illustrate how these genome sequences can be utilized to aid our understanding of strain-specific functional and phenotypic differences.

Database URL:www.yeastgenome.org

AGAPE (Automated Genome Analysis PipelinE) for pan-genome analysis of Saccharomyces cerevisiae

The characterization and public release of genome sequences from thousands of organisms is expanding the scope for genetic variation studies. However, understanding the phenotypic consequences of genetic variation remains a challenge in eukaryotes due to the complexity of the genotype-phenotype map. One approach to this is the intensive study of model systems for which diverse sources of information can be accumulated and integrated. Saccharomyces cerevisiae is an extensively studied model organism, with well-known protein functions and thoroughly curated phenotype data. To develop and expand the available resources linking genomic variation with function in yeast, we aim to model the pan-genome of S. cerevisiae. To initiate the yeast pan-genome, we newly sequenced or re-sequenced the genomes of 25 strains that are commonly used in the yeast research community using advanced sequencing technology at high quality. We also developed a pipeline for automated pan-genome analysis, which integrates the steps of assembly, annotation, and variation calling. To assign strain-specific functional annotations, we identified genes that were not present in the reference genome. We classified these according to their presence or absence across strains and characterized each group of genes with known functional and phenotypic features. The functional roles of novel genes not found in the reference genome and associated with strains or groups of strains appear to be consistent with anticipated adaptations in specific lineages. As more S. cerevisiae strain genomes are released, our analysis can be used to collate genome data and relate it to lineage-specific patterns of genome evolution. Our new tool set will enhance our understanding of genomic and functional evolution in S. cerevisiae, and will be available to the yeast genetics and molecular biology community.

Making the difference: integrating structural variation detection tools

From prokaryotes to eukaryotes, phenotypic variation, adaptation and speciation has been associated with structural variation between genomes of individuals within the same species. Many computer algorithms detecting such variations (callers) have recently been developed, spurred by the advent of the next-generation sequencing technology. Such callers mainly exploit split-read mapping or paired-end read mapping. However, as different callers are geared towards different types of structural variation, there is still no single caller that can be considered a community standard; instead, increasingly the various callers are combined in integrated pipelines. In this article, we review a wide range of callers, discuss challenges in the integration step and present a survey of pipelines used in population genomics studies. Based on our findings, we provide general recommendations on how to set-up such pipelines. Finally, we present an outlook on future challenges in structural variation detection.