Bioactive glycans in a microbiome-directed food for children with malnutrition.

Evidence is accumulating that perturbed postnatal development of the gut microbiome contributes to childhood malnutrition1-4. Here we analyse biospecimens from a randomized, controlled trial of a microbiome-directed complementary food (MDCF-2) that produced superior rates of weight gain compared with a calorically more dense conventional ready-to-use supplementary food in 12-18-month-old Bangladeshi children with moderate acute malnutrition4. We reconstructed 1,000 bacterial genomes (metagenome-assembled genomes (MAGs)) from the faecal microbiomes of trial participants, identified 75 MAGs of which the abundances were positively associated with ponderal growth (change in weight-for-length Z score (WLZ)), characterized changes in MAG gene expression as a function of treatment type and WLZ response, and quantified carbohydrate structures in MDCF-2 and faeces. The results reveal that two Prevotella copri MAGs that are positively associated with WLZ are the principal contributors to MDCF-2-induced expression of metabolic pathways involved in utilizing the component glycans of MDCF-2. The predicted specificities of carbohydrate-active enzymes expressed by their polysaccharide-utilization loci are correlated with (1) the in vitro growth of Bangladeshi P. copri strains, possessing varying degrees of polysaccharide-utilization loci and genomic conservation with these MAGs, in defined medium containing different purified glycans representative of those in MDCF-2, and (2) the levels of faecal carbohydrate structures in the trial participants. These associations suggest that identifying bioactive glycan structures in MDCFs metabolized by growth-associated bacterial taxa will help to guide recommendations about their use in children with acute malnutrition and enable the development of additional formulations.

Glycidamide-induced hypermutation in yeast single-stranded DNA reveals a ubiquitous clock-like mutational motif in humans.

Mutagens often prefer specific nucleotides or oligonucleotide motifs that can be revealed by studying the hypermutation spectra in single-stranded (ss) DNA. We utilized a yeast model to explore mutagenesis by glycidamide, a simple epoxide formed endogenously in humans from the environmental toxicant acrylamide. Glycidamide caused ssDNA hypermutation in yeast predominantly in cytosines and adenines. The most frequent mutations in adenines occurred in the nAt→nGt trinucleotide motif. Base substitutions A→G in this motif relied on Rev1 translesion polymerase activity. Inactivating Rev1 did not alter the nAt trinucleotide preference, suggesting it may be an intrinsic specificity of the chemical reaction between glycidamide and adenine in the ssDNA. We found this mutational motif enriched in published sequencing data from glycidamide-treated mouse cells and ubiquitous in human cancers. In cancers, this motif was positively correlated with the single base substitution (SBS) smoking-associated SBS4 signature, with the clock-like signatures SBS1, SBS5, and was strongly correlated with smoking history and with age of tumor donors. Clock-like feature of the motif was also revealed in cells of human skin and brain. Given its pervasiveness, we propose that this mutational motif reflects mutagenic lesions to adenines in ssDNA from a potentially broad range of endogenous and exogenous agents.

Predicting Structural Susceptibility of Proteins to Proteolytic Processing.

The importance of 3D protein structure in proteolytic processing is well known. However, despite the plethora of existing methods for predicting proteolytic sites, only a few of them utilize the structural features of potential substrates as predictors. Moreover, to our knowledge, there is currently no method available for predicting the structural susceptibility of protein regions to proteolysis. We developed such a method using data from CutDB, a database that contains experimentally verified proteolytic events. For prediction, we utilized structural features that have been shown to influence proteolysis in earlier studies, such as solvent accessibility, secondary structure, and temperature factor. Additionally, we introduced new structural features, including length of protruded loops and flexibility of protein termini. To maximize the prediction quality of the method, we carefully curated the training set, selected an appropriate machine learning method, and sampled negative examples to determine the optimal positive-to-negative class size ratio. We demonstrated that combining our method with models of protease primary specificity can outperform existing bioinformatics methods for the prediction of proteolytic sites. We also discussed the possibility of utilizing this method for bioinformatics prediction of other post-translational modifications.

Lineage Conversion in Pediatric B-Cell Precursor Acute Leukemia under Blinatumomab Therapy.

We report incidence and deep molecular characteristics of lineage switch in 182 pediatric patients affected by B-cell precursor acute lymphoblastic leukemia (BCP-ALL), who were treated with blinatumomab. We documented six cases of lineage switch that occurred after or during blinatumomab exposure. Therefore, lineage conversion was found in 17.4% of all resistance cases (4/27) and 3.2% of relapses (2/63). Half of patients switched completely from BCP-ALL to CD19-negative acute myeloid leukemia, others retained CD19-positive B-blasts and acquired an additional CD19-negative blast population: myeloid or unclassifiable. Five patients had KMT2A gene rearrangements; one had TCF3::ZNF384 translocation. The presented cases showed consistency of gene rearrangements and fusion transcripts across initially diagnosed leukemia and lineage switch. In two of six patients, the clonal architecture assessed by IG/TR gene rearrangements was stable, while in others, loss of clones or gain of new clones was noted. KMT2A-r patients demonstrated very few additional mutations, while in the TCF3::ZNF384 case, lineage switch was accompanied by a large set of additional mutations. The immunophenotype of an existing leukemia sometimes changes via different mechanisms and with different additional molecular changes. Careful investigation of all BM compartments together with all molecular -minimal residual disease studies can lead to reliable identification of lineage switch.

Basis for substrate recognition and distinction by matrix metalloproteinases.

Genomic sequencing and structural genomics produced a vast amount of sequence and structural data, creating an opportunity for structure-function analysis in silico [Radivojac P, et al. (2013) Nat Methods 10(3):221-227]. Unfortunately, only a few large experimental datasets exist to serve as benchmarks for function-related predictions. Furthermore, currently there are no reliable means to predict the extent of functional similarity among proteins. Here, we quantify structure-function relationships among three phylogenetic branches of the matrix metalloproteinase (MMP) family by comparing their cleavage efficiencies toward an extended set of phage peptide substrates that were selected from ∼ 64 million peptide sequences (i.e., a large unbiased representation of substrate space). The observed second-order rate constants [k(obs)] across the substrate space provide a distance measure of functional similarity among the MMPs. These functional distances directly correlate with MMP phylogenetic distance. There is also a remarkable and near-perfect correlation between the MMP substrate preference and sequence identity of 50-57 discontinuous residues surrounding the catalytic groove. We conclude that these residues represent the specificity-determining positions (SDPs) that allowed for the expansion of MMP proteolytic function during evolution. A transmutation of only a few selected SDPs proximal to the bound substrate peptide, and contributing the most to selectivity among the MMPs, is sufficient to enact a global change in the substrate preference of one MMP to that of another, indicating the potential for the rational and focused redesign of cleavage specificity in MMPs.

Mycobacterial nicotinate mononucleotide adenylyltransferase: structure, mechanism, and implications for drug discovery.

Nicotinate mononucleotide adenylyltransferase NadD is an essential enzyme in the biosynthesis of the NAD cofactor, which has been implicated as a target for developing new antimycobacterial therapies. Here we report the crystal structure of Mycobacterium tuberculosis NadD (MtNadD) at a resolution of 2.4 Å. A remarkable new feature of the MtNadD structure, compared with other members of this enzyme family, is a 310 helix that locks the active site in an over-closed conformation. As a result, MtNadD is rendered inactive as it is topologically incompatible with substrate binding and catalysis. Directed mutagenesis was also used to further dissect the structural elements that contribute to the interactions of the two MtNadD substrates, i.e. ATP and nicotinic acid mononucleotide (NaMN). For inhibitory profiling of partially active mutants and wild type MtNadD, we used a small molecule inhibitor of MtNadD with moderate affinity (Ki ∼ 25 µM) and antimycobacterial activity (MIC80) ∼ 40-80 µM). This analysis revealed interferences with some of the residues in the NaMN binding subsite consistent with the competitive inhibition observed for the NaMN substrate (but not ATP). A detailed steady-state kinetic analysis of MtNadD suggests that ATP must first bind to allow efficient NaMN binding and catalysis. This sequential mechanism is consistent with the requirement of transition to catalytically competent (open) conformation hypothesized from structural modeling. A possible physiological significance of this mechanism is to enable the down-regulation of NAD synthesis under ATP-limiting dormancy conditions. These findings point to a possible new strategy for designing inhibitors that lock the enzyme in the inactive over-closed conformation.

Comparative immunogenicity and structural analysis of epitopes of different bacterial L-asparaginases.

BACKGROUND: E.coli type II L-asparaginase is widely used for treatment of acute lymphoblastic leukemia. However, serious side effects such as allergic or hypersensitivity reactions are common for L-asparaginase treatment. Methods for minimizing immune response on L-asparaginase treatment in human include bioengeneering of less immunogenic version of the enzyme or utilizing the homologous enzymes of different origin. To rationalize these approaches we compared immunogenicity of L-asparaginases from five bacterial organisms and performed sequence-structure analysis of the presumable epitope regions. METHODS: IgG and IgM immune response in C57B16 mice after immunization with Wollinella succinogenes type II (WsA), Yersinia pseudotuberculosis type II (YpA), Erwinia carotovora type II (EwA), and Rhodospirillum rubrum type I (RrA) and Escherichia coli type II (EcA) L-asparaginases was evaluated using standard ELISA method. The comparative bioinformatics analysis of structure and sequence of the bacterial L-asparaginases presumable epitope regions was performed. RESULTS: We showed different immunogenic properties of five studied L-asparaginases and confirmed the possibility of replacement of EcA with L-asparaginase from different origin as a second-line treatment. Studied L-asparaginases might be placed in the following order based on the immunogenicity level: YpA > RrA, WsA ≥ EwA > EcA. Most significant cross-immunogenicity was shown between EcA and YpA. We propose that a long N-terminus of YpA enzyme enriched with charged aminoacids and tryptophan could be a reason of higher immunogenicity of YpA in comparison with other considered enzymes. Although the recognized structural and sequence differences in putative epitope regions among five considered L-asparaginases does not fully explain experimental observation of the immunogenicity of the enzymes, the performed analysis set the foundation for further research in this direction. CONCLUSIONS: The performed studies showed different immunogenic properties of L-asparaginases and confirmed the possibility of replacement of EcA with L-asparaginase from different origin. The preferable enzymes for the second line treatment are WsA, RrA, or EwA.

Sequence and conformational specificity in substrate recognition: several human Kunitz protease inhibitor domains are specific substrates of mesotrypsin.

Mesotrypsin is an isoform of trypsin that is uniquely resistant to polypeptide trypsin inhibitors and can cleave some inhibitors rapidly. Previous studies have shown that the amyloid precursor protein Kunitz protease inhibitor domain (APPI) is a specific substrate of mesotrypsin and that stabilization of the APPI cleavage site in a canonical conformation contributes to recognition by mesotrypsin. We hypothesized that other proteins possessing potential cleavage sites stabilized in a similar conformation might also be mesotrypsin substrates. Here we evaluated a series of candidate substrates, including human Kunitz protease inhibitor domains from amyloid precursor-like protein 2 (APLP2), bikunin, hepatocyte growth factor activator inhibitor type 2 (HAI2), tissue factor pathway inhibitor-1 (TFPI1), and tissue factor pathway inhibitor-2 (TFPI2), as well as E-selectin, an unrelated protein possessing a potential cleavage site displaying canonical conformation. We find that Kunitz domains within APLP2, bikunin, and HAI2 are cleaved by mesotrypsin with kinetic profiles of specific substrates. TFPI1 and TFPI2 Kunitz domains are cleaved less efficiently by mesotrypsin, and E-selectin is not cleaved at the anticipated site. Cocrystal structures of mesotrypsin with HAI2 and bikunin Kunitz domains reveal the mode of mesotrypsin interaction with its canonical substrates. Our data suggest that major determinants of mesotrypsin substrate specificity include sequence preferences at the P1 and P'2 positions along with conformational stabilization of the cleavage site in the canonical conformation. Mesotrypsin up-regulation has been implicated previously in cancer progression, and proteolytic clearance of Kunitz protease inhibitors offers potential mechanisms by which mesotrypsin may mediate pathological effects in cancer.

Horizontal gene transfer and genome evolution in Methanosarcina.

BACKGROUND: Genomes of Methanosarcina spp. are among the largest archaeal genomes. One suggested reason for that is massive horizontal gene transfer (HGT) from bacteria. Genes of bacterial origin may be involved in the central metabolism and solute transport, in particular sugar synthesis, sulfur metabolism, phosphate metabolism, DNA repair, transport of small molecules etc. Horizontally transferred (HT) genes are considered to play the key role in the ability of Methanosarcina spp. to inhabit diverse environments. At the moment, genomes of three Methanosarcina spp. have been sequenced, and while these genomes vary in length and number of protein-coding genes, they all have been shown to accumulate HT genes. However, previous estimates had been made when fewer archaeal genomes were known. Moreover, several Methanosarcinaceae genomes from other genera have been sequenced recently. Here, we revise the census of genes of bacterial origin in Methanosarcinaceae. RESULTS: About 5% of Methanosarcina genes have been shown to be horizontally transferred from various bacterial groups to the last common ancestor either of Methanosarcinaceae, or Methanosarcina, or later in the evolution. Simulation of the composition of the NCBI protein non-redundant database for different years demonstrates that the estimates of the HGT rate have decreased drastically since 2002, the year of publication of the first Methanosarcina genome. The phylogenetic distribution of HT gene donors is non-uniform. Most HT genes were transferred from Firmicutes and Proteobacteria, while no HGT events from Actinobacteria to the common ancestor of Methanosarcinaceae were found. About 50% of HT genes are involved in metabolism. Horizontal transfer of transcription factors is not common, while 46% of horizontally transferred genes have demonstrated differential expression in a variety of conditions. HGT of complete operons is relatively infrequent and half of HT genes do not belong to operons. CONCLUSIONS: While genes of bacterial origin are still more frequent in Methanosarcinaceae than in other Archaea, most HGT events described earlier as Methanosarcina-specific seem to have occurred before the divergence of Methanosarcinaceae. Genes horizontally transferred from bacteria to archaea neither tend to be transferred with their regulators, nor in long operons.

Functional diversification of ROK-family transcriptional regulators of sugar catabolism in the Thermotogae phylum.

Large and functionally heterogeneous families of transcription factors have complex evolutionary histories. What shapes specificities toward effectors and DNA sites in paralogous regulators is a fundamental question in biology. Bacteria from the deep-branching lineage Thermotogae possess multiple paralogs of the repressor, open reading frame, kinase (ROK) family regulators that are characterized by carbohydrate-sensing domains shared with sugar kinases. We applied an integrated genomic approach to study functions and specificities of regulators from this family. A comparative analysis of 11 Thermotogae genomes revealed novel mechanisms of transcriptional regulation of the sugar utilization networks, DNA-binding motifs and specific functions. Reconstructed regulons for seven groups of ROK regulators were validated by DNA-binding assays using purified recombinant proteins from the model bacterium Thermotoga maritima. All tested regulators demonstrated specific binding to their predicted cognate DNA sites, and this binding was inhibited by specific effectors, mono- or disaccharides from their respective sugar catabolic pathways. By comparing ligand-binding domains of regulators with structurally characterized kinases from the ROK family, we elucidated signature amino acid residues determining sugar-ligand regulator specificity. Observed correlations between signature residues and the sugar-ligand specificities provide the framework for structure functional classification of the entire ROK family.

Sequence-derived structural features driving proteolytic processing.

Proteolytic signaling, or regulated proteolysis, is an essential part of many important pathways such as Notch, Wnt, and Hedgehog. How the structure of the cleaved substrate regions influences the efficacy of proteolytic processing remains underexplored. Here, we analyzed the relative importance in proteolysis of various structural features derived from substrate sequences using a dataset of more than 5000 experimentally verified proteolytic events captured in CutDB. Accessibility to the solvent was recognized as an essential property of a proteolytically processed polypeptide chain. Proteolytic events were found nearly uniformly distributed among three types of secondary structure, although with some enrichment in loops. Cleavages in α-helices were found to be relatively abundant in regions apparently prone to unfolding, while cleavages in ß-structures tended to be located at the periphery of ß-sheets. Application of the same statistical procedures to proteolytic events divided into separate sets according to the catalytic classes of proteases proved consistency of the results and confirmed that the structural mechanisms of proteolysis are universal. The estimated prediction power of sequence-derived structural features, which turned out to be sufficiently high, presents a rationale for their use in bioinformatic prediction of proteolytic events.

Genome-wide bioinformatics analysis of human protease capacity for proteolytic cleavage of the SARS-CoV-2 spike glycoprotein.

Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) primarily enters the cell by binding the virus's spike (S) glycoprotein to the angiotensin-converting enzyme 2 receptor on the cell surface, followed by proteolytic cleavage by host proteases. Studies have identified furin and transmembrane protease serine 2 proteases in priming and triggering cleavages of the S glycoprotein, converting it into a fusion-competent form and initiating membrane fusion, respectively. Alternatively, SARS-CoV-2 can enter the cell through the endocytic pathway, where activation is triggered by lysosomal cathepsin L. However, other proteases are also suspected to be involved in both entry routes. In this study, we conducted a genome-wide bioinformatics analysis to explore the capacity of human proteases in hydrolyzing peptide bonds of the S glycoprotein. Predictive models of sequence specificity for 169 human proteases were constructed and applied to the S glycoprotein together with the method for predicting structural susceptibility to proteolysis of protein regions. After validating our approach on extensively studied S2' and S1/S2 cleavage sites, we applied our method to each peptide bond of the S glycoprotein across all 169 proteases. Our results indicate that various members of the proprotein convertase subtilisin/kexin type, type II transmembrane family serine protease, and kallikrein families, as well as specific coagulation factors, are capable of cleaving S2' or S1/S2 sites. We have also identified a potential cleavage site of cathepsin L at the K790 position within the S2' loop. Structural analysis suggests that cleavage of this site induces conformational changes similar to the cleavage at the R815 (S2') position, leading to the exposure of the fusion peptide and subsequent fusion with the membrane. Other potential cleavage sites and the influence of mutations in common SARS-CoV-2 variants on proteolytic efficiency are discussed.IMPORTANCEThe entry of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) into the cell, activated by host proteases, is considerably more complex in coronaviruses than in most other viruses and is not fully understood. There is evidence that other proteases beyond the known furin and transmembrane protease serine 2 can activate the spike protein. Another example of uncertainty is the cleavage site for the alternative endocytic route of SARS-CoV-2 entrance, which is still unknown. Bioinformatics methods, modeling protease specificity and estimating the structural susceptibility of protein regions to proteolysis, can aid in studying this topic by predicting the involved proteases and their cleavage sites, thereby substantially reducing the amount of experimental work. Elucidating the mechanisms of spike protein activation is crucial for preventing possible future coronavirus pandemics and developing antiviral drugs.

Evolution of pan-genomes of Escherichia coli, Shigella spp., and Salmonella enterica.

Multiple sequencing of genomes belonging to a bacterial species allows one to analyze and compare statistics and dynamics of the gene complements of species, their pan-genomes. Here, we analyzed multiple genomes of Escherichia coli, Shigella spp., and Salmonella enterica. We demonstrate that the distribution of the number of genomes harboring a gene is well approximated by a sum of two power functions, describing frequent genes (present in many strains) and rare genes (present in few strains). The virtual absence of Shigella-specific genes not present in E. coli genomes confirms previous observations that Shigella is not an independent genus. While the pan-genome size is increasing with each new strain, the number of genes present in a fixed fraction of strains stabilizes quickly. For instance, slightly fewer than 4,000 genes are present in at least half of any group of E. coli genomes. Comparison of S. enterica and E. coli pan-genomes revealed the existence of a common periphery, that is, genes present in some but not all strains of both species. Analysis of phylogenetic trees demonstrates that rare genes from the periphery likely evolve under horizontal transfer, whereas frequent periphery genes may have been inherited from the periphery genome of the common ancestor.

Genomic reconstruction of the transcriptional regulatory network in Bacillus subtilis.

The adaptation of microorganisms to their environment is controlled by complex transcriptional regulatory networks (TRNs), which are still only partially understood even for model species. Genome scale annotation of regulatory features of genes and TRN reconstruction are challenging tasks of microbial genomics. We used the knowledge-driven comparative-genomics approach implemented in the RegPredict Web server to infer TRN in the model Gram-positive bacterium Bacillus subtilis and 10 related Bacillales species. For transcription factor (TF) regulons, we combined the available information from the DBTBS database and the literature with bioinformatics tools, allowing inference of TF binding sites (TFBSs), comparative analysis of the genomic context of predicted TFBSs, functional assignment of target genes, and effector prediction. For RNA regulons, we used known RNA regulatory motifs collected in the Rfam database to scan genomes and analyze the genomic context of new RNA sites. The inferred TRN in B. subtilis comprises regulons for 129 TFs and 24 regulatory RNA families. First, we analyzed 66 TF regulons with previously known TFBSs in B. subtilis and projected them to other Bacillales genomes, resulting in refinement of TFBS motifs and identification of novel regulon members. Second, we inferred motifs and described regulons for 28 experimentally studied TFs with previously unknown TFBSs. Third, we discovered novel motifs and reconstructed regulons for 36 previously uncharacterized TFs. The inferred collection of regulons is available in the RegPrecise database (http://regprecise.lbl.gov/) and can be used in genetic experiments, metabolic modeling, and evolutionary analysis.

Genomic reconstruction of transcriptional regulatory networks in lactic acid bacteria.

BACKGROUND: Genome scale annotation of regulatory interactions and reconstruction of regulatory networks are the crucial problems in bacterial genomics. The Lactobacillales order of bacteria collates various microorganisms having a large economic impact, including both human and animal pathogens and strains used in the food industry. Nonetheless, no systematic genome-wide analysis of transcriptional regulation has been previously made for this taxonomic group. RESULTS: A comparative genomics approach was used for reconstruction of transcriptional regulatory networks in 30 selected genomes of lactic acid bacteria. The inferred networks comprise regulons for 102 orthologous transcription factors (TFs), including 47 novel regulons for previously uncharacterized TFs. Numerous differences between regulatory networks of the Streptococcaceae and Lactobacillaceae groups were described on several levels. The two groups are characterized by substantially different sets of TFs encoded in their genomes. Content of the inferred regulons and structure of their cognate TF binding motifs differ for many orthologous TFs between the two groups. Multiple cases of non-orthologous displacements of TFs that control specific metabolic pathways were reported. CONCLUSIONS: The reconstructed regulatory networks substantially expand the existing knowledge of transcriptional regulation in lactic acid bacteria. In each of 30 studied genomes the obtained regulatory network contains on average 36 TFs and 250 target genes that are mostly involved in carbohydrate metabolism, stress response, metal homeostasis and amino acids biosynthesis. The inferred networks can be used for genetic experiments, functional annotations of genes, metabolic reconstruction and evolutionary analysis. All reconstructed regulons are captured within the Streptococcaceae and Lactobacillaceae collections in the RegPrecise database (http://regprecise.lbl.gov).

Comparative genomics of metabolic capacities of regulons controlled by cis-regulatory RNA motifs in bacteria.

BACKGROUND: In silico comparative genomics approaches have been efficiently used for functional prediction and reconstruction of metabolic and regulatory networks. Riboswitches are metabolite-sensing structures often found in bacterial mRNA leaders controlling gene expression on transcriptional or translational levels.An increasing number of riboswitches and other cis-regulatory RNAs have been recently classified into numerous RNA families in the Rfam database. High conservation of these RNA motifs provides a unique advantage for their genomic identification and comparative analysis. RESULTS: A comparative genomics approach implemented in the RegPredict tool was used for reconstruction and functional annotation of regulons controlled by RNAs from 43 Rfam families in diverse taxonomic groups of Bacteria. The inferred regulons include ~5200 cis-regulatory RNAs and more than 12000 target genes in 255 microbial genomes. All predicted RNA-regulated genes were classified into specific and overall functional categories. Analysis of taxonomic distribution of these categories allowed us to establish major functional preferences for each analyzed cis-regulatory RNA motif family. Overall, most RNA motif regulons showed predictable functional content in accordance with their experimentally established effector ligands. Our results suggest that some RNA motifs (including thiamin pyrophosphate and cobalamin riboswitches that control the cofactor metabolism) are widespread and likely originated from the last common ancestor of all bacteria. However, many more analyzed RNA motifs are restricted to a narrow taxonomic group of bacteria and likely represent more recent evolutionary innovations. CONCLUSIONS: The reconstructed regulatory networks for major known RNA motifs substantially expand the existing knowledge of transcriptional regulation in bacteria. The inferred regulons can be used for genetic experiments, functional annotations of genes, metabolic reconstruction and evolutionary analysis. The obtained genome-wide collection of reference RNA motif regulons is available in the RegPrecise database (http://regprecise.lbl.gov/).

RegPrecise 3.0--a resource for genome-scale exploration of transcriptional regulation in bacteria.

BACKGROUND: Genome-scale prediction of gene regulation and reconstruction of transcriptional regulatory networks in prokaryotes is one of the critical tasks of modern genomics. Bacteria from different taxonomic groups, whose lifestyles and natural environments are substantially different, possess highly diverged transcriptional regulatory networks. The comparative genomics approaches are useful for in silico reconstruction of bacterial regulons and networks operated by both transcription factors (TFs) and RNA regulatory elements (riboswitches). DESCRIPTION: RegPrecise (http://regprecise.lbl.gov) is a web resource for collection, visualization and analysis of transcriptional regulons reconstructed by comparative genomics. We significantly expanded a reference collection of manually curated regulons we introduced earlier. RegPrecise 3.0 provides access to inferred regulatory interactions organized by phylogenetic, structural and functional properties. Taxonomy-specific collections include 781 TF regulogs inferred in more than 160 genomes representing 14 taxonomic groups of Bacteria. TF-specific collections include regulogs for a selected subset of 40 TFs reconstructed across more than 30 taxonomic lineages. Novel collections of regulons operated by RNA regulatory elements (riboswitches) include near 400 regulogs inferred in 24 bacterial lineages. RegPrecise 3.0 provides four classifications of the reference regulons implemented as controlled vocabularies: 55 TF protein families; 43 RNA motif families; ~150 biological processes or metabolic pathways; and ~200 effectors or environmental signals. Genome-wide visualization of regulatory networks and metabolic pathways covered by the reference regulons are available for all studied genomes. A separate section of RegPrecise 3.0 contains draft regulatory networks in 640 genomes obtained by an conservative propagation of the reference regulons to closely related genomes. CONCLUSIONS: RegPrecise 3.0 gives access to the transcriptional regulons reconstructed in bacterial genomes. Analytical capabilities include exploration of: regulon content, structure and function; TF binding site motifs; conservation and variations in genome-wide regulatory networks across all taxonomic groups of Bacteria. RegPrecise 3.0 was selected as a core resource on transcriptional regulation of the Department of Energy Systems Biology Knowledgebase, an emerging software and data environment designed to enable researchers to collaboratively generate, test and share new hypotheses about gene and protein functions, perform large-scale analyses, and model interactions in microbes, plants, and their communities.

Molecular Heterogeneity of Pediatric AML with Atypical Promyelocytes Accumulation in Children-A Single Center Experience.

Acute promyelocytic leukemia (APL) pathogenesis is based on RARA gene translocations, which are of high importance in the diagnosis of and proper therapy selection for APL. However, in some cases acute myeloid leukemia (AML) demonstrates APL-like morphological features such as atypical promyelocytes accumulation. This type of AML is characterized by the involvement of other RAR family members or completely different genes. In the present study, we used conventional karyotyping, FISH and high-throughput sequencing in a group of 271 de novo AML with atypical promyelocytes accumulation. Of those, 255 cases were shown to carry a typical chromosomal translocation t(15;17)(q24;q21) with PML::RARA chimeric gene formation (94.1%). Other RARA-positive cases exhibited cryptic PML::RARA fusion without t(15;17)(q24;q21) (1.8%, n = 5) and variant t(5;17)(q35;q21) translocation with NPM1::RARA chimeric gene formation (1.5%, n = 4). However, 7 RARA-negative AMLs with atypical promyelocytes accumulation were also discovered. These cases exhibited TBL1XR1::RARB and KMT2A::SEPT6 fusions as well as mutations, e.g., NPM1 insertion and non-recurrent chromosomal aberrations. Our findings demonstrate the genetic diversity of AML with APL-like morphological features, which is of high importance for successful therapy implementation.

Bioactive glycans in a microbiome-directed food for malnourished children.

Evidence is accumulating that perturbed postnatal development of the gut microbiome contributes to childhood malnutrition1-4. Designing effective microbiome-directed therapeutic foods to repair these perturbations requires knowledge about how food components interact with the microbiome to alter its expressed functions. Here we use biospecimens from a randomized, controlled trial of a microbiome-directed complementary food prototype (MDCF-2) that produced superior rates of weight gain compared to a conventional ready-to-use supplementary food (RUSF) in 12-18-month-old Bangladeshi children with moderate acute malnutrition (MAM)4. We reconstructed 1000 bacterial genomes (metagenome-assembled genomes, MAGs) present in their fecal microbiomes, identified 75 whose abundances were positively associated with weight gain (change in weight-for-length Z score, WLZ), characterized gene expression changes in these MAGs as a function of treatment type and WLZ response, and used mass spectrometry to quantify carbohydrate structures in MDCF-2 and feces. The results reveal treatment-induced changes in expression of carbohydrate metabolic pathways in WLZ-associated MAGs. Comparing participants consuming MDCF-2 versus RUSF, and MDCF-2-treated children in the upper versus lower quartiles of WLZ responses revealed that two Prevotella copri MAGs positively associated with WLZ were principal contributors to MDCF-2-induced expression of metabolic pathways involved in utilization of its component glycans. Moreover, the predicted specificities of carbohydrate active enzymes expressed by polysaccharide utilization loci (PULs) in these two MAGs correlate with the (i) in vitro growth of Bangladeshi P. copri strains, possessing differing degrees of PUL and overall genomic content similarity to these MAGs, cultured in defined medium containing different purified glycans representative of those in MDCF-2, and (ii) levels of carbohydrate structures identified in feces from clinical trial participants. In the accompanying paper5, we use a gnotobiotic mouse model colonized with age- and WLZ-associated bacterial taxa cultured from this study population, and fed diets resembling those consumed by study participants, to directly test the relationship between P. copri, MDCF-2 glycan metabolism, host ponderal growth responses, and intestinal gene expression and metabolism. The ability to identify bioactive glycan structures in MDCFs that are metabolized by growth-associated bacterial taxa will help guide recommendations about use of this MDCF for children with acute malnutrition representing different geographic locales and ages, as well as enable development of bioequivalent, or more efficacious, formulations composed of culturally acceptable and affordable ingredients.

Evolution of transcriptional regulation in closely related bacteria.

BACKGROUND: The exponential growth of the number of fully sequenced genomes at varying taxonomic closeness allows one to characterize transcriptional regulation using comparative-genomics analysis instead of time-consuming experimental methods. A transcriptional regulatory unit consists of a transcription factor, its binding site and a regulated gene. These units constitute a graph which contains so-called "network motifs", subgraphs of a given structure. Here we consider genomes of closely related Enterobacteriales and estimate the fraction of conserved network motifs and sites as well as positions under selection in various types of non-coding regions. RESULTS: Using a newly developed technique, we found that the highest fraction of positions under selection, approximately 50%, was observed in synvergon spacers (between consecutive genes from the same strand), followed by ~45% in divergon spacers (common 5'-regions), and ~10% in convergon spacers (common 3'-regions). The fraction of selected positions in functional regions was higher, 60% in transcription factor-binding sites and ~45% in terminators and promoters. Small, but significant differences were observed between Escherichia coli and Salmonella enterica. This fraction is similar to the one observed in eukaryotes.The conservation of binding sites demonstrated some differences between types of regulatory units. In E. coli, strains the interactions of the type "local transcriptional factor gene" turned out to be more conserved in feed-forward loops (FFLs) compared to non-motif interactions. The coherent FFLs tend to be less conserved than the incoherent FFLs. A natural explanation is that the former imply functional redundancy. CONCLUSIONS: A naïve hypothesis that FFL would be highly conserved turned out to be not entirely true: its conservation depends on its status in the transcriptional network and also from its usage. The fraction of positions under selection in intergenic regions of bacterial genomes is roughly similar to that of eukaryotes. Known regulatory sites explain 20±5% of selected positions.