MultiMetEval: comparative and multi-objective analysis of genome-scale metabolic models

Comparative metabolic modelling is emerging as a novel field, supported by the development of reliable and standardized approaches for constructing genome-scale metabolic models in high throughput. New software solutions are needed to allow efficient comparative analysis of multiple models in the context of multiple cellular objectives. Here, we present the user-friendly software framework Multi-Metabolic Evaluator (MultiMetEval), built upon SurreyFBA, which allows the user to compose collections of metabolic models that together can be subjected to flux balance analysis. Additionally, MultiMetEval implements functionalities for multi-objective analysis by calculating the Pareto front between two cellular objectives. Using a previously generated dataset of 38 actinobacterial genome-scale metabolic models, we show how these approaches can lead to exciting novel insights. Firstly, after incorporating several pathways for the biosynthesis of natural products into each of these models, comparative flux balance analysis predicted that species like Streptomyces that harbour the highest diversity of secondary metabolite biosynthetic gene clusters in their genomes do not necessarily have the metabolic network topology most suitable for compound overproduction. Secondly, multi-objective analysis of biomass production and natural product biosynthesis in these actinobacteria shows that the well-studied occurrence of discrete metabolic switches during the change of cellular objectives is inherent to their metabolic network architecture. Comparative and multi-objective modelling can lead to insights that could not be obtained by normal flux balance analyses. MultiMetEval provides a powerful platform that makes these analyses straightforward for biologists. Sources and binaries of MultiMetEval are freely available from https://github.com/PiotrZakrzewski/MetEv​al/downloads

Meta-analysis of GWAS of over 16,000 individuals with autism spectrum disorder highlights a novel locus at 10q24.32 and a significant overlap with schizophrenia.

To access publisher's full text version of this article, please click on the hyperlink in Additional Links field or click on the hyperlink at the top of the page marked FilesOver the past decade genome-wide association studies (GWAS) have been applied to aid in the understanding of the biology of traits. The success of this approach is governed by the underlying effect sizes carried by the true risk variants and the corresponding statistical power to observe such effects given the study design and sample size under investigation. Previous ASD GWAS have identified genome-wide significant (GWS) risk loci; however, these studies were of only of low statistical power to identify GWS loci at the lower effect sizes (odds ratio (OR) <1.15).We conducted a large-scale coordinated international collaboration to combine independent genotyping data to improve the statistical power and aid in robust discovery of GWS loci. This study uses genome-wide genotyping data from a discovery sample (7387 ASD cases and 8567 controls) followed by meta-analysis of summary statistics from two replication sets (7783 ASD cases and 11359 controls; and 1369 ASD cases and 137308 controls).We observe a GWS locus at 10q24.32 that overlaps several genes including PITX3, which encodes a transcription factor identified as playing a role in neuronal differentiation and CUEDC2 previously reported to be associated with social skills in an independent population cohort. We also observe overlap with regions previously implicated in schizophrenia which was further supported by a strong genetic correlation between these disorders (Rg = 0.23; P = 9 × 10(-6)). We further combined these Psychiatric Genomics Consortium (PGC) ASD GWAS data with the recent PGC schizophrenia GWAS to identify additional regions which may be important in a common neurodevelopmental phenotype and identified 12 novel GWS loci. These include loci previously implicated in ASD such as FOXP1 at 3p13, ATP2B2 at 3p25.3, and a 'neurodevelopmental hub' on chromosome 8p11.23.This study is an important step in the ongoing endeavour to identify the loci which underpin the common variant signal in ASD. In addition to novel GWS loci, we have identified a significant genetic correlation with schizophrenia and association of ASD with several neurodevelopmental-related genes such as EXT1, ASTN2, MACROD2, and HDAC4.National Institutes of Mental Health (NIMH, USA) ACE Network Autism Genetic Resource Exchange (AGRE) is a program of Autism Speaks (USA) The Autism Genome Project (AGP) from Autism Speaks (USA) Canadian Institutes of Health Research (CIHR), Genome Canada Health Research Board (Ireland) Hilibrand Foundation (USA) Medical Research Council (UK) National Institutes of Health (USA) Ontario Genomics Institute University of Toronto McLaughlin Centre Simons Foundation Johns Hopkins Autism Consortium of Boston NLM Family foundation National Institute of Health grants National Health Medical Research Council Scottish Rite Spunk Fund, Inc. Rebecca and Solomon Baker Fund APEX Foundation National Alliance for Research in Schizophrenia and Affective Disorders (NARSAD) endowment fund of the Nancy Pritzker Laboratory (Stanford) Autism Society of America Janet M. Grace Pervasive Developmental Disorders Fund The Lundbeck Foundation universities and university hospitals of Aarhus and Copenhagen Stanley Foundation Centers for Disease Control and Prevention (CDC) Netherlands Scientific Organization Dutch Brain Foundation VU University Amsterdam Trinity Centre for High Performance Computing through Science Foundation Ireland Autism Genome Project (AGP) from Autism Speak

Polygenic transmission disequilibrium confirms that common and rare variation act additively to create risk for autism spectrum disorders

Autism spectrum disorder (ASD) risk is influenced by common polygenic and de novo variation. We aimed to clarify the influence of polygenic risk for ASD and to identify subgroups of ASD cases, including those with strongly acting de novo variants, in which polygenic risk is relevant. Using a novel approach called the polygenic transmission disequilibrium test and data from 6,454 families with a child with ASD, we show that polygenic risk for ASD, schizophrenia, and greater educational attainment is over-transmitted to children with ASD. These findings hold independent of proband IQ. We find that polygenic variation contributes additively to risk in ASD cases who carry a strongly acting de novo variant. Lastly, we show that elements of polygenic risk are independent and differ in their relationship with phenotype. These results confirm that the genetic influences on ASD are additive and suggest that they create risk through at least partially distinct etiologic pathways

Longitudinal analysis of vertebral fracture and BMD in a Canadian cohort of adult cystic fibrosis patients

<p>Abstract</p> <p>Background</p> <p>Vertebral fractures in patients with cystic fibrosis (CF) may contribute to an accelerated decline in lung function and can be a contraindication to lung transplantation. In this study, we examined longitudinal change in bone mineral density (BMD) and the prevalence of vertebral fractures in adult CF patients, without lung-transplant, attending a Canadian specialty clinic.</p> <p>Methods</p> <p>Retrospective chart review of all patients attending an Adult Cystic Fibrosis Clinic at Hamilton Health Sciences in Hamilton, Canada. Forty-nine of 56 adults met inclusion criteria. Chest radiographs were graded by consensus approach using Genant's semi-quantitative method to identify and grade fractured vertebrae. Dual x-ray absorptiometry (DXA) scans were also reviewed.</p> <p>Results</p> <p>The mean age of the cohort was 25.2 years (SD 9.4), 43% were male. The mean body mass index (BMI) was 19.8 (2.8) for males and 21.7 (5.1) for females. At baseline, the rate of at least one vertebral fracture was 16.3%; rising to 21.3% (prevalent and incident) after a 3-year follow-up. The mean BMD T-or Z-scores at baseline were -0.80 (SD 1.1) at the lumbar spine, -0.57 (SD 0.97) at the proximal femur, and -0.71 (SD 1.1) at the whole body. Over approximately 4-years, the mean percent change in BMD was -1.93% at the proximal femur and -0.73% at the lumbar spine.</p> <p>Conclusion</p> <p>Approximately one in five CF patients demonstrated at least one or more vertebral fractures. Moderate declines in BMD were observed. Given the high rate of vertebral fractures noted in this cohort of adult CF patients, and the negative impact they have on compromised lung functioning, regular screening for vertebral fractures should be considered on routine chest radiographs.</p

Genome sequencing and comparative genomics of the broad host-range pathogen Rhizoctonia solani AG8

Rhizoctonia solani is a soil-borne basidiomycete fungus with a necrotrophic lifestyle which is classified into fourteen reproductively incompatible anastomosis groups (AGs). One of these, AG8, is a devastating pathogen causing bare patch of cereals, brassicas and legumes. R. solani is a multinucleate heterokaryon containing significant heterozygosity within a single cell. This complexity posed significant challenges for the assembly of its genome. We present a high quality genome assembly of R. solani AG8 and a manually curated set of 13,964 genes supported by RNA-seq. The AG8 genome assembly used novel methods to produce a haploid representation of its heterokaryotic state. The whole-genomes of AG8, the rice pathogen AG1-IA and the potato pathogen AG3 were observed to be syntenic and co-linear. Genes and functions putatively relevant to pathogenicity were highlighted by comparing AG8 to known pathogenicity genes, orthology databases spanning 197 phytopathogenic taxa and AG1-IA.We also observed SNP-level “hypermutation” of CpG dinucleotides to TpG between AG8 nuclei, with similarities to repeat-induced point mutation (RIP). Interestingly, gene-coding regions were widely affected along with repetitive DNA, which has not been previously observed for RIP in mononuclear fungi of the Pezizomycotina. The rate of heterozygous SNP mutations within this single isolate of AG8 was observed to be higher than SNP mutation rates observed across populations of most fungal species compared. Comparative analyses were combined to predict biological processes relevant to AG8 and 308 proteins with effector-like characteristics, forming a valuable resource for further study of this pathosystem. Predicted effector-like proteins had elevated levels of non-synonymous point mutations relative to synonymous mutations (dN/dS), suggesting that they may be under diversifying selection pressures. In addition, the distant relationship to sequenced necrotrophs of the Ascomycota suggests the R. solani genome sequence may prove to be a useful resource in future comparative analysis of plant pathogens

XIAP Regulates Cytosol-Specific Innate Immunity to Listeria Infection

The inhibitor of apoptosis protein (IAP) family has been implicated in immune regulation, but the mechanisms by which IAP proteins contribute to immunity are incompletely understood. We show here that X-linked IAP (XIAP) is required for innate immune control of Listeria monocytogenes infection. Mice deficient in XIAP had a higher bacterial burden 48 h after infection than wild-type littermates, and exhibited substantially decreased survival. XIAP enhanced NF-κB activation upon L. monocytogenes infection of activated macrophages, and prolonged phosphorylation of Jun N-terminal kinase (JNK) specifically in response to cytosolic bacteria. Additionally, XIAP promoted maximal production of pro-inflammatory cytokines upon bacterial infection in vitro or in vivo, or in response to combined treatment with NOD2 and TLR2 ligands. Together, our data suggest that XIAP regulates innate immune responses to L. monocytogenes infection by potentiating synergy between Toll-like receptors (TLRs) and Nod-like receptors (NLRs) through activation of JNK- and NF-κB–dependent signaling

Potentiation of Epithelial Innate Host Responses by Intercellular Communication

The epithelium efficiently attracts immune cells upon infection despite the low number of pathogenic microbes and moderate levels of secreted chemokines per cell. Here we examined whether horizontal intercellular communication between cells may contribute to a coordinated response of the epithelium. Listeria monocytogenes infection, transfection, and microinjection of individual cells within a polarized intestinal epithelial cell layer were performed and activation was determined at the single cell level by fluorescence microscopy and flow cytometry. Surprisingly, chemokine production after L. monocytogenes infection was primarily observed in non-infected epithelial cells despite invasion-dependent cell activation. Whereas horizontal communication was independent of gap junction formation, cytokine secretion, ion fluxes, or nitric oxide synthesis, NADPH oxidase (Nox) 4-dependent oxygen radical formation was required and sufficient to induce indirect epithelial cell activation. This is the first report to describe epithelial cell-cell communication in response to innate immune activation. Epithelial communication facilitates a coordinated infectious host defence at the very early stage of microbial infection

Microstructural evolution and trace element mobility in Witwatersrand pyrite

Microstructural analysis of pyrite from a single sample of Witwatersrand conglomerate indicates a complex deformation history involving components of both plastic and brittle deformation. Internal deformation associated with dislocation creep is heterogeneously developed within grains, shows no systematic relationship to bulk rock strain or the location of grain boundaries and is interpreted to represent an episode of pyrite deformation that predates the incorporation of detrital pyrite grains into the Central Rand conglomerates. In contrast, brittle deformation, manifest by grain fragmentation that transects dislocation-related microstructures, is spatially related to grain contacts and is interpreted to represent post-depositional deformation of the Central Rand conglomerates. Analysis of the low-angle boundaries associated with the early dislocation creep phase of deformation indicates the operation of {100} slip systems. However, some orientation boundaries have geometrical characteristics that are not consistent with simple {100} deformation.These boundaries may represent the combination of multiple slip systems or the operation of the previously unrecognized {120} slip system. These boundaries are associated with order of magnitude enrichments in As, Ni and Co that indicate a deformation control on the remobilization of trace elements within pyrite and a potential slip system control on the effectiveness of fast-diffusion pathways. The results confirm the importance of grain-scale elemental remobilization within pyrite prior to their incorporation into the Witwatersrand gold-bearing conglomerates. Since the relationship between gold and pyrite is intimately related to the trace element geochemistry of pyrite, the results have implications for the application of minor element geochemistry to ore deposit formation, suggest a reason for heterogeneous conductivity and localized gold precipitation in natural pyrite and provide a framework for improving mineral processing

Type 1 Fimbriae, a Colonization Factor of Uropathogenic Escherichia coli, Are Controlled by the Metabolic Sensor CRP-cAMP

Type 1 fimbriae are a crucial factor for the virulence of uropathogenic Escherichia coli during the first steps of infection by mediating adhesion to epithelial cells. They are also required for the consequent colonization of the tissues and for invasion of the uroepithelium. Here, we studied the role of the specialized signal transduction system CRP-cAMP in the regulation of type 1 fimbriation. Although initially discovered by regulating carbohydrate metabolism, the CRP-cAMP complex controls a major regulatory network in Gram-negative bacteria, including a broad subset of genes spread into different functional categories of the cell. Our results indicate that CRP-cAMP plays a dual role in type 1 fimbriation, affecting both the phase variation process and fimA promoter activity, with an overall repressive outcome on fimbriation. The dissection of the regulatory pathway let us conclude that CRP-cAMP negatively affects FimB-mediated recombination by an indirect mechanism that requires DNA gyrase activity. Moreover, the underlying studies revealed that CRP-cAMP controls the expression of another global regulator in Gram-negative bacteria, the leucine-responsive protein Lrp. CRP-cAMP-mediated repression is limiting the switch from the non-fimbriated to the fimbriated state. Consistently, a drop in the intracellular concentration of cAMP due to altered physiological conditions (e.g. growth in presence of glucose) increases the percentage of fimbriated cells in the bacterial population. We also provide evidence that the repression of type 1 fimbriae by CRP-cAMP occurs during fast growth conditions (logarithmic phase) and is alleviated during slow growth (stationary phase), which is consistent with an involvement of type 1 fimbriae in the adaptation to stress conditions by promoting biofilm growth or entry into host cells. Our work suggests that the metabolic sensor CRP-cAMP plays a role in coupling the expression of type 1 fimbriae to environmental conditions, thereby also affecting subsequent attachment and colonization of host tissues

Meta-analysis of GWAS of over 16,000 individuals with autism spectrum disorder highlights a novel locus at 10q24.32 and a significant overlap with schizophrenia

Background Over the past decade genome-wide association studies (GWAS) have been applied to aid in the understanding of the biology of traits. The success of this approach is governed by the underlying effect sizes carried by the true risk variants and the corresponding statistical power to observe such effects given the study design and sample size under investigation. Previous ASD GWAS have identified genome-wide significant (GWS) risk loci; however, these studies were of only of low statistical power to identify GWS loci at the lower effect sizes (odds ratio (OR) <1.15). Methods We conducted a large-scale coordinated international collaboration to combine independent genotyping data to improve the statistical power and aid in robust discovery of GWS loci. This study uses genome-wide genotyping data from a discovery sample (7387 ASD cases and 8567 controls) followed by meta-analysis of summary statistics from two replication sets (7783 ASD cases and 11359 controls; and 1369 ASD cases and 137308 controls). Results We observe a GWS locus at 10q24.32 that overlaps several genes including PITX3, which encodes a transcription factor identified as playing a role in neuronal differentiation and CUEDC2 previously reported to be associated with social skills in an independent population cohort. We also observe overlap with regions previously implicated in schizophrenia which was further supported by a strong genetic correlation between these disorders (Rg = 0.23; P=9 ×10−6). We further combined these Psychiatric Genomics Consortium (PGC) ASD GWAS data with the recent PGC schizophrenia GWAS to identify additional regions which may be important in a common neurodevelopmental phenotype and identified 12 novel GWS loci. These include loci previously implicated in ASD such as FOXP1 at 3p13, ATP2B2 at 3p25.3, and a ‘neurodevelopmental hub’ on chromosome 8p11.23. Conclusions This study is an important step in the ongoing endeavour to identify the loci which underpin the common variant signal in ASD. In addition to novel GWS loci, we have identified a significant genetic correlation with schizophrenia and association of ASD with several neurodevelopmental-related genes such as EXT1, ASTN2, MACROD2, and HDAC4