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

Evidence for a Grooming Claw in a North American Adapiform Primate: Implications for Anthropoid Origins

Among fossil primates, the Eocene adapiforms have been suggested as the closest relatives of living anthropoids (monkeys, apes, and humans). Central to this argument is the form of the second pedal digit. Extant strepsirrhines and tarsiers possess a grooming claw on this digit, while most anthropoids have a nail. While controversial, the possible presence of a nail in certain European adapiforms has been considered evidence for anthropoid affinities. Skeletons preserved well enough to test this idea have been lacking for North American adapiforms. Here, we document and quantitatively analyze, for the first time, a dentally associated skeleton of Notharctus tenebrosus from the early Eocene of Wyoming that preserves the complete bones of digit II in semi-articulation. Utilizing twelve shape variables, we compare the distal phalanges of Notharctus tenebrosus to those of extant primates that bear nails (n = 21), tegulae (n = 4), and grooming claws (n = 10), and those of non-primates that bear claws (n = 7). Quantitative analyses demonstrate that Notharctus tenebrosus possessed a grooming claw with a surprisingly well-developed apical tuft on its second pedal digit. The presence of a wide apical tuft on the pedal digit II of Notharctus tenebrosus may reflect intermediate morphology between a typical grooming claw and a nail, which is consistent with the recent hypothesis that loss of a grooming claw occurred in a clade containing adapiforms (e.g. Darwinius masillae) and anthropoids. However, a cladistic analysis including newly documented morphologies and thorough representation of characters acknowledged to have states constituting strepsirrhine, haplorhine, and anthropoid synapomorphies groups Notharctus tenebrosus and Darwinius masillae with extant strepsirrhines rather than haplorhines suggesting that the form of pedal digit II reflects substantial homoplasy during the course of early primate evolution

Summary of the DREAM8 Parameter Estimation Challenge: Toward Parameter Identification for Whole-Cell Models

Whole-cell models that explicitly represent all cellular components at the molecular level have the potential to predict phenotype from genotype. However, even for simple bacteria, whole-cell models will contain thousands of parameters, many of which are poorly characterized or unknown. New algorithms are needed to estimate these parameters and enable researchers to build increasingly comprehensive models. We organized the Dialogue for Reverse Engineering Assessments and Methods (DREAM) 8 Whole-Cell Parameter Estimation Challenge to develop new parameter estimation algorithms for whole-cell models. We asked participants to identify a subset of parameters of a whole-cell model given the model’s structure and in silico “experimental” data. Here we describe the challenge, the best performing methods, and new insights into the identifiability of whole-cell models. We also describe several valuable lessons we learned toward improving future challenges. Going forward, we believe that collaborative efforts supported by inexpensive cloud computing have the potential to solve whole-cell model parameter estimation

Separate loci underlie resistance to root infection and leaf scorch during soybean sudden death syndrome

Soybean [Glycine max (L.) Merr.] cultivars show differences in their resistance to both the leaf scorch and root rot of sudden death syndrome (SDS). The syndrome is caused by root colonization by Fusarium virguliforme (ex. F. solani f. sp. glycines). Root susceptibility combined with reduced leaf scorch resistance has been associated with resistance to Heterodera glycines HG Type 1.3.6.7 (race 14) of the soybean cyst nematode (SCN). In contrast, the rhg1 locus underlying resistance to Hg Type 0 was found clustered with three loci for resistance to SDS leaf scorch and one for root infection. The aims of this study were to compare the inheritance of resistance to leaf scorch and root infection in a population that segregated for resistance to SCN and to identify the underlying quantitative trait loci (QTL). “Hartwig”, a cultivar partially resistant to SDS leaf scorch, F. virguliforme root infection and SCN HG Type 1.3.6.7 was crossed with the partially susceptible cultivar “Flyer”. Ninety-two F5-derived recombinant inbred lines and 144 markers were used for map development. Four QTL found in earlier studies were confirmed. One contributed resistance to leaf scorch on linkage group (LG) C2 (Satt277; P = 0.004, R 2 = 15%). Two on LG G underlay root infection at R8 (Satt038; P = 0.0001 R 2 = 28.1%; Satt115; P = 0.003, R 2 = 12.9%). The marker Satt038 was linked to rhg1 underlying resistance to SCN Hg Type 0. The fourth QTL was on LG D2 underlying resistance to root infection at R6 (Satt574; P = 0.001, R 2 = 10%). That QTL was in an interval previously associated with resistance to both SDS leaf scorch and SCN Hg Type 1.3.6.7. The QTL showed repulsion linkage with resistance to SCN that may explain the relative susceptibility to SDS of some SCN resistant cultivars. One additional QTL was discovered on LG G underlying resistance to SDS leaf scorch measured by disease index (Satt130; P = 0.003, R 2 = 13%). The loci and markers will provide tagged alleles with which to improve the breeding of cultivars combining resistances to SDS leaf scorch, root infection and SCN HG Type 1.3.6.7

Socio-demographic factors associated with pet ownership amongst adolescents from a UK birth cohort

Background: In developed nations, pet ownership is common within families. Both physical and psychological health benefits may result from owning a pet during childhood and adolescence. However, it is difficult to determine whether these benefits are due to pet ownership directly or to factors linked to both pet ownership and health. Previous research found associations between a range of socio-demographic factors and pet ownership in seven-year-old children from a UK cohort. The current study extends this research to adolescence, considering that these factors may be important to consider in future Human-Animal Interaction (HAI) research across childhood.Results:The Avon Longitudinal Study of Parents and Children (ALSPAC) collected pet ownership data prospectively via maternal reports from gestation up to age 10 years old and via self-report retrospectively at age 18 for ages 11(n= 3063) to 18 years old (n= 3098) on cats, dogs, rabbits, rodents, birds, fish, tortoise/turtles and horses. The dataset also contains a wide range of potential confounders, including demographic and socio-economic variables.The ownership of all pet types peaked at age 11 (80%) and then decreased during adolescence, with the exclusion of cats which remained constant (around 30%), and dogs which increased through 11–18 years (26–37%). Logistic regression was used to build multivariable models for ownership of each pet type at age 13 years, and the factors identified in these models were compared to previously published data for 7 year-olds in the same cohort. There was some consistency with predictors reported at age 7. Generally sex, birth order, maternal age, maternal education, number of people in the household, house type, and concurrent ownership of other pets were associated with pet ownership at both 7 and 13 years (the direction of association varied according to pet type).Factors that were no longer associated with adolescent pet ownership included child ethnicity, paternal education,and parental social class.Conclusions:A number of socio-demographic factors are associated with pet ownership in childhood and adolescence and they differ according to the type of pet, and age of child. These factors are potential confounders that must be considered in future HAI studies

Multiple Translocation of the AVR-Pita Effector Gene among Chromosomes of the Rice Blast Fungus Magnaporthe oryzae and Related Species

Magnaporthe oryzae is the causal agent of rice blast disease, a devastating problem worldwide. This fungus has caused breakdown of resistance conferred by newly developed commercial cultivars. To address how the rice blast fungus adapts itself to new resistance genes so quickly, we examined chromosomal locations of AVR-Pita, a subtelomeric gene family corresponding to the Pita resistance gene, in various isolates of M. oryzae (including wheat and millet pathogens) and its related species. We found that AVR-Pita (AVR-Pita1 and AVR-Pita2) is highly variable in its genome location, occurring in chromosomes 1, 3, 4, 5, 6, 7, and supernumerary chromosomes, particularly in rice-infecting isolates. When expressed in M. oryzae, most of the AVR-Pita homologs could elicit Pita-mediated resistance, even those from non-rice isolates. AVR-Pita was flanked by a retrotransposon, which presumably contributed to its multiple translocation across the genome. On the other hand, family member AVR-Pita3, which lacks avirulence activity, was stably located on chromosome 7 in a vast majority of isolates. These results suggest that the diversification in genome location of AVR-Pita in the rice isolates is a consequence of recognition by Pita in rice. We propose a model that the multiple translocation of AVR-Pita may be associated with its frequent loss and recovery mediated by its transfer among individuals in asexual populations. This model implies that the high mobility of AVR-Pita is a key mechanism accounting for the rapid adaptation toward Pita. Dynamic adaptation of some fungal plant pathogens may be achieved by deletion and recovery of avirulence genes using a population as a unit of adaptation

Whole genome identification of Mycobacterium tuberculosis vaccine candidates by comprehensive data mining and bioinformatic analyses

<p>Abstract</p> <p>Background</p> <p><it>Mycobacterium tuberculosis</it>, the causative agent of tuberculosis (TB), infects ~8 million annually culminating in ~2 million deaths. Moreover, about one third of the population is latently infected, 10% of which develop disease during lifetime. Current approved prophylactic TB vaccines (BCG and derivatives thereof) are of variable efficiency in adult protection against pulmonary TB (0%–80%), and directed essentially against early phase infection.</p> <p>Methods</p> <p>A genome-scale dataset was constructed by analyzing published data of: (1) global gene expression studies under conditions which simulate intra-macrophage stress, dormancy, persistence and/or reactivation; (2) cellular and humoral immunity, and vaccine potential. This information was compiled along with revised annotation/bioinformatic characterization of selected gene products and <it>in silico </it>mapping of T-cell epitopes. Protocols for scoring, ranking and prioritization of the antigens were developed and applied.</p> <p>Results</p> <p>Cross-matching of literature and <it>in silico</it>-derived data, in conjunction with the prioritization scheme and biological rationale, allowed for selection of 189 putative vaccine candidates from the entire genome. Within the 189 set, the relative distribution of antigens in 3 functional categories differs significantly from their distribution in the whole genome, with reduction in the Conserved hypothetical category (due to improved annotation) and enrichment in Lipid and in Virulence categories. Other prominent representatives in the 189 set are the PE/PPE proteins; iron sequestration, nitroreductases and proteases, all within the Intermediary metabolism and respiration category; ESX secretion systems, resuscitation promoting factors and lipoproteins, all within the Cell wall category. Application of a ranking scheme based on qualitative and quantitative scores, resulted in a list of 45 best-scoring antigens, of which: 74% belong to the dormancy/reactivation/resuscitation classes; 30% belong to the Cell wall category; 13% are classical vaccine candidates; 9% are categorized Conserved hypotheticals, all potentially very potent T-cell antigens.</p> <p>Conclusion</p> <p>The comprehensive literature and <it>in silico</it>-based analyses allowed for the selection of a repertoire of 189 vaccine candidates, out of the whole-genome 3989 ORF products. This repertoire, which was ranked to generate a list of 45 top-hits antigens, is a platform for selection of genes covering all stages of <it>M. tuberculosis </it>infection, to be incorporated in rBCG or subunit-based vaccines.</p

Antiinflammatory Therapy with Canakinumab for Atherosclerotic Disease

Background: Experimental and clinical data suggest that reducing inflammation without affecting lipid levels may reduce the risk of cardiovascular disease. Yet, the inflammatory hypothesis of atherothrombosis has remained unproved. Methods: We conducted a randomized, double-blind trial of canakinumab, a therapeutic monoclonal antibody targeting interleukin-1β, involving 10,061 patients with previous myocardial infarction and a high-sensitivity C-reactive protein level of 2 mg or more per liter. The trial compared three doses of canakinumab (50 mg, 150 mg, and 300 mg, administered subcutaneously every 3 months) with placebo. The primary efficacy end point was nonfatal myocardial infarction, nonfatal stroke, or cardiovascular death. RESULTS: At 48 months, the median reduction from baseline in the high-sensitivity C-reactive protein level was 26 percentage points greater in the group that received the 50-mg dose of canakinumab, 37 percentage points greater in the 150-mg group, and 41 percentage points greater in the 300-mg group than in the placebo group. Canakinumab did not reduce lipid levels from baseline. At a median follow-up of 3.7 years, the incidence rate for the primary end point was 4.50 events per 100 person-years in the placebo group, 4.11 events per 100 person-years in the 50-mg group, 3.86 events per 100 person-years in the 150-mg group, and 3.90 events per 100 person-years in the 300-mg group. The hazard ratios as compared with placebo were as follows: in the 50-mg group, 0.93 (95% confidence interval [CI], 0.80 to 1.07; P = 0.30); in the 150-mg group, 0.85 (95% CI, 0.74 to 0.98; P = 0.021); and in the 300-mg group, 0.86 (95% CI, 0.75 to 0.99; P = 0.031). The 150-mg dose, but not the other doses, met the prespecified multiplicity-adjusted threshold for statistical significance for the primary end point and the secondary end point that additionally included hospitalization for unstable angina that led to urgent revascularization (hazard ratio vs. placebo, 0.83; 95% CI, 0.73 to 0.95; P = 0.005). Canakinumab was associated with a higher incidence of fatal infection than was placebo. There was no significant difference in all-cause mortality (hazard ratio for all canakinumab doses vs. placebo, 0.94; 95% CI, 0.83 to 1.06; P = 0.31). Conclusions: Antiinflammatory therapy targeting the interleukin-1β innate immunity pathway with canakinumab at a dose of 150 mg every 3 months led to a significantly lower rate of recurrent cardiovascular events than placebo, independent of lipid-level lowering. (Funded by Novartis; CANTOS ClinicalTrials.gov number, NCT01327846.