Integrative Modelling of Signalling Network Dynamics Identifies Cell Type-selective Therapeutic Strategies for FGFR4-driven Cancers Oncogenic FGFR4 signalling represents a potential therapeutic target in many cancer types, including triple negative breast cancer (TNBC) and hepatocellular carcinoma (HCC). However, resistance to single-agent therapy directed at FGFR4 remains a major challenge, prompting the need to identify more effective combinatorial therapeutic strategies. Here, we integrated computational network modelling and experimental validation to characterise dynamic reprogramming of the FGFR4 signalling network in TNBC following FGFR4 kinase inhibition. We found that AKT, which signals downstream of FGFR4, displayed a rapid and potent reactivation following FGFR4 targeting. Through model-based simulation and systematic prediction of the effect of co-targeting specific network nodes, we predicted, and validated experimentally, strong synergism of co-targeting FGFR4 and particular ErbB kinases or AKT, but not the upstream kinase PI3K. Further, incorporation of protein expression data from hundreds of cancer cell types enabled us to adapt our model to other diverse cellular contexts, leading to the prediction that while AKT rebound occurs frequently, it is not a general phenomenon. Instead, ERK is reactivated in a subset of cell types, including the FGFR4-driven HCC cell line Hep3B. This was subsequently corroborated, and moreover, co-targeting FGFR4 and MEK in Hep3B cells markedly enhanced inhibition of cell proliferation. Overall, these findings provide novel insights into the dynamics of drug-induced network remodelling in cancer cells, highlight the impact of protein expression heterogeneity on network response to targeted therapy and identify candidate cell type-selective combination treatments for FGFR4-driven cancer.
A flux-based machine learning model to simulate the impact of pathogen metabolic heterogeneity on drug interactions AO_SCPLOWBSTRACTC_SCPLOWDrug combinations are a promising strategy to counter antibiotic resistance. However, current experimental and computational approaches do not account for the entire complexity involved in combination therapy design, such as the effect of pathogen metabolic heterogeneity, changes in the growth environment, drug treatment order and time interval. To address these limitations, we present a comprehensive approach that uses genome-scale metabolic modeling and machine learning to guide combination therapy design. Our mechanistic approach (a) accommodates diverse data types, (b) accounts for time- and order-specific interactions, and (c) accurately predicts drug interactions in various growth conditions and their robustness to pathogen metabolic heterogeneity. Our approach achieved high accuracy (AUC = 0.83 for synergy, AUC = 0.98 for antagonism) in predicting drug interactions for E. coli cultured in 57 metabolic conditions based on experimental validation. The entropy in bacterial metabolic response was predictive of combination therapy outcomes across time scales and growth conditions. Simulation of metabolic heterogeneity using population-FBA identified two sub-populations of E. coli cells defined by regulation of cis- to trans-isomerization of unsaturated fatty acids with distinct effects on drug interactions. Analysis of the vast landscape of condition-specific drug interactions revealed a set of 24 robustly synergistic drug combinations with potential for clinical use.

SO_SCPLOWIGNIFICANCEC_SCPLOWWorldwide, 700,000 people die each year from drug-resistant infections. Drug combinations have great potential to reduce the spread of drug-resistant bacteria. However, their potency is impacted by both the pathogen growth environment and the heterogeneity in pathogen metabolism. The metabolic heterogeneity in a pathogen population allows them to survive antibiotic treatment. Here we present a flexible machine-learning framework that utilizes diverse data types to effectively search through the large design space of both sequential and simultaneous combination therapies across hundreds of simulated growth conditions and pathogen metabolic states. Our approach can serve as a useful guide for the selection of robustly synergistic drug combinations.
N-glycan chitobiose core biosynthesis by Agl24 strengthens the hypothesis of an archaeal origin of the eukaryal N-glycosylation Protein N-glycosylation is the most common posttranslational modifications found in all three domains of life. The crenarchaeal N-glycosylation begins with the synthesis of a lipid-linked chitobiose core structure, identical to that in eukaryotes. Here, we report the identification of a thermostable archaeal beta-1,4-N-acetylglucosaminyltransferase, named archaeal glycosylation enzyme 24 (Agl24), responsible for the synthesis of the N-glycan chitobiose core. Biochemical characterization confirmed the function as an inverting {beta}-D-GlcNAc-(1[-&gt;]4)--D-GlcNAc-diphosphodolichol glycosyltransferase. Substitution of a conserved histidine residue, found also in the eukaryotic and bacterial homologs, demonstrated its functional importance for Agl24. Furthermore, bioinformatics and structural modeling revealed strong similarities between Agl24 and both the eukaryotic Alg14/13 and a distant relation to the bacterial MurG, which catalyze the identical or a similar process, respectively. Our data, complemented by phylogenetic analysis of Alg13 and Alg14, revealed similar sequences in Asgardarchaeota, further supporting the hypothesis that the Alg13/14 homologs in eukaryotes have been acquired during eukaryogenesis.

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=95 SRC="FIGDIR/small/427365v2_ufig1.gif" ALT="Figure 1">
View larger version (23K):
org.highwire.dtl.DTLVardef@13d1ffdorg.highwire.dtl.DTLVardef@1ff48d0org.highwire.dtl.DTLVardef@1056e75org.highwire.dtl.DTLVardef@196ff58_HPS_FORMAT_FIGEXP  M_FIG C_FIG HighlightsO_LIFirst identification and characterization of a thermostable {beta}-D-GlcNAc-(1[-&gt;]4)--D-GlcNAc-diphosphodolichol glycosyltransferase (GT family 28) in Archaea.
C_LIO_LIA highly conserved histidine, within a GGH motif in Agl24, Alg14, and MurG, is essential for function of Agl24.
C_LIO_LIAgl24-like homologs are broadly distributed among Archaea.
C_LIO_LIThe eukaryotic Alg13 and Alg14 are closely related to the Asgard homologs, suggesting their acquisition during eukaryogenesis.
C_LI Microbial Communities Performing Hydrogen Solventogenic Metabolism of Volatile Fatty Acids 1This work evaluates four different physicochemical pretreatments (acidic, thermal, acidic-thermal and thermal-acidic) on an anaerobic inoculum used for alcohol production from acetate and butyrate. All experiments were conducted in single batches using acetate and butyrate as substrates at 30 {degrees}C and with a pressurized headspace of pure H2 at 2.15 atm (218.2 MPa). Thermal and acidic-thermal pretreatments lead to higher production of both ethanol and butanol. Mathematical modelling shows that the highest attainable concentrations of ethanol and butanol produced were 122 mg L-1 and 97 mg L-1 for the thermal pretreatment (after 17.5 days) and 87 mg L-1 and 143 mg L-1 for the acidic-thermal pretreatment (after 18.9 days). Acetate was produced in all assays. Thermodynamic data indicated that a high H2 partial pressure favoured solventogenic metabolic pathways. Finally, sequencing data showed that both thermal and acidic-thermal pretreatments selected mainly the bacterial genera Pseudomonas, Brevundimonas and Clostridium. The acidic-thermal pretreatment selected a bacterial community more adapted to the conversion of acetate and butyrate into ethanol and butanol, respectively. Thermal-acidic pretreatment was unstable, showing significant variability between replicates. Acidic pretreatment showed the lowest alcohol production. CCL28 modulates neutrophil responses during infection with mucosal pathogens The chemokine CCL28 is highly expressed in mucosal tissues, but its role during infection is not well understood. Here we show that CCL28 promotes neutrophil accumulation to the gut of mice infected with Salmonella and to the lung of mice infected with Acinetobacter. Neutrophils isolated from the infected mucosa expressed the CCL28 receptors CCR3 and, to a lesser extent, CCR10, on their surface. The functional consequences of CCL28 deficiency were different between the two infections, as Ccl28-/- mice were highly susceptible to Salmonella gut infection, but highly resistant to otherwise lethal Acinetobacter lung infection. In vitro, unstimulated neutrophils harbored pre-formed intracellular CCR3 that was rapidly mobilized to the cell surface following phagocytosis or inflammatory stimuli. Moreover, CCL28 stimulation enhanced neutrophil antimicrobial activity, production of reactive oxygen species, and formation of extracellular traps, all processes that were largely dependent on CCR3. Consistent with the different outcomes in the two infection models, neutrophil stimulation with CCL28 boosted the killing of Salmonella but not of Acinetobacter. CCL28 thus plays a critical role in the immune response to mucosal pathogens by increasing neutrophil accumulation and activation, which can enhance pathogen clearance but also exacerbate disease depending on the mucosal site and the infectious agent. Gene loss during a transition to multicellularity Multicellular evolution is a major transition associated with momentous diversification of multiple lineages and increased developmental complexity. The volvocine algae comprise a valuable system for the study of this transition, as they span from unicellular to undifferentiated and differentiated multicellular morphologies despite their genomes being highly similar, suggesting multicellular evolution requires few genetic changes to undergo dramatic shifts in developmental complexity. Here, the evolutionary dynamics of six volvocine genomes were examined, where a gradual loss of genes was observed in parallel to the co-option of a few key genes. Protein complexes in the six species exhibited a high degree of novel interactions, suggesting that gene loss plays a role in evolutionary novelty. This finding was supported by gene network modeling, where gene loss outpaces gene gain in generating novel stable network states. These results suggest developmental complexity may be driven by gene loss, and not just by gene gain.

Significance StatementIncreased developmental complexity is thought to evolve mainly through genetic innovation and co-option. Comparison of the genomes of the closely related volvocines revealed that despite large changes in developmental complexity their genomes are exhibiting significant gene loss. Further, a burst of gene loss is shown to occur at the transition to undifferentiated multicellularity by gene inactivation and decay. This likely results in changes to protein-protein interactions within the cell, suggesting evolutionary novelty does not always require gene gain. Using empirical and modelling approaches we demonstrate gene loss can more easily produce protein network novelty than gene gain. We propose that gene loss can be a driver of biological innovation, allowing for reconfiguration of gene networks and differential use of existing functional repertoires. First evidence of in vitro cytotoxic effects of marine microlitter on Merluccius merluccius and Mullus barbatus, two Mediterranean commercial fish species Marine litter is composed mainly of plastics and is recognized as a serious threats to marine ecosystems. Ecotoxicological approaches have started elucidating the potential severity of microplastics (MPs) in controlled laboratory studies with pristine materials but no information exist on marine environmental microlitter as a whole. Here, we characterized the litter in the coastal Northern Tyrrhenian sea and in the stomach of two fish species of socio-economic importance, and exposed primary cell cultures of mucosal and lymphoid organs to marine microlitter for evaluating possible cytotoxic effects. An average of 0.30 {+/-} 0.02 microlitter items m-3 was found in water samples. FT-IR analysis revealed that plastic particles, namely HDPE, polyamide and polypropylene were present in 100% and 83.3% of Merluccius merluccius and Mullus barbatus analyzed, which overall ingested 14.67 {+/-} 4.10 and 5.50 {+/-} 1.97 items/individual, respectively. Moreover, microlitter was confirmed as a vector of microorganisms. Lastly, the apical end-point of viability was found to be significantly reduced in splenic cells exposed in vitro to two microlitter conditions. Considering the role of the spleen in the mounting of adaptive immune responses, our results warrant more in-depth investigations for clarifying the actual susceptibility of these two species to anthropogenic microlitter.

HighlightsO_LI0.30 {+/-} 0.02 microlitter items m-3 were found at the surface of coastal Northern Tyrrhenian sea
C_LIO_LI14.67 {+/-} 4.10 and 5.50 {+/-} 1.97 items/individual were retrieved from the stomach of hakes and mullets
C_LIO_LIThe ingested microlitter contained plastic items
C_LIO_LIMicrolitter was validated as a carrier of bacteria, fungi and flagellates
C_LIO_LISplenic cells exposed to two microlitter conditions for 72 hours suffered cytotoxicity
C_LI Climate warming changes synchrony of plants and pollinators Climate warming changes the phenology of many species. When interacting organisms respond differently, climate change may disrupt their interactions and affect the stability of ecosystems. Here, we used GBIF occurrence records to examine phenology trends in plants and their associated insect pollinators in Germany since the 1980s. We found strong phenological advances in plants, but differences in the extent of shifts among pollinator groups. The temporal trends in plant and insect phenologies were generally associated with interannual temperature variation, and thus likely driven by climate change. When examining the synchrony of species-level plant-pollinator interactions, their temporal trends differed among pollinator groups. Overall, plant-pollinator interactions become more synchronized, mainly since the phenology of plants responded more strongly to climate change than that of the pollinators. However, if the observed trends continue, many interactions may become more asynchronous again in the future. Our study suggests that climate change affects the phenologies of both plants and insects, and that it also influences the synchrony of plant-pollinator interactions. Reuse of a wing venation gene-regulatory network in patterning the eyespot rings of butterflies Novel organismal traits might reuse ancestral gene-regulatory networks (GRNs) in their development, but data supporting this mechanism are still sparse. Here we propose that a conserved venation gene network, used in patterning insect wings, has been recruited to pattern the sharp and distinct rings of color in butterfly eyespot development. This network involves Decapentaplegic (Dpp), acting as a central morphogen, activating optix and spalt at different concentration thresholds. It also involves Spalt cross-regulating optix resulting in the formation of a sharp boundary between the two genes. First, using in situ hybridizations, we show that dpp is expressed in the eyespot centers during the early pupal stage and that CRISPR-Cas9 embryonic knockouts of dpp, optix, and spalt result in the loss of eyespots, loss of the orange ring, and loss of the black scales, respectively. Furthermore, using RNAi during the pupal stage, we show that dpp down-regulation results in the complete or partial loss of eyespots. Using CRISPR disruptions of optix or spalt, followed by immunostainings, we show that Spalt represses optix in the central region of the eyespot, limiting optix expression to a more peripheral ring. We further illustrate that optix disruptions alter scale pigmentation and scale ultrastructure as they transform orange to brown scales. Altogether, the regulatory interactions between dpp, spalt and optix mimic those found in earlier larval developmental stages, in the anterior compartment of an insects wing, where these genes play a role in venation patterning. The network similarities suggest that the venation GRN was co-opted to aid in the differentiation of the eyespot rings. Text has been revised, image definition has been revised. Maize is a staple food all over the world. Models for climate change suggest that, in the future, cloud formation might be reduced in the tropics increasing the exposure to Ultraviolet-B (UV-B) radiation, a DNA-damaging agent. UV-B (290 to 320 nm) has been shown to affect yield in maize. In this project we have determined the differences in DNA repair efficiencies between U.S. inbred lines B73 and Mo17, and Central American purple landraces from Guatemala and Costa Rica. Our results from single cell electrophoresis experiments (Comet Assay) suggest that the landrace Pujagua Santa Cruz (P1, Costa Rican) was resistant to damage caused by the radiomimetic agent zeocin (24h/100 micrograms per mL), while landrace Pujagua La Cruz (P2, Costa Rican) was able to repair DNA damage after one hour. On the other hand, line Mo17 (Missouri, USA) was unable to repair the damage, while B73 (Iowa, USA) and the landraces Jocopilas (Guatemalan), Orotina Congo (Costa Rican) and Talamanca (Costa Rican) were partially able to repair DNA damage. High Resolution Melting (HRM) curve analysis of putative homologous DNA repair gene ZeaATM1 showed that both P1 and P2 had differences in the melting temperatures for this gene compared to B73 and Mo17, while P1 showed additional differences in ZeaSOG1, ZeaRAD51 and ZeaBRCA1, suggesting that in this landrace the presence of polymorphisms may be common among key genes for this pathway. Taken together our results suggest that key adaptive differences in DNA repair efficiencies exist between inbred lines and landraces of maize and that some Central America landraces could be used as a valuable pool of alleles for plant breeding aiming to increase tolerance to radiation. SUMM4 biochemically complex couples insulator function and DNA replication timingcontrol The mechanisms that establish DNA replication timing programs in eukaryotes remain incompletely understood. Drosophila SNF2-related factor SUUR imparts under-replication (UR) of late-replicating intercalary heterochromatin (IH) in polytene chromosomes. We developed a proteomics technique termed MERCI to isolate a native complex SUMM4 comprising SUUR and chromatin boundary protein Mod(Mdg4)-67.2. Mod(Mdg4) stimulates the ATPase activity of SUUR and is required for its normal spatiotemporal distribution in vivo. Both SUMM4 subunits mediate the activities of gypsy insulator disrupting enhancer-promoter interactions and establishing chromatin barriers. Furthermore, SuUR or mod(mdg4) mutations reverse UR of IH. Our findings uncover a critical role for architectural proteins in attenuating replication fork progression and suggest an alternative mechanism for DNA replication timing that does not depend on an asynchronous firing of replication origins.

One-Sentence SummaryA stable protein complex comprising an insulator factor and a SNF2-like ATPase imparts late replication of heterochromatin. 3D Reconstructions of Mouse Skeletal Muscle Reveal a Decrease in the MICOS Complex and Altered Mitochondrial Networks Background: Skeletal muscle gradually loses mass, strength, endurance, and oxidative capacity during aging. Studies of bioenergetics and protein turnover show that mitochondria mediate this decline in function. Mitochondria cristae are essential for the production of ATP. While mitochondrial aging is associated with endoplasmic reticulum stress, fragmented mitochondria, and decreased mitochondrial capacity, the genes associated with morphological changes in mitochondria during aging still requires further elucidation Further, it is not completely understood how 3D mitochondrial networks and the specialization of mitochondria are altered during aging. Methods: We measured changes in mitochondrial morphology and mitochondrial connectivity during the aging of the mouse gastrocnemius muscle through serial block facing-scanning electron microscopy and 3D reconstruction. Nanotunnels were also measured through 3D reconstruction. CRISPR/Cas9 KO was performed to examine changes in mitochondria upon loss of the mitochondrial contact site and cristae organizing system (MICOS) complex and optic atrophy 1 (OPA-1). Metabolomics were used to find key metabolite pathways changed upon off of the MICOS complex. Results: We found changes in mitochondrial network configuration, nanotunneling, size, shape, number, contact sites, cristae organizing system (MICOS) dynamics and gene expression in skeletal muscle during aging. We also found an association of OPA-1 and the MICOS complex in the gastrocnemius with mitochondrial aging. Further, the loss of the MICOS complex was linked with decreased oxidative capacity, and altered mitochondrial metabolism. Conclusions: To our knowledge, studies of mitochondrial changes in skeletal muscle throughout aging remains limited. 3D reconstructions of nanotunnels elucidated novel patterns of mitochondrial aging in skeletal muscle. MICOS proteins decreased with age and mitochondrial morphology was similar between aged skeletal muscle and that of young mice with MICOS protein loss. In tandem, our data suggest a relationship between the MICOS complex and aging, which could be potentially linked to disease states with additional 3D reconstruction. Kre28-Spc105 interaction is essential for Spc105 loading at the kinetochore. Kinetochores are macromolecular protein assemblies that attach sister chromatids to spindle microtubules and mediate accurate chromosome segregation during mitosis. The outer kinetochore consists of the KMN network, a protein super-complex comprising Knl1 (yeast Spc105), Mis12 (yeast Mtw1), and Ndc80 (yeast Ndc80), which harbors sites for microtubule binding. Within the KMN network, Spc105 acts as an interaction hub of components involved in spindle assembly checkpoint (SAC) signaling. It is known that Spc105 forms a complex with kinetochore component Kre28. However, where Kre28 physically localizes in the budding yeast kinetochore is not clear. The exact function of Kre28 at the kinetochore is also unknown. Here, we investigate how Spc105 and Kre28 interact and how they are organized within bioriented yeast kinetochores using genetics and cell biological experiments. Our microscopy data show that Spc105 and Kre28 localize at the kinetochore with a 1:1 stoichiometry. We also show that the Kre28-Spc105 interaction is important for Spc105 protein turn-over and essential for their mutual recruitment at the kinetochores. We created several truncation mutants of kre28 that affect Spc105 loading at the kinetochores. When over-expressed, these mutants sustain the cell viability, but SAC signaling and kinetochore biorientation are impaired. Therefore, we conclude that Kre28 contributes to chromosome biorientation and high-fidelity segregation at least indirectly by regulating Spc105 localization at the kinetochores. The relative speed of the glucorticoid stress response varies independently of scope and is predicted by environmental variability and longevity across birds. The acute glucocorticoid response is a key mediator of the coordinated vertebrate response to unpredictable challenges. Rapid glucocorticoid increases initiate changes that allow animals to cope with stressors. The scope of the glucocorticoid response -- defined here as the absolute increase in glucocorticoids -- is associated with individual differences in performance and varies across species with environment and life history. In addition to varying in scope, responses can differ enormously in speed; however, relatively little is known about whether speed and absolute glucocorticoid levels covary, how selection shapes speed, or what aspects of speed are important. We used corticosterone samples collected at 5 time points from 1,750 individuals of 60 species of birds to ask i) how the speed and scope of the glucocorticoid response covary and ii) whether variation in absolute or relative speed is predicted by environmental context or life history. As predicted by a recent model, faster absolute glucocorticoid responses were strongly associated with a larger scope. Despite this covariation, the relative speed of the glucocorticoid response (standardized within species) varied independently of absolute scope, suggesting that selection could operate on both features independently. Species with faster relative glucocorticoid responses lived in locations with more variable temperature and had shorter lifespans. Our results suggest that rapid changes associated with the speed of the glucocorticoid response, such as those occurring through non-genomic receptors, might be an important determinant of coping ability and we emphasize the need for studies designed to measure speed independently of absolute glucocorticoid levels. Polycystin promotes renal fibrosis through EZH2 Autosomal polycystic kidney disease (ADPKD) is caused by mutations in PKD1or PKD2 gene, which encodes polycystin-1 or polycystin-2 respectively. Paradoxically, polycystins are up-regulated in human ADPKD kidneys and overexpression of polycystin-1 even induced renal cyst formation in mice, suggesting that up-regulation of polycystin is pathogenic. In this study, we determined the role of polycystin in renal fibrosis and explored its underlying mechanisms. We found that the expression of polycystin-1 or polycystin-2 was positively correlated with the expression of collagen-I in the kidney of three different fibrotic mouse models. Pharmaceutical inhibition or genetic knockout of polycystin-2 reduced the expression of pro-fibrotic markers in fibrotic kidneys of two animal models. Moreover, conditional knockout of Pkd1 gene attenuated renal fibrosis in two mouse models. Enhancer of Zeste homolog 2 (EZH2) is a methyltransferase which plays important role in renal fibrosis. We found that EZH2 was positively correlated with the expression of polycystins in fibrotic kidneys. Conditional knockout of the Ezh2 attenuated the anti-fibrotic responses induced by Pkd1 deletion in fibrotic kidneys.

In conclusion, polycystins are up-regulated in injured kidneys and polycystins promote renal fibrosis through up-regulation of EZH2 in fibrotic kidneys, suggesting that renal fibrosis in polycystic kidneys is not the direct effect of PKD gene mutations.

SignificanceMutations of polycystins cause autosomal polycystic kidney disease (ADPKD) which is associated with increased renal fibrosis in cystic kidneys. However, the role of polycystin in renal fibrosis remains unclear. Here, we showed that the expression of polycystins is positively correlated with the development of renal fibrosis in injured kidneys. Pharmaceutical inhibition or genetic deletion of polycystin-2 attenuated renal fibrosis in fibrotic kidneys. Moreover, genetic deletion of polycystin-1 inhibited renal fibrosis in two different animal models. Finally, we showed that ablation of the methyltransferase EZH2 mitigated the anti-fibrotic effect of polycystin-1 deletion in fibrotic kidneys. We conclude that polycystin promotes through EZH2 and renal fibrosis in cystic kidney is not the direct effect of PKD gene mutations. Diverse stem-chondrichthyan oral structures and evidence for an independently acquired acanthodid dentition The teeth of sharks famously form a series of transversely-organised files with conveyor-belt replacement that are borne directly on the jaw cartilages, in contrast to the dermal plate-borne dentition of bony fishes that undergoes site-specific replacement. A major obstacle in understanding how this system evolved is the poorly understood relationships of the earliest chondrichthyans and the profusion of morphologically and terminologically diverse bones, cartilages, splints and whorls that they possess. Here we use tomographic methods to investigate mandibular structures in several early branching  acanthodian-grade stem-chondrichthyans. We show that the dentigerous jaw bones of disparate genera of ischnacanthids are united by a common construction, being growing bones with non-shedding dentition. Mandibular splints, which support the ventro-lateral edge of the Meckels cartilage in some taxa, are formed from dermal bone and may be an acanthodid synapomorphy. We demonstrate that the teeth of Acanthodopsis are borne directly on the mandibular cartilage and that this taxon is deeply nested within an edentulous radiation, representing an unexpected independent origin of teeth. Many or even all of the range of unusual oral structures may be apomorphic, but they should nonetheless be considered when building hypotheses of tooth and jaw evolution, both in chondrichthyans and more broadly. From fitting the average to fitting the individual: A cautionary tale for mathematicalmodelers An outstanding challenge in the clinical care of cancer is moving from a one-size-fits-all approach that relies on population-level statistics towards personalized therapeutic design. Mathematical modeling is a powerful tool in treatment personalization, as it allows for the incorporation of patient-specific data so that treatment can be tailor-designed to the individual. Herein, we work with a mathematical model of murine cancer immunotherapy that has been previously-validated against the average of an experimental dataset. We ask the question: what happens if we try to use this same model to perform personalized fits, and therefore make individualized treatment recommendations? Typically, this would be done by choosing a single fitting methodology, and a single cost function, identifying the individualized best-fit parameters, and extrapolating from there to make personalized treatment recommendations. Our analyses show the potentially problematic nature of this approach, as predicted personalized treatment response proved to be sensitive to the fitting methodology utilized. We also demonstrate how a small amount of the right additional experimental measurements could go a long way to improve consistency in personalized fits. Finally, we show how quantifying the robustness of the average response could also help improve confidence in personalized treatment recommendations.

Author summaryAs we enter the era of healthcare where personalized medicine becomes a more common approach to treating cancer patients, harnessing the power of mathematical models will only become more essential. Using a preclinical dataset on cancer immunotherapy, we explore the challenges and limitations that arise when trying to move from fitting and making predictions for the population-level average, to fitting and making predictions for an individual. We find that the standard of approach of picking a single fitting methodology and a single cost function may end up having limited predictive value when applied to individual data. We also show how having a small amount of the right additional experimental data, and establishing the robustness of average treatment response, can help improve confidence in personalized model predictions. Ultrahigh-Throughput Ambient MS: Direct Analysis at 22 Samples per Second by Infrared Matrix-Assisted Laser Desorption Electrospray Ionization Mass Spectrometry Infrared Matrix-Assisted Laser Desorption Electrospray Ionization (IR-MALDESI) mass spectrometry is an ambient-direct sampling method being developed for high-throughput, label-free, biochemical screening of large-scale compound libraries. Here, we report the development of an ultrahigh-throughput continuous motion IR-MALDESI sampling approach capable of acquiring data at rates up to 22.7 samples per second in a 384-well microtiter plate. At top speed, less than 1% analyte carryover is observed from well-to-well and signal intensity relative standard deviations (RSD) of 11.5% and 20.9% for 3 M 1-hydroxymidazolam and 12 M dextrorphan, respectively, are achieved. The ability to perform parallel kinetics studies on 384 samples with ~30s time resolution using an isocitrate dehydrogenase 1 (IDH1) enzyme assay is shown. Finally, we demonstrate the repeatability and throughput of our approach by measuring 115,200 samples from 300 microtiter plate reads consecutively over 5.54 hours with RSDs under 8.14% for each freshly introduced plate. Taken together, these results demonstrate the use of IR-MALDESI at sample acquisition rates that surpass other currently reported direct sampling mass spectrometry approaches used for high throughput compound screening.

For Table of Contents Only

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=125 SRC="FIGDIR/small/465730v3_ufig1.gif" ALT="Figure 1">
View larger version (49K):
org.highwire.dtl.DTLVardef@5f916corg.highwire.dtl.DTLVardef@4473b4org.highwire.dtl.DTLVardef@82e8forg.highwire.dtl.DTLVardef@128b57b_HPS_FORMAT_FIGEXP  M_FIG C_FIG Embracing curiosity eliminates the exploration-exploitation dilemma Balancing exploration with exploitation is seen as a mathematically intractable dilemma that all animals face. In this paper, we provide an alternative view of this classic problem that does not depend on exploring to optimize for reward. We argue that the goal of exploration should be pure curiosity, or learning for learnings sake. Through theory and simulations we prove that explore-exploit problems based on this can be solved by a simple rule that yields optimal solutions: when information is more valuable than rewards, be curious, otherwise seek rewards. We show that this rule performs well and robustly under naturalistic constraints. We suggest three criteria can be used to distinguish our approach from other theories. Optogenetic stimulation of lateral hypothalamic orexin/dynorphin inputs in the ventral tegmental area potentiates mesolimbic dopamine neurotransmission and promotes reward-seeking behaviours. Reward and reinforcement processes are critical for survival and propagation of genes. While numerous brain systems underlie these processes, a cardinal role is ascribed to mesolimbic dopamine. However, ventral tegmental area (VTA) dopamine neurons receive complex innervation and various neuromodulatory factors, including input from lateral hypothalamic (LH) orexin/hypocretin neurons which also express and co-release the neuropeptide, dynorphin. Dynorphin in the VTA induces aversive conditioning through the Kappa opioid receptor (KOR) and decreases dopamine when administered intra-VTA. Exogenous application of orexin or orexin 1 receptor (oxR1) antagonists in the VTA bidirectionally modulates dopamine-driven motivation and reward-seeking behaviours, including the attribution of motivational value to primary rewards and associated conditioned stimuli. However, the effect of endogenous stimulation of LH orexin/dynorphin-containing projections to the VTA and the potential contribution of co-released dynorphin on mesolimbic dopamine and reward related processes remains uncharacterised. We combined optogenetic, electrochemical, and behavioural approaches to examine this. We found that optical stimulation of LH orexin/dynorphin inputs in the VTA potentiates mesolimbic dopamine neurotransmission in the nucleus accumbens (NAc) core, produces real time and conditioned place preference, and increases the food cue-directed orientation in a Pavlovian conditioning procedure. LH orexin/dynorphin potentiation of NAc dopamine release and real time place preference was blocked by an oxR1, but not KOR antagonist. Thus, rewarding effects associated with optical stimulation of LH orexin/dynorphin inputs in the VTA are predominantly driven by orexin rather than dynorphin. Trade-offs in biosensor optimization for dynamic pathway engineering Recent progress in synthetic biology allows the construction of dynamic control circuits for metabolic engineering. This technology promises to overcome many challenges encountered in traditional pathway engineering, thanks to their ability to self-regulate gene expression in response to bioreactor perturbations. The central components in these control circuits are metabolite biosensors that read out pathway signals and actuate enzyme expression. However, the construction of metabolite biosensors is a major bottleneck for strain design, and a key challenge is to understand the relation between biosensor dose-response curves and pathway performance. Here we employ multiobjective optimization to quantify performance trade-offs that arise in the design and calibration of metabolite biosensors. Our approach reveals strategies for tuning dose-response curves along an optimal trade-off between production flux and the cost of an increased expression burden on the host. We explore properties of control architectures built in the literature, and identify their advantages and caveats in terms of performance and robustness to growth conditions and leaky promoters. We demonstrate the optimality of a control circuit for glucaric acid production in Escherichia coli, which has been shown to increase titer by 2.5-fold as compared to static designs. Our results lay the groundwork for the automated design of control circuits for pathway engineering, with applications in the food, energy and pharmaceutical sectors. Prediction of the sex-determination gene in tunas (Thunnus fishes) Fish species have a variety of sex determination systems. Tunas (genus Thunnus) have an XY genetic sex-determination system. However, the Y chromosome or responsible locus has not yet been identified in males. In a previous study, a female genome of Pacific bluefin tuna (T. orientalis) was sequenced, and candidates for sex-associated DNA polymorphisms were identified by a genome-wide association study using resequencing data. In the present study, we sequenced a male genome of Pacific bluefin tuna by long-read and linked-read sequencing technologies, and explored male-specific loci through a comparison with the female genome. As a result, we found a unique region carrying the male-specific haplotype, where a homolog of estrogen sulfotransferase gene was predicted to be encoded. The genome-wide mapping of previously resequenced data indicated that, among the functionally annotated genes, only this gene, named sult1st6y, was paternally inherited in the males of Pacific bluefin tuna. We reviewed the RNA-seq data of southern bluefin tuna (T. maccoyii) in the public database and found that sult1st6y of southern bluefin tuna was expressed in all male testes, but absent or suppressed in the female ovary. Since estrogen sulfotransferase is responsible for the inactivation of estrogens, it is reasonable to assume that the expression of sult1st6y in gonad cells may inhibit female development, thereby inducing the individuals to become males. Thus, our results raise a promising hypothesis that sult1st6y is the sex-determination gene in Thunnus fishes, or at least functions at a crucial point in the sex-differentiation cascade. Computational study and design of effective siRNAs to silence structural proteins associated genes of Indian SARS-CoV-2 strains SARS-CoV-2 is a highly transmissible and pathogenic coronavirus that first emerged in late 2019 and has since triggered a pandemic of acute respiratory disease named COVID-19 which poses a significant threat to all public health institutions in the absence of specific antiviral treatment. Since the outbreak began in March 2020, India has reported 4.77 lakh Coronavirus deaths, according to the World Health Organization (WHO). The innate RNA interference (RNAi) pathway, on the other hand, allows for the development of nucleic acid-based antiviral drugs in which complementary small interfering RNAs (siRNAs) mediate the post-transcriptional gene silencing (PTGS) of target mRNA. Therefore, in this current study, the potential of RNAi was harnessed to construct siRNA molecules that target the consensus regions of specific structural proteins associated genes of SARS-CoV-2, such as the envelope protein gene (E), membrane protein gene (M), nucleocapsid phosphoprotein gene (N), and surface glycoprotein gene (S) which are important for the viral pathogenesis. Conserved sequences of 811 SARS-CoV-2 strains from around India were collected to design 21 nucleotides long siRNA duplex based on various computational algorithms and parameters targeting E, M, N and S genes. The proposed siRNA molecules possessed sufficient nucleotide-based and other features for effective gene silencing and BLAST results revealed that siRNAs targets have no significant matches across the whole human genome and hence, siRNAs were found to have no off-target effects on the genome, ruling out the possibility of off-target silencing. Finally, out of 157 computationally identified siRNAs, only 4 effective siRNA molecules were selected for each target gene which is proposed to exert the best action based on GC content, free energy of folding, free energy of binding, melting temperature, heat capacity and molecular docking analysis with Human AGO2 protein. Our engineered siRNA candidates could be used as a genome-level therapeutic treatment against various sequenced SARS-CoV-2 strains in India. However, future applications will necessitate additional validations in vitro and in vivo animal models. Changes to the mtDNA copy number during yeast culture growth We show that the mtDNA copy number in growing cultures of the yeast Saccharomyces cerevisiae increases by a factor of up to 4, being lowest ([~]10 per haploid genome) and stable during rapid fermentative growth, and highest at the end of the respiratory phase. When yeast are grown on glucose, the onset of the mtDNA copy number increase coincides with the early stages of the diauxic shift, and the increase continues through respiration. A lesser yet still substantial copy number increase occurs when yeast are grown on a nonfermentable carbon source, i.e. when there is no diauxic shift. The mtDNA copy number increase during and for some time after the diauxic shift is not driven by an increase in cell size. The copy number increase occurs in both haploid and diploid strains, but is markedly attenuated in a diploid wild isolate that is a ready sporulator. Strain-to-strain differences in mtDNA copy number are least apparent in fermentation and most apparent in late respiration or stationary phase. While changes in mitochondrial morphology and function were previously known to accompany changes in physiological state, it had not been previously shown that the mtDNA copy number changes substantially over time in a clonal growing culture. The mtDNA copy number in yeast is therefore a highly dynamic phenotype. IRGM/Irgm1 Aggravates Progression of Atherosclerosis by Inducing Macrophage Apoptosis through the MAPK Signaling Pathway AimsAtherosclerosis underlies most cardiovascular diseases, among which acute coronary syndrome (ACS) caused by plaque rupture (PR) often leads to death. Immune-related GTPases (IRGM/Irgm1) have been extensively studied in inflammatory diseases, but their role in atherosclerosis is unclear. Determining how IRGM/Irgm1 promotes atherosclerotic plaque vulnerability will provide information for new biomarkers and/or therapeutic targets.

Methods and resultsWe identified ruptured and unruptured plaques by optical coherence tomography, and found that serum IRGM was highly expressed in patients with ST-segment elevation myocardial infarction in PR. We used ApoE-/-Irgm1+/+, ApoE-/-Irgm1+/- mice and chimeric mice to establish a model of advanced atherosclerosis. The results of pathological experiments showed that Irgm1 caused plaque necrosis. The ratio of neutral lipids and cholesterol crystals increases, while the content of collagen fibers decreases, aggravating the destabilization of atherosclerotic plaques. In vitro, we used multiple approaches to confirm that Irgm1 promotes macrophage apoptosis by promoting the production of reactive oxygen species and activating the MAPK signaling pathway.

ConclusionsIRGM may be a potential risk factor for PR. Mechanistic studies have shown that IRGM/Irgm1 contributes to the formation and rupture of fragile plaques. This is partly mediated by the induction of macrophage apoptosis via the MAPK signaling pathway. IRGM may offer new strategies for early treatment of ACS.

Translation viewOur findings indicate that IRGM/Irgm1 contributes to formation and rupture of vulnerable plaques. It suggests that IRGM may provide a new target for the early treatment of ACS. A tool for assessing alignment of biomedical data repositories with open, FAIR, citation and trustworthy principles Increasing attention is being paid to the operation of biomedical data repositories in light of efforts to improve how scientific data is handled and made available for the long term. Multiple groups have produced recommendations for functions that biomedical repositories should support, with many using requirements of the FAIR data principles as guidelines. However, FAIR is but one set of principles that has arisen out of the open science community. They are joined by principles governing open science, data citation and trustworthiness, all of which are important aspects for biomedical data repositories to support. Together, these define a framework for data repositories that we call OFCT: Open, FAIR, Citable and Trustworthy. Here we developed an instrument using the open source PolicyModels toolkit that attempts to operationalize key aspects of OFCT principles and piloted the instrument by evaluating eight biomedical community repositories listed by the NIDDK Information Network (dkNET.org). Repositories included both specialist repositories that focused on a particular data type or domain, in this case diabetes and metabolomics, and generalist repositories that accept all data types and domains. The goal of this work was both to obtain a sense of how much the design of current biomedical data repositories align with these principles and to augment the dkNET listing with additional information that may be important to investigators trying to choose a repository, e.g., does the repository fully support data citation? The evaluation was performed from March to November 2020 through inspection of documentation and interaction with the sites by the authors. Overall, although there was little explicit acknowledgement of any of the OFCT principles in our sample, the majority of repositories provided at least some support for their tenets.