Genome-wide CRISPRi screen identifies enhanced autolithotrophic phenotypes in acetogenic bacterium Eubacterium limosum.

Acetogenic bacteria are a unique biocatalyst that highly promises to develop the sustainable bioconversion of carbon oxides (e.g., CO and CO2) into multicarbon biochemicals. Genotype-phenotype relationships are important for engineering their metabolic capability to enhance their biocatalytic performance; however, systemic investigation on the fitness contribution of individual gene has been limited. Here, we report genome-scale CRISPR interference screening using 41,939 guide RNAs designed from the E. limosum genome, one of the model acetogenic species, where all genes were targeted for transcriptional suppression. We investigated the fitness contributions of 96% of the total genes identified, revealing the gene fitness and essentiality for heterotrophic and autotrophic metabolisms. Our data show that the Wood-Ljungdahl pathway, membrane regeneration, membrane protein biosynthesis, and butyrate synthesis are essential for autotrophic acetogenesis in E. limosum. Furthermore, we discovered genes that are repression targets that unbiasedly increased autotrophic growth rates fourfold and acetoin production 1.5-fold compared to the wild-type strain under CO2-H2 conditions. These results provide insight for understanding acetogenic metabolism and genome engineering in acetogenic bacteria.

Serial adaptive laboratory evolution enhances mixed carbon metabolic capacity of Escherichia coli.

Microbes have inherent capacities for utilizing various carbon sources, however they often exhibit sub-par fitness due to low metabolic efficiency. To test whether a bacterial strain can optimally utilize multiple carbon sources, Escherichia coli was serially evolved in L-lactate and glycerol. This yielded two end-point strains that evolved first in L-lactate then in glycerol, and vice versa. The end-point strains displayed a universal growth advantage on single and a mixture of adaptive carbon sources, enabled by a concerted action of carbon source-specialists and generalist mutants. The combination of just four variants of glpK, ppsA, ydcI, and rph-pyrE, accounted for more than 80% of end-point strain fitness. In addition, machine learning analysis revealed a coordinated activity of transcriptional regulators imparting condition-specific regulation of gene expression. The effectiveness of the serial adaptive laboratory evolution (ALE) scheme in bioproduction applications was assessed under single and mixed-carbon culture conditions, in which serial ALE strain exhibited superior productivity of acetoin compared to ancestral strains. Together, systems-level analysis elucidated the molecular basis of serial evolution, which hold potential utility in bioproduction applications.

Synthetic 3'-UTR valves for optimal metabolic flux control in Escherichia coli.

As the design of genetic circuitry for synthetic biology becomes more sophisticated, diverse regulatory bioparts are required. Despite their importance, well-characterized 3'-untranslated region (3'-UTR) bioparts are limited. Thus, transcript 3'-ends require further investigation to understand the underlying regulatory role and applications of the 3'-UTR. Here, we revisited the use of Term-Seq in the Escherichia coli strain K-12 MG1655 to enhance our understanding of 3'-UTR regulatory functions and to provide a diverse collection of tunable 3'-UTR bioparts with a wide termination strength range. Comprehensive analysis of 1,629 transcript 3'-end positions revealed multiple 3'-termini classes generated through transcription termination and RNA processing. The examination of individual Rho-independent terminators revealed a reduction in downstream gene expression over a wide range, which led to the design of novel synthetic metabolic valves that control metabolic fluxes in branched pathways. These synthetic metabolic valves determine the optimal balance of heterologous pathways for maximum target biochemical productivity. The regulatory strategy using 3'-UTR bioparts is advantageous over promoter- or 5'-UTR-based transcriptional control as it modulates gene expression at transcription levels without trans-acting element requirements (e.g. transcription factors). Our results provide a foundational platform for 3'-UTR engineering in synthetic biology applications.

RiboRid: A low cost, advanced, and ultra-efficient method to remove ribosomal RNA for bacterial transcriptomics.

RNA sequencing techniques have enabled the systematic elucidation of gene expression (RNA-Seq), transcription start sites (differential RNA-Seq), transcript 3' ends (Term-Seq), and post-transcriptional processes (ribosome profiling). The main challenge of transcriptomic studies is to remove ribosomal RNAs (rRNAs), which comprise more than 90% of the total RNA in a cell. Here, we report a low-cost and robust bacterial rRNA depletion method, RiboRid, based on the enzymatic degradation of rRNA by thermostable RNase H. This method implemented experimental considerations to minimize nonspecific degradation of mRNA and is capable of depleting pre-rRNAs that often comprise a large portion of RNA, even after rRNA depletion. We demonstrated the highly efficient removal of rRNA up to a removal efficiency of 99.99% for various transcriptome studies, including RNA-Seq, Term-Seq, and ribosome profiling, with a cost of approximately $10 per sample. This method is expected to be a robust method for large-scale high-throughput bacterial transcriptomic studies.

Acetogenic bacteria utilize light-driven electrons as an energy source for autotrophic growth.

Acetogenic bacteria use cellular redox energy to convert CO2 to acetate using the Wood-Ljungdahl (WL) pathway. Such redox energy can be derived from electrons generated from H2 as well as from inorganic materials, such as photoresponsive semiconductors. We have developed a nanoparticle-microbe hybrid system in which chemically synthesized cadmium sulfide nanoparticles (CdS-NPs) are displayed on the cell surface of the industrial acetogen Clostridium autoethanogenum The hybrid system converts CO2 into acetate without the need for additional energy sources, such as H2, and uses only light-induced electrons from CdS-NPs. To elucidate the underlying mechanism by which C. autoethanogenum uses electrons generated from external energy sources to reduce CO2, we performed transcriptional analysis. Our results indicate that genes encoding the metal ion or flavin-binding proteins were highly up-regulated under CdS-driven autotrophic conditions along with the activation of genes associated with the WL pathway and energy conservation system. Furthermore, the addition of these cofactors increased the CO2 fixation rate under light-exposure conditions. Our results demonstrate the potential to improve the efficiency of artificial photosynthesis systems based on acetogenic bacteria integrated with photoresponsive nanoparticles.

Functional cooperation of the glycine synthase-reductase and Wood-Ljungdahl pathways for autotrophic growth of Clostridium drakei.

Among CO2-fixing metabolic pathways in nature, the linear Wood-Ljungdahl pathway (WLP) in phylogenetically diverse acetate-forming acetogens comprises the most energetically efficient pathway, requires the least number of reactions, and converts CO2 to formate and then into acetyl-CoA. Despite two genes encoding glycine synthase being well-conserved in WLP gene clusters, the functional role of glycine synthase under autotrophic growth conditions has remained uncertain. Here, using the reconstructed genome-scale metabolic model iSL771 based on the completed genome sequence, transcriptomics, 13C isotope-based metabolite-tracing experiments, biochemical assays, and heterologous expression of the pathway in another acetogen, we discovered that the WLP and the glycine synthase pathway are functionally interconnected to fix CO2, subsequently converting CO2 into acetyl-CoA, acetyl-phosphate, and serine. Moreover, the functional cooperation of the pathways enhances CO2 consumption and cellular growth rates via bypassing reducing power required reactions for cellular metabolism during autotrophic growth of acetogens.

Genome-scale analysis of genetic regulatory elements in Streptomyces avermitilis MA-4680 using transcript boundary information.

BACKGROUND: The gram-positive bacterium, Streptomyces avermitilis, holds industrial importance as the producer of avermectin, a widely used anthelmintic agent, and a heterologous expression host of secondary metabolite-biosynthetic gene clusters. Despite its industrial importance, S. avermitilis' genome organization and regulation of gene expression remain poorly understood. In this study, four different types of Next-Generation Sequencing techniques, including dRNA-Seq, Term-Seq, RNA-Seq and ribosome profiling, were applied to S. avermitilis to determine transcription units of S. avermitilis at a genome-wide level and elucidate regulatory elements for transcriptional and translational control of individual transcription units. RESULT: By applying dRNA-Seq and Term-Seq to S. avermitilis MA-4680, a total of 2361 transcription start sites and 2017 transcript 3'-end positions were identified, respectively, leading to determination of 1601 transcription units encoded in S. avermitilis' genome. Cataloguing the transcription units and integrated analysis of multiple high-throughput data types revealed the presence of diverse regulatory elements for gene expression, such as promoters, 5'-UTRs, terminators, 3'-UTRs and riboswitches. The conserved promoter motifs were identified from 2361 transcription start sites as 5'-TANNNT and 5'-BTGACN for the - 10 and - 35 elements, respectively. The - 35 element and spacer lengths between - 10 and - 35 elements were critical for transcriptional regulation of functionally distinct genes, suggesting the involvement of unique sigma factors. In addition, regulatory sequences recognized by antibiotic regulatory proteins were identified from the transcription start site information. Analysis of the 3'-end of RNA transcript revealed that stem structure formation is a major determinant for transcription termination of most transcription units. CONCLUSIONS: The transcription unit architecture elucidated from the transcripts' boundary information provides insights for unique genetic regulatory mechanisms of S. avermitilis. Our findings will elevate S. avermitilis' potential as a production host for a diverse set of secondary metabolites.

Development of highly characterized genetic bioparts for efficient gene expression in CO2-fixing Eubacterium limosum.

Acetogenic bacteria demonstrate industrial potential for utilizing carbon dioxide (CO2) for biochemical production using the Wood-Ljungdahl pathway. However, the metabolic engineering of acetogenic bacteria has been hampered by the limited number of available genetic bioparts for gene expression. Here, we integrated RNA sequencing, ribosome profiling, differential RNA sequencing, and RNA 3'-end sequencing results of Eubacterium limosum to establish genetic bioparts, such as promoters, 5' untranslated regions, and transcript terminators, to regulate transcriptional and translational expression of genes composing of biosynthetic pathways. In addition, a transformation method for the strain was developed to efficiently deliver the obtained genetic bioparts into cells, resulting in a transformation efficiency of 2.5 × 105 CFU/µg DNA. Using this method, the genetic bioparts were efficiently introduced, and their strengths were measured, which were then applied to optimize the heterologous expression of acetolactate synthase and acetolactate decarboxylase for non-native biochemical acetoin production. The strategy developed in this study is the first report on integrating multi-omics data for biopart development of CO2 or syngas utilizing acetogenic bacteria, which lays a foundation for the efficient production of biochemicals from CO2 or syngas as a carbon feedstock under autotrophic growth conditions.

Adaptive laboratory evolution of Escherichia coli W enhances gamma-aminobutyric acid production using glycerol as the carbon source.

The microbial conversion of glycerol into value-added commodity products has emerged as an attractive means to meet the demands of biosustainability. However, glycerol is a non-preferential carbon source for productive fermentation because of its low energy density. We employed evolutionary and metabolic engineering in tandem to construct an Escherichia coli strain with improved GABA production using glycerol as the feedstock carbon. Adaptive evolution of E. coli W under glycerol-limited conditions for 1300 generations harnessed an adapted strain with a metabolic system optimized for glycerol utilization. Mutation profiling, enzyme kinetic assays, and transcriptome analysis of the adapted strain allowed us to decipher the basis of glycerol adaptation at the molecular level. Importantly, increased substrate influx mediated by the mutant glpK and modulation of intracellular cAMP levels were the key drivers of improved fitness in the glycerol-limited condition. Leveraging the enhanced capability of glycerol utilization in the strain, we constructed a GABA-producing E. coli W-derivative with superior GABA production compared to the wild-type. Furthermore, rationally designed inactivation of the non-essential metabolic genes, including ackA, mgsA, and gabT, in the glycerol-adapted strain improved the final GABA titer and specific productivity by 3.9- and 4.3-fold, respectively, compared with the wild-type.

Systems and synthetic biology to elucidate secondary metabolite biosynthetic gene clusters encoded in Streptomyces genomes.

Covering: 2010 to 2020 Over the last few decades, Streptomyces have been extensively investigated for their ability to produce diverse bioactive secondary metabolites. Recent advances in Streptomyces research have been largely supported by improvements in high-throughput technology 'omics'. From genomics, numerous secondary metabolite biosynthetic gene clusters were predicted, increasing their genomic potential for novel bioactive compound discovery. Additional omics, including transcriptomics, translatomics, interactomics, proteomics and metabolomics, have been applied to obtain a system-level understanding spanning entire bioprocesses of Streptomyces, revealing highly interconnected and multi-layered regulatory networks for secondary metabolism. The comprehensive understanding derived from this systematic information accelerates the rational engineering of Streptomyces to enhance secondary metabolite production, integrated with the exploitation of the highly efficient 'Design-Build-Test-Learn' cycle in synthetic biology. In this review, we describe the current status of omics applications in Streptomyces research to better understand the organism and exploit its genetic potential for higher production of valuable secondary metabolites and novel secondary metabolite discovery.

Primary transcriptome and translatome analysis determines transcriptional and translational regulatory elements encoded in the Streptomyces clavuligerus genome.

Determining transcriptional and translational regulatory elements in GC-rich Streptomyces genomes is essential to elucidating the complex regulatory networks that govern secondary metabolite biosynthetic gene cluster (BGC) expression. However, information about such regulatory elements has been limited for Streptomyces genomes. To address this limitation, a high-quality genome sequence of ß-lactam antibiotic-producing Streptomyces clavuligerus ATCC 27 064 is completed, which contains 7163 newly annotated genes. This provides a fundamental reference genome sequence to integrate multiple genome-scale data types, including dRNA-Seq, RNA-Seq and ribosome profiling. Data integration results in the precise determination of 2659 transcription start sites which reveal transcriptional and translational regulatory elements, including -10 and -35 promoter components specific to sigma (σ) factors, and 5'-untranslated region as a determinant for translation efficiency regulation. Particularly, sequence analysis of a wide diversity of the -35 components enables us to predict potential σ-factor regulons, along with various spacer lengths between the -10 and -35 elements. At last, the primary transcriptome landscape of the ß-lactam biosynthetic pathway is analyzed, suggesting temporal changes in metabolism for the synthesis of secondary metabolites driven by transcriptional regulation. This comprehensive genetic information provides a versatile genetic resource for rational engineering of secondary metabolite BGCs in Streptomyces.

Discovery of novel secondary metabolites encoded in actinomycete genomes through coculture.

Actinomycetes are a rich source of bioactive natural products important for novel drug leads. Recent genome mining approaches have revealed an enormous number of secondary metabolite biosynthetic gene clusters (smBGCs) in actinomycetes. However, under standard laboratory culture conditions, many smBGCs are silent or cryptic. To activate these dormant smBGCs, several approaches, including culture-based or genetic engineering-based strategies, have been developed. Above all, coculture is a promising approach to induce novel secondary metabolite production from actinomycetes by mimicking an ecological habitat where cryptic smBGCs may be activated. In this review, we introduce coculture studies that aim to expand the chemical diversity of actinomycetes, by categorizing the cases by the type of coculture partner. Furthermore, we discuss the current challenges that need to be overcome to support the elicitation of novel bioactive compounds from actinomycetes.

Genome-scale analysis of Acetobacterium bakii reveals the cold adaptation of psychrotolerant acetogens by post-transcriptional regulation.

Acetogens synthesize acetyl-CoA via CO2 or CO fixation, producing organic compounds. Despite their ecological and industrial importance, their transcriptional and post-transcriptional regulation has not been systematically studied. With completion of the genome sequence of Acetobacterium bakii (4.28-Mb), we measured changes in the transcriptome of this psychrotolerant acetogen in response to temperature variations under autotrophic and heterotrophic growth conditions. Unexpectedly, acetogenesis genes were highly up-regulated at low temperatures under heterotrophic, as well as autotrophic, growth conditions. To mechanistically understand the transcriptional regulation of acetogenesis genes via changes in RNA secondary structures of 5'-untranslated regions (5'-UTR), the primary transcriptome was experimentally determined, and 1379 transcription start sites (TSS) and 1100 5'-UTR were found. Interestingly, acetogenesis genes contained longer 5'-UTR with lower RNA-folding free energy than other genes, revealing that the 5'-UTRs control the RNA abundance of the acetogenesis genes under low temperature conditions. Our findings suggest that post-transcriptional regulation via RNA conformational changes of 5'-UTRs is necessary for cold-adaptive acetogenesis.

System-level understanding of gene expression and regulation for engineering secondary metabolite production in Streptomyces.

The gram-positive bacterium, Streptomyces, is noticed for its ability to produce a wide array of pharmaceutically active compounds through secondary metabolism. To discover novel bioactive secondary metabolites and increase the production, Streptomyces species have been extensively studied for the past decades. Among the cellular components, RNA molecules play important roles as the messengers for gene expression and diverse regulations taking place at the RNA level. Thus, the analysis of RNA-level regulation is critical to understanding the regulation of Streptomyces' metabolism and secondary metabolite production. A dramatic advance in Streptomyces research was made recently, by exploiting high-throughput technology to systematically understand RNA levels. In this review, we describe the current status of the system-wide investigation of Streptomyces in terms of RNA, toward expansion of its genetic potential for secondary metabolite synthesis.

Synthetic Biology Approaches in The Development of Engineered Therapeutic Microbes.

Since the intimate relationship between microbes and human health has been uncovered, microbes have been in the spotlight as therapeutic targets for several diseases. Microbes contribute to a wide range of diseases, such as gastrointestinal disorders, diabetes and cancer. However, as host-microbiome interactions have not been fully elucidated, treatments such as probiotic administration and fecal transplantations that are used to modulate the microbial community often cause nonspecific results with serious safety concerns. As an alternative, synthetic biology can be used to rewire microbial networks such that the microbes can function as therapeutic agents. Genetic sensors can be transformed to detect biomarkers associated with disease occurrence and progression. Moreover, microbes can be reprogrammed to produce various therapeutic molecules from the host and bacterial proteins, such as cytokines, enzymes and signaling molecules, in response to a disturbed physiological state of the host. These therapeutic treatment systems are composed of several genetic parts, either identified in bacterial endogenous regulation systems or developed through synthetic design. Such genetic components are connected to form complex genetic logic circuits for sophisticated therapy. In this review, we discussed the synthetic biology strategies that can be used to construct engineered therapeutic microbes for improved microbiome-based treatment.

Peptide Transporter CstA Imports Pyruvate in Escherichia coli K-12.

Pyruvate is an important intermediate of central carbon metabolism and connects a variety of metabolic pathways in Escherichia coli Although the intracellular pyruvate concentration is dynamically altered and tightly balanced during cell growth, the pyruvate transport system remains unclear. Here, we identified a pyruvate transporter in E. coli using high-throughput transposon sequencing. The transposon mutant library (a total of 5 × 105 mutants) was serially grown with a toxic pyruvate analog (3-fluoropyruvate [3FP]) to enrich for transposon mutants lacking pyruvate transport function. A total of 52 candidates were selected on the basis of a stringent enrichment level of transposon insertion frequency in response to 3FP treatment. Subsequently, their pyruvate transporter function was examined by conventional functional assays, such as those measuring growth inhibition by the toxic pyruvate analog and pyruvate uptake activity. The pyruvate transporter system comprises CstA and YbdD, which are known as a peptide transporter and a conserved protein, respectively, whose functions are associated with carbon starvation conditions. In addition to the presence of more than one endogenous pyruvate importer, it has been suggested that the E. coli genome encodes constitutive and inducible pyruvate transporters. Our results demonstrated that CstA and YbdD comprise the constitutive pyruvate transporter system in E. coli, which is consistent with the tentative genomic locus previously suggested and the functional relationship with the extracellular pyruvate sensing system. The identification of this pyruvate transporter system provides valuable genetic information for understanding the complex process of pyruvate metabolism in E. coliIMPORTANCE Pyruvate is an important metabolite as a central node in bacterial metabolism, and its intracellular levels are tightly regulated to maintain its functional roles in highly interconnected metabolic pathways. However, an understanding of the mechanism of how bacterial cells excrete and transport pyruvate remains elusive. Using high-throughput transposon sequencing followed by pyruvate uptake activity testing of the selected candidate genes, we found that a pyruvate transporter system comprising CstA and YbdD, currently annotated as a peptide transporter and a conserved protein, respectively, constitutively transports pyruvate. The identification of the physiological role of the pyruvate transporter system provides valuable genetic information for understanding the complex pyruvate metabolism in Escherichia coli.

Genome-scale analysis of syngas fermenting acetogenic bacteria reveals the translational regulation for its autotrophic growth.

BACKGROUND: Acetogenic bacteria constitute promising biocatalysts for the conversion of CO2/H2 or synthesis gas (H2/CO/CO2) into biofuels and value-added biochemicals. These microorganisms are naturally capable of autotrophic growth via unique acetogenesis metabolism. Despite their biosynthetic potential for commercial applications, a systemic understanding of the transcriptional and translational regulation of the acetogenesis metabolism remains unclear. RESULTS: By integrating genome-scale transcriptomic and translatomic data, we explored the regulatory logic of the acetogenesis to convert CO2 into biomass and metabolites in Eubacterium limosum. The results indicate that majority of genes associated with autotrophic growth including the Wood-Ljungdahl pathway, the reduction of electron carriers, the energy conservation system, and gluconeogenesis were transcriptionally upregulated. The translation efficiency of genes in cellular respiration and electron bifurcation was also highly enhanced. In contrast, the transcriptionally abundant genes involved in the carbonyl branch of the Wood-Ljungdahl pathway, as well as the ion-translocating complex and ATP synthase complex in the energy conservation system, showed decreased translation efficiency. The translation efficiencies of genes were regulated by 5'UTR secondary structure under the autotrophic growth condition. CONCLUSIONS: The results illustrated that the acetogenic bacteria reallocate protein synthesis, focusing more on the translation of genes for the generation of reduced electron carriers via electron bifurcation, rather than on those for carbon metabolism under autotrophic growth.

Applications of CRISPR/Cas System to Bacterial Metabolic Engineering.

The clustered regularly interspaced short palindromic repeats/CRISPR-associated (CRISPR/Cas) adaptive immune system has been extensively used for gene editing, including gene deletion, insertion, and replacement in bacterial and eukaryotic cells owing to its simple, rapid, and efficient activities in unprecedented resolution. Furthermore, the CRISPR interference (CRISPRi) system including deactivated Cas9 (dCas9) with inactivated endonuclease activity has been further investigated for regulation of the target gene transiently or constitutively, avoiding cell death by disruption of genome. This review discusses the applications of CRISPR/Cas for genome editing in various bacterial systems and their applications. In particular, CRISPR technology has been used for the production of metabolites of high industrial significance, including biochemical, biofuel, and pharmaceutical products/precursors in bacteria. Here, we focus on methods to increase the productivity and yield/titer scan by controlling metabolic flux through individual or combinatorial use of CRISPR/Cas and CRISPRi systems with introduction of synthetic pathway in industrially common bacteria including Escherichia coli. Further, we discuss additional useful applications of the CRISPR/Cas system, including its use in functional genomics.

The architecture of ArgR-DNA complexes at the genome-scale in Escherichia coli.

DNA-binding motifs that are recognized by transcription factors (TFs) have been well studied; however, challenges remain in determining the in vivo architecture of TF-DNA complexes on a genome-scale. Here, we determined the in vivo architecture of Escherichia coli arginine repressor (ArgR)-DNA complexes using high-throughput sequencing of exonuclease-treated chromatin-immunoprecipitated DNA (ChIP-exo). The ChIP-exo has a unique peak-pair pattern indicating 5' and 3' ends of ArgR-binding region. We identified 62 ArgR-binding loci, which were classified into three groups, comprising single, double and triple peak-pairs. Each peak-pair has a unique 93 base pair (bp)-long (±2 bp) ArgR-binding sequence containing two ARG boxes (39 bp) and residual sequences. Moreover, the three ArgR-binding modes defined by the position of the two ARG boxes indicate that DNA bends centered between the pair of ARG boxes facilitate the non-specific contacts between ArgR subunits and the residual sequences. Additionally, our approach may also reveal other fundamental structural features of TF-DNA interactions that have implications for studying genome-scale transcriptional regulatory networks.

Functional elucidation of the non-coding RNAs of Kluyveromyces marxianus in the exponential growth phase.

BACKGROUND: Non-coding RNAs (ncRNAs), which perform diverse regulatory roles, have been found in organisms from all superkingdoms of life. However, there have been limited numbers of studies on the functions of ncRNAs, especially in nonmodel organisms such as Kluyveromyces marxianus that is widely used in the field of industrial biotechnology. RESULTS: In this study, we measured changes in transcriptome at three time points during the exponential growth phase of K. marxianus by using strand-specific RNA-seq. We found that approximately 60% of the transcriptome consists of ncRNAs transcribed from antisense and intergenic regions of the genome that were transcribed at lower levels than mRNA. In the transcriptome, a substantial number of long antisense ncRNAs (lancRNAs) are differentially expressed and enriched in carbohydrate and energy metabolism pathways. Furthermore, this enrichment is evolutionarily conserved, at least in yeast. Particularly, the mode of regulation of mRNA/lancRNA pairs is associated with mRNA transcription levels; the correlation between the pairs is positive at high mRNA transcriptional levels and negative at low levels. In addition, significant induction of mRNA and coverage of more than half of the mRNA sequence by a lancRNA strengthens the positive correlation between mRNA/lancRNA pairs. CONCLUSIONS: Transcriptome sequencing of K. marxianus in the exponential growth phase reveals pervasive transcription of ncRNAs with evolutionarily conserved functions. Studies of the mode of regulation of mRNA/lancRNA pairs suggest that induction of lancRNA may be associated with switch-like behavior of mRNA/lancRNA pairs and efficient regulation of the carbohydrate and energy metabolism pathways in the exponential growth phase of K. marxianus being used in industrial applications.