Metagenomic Insights Into Functional and Taxonomic Compositions of an Activated Sludge Microbial Community Treating Leachate of a Completed Landfill: A Pathway-Based Analysis

Abstract

Upcycling wastes into valuable products by mixed microbial communities has recently received considerable attention. Sustainable production of high-value substances from one-carbon (C1) compounds, e.g., methanol supplemented as an external electron donor in bioreactors for wastewater treatment, is a promising application of upcycling. This study undertook a gene-centric approach to screen valuable production potentials from mixed culture biomass, removing organic carbon and nitrogen from landfill leachate. To this end, the microbial community of the activated sludge from a landfill leachate treatment plant and its metabolic potential for the production of seven valuable products were investigated. The DNA extracted from the activated sludge was subjected to shotgun metagenome sequencing to analyze the microbial taxonomy and functions associated with producing the seven products. The functional analysis confirmed that the activated sludge could produce six of the valuable products, ectoine, polyhydroxybutyrate (PHB), zeaxanthin, astaxanthin, acetoin, and 2,3-butanediol. Quantification of the detected functional gene hit numbers for these valuable products as a primary trial identified a potential rate-limiting metabolic pathway, e.g., conversion of L-2,4-diaminobutyrate into N-γ-acetyl-L2,4,-diaminobutyrate during the ectoine biosynthesis. Overall, this study demonstrated that primary screening by the proposed gene-centric approach can be used to evaluate the potential for the production of valuable products using mixed culture or single microbe in engineered systems. The proposed approach can be expanded to sites where water purification is highly required, but resource recovery, or upcycling has not been implemented.

Introduction

Microbial resource management (MRM), which enables a better understanding and harnessing of microorganisms in ecosystems for sustainable development, has been advocated and implemented in engineered systems (Verstraete, 2007). One application of MRM is recovery or upcycling of valuable resources from wastes and wastewater (McCarty et al., 2011). MRM does not hinge upon the use of one species, but rather a myriad of uncultured microorganisms with various functions that are not fully understood. Accordingly, exploring and comprehending microbial community compositions and functions could pave the way for the innovation of novel bioprocesses, contributing to resource recovery, production of valuable materials, and energy-saving by wastewater treatment systems.

Among targeted sites amenable to MRM applications, bioreactors treating industrial water and wastewater, e.g., waste produced during food processing and by livestock industries (Matassa et al., 2015; Spiller et al., 2020) and landfill leachate, are promising. From a societal perspective, the demand for appropriate leachate treatment is expected to surge as the amount of waste increases from 2.0 billion tons in 2016 to an expected 3.4 billion tons by 2050 (Kaza et al., 2018). Landfill leachate consists of small molecular organic substances and biodegradable volatile fatty acids in the early stages of the landfill process, while ammonia and persistent organic substances are the main components in the later stages following microbial decomposition of major organic waste constituents (Costa et al., 2019). Methanol is often used as a C1 electron donor to replenish insufficient supplies of electron donors for nitrogen oxides produced during ammonia oxidation in landfill leachate treatment after long-term implementation to facilitate denitrification (Ilies and Mavinic, 2001). Recent studies have shown that methanotrophs and methylotrophs, which are both C1-assimilating bacteria, synthesize valuable products including biodegradable polymers and active ingredients in cosmetics (Zamakhaeva et al., 2017). Given that leachate treatment plants likely harbor microbial consortia that assimilate C1 compounds and ammonia, leachate can be expected to function as a site at which these simple wastewater constituents can be converted into valuable products.

Methanol supplied to activated sludge of a wastewater treatment plant and indigenously high salinity (1–3%) triggers changes in microbial community compositions (Osaka et al., 2006). Moreover, microbial communities enriched by high methanol concentrations are reportedly capable of producing high-value products such as ectoine (Doronina et al., 2010), polyhydroxybutyrate (PHB; Zamakhaeva et al., 2017), acetoin (3-hydroxy-2-butanone; Schendel et al., 1990; Gogleva et al., 2010), and carotenoids (Berry et al., 2003; Lee et al., 2004; Guo et al., 2017). Despite their benefits (Supplementary Table 1), chemical industries face challenges with cost-effective and high-yield production of these products (Pastor et al., 2010; Xiao and Lu, 2014; Levett et al., 2016; Zhang, 2018). However, the development of a biorefinery platform based on the MRM concept has the potential to overcome the challenges. Therefore, bioprocesses based on mixed microbial communities are potential surrogates to commercial chemical processes. Accordingly, it is important to explore and harness effective microbes harboring the genetic potential to synthesize valuable products. However, the phylogeny and physiology of such microbes and varieties of targeted valuable products remain mostly unexploited.

The design and operation of bioreactors for treating wastewater has been reliant on an empirical top-down approach, in which macro-scale variables, e.g., pH and dissolved oxygen, are precisely controlled because of limited understanding of the intricate and diverse compositions of microbial communities in wastewater (Lawson et al., 2019). However, recent advances in high-throughput sequencing technologies, genome-centric with binning, and multi-omics approaches have enabled analysis of micro-scale microbiological processes, including metabolic networks, and interactions among microbes (Nobu et al., 2015; Lawson et al., 2017). Despite massive and invaluable information on functions and taxonomies of microbial communities, genome-centric and multi-omics approaches, however, are labor-intensive and time-consuming, and require strong expertise on microbial ecology and bioinformatics, which does not allow engineers to fully harness the approaches. Therefore, a simple platform to rapidly explore the possibility for a microbial consortium to produce valuable products is demanded.

This study advocates a simple gene-centric method without binning to screen the capability of synthesizing valuable products. Taking a landfill as an MRM site, the analysis was conducted to elucidate the microbial community composition and functions of activated sludge in a landfill leachate treatment system in which methanol was supplemented as an electron donor for nitrogen removal. Using shotgun metagenomics analysis, the metabolic potential for the production of valuable products was targeted. The specific goals of this study were: (i) to determine if an individual species of bacteria or microbial community in the activated sludge possesses the metabolic potential for the synthesis of targeted products at the gene level; and (ii) to infer potential pathways that could act as rate-limiting steps during production.

Materials and Methods

Sampling

Activated sludge of a landfill leachate treatment plant at a completed landfill (Tokyo, Japan) was collected on 20 July 2018. The landfill was in operation for several decades, during which time the leachate was treated by an activated sludge system in a sequencing batch mode. The sequencing batch reactor (SBR) had an effective volume of 740 m³ and received 74.2 m³/day of leachate on average. One cycle of the SBR operation lasted 2 h and consisted of 10 min for inflow, 45 min for mixing to ensure anoxic conditions, 15 min for aeration, and 50 min for settling and outflow of the treated water. Prior to entry of the influent to the SBR, the leachate was mixed with groundwater at a volumetric ratio of 3 (leachate): 1 (groundwater). Tables 1, 2 summarize the operating conditions, and the leachate characteristics and SBR performance, respectively. The organic carbon concentration was low (BOD < 15 mg/L and COD < 50 mg/L), whereas the ammonium concentration in the leachate was as high as 150 mg-N/L. Because of the imbalance of the carbon/nitrogen ratio, about 5.3 L of 50% v/v methanol was added during each inflow period as an external electron donor for denitrification. The methanol amount allowed a relative loading of COD/T-N (kg-COD/kg-N) of 4.38 ± 1.80, which is higher than the denitrification stoichiometry with methanol.

DNA Extraction

Twenty liters of activated sludge were taken from the SBR and immediately transferred to the laboratory. After thoroughly mixing the activated sludge, 15 mL of wastewater containing 22.2 mg-wet weight-activated sludge were used for DNA extraction by a DNA extraction kit (Fast DNA^TM Spin Kit for Soil, MP-Biomedicals, Santa Ana, CA, United States) according to the manufacturer’s protocol. The DNA concentration was measured using a spectrophotometer (NanoDrop 2000c, Thermo Scientific, Wilmington, DE, United States).

Sequence Processing and Metagenomic Assembly

The sequencing and metagenomic assembly have been described in detail in a previous study (Yasuda et al., 2020). Briefly, sequencing was performed with an Ion Torrent (ION Torrent Ion S5) using an Ion 530^TM Chip. Torrent Suite v 4.4.2 (Thermo Fisher Scientific, Germany) was used for base calling and to run multiplexing according to the manufacturer’s protocol. After quality filtering of metagenomic reads and trimming both sides of the reads to ensure a Q20 quality score, 20.24 million reads were processed, giving 13,005,150 reads. Megahit v 1.1.3 (Li et al., 2015) was used at the minimum contig length of 1,000 bp for assembly, resulting in acquisition of 104.86 Mbp contigs. Mapping the reads to the assembly was accomplished with Bowtie2 v 2.3.4.2 (Langmead and Salzberg, 2013), and the subsequent identification of single-copy bacterial and archaeal genes was achieved using HMMER v 3.2 (Finn et al., 2011). Taxonomy was assigned to contigs using Centrifuge v 1.0.3 (Kim et al., 2016). The latest taxonomic information was used for classification (Waite et al., 2017, 2020; Parks et al., 2018). Default parameters were used for all software unless otherwise specified.

Functional Genes and Metabolic Pathways of Valuable Products

Amino-acid sequence data obtained by the HMMER procedure were used to query the Kyoto Encyclopedia of Genes and Genomes (KEGG; Kanehisa and Goto, 2000) with GhostKOALA v 2.1 (Kanehisa et al., 2016). The existence of functional genes and metabolic pathways associated with the production of valuable products by either microbial communities or a single species of indigenous bacteria were confirmed. Based on a literature survey for valuable products microbiologically synthesized from C1 compounds, the metabolic pathways and related functional genes for the production of ectoine (Pastor et al., 2010), PHB (Orita et al., 2014), carotenoids (astaxanthin and zeaxanthin; Ye and Bhatia, 2012; Zhang, 2018), and acetoin (Nielsen et al., 2010) were selected. Carotenoids are synthesized from acetyl-CoA to terpenoid backbone synthesis via either a mevalonate pathway or a non-mevalonate pathway (Supplementary Figure 1 and Supplementary Table 1). The majority of bacteria possess the latter, whereas some use the former pathway (Hoshino and Gaucher, 2018; Rico et al., 2019). Therefore, this study investigated both pathways.

This research assumed that all eubacteria possess central metabolisms, e.g., glycolysis and the TCA cycle. Essential functional genes and enzymes regulating biosynthesis of each valuable product are listed in Table 3. In addition to these compounds that C1-assimilating bacteria are potentially able to synthesize, 2,3-butanediol, which is a promising fuel that is chemically produced from acetoin (Nielsen et al., 2010; Ji et al., 2011), and lutein, which is a carotenoid mainly produced by plants and is known to play a role in vision enhancement (Ye and Bhatia, 2012; Nwachukwu et al., 2016; Zhang, 2018), were also selected.

Visualization of Figures

Taxonomic classification results (Figures 1, 2 and Supplementary Figure 1) and repertories of metabolic pathways (Figure 3) were visualized with R (R Core Team, 2020), Rstudio (RStudio_Team, 2020), and Tidyverse package (Wickham et al., 2019).

FIGURE 1

FIGURE 2

FIGURE 3

Nucleotide Accession Numbers

The sequencing data generated by metagenomic analysis have been deposited in DDBJ with the following accession numbers: DRA009549, BioProject: PRJDB9252, and BioSample: SAMD00202538.

Results

Landfill Leachate Compositions and Reactor Performances

The SBR-based activated sludge system in this study intermittently received methanol as an external electron donor to replenish the limited amount of organic carbon constituents (average COD concentration of 43.9 ± 9.2 mg/L) from the closed landfill. The leachate contained a chloride concentration of 7,693 ± 1,767 mg/L, and the salinity was calculated at 1.23% by the correlation with electron conductivity and temperature (APHA et al., 1992), which is classified as brackish (Table 2). The average removal efficiencies of COD (including the added methanol) and nitrogen were 95.5 ± 1.6 and 97.2 ± 2.1%, and the volumetric removal rates of COD and T-N were 26.1 ± 8.8 g-COD/m³/day and 16.6 ± 3.1 g-N/m³/day, respectively. Nitrification was nearly complete in the SBR, resulting in an effluent ammonium concentration of below 0.05 mg N/L on average.

Microbial Community of Mixed Culture Biomass in Activated Sludge

The class-level taxonomic composition of the activated sludge biomass is summarized in Figure 1. The class Gamma-proteobacteria was most abundant, accounting for up to 35% of the total population. This was followed by Alpha-proteobacteria (23%), and Actinobacteria (19%). The top 50 most abundant microbes at the genus and species levels and the corresponding relative abundances are shown in Figure 2 and Supplementary Figure 1, respectively. Methyloceanibacter caenitepidi, which is a marine methanol-assimilating bacterium, was the most abundant in the biomass (Supplementary Figure 1). Based on read numbers, its relative abundance was 1.3 times (1.2%) higher than that of the second most predominant species, Rhodopseudomonas palustris (0.9%). As for major guilds for ammonia oxidation, the genera Candidatus Nitrosopumilus and Nitrospira, known as ammonia-oxidizing archaea and complete ammonia-oxidizing bacteria, respectively, were detected (Figure 2 and Supplementary Figure 1).

Metabolic Pathways and Genetic Production Potential of the Mixed Culture

The results of mapping our sample to reference metabolic pathways in the KEGG are summarized in Figure 3. The highest relative gene abundance was associated with carbohydrate metabolism (12%) and genetic information processing (12%), followed by protein families: genetic information processing (11%), and protein families: signaling and cellular processes (10%). The biosynthesis pathways for each valuable product after mapping with mixed culture data are illustrated in Figure 4. The mixed culture biomass in the activated sludge possessed a full set of functional genes to synthesize six of the valuable products, ectoine, PHB, zeaxanthin, astaxanthin, acetoin, and 2,3-butanediol (Figure 4 and Table 4). Comparing the hit numbers of the functional genes for ectoine, acetoin, and 2,3-butanediol for each enzymatic reaction step (Figures 4A,D,E) revealed a substantial difference in hit numbers. Specifically, the lowest and highest gene hit numbers for ectoine biosynthesis were ectA (177) and lysC (12,563) for conversion of L-2,4-diaminobutyrate into N-γ-acetyl-L2,4,-diaminobutyrate and of aspartate into aspartyl phosphate, respectively. These hit numbers were used to designate the degree of the rate-limiting reaction for biosynthesis of a valuable product (e.g., 0.014 = 177/12,563). Acetoin and 2,3-butanediol had the common functional gene alsD, which had a considerably lower hit number (115) than the others, leading to a relative abundance of the lowest/highest gene hit number being 0.0045. In contrast, PHB, zeaxanthin, and astaxanthin had relatively evenly-distributed gene hit numbers for their enzymatic reactions. For PHB biosynthesis, the hit numbers of phaA (2,532) and phaB (2,476), which encode acetyl-CoA acetyltransferase and acetoacetyl-CoA reductase, respectively, were lower than the one of phaC (14,491), which encodes poly (3-hydroxyalkanoate) polymerase. These hit numbers resulted in relative abundances of the lowest/highest gene hit number of 0.175 (phaA/phaC) and 0.171 (phaB/phaC), respectively. The same trend appeared with zeaxanthin and astaxanthin syntheses, resulting in a relative abundance of the lowest/highest hit number of 0.11. Terpenoid backbone biosynthesis pathways after gene-mapping are shown in Supplementary Figure 2. In addition, the results of gene-mapping of each valuable product to the KEGG reference pathways are shown in Supplementary Figures 3–5 and those of the terpenoid backbone biosynthesis are shown in Supplementary Figure 6. The mapping results shown in the figures indicate that the pathways for the production of the compounds are connected.

FIGURE 4

Genetic Potential of a Single Bacterium to Produce Valuable Compounds

Linking the phylogeny with functions revealed that PHB, acetoin, and 2,3-butanediol can be genetically produced by a single bacterial species (Table 4). However, individual species of bacteria lack certain functional genes for the production of ectoine, lutein, zeaxanthin, and astaxanthin, implying that syntrophic interactions between different microbial species are needed for synthesis of these compounds. M. caenitepidi, Janthinobacterium sp. Marseille, and Ca. Filomicrobium marinum W harbored phaABC genes encoding enzymes for PHB biosynthesis, while Ca. promineofilum breve, Caldilinea aerophila, Lonsdalea britannica, and Orrella dioscoreae, possessed genes essential to acetoin production, i.e., ilvM, alsD, and butA (Figure 4). Among the bacteria capable of producing acetoin, Caldilinea aerophila, Lonsdalea britannica, and Orrella dioscoreae harbor butA and BDH, which genetically produce 2,3-butanediol from acetoin.

Discussion

Gene-Centric Approach That Facilitates the Development of Upcycling Bioprocesses From Landfill Leachates

The MRM concept (Verstraete, 2007) paves the way for a paradigm shift in which wastewater can be regarded as a resource for valuable products. Applying a bottom-up approach based on gene-centric analysis to the MRM concept is essential for the feasibility of valuable product biosynthesis based on metabolic network information. Therefore, feasibility assessments based on the metagenome will help to establish innovative bioprocesses for valuable product biosynthesis. However, this approach has not been applied to a bioreactor for nitrification and denitrification in landfill leachate containing high concentrations of ammonia (e.g., several hundred mg N/L) and salt (above > 1%; Díaz et al., 2019; Shu et al., 2019). This study, therefore, advocates a rapid screening procedure based on the gene-centric analysis from metagenomic data in a simple manner. To the best of our knowledge, this is the first metagenomics study that investigates the metabolic potential for valuable product biosynthesis by activated sludge used for leachate treatment at a closed landfill site. The application of gene-centric metagenomic analysis confirmed that the activated sludge possessed the metabolic potential for the biosynthesis of six targeted products. The holistic and exclusive screening of metabolic pathway networks for valuable products based on gene-centric metagenomic analysis revealed (i) the potential for valuable product biosynthesis (Table 4) and (ii) inferred potential rate-limiting steps within the metabolic pathways (Figure 4). Although the taxonomic and functional information attained by the gene-centric analysis is not as profound as that by the genome-centric counterpart, the information provides the implication if a biomass sample possesses metabolic potentials for valuable product production or not. The gene-centric approach implemented in this study will facilitate the development of upcycling bioprocesses from landfill leachates in which organic carbon and nitrogen are abundant.

The gene hit number analysis suggests a potential rate-limiting step during the bio-transformations for value-added compounds, e.g., the conversion from L-2,4-diaminobutylate to N-γ-acetyl-L2,4-diaminobutyrate in a series of ectoine production from aspartate (Figure 4). A challenge in this study stems from the authenticity of the gene hit number as an indicator for rate-limiting steps in the metabolic pathways because the differences in gene expression affecting RNA and protein abundances, and the existence of pseudogenes may also determine enzymatic reactions. Although this challenge is under debate, comparing the gene hit numbers may provide rough measures of bottleneck enzymes in metabolic reactions. The credibility of quantifying gene hit numbers warrants future investigation.

Metabolic Potentials of Activated Sludge Biomass at a Landfill Leachate Treatment Plant

This study demonstrated that the mixed microbial community in a landfill leachate treatment site harbored functional genes associated with the syntheses of ectoine, PHB, zeaxanthin, astaxanthin, acetoin, and 2,3-butanediol (Table 4). The presence of these genes suggests the potential for use of this activated sludge biomass in biosynthesis. As previously demonstrated, microbial consortia in activated sludge have been shown to produce short-chain fatty acids (Li et al., 2014), enzymes (Yu et al., 2009), polyhydroxyalkanoate (PHA; Lee and Yu, 1997), bio-flocculants (Wang et al., 2007), and fuels (Salama et al., 2017). As for ectoine synthesis, the analysis of the hit numbers of each functional gene suggested the conversion of L-2,4-diaminobutylate into N-γ-acetyl-L2,4-diaminobutyrate mediated by L-2,4-diaminobutyric acid acetyltransferase (EctA) as a rate-limiting step. Previous studies of EctA purified from halophilic bacteria revealed that suitable conditions for this process were pH 8.2–9.5, a temperature of 20°C, and a salt concentration of 2.3% (Ono et al., 1999; Richter et al., 2020), which do not coincide with the SBR conditions in this study (pH 7.35 ± 0.08, temperature of 28.8°C ± 1.4°C, and salinity of 1.23%). Moreover, caution should be taken when considering these findings because the conditions required to activate EctA in vitro and in vivo may not be the same. Accordingly, future studies to identify suitable environmental conditions for bacteria possessing EctA are needed.

Potential Valuable Products Produced by a Single Microorganism

The metagenomic analysis of the activated sludge used to treat landfill leachate allowed classification of the production potential for valuable products by (i) individual species of bacteria or (ii) mixed microbial communities (Table 4). The results demonstrated that individual species of bacteria present in the activated sludge of the landfill leachate treatment had the genetic potential to synthesize PHB, acetoin, and 2,3-butanediol among the seven products listed in Table 4.

Our metagenomic analysis unraveled that the species proximal to M. caenitepidi, Ca. Filomicrobium marinum W, and Janthinobacterium sp. Marseille were all found to possess the entire set of genes encoding for PHB synthesis. The metagenomic analysis of the taxonomy and functions revealed for the first time that M. caenitepidi, a marine methanol-assimilating bacterium, was dominant in the activated sludge of the landfill leachate treatment facility (Supplementary Figure 1). The presence of M. caenitepidi aside from the ocean has not been reported, except for one study (Yasuda et al., 2020). M. caenitepidi Gela4 was isolated from marine sediment near a hydrothermal vent containing methane in Kagoshima Bay, Japan (Takeuchi et al., 2014) and from a Belgian North Sea sediment (Vekeman et al., 2016). The genus Filomicrobium was detected in a methanol-fed denitrifying bioreactor from saline water (Rissanen et al., 2016). Ca. Filomicrobium marinum W (Henriques and De Marco, 2015), which was first isolated from North Atlantic surface seawater, is a methylotroph capable of potentially synthesizing PHB. Janthinobacterium isolates acquired from the Antarctic possess functional genes for PHB biosynthesis. Although they thrive in marine environments (Tan et al., 2020), the detection of Janthinobacterium sp. as a PHB producer in the landfill leachate treatment site makes sense as the presence of methanol and salinity provide a comparable environment. These species require greater ecophysiological descriptions to leverage their PHB synthesis potentials, which will be a follow-up aim.

Bacteria harboring phaABC store PHB intracellularly under feast-famine conditions (Inoue et al., 2018). Because the SBR in this study may induce feast–famine conditions over one cycle, these three bacteria accumulated PHB to attain energy for their preferential growth. Accordingly, the conditions likely exerted high selection pressure to preferentially enrich the sludge with these three potential PHB producers. Identification of their presence and functions paves the way for developing a PHB-producing process in a water/wastewater treatment system. However, this preliminary study was designed to identify the functions and compositions of an activated sludge microbiome treating landfill leachate, not to quantify water and intracellular constituents in the SBR. Therefore, future investigations to evaluate long-term SBR operation and the associated rates of production of valuable products in accordance with microbial population dynamics in the activated sludge are warranted.

This study also revealed that the bacteria listed in Table 4 harbor genes capable of producing acetoin and 2,3-butanediol. However, the ecological niches of taxonomic proxies for these bacteria are not consistent. Ca. Promineofilum and Caldilinea aerophile, which belong to the phylum Chloroflexi, are commonly found in municipal wastewater treatment plants (McIlroy et al., 2016) and hot spring sulfur turf mats (Sekiguchi et al., 2003), while Lonsdalea Britannica and Orrella dioscoreae are Gamma-proteobacteria isolated from the leaves of Quercus robur (Li et al., 2017) and Dioscorea sansibarensis (Carlier et al., 2017), respectively. Despite their genetic potentials to produce valuable products, in-depth physiological activity tests for actual production have not been conducted (Sekiguchi et al., 2003; Brady et al., 2012; McIlroy et al., 2016; Carlier et al., 2017; Li et al., 2017). Future research will be undertaken to confirm the product availability, yield, and production rate by bacterial species harboring genes related to PHB, acetoin, and 2,3-butanediol production.

Metagenomic Screening of Promising Bacteria Capable of Producing Valuable Products

In-depth analysis of mixed microbial community compositions and functions by gene-centric metagenomic analysis, proposed in this study, could facilitate a primary screening strategy for exploration and leveraging of promising microbial communities or species. Our gene-centric analysis confirmed that the activated sludge investigated did not harbor the functional genes encoding enzymes to synthesize lutein (Supplementary Figure 5). Moreover, inoculating the activated sludge would not succeed in enrichment with a microbial guild producing lutein. Pre-screening by gene-centric metagenomic analysis allows determination of whether a single microbe (I), a mixed microbial community (II), or no microbes (III) synthesize a valuable product (Figure 5). Such classifications may enable identification of the conditions required for a key mixed microbial community or an individual microbe to be enriched (Cases I and II Figure 5), ultimately allowing the development of an innovative microbial process. Accordingly, implementing bioaugmentation strategies (Laocharoen et al., 2015; Lin et al., 2018) or retrofitting treatment processes for valuable product biosynthesis, analogous to the concept of bioaugmentation technology where nitrifying bacteria are enriched in a tank installed in a sludge recycle line (Salem et al., 2003), are potential applications of this screening process (Figure 5). These findings will also enable research toward isolation or selective enrichment of promising microbes using suitable environmental conditions identified from metagenomic data. Moreover, if screening reveals no microbes have the desired properties, the method developed here will save a great deal of time and labor. The advocated scheme could be beneficial especially for engineers to bridge fundamental and applied research toward resource recovery and upcycling.

FIGURE 5

Conclusion

In this study, gene-centric metagenomic analysis of an activated sludge used to treat landfill leachate with a high nitrogen concentration and methanol supplementation revealed the genetic potential for biosynthesis of valuable products by microorganisms in the activated sludge. Specifically, it was confirmed that PHB, ectoine, astaxanthin, zeaxanthin, acetoin, and 2,3-butanediol could be biosynthesized by the mixed microbial community in the activated sludge. Moreover, the pathway-based analysis implied a rate-limiting step for each product. PHB, acetoin and 2,3-butanediol also had the biosynthetic potential to be produced by individual species of microorganisms in the activated sludge. Furthermore, a marine methanol-assimilating bacterium with the genetic potential for PHB biosynthesis, M. caenitepidi, was shown to be dominant in the activated sludge. Our results demonstrate the potential for application of this method as tool for primary screening of promising microbial species in development of bioprocesses for valuable product biosynthesis.

Overall, primary screening by gene-centric metagenomic analysis can be used to evaluate the potential for valuable compound production by microbial communities or individual species in engineered systems was demonstrated. The analysis specifically provides engineers with an opportunity to preliminarily test a possibility of mixed culture biomass to produce valuable products. Future studies to link the screening by metagenomic analysis as described in this study with enrichment and isolation of promising microbial species are warranted.

References

• 1

  APHA, AWWA, and WEF (1992). Standard Methods for the Examination of Water and Wastewater, 18th Edn. Washington, DC: American Public Health Association.

  • Google Scholar
  • View reference in article

• 2

  BerryA.JanssensD.HuM.JoreJ. P. M.HosteB.CleenwerckI.et al (2003). Paracoccus zeaxanthinifaciens sp. nov., a zeaxanthin-producing bacterium.Int. J. Syst. Evol. Microbiol.53231–238. 10.1099/ijs.0.02368-0

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 3

  BradyC. L.CleenwerckI.DenmanS.VenterS. N.Rodríguez-PalenzuelaP.CoutinhoT. A.et al (2012). Proposal to reclassify Brenneria quercina (Hildebrand and Schroth 1967) Hauben et al. 1999 into a new genus, Lonsdalea gen. nov., as Lonsdalea quercina comb. nov., descriptions of Lonsdalea quercina subsp. quercina comb. nov., Lonsdalea quercina subsp. ib.Int. J. Syst. Evol. Microbiol.621592–1602. 10.1099/ijs.0.035055-0

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 4

  CarlierA.CnockaertM.FehrL.VandammeP.EberlL. (2017). Draft genome and description of Orrella dioscoreae gen. nov. sp. nov., a new species of Alcaligenaceae isolated from leaf acumens of Dioscorea sansibarensis.Syst. Appl. Microbiol.4011–21. 10.1016/j.syapm.2016.10.002

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 5

  CostaA. M.AlfaiaR. G.deS. M.CamposJ. C. (2019). Landfill leachate treatment in Brazil – an overview.J. Environ. Manage.232110–116. 10.1016/j.jenvman.2018.11.006

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 6

  DíazA. I.OulegoP.LacaA.GonzálezJ. M.DíazM. (2019). Metagenomic analysis of bacterial communities from a nitrification–denitrification treatment of landfill leachates.Clean Soil Air Water471–8. 10.1002/clen.201900156

  • CrossRef
  • Google Scholar
  • View reference in article

• 7

  DoroninaN. V.EzhovV. A.BeschastnyiA. P.TrotsenkoY. A. (2010). Biosynthesis of the bioprotectant ectoin by aerobic methylotrophic bacteria from methanol.Appl. Biochem. Microbiol.46173–176. 10.1134/S0003683810020080

  • CrossRef
  • Google Scholar
  • View reference in article

• 8

  FinnR. D.ClementsJ.EddyS. R. (2011). HMMER web server: interactive sequence similarity searching.Nucleic Acids Res.3929–37. 10.1093/nar/gkr367

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 9

  GoglevaA. A.KaparullinaE. N.DoroninaN. V.TrotsenkoY. A. (2010). Methylophilus flavus sp. nov. and Methylophilus luteus sp. nov., aerobic, methylotrophic bacteria associated with plants.Int. J. Syst. Evol. Microbiol.602623–2628. 10.1099/ijs.0.019455-0

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 10

  GuoW.LiD.HeR.WuM.ChenW.GaoF.et al (2017). Synthesizing value-added products from methane by a new Methylomonas.J. Appl. Microbiol.1231214–1227. 10.1111/jam.13581

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 11

  HenriquesA. C.De MarcoP. (2015). Complete genome sequences of two strains of “Candidatus Filomicrobium marinum,” a methanesulfonate-degrading species.Genome Announc.31–2. 10.1128/genomeA.00160-15

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 12

  HoshinoY.GaucherE. A. (2018). On the origin of isoprenoid biosynthesis.Mol. Biol. Evol.352185–2197. 10.1093/molbev/msy120

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 13

  IliesP.MavinicD. S. (2001). The effect of decreased ambient temperature on the biological nitrification and denitrification of a high ammonia landfill leachate.Water Res.352065–2072. 10.1016/S0043-1354(00)00477-2

  • CrossRef
  • Google Scholar
  • View reference in article

• 14

  InoueD.SuzukiY.SawadaK.SeiK. (2018). Polyhydroxyalkanoate accumulation ability and associated microbial community in activated sludge-derived acetate-fed microbial cultures enriched under different temperature and pH conditions.J. Biosci. Bioeng.125339–345. 10.1016/j.jbiosc.2017.09.008

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 15

  JiX.HuangH.OuyangP. (2011). Microbial 2,3-butanediol production: a state-of-the-art review.Biotechnol. Adv.29351–364. 10.1016/j.biotechadv.2011.01.007

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 16

  KanehisaM.GotoS. (2000). KEGG: kyoto encyclopedia of genes and genomes.Nucleic Acids Res.2827–30.

  • Google Scholar

• 17

  KanehisaM.SatoY.MorishimaK. (2016). BlastKOALA and GhostKOALA: KEGG tools for functional characterization of genome and metagenome sequences.J. Mol. Biol.428726–731. 10.1016/j.jmb.2015.11.006

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 18

  KazaS.YaoL.Bhada-TataP.Van WoerdenF. (2018). What a Waste 2.0: A Global Snapshot of Solid Waste Management to 2050.Washington, DC: The World Bank.

  • Google Scholar
  • View reference in article

• 19

  KimD.SongL.BreitwieserF. P.SalzbergS. L. (2016). Centrifuge: rapid and sensitive classification of metagenomic sequences.Genome Res.261721–1729. 10.1101/gr.210641.116

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 20

  LangmeadB.SalzbergS. (2013). Fast gapped-read alignment with Bowtie 2.Nat. Methods9357–359. 10.1038/nmeth.1923.Fast

  • CrossRef
  • Google Scholar
  • View reference in article

• 21

  LaocharoenS.ReungsangA.PlangklangP. (2015). Bioaugmentation of Lactobacillus delbrueckii ssp. bulgaricus TISTR 895 to enhance bio-hydrogen production of Rhodobacter sphaeroides KKU-PS5.Biotechnol. Biofuels81–16. 10.1186/s13068-015-0375-z

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 22

  LawsonC. E.HarcombeW. R.HatzenpichlerR.LindemannS. R.LöfflerF. E.O’MalleyM. A.et al (2019). Common principles and best practices for engineering microbiomes.Nat. Rev. Microbiol.17725–741. 10.1038/s41579-019-0255-9

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 23

  LawsonC. E.WuS.BhattacharjeeA. S.HamiltonJ. J.McMahonK. D.GoelR.et al (2017). Metabolic network analysis reveals microbial community interactions in anammox granules.Nat. Commun.81–12. 10.1038/ncomms15416

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 24

  LeeJ. H.KimY. S.ChoiT. J.LeeW. J.KimY. T. (2004). Paracoccus haeundaensis sp. nov., a Gram-negative, halophilic, astaxanthin-producing bacterium.Int. J. Syst. Evol. Microbiol.541699–1702. 10.1099/ijs.0.63146-0

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 25

  LeeS.YuJ. (1997). Production of biodegradable thermoplastics from municipal sludge by a two-stage bioprocess.Resour. Conserv. Recycl.19151–164. 10.1016/S0921-3449(96)01157-3

  • CrossRef
  • Google Scholar

• 26

  LevettI.BirkettG.DaviesN.BellA.LangfordA.LaycockB.et al (2016). Techno-economic assessment of poly-3-hydroxybutyrate (PHB) production from methane - the case for thermophilic bioprocessing.J. Environ. Chem. Eng.43724–3733. 10.1016/j.jece.2016.07.033

  • CrossRef
  • Google Scholar

• 27

  LiD.LiuC.LuoR.SadakaneK.LamT. (2015). MEGAHIT: an ultra-fast single-node solution for large and complex metagenomics assembly via succinct de Bruijn graph.Bioinformatics311674–1676. 10.1093/bioinformatics/btv033

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 28

  LiX. L.PengY. Z.RenN. Q.LiB. K.ChaiT. Z.ZhangL. (2014). Effect of temperature on short chain fatty acids (SCFAs) accumulation and microbiological transformation in sludge alkaline fermentation with Ca(OH)2 adjustment.Water Res.6134–45. 10.1016/j.watres.2014.03.030

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 29

  LiY.XueH.GuoL. M.KoltayA.Palacio-BielsaA.ChangJ.et al (2017). Elevation of three subspecies of Lonsdalea quercina to species level: Lonsdalea britannica sp. nov., Lonsdalea iberica sp. nov. and Lonsdalea populi sp. nov.Int. J. Syst. Evol. Microbiol.674680–4684. 10.1099/ijsem.0.002353

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 30

  LinZ.SunD.DangY.HolmesD. E. (2018). Significant enhancement of nitrous oxide energy yields from wastewater achieved by bioaugmentation with a recombinant strain of Pseudomonas aeruginosa.Sci. Rep.8:11916. 10.1038/s41598-018-30326-8

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 31

  MatassaS.BatstoneD. J.HülsenT.SchnoorJ.VerstraeteW. (2015). Can direct conversion of used nitrogen to new feed and protein help feed the world?Environ. Sci. Technol.495247–5254. 10.1021/es505432w

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 32

  McCartyP. L.BaeJ.KimJ. (2011). Domestic wastewater treatment as a net energy producer-can this be achieved?Environ. Sci. Technol.457100–7106. 10.1021/es2014264

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 33

  McIlroyS. J.KarstS. M.NierychloM.DueholmM. S.AlbertsenM.KirkegaardR. H.et al (2016). Genomic and in situ investigations of the novel uncultured Chloroflexi associated with 0092 morphotype filamentous bulking in activated sludge.ISME J.102223–2234. 10.1038/ismej.2016.14

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 34

  MisawaN.SatomiY.KondoK.YokoyamaA.KajiwaraS.SaitoT.et al (1995). Structure and functional analysis of a marine bacterial carotenoid biosynthesis gene cluster and astaxanthin biosynthetic pathway proposed at the gene level.J. Bacteriol.177, 6575–6584. 10.1128/jb.177.22.6575-6584.1995

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 35

  NielsenD. R.YoonS.YuanC. J.PratherK. L. J. (2010). Metabolic engineering of acetoin and meso-2,3-butanediol biosynthesis in E. coli.Biotechnol. J.5274–284. 10.1002/biot.200900279

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 36

  NobuM. K.NarihiroT.RinkeC.KamagataY.TringeS. G.WoykeT.et al (2015). Microbial dark matter ecogenomics reveals complex synergistic networks in a methanogenic bioreactor.ISME J.91710–1722. 10.1038/ismej.2014.256

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 37

  NwachukwuI. D.UdenigweC. C.AlukoR. E. (2016). Lutein and zeaxanthin: production technology, bioavailability, mechanisms of action, visual function, and health claim status.Trends Food Sci. Technol.4974–84. 10.1016/j.tifs.2015.12.005

  • CrossRef
  • Google Scholar
  • View reference in article

• 38

  OnoH.SawadaK.KhunajakrN.TaoT.YamamotoM.HiramotoM.et al (1999). Characterization of biosynthetic enzymes for ectoine as a compatible solute in a moderately halophilic eubacterium, Halomonas elongata.J. Bacteriol.18191–99. 10.1128/jb.181.1.91-99.1999

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 39

  OritaI.NishikawaK.NakamuraS.FukuiT. (2014). Biosynthesis of polyhydroxyalkanoate copolymers from methanol by Methylobacterium extorquens AM1 and the engineered strains under cobalt-deficient conditions.Appl. Microbiol. Biotechnol.983715–3725. 10.1007/s00253-013-5490-9

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 40

  OsakaT.YoshieS.TsunedaS.HirataA.IwamiN.InamoriY. (2006). Identification of acetate- or methanol-assimilating bacteria under nitrate-reducing conditions by stable-isotope probing.Microb. Ecol.52253–266. 10.1007/s00248-006-9071-7

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 41

  ParksD. H.ChuvochinaM.WaiteD. W.RinkeC.SkarshewskiA.ChaumeilP. A.et al (2018). A standardized bacterial taxonomy based on genome phylogeny substantially revises the tree of life.Nat. Biotechnol.36:996. 10.1038/nbt.4229

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 42

  PastorJ. M.SalvadorM.ArgandoñaM.BernalV.Reina-BuenoM.CsonkaL. N.et al (2010). Ectoines in cell stress protection: uses and biotechnological production.Biotechnol. Adv.28782–801. 10.1016/j.biotechadv.2010.06.005

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 43

  R Core Team (2020). R: A Language and Environment for Statistical Computing.Vienna: R Foundation Statistical Computing.

  • Google Scholar
  • View reference in article

• 44

  RichterA. A.KobusS.CzechL.HoeppnerA.ZarzyckiJ.ErbT. J.et al (2020). The architecture of the diaminobutyrate acetyltransferase active site provides mechanistic insight into the biosynthesis of the chemical chaperone ectoine.J. Biol. Chem.2952822–2838. 10.1074/jbc.RA119.011277

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 45

  RicoJ.DuquesneK.PetitJ. L.MariageA.DariiE.PeruchF.et al (2019). Exploring natural biodiversity to expand access to microbial terpene synthesis.Microb. Cell Fact.181–10. 10.1186/s12934-019-1074-4

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 46

  RissanenA. J.OjalaA.DernjatinM.JaakkolaJ.TiirolaM. (2016). Methylophaga and Hyphomicrobium can be used as target genera in monitoring saline water methanol-utilizing denitrification.J. Ind. Microbiol. Biotechnol.43, 1647–1657. 10.1007/s10295-016-1839-2

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 47

  RStudio_Team (2020). RStudio: Integrated Development for R.Boston, MA: RStudio.

  • Google Scholar
  • View reference in article

• 48

  SalamaE. S.KuradeM. B.Abou-ShanabR. A. I.El-DalatonyM. M.YangI. S.MinB.et al (2017). Recent progress in microalgal biomass production coupled with wastewater treatment for biofuel generation.Renew. Sustain. Energy Rev.791189–1211. 10.1016/j.rser.2017.05.091

  • CrossRef
  • Google Scholar
  • View reference in article

• 49

  SalemS.BerendsD. H. J. G.HeijnenJ. J.Van LoosdrechtM. C. M. (2003). Bio-augmentation by nitrification with return sludge.Water Res.371794–1804. 10.1016/S0043-1354(02)00550-X

  • CrossRef
  • Google Scholar
  • View reference in article

• 50

  SchendelF. J.BremmonC. E.FlickingerM. C.GuettlerM.HansonR. S. (1990). L-lysine production at 50°C by mutants of a newly isolated and characterized methylotrophic Bacillus sp.Appl. Environ. Microbiol.56963–970. 10.1128/aem.56.4.963-970.1990

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 51

  SekiguchiY.YamadaT.HanadaS.OhashiA.HaradaH.KamagataY. (2003). Anaerolinea thermophila gen. nov., sp. nov. and Caldilinea aerophila gen. nov., sp. nov., novel filamentous thermophiles that represent a previously uncultured lineage of the domain bacteria at the subphylum level.Int. J. Syst. Evol. Microbiol.531843–1851. 10.1099/ijs.0.02699-0

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 52

  ShuD.GuoJ.ZhangB.HeY.WeiG. (2019). rDNA- and rRNA-derived communities present divergent assemblage patterns and functional traits throughout full-scale landfill leachate treatment process trains.Sci. Total Environ.6461069–1079. 10.1016/j.scitotenv.2018.07.388

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 53

  SpillerM.MuysM.PapiniG.SakarikaM.BuyleM.VlaeminckS. E. (2020). Environmental impact of microbial protein from potato wastewater as feed ingredient: comparative consequential life cycle assessment of three production systems and soybean meal.Water Res.1711–12. 10.1016/j.watres.2019.115406

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 54

  TakeuchiM.KatayamaT.YamagishiT.HanadaS.TamakiH.KamagataY.et al (2014). Methyloceanibacter caenitepidi gen. nov., sp. nov., a facultatively methylotrophic bacterium isolated from marine sediments near a hydrothermal vent.Int. J. Syst. Evol. Microbiol.64462–468. 10.1099/ijs.0.053397-0

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 55

  TanI. K. P.FoongC. P.TanH. T.LimH.ZainN. A. A.TanY. C.et al (2020). Polyhydroxyalkanoate (PHA) synthase genes and PHA-associated gene clusters in Pseudomonas spp. and Janthinobacterium spp. isolated from Antarctica.J. Biotechnol.31318–28. 10.1016/j.jbiotec.2020.03.006

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 56

  VekemanB.KerckhofF.CremersG.VosP.De, VandammeP.et al (2016). New methyloceanibacter diversity from North Sea sediments includes methanotroph containing solely the soluble methane monooxygenase.Environ. Microbiol.184523–4536. 10.1111/1462-2920.13485

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 57

  VerstraeteW. (2007). Microbial ecology and environmental biotechnology.ISME J.14–8. 10.1038/ismej.2007.7

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 58

  WaiteD. W.ChuvochinaM.PelikanC.ParksD. H.YilmazP.WagnerM.et al (2020). Proposal to reclassify the proteobacterial classes deltaproteobacteria and oligoflexia, and the phylum thermodesulfobacteria into four phyla reflecting major functional capabilities.Int. J. Syst. Evol. Microbiol.705972–6016. 10.1099/ijsem.0.004213

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 59

  WaiteD. W.VanwonterghemI.RinkeC.ParksD. H.ZhangY.TakaiK.et al (2017). Comparative genomic analysis of the class Epsilonproteobacteria and proposed reclassification to epsilonbacteraeota (phyl. nov.).Front. Microbiol.8:682. 10.3389/fmicb.2017.00682

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 60

  WangS. G.GongW. X.LiuX. W.TianL.YueQ. Y.GaoB. Y. (2007). Production of a novel bioflocculant by culture of Klebsiella mobilis using dairy wastewater.Biochem. Eng. J.3681–86. 10.1016/j.bej.2007.02.003

  • CrossRef
  • Google Scholar
  • View reference in article

• 61

  WickhamH.AverickM.BryanJ.ChangW.McGowanL.FrançoisR.et al (2019). Welcome to the tidyverse.J. Open Source Softw.4:1686. 10.21105/joss.01686

  • CrossRef
  • Google Scholar
  • View reference in article

• 62

  XiaoZ.LuJ. R. (2014). Strategies for enhancing fermentative production of acetoin: a review.Biotechnol. Adv.32492–503. 10.1016/j.biotechadv.2014.01.002

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 63

  YasudaS.SuenagaT.OrschlerL.AgrawalS.LacknerS.TeradaA. (2020). Identification of a metagenome-assembled genome of an uncultured Methyloceanibacter sp. strain acquired from an activated sludge system used for landfill leachate treatment.Microbiol. Resour. Announc.91–2. 10.1128/MRA.00771-20

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 64

  YeV. M.BhatiaS. K. (2012). Pathway engineering strategies for production of beneficial carotenoids in microbial hosts.Biotechnol. Lett.341405–1414. 10.1007/s10529-012-0921-8

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 65

  YuG.HeP.ShaoL.ZhuY. (2009). Enzyme extraction by ultrasound from sludge flocs.J. Environ. Sci.21204–210. 10.1016/S1001-0742(08)62252-4

  • CrossRef
  • Google Scholar
  • View reference in article

• 66

  ZamakhaevaS. A.FedorovD. N.TrotsenkoY. A. (2017). Methylotrophic producers of bioplastics (Review).Appl. Biochem. Microbiol.53389–400. 10.1134/S0003683817040147

  • CrossRef
  • Google Scholar
  • View reference in article

• 67

  ZhangC. (2018). Biosynthesis of Carotenoids and Apocarotenoids by Microorganisms and Their Industrial Potential.London: IntechOpen.

  • Google Scholar
  • View reference in article