Predicting functional associations from metabolism using bi-partite network algorithms

<p>Abstract</p> <p>Background</p> <p>Metabolic reconstructions contain detailed information about metabolic enzymes and their reactants and products. These networks can be used to infer functional associations between metabolic enzymes. Many methods are based on the number of metabolites shared by two enzymes, or the shortest path between two enzymes. Metabolite sharing can miss associations between non-consecutive enzymes in a serial pathway, and shortest-path algorithms are sensitive to high-degree metabolites such as water and ATP that create connections between enzymes with little functional similarity.</p> <p>Results</p> <p>We present new, fast methods to infer functional associations in metabolic networks. A local method, the degree-corrected Poisson score, is based only on the metabolites shared by two enzymes, but uses the known metabolite degree distribution. A global method, based on graph diffusion kernels, predicts associations between enzymes that do not share metabolites. Both methods are robust to high-degree metabolites. They out-perform previous methods in predicting shared Gene Ontology (GO) annotations and in predicting experimentally observed synthetic lethal genetic interactions. Including cellular compartment information improves GO annotation predictions but degrades synthetic lethal interaction prediction. These new methods perform nearly as well as computationally demanding methods based on flux balance analysis.</p> <p>Conclusions</p> <p>We present fast, accurate methods to predict functional associations from metabolic networks. Biological significance is demonstrated by identifying enzymes whose strong metabolic correlations are missed by conventional annotations in GO, most often enzymes involved in transport vs. synthesis of the same metabolite or other enzyme pairs that share a metabolite but are separated by conventional pathway boundaries. More generally, the methods described here may be valuable for analyzing other types of networks with long-tailed degree distributions and high-degree hubs.</p

Deciphering the Catalytic Machinery in 30S Ribosome Assembly GTPase YqeH

YqeH, a circularly permuted GTPase (cpGTPase), which is conserved across bacteria and eukaryotes including humans is important for the maturation of small (30S) ribosomal subunit in Bacillus subtilis. Recently, we have shown that it binds 30S in a GTP/GDP dependent fashion. However, the catalytic machinery employed to hydrolyze GTP is not recognized for any of the cpGTPases, including YqeH. This is because they possess a hydrophobic substitution in place of a catalytic glutamine (present in Ras-like GTPases). Such GTPases were categorized as HAS-GTPases and were proposed to follow a catalytic mechanism, different from the Ras-like proteins.MnmE, another HAS-GTPase, but not circularly permuted, utilizes a potassium ion and water mediated interactions to drive GTP hydrolysis. Though the G-domain of MnmE and YqeH share only approximately 25% sequence identity, the conservation of characteristic sequence motifs between them prompted us to probe GTP hydrolysis machinery in YqeH, by employing homology modeling in conjunction with biochemical experiments. Here, we show that YqeH too, uses a potassium ion to drive GTP hydrolysis and stabilize the transition state. However, unlike MnmE, it does not dimerize in the transition state, suggesting alternative ways to stabilize switches I and II. Furthermore, we identify a potential catalytic residue in Asp-57, whose recognition, in the absence of structural information, was non-trivial due to the circular permutation in YqeH. Interestingly, when compared with MnmE, helix alpha2 that presents Asp-57 is relocated towards the N-terminus in YqeH. An analysis of the YqeH homology model, suggests that despite such relocation, Asp-57 may facilitate water mediated catalysis, similarly as the catalytic Glu-282 of MnmE. Indeed, an abolished catalysis by D57I mutant supports this inference.An uncommon means to achieve GTP hydrolysis utilizing a K(+) ion has so far been demonstrated only for MnmE. Here, we show that YqeH also utilizes a similar mechanism. While the catalytic machinery is similar in both, mechanistic differences may arise based on the way they are deployed. It appears that K(+) driven mechanism emerges as an alternative theme to stabilize the transition state and hydrolyze GTP in a subset of GTPases, such as the HAS-GTPases

Search for new phenomena in final states with an energetic jet and large missing transverse momentum in pp collisions at √ s = 8 TeV with the ATLAS detector

Results of a search for new phenomena in final states with an energetic jet and large missing transverse momentum are reported. The search uses 20.3 fb−1 of √ s = 8 TeV data collected in 2012 with the ATLAS detector at the LHC. Events are required to have at least one jet with pT > 120 GeV and no leptons. Nine signal regions are considered with increasing missing transverse momentum requirements between Emiss T > 150 GeV and Emiss T > 700 GeV. Good agreement is observed between the number of events in data and Standard Model expectations. The results are translated into exclusion limits on models with either large extra spatial dimensions, pair production of weakly interacting dark matter candidates, or production of very light gravitinos in a gauge-mediated supersymmetric model. In addition, limits on the production of an invisibly decaying Higgs-like boson leading to similar topologies in the final state are presente

Visualizing variation within Global Pneumococcal Sequence Clusters (GPSCs) and country population snapshots to contextualize pneumococcal isolates

Knowledge of pneumococcal lineages, their geographic distribution and antibiotic resistance patterns, can give insights into global pneumococcal disease. We provide interactive bioinformatic outputs to explore such topics, aiming to increase dissemination of genomic insights to the wider community, without the need for specialist training. We prepared 12 country-specific phylogenetic snapshots, and international phylogenetic snapshots of 73 common Global Pneumococcal Sequence Clusters (GPSCs) previously defined using PopPUNK, and present them in Microreact. Gene presence and absence defined using Roary, and recombination profiles derived from Gubbins are presented in Phandango for each GPSC. Temporal phylogenetic signal was assessed for each GPSC using BactDating. We provide examples of how such resources can be used. In our example use of a country-specific phylogenetic snapshot we determined that serotype 14 was observed in nine unrelated genetic backgrounds in South Africa. The international phylogenetic snapshot of GPSC9, in which most serotype 14 isolates from South Africa were observed, highlights that there were three independent sub-clusters represented by South African serotype 14 isolates. We estimated from the GPSC9-dated tree that the sub-clusters were each established in South Africa during the 1980s. We show how recombination plots allowed the identification of a 20kb recombination spanning the capsular polysaccharide locus within GPSC97. This was consistent with a switch from serotype 6A to 19A estimated to have occured in the 1990s from the GPSC97-dated tree. Plots of gene presence/absence of resistance genes (tet, erm, cat) across the GPSC23 phylogeny were consistent with acquisition of a composite transposon. We estimated from the GPSC23-dated tree that the acquisition occurred between 1953 and 1975. Finally, we demonstrate the assignment of GPSC31 to 17 externally generated pneumococcal serotype 1 assemblies from Utah via Pathogenwatch. Most of the Utah isolates clustered within GPSC31 in a USA-specific clade with the most recent common ancestor estimated between 1958 and 1981. The resources we have provided can be used to explore to data, test hypothesis and generate new hypotheses. The accessible assignment of GPSCs allows others to contextualize their own collections beyond the data presented here

Protection of Mice against Lethal Challenge with 2009 H1N1 Influenza A Virus by 1918-Like and Classical Swine H1N1 Based Vaccines

The recent 2009 pandemic H1N1 virus infection in humans has resulted in nearly 5,000 deaths worldwide. Early epidemiological findings indicated a low level of infection in the older population (>65 years) with the pandemic virus, and a greater susceptibility in people younger than 35 years of age, a phenomenon correlated with the presence of cross-reactive immunity in the older population. It is unclear what virus(es) might be responsible for this apparent cross-protection against the 2009 pandemic H1N1 virus. We describe a mouse lethal challenge model for the 2009 pandemic H1N1 strain, used together with a panel of inactivated H1N1 virus vaccines and hemagglutinin (HA) monoclonal antibodies to dissect the possible humoral antigenic determinants of pre-existing immunity against this virus in the human population. By hemagglutinination inhibition (HI) assays and vaccination/challenge studies, we demonstrate that the 2009 pandemic H1N1 virus is antigenically similar to human H1N1 viruses that circulated from 1918–1943 and to classical swine H1N1 viruses. Antibodies elicited against 1918-like or classical swine H1N1 vaccines completely protect C57B/6 mice from lethal challenge with the influenza A/Netherlands/602/2009 virus isolate. In contrast, contemporary H1N1 vaccines afforded only partial protection. Passive immunization with cross-reactive monoclonal antibodies (mAbs) raised against either 1918 or A/California/04/2009 HA proteins offered full protection from death. Analysis of mAb antibody escape mutants, generated by selection of 2009 H1N1 virus with these mAbs, indicate that antigenic site Sa is one of the conserved cross-protective epitopes. Our findings in mice agree with serological data showing high prevalence of 2009 H1N1 cross-reactive antibodies only in the older population, indicating that prior infection with 1918-like viruses or vaccination against the 1976 swine H1N1 virus in the USA are likely to provide protection against the 2009 pandemic H1N1 virus. This data provides a mechanistic basis for the protection seen in the older population, and emphasizes a rationale for including vaccination of the younger, naïve population. Our results also support the notion that pigs can act as an animal reservoir where influenza virus HAs become antigenically frozen for long periods of time, facilitating the generation of human pandemic viruses

Primula vulgaris (primrose) genome assembly, annotation and gene expression, with comparative genomics on the heterostyly supergene

Primula vulgaris (primrose) exhibits heterostyly: plants produce self-incompatible pin- or thrum-form flowers, with anthers and stigma at reciprocal heights. Darwin concluded that this arrangement promotes insect-mediated cross-pollination; later studies revealed control by a cluster of genes, or supergene, known as the S (Style length) locus. The P. vulgaris S locus is absent from pin plants and hemizygous in thrum plants (thrum-specific); mutation of S locus genes produces self-fertile homostyle flowers with anthers and stigma at equal heights. Here, we present a 411 Mb P. vulgaris genome assembly of a homozygous inbred long homostyle, representing ~87% of the genome. We annotate over 24,000 P. vulgaris genes, and reveal more genes up-regulated in thrum than pin flowers. We show reduced genomic read coverage across the S locus in other Primula species, including P. veris, where we define the conserved structure and expression of the S locus genes in thrum. Further analysis reveals the S locus has elevated repeat content (64%) compared to the wider genome (37%). Our studies suggest conservation of S locus genetic architecture in Primula, and provide a platform for identification and evolutionary analysis of the S locus and downstream targets that regulate heterostyly in diverse heterostylous species

International links between Streptococcus pneumoniae vaccine serotype 4 sequence type (ST) 801 in Northern European shipyard outbreaks of invasive pneumococcal disease

Background: Pneumococcal disease outbreaks of vaccine preventable serotype 4 sequence type (ST)801 in shipyards have been reported in several countries. We aimed to use genomics to establish any international links between them. Methods: Sequence data from ST801-related outbreak isolates from Norway (n = 17), Finland (n = 11) and Northern Ireland (n = 2) were combined with invasive pneumococcal disease surveillance from the respective countries, and ST801-related genomes from an international collection (n = 41 of > 40,000), totalling 106 genomes. Raw data were mapped and recombination excluded before phylogenetic dating. Results: Outbreak isolates were relatively diverse, with up to 100 SNPs (single nucleotide polymorphisms) and a common ancestor estimated around the year 2000. However, 19 Norwegian and Finnish isolates were nearly indistinguishable (0–2 SNPs) with the common ancestor dated around 2017. Conclusion: The total diversity of ST801 within the outbreaks could not be explained by recent transmission alone, suggesting that harsh environmental and associated living conditions reported in the shipyards may facilitate invasion of colonising pneumococci. However, near identical strains in the Norwegian and Finnish outbreaks does suggest that transmission between international shipyards also contributed to those outbreaks. This indicates the need for improved preventative measures in this working population including pneumococcal vaccination

Pneumococcal lineages associated with serotype replacement and antibiotic resistance in childhood invasive pneumococcal disease in the post-PCV13 era: an international whole-genome sequencing study

Background: Invasive pneumococcal disease remains an important health priority owing to increasing disease incidence caused by pneumococci expressing non-vaccine serotypes. We previously defined 621 Global Pneumococcal Sequence Clusters (GPSCs) by analysing 20 027 pneumococcal isolates collected worldwide and from previously published genomic data. In this study, we aimed to investigate the pneumococcal lineages behind the predominant serotypes, the mechanism of serotype replacement in disease, as well as the major pneumococcal lineages contributing to invasive pneumococcal disease in the post-vaccine era and their antibiotic resistant traits. / Methods: We whole-genome sequenced 3233 invasive pneumococcal disease isolates from laboratory-based surveillance programmes in Hong Kong (n=78), Israel (n=701), Malawi (n=226), South Africa (n=1351), The Gambia (n=203), and the USA (n=674). The genomes represented pneumococci from before and after pneumococcal conjugate vaccine (PCV) introductions and were from children younger than 3 years. We identified predominant serotypes by prevalence and their major contributing lineages in each country, and assessed any serotype replacement by comparing the incidence rate between the pre-PCV and PCV periods for Israel, South Africa, and the USA. We defined the status of a lineage as vaccine-type GPSC (≥50% 13-valent PCV [PCV13] serotypes) or non-vaccine-type GPSC (>50% non-PCV13 serotypes) on the basis of its initial serotype composition detected in the earliest vaccine period to measure their individual contribution toward serotype replacement in each country. Major pneumococcal lineages in the PCV period were identified by pooled incidence rate using a random effects model. / Findings: The five most prevalent serotypes in the PCV13 period varied between countries, with only serotypes 5, 12F, 15B/C, 19A, 33F, and 35B/D common to two or more countries. The five most prevalent serotypes in the PCV13 period varied between countries, with only serotypes 5, 12F, 15B/C, 19A, 33F, and 35B/D common to two or more countries. These serotypes were associated with more than one lineage, except for serotype 5 (GPSC8). Serotype replacement was mainly mediated by expansion of non-vaccine serotypes within vaccine-type GPSCs and, to a lesser extent, by increases in non-vaccine-type GPSCs. A globally spreading lineage, GPSC3, expressing invasive serotypes 8 in South Africa and 33F in the USA and Israel, was the most common lineage causing non-vaccine serotype invasive pneumococcal disease in the PCV13 period. We observed that same prevalent non-vaccine serotypes could be associated with distinctive lineages in different countries, which exhibited dissimilar antibiotic resistance profiles. In non-vaccine serotype isolates, we detected significant increases in the prevalence of resistance to penicillin (52 [21%] of 249 vs 169 [29%] of 575, p=0·0016) and erythromycin (three [1%] of 249 vs 65 [11%] of 575, p=0·0031) in the PCV13 period compared with the pre-PCV period. / Interpretation: Globally spreading lineages expressing invasive serotypes have an important role in serotype replacement, and emerging non-vaccine serotypes associated with different pneumococcal lineages in different countries might be explained by local antibiotic-selective pressures. Continued genomic surveillance of the dynamics of the pneumococcal population with increased geographical representation in the post-vaccine period will generate further knowledge for optimising future vaccine design. / Funding: Bill & Melinda Gates Foundation, Wellcome Sanger Institute, and the US Centers for Disease Control

Built Shallow to Maintain Homeostasis and Persistent Infection: Insight into the Transcriptional Regulatory Network of the Gastric Human Pathogen Helicobacter pylori

Transcriptional regulatory networks (TRNs) transduce environmental signals into coordinated output expression of the genome. Accordingly, they are central for the adaptation of bacteria to their living environments and in host–pathogen interactions. Few attempts have been made to describe a TRN for a human pathogen, because even in model organisms, such as Escherichia coli, the analysis is hindered by the large number of transcription factors involved. In light of the paucity of regulators, the gastric human pathogen Helicobacter pylori represents a very appealing system for understanding how bacterial TRNs are wired up to support infection in the host. Herein, we review and analyze the available molecular and “-omic” data in a coherent ensemble, including protein–DNA and protein–protein interactions relevant for transcriptional control of pathogenic responses. The analysis covers ∼80% of the annotated H. pylori regulators, and provides to our knowledge the first in-depth description of a TRN for an important pathogen. The emerging picture indicates a shallow TRN, made of four main modules (origons) that process the physiological responses needed to colonize the gastric niche. Specific network motifs confer distinct transcriptional response dynamics to the TRN, while long regulatory cascades are absent. Rather than having a plethora of specialized regulators, the TRN of H. pylori appears to transduce separate environmental inputs by using different combinations of a small set of regulators

Measurement of VH, H → b b ¯ production as a function of the vector-boson transverse momentum in 13 TeV pp collisions with the ATLAS detector

Cross-sections of associated production of a Higgs boson decaying into bottom-quark pairs and an electroweak gauge boson, W or Z, decaying into leptons are measured as a function of the gauge boson transverse momentum. The measurements are performed in kinematic fiducial volumes defined in the `simplified template cross-section' framework. The results are obtained using 79.8 fb-1 of proton-proton collisions recorded by the ATLAS detector at the Large Hadron Collider at a centre-of-mass energy of 13 TeV. All measurements are found to be in agreement with the Standard Model predictions, and limits are set on the parameters of an effective Lagrangian sensitive to modifications of the Higgs boson couplings to the electroweak gauge bosons