PaKman: A Scalable Algorithm for Generating Genomic Contigs on Distributed Memory Machines
Not Available
SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Not Available
Efficient redirection of microbial metabolism into the abundant production of desired bioproducts remains non-trivial. Here, we used flux-based modeling approaches to improve yields of fatty acids in Saccharomyces cerevisiae. We combined 13 C labeling data with comprehensive genome-scale models to shed light onto microbial metabolism and improve metabolic engineering efforts. We concentrated on studying the balance of acetyl-CoA, a precursor metabolite for the biosynthesis of fatty acids. A genome-wide acetyl-CoA balance study showed ATP citrate lyase from Yarrowia lipolytica as a robust source of cytoplasmic acetyl-CoA and malate synthase as a desirable target for downregulation in terms of acetyl-CoA consumption. These genetic modifications were applied to S. cerevisiae WRY2, a strain that is capable of producing 460 mg/L of free fatty acids. With the addition of ATP citrate lyase and downregulation of malate synthase, the engineered strain produced 26% more free fatty acids. Further increases in free fatty acid production of 33% were obtained by knocking out the cytoplasmic glycerol-3-phosphate dehydrogenase, which flux analysis had shown was competing for carbon flux upstream with the carbon flux through the acetyl-CoA production pathway in the cytoplasm. In total, the genetic interventions applied in this work increased fatty acid production by ~70%.
Here, we present the first set of benchmark problems for phase field models that are being developed by the Center for Hierarchical Materials Design (CHiMaD) and the National Institute of Standards and Technology (NIST). While many scientific research areas use a limited set of well-established software, the growing phase field community continues to develop a wide variety of codes and lacks benchmark problems to consistently evaluate the numerical performance of new implementations. Phase field modeling has become significantly more popular as computational power has increased and is now becoming mainstream, driving the need for benchmark problems to validate and verify new implementations. We follow the example set by the micromagnetics community to develop an evolving set of benchmark problems that test the usability, computational resources, numerical capabilities and physical scope of phase field simulation codes. In this paper, we propose two benchmark problems that cover the physics of solute diffusion and growth and coarsening of a second phase via a simple spinodal decomposition model and a more complex Ostwald ripening model. We demonstrate the utility of benchmark problems by comparing the results of simulations performed with two different adaptive time stepping techniques, and we discuss the needs of future benchmark problems. The development of benchmark problems will enable the results of quantitative phase field models to be confidently incorporated into integrated computational materials science and engineering (ICME), an important goal of the Materials Genome Initiative.
DNA ligase 1 (LIG1, Cdc9 in yeast) finalizes eukaryotic nuclear DNA replication by sealing Okazaki fragments using DNA end-joining reactions that strongly discriminate against incorrectly paired DNA substrates. Whether intrinsic ligation fidelity contributes to the accuracy of replication of the nuclear genome is unknown. Here, we show that an engineered low-fidelity LIG1 Cdc9 variant confers a novel mutator phenotype in yeast typified by the accumulation of single base insertion mutations in homonucleotide runs. The rate at which these additions are generated increases upon concomitant inactivation of DNA mismatch repair, or by inactivation of the Fen1 Rad27 Okazaki fragment maturation (OFM) nuclease. Biochemical and structural data establish that LIG1 Cdc9 normally avoids erroneous ligation of DNA polymerase slippage products, and this protection is compromised by mutation of a LIG1 Cdc9 high-fidelity metal binding site. Collectively, our data indicate that high-fidelity DNA ligation is required to prevent insertion mutations, and that this may be particularly critical following strand displacement synthesis during the completion of OFM.
Inspired by the evolution of biological systems, genetic algorithms have been applied to generate solutions for optimization problems in a variety of scientific and engineering disciplines. For a given problem, a suitable genome representation must be defined along with a mutation operator to generate subsequent generations. Unlike natural systems which display a variety of complex rearrangements (e.g. mobile genetic elements), mutation for genetic algorithms commonly utilizes only random point-wise changes. Furthermore, generalizing beyond point-wise mutations poses a key difficulty as useful genome rearrangements depend on the representation and problem domain. To move beyond the limitations of manually defined point-wise changes, here we propose the use of techniques from masked language models to automatically generate mutations. As a first step, common subsequences within a given population are used to generate a vocabulary. The vocabulary is then used to tokenize each genome. A masked language model is trained on the tokenized data in order to generate possible rearrangements (i.e. mutations). In order to illustrate the proposed strategy, we use string representations of molecules and use a genetic algorithm to optimize for drug-likeness and synthesizability. Finally, our results show that moving beyond random point-wise mutations accelerates genetic algorithm optimization.
Targeted insertion of transgenes at pre-determined plant genomic safe harbors provides a desirable alternative to insertions at random sites achieved through conventional methods. Most existing cases of targeted gene insertion in plants have either relied on the presence of a selectable marker gene in the insertion cassette or occurred at low frequency with relatively small DNA fragments (<1.8 kb). Here, we report the use of an optimized CRISPR-Cas9-based method to achieve the targeted insertion of a 5.2 kb carotenoid biosynthesis cassette at two genomic safe harbors in rice. We obtain marker-free rice plants with high carotenoid content in the seeds and no detectable penalty in morphology or yield. Whole-genome sequencing reveals the absence of off-target mutations by Cas9 in the engineered plants. These results demonstrate targeted gene insertion of marker-free DNA in rice using CRISPR-Cas9 genome editing, and offer a promising strategy for genetic improvement of rice and other crops.
Conversion of CO 2 to value-added products presents an opportunity to reduce GHG emissions while generating revenue. Formate, which can be generated by the electrochemical reduction of CO 2 , has been proposed as a promising intermediate compound for microbial upgrading. Here we present progress towards improving the soil bacterium Cupriavidus necator H16, which is capable of growing on formate as its sole source of carbon and energy using the Calvin-Benson-Bassham (CBB) cycle, as a host for formate utilization. Using adaptive laboratory evolution, we generated several isolates that exhibited faster growth rates on formate. The genomes of these isolates were sequenced, and resulting mutations were systematically reintroduced by metabolic engineering, to identify those that improved growth. Further, the metabolic impact of several mutations was investigated further using RNA-seq transcriptomics. We found that deletion of a transcriptional regulator implicated in quorum sensing, PhcA, reduced expression of several operons and led to improved growth on formate. Growth was also improved by deleting large genomic regions present on the extrachromosomal megaplasmid pHG1, particularly two hydrogenase operons and the megaplasmid CBB operon, one of two copies present in the genome. Based on these findings, we generated a rationally engineered ..delta..phcA and megaplasmid-deficient strain that exhibited a 24% faster maximum growth rate on formate. Moreover, this strain achieved a 7% growth rate improvement on succinate and a 19% increase on fructose, demonstrating the broad utility of microbial genome reduction. This strain has the potential to serve as an improved microbial chassis for biological conversion of formate to value-added products.
Poly(ethylene terephthalate) (PET) is one of the most ubiquitous plastics and can be depolymerized through biological and chemo-catalytic routes to its constituent monomers, terephthalic acid (TPA) and ethylene glycol (EG). TPA and EG can be re-synthesized into PET for closed-loop recycling or microbially converted into higher-value products for open-loop recycling. Here, in this study, we expand on our previous efforts engineering and applying Pseudomonas putida KT2440 for PET conversion by employing adaptive laboratory evolution (ALE) to improve TPA catabolism. Three P. putida strains with varying degrees of metabolic engineering for EG catabolism underwent an automation-enabled ALE campaign on TPA, a TPA and EG mixture, and glucose as a control. ALE increased the growth rate on TPA and TPA-EG mixtures by 4.1- and 3.5-fold, respectively, in approximately 350 generations. Evolved isolates were collected at the midpoints and endpoints of 39 independent ALE experiments, and growth rates were increased by 0.15 and 0.20 h -1 on TPA and a TPA-EG, respectively, in the best performing isolates. Whole-genome re-sequencing identified multiple converged mutations, including loss-of-function mutations to global regulators gacS, gacA, and turA along with large duplication and intergenic deletion events that impacted the heterologously-expressed tphAB II catabolic genes. Reverse engineering of these targets confirmed causality, and a strain with all three regulators deleted and second copies of tphAB II and tpaK displayed improved TPA utilization compared to the base strain. Taken together, an iterative strain engineering process involving heterologous pathway engineering, ALE, whole genome sequencing, and genome editing identified five genetic interventions that improve P. putida growth on TPA, aimed at developing enhanced whole-cell biocatalysts for PET upcycling.
Abstract As genomic databases expand and artificial intelligence tools advance, there is a growing demand for efficient characterization of large numbers of proteins. To this end, here we describe a generalizable pipeline for high-throughput protein purification using small-scale expression in E. coli and an affordable liquid-handling robot. This low-cost platform enables the purification of 96 proteins in parallel with minimal waste and is scalable for processing hundreds of proteins weekly per user. We demonstrate the performance of this method with the expression and purification of the leading poly(ethylene terephthalate) hydrolases reported in the literature. Replicate experiments demonstrated reproducibility and enzyme purity and yields (up to 400 µg) sufficient for comprehensive analyses of both thermostability and activity, generating a standardized benchmark dataset for comparing these plastic-degrading enzymes. The cost-effectiveness and ease of implementation of this platform render it broadly applicable to diverse protein characterization challenges in the biological sciences.
Thermoanaerobacterium saccharolyticum is a hemicellulose-degrading thermophilic anaerobe that was previously engineered to produce ethanol at high yield. For this research, a major project was undertaken to develop this organism into an industrial biocatalyst, but the lack of genome information and resources were recognized early on as a key limitation.
Waste carbon from industrial point sources can be captured, stored, and/or transformed using electrochemical conversion or "electrons to molecules" technologies using low-cost renewable electricity. One such process involves electrocatalytic reduction of CO2 to generate formate/formic acid, a C1 carboxylic acid. Formate is a promising potential feedstock for microbial upgrading, as it is water soluble and can be consumed as the sole source of carbon and energy by some microbial species, such as the soil bacterium Cupriavidus necator. Here we will present progress toward improving C. necator as a host for biological conversion of formate to value-added products. Using the power of adaptive laboratory evolution, we were able to isolate mutants of C. necator with significantly faster growth rates on formate. We then sequenced the genomes of these strains, elucidated the metabolic role of the mutations we found, and then used these insights to build rationally engineered strains that outperform even the best evolved isolates. These results highlight the utility of "genome streamlining" as a route for generating platform strains with potential industrial applications.
The development of microbial chassis for the production of a variety of biochemicals and biofuels is a growing area of research. How to efficiently manipulate genetic information to achieve optimal production of these compounds is a key area of focus in the field. In recent years, clustered regularly interspaced palindromic repeats (CRISPR) and its associated proteins (Cas) have become a popular strategy for gene editing and regulation in many organisms due to its versatility and efficacy. Here, we describe methods developed utilizing CRISPR-Cas systems for engineering microbial cell factories (e.g., bacteria and yeast) at the single gene to genome scale for mutagenesis and transcriptional regulation of target genes. Finally, we provide a perspective on the challenges and opportunities for the applications of advanced CRISPR-Cas-based tools for engineering microbial cell factories.
Biorefinery processes for converting lignocellulosic biomass to fuels and chemicals proceed via an integrated series of steps. Biomass is first pretreated and deconstructed using chemical catalysts and/or enzymes to liberate sugar monomers and lignin fragments. Deconstruction is followed by a conversion step in which engineered host organisms assimilate the released sugar monomers and lignin fragments, and produce value-added fuels and chemicals. Over the past couple of decades, a significant amount of work has been done to develop innovative biomass deconstruction and conversion processes that efficiently solubilize biomass, separate lignin from the biomass, maximize yields of bioavailable sugars and lignin fragments and convert the majority of these carbon sources into fuels, commodity chemicals, and materials. Herein, we advocate that advanced in silico approaches provide a theoretical framework for developing efficient processes for lignocellulosic biomass valorization and maximizing yields of sugars and lignin fragments during deconstruction and fuel and chemical titers during conversion. This manuscript surveys the latest developments in lignocellulosic biomass valorization with special attention given to highlighting computational approaches used in process optimization for lignocellulose pretreatment; enzyme engineering for enhanced saccharification and delignification; and prediction of the genome modification necessary for desired pathway fine-tuning to upgrade products from biomass deconstruction into value-added products. Physics-based modeling approaches such as density functional theory calculations and molecular dynamics simulations have been most impactful in studies aimed at exploring the molecular level details of solvent-biomass interactions, reaction mechanisms occurring in biomass-solvent systems, and the catalytic mechanisms and engineering of enzymes involved in biomass degradation. More recently, with ever increasing amounts of data from, for example, advanced mutli-omics experiments, machine learning approaches have begun to make important contributions in synthetic biology and optimization of metabolic pathways for production of biofuels and chemicals.
Production of industrially relevant compounds in microbial cell factories can employ either genomes or plasmids as an expression platform. Selection of plasmids as pathway carriers is advantageous for rapid demonstration but poses a challenge of stability. Yarrowia lipolytica has attracted great attention in the past decade for the biosynthesis of chemicals related to fatty acids at titers attractive to industry, and many genetic tools have been developed to explore its oleaginous potential. Here, our recent studies on the autonomously replicating sequences (ARSs) of nonconventional yeasts revealed that the ARSs from Y. lipolytica showcase a unique structure that includes a previously unannotated sequence (spacer) linking the origin of replication (ORI) and the centromeric (CEN) element and plays a critical role in modulating plasmid behavior. Maintaining a native 645-bp spacer yielded a 2.2-fold increase in gene expression and 1.7-fold higher plasmid stability compared to a more universally employed minimized ARS. Testing the modularity of the ARS sub-elements indicated that plasmid stability exhibits a pronounced cargo dependency. Instability caused both plasmid loss and intramolecular rearrangements. Altogether, our work clarifies the appropriate application of various ARSs for the scientific community and sheds light on a previously unexplored DNA element as a potential target for engineering Y. lipolytica.
Support from this award enabled us to expand critical genomic and epigenomic resources available for the study and engineering of duckweeds, and to complete design-build-test-learn cycles that culminated in the engineering of high oil duckweeds. We shared our findings with the public and scientific community in a series of peer-reviewed publications highlighted within. We successfully reengineered metabolism, dramatically increasing oil content by diverting a substantial portion of accumulated carbon from starch. This sets the stage for the use of fast growing aquatic duckweeds as robust, highly productive sources of biodiesel, with the capacity to outperform soybean oil crops.
Abstract In vivo genetic diversifiers have previously enabled efficient searches of genetic variant fitness landscapes for continuous directed evolution. However, existing genomic diversification modalities for mammalian genomic loci exclusively rely on deaminases to generate transition mutations within target loci, forfeiting access to most missense mutations. Here, we engineer CRISPR-guided error-prone DNA polymerases (EvolvR) to diversify all four nucleotides within genomic loci in mammalian cells. We demonstrate that EvolvR generates both transition and transversion mutations throughout a mutation window of at least 40 bp and implement EvolvR to evolve previously unreported drug-resistantMAP2K1variants via substitutions not achievable with deaminases. Moreover, we discover that the nickase’s mismatch tolerance limits EvolvR’s mutation window and substitution biases in a gRNA-specific fashion. To compensate for gRNA-to-gRNA variability in mutagenesis, we maximize the number of gRNA target sequences by incorporating a PAM-flexible nickase into EvolvR. Finally, we find a strong correlation between predicted free energy changes underlying R-loop formation and EvolvR’s performance using a given gRNA. The EvolvR system diversifies all four nucleotides to enable the evolution of mammalian cells, while nuclease and gRNA-specific properties underlying nickase fidelity can be engineered to further enhance EvolvR’s mutation rates.
Extensive studies in prostate cancer and other malignancies have revealed that L -methionine ( L -Met) and its metabolites play a critical role in tumorigenesis. Preclinical and clinical studies have demonstrated that systemic restriction of serum L -Met, either via partial dietary restriction or with bacterial L -Met–degrading enzymes exerts potent antitumor effects. However, administration of bacterial L -Met–degrading enzymes has not proven practical for human therapy because of problems with immunogenicity. As the human genome does not encode L -Met–degrading enzymes, we engineered the human cystathionine-γ-lyase (hMGL-4.0) to catalyze the selective degradation of L -Met. At therapeutically relevant dosing, hMGL-4.0 reduces serum L -Met levels to >75% for >72 h and significantly inhibits the growth of multiple prostate cancer allografts/xenografts without weight loss or toxicity. We demonstrate that in vitro, hMGL-4.0 causes tumor cell death, associated with increased reactive oxygen species, S-adenosyl-methionine depletion, global hypomethylation, induction of autophagy, and robust poly(ADP-ribose) polymerase (PARP) cleavage indicative of DNA damage and apoptosis.
The variability in phenotypic outcomes among biological replicates in engineered microbial factories presents a captivating mystery. Establishing the association between phenotypic variability and genetic drivers is important to solve this intricate puzzle. Here, we applied a previously developed auxin-inducible depletion of hexokinase 2 as a metabolic engineering strategy for improved nerolidol production in Saccharomyces cerevisiae, and biological replicates exhibit a dichotomy in nerolidol production of either 3.5 or 2.5 g L –1 nerolidol. Harnessing Oxford Nanopore’s long-read genomic sequencing, we reveal a potential genetic cause—the chromosome integration of a 2μ sequence-based yeast episomal plasmid, encoding the expression cassettes for nerolidol synthetic enzymes. This finding was reinforced through chromosome integration revalidation, engineering nerolidol and valencene production strains, and generating a diverse pool of yeast clones, each uniquely fingerprinted by gene copy numbers, plasmid integrations, other genomic rearrangements, protein expression levels, growth rate, and target product productivities. Τhe best clone in two strains produced 3.5 g L –1 nerolidol and ~0.96 g L –1 valencene. Comparable genotypic and phenotypic variations were also generated through the integration of a yeast integrative plasmid lacking 2μ sequences. Our work shows that multiple factors, including plasmid integration status, subchromosomal location, gene copy number, sesquiterpene synthase expression level, and genome rearrangement, together play a complicated determinant role on the productivities of sesquiterpene product. Integration of yeast episomal/integrative plasmids may be used as a versatile method for increasing the diversity and optimizing the efficiency of yeast cell factories, thereby uncovering metabolic control mechanisms.