From breeding to genome design: A genomic makeover for potatoes.

Potato breeding efforts have long been hindered by the genetic consequences of millennia of clonal propagation. To mitigate genomic constraints, Zhang et al. leverage an unprecedented scale of sequencing and marker-assisted breeding to unlock traits that have not been possible through classical breeding, providing a blueprint for plant genome design.

How to engineer the unknown: Advancing a quantitative and predictive understanding of plant and soil biology to address climate change.

Our basic understanding of carbon cycling in the biosphere remains qualitative and incomplete, precluding our ability to effectively engineer novel solutions to climate change. How can we attempt to engineer the unknown? This challenge has been faced before in plant biology, providing a roadmap to guide future efforts. We use examples from over a century of photosynthesis research to illustrate the key principles that will set future plant engineering on a solid footing, namely, an effort to identify the key control variables, quantify the effects of systematically tuning these variables, and use theory to account for these observations. The main contributions of plant synthetic biology will stem not from delivering desired genotypes but from enabling the kind of predictive understanding necessary to rationally design these genotypes in the first place. Only then will synthetic plant biology be able to live up to its promise.

Carbon isotope fractionation by an ancestral rubisco suggests that biological proxies for CO2 through geologic time should be reevaluated.

The history of Earth's carbon cycle reflects trends in atmospheric composition convolved with the evolution of photosynthesis. Fortunately, key parts of the carbon cycle have been recorded in the carbon isotope ratios of sedimentary rocks. The dominant model used to interpret this record as a proxy for ancient atmospheric CO2 is based on carbon isotope fractionations of modern photoautotrophs, and longstanding questions remain about how their evolution might have impacted the record. Therefore, we measured both biomass (εp) and enzymatic (εRubisco) carbon isotope fractionations of a cyanobacterial strain (Synechococcus elongatus PCC 7942) solely expressing a putative ancestral Form 1B rubisco dating to ≫1 Ga. This strain, nicknamed ANC, grows in ambient pCO2 and displays larger εp values than WT, despite having a much smaller εRubisco (17.23 ± 0.61‰ vs. 25.18 ± 0.31‰, respectively). Surprisingly, ANC εp exceeded ANC εRubisco in all conditions tested, contradicting prevailing models of cyanobacterial carbon isotope fractionation. Such models can be rectified by introducing additional isotopic fractionation associated with powered inorganic carbon uptake mechanisms present in Cyanobacteria, but this amendment hinders the ability to accurately estimate historical pCO2 from geological data. Understanding the evolution of rubisco and the CO2 concentrating mechanism is therefore critical for interpreting the carbon isotope record, and fluctuations in the record may reflect the evolving efficiency of carbon fixing metabolisms in addition to changes in atmospheric CO2.

Structural and biochemical basis for regiospecificity of the flavonoid glycosyltransferase UGT95A1.

Glycosylation is a predominant strategy plants use to fine-tune the properties of small molecule metabolites to affect their bioactivity, transport, and storage. It is also important in biotechnology and medicine as many glycosides are utilized in human health. Small molecule glycosylation is largely carried out by family 1 glycosyltransferases. Here, we report a structural and biochemical investigation of UGT95A1, a family 1 GT enzyme from Pilosella officinarum that exhibits a strong, unusual regiospecificity for the 3'-O position of flavonoid acceptor substrate luteolin. We obtained an apo crystal structure to help drive the analyses of a series of binding site mutants, revealing that while most residues are tolerant to mutations, key residues M145 and D464 are important for overall glycosylation activity. Interestingly, E347 is crucial for maintaining the strong preference for 3'-O glycosylation, while R462 can be mutated to increase regioselectivity. The structural determinants of regioselectivity were further confirmed in homologous enzymes. Our study also suggests that the enzyme contains large, highly dynamic, disordered regions. We showed that while most disordered regions of the protein have little to no implication in catalysis, the disordered regions conserved among investigated homologs are important to both the overall efficiency and regiospecificity of the enzyme. This report represents a comprehensive in-depth analysis of a family 1 GT enzyme with a unique substrate regiospecificity and may provide a basis for enzyme functional prediction and engineering.

Cynipid wasps systematically reprogram host metabolism and restructure cell walls in developing galls.

Many insects have evolved the ability to manipulate plant growth to generate extraordinary structures called galls, in which insect larva can develop while being sheltered and feeding on the plant. In particular, cynipid (Hymenoptera: Cynipidae) wasps have evolved to form morphologically complex galls and generate an astonishing array of gall shapes, colors, and sizes. However, the biochemical basis underlying these remarkable cellular and developmental transformations remains poorly understood. A key determinant in plant cellular development is cell wall deposition that dictates the physical form and physiological function of newly developing cells, tissues, and organs. However, it is unclear to what degree cell walls are restructured to initiate and support the formation of new gall tissue. Here, we characterize the molecular alterations underlying gall development using a combination of metabolomic, histological, and biochemical techniques to elucidate how valley oak (Quercus lobata) leaf cells are reprogrammed to form galls. Strikingly, gall development involves an exceptionally coordinated spatial deposition of lignin and xylan to form de novo gall vasculature. Our results highlight how cynipid wasps can radically change the metabolite profile and restructure the cell wall to enable the formation of galls, providing insights into the mechanism of gall induction and the extent to which plants can be entirely reprogrammed to form unique structures and organs.

A dolabralexin-deficient mutant provides insight into specialized diterpenoid metabolism in maize.

Two major groups of specialized metabolites in maize (Zea mays), termed kauralexins and dolabralexins, serve as known or predicted diterpenoid defenses against pathogens, herbivores, and other environmental stressors. To consider the physiological roles of the recently discovered dolabralexin pathway, we examined dolabralexin structural diversity, tissue-specificity, and stress-elicited production in a defined biosynthetic pathway mutant. Metabolomics analyses support a larger number of dolabralexin pathway products than previously known. We identified dolabradienol as a previously undetected pathway metabolite and characterized its enzymatic production. Transcript and metabolite profiling showed that dolabralexin biosynthesis and accumulation predominantly occur in primary roots and show quantitative variation across genetically diverse inbred lines. Generation and analysis of CRISPR-Cas9-derived loss-of-function Kaurene Synthase-Like 4 (Zmksl4) mutants demonstrated dolabralexin production deficiency, thus supporting ZmKSL4 as the diterpene synthase responsible for the conversion of geranylgeranyl pyrophosphate precursors into dolabradiene and downstream pathway products. Zmksl4 mutants further display altered root-to-shoot ratios and root architecture in response to water deficit. Collectively, these results demonstrate dolabralexin biosynthesis via ZmKSL4 as a committed pathway node biochemically separating kauralexin and dolabralexin metabolism, and suggest an interactive role of maize dolabralexins in plant vigor during abiotic stress.

Discovery of photosynthesis genes through whole-genome sequencing of acetate-requiring mutants of Chlamydomonas reinhardtii.

Large-scale mutant libraries have been indispensable for genetic studies, and the development of next-generation genome sequencing technologies has greatly advanced efforts to analyze mutants. In this work, we sequenced the genomes of 660 Chlamydomonas reinhardtii acetate-requiring mutants, part of a larger photosynthesis mutant collection previously generated by insertional mutagenesis with a linearized plasmid. We identified 554 insertion events from 509 mutants by mapping the plasmid insertion sites through paired-end sequences, in which one end aligned to the plasmid and the other to a chromosomal location. Nearly all (96%) of the events were associated with deletions, duplications, or more complex rearrangements of genomic DNA at the sites of plasmid insertion, and together with deletions that were unassociated with a plasmid insertion, 1470 genes were identified to be affected. Functional annotations of these genes were enriched in those related to photosynthesis, signaling, and tetrapyrrole synthesis as would be expected from a library enriched for photosynthesis mutants. Systematic manual analysis of the disrupted genes for each mutant generated a list of 253 higher-confidence candidate photosynthesis genes, and we experimentally validated two genes that are essential for photoautotrophic growth, CrLPA3 and CrPSBP4. The inventory of candidate genes includes 53 genes from a phylogenomically defined set of conserved genes in green algae and plants. Altogether, 70 candidate genes encode proteins with previously characterized functions in photosynthesis in Chlamydomonas, land plants, and/or cyanobacteria; 14 genes encode proteins previously shown to have functions unrelated to photosynthesis. Among the remaining 169 uncharacterized genes, 38 genes encode proteins without any functional annotation, signifying that our results connect a function related to photosynthesis to these previously unknown proteins. This mutant library, with genome sequences that reveal the molecular extent of the chromosomal lesions and resulting higher-confidence candidate genes, will aid in advancing gene discovery and protein functional analysis in photosynthesis.

Plant glycosyltransferases for expanding bioactive glycoside diversity.

Glycosylation is a successful strategy to alter the pharmacological properties of small molecules, and it has emerged as a unique approach to expand the chemical space of natural products that can be explored in drug discovery. Traditionally, most glycosylation events have been carried out chemically, often requiring many protection and deprotection steps to achieve a target molecule. Enzymatic glycosylation by glycosyltransferases could provide an alternative strategy for producing new glycosides. In particular, the glycosyltransferase family has greatly expanded in plants, representing a rich enzymatic resource to mine and expand the diversity of glycosides with novel bioactive properties. This article highlights previous and prospective uses for plant glycosyltransferases in generating bioactive glycosides and altering their pharmacological properties.

Alternative Identification of Glycosides Using MS/MS Matching with an In Silico-Modified Aglycone Mass Spectra Library.

Glycosylation of metabolites serves multiple purposes. Adding sugars makes metabolites more water soluble and improves their biodistribution, stability, and detoxification. In plants, the increase in melting points enables storing otherwise volatile compounds that are released by hydrolysis when needed. Classically, glycosylated metabolites were identified by mass spectrometry (MS/MS) using [M-sugar] neutral losses. Herein, we studied 71 pairs of glycosides with their respective aglycones, including hexose, pentose, and glucuronide moieties. Using liquid chromatography (LC) coupled to electrospray ionization high-resolution mass spectrometry, we detected the classic [M-sugar] product ions for only 68% of glycosides. Instead, we found that most aglycone MS/MS product ions were conserved in the MS/MS spectra of their corresponding glycosides, even when no [M-sugar] neutral losses were observed. We added pentose and hexose units to the precursor masses of an MS/MS library of 3057 aglycones to enable rapid identification of glycosylated natural products with standard MS/MS search algorithms. When searching unknown compounds in untargeted LC-MS/MS metabolomics data of chocolate and tea, we structurally annotated 108 novel glycosides in standard MS-DIAL data processing. We uploaded this new in silico-glycosylated product MS/MS library to GitHub to enable users to detect natural product glycosides without authentic chemical standards.

Heinz-resistant tomato cultivars exhibit a lignin-based resistance to field dodder (Cuscuta campestris) parasitism.

Cuscuta species (dodders) are agriculturally destructive, parasitic angiosperms. These parasitic plants use haustoria as physiological bridges to extract nutrients and water from hosts. Cuscuta campestris has a broad host range and wide geographical distribution. While some wild tomato relatives are resistant, cultivated tomatoes are generally susceptible to C. campestris infestations. However, some specific Heinz tomato (Solanum lycopersicum) hybrid cultivars exhibit resistance to dodders in the field, but their defense mechanism was previously unknown. Here, we discovered that the stem cortex in these resistant lines responds with local lignification upon C. campestris attachment, preventing parasite entry into the host. Lignin Induction Factor 1 (LIF1, an AP2-like transcription factor), SlMYB55, and Cuscuta R-gene for Lignin-based Resistance 1, a CC-NBS-LRR (CuRLR1) are identified as factors that confer host resistance by regulating lignification. SlWRKY16 is upregulated upon C. campestris infestation and potentially negatively regulates LIF1 function. Intriguingly, CuRLR1 may play a role in signaling or function as an intracellular receptor for receiving Cuscuta signals or effectors, thereby regulating lignification-based resistance. In summary, these four regulators control the lignin-based resistance response in specific Heinz tomato cultivars, preventing C. campestris from parasitizing resistant tomatoes. This discovery provides a foundation for investigating multilayer resistance against Cuscuta species and has potential for application in other essential crops attacked by parasitic plants.

Accumulation of high-value bioproducts in planta can improve the economics of advanced biofuels.

Coproduction of high-value bioproducts at biorefineries is a key factor in making biofuels more cost-competitive. One strategy for generating coproducts is to directly engineer bioenergy crops to accumulate bioproducts in planta that can be fractionated and recovered at biorefineries. Here, we develop quantitative insights into the relationship between bioproduct market value and target accumulation rates by investigating a set of industrially relevant compounds already extracted from plant sources with a wide range of market prices and applications, including <$10/kg (limonene, latex, and polyhydroxybutyrate [PHB]), $10 to $100/kg (cannabidiol), and >$100/kg (artemisinin). These compounds are used to identify a range of mass fraction thresholds required to achieve net economic benefits for biorefineries and the additional amounts needed to reach a target $2.50/gal biofuel selling price, using cellulosic ethanol production as a test case. Bioproduct market prices and recovery costs determine the accumulation threshold; we find that moderate- to high-value compounds (i.e., cannabidiol and artemisinin) offer net economic benefits at accumulation rates of just 0.01% dry weight (dwt) to 0.02 dwt%. Lower-value compounds, including limonene, latex, and PHB, require at least an order-of-magnitude greater accumulation to overcome additional extraction and recovery costs (0.3 to 1.2 dwt%). We also find that a diversified approach is critical. For example, global artemisinin demand could be met with fewer than 10 biorefineries, while global demand for latex is equivalent to nearly 180 facilities. Our results provide a roadmap for future plant metabolic engineering efforts aimed at increasing the value derived from bioenergy crops.

Plant-based engineering for production of high-valued natural products.

Covering: up to March 2022Plants are a unique source of complex specialized metabolites, many of which play significant roles in human society. In many cases, however, the availability of these metabolites from naturally occurring sources fails to meet current demands. Thus, there is much interest in expanding the production capacity of target plant molecules. Traditionally, plant breeding, chemical synthesis, and microbial fermentation are considered the primary routes towards large scale production of natural products. Here, we explore the advances, challenges, and future of plant engineering as a complementary path. Although plants are an integral part of our food and agricultural systems and sustain an extensive array of chemical constituents, their complex genetics and physiology have prevented the optimal exploitation of plants as a production chassis. We highlight emerging engineering tools and scientific advances developed in recent years that have improved the prospects of using plants as a sustainable and scalable production platform. We also discuss technological limitations and overall economic outlook of plant-based production of natural products.

Nitrogen Metabolism in Pseudomonas putida: Functional Analysis Using Random Barcode Transposon Sequencing.

Pseudomonas putida KT2440 has long been studied for its diverse and robust metabolisms, yet many genes and proteins imparting these growth capacities remain uncharacterized. Using pooled mutant fitness assays, we identified genes and proteins involved in the assimilation of 52 different nitrogen containing compounds. To assay amino acid biosynthesis, 19 amino acid drop-out conditions were also tested. From these 71 conditions, significant fitness phenotypes were elicited in 672 different genes including 100 transcriptional regulators and 112 transport-related proteins. We divide these conditions into 6 classes, and propose assimilatory pathways for the compounds based on this wealth of genetic data. To complement these data, we characterize the substrate range of three promiscuous aminotransferases relevant to metabolic engineering efforts in vitro. Furthermore, we examine the specificity of five transcriptional regulators, explaining some fitness data results and exploring their potential to be developed into useful synthetic biology tools. In addition, we use manifold learning to create an interactive visualization tool for interpreting our BarSeq data, which will improve the accessibility and utility of this work to other researchers. IMPORTANCE Understanding the genetic basis of P. putida's diverse metabolism is imperative for us to reach its full potential as a host for metabolic engineering. Many target molecules of the bioeconomy and their precursors contain nitrogen. This study provides functional evidence linking hundreds of genes to their roles in the metabolism of nitrogenous compounds, and provides an interactive tool for visualizing these data. We further characterize several aminotransferases, lactamases, and regulators, which are of particular interest for metabolic engineering.

Design of orthogonal regulatory systems for modulating gene expression in plants.

Agricultural biotechnology strategies often require the precise regulation of multiple genes to effectively modify complex plant traits. However, most efforts are hindered by a lack of characterized tools that allow for reliable and targeted expression of transgenes. We have successfully engineered a library of synthetic transcriptional regulators that modulate expression strength in planta. By leveraging orthogonal regulatory systems from Saccharomyces spp., we have developed a strategy for the design of synthetic activators, synthetic repressors, and synthetic promoters and have validated their use in Nicotiana benthamiana and Arabidopsis thaliana. This characterization of contributing genetic elements that dictate gene expression represents a foundation for the rational design of refined synthetic regulators. Our findings demonstrate that these tools provide variation in transcriptional output while enabling the concerted expression of multiple genes in a tissue-specific and environmentally responsive manner, providing a basis for generating complex genetic circuits that process endogenous and environmental stimuli.

In-planta production of the biodegradable polyester precursor 2-pyrone-4,6-dicarboxylic acid (PDC): Stacking reduced biomass recalcitrance with value-added co-product.

2-Pyrone-4,6-dicarboxylic acid (PDC), a chemically stable intermediate that naturally occurs during microbial degradation of lignin by bacteria, represents a promising building block for diverse biomaterials and polyesters such as biodegradable plastics. The lack of a chemical synthesis method has hindered large-scale utilization of PDC and metabolic engineering approaches for its biosynthesis have recently emerged. In this study, we demonstrate a strategy for the production of PDC via manipulation of the shikimate pathway using plants as green factories. In tobacco leaves, we first showed that transient expression of bacterial feedback-resistant 3-deoxy-D-arabinoheptulosonate 7-phosphate synthase (AroG) and 3-dehydroshikimate dehydratase (QsuB) produced high titers of protocatechuate (PCA), which was in turn efficiently converted into PDC upon co-expression of PCA 4,5-dioxygenase (PmdAB) and 4-carboxy-2-hydroxymuconate-6-semialdehyde dehydrogenase (PmdC) derived from Comamonas testosteroni. We validated that stable expression of AroG in Arabidopsis in a genetic background containing the QsuB gene enhanced PCA content in plant biomass, presumably via an increase of the carbon flux through the shikimate pathway. Further, introducing AroG and the PDC biosynthetic genes (PmdA, PmdB, and PmdC) into the Arabidopsis QsuB background, or introducing the five genes (AroG, QsuB, PmdA, PmdB, and PmdC) stacked on a single construct into wild-type plants, resulted in PDC titers of ~1% and ~3% dry weight in plant biomass, respectively. Consistent with previous studies of plants expressing QsuB, all PDC producing lines showed strong reduction in lignin content in stems. This low lignin trait was accompanied with improvements of biomass saccharification efficiency due to reduced cell wall recalcitrance to enzymatic degradation. Importantly, most transgenic lines showed no reduction in biomass yields. Therefore, we conclude that engineering plants with the proposed de-novo PDC pathway provides an avenue to enrich biomass with a value-added co-product while simultaneously improving biomass quality for the supply of fermentable sugars. Implementing this strategy into bioenergy crops has the potential to support existing microbial fermentation approaches that exploit lignocellulosic biomass feedstocks for PDC production.

The More, the Better? Improving VR Firefighting Training System with Realistic Firefighter Tools as Controllers.

A virtual reality (VR) controller plays a key role in supporting interactions between users and the virtual environment. This paper investigates the relationship between the user experience and VR control device modality. We developed a VR firefighting training system integrated with four control devices adapted from real firefighting tools. We iteratively improved the controllers and VR system through a pilot study with six participants and conducted a user study with 30 participants to assess two salient human factor constructs-perceived presence and cognitive load-with three device modality conditions (two standard VR controllers, four real tools, and a hybrid of one real tool and one standard VR controller). We found that having more realistic devices that simulate real tools does not necessarily guarantee a higher level of user experience, highlighting a strategic approach to the development and utilization of VR control devices. Our study gives empirical insights on establishing appropriate combinations of VR control device modality in the context of field-based VR simulation and training.

Fatty Acid and Alcohol Metabolism in Pseudomonas putida: Functional Analysis Using Random Barcode Transposon Sequencing.

With its ability to catabolize a wide variety of carbon sources and a growing engineering toolkit, Pseudomonas putida KT2440 is emerging as an important chassis organism for metabolic engineering. Despite advances in our understanding of the organism, many gaps remain in our knowledge of the genetic basis of its metabolic capabilities. The gaps are particularly noticeable in our understanding of both fatty acid and alcohol catabolism, where many paralogs putatively coding for similar enzymes coexist, making biochemical assignment via sequence homology difficult. To rapidly assign function to the enzymes responsible for these metabolisms, we leveraged random barcode transposon sequencing (RB-Tn-Seq). Global fitness analyses of transposon libraries grown on 13 fatty acids and 10 alcohols produced strong phenotypes for hundreds of genes. Fitness data from mutant pools grown on fatty acids of varying chain lengths indicated specific enzyme substrate preferences and enabled us to hypothesize that DUF1302/DUF1329 family proteins potentially function as esterases. From the data, we also postulate catabolic routes for the two biogasoline molecules isoprenol and isopentanol, which are catabolized via leucine metabolism after initial oxidation and activation with coenzyme A (CoA). Because fatty acids and alcohols may serve as both feedstocks and final products of metabolic-engineering efforts, the fitness data presented here will help guide future genomic modifications toward higher titers, rates, and yields.IMPORTANCE To engineer novel metabolic pathways into P. putida, a comprehensive understanding of the genetic basis of its versatile metabolism is essential. Here, we provide functional evidence for the putative roles of hundreds of genes involved in the fatty acid and alcohol metabolism of the bacterium. These data provide a framework facilitating precise genetic changes to prevent product degradation and to channel the flux of specific pathway intermediates as desired.

Investigation of Bar-seq as a method to study population dynamics of Saccharomyces cerevisiae deletion library during bioreactor cultivation.

BACKGROUND: Despite the latest advancements in metabolic engineering for genome editing and characterization of host performance, the successful development of robust cell factories used for industrial bioprocesses and accurate prediction of the behavior of microbial systems, especially when shifting from laboratory-scale to industrial conditions, remains challenging. To increase the probability of success of a scale-up process, data obtained from thoroughly performed studies mirroring cellular responses to typical large-scale stimuli may be used to derive crucial information to better understand potential implications of large-scale cultivation on strain performance. This study assesses the feasibility to employ a barcoded yeast deletion library to assess genome-wide strain fitness across a simulated industrial fermentation regime and aims to understand the genetic basis of changes in strain physiology during industrial fermentation, and the corresponding roles these genes play in strain performance. RESULTS: We find that mutant population diversity is maintained through multiple seed trains, enabling large scale fermentation selective pressures to act upon the community. We identify specific deletion mutants that were enriched in all processes tested in this study, independent of the cultivation conditions, which include MCK1, RIM11, MRK1, and YGK3 that all encode homologues of mammalian glycogen synthase kinase 3 (GSK-3). Ecological analysis of beta diversity between all samples revealed significant population divergence over time and showed feed specific consequences of population structure. Further, we show that significant changes in the population diversity during fed-batch cultivations reflect the presence of significant stresses. Our observations indicate that, for this yeast deletion collection, the selection of the feeding scheme which affects the accumulation of the fermentative by-product ethanol impacts the diversity of the mutant pool to a higher degree as compared to the pH of the culture broth. The mutants that were lost during the time of most extreme population selection suggest that specific biological processes may be required to cope with these specific stresses. CONCLUSIONS: Our results demonstrate the feasibility of Bar-seq to assess fermentation associated stresses in yeast populations under industrial conditions and to understand critical stages of a scale-up process where variability emerges, and selection pressure gets imposed. Overall our work highlights a promising avenue to identify genetic loci and biological stress responses required for fitness under industrial conditions.

Evolution of the 3-hydroxypropionate bicycle and recent transfer of anoxygenic photosynthesis into the Chloroflexi.

Various lines of evidence from both comparative biology and the geologic record make it clear that the biochemical machinery for anoxygenic photosynthesis was present on early Earth and provided the evolutionary stock from which oxygenic photosynthesis evolved ca. 2.3 billion years ago. However, the taxonomic identity of these early anoxygenic phototrophs is uncertain, including whether or not they remain extant. Several phototrophic bacterial clades are thought to have evolved before oxygenic photosynthesis emerged, including the Chloroflexi, a phylum common across a wide range of modern environments. Although Chloroflexi have traditionally been thought to be an ancient phototrophic lineage, genomics has revealed a much greater metabolic diversity than previously appreciated. Here, using a combination of comparative genomics and molecular clock analyses, we show that phototrophic members of the Chloroflexi phylum are not particularly ancient, having evolved well after the rise of oxygen (ca. 867 million years ago), and thus cannot be progenitors of oxygenic photosynthesis. Similarly, results show that the carbon fixation pathway that defines this clade-the 3-hydroxypropionate bicycle-evolved late in Earth history as a result of a series of horizontal gene transfer events, explaining the lack of geological evidence for this pathway based on the carbon isotope record. These results demonstrate the role of horizontal gene transfer in the recent metabolic innovations expressed within this phylum, including its importance in the development of a novel carbon fixation pathway.

MetaPOAP: presence or absence of metabolic pathways in metagenome-assembled genomes.

Summary: Genome-resolved metagenomics allows the construction of draft microbial genomes from short-read shotgun metagenomics (Metagenome-Assembled Genomes, or MAGs); however, even high-quality MAGs are typically somewhat incomplete and contain a small amount of contaminant sequence, making accurate prediction of metabolic potential challenging. Here, we describe MetaPOAP, an algorithm for probabalistic assessment of the statistical likelihoods for the presence or absence of metabolic pathways in MAGs. Availability and implementation: MetaPOAP is available as Python scripts on GitHub or from the Fischer lab webpage, https://github.com/lmward/MetaPOAP. Supplementary information: Supplementary data are available at Bioinformatics online.