Search

Estimating causal queries, such as changes in protein abundance in response to a perturbation, is a fundamental task in the analysis of biomolecular pathways. The estimation requires experimental measurements on the pathway components. However, in practice many pathway components are left unobserved (latent) because they are either unknown, or difficult to measure. Latent variable models (LVMs) are well-suited for such estimation. Unfortunately, LVM-based estimation of causal queries can be inaccurate when parameters of the latent variables are not uniquely identified, or when the number of latent variables is misspecified. This has limited the use of LVMs for causal inference in biomolecular pathways. In this manuscript, we propose a general and practical approach for LVM-based estimation of causal queries. We prove that, despite the challenges above, LVM-based estimators of causal queries are accurate if the queries are identifiable according to Pearl's do-calculus, and describe an algorithm for its estimation. We illustrate the breadth and the practical utility of this approach for estimating causal queries in four synthetic and two experimental case studies, where structures of biomolecular pathways challenge the existing methods for causal query estimation. The code and the data documenting all the case studies are available at \url{https://github.com/srtaheri/LVMwithDoCalculus}, Comment: https://academic.oup.com/bioinformatics/article/38/Supplement_1/i350/6617530

Optimizing the metabolism of microbial cell factories for yields and titers is a critical step for economically viable production of bioproducts and biofuels. In this process, tuning the expression of individual enzymes to obtain the desired pathway flux is a challenging step, in which data from separate multiomics techniques must be integrated with existing biological knowledge to determine where changes should be made. Following a design-build-test-learn strategy, building on recent advances in Bayesian metabolic control analysis, we identify key enzymes in the oleaginous yeast Yarrowia lipolytica that correlate with the production of itaconate by integrating a metabolic model with multiomics measurements. To this extent, we quantify the uncertainty for a variety of key parameters, known as flux control coefficients (FCCs), needed to improve the bioproduction of target metabolites and statistically obtain key correlations between the measured enzymes and boundary flux. Based on the top five significant FCCs and five correlated enzymes, our results show phosphoglycerate mutase, acetyl-CoA synthetase (ACSm), carbonic anhydrase (HCO3E), pyrophosphatase (PPAm), and homoserine dehydrogenase (HSDxi) enzymes in rate-limiting reactions that can lead to increased itaconic acid production.

Self-driving labs (SDLs) combine fully automated experiments with artificial intelligence (AI) that decides the next set of experiments. Taken to their ultimate expression, SDLs could usher a new paradigm of scientific research, where the world is probed, interpreted, and explained by machines for human benefit. While there are functioning SDLs in the fields of chemistry and materials science, we contend that synthetic biology provides a unique opportunity since the genome provides a single target for affecting the incredibly wide repertoire of biological cell behavior. However, the level of investment required for the creation of biological SDLs is only warranted if directed towards solving difficult and enabling biological questions. Here, we discuss challenges and opportunities in creating SDLs for synthetic biology., Comment: 17 pages, 3 figures. Submitted for publication in Current Opinion in Biotechnology. Updated figure 3 in this version

Lipids play a fundamental role in fungal cell biology, being essential cell membrane components and major targets of antifungal drugs. A deeper knowledge of lipid metabolism is key for developing new drugs and a better understanding of fungal pathogenesis. Here, we built a comprehensive map of the Histoplasma capsulatum lipid metabolic pathway by incorporating proteomic and lipidomic analyses. We performed genetic complementation and overexpression of H. capsulatum genes in Saccharomyces cerevisiae to validate reactions identified in the map and to determine enzymes responsible for catalyzing orphan reactions. The map led to the identification of both the fatty acid desaturation and the sphingolipid biosynthesis pathways as targets for drug development. We found that the sphingolipid biosynthesis inhibitor myriocin, the fatty acid desaturase inhibitor thiocarlide, and the fatty acid analog 10-thiastearic acid inhibit H. capsulatum growth in nanomolar to low-micromolar concentrations. These compounds also reduced the intracellular infection in an alveolar macrophage cell line. Overall, this lipid metabolic map revealed pathways that can be targeted for drug development. IMPORTANCE It is estimated that 150 people die per hour due to the insufficient therapeutic treatments to combat fungal infections. A major hurdle to developing antifungal therapies is the scarce knowledge on the fungal metabolic pathways and mechanisms of virulence. In this context, fungal lipid metabolism is an excellent candidate for developing drugs due to its essential roles in cellular scaffolds, energy storage, and signaling transductors. Here, we provide a detailed map of Histoplasma capsulatum lipid metabolism. The map revealed points of this fungus lipid metabolism that can be targeted for developing antifungal drugs.

Optimizing the metabolism of microbial cell factories for yields and titers is a critical step for economically viable production of bioproducts and biofuels. In this process, tuning the expression of individual enzymes to obtain the desired pathway flux is a challenging step, in which data from separate multiomics techniques must be integrated with existing biological knowledge to determine where changes should be made. Following a design-build-test-learn strategy, building on recent advances in Bayesian metabolic control analysis, we identify key enzymes in the oleaginous yeast

Discovery of a biomarker for patients at high risk of progression to severe Coronavirus Disease 2019 (COVID-19) is critical for clinical management, particularly in areas of the world where widespread vaccine distribution and herd immunity have yet to be achieved. Herein, we characterize peripheral blood from 218 COVID-19 patients with flow cytometry and provide evidence that megakaryocytes are a target for infection by Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2). We demonstrate a positive correlation between infected megakaryocytes expressing the protein calprotectin (also called S100A8/A9), a known marker of COVID-19 severity. Additionally, we show that infected, calprotectin expressing megakaryocytes are correlated with COVID-19 severity and are a prognostic indicator of 30-day clinical outcomes including respiratory failure, thrombotic events, acute kidney injury (AKI), ICU admission, and mechanical ventilation. These findings represent a novel SARS-CoV-2 infection target with significant clinical implications as a biomarker for clinical outcomes associated with severe COVID-19.

We describe a large-scale community effort to build an open-access, interoperable, and computable repository of COVID-19 molecular mechanisms - the COVID-19 Disease Map. We discuss the tools, platforms, and guidelines necessary for the distributed development of its contents by a multi-faceted community of biocurators, domain experts, bioinformaticians, and computational biologists. We highlight the role of relevant databases and text mining approaches in the enrichment and validation of the curated mechanisms. We describe the contents of the Map and their relevance to the molecular pathophysiology of COVID-19 and the analytical and computational modelling approaches that can be applied for mechanistic data interpretation and predictions. We conclude by demonstrating concrete applications of our work through several use cases and highlight new testable hypotheses.; We describe a large-scale community effort to build an open-access, interoperable, and computable repository of COVID-19 molecular mechanisms - the COVID-19 Disease Map. We discuss the tools, platforms, and guidelines necessary for the distributed development of its contents by a multi-faceted community of biocurators, domain experts, bioinformaticians, and computational biologists. We highlight the role of relevant databases and text mining approaches in the enrichment and validation of the curated mechanisms. We describe the contents of the Map and their relevance to the molecular pathophysiology of COVID-19 and the analytical and computational modelling approaches that can be applied for mechanistic data interpretation and predictions. We conclude by demonstrating concrete applications of our work through several use cases and highlight new testable hypotheses.

Counterfactual inference is a useful tool for comparing outcomes of interventions on complex systems. It requires us to represent the system in form of a structural causal model, complete with a causal diagram, probabilistic assumptions on exogenous variables, and functional assignments. Specifying such models can be extremely difficult in practice. The process requires substantial domain expertise, and does not scale easily to large systems, multiple systems, or novel system modifications. At the same time, many application domains, such as molecular biology, are rich in structured causal knowledge that is qualitative in nature. This manuscript proposes a general approach for querying a causal biological knowledge graph, and converting the qualitative result into a quantitative structural causal model that can learn from data to answer the question. We demonstrate the feasibility, accuracy and versatility of this approach using two case studies in systems biology. The first demonstrates the appropriateness of the underlying assumptions and the accuracy of the results. The second demonstrates the versatility of the approach by querying a knowledge base for the molecular determinants of a severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)-induced cytokine storm, and performing counterfactual inference to estimate the causal effect of medical countermeasures for severely ill patients., In proceeding of IEEE, Transactions on Big Data

Biological engineering of microorganisms to produce value-added chemicals is a promising route to sustainable manufacturing. However, overproduction of metabolic intermediates at high titer, rate, and yield from inexpensive substrates is challenging in non-model systems where limited information is available regarding metabolic flux and its control in production conditions. Integrated multi-omic analyses of engineered strains offers an in-depth look at metabolites and proteins directly involved in growth and production of target and non-target bioproducts. Here we applied multi-omic analyses to overproduction of the polymer precursor 3-hydroxypropionic acid (3HP) in the filamentous fungus Aspergillus pseudoterreus. A synthetic pathway consisting of aspartate decarboxylase, beta-alanine pyruvate transaminase, and 3HP dehydrogenase was designed and built for A. pseudoterreus. Strains with single- and multi-copy integration events were isolated and multi-omics analysis consisting of intracellular and extracellular metabolomics and targeted and global proteomics was used to interrogate the strains in shake-flask and bioreactor conditions. Production of a variety of co-products (organic acids and glycerol) and oxidative degradation of 3HP were identified as metabolic pathways competing with 3HP production. Intracellular accumulation of nitrogen as 2,4-diaminobutanoate was identified as an off-target nitrogen sink that may also limit flux through the engineered 3HP pathway. Elimination of the high-expression oxidative 3HP degradation pathway by deletion of a putative malonate semialdehyde dehydrogenase improved the yield of 3HP by 3.4 × after 10 days in shake-flask culture. This is the first report of 3HP production in a filamentous fungus amenable to industrial scale biomanufacturing of organic acids at high titer and low pH.

Systems biology has experienced dramatic growth in the number, size, and complexity of computational models. To reproduce simulation results and reuse models, researchers must exchange unambiguous model descriptions. We review the latest edition of the Systems Biology Markup Language (SBML), a format designed for this purpose. A community of modelers and software authors developed SBML Level 3 over the past decade. Its modular form consists of a core suited to representing reaction‐based models and packages that extend the core with features suited to other model types including constraint‐based models, reaction‐diffusion models, logical network models, and rule‐based models. The format leverages two decades of SBML and a rich software ecosystem that transformed how systems biologists build and interact with models. More recently, the rise of multiscale models of whole cells and organs, and new data sources such as single‐cell measurements and live imaging, has precipitated new ways of integrating data with models. We provide our perspectives on the challenges presented by these developments and how SBML Level 3 provides the foundation needed to support this evolution., Over the past two decades, scientists from different fields have been developing SBML, a standard format for encoding computational models in biology and medicine. This article summarizes recent progress and gives perspectives on emerging challenges.

Climate change is predicted to result in increased drought extent and intensity in the highly productive, former tallgrass prairie region of the continental United States. These soils store large reserves of carbon. The decrease in soil moisture due to drought has largely unknown consequences on soil carbon cycling and other key biogeochemical cycles carried out by soil microbiomes. In this study, we found that soil drying had a significant impact on the structure and function of soil microbial communities, including shifts in expression of specific metabolic pathways, such as those leading toward production of osmoprotectant compounds. This study demonstrates the application of an untargeted multi-omics approach to decipher details of the soil microbial community’s metaphenotypic response to environmental perturbations and should be applicable to studies of other complex microbial systems as well., Climate change is causing shifts in precipitation patterns in the central grasslands of the United States, with largely unknown consequences on the collective physiological responses of the soil microbial community, i.e., the metaphenome. Here, we used an untargeted omics approach to determine the soil microbial community’s metaphenomic response to soil moisture and to define specific metabolic signatures of the response. Specifically, we aimed to develop the technical approaches and metabolic mapping framework necessary for future systematic ecological studies. We collected soil from three locations at the Konza Long-Term Ecological Research (LTER) field station in Kansas, and the soils were incubated for 15 days under dry or wet conditions and compared to field-moist controls. The microbiome response to wetting or drying was determined by 16S rRNA amplicon sequencing, metatranscriptomics, and metabolomics, and the resulting shifts in taxa, gene expression, and metabolites were assessed. Soil drying resulted in significant shifts in both the composition and function of the soil microbiome. In contrast, there were few changes following wetting. The combined metabolic and metatranscriptomic data were used to generate reaction networks to determine the metaphenomic response to soil moisture transitions. Site location was a strong determinant of the response of the soil microbiome to moisture perturbations. However, some specific metabolic pathways changed consistently across sites, including an increase in pathways and metabolites for production of sugars and other osmolytes as a response to drying. Using this approach, we demonstrate that despite the high complexity of the soil habitat, it is possible to generate insight into the effect of environmental change on the soil microbiome and its physiology and functions, thus laying the groundwork for future, targeted studies. IMPORTANCE Climate change is predicted to result in increased drought extent and intensity in the highly productive, former tallgrass prairie region of the continental United States. These soils store large reserves of carbon. The decrease in soil moisture due to drought has largely unknown consequences on soil carbon cycling and other key biogeochemical cycles carried out by soil microbiomes. In this study, we found that soil drying had a significant impact on the structure and function of soil microbial communities, including shifts in expression of specific metabolic pathways, such as those leading toward production of osmoprotectant compounds. This study demonstrates the application of an untargeted multi-omics approach to decipher details of the soil microbial community’s metaphenotypic response to environmental perturbations and should be applicable to studies of other complex microbial systems as well.

Metabolic network modeling of microbial communities provides an in‐depth understanding of community‐wide metabolic and regulatory processes. Compared to single organism analyses, community metabolic network modeling is more complex because it needs to account for interspecies interactions. To date, most approaches focus on reconstruction of high‐quality individual networks so that, when combined, they can predict community behaviors as a result of interspecies interactions. However, this conventional method becomes ineffective for communities whose members are not well characterized and cannot be experimentally interrogated in isolation. Here, we tested a new approach that uses community‐level data as a critical input for the network reconstruction process. This method focuses on directly predicting interspecies metabolic interactions in a community, when axenic information is insufficient. We validated our method through the case study of a bacterial photoautotroph–heterotroph consortium that was used to provide data needed for a community‐level metabolic network reconstruction. Resulting simulations provided experimentally validated predictions of how a photoautotrophic cyanobacterium supports the growth of an obligate heterotrophic species by providing organic carbon and nitrogen sources. J. Cell. Physiol. 231: 2339–2345, 2016. © 2016 The Authors. Journal of Cellular Physiology Published by Wiley Periodicals, Inc.

We report the application of a recently proposed approach for modeling biological systems using a maximum entropy production rate principle in lieu of having in vivo rate constants. The method is applied in four steps: (1) a new ordinary differential equation (ODE) based optimization approach based on Marcelin’s 1910 mass action equation is used to obtain the maximum entropy distribution; (2) the predicted metabolite concentrations are compared to those generally expected from experiments using a loss function from which post-translational regulation of enzymes is inferred; (3) the system is re-optimized with the inferred regulation from which rate constants are determined from the metabolite concentrations and reaction fluxes; and finally (4) a full ODE-based, mass action simulation with rate parameters and allosteric regulation is obtained. From the last step, the power characteristics and resistance of each reaction can be determined. The method is applied to the central metabolism of Neurospora crassa and the flow of material through the three competing pathways of upper glycolysis, the non-oxidative pentose phosphate pathway, and the oxidative pentose phosphate pathway are evaluated as a function of the NADP/NADPH ratio. It is predicted that regulation of phosphofructokinase (PFK) and flow through the pentose phosphate pathway are essential for preventing an extreme level of fructose 1,6-bisphophate accumulation. Such an extreme level of fructose 1,6-bisphophate would otherwise result in a glassy cytoplasm with limited diffusion, dramatically decreasing the entropy and energy production rate and, consequently, biological competitiveness.

Transcriptional and translational feedback loops in fungi and animals drive circadian rhythms in transcript levels that provide output from the clock, but post-transcriptional mechanisms also contribute. To determine the extent and underlying source of this regulation, we applied newly developed analytical tools to a long-duration, deeply sampled, circadian proteomics time course comprising half of the proteome. We found a quarter of expressed proteins are clock regulated, but >40% of these do not arise from clock-regulated transcripts, and our analysis predicts that these protein rhythms arise from oscillations in translational rates. Our data highlighted the impact of the clock on metabolic regulation, with central carbon metabolism reflecting both transcriptional and post-transcriptional control and opposing metabolic pathways showing peak activities at different times of day. The transcription factor CSP-1 plays a role in this metabolic regulation, contributing to the rhythmicity and phase of clock-regulated proteins.

An unexpected July 2014 decision of the U.S. Court of Appeals for the District of Columbia Circuit pertaining to the Committee on Foreign Investment in the United States (CFIUS) merits consideration by both foreign investors and the U.S. national security community. Reversing a decision of the District Court below, the Circuit Court found that CFIUS and the President denied Ralls Corporation (Ralls), a foreign investor, its constitutionally protected due process rights in connection with requiring Ralls to divest acquired property on national security grounds. The case centered on a July 2012 finding by CFIUS that Ralls' acquisition of several windfarm projects located near a U.S. Naval facility in Oregon posed a national security risk to the United States. The President then issued an historic order-just the second of its kind-mandating the unwinding of the completed transaction pursuant to a national security review. Ralls then took CFIUS and the President to court, arguing CFIUS and the President denied its due process. Between two hearings in February and October 2013, the District Court dismissed all of Ralls' claims, finding national security interests outweighed Ralls' interest in learning the specific grounds for the President's divestment order. In the recent decision, the Circuit Court reversed and remanded to the District Court, finding CFIUS and the President failed to give constitutional due process to Ralls. The Circuit Court also ruled that CFIUS must share all unclassified information on which it relied and provide Ralls an opportunity to rebut. As we explore below, the Circuit Court's decision could impact the CFIUS review process and considerations for foreign investors. Even if the decision ultimately does not stand, as we have discussed at length in companion articles, the Ralls case underscores the importance for potential foreign investors in the United States to carefully consider CFIUS issues early in the investment process.

The filamentous fungus Neurospora crassa responds to light in complex ways. To thoroughly study the transcriptional response of this organism to light, RNA-seq was used to analyze capped and polyadenylated mRNA prepared from mycelium grown for 24 hr in the dark and then exposed to light for 0 (control) 15, 60, 120, and 240 min. More than three-quarters of all defined protein coding genes (79%) were expressed in these cells. The increased sensitivity of RNA-seq compared with previous microarray studies revealed that the RNA levels for 31% of expressed genes were affected two-fold or more by exposure to light. Additionally, a large class of mRNAs, enriched for transcripts specifying products involved in rRNA metabolism, showed decreased expression in response to light, indicating a heretofore undocumented effect of light on this pathway. Based on measured changes in mRNA levels, light generally increases cellular metabolism and at the same time causes significant oxidative stress to the organism. To deal with this stress, protective photopigments are made, antioxidants are produced, and genes involved in ribosome biogenesis are transiently repressed.

In 2012, pursuant to a national security review conducted by the Committee on Foreign Investment in the United States (CFIUS), President Obama ordered the unwinding of Ralls Corporation's ("Ralls") acquisition of certain US wind farm companies. Ralls' legal challenge of the President's Order, and the court's opinions to date, highlight important considerations for both foreign investors and the national security community. This article contextualizes the case within a broader debate about how well CFIUS balances openness to foreign investment and protection of national security, points out important gaps in the national security review process, and considers opportunities for greater transparency.

A variety of microbial communities and their genes (the microbiome) exist throughout the human body, with fundamental roles in human health and disease. The National Institutes of Health (NIH)-funded Human Microbiome Project Consortium has established a population-scale framework to develop metagenomic protocols, resulting in a broad range of quality-controlled resources and data including standardized methods for creating, processing and interpreting distinct types of high-throughput metagenomic data available to the scientific community. Here we present resources from a population of 242 healthy adults sampled at 15 or 18 body sites up to three times, which have generated 5,177 microbial taxonomic profiles from 16S ribosomal RNA genes and over 3.5 terabases of metagenomic sequence so far. In parallel, approximately 800 reference strains isolated from the human body have been sequenced. Collectively, these data represent the largest resource describing the abundance and variety of the human microbiome, while providing a framework for current and future studies.

A recent decision of the US District Court for the District of Columbia (District Court) merits consideration by both foreign investors and the US national security community. The case involved a dispute between a foreign investor, Ralls Corporation (Ralls), and the US Government. In July 2012, the Committee on Foreign Investment in the United States (CFIUS) found that Ralls' acquisition of several windfarm projects located near a US Naval facility in Oregon posed a national security risk to the United States. As we discussed at length in a prior article in this Journal, the President then issued an historic order-just the second of its kind-mandating the unwinding of a completed transaction pursuant to a national security review. Ralls then took CFIUS and the President to court, marking the first time an aggrieved foreign investor had done so. After an initial hearing, in February 2013 the District Court dismissed all but one of Ralls' claims, granting an additional hearing only to determine whether Ralls' due process claim merited further consideration. On October 9, 2013, the District Court rejected Ralls' due process argument as well, finding national security interests outweighed Ralls' interest in learning the specific grounds for President Obama's divestment order. For prospective foreign investors, this ruling offers an important lesson, underscoring the importance of seeking prior approval from CFIUS before entering into a "covered transaction."

Significant advances in system-level modeling of cellular behavior can be achieved based on constraints derived from genomic information and on optimality hypotheses. For steady-state models of metabolic networks, mass conservation and reaction stoichiometry impose linear constraints on metabolic fluxes. Different objectives, such as maximization of growth rate or minimization of flux distance from a reference state, can be tested in different organisms and conditions. In particular, we have suggested that the metabolic properties of mutant bacterial strains are best described by an algorithm that performs a minimization of metabolic adjustment (MOMA) upon gene deletion. The increasing availability of many annotated genomes paves the way for a systematic application of these flux balance methods to a large variety of organisms. However, such a high throughput goal crucially depends on our capacity to build metabolic flux models in a fully automated fashion. Here we describe a pipeline for generating models from annotated genomes and discuss the current obstacles to full automation. In addition, we propose a framework for the integration of flux modeling results and high throughput proteomic data, which can potentially help in the inference of whole-cell kinetic parameters.

The filamentous fungus Neurospora crassa played a central role in the development of twentieth-century genetics, biochemistry and molecular biology, and continues to serve as a model organism for eukaryotic biology. Here, we have reconstructed a genome-scale model of its metabolism. This model consists of 836 metabolic genes, 257 pathways, 6 cellular compartments, and is supported by extensive manual curation of 491 literature citations. To aid our reconstruction, we developed three optimization-based algorithms, which together comprise Fast Automated Reconstruction of Metabolism (FARM). These algorithms are: LInear MEtabolite Dilution Flux Balance Analysis (limed-FBA), which predicts flux while linearly accounting for metabolite dilution; One-step functional Pruning (OnePrune), which removes blocked reactions with a single compact linear program; and Consistent Reproduction Of growth/no-growth Phenotype (CROP), which reconciles differences between in silico and experimental gene essentiality faster than previous approaches. Against an independent test set of more than 300 essential/non-essential genes that were not used to train the model, the model displays 93% sensitivity and specificity. We also used the model to simulate the biochemical genetics experiments originally performed on Neurospora by comprehensively predicting nutrient rescue of essential genes and synthetic lethal interactions, and we provide detailed pathway-based mechanistic explanations of our predictions. Our model provides a reliable computational framework for the integration and interpretation of ongoing experimental efforts in Neurospora, and we anticipate that our methods will substantially reduce the manual effort required to develop high-quality genome-scale metabolic models for other organisms., Author Summary Few organisms have been as foundational to the development of modern genetics and cellular metabolism as Neurospora crassa. Given the wealth of knowledge available for this filamentous fungus, the effort required to manually curate a high-quality genome-scale metabolic reconstruction would be daunting. To aid the reconstruction process, we developed three optimization-based algorithms. The first algorithm predicts flux while linearly accounting for metabolite dilution; the second algorithm removes blocked reactions with one compact linear program; and the third algorithm reconciles differences between in silico predictions and experimental observations of mutant viability. We have used these algorithms to develop the first genome-scale metabolic model for Neurospora. We have validated the accuracy of our model against an independent test set of more than 300 growth/no-growth phenotypes, and our model displays 93% sensitivity and specificity. Simulating the biochemical genetics experiments originally performed on Neurospora, we comprehensively predicted essential genes, nutrient rescues of auxotroph mutants and synthetic lethal interactions. With these predictions, we provide potential mechanistic insight into known mutant phenotypes, and testable hypotheses for novel mutant phenotypes. The model, the algorithms and the testable hypotheses provide a computational foundation for the study of Neurospora crassa metabolism.

The sequencing of complete genomes has accelerated biomedical research by providing information about the overall coding capacity of bacterial chromosomes. The original TB annotation resulted in putative functional assignment of similar to 60% of the genes to specific metabolic functions, however, the other 40% of the encoded ORFs where annotated as conserved hypothetical proteins, hypothetical proteins or encoding proteins of unknown function. The TB research community is now at the beginning of the next phases of post-genomics; namely reannotation and functional characterization by targeted experimentation. Arguably, this is the most significant time for basic microbiology in recent history. To foster basic TB research, the Tuberculosis Community Annotation Project (TBCAP) jamboree exercise began the reannotation effort by providing additional information for previous annotations, and refining and substantiating the functional assignment of ORFs and genes within metabolic pathways. The overall goal of the TBCAP 2012 exercise was to gather and compile various data types and use this information with oversight from the scientific community to provide additional information to support the functional annotations of encoding genes. Another objective of this effort was to standardize the publicly accessible Mycobacterium tuberculosis reference sequence and its annotation. The greatest benefit of functional annotation information of genome sequence is that it fuels TB research for drug discovery, diagnostics, vaccine development and epidemiology. (C) 2012 Elsevier Ltd. All rights reserved.

Loa loa, the African eyeworm, is a major filarial pathogen of humans. Unlike most filariae, Loa loa does not contain the obligate intracellular Wolbachia endosymbiont. We describe the 91.4 Mb genome of Loa loa, and the genome of the related filarial parasite Wuchereria bancrofti, and predict 14,907 Loa loa genes based on microfilarial RNA sequencing. By comparing these genomes to that of another filarial parasite, Brugia malayi, and to several other nematode genomes, we demonstrate synteny among filariae but not with non-parasitic nematodes. The Loa loa genome encodes many immunologically relevant genes, as well as protein kinases targeted by drugs currently approved for humans. Despite lacking Wolbachia, Loa loa shows no new metabolic synthesis or transport capabilities compared to other filariae. These results suggest that the role played by Wolbachia in filarial biology is more subtle than previously thought and reveal marked differences between parasitic and non-parasitic nematodes.

The enterococci are Gram-positive lactic acid bacteria that inhabit the gastrointestinal tracts of diverse hosts. However, Enterococcus faecium and E. faecalis have emerged as leading causes of multidrug-resistant hospital-acquired infections. The mechanism by which a well-adapted commensal evolved into a hospital pathogen is poorly understood. In this study, we examined high-quality draft genome data for evidence of key events in the evolution of the leading causes of enterococcal infections, including E. faecalis, E. faecium, E. casseliflavus, and E. gallinarum. We characterized two clades within what is currently classified as E. faecium and identified traits characteristic of each, including variation in operons for cell wall carbohydrate and putative capsule biosynthesis. We examined the extent of recombination between the two E. faecium clades and identified two strains with mosaic genomes. We determined the underlying genetics for the defining characteristics of the motile enterococci E. casseliflavus and E. gallinarum. Further, we identified species-specific traits that could be used to advance the detection of medically relevant enterococci and their identification to the species level., IMPORTANCE The enterococci, in particular, vancomycin-resistant enterococci, have emerged as leading causes of multidrug-resistant hospital-acquired infections. In this study, we examined genome sequence data to define traits with the potential to influence host-microbe interactions and to identify sequences and biochemical functions that could form the basis for the rapid identification of enterococcal species or lineages of importance in clinical and environmental samples.

There was an error in the placement of Figures 1 and 2. The figures were inadvertently swapped. The figure legends are in the correct locations.

Microbial communities carry out the majority of the biochemical activity on the planet, and they play integral roles in processes including metabolism and immune homeostasis in the human microbiome. Shotgun sequencing of such communities' metagenomes provides information complementary to organismal abundances from taxonomic markers, but the resulting data typically comprise short reads from hundreds of different organisms and are at best challenging to assemble comparably to single-organism genomes. Here, we describe an alternative approach to infer the functional and metabolic potential of a microbial community metagenome. We determined the gene families and pathways present or absent within a community, as well as their relative abundances, directly from short sequence reads. We validated this methodology using a collection of synthetic metagenomes, recovering the presence and abundance both of large pathways and of small functional modules with high accuracy. We subsequently applied this method, HUMAnN, to the microbial communities of 649 metagenomes drawn from seven primary body sites on 102 individuals as part of the Human Microbiome Project (HMP). This provided a means to compare functional diversity and organismal ecology in the human microbiome, and we determined a core of 24 ubiquitously present modules. Core pathways were often implemented by different enzyme families within different body sites, and 168 functional modules and 196 metabolic pathways varied in metagenomic abundance specifically to one or more niches within the microbiome. These included glycosaminoglycan degradation in the gut, as well as phosphate and amino acid transport linked to host phenotype (vaginal pH) in the posterior fornix. An implementation of our methodology is available at http://huttenhower.sph.harvard.edu/human​n. This provides a means to accurately and efficiently characterize microbial metabolic pathways and functional modules directly from high-throughput sequencing reads, enabling the determination of community roles in the HMP cohort and in future metagenomic studies., National Institutes of Health (U.S.) (U54HG004968)