MINNO: An Open Source Software for Refining Metabolic Networks and Investigating Complex Network Activity Using Empirical Metabolomics Data.

Metabolomics is a powerful tool for uncovering biochemical diversity in a wide range of organisms. Metabolic network modeling is commonly used to frame metabolomics data in the context of a broader biological system. However, network modeling of poorly characterized nonmodel organisms remains challenging due to gene homology mismatches which lead to network architecture errors. To address this, we developed the Metabolic Interactive Nodular Network for Omics (MINNO), a web-based mapping tool that uses empirical metabolomics data to refine metabolic networks. MINNO allows users to create, modify, and interact with metabolic pathway visualizations for thousands of organisms, in both individual and multispecies contexts. Herein, we illustrate the use of MINNO in elucidating the metabolic networks of understudied species, such as those of the Borrelia genus, which cause Lyme and relapsing fever diseases. Using a hybrid genomics-metabolomics modeling approach, we constructed species-specific metabolic networks for threeBorrelia species. Using these empirically refined networks, we were able to metabolically differentiate these species via their nucleotide metabolism, which cannot be predicted from genomic networks. Additionally, using MINNO, we identified 18 missing reactions from the KEGG database, of which nine were supported by the primary literature. These examples illustrate the use of metabolomics for the empirical refining of genetically constructed networks and show how MINNO can be used to study nonmodel organisms.

Metabolic Interactive Nodular Network for Omics (MINNO): Refining and investigating metabolic networks based on empirical metabolomics data.

Metabolomics is a powerful tool for uncovering biochemical diversity in a wide range of organisms, and metabolic network modeling is commonly used to frame results in the context of a broader homeostatic system. However, network modeling of poorly characterized, non-model organisms remains challenging due to gene homology mismatches. To address this challenge, we developed Metabolic Interactive Nodular Network for Omics (MINNO), a web-based mapping tool that takes in empirical metabolomics data to refine metabolic networks for both model and unusual organisms. MINNO allows users to create and modify interactive metabolic pathway visualizations for thousands of organisms, in both individual and multi-species contexts. Herein, we demonstrate an important application of MINNO in elucidating the metabolic networks of understudied species, such as those of the Borrelia genus, which cause Lyme disease and relapsing fever. Using a hybrid genomics-metabolomics modeling approach, we constructed species-specific metabolic networks for three Borrelia species. Using these empirically refined networks, we were able to metabolically differentiate these genetically similar species via their nucleotide and nicotinate metabolic pathways that cannot be predicted from genomic networks. These examples illustrate the use of metabolomics for the empirical refining of genetically constructed networks and show how MINNO can be used to study non-model organisms.

Choreography of Lyme Disease Spirochete Adhesins To Promote Vascular Escape.

The Lyme disease spirochete Borrelia burgdorferi sensu lato can cause a multitude of clinical manifestations because of its ability to disseminate into any organ system via migration through soft tissue, the lymphatic system, and the circulatory system. The latter is believed to constitute the predominant pathway for dissemination to distal sites from the inoculating tick bite. In spite of its importance, the hematogenous dissemination process remains largely uncharacterized, particularly due to difficulties studying this process in a living host and the lack of an in vitro system that recapitulates animal infection. In the current work, we provide the first information regarding the stage of the vascular transmigration pathway where three important adhesins function during invasion of mouse knee joint peripheral tissue from postcapillary venules. Using intravital imaging coupled with genetic experiments employing sequential double infection, we show a complex temporal choreography of P66, decorin binding proteins (DbpA/B), and outer surface protein C (OspC) at discrete steps along the pathway of vascular escape, underscoring the importance of B. burgdorferi adhesins in hematogenous dissemination in the mouse knee joint and the complexity of vascular transmigration by a disseminating pathogen. IMPORTANCE Lyme disease is caused by the spirochete Borrelia burgdorferi, which is transmitted by a bite from an infected tick. Disease development involves a complex series of host-pathogen interactions as well as dissemination of the infecting organisms to sites distal to the original tick bite. The predominant pathway for this is believed to be hematogenous dissemination. The mechanism by which the spirochetes escape circulation is unknown. Here, using intravital microscopy, where the Lyme spirochete can be observed in a living mouse, we have studied the stage in the vascular escape process where each of three surface adhesins functions to facilitate escape of the spirochete from postcapillary venules to invade mouse knee joint peripheral tissue. A complex pattern of involvement at various locations in the multistage process is described using a unique experimental approach that is applicable to other disseminating pathogens.

The Putative Endonuclease Activity of MutL Is Required for the Segmental Gene Conversion Events That Drive Antigenic Variation of the Lyme Disease Spirochete.

The Lyme disease spirochete Borrelia burgdorferi, encodes an elaborate antigenic variation system that promotes the ongoing variation of a major surface lipoprotein, VlsE. Changes in VlsE are continual and always one step ahead of the host acquired immune system, which requires 1-2 weeks to generate specific antibodies. By the time this happens, new VlsE variants have arisen that escape immunosurveillance, providing an avenue for persistent infection. This antigenic variation system is driven by segmental gene conversion events that transfer information from a series of silent cassettes (vls2-16) to the expression locus, vlsE. The molecular details of this process remain elusive. Recombinational switching at vlsE is RecA-independent and the only required factor identified to date is the RuvAB branch migrase. In this work we have used next generation long-read sequencing to analyze the effect of several DNA replication/recombination/repair gene disruptions on the frequency of gene conversions at vlsE and report a requirement for the mismatch repair protein MutL. Site directed mutagenesis of mutL suggests that the putative MutL endonuclease activity is required for recombinational switching at vlsE. This is the first report of an unexpected essential role for MutL in a bacterial recombination system and expands the known function of this protein as well as our knowledge of the details of the novel recombinational switching mechanism for vlsE variation.

VlsE, the nexus for antigenic variation of the Lyme disease spirochete, also mediates early bacterial attachment to the host microvasculature under shear force.

Hematogenous dissemination is a critical step in the evolution of local infection to systemic disease. The Lyme disease (LD) spirochete, which efficiently disseminates to multiple tissues, has provided a model for this process, in particular for the key early event of pathogen adhesion to the host vasculature. This occurs under shear force mediated by interactions between bacterial adhesins and mammalian cell-surface proteins or extracellular matrix (ECM). Using real-time intravital imaging of the Lyme spirochete in living mice, we previously identified BBK32 as the first LD spirochetal adhesin demonstrated to mediate early vascular adhesion in a living mouse; however, deletion of bbk32 resulted in loss of only about half of the early interactions, suggesting the existence of at least one other adhesin (adhesin-X) that promotes early vascular interactions. VlsE, a surface lipoprotein, was identified long ago by its capacity to undergo rapid antigenic variation, is upregulated in the mammalian host and required for persistent infection in immunocompetent mice. In immunodeficient mice, VlsE shares functional overlap with OspC, a multi-functional protein that displays dermatan sulfate-binding activity and is required for joint invasion and colonization. In this research, using biochemical and genetic approaches as well as intravital imaging, we have identified VlsE as adhesin-X; it is a dermatan sulfate (DS) adhesin that efficiently promotes transient adhesion to the microvasculature under shear force via its DS binding pocket. Intravenous inoculation of mice with a low-passage infectious B. burgdorferi strain lacking both bbk32 and vlsE almost completely eliminated transient microvascular interactions. Comparative analysis of binding parameters of VlsE, BBK32 and OspC provides a possible explanation why these three DS adhesins display different functionality in terms of their ability to promote early microvascular interactions.

Strain-specific joint invasion and colonization by Lyme disease spirochetes is promoted by outer surface protein C.

Lyme disease, caused by Borrelia burgdorferi, B. afzelii and B. garinii, is a chronic, multi-systemic infection and the spectrum of tissues affected can vary with the Lyme disease strain. For example, whereas B. garinii infection is associated with neurologic manifestations, B. burgdorferi infection is associated with arthritis. The basis for tissue tropism is poorly understood, but has been long hypothesized to involve strain-specific interactions with host components in the target tissue. OspC (outer surface protein C) is a highly variable outer surface protein required for infectivity, and sequence differences in OspC are associated with variation in tissue invasiveness, but whether OspC directly influences tropism is unknown. We found that OspC binds to the extracellular matrix (ECM) components fibronectin and/or dermatan sulfate in an OspC variant-dependent manner. Murine infection by isogenic B. burgdorferi strains differing only in their ospC coding region revealed that two OspC variants capable of binding dermatan sulfate promoted colonization of all tissues tested, including joints. However, an isogenic strain producing OspC from B. garinii strain PBr, which binds fibronectin but not dermatan sulfate, colonized the skin, heart and bladder, but not joints. Moreover, a strain producing an OspC altered to recognize neither fibronectin nor dermatan sulfate displayed dramatically reduced levels of tissue colonization that were indistinguishable from a strain entirely deficient in OspC. Finally, intravital microscopy revealed that this OspC mutant, in contrast to a strain producing wild type OspC, was defective in promoting joint invasion by B. burgdorferi in living mice. We conclude that OspC functions as an ECM-binding adhesin that is required for joint invasion, and that variation in OspC sequence contributes to strain-specific differences in tissue tropism displayed among Lyme disease spirochetes.

Changing of the guard: How the Lyme disease spirochete subverts the host immune response.

Lyme disease, also known as Lyme borreliosis, is the most common tick-transmitted disease in the Northern Hemisphere. The disease is caused by the bacterial spirochete Borrelia burgdorferi and other related Borrelia species. One of the many fascinating features of this unique pathogen is an elaborate system for antigenic variation, whereby the sequence of the surface-bound lipoprotein VlsE is continually modified through segmental gene conversion events. This perpetual changing of the guard allows the pathogen to remain one step ahead of the acquired immune response, enabling persistent infection. Accordingly, the vls locus is the most evolutionarily diverse genetic element in Lyme disease-causing borreliae. Small stretches of information are transferred from a series of silent cassettes in the vls locus to generate an expressed mosaic vlsE gene version that contains genetic information from several different silent cassettes, resulting in ∼1040 possible vlsE sequences. Yet, despite its extreme evolutionary flexibility, the locus has rigidly conserved structural features. These include a telomeric location of the vlsE gene, an inverse orientation of vlsE and the silent cassettes, the presence of nearly perfect inverted repeats of ∼100 bp near the 5' end of vlsE, and an exceedingly high concentration of G runs in vlsE and the silent cassettes. We discuss the possible roles of these evolutionarily conserved features, highlight recent findings from several studies that have used next-generation DNA sequencing to unravel the switching process, and review advances in the development of a mini-vls system for genetic manipulation of the locus.

Outer membrane protein I is associated with poly-ß-hydroxybutyrate granules and is necessary for optimal polymer accumulation in Azotobacter vinelandii on solid medium.

Azotobacter vinelandii is a soil bacterium that is able to synthesize poly-ß-hydroxybutyrate (PHB), a polymer used to produce biodegradable plastic. PHB is stored in the cytoplasm as granules surrounded by several proteins such as the major phasin PhbP, PHB synthase and PHB depolymerase, among others. Many studies have reported the presence of membrane proteins on PHB granules due to contamination during the polymer extraction procedures. Previously, the outer membrane protein I (OprI) was detected on the polymer granules in A. vinelandii. In this study, by using random transposon mutagenesis, we identified that a mutation in the oprI gene diminished PHB accumulation in A. vinelandii on solid medium. Electron microscopy confirmed the low polymer production by the oprI mutant. Analysis of PHB granules by Tricine-SDS-PAGE revealed that the absence of OprI affected the protein profile of the granules, suggesting that OprI could have a structural role in A. vinelandii. Thus, some membrane proteins on PHB granules may not be artefacts as previously described.

Antigenic variation in the Lyme spirochete: detailed functional assessment of recombinational switching at vlsE in the JD1 strain of Borrelia burgdorferi.

Borrelia burgdorferi is a causative agent of Lyme disease and establishes long-term infection in mammalian hosts. Persistence is promoted by the VlsE antigenic variation system, which generates combinatorial diversity of VlsE through unidirectional, segmental gene conversion from an array of silent cassettes. Here we explore the variants generated by the vls system of strain JD1, which has divergent sequence and structural elements from the type strain B31, the only B. burgdorferi strain in which recombinational switching at vlsE has been studied in detail. We first completed the sequencing of the vls region in JD1, uncovering a previously unreported 114 bp inverted repeat sequence upstream of vlsE. A five-week infection of WT and SCID mice was used for PacBio long read sequencing along with our recently developed VAST pipeline to analyze recombinational switching at vlsE from 40,000 sequences comprising 226,000 inferred recombination events. We show that antigenic variation in B31 and JD1 is highly similar, despite the lack of 17 bp direct repeats in JD1, a somewhat different arrangement of the silent cassettes, divergent inverted repeat sequences and general divergence in the vls sequences. We also present data that strongly suggest that dimerization is required for in vivo functionality of VlsE.

A Borrelia burgdorferi mini-vls system that undergoes antigenic switching in mice: investigation of the role of plasmid topology and the long inverted repeat.

Borrelia burgdorferi evades the host immune system by switching the surface antigen. VlsE, in a process known as antigenic variation. The DNA mechanisms and genetic elements present on the vls locus that participate in the switching process remain to be elucidated. Manipulating the vls locus has been difficult due to its instability on Escherichia coli plasmids. In this study, we generated for the first time a mini-vls system composed of a single silent vlsE variable region (silent cassette 2) through the vlsE gene by performing some cloning steps directly in a highly transformable B. burgdorferi strain. Variants of the mini system were constructed with or without the long inverted repeat (IR) located upstream of vlsE and on both circular and linear plasmids to investigate the importance of the IR and plasmid topology on recombinational switching at vlsE. Amplicon sequencing using PacBio long read technology and analysis of the data with our recently reported pipeline and VAST software showed that the system undergoes switching in mice in both linear and circular versions and that the presence of the hairpin does not seem to be crucial in the linear version, however it is required when the topology is circular.

Antigenic Variation in the Lyme Spirochete: Insights into Recombinational Switching with a Suggested Role for Error-Prone Repair.

The Lyme disease spirochete, Borrelia burgdorferi, uses antigenic variation as a strategy to evade the host's acquired immune response. New variants of surface-localized VlsE are generated efficiently by unidirectional recombination from 15 unexpressed vls cassettes into the vlsE locus. Using algorithms to analyze switching from vlsE sequencing data, we characterize a population of over 45,000 inferred recombination events generated during mouse infection. We present evidence for clustering of these recombination events within the population and along the vlsE gene, a role for the direct repeats flanking the variable region in vlsE, and the importance of sequence homology in determining the location of recombination, despite RecA's dispensability. Finally, we report that non-templated sequence variation is strongly associated with recombinational switching and occurs predominantly at the 5' end of conversion tracts. This likely results from an error-prone repair mechanism operational during recombinational switching that elevates the mutation rate > 5,000-fold in switched regions.

Analysis of recombinational switching at the antigenic variation locus of the Lyme spirochete using a novel PacBio sequencing pipeline.

The Lyme disease spirochete evades the host immune system by combinatorial variation of VlsE, a surface antigen. Antigenic variation occurs via segmental gene conversion from contiguous silent cassettes into the vlsE locus. Because of the high degree of similarity between switch variants and the size of vlsE, short-read NGS technologies have been unsuitable for sequencing vlsE populations. Here we use PacBio sequencing technology coupled with the first fully-automated software pipeline (VAST) to accurately process NGS data by minimizing error frequency, eliminating heteroduplex errors and accurately aligning switch variants. We extend earlier studies by showing use of almost all of the vlsE SNP repertoire. In different tissues of the same mouse, 99.6% of the variants were unique, suggesting that dissemination of Borrelia burgdorferi is predominantly unidirectional with little tissue-to-tissue hematogenous dissemination. We also observed a similar number of variants in SCID and wild-type mice, a heatmap of location and frequency of amino acid changes on the 3D structure and note differences observed in SCID versus wild type mice that hint at possible amino acid function. Our observed selection against diversification of residues at the dimer interface in wild-type mice strongly suggests that dimerization is required for in vivo functionality of vlsE.

Biased Gene Conversion in Rhizobium etli Is Caused by Preferential Double-Strand Breaks on One of the Recombining Homologs.

UNLABELLED: Gene conversion, the nonreciprocal transfer of information during homologous recombination, is the main process that maintains identity between members of multigene families. Gene conversion in the nitrogenase (nifH) multigene family of Rhizobium etli was analyzed by using a two-plasmid system, where each plasmid carried a copy of nifH. One of the nifH copies was modified, creating restriction fragment length polymorphisms (RFLPs) spaced along the gene. Once the modified plasmid was introduced into R. etli, selection was done for cointegration with a resident plasmid lacking the RFLPs. Most of the cointegrate molecules harbor gene conversion events, biased toward a gain of RFLPs. This bias may be explained under the double-strand break repair model by proposing that the nifH gene lacking the RFLPs suffers a DNA double-strand break, so the incoming plasmid functions as a template for repairing the homolog on the resident plasmid. To support this proposal, we cloned an SceI site into the nifH homolog that had the RFLPs used for scoring gene conversion. In vivo expression of the meganuclease I-SceI allowed the generation of a double-strand break on this homolog. Upon introduction of this modified plasmid into an R. etli strain lacking I-SceI, biased gene conversion still favored the retention of markers on the incoming plasmid. In contrast, when the recipient strain ectopically expressed I-SceI, a dramatic reversal in gene conversion bias was seen, favoring the preservation of resident sequences. These results show that biased gene conversion is caused by preferential double-strand breaks on one of the recombining homologs. IMPORTANCE: In this work, we analyzed gene conversion by using a system that entails horizontal gene transfer followed by homologous recombination in the recipient cell. Most gene conversion events are biased toward the acquisition of the incoming sequences, ranging in size from 120 bp to 800 bp. This bias is due to preferential cutting of resident DNA and can be reversed upon introduction of a double-strand break on the incoming sequence. Since conditions used in this work are similar to those in horizontal gene transfer, it provides evidence that, upon transfer, the resident DNA preferentially acquires gene variants.

High polyhydroxybutyrate production in Pseudomonas extremaustralis is associated with differential expression of horizontally acquired and core genome polyhydroxyalkanoate synthase genes.

Pseudomonas extremaustralis produces mainly polyhydroxybutyrate (PHB), a short chain length polyhydroxyalkanoate (sclPHA) infrequently found in Pseudomonas species. Previous studies with this strain demonstrated that PHB genes are located in a genomic island. In this work, the analysis of the genome of P. extremaustralis revealed the presence of another PHB cluster phbFPX, with high similarity to genes belonging to Burkholderiales, and also a cluster, phaC1ZC2D, coding for medium chain length PHA production (mclPHA). All mclPHA genes showed high similarity to genes from Pseudomonas species and interestingly, this cluster also showed a natural insertion of seven ORFs not related to mclPHA metabolism. Besides PHB, P. extremaustralis is able to produce mclPHA although in minor amounts. Complementation analysis demonstrated that both mclPHA synthases, PhaC1 and PhaC2, were functional. RT-qPCR analysis showed different levels of expression for the PHB synthase, phbC, and the mclPHA synthases. The expression level of phbC, was significantly higher than the obtained for phaC1 and phaC2, in late exponential phase cultures. The analysis of the proteins bound to the PHA granules showed the presence of PhbC and PhaC1, whilst PhaC2 could not be detected. In addition, two phasin like proteins (PhbP and PhaI) associated with the production of scl and mcl PHAs, respectively, were detected. The results of this work show the high efficiency of a foreign gene (phbC) in comparison with the mclPHA core genome genes (phaC1 and phaC2) indicating that the ability of P. extremaustralis to produce high amounts of PHB could be explained by the different expression levels of the genes encoding the scl and mcl PHA synthases.

Posttranscriptional regulation of PhbR, the transcriptional activator of polyhydroxybutyrate synthesis, by iron and the sRNA ArrF in Azotobacter vinelandii.

Azotobacter vinelandii is a Gram-negative bacterium able to synthesize poly-ß-hydroxybutyrate (PHB), a biodegradable plastic of industrial interest. The phbBAC operon encodes the enzymes of PHB synthesis and is activated by the transcriptional regulator PhbR and the sigma factor RpoS. Iron limitation has been previously reported to increase PHB accumulation in A. vinelandii; however, the mechanism by which iron controls PHB synthesis is unknown. Under iron starvation in Escherichia coli, the RyhB sRNA modulates the translation of genes involved in iron homeostasis. ArrF is the RyhB analogue in A. vinelandii and similarly increases in quantity during Fe(2+) depletion. In this study, we evaluate the effect of iron and ArrF on PHB accumulation, and on phbR and phbBAC expression in A. vinelandii strain UW136. Using transcriptional and translational fusions of phbR and phbB with gusA reporter gene, we found that iron limitation increased the expression of phbBAC at the transcriptional level and posttranscriptionally increased the expression of phbR. We also found that the ArrF sRNA is a positive regulator of phbR expression at the posttranscriptional level. Collectively, these data suggest that iron limitation increases the translation of phbR through ArrF.

RsmA post-transcriptionally controls PhbR expression and polyhydroxybutyrate biosynthesis in Azotobacter vinelandii.

In Azotobacter vinelandii the two-component GacS/GacA system is required for synthesis of polyhydroxybutyrate (PHB) and of the exopolysaccharide alginate. The RsmA protein was shown to interact with the alginate biosynthetic algD mRNA, acting as a translational repressor, and GacA was found to activate transcription of the rsmZ1 and rsmZ2 genes that encode small RNAs interacting with RsmA to counteract its repressor activity. The phbBAC operon encodes the enzymes of PHB synthesis and is activated by the transcriptional regulator PhbR. This study shows that GacA is required for transcription of one rsmY and seven rsmZ1-rsmZ7 genes present in the A. vinelandii genome, and that inactivation of rsmA results in increased PHB production. Transcriptional and translational phbR-gusA gene fusions were used to show that the gacA mutation negatively affected the expression of the phbR gene at the translational level. We also demonstrated an in vitro interaction of RsmA with RNAs corresponding to phbB and phbR mRNA leaders, and showed that the stability of phbR and phbB mRNAs is increased in the rsmA mutant. Taken together these results indicate that in A. vinelandii, RsmA post-transcriptionally represses the expression of PhbR.

Transcriptional activation of the Azotobacter vinelandii polyhydroxybutyrate biosynthetic genes phbBAC by PhbR and RpoS.

We previously showed that in Azotobacter vinelandii, accumulation of polyhydroxybutyrate (PHB) occurs mainly during the stationary phase, and that a mutation in phbR, encoding a transcriptional regulator of the AraC family, reduces PHB accumulation. In this study, we characterized the roles of PhbR and RpoS, a central regulator during stationary phase in bacteria, in the regulation of expression of the PHB biosynthetic operon phbBAC and phbR. We showed that inactivation of rpoS reduced PHB accumulation, similar to the phbR mutation, and inactivation of both rpoS and phbR resulted in an inability to produce PHB. We carried out expression studies with the wild-type, and the rpoS, phbR and double rpoS-phbR mutant strains, using quantitative RT-PCR, as well as phbB : : gusA and phbR : : gusA gene fusions. These studies showed that both PhbR and RpoS act as activators of phbB and phbR, and revealed a role for PhbR as an autoactivator. We also demonstrated that PhbR binds specifically to two almost identical 18 bp sites, TGTCACCAA-N(4)-CACTA and TGTCACCAA-N(4)-CAGTA, present in the phbB promoter region. The activation of phbB and phbR transcription by RpoS reported here is in agreement with the observation that accumulation of PHB in A. vinelandii occurs mainly during the stationary phase.

The extent of migration of the Holliday junction is a crucial factor for gene conversion in Rhizobium etli.

Gene conversion, defined as the nonreciprocal transfer of DNA, is one result of homologous recombination. Three steps in recombination could give rise to gene conversion: (i) DNA synthesis for repair of the degraded segment, (ii) Holliday junction migration, leading to heteroduplex formation, and (iii) repair of mismatches in the heteroduplex. There are at least three proteins (RuvAB, RecG, and RadA) that participate in the second step. Their roles have been studied for homologous recombination, but evidence of their relative role in gene conversion is lacking. In this work, we showed the effect on gene conversion of mutations in ruvB, recG, and radA in Rhizobium etli, either alone or in combination, using a cointegration strategy previously developed in our laboratory. The results indicate that the RuvAB system is highly efficient for gene conversion, since its absence provokes smaller gene conversion segments than those in the wild type as well as a shift in the preferred position of conversion tracts. The RecG system possesses a dual role for gene conversion. Inactivation of recG leads to longer gene conversion tracts than those in the wild type, indicating that its activity may hinder heteroduplex extension. However, under circumstances where it is the only migration activity present (as in the ruvB radA double mutant), conversion segments can still be seen, indicating that RecG can also promote gene conversion. RadA is the least efficient system in R. etli but is still needed for the production of detectable gene conversion tracts.