Genomic analysis of Caldalkalibacillus thermarum TA2.A1 reveals aerobic alkaliphilic metabolism and evolutionary hallmarks linking alkaliphilic bacteria and plant life.

The aerobic thermoalkaliphile Caldalkalibacillus thermarum strain TA2.A1 is a member of a separate order of alkaliphilic bacteria closely related to the Bacillales order. Efforts to relate the genomic information of this evolutionary ancient organism to environmental adaptation have been thwarted by the inability to construct a complete genome. The existing draft genome is highly fragmented due to repetitive regions, and gaps between and over repetitive regions were unbridgeable. To address this, Oxford Nanopore Technology's MinION allowed us to span these repeats through long reads, with over 6000-fold coverage. This resulted in a single 3.34 Mb circular chromosome. The profile of transporters and central metabolism gives insight into why the organism prefers glutamate over sucrose as carbon source. We propose that the deamination of glutamate allows alkalization of the immediate environment, an excellent example of how an extremophile modulates environmental conditions to suit its own requirements. Curiously, plant-like hallmark electron transfer enzymes and transporters are found throughout the genome, such as a cytochrome b6c1 complex and a CO2-concentrating transporter. In addition, multiple self-splicing group II intron-encoded proteins closely aligning to those of a telomerase reverse transcriptase in Arabidopsis thaliana were revealed. Collectively, these features suggest an evolutionary relationship to plant life.

Improving Industrially Relevant Phenotypic Traits by Engineering Chromosome Copy Number in Saccharomyces pastorianus.

The lager-brewing yeast Saccharomyces pastorianus is a hybrid between S. cerevisiae and S. eubayanus with an exceptional degree of aneuploidy. While chromosome copy number variation (CCNV) is present in many industrial Saccharomyces strains and has been linked to various industrially-relevant traits, its impact on the brewing performance of S. pastorianus remains elusive. Here we attempt to delete single copies of chromosomes which are relevant for the production of off-flavor compound diacetyl by centromere silencing. However, the engineered strains display CNV of multiple non-targeted chromosomes. We attribute this unintended CCNV to inherent instability and to a mutagenic effect of electroporation and of centromere-silencing. Regardless, the resulting strains displayed large phenotypic diversity. By growing centromere-silenced cells in repeated sequential batches in medium containing 10% ethanol, mutants with increased ethanol tolerance were obtained. By using CCNV mutagenesis by exposure to the mitotic inhibitor MBC, selection in the same set-up yielded even more tolerant mutants that would not classify as genetically modified organisms. These results show that CCNV of alloaneuploid S. pastorianus genomes is highly unstable, and that CCNV mutagenesis can generate broad diversity. Coupled to effective selection or screening, CCNV mutagenesis presents a potent tool for strain improvement.

The complete genome sequence of the nitrile biocatalyst Rhodocccus rhodochrous ATCC BAA-870.

BACKGROUND: Rhodococci are industrially important soil-dwelling Gram-positive bacteria that are well known for both nitrile hydrolysis and oxidative metabolism of aromatics. Rhodococcus rhodochrous ATCC BAA-870 is capable of metabolising a wide range of aliphatic and aromatic nitriles and amides. The genome of the organism was sequenced and analysed in order to better understand this whole cell biocatalyst. RESULTS: The genome of R. rhodochrous ATCC BAA-870 is the first Rhodococcus genome fully sequenced using Nanopore sequencing. The circular genome contains 5.9 megabase pairs (Mbp) and includes a 0.53 Mbp linear plasmid, that together encode 7548 predicted protein sequences according to BASys annotation, and 5535 predicted protein sequences according to RAST annotation. The genome contains numerous oxidoreductases, 15 identified antibiotic and secondary metabolite gene clusters, several terpene and nonribosomal peptide synthetase clusters, as well as 6 putative clusters of unknown type. The 0.53 Mbp plasmid encodes 677 predicted genes and contains the nitrile converting gene cluster, including a nitrilase, a low molecular weight nitrile hydratase, and an enantioselective amidase. Although there are fewer biotechnologically relevant enzymes compared to those found in rhodococci with larger genomes, such as the well-known Rhodococcus jostii RHA1, the abundance of transporters in combination with the myriad of enzymes found in strain BAA-870 might make it more suitable for use in industrially relevant processes than other rhodococci. CONCLUSIONS: The sequence and comprehensive description of the R. rhodochrous ATCC BAA-870 genome will facilitate the additional exploitation of rhodococci for biotechnological applications, as well as enable further characterisation of this model organism. The genome encodes a wide range of enzymes, many with unknown substrate specificities supporting potential applications in biotechnology, including nitrilases, nitrile hydratase, monooxygenases, cytochrome P450s, reductases, proteases, lipases, and transaminases.

Chromosome level assembly and comparative genome analysis confirm lager-brewing yeasts originated from a single hybridization.

BACKGROUND: The lager brewing yeast, S. pastorianus, is a hybrid between S. cerevisiae and S. eubayanus with extensive chromosome aneuploidy. S. pastorianus is subdivided into Group 1 and Group 2 strains, where Group 2 strains have higher copy number and a larger degree of heterozygosity for S. cerevisiae chromosomes. As a result, Group 2 strains were hypothesized to have emerged from a hybridization event distinct from Group 1 strains. Current genome assemblies of S. pastorianus strains are incomplete and highly fragmented, limiting our ability to investigate their evolutionary history. RESULTS: To fill this gap, we generated a chromosome-level genome assembly of the S. pastorianus strain CBS 1483 from Oxford Nanopore MinION DNA sequencing data and analysed the newly assembled subtelomeric regions and chromosome heterozygosity. To analyse the evolutionary history of S. pastorianus strains, we developed Alpaca: a method to compute sequence similarity between genomes without assuming linear evolution. Alpaca revealed high similarities between the S. cerevisiae subgenomes of Group 1 and 2 strains, and marked differences from sequenced S. cerevisiae strains. CONCLUSIONS: Our findings suggest that Group 1 and Group 2 strains originated from a single hybridization involving a heterozygous S. cerevisiae strain, followed by different evolutionary trajectories. The clear differences between both groups may originate from a severe population bottleneck caused by the isolation of the first pure cultures. Alpaca provides a computationally inexpensive method to analyse evolutionary relationships while considering non-linear evolution such as horizontal gene transfer and sexual reproduction, providing a complementary viewpoint beyond traditional phylogenetic approaches.

Allele-specific genome editing using CRISPR-Cas9 is associated with loss of heterozygosity in diploid yeast.

Targeted DNA double-strand breaks (DSBs) with CRISPR-Cas9 have revolutionized genetic modification by enabling efficient genome editing in a broad range of eukaryotic systems. Accurate gene editing is possible with near-perfect efficiency in haploid or (predominantly) homozygous genomes. However, genomes exhibiting polyploidy and/or high degrees of heterozygosity are less amenable to genetic modification. Here, we report an up to 99-fold lower gene editing efficiency when editing individual heterozygous loci in the yeast genome. Moreover, Cas9-mediated introduction of a DSB resulted in large scale loss of heterozygosity affecting DNA regions up to 360 kb and up to 1700 heterozygous nucleotides, due to replacement of sequences on the targeted chromosome by corresponding sequences from its non-targeted homolog. The observed patterns of loss of heterozygosity were consistent with homology directed repair. The extent and frequency of loss of heterozygosity represent a novel mutagenic side-effect of Cas9-mediated genome editing, which would have to be taken into account in eukaryotic gene editing. In addition to contributing to the limited genetic amenability of heterozygous yeasts, Cas9-mediated loss of heterozygosity could be particularly deleterious for human gene therapy, as loss of heterozygous functional copies of anti-proliferative and pro-apoptotic genes is a known path to cancer.

Galacturonate Metabolism in Anaerobic Chemostat Enrichment Cultures: Combined Fermentation and Acetogenesis by the Dominant sp. nov. "Candidatus Galacturonibacter soehngenii".

Agricultural residues such as sugar beet pulp and citrus peel are rich in pectin, which contains galacturonic acid as a main monomer. Pectin-rich residues are underexploited as feedstocks for production of bulk chemicals or biofuels. The anaerobic, fermentative conversion of d-galacturonate in anaerobic chemostat enrichment cultures provides valuable information toward valorization of these pectin-rich feedstocks. Replicate anaerobic chemostat enrichments, with d-galacturonate as the sole limiting carbon source and inoculum from cow rumen content and rotting orange peels, yielded stable microbial communities, which were dominated by a novel Lachnospiraceae species, for which the name "Candidatus Galacturonibacter soehngenii" was proposed. Acetate was the dominant catabolic product, with formate and H2 as coproducts. The observed molar ratio of acetate and the combined amounts of H2 and formate deviated significantly from 1, which suggested that some of the hydrogen and CO2 formed during d-galacturonate fermentation was converted into acetate via the Wood-Ljungdahl acetogenesis pathway. Indeed, metagenomic analysis of the enrichment cultures indicated that the genome of "Candidatus G. soehngenii" encoded enzymes of the adapted Entner-Doudoroff pathway for d-galacturonate metabolism as well as enzymes of the Wood-Ljungdahl pathway. The simultaneous operation of these pathways may provide a selective advantage under d-galacturonate-limited conditions by enabling a higher specific ATP production rate and lower residual d-galacturonate concentration than would be possible with a strictly fermentative metabolism of this carbon and energy source.IMPORTANCE This study on d-galacturonate metabolism by open, mixed-culture enrichments under anaerobic, d-galacturonate-limited chemostat conditions shows a stable and efficient fermentation of d-galacturonate into acetate as the dominant organic fermentation product. This fermentation stoichiometry and population analyses provide a valuable baseline for interpretation of the conversion of pectin-rich agricultural feedstocks by mixed microbial cultures. Moreover, the results of this study provide a reference for studies on the microbial metabolism of d-galacturonate under different cultivation regimes.

Low affinity uniporter carrier proteins can increase net substrate uptake rate by reducing efflux.

Many organisms have several similar transporters with different affinities for the same substrate. Typically, high-affinity transporters are expressed when substrate is scarce and low-affinity ones when it is abundant. The benefit of using low instead of high-affinity transporters remains unclear, especially when additional nutrient sensors are present. Here, we investigate two hypotheses. It was previously hypothesized that there is a trade-off between the affinity and the catalytic efficiency of transporters, and we find some but no definitive support for it. Additionally, we propose that for uptake by facilitated diffusion, at saturating substrate concentrations, lowering the affinity enhances the net uptake rate by reducing substrate efflux. As a consequence, there exists an optimal, external-substrate-concentration dependent transporter affinity. A computational model of Saccharomyces cerevisiae glycolysis shows that using the low affinity HXT3 transporter instead of the high affinity HXT6 enhances the steady-state flux by 36%. We tried to test this hypothesis with yeast strains expressing a single glucose transporter modified to have either a high or a low affinity. However, due to the intimate link between glucose perception and metabolism, direct experimental proof for this hypothesis remained inconclusive. Still, our theoretical results provide a novel reason for the presence of low-affinity transport systems.

Nanopore sequencing enables near-complete de novo assembly of Saccharomyces cerevisiae reference strain CEN.PK113-7D.

The haploid Saccharomyces cerevisiae strain CEN.PK113-7D is a popular model system for metabolic engineering and systems biology research. Current genome assemblies are based on short-read sequencing data scaffolded based on homology to strain S288C. However, these assemblies contain large sequence gaps, particularly in subtelomeric regions, and the assumption of perfect homology to S288C for scaffolding introduces bias. In this study, we obtained a near-complete genome assembly of CEN.PK113-7D using only Oxford Nanopore Technology's MinION sequencing platform. Fifteen of the 16 chromosomes, the mitochondrial genome and the 2-µm plasmid are assembled in single contigs and all but one chromosome starts or ends in a telomere repeat. This improved genome assembly contains 770 Kbp of added sequence containing 248 gene annotations in comparison to the previous assembly of CEN.PK113-7D. Many of these genes encode functions determining fitness in specific growth conditions and are therefore highly relevant for various industrial applications. Furthermore, we discovered a translocation between chromosomes III and VIII that caused misidentification of a MAL locus in the previous CEN.PK113-7D assembly. This study demonstrates the power of long-read sequencing by providing a high-quality reference assembly and annotation of CEN.PK113-7D and places a caveat on assumed genome stability of microorganisms.

Growth-rate dependency of de novo resveratrol production in chemostat cultures of an engineered Saccharomyces cerevisiae strain.

INTRODUCTION: Saccharomyces cerevisiae has become a popular host for production of non-native compounds. The metabolic pathways involved generally require a net input of energy. To maximize the ATP yield on sugar in S. cerevisiae, industrial cultivation is typically performed in aerobic, sugar-limited fed-batch reactors which, due to constraints in oxygen transfer and cooling capacities, have to be operated at low specific growth rates. Because intracellular levels of key metabolites are growth-rate dependent, slow growth can significantly affect biomass-specific productivity. Using an engineered Saccharomyces cerevisiae strain expressing a heterologous pathway for resveratrol production as a model energy-requiring product, the impact of specific growth rate on yeast physiology and productivity was investigated in aerobic, glucose-limited chemostat cultures. RESULTS: Stoichiometric analysis revealed that de novo resveratrol production from glucose requires 13 moles of ATP per mole of produced resveratrol. The biomass-specific production rate of resveratrol showed a strong positive correlation with the specific growth rate. At low growth rates a substantial fraction of the carbon source was invested in cellular maintenance-energy requirements (e.g. 27 % at 0.03 h(-1)). This distribution of resources was unaffected by resveratrol production. Formation of the by-products coumaric, phloretic and cinnamic acid had no detectable effect on maintenance energy requirement and yeast physiology in chemostat. Expression of the heterologous pathway led to marked differences in transcript levels in the resveratrol-producing strain, including increased expression levels of genes involved in pathways for precursor supply (e.g. ARO7 and ARO9 involved in phenylalanine biosynthesis). The observed strong differential expression of many glucose-responsive genes in the resveratrol producer as compared to a congenic reference strain could be explained from higher residual glucose concentrations and higher relative growth rates in cultures of the resveratrol producer. CONCLUSIONS: De novo resveratrol production by engineered S. cerevisiae is an energy demanding process. Resveratrol production by an engineered strain exhibited a strong correlation with specific growth rate. Since industrial production in fed-batch reactors typically involves low specific growth rates, this study emphasizes the need for uncoupling growth and product formation via energy-requiring pathways.

Physiological and transcriptional responses of anaerobic chemostat cultures of Saccharomyces cerevisiae subjected to diurnal temperature cycles.

Diurnal temperature cycling is an intrinsic characteristic of many exposed microbial ecosystems. However, its influence on yeast physiology and the yeast transcriptome has not been studied in detail. In this study, 24-h sinusoidal temperature cycles, oscillating between 12°C and 30°C, were imposed on anaerobic, glucose-limited chemostat cultures of Saccharomyces cerevisiae. After three diurnal temperature cycles (DTC), concentrations of glucose and extracellular metabolites as well as CO2 production rates showed regular, reproducible circadian rhythms. DTC also led to waves of transcriptional activation and repression, which involved one-sixth of the yeast genome. A substantial fraction of these DTC-responsive genes appeared to respond primarily to changes in the glucose concentration. Elimination of known glucose-responsive genes revealed an overrepresentation of previously identified temperature-responsive genes as well as genes involved in the cell cycle and de novo purine biosynthesis. In-depth analysis demonstrated that DTC led to a partial synchronization of the cell cycle of the yeast populations in chemostat cultures, which was lost upon release from DTC. Comparison of DTC results with data from steady-state cultures showed that the 24-h DTC was sufficiently slow to allow S. cerevisiae chemostat cultures to acclimate their transcriptome and physiology at the DTC temperature maximum and to approach acclimation at the DTC temperature minimum. Furthermore, this comparison and literature data on growth rate-dependent cell cycle phase distribution indicated that cell cycle synchronization was most likely an effect of imposed fluctuations of the relative growth rate (µ/µmax) rather than a direct effect of temperature.