The genome of the cereal pest Sitophilus oryzae: a transposable element haven BackgroundThe rice weevil Sitophilus oryzae is one of the most important agricultural pests, causing extensive damage to cereal in fields and to stored grains. S. oryzae has an intracellular symbiotic relationship (endosymbiosis) with the Gram-negative bacterium Sodalis pierantonius and is a valuable model to decipher host-symbiont molecular interactions.

ResultsWe sequenced the Sitophilus oryzae genome using a combination of short and long reads to produce the best assembly for a Curculionidae species to date. We show that S. oryzae has undergone successive bursts of transposable element (TE) amplification, representing 72% of the genome. In addition, we show that many TE families are transcriptionally active, and changes in their expression are associated with insect endosymbiotic state. S. oryzae has undergone a high gene expansion rate, when compared to other beetles. Reconstruction of host-symbiont metabolic networks revealed that, despite its recent association with cereal weevils (30 Kyear), S. pierantonius relies on the host for several amino acids and nucleotides to survive and to produce vitamins and essential amino-acids required for insect development and cuticle biosynthesis.

ConclusionsHere we present the genome of an agricultural pest beetle, which may act as a foundation for pest control. In addition, S. oryzae may be a useful model for endosymbiosis, and studying TE evolution and regulation, along with the impact of TEs on eukaryotic genomes.
Epigenetic scores for the circulating proteome as tools for disease prediction Protein biomarkers have been identified across many age-related morbidities. However, characterising epigenetic influences could further inform disease predictions. Here, we leverage epigenome-wide data to study links between the DNAm signatures of the circulating proteome and incident diseases. Using data from four cohorts, we trained and tested epigenetic scores (EpiScores) for 953 plasma proteins, identifying 109 scores that explained between 1% and 58% of the variance in protein levels after adjusting for known protein quantitative trait loci (pQTL) genetic effects. By projecting these EpiScores into an independent sample, (Generation Scotland; n=9,537) and relating them to incident morbidities over a follow-up of 14 years, we uncovered 137 EpiScore - disease associations. These associations were largely independent of immune cell proportions, common lifestyle and health factors and biological aging. Notably, we found that our diabetes-associated EpiScores highlighted previous top biomarker associations from proteome-wide assessments of diabetes. These EpiScores for protein levels can therefore be a valuable resource for disease prediction and risk stratification.
NUDT18 catalyzes the hydrolysis of active metabolites of the antivirals Remdesivir, Ribavirin and Molnupiravir Remdesivir and Molnupiravir have gained considerable interest due to their activity against SARS-CoV-2. Cellular hydrolysis of their active triphosphate forms, Remdesivir-TP and Molnupiravir-TP, would decrease drug efficiency. We therefore tested Remdesivir-TP as a substrate against a panel of human hydrolases and found that NUDT18 catalyzes the hydrolysis of Remdesivir-TP. The kcat value of NUDT18 for Remdesivir-TP was determined to 2.6 s-1 and the Km value was 156 M, suggesting that NUDT18 catalyzed hydrolysis of Remdesivir-TP occurs in cells. We demonstrate that the triphosphates of the antivirals Ribavirin and Molnupiravir are hydrolyzed by NUDT18, albeit with a lower efficiency compared to Remdesivir-TP. NUDT18 also hydrolyses the triphosphates of Sofosbuvir and Aciclovir although with significantly lower activity. These results suggest that NUDT18 can act as a cellular sanitizer of modified nucleotides and may influence the antiviral efficacy of Remdesivir, Molnupiravir and Ribavirin. NUDT18 is expressed in respiratory epithelial cells and may limit the antiviral efficacy of Remdesivir and Molnupiravir against SARS-CoV2 replication by decreasing the intracellular concentration of their active metabolites at their intended site of action. Alkaline phosphatase activity of serum affects osteogenic differentiation cultures Fetal bovine serum (FBS) is a widely used supplement in cell culture medium, despite its known variability in composition which greatly affects cellular function and consequently the outcome of studies. In bone tissue engineering, the deposited mineralized matrix is one of the main outcome parameters, but using different brands of FBS can result in large variations. Alkaline phosphatase (ALP) is present in FBS. Not only is ALP used to judge the osteogenic differentiation of bone cells, it may affect deposition of mineralized matrix. The present study focused on the enzymatic activity of ALP in FBS of different suppliers and its contribution to mineralization in osteogenic differentiation cultures. It was hypothesized that culturing cells in a medium with high intrinsic ALP activity of FBS will lead to higher mineral deposition compared to media with lower ALP activity. The used FBS types were shown to have significant differences in enzymatic ALP activity. Our results indicate that the ALP activity of the medium not only affected the deposited mineralized matrix but also the osteogenic differentiation of cells as measured by a changed cellular ALP activity of human bone marrow derived mesenchymal stromal cells (hBMSC). In media with low inherent ALP activity, the cellular ALP activity was increased and played the major role in the mineralization process; while, in media with high intrinsic ALP activity contribution from the serum, less cellular ALP activity was measured and the ALP activity of the medium also contributed to mineral formation substantially. Our results highlight the diverse effects of ALP activity intrinsic to FBS on osteogenic differentiation and matrix mineralization and how FBS can determine the experimental outcomes, in particular for studies investigating matrix mineralization. Once again, the need to replace FBS with more controlled and known additives is highlighted.

Graphical abstract

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=55 SRC="FIGDIR/small/462693v2_ufig1.gif" ALT="Figure 1">
View larger version (17K):
org.highwire.dtl.DTLVardef@902078org.highwire.dtl.DTLVardef@43b6fborg.highwire.dtl.DTLVardef@1a244fborg.highwire.dtl.DTLVardef@5ef819_HPS_FORMAT_FIGEXP  M_FIG C_FIG Tunable Cellular Localization and Extensive Cytoskeleton-Interplay of Reflectins Reflectin proteins are nature copolymers consist of repeated canonical domains, locating in a biophotonic system called Bragg lamellae and manipulating the dynamic structural coloration of iridocytes. Their biological functions are intriguing, but the underlying mechanism is not fully understood. In this study, reflectin A1, A2, B1 and C were found to present distinguished cyto-/nucleoplasmic localization preferences. Comparable intracellular localization can be reproduced by truncated reflectin variants, suggesting a conceivable evolutionary order among reflectin proteins. Secondly, the size-dependent access of reflectin variants into nucleus demonstrate a potential model of how reflectins get into Bragg lamellae. Thirdly, RfA1 was found to extensively interact with cytoskeleton, including its binding to actin and microtubule minus-end-directed movement. This implies cytoskeleton system plays fundamental role during the organization and transportation of reflectin proteins. Findings presented here recommend reflectins as programmable biomaterials which can be used to decipher their evolution processes, to delineate their biological mechanism, and to achieve tunable intracellular targeting as editable tags. How do stomata evolve on the road to C4 photosynthesis? C4 photosynthesis optimizes plants carbon and water relations, allowing high photosynthetic rate with low stomatal conductance. Stomata have long been believed as a part of C4 syndrome. However, it remains unclear how stomata traits evolved along the path from C3 to C4. Stomatal patterning was examined in Flaveria genus, a model for studying C4 evolution. Comparative, transgenic and semi-in-vitro experiments were used to study molecular basis that underlies stomatal traits along C4 evolution. Novel results: the evolution from C3 to C4 species through intermediate species is accompanied by a stepwise rather than an abrupt change in the stomatal traits. The initial change occurs near Type II and dramatic change occurs at the C4-like species. On the road to C4, stomata become less in number but bigger in size and changes in stomatal density dominates changes in maximum stomatal conductance (gsmax). The reduction of FSTOMAGEN expression underlies altered gsmax between Flaveria species with different photosynthetic pathways and likely occurs in other C4 lineages. Our study provides insight into the pattern, mechanism and role of stomatal evolution along the road towards C4. This work highlights the stomatal traits in the current C4 evolutionary model and the co-evolution of photosynthetic pathway and stomata. RNA-guided Cas12a and Cas9 knockout and knockin of the omega-1 locus of Schistosoma mansoni We compared the efficiency of gene knockout (KO) and precision of insertion (knock-in, KI) of the RNA-guided AsCas12a nuclease of Acidaminococcus sp. with that of SpCas9 from Streptococcus pyogenes, aiming to enhance the functional genomics toolkit for Schistosoma mansoni. Programmed DNA cleavages catalyzed by Cas12a and Cas9 result in staggered and blunt ended strand breaks, respectively. TTTV, the optimal protospacer adjacent motif for AsCas12a would occur frequently within the AT-rich genome of this platyhelminth. We deployed optimized conditions for the ratio of guide RNAs to the nuclease, donor templates, and electroporation parameters, to target a key enzyme termed omega-1 that is secreted by the schistosome egg. AsCas12a was more efficient than SpCas9 for gene knockout of omega-1 as determined by tracking of indels by decomposition (P < 0.001). Resulting from CRISPREsso2 analysis, most mutations were deletions; SpCas9 induced short deletions of 3 nt in length whereas AsCas12a induced deletions of 2 to 26 nt. Knockout efficiency of both nucleases markedly increased in the presence of short, single stranded oligodeoxynucleotide (ssODN) donor templates. With AsCas12a, ssODNs representative of both the non-CRISPR target (NT) and target (T) strands of the targeted gene were tested, resulting in KO efficiencies of 15.67, 28.71 and 21.43% in the SpCas9 plus donor ssODN, AsCas12a plus NT-ssODN, and AsCas12a plus T-ssODN groups, respectively. Trans cleavage activity against the ssODNs by activated AsCas12a was not apparent in vitro. Programmed SpCas9 editing led to more precise transgene insertion than AsCas12a, with KI efficiencies of 17.07% for the KI_SpCas9 group, 14.58% for KI_AsCas12a-NT-ssODN and 12.37% for KI_AsCas12a-T-ssODN. Although AsCas12a induced fewer mutations per genome than SpCas9, the phenotypic impact on transcription and expression of omega-1 was similar for both nucleases. These findings revealed that AsCas12a and SpCas9 both provide tractable routes for RNA-guided programmed mutation of the genome of the schistosome egg. Are microbes growing on flowers evil? Effects of old flower microbes on fruit set in a wild ginger with one-day flowers, Alpinia japonica (Zingiberaceae) Flowers are colonized and inhabited by diverse microbes. Plants rapidly replace flowers with short lifespan, and old flowers senesce. This may contribute to avoiding adverse effects of the microbes. In this study, we investigate if the flower microbial community on old flowers impedes fruit and seed production in a wild ginger with one-day flowers. We inoculated newly opened flowers with old flower microbes, and monitored the effects on fruit and seed set. We also assessed prokaryotic communities on the flowers using amplicon sequencing. We found six bacterial amplicon sequence variants (ASVs) whose proportions were increased on the inoculated flowers. These ASVs were also found on flower buds and flowers that were bagged by net or paper during anthesis, suggesting that they had been present in small numbers prior to flowering. Fruit set was negatively associated with the proportions of these ASVs, while seed set was not. The results suggest that old flowers harbor microbial communities different from those at anthesis, and that the microbes abundant on old flowers negatively affect plant reproduction. Though the short lifespan of flowers has gotten little attention, it might be an essential defense mechanism to cope with antagonistic microbes that rapidly proliferate on the flowers. Using machine learning to design adeno-associated virus capsids with high likelihood of viral assembly We study the application of machine learning in designing adeno-associated virus (AAV) capsid sequences with high likelihood of viral assembly, i.e. capsid viability. Specifically, we design and implement Origami, a model-based optimization algorithm, to identify highly viable capsid sequences within the vast space of 2033 possibilities. Our evaluation shows that Origami performs well in terms of optimality and diversity of model-designed sequences. Moreover, these sequences are ranked according to their viability score. This helps designing experiments given budget constraint. Faster-haplodiploid evolution under divergence-with-gene-flow: simulations and empirical data from pine-feeding hymenopterans Although haplodiploidy is widespread in nature, the evolutionary consequences of this sex determination mechanism are not well characterized. Here, we examine how genome-wide hemizygosity and a lack of recombination in haploid males affects genomic differentiation in populations that diverge via natural selection while experiencing gene flow. First, we simulated diploid and haplodiploid "genomes" (500-kb loci) evolving under an isolation-with-migration model with mutation, drift, selection, migration, and recombination; and examined differentiation at neutral sites both tightly and loosely linked to a divergently selected site. So long as there is divergent selection and migration, sex-limited hemizygosity and recombination cause elevated differentiation (i.e., produce a "faster-haplodiploid effect") in haplodiploid populations relative to otherwise equivalent diploid populations, for both recessive and codominant mutations. Second, we used genome-wide SNP data to model divergence history and describe patterns of genomic differentiation between sympatric populations of Neodiprion lecontei and N. pinetum, a pair of pine sawfly species (order: Hymenoptera; family: Diprionidae) that are specialized on different pine hosts. These analyses support a history of continuous gene exchange throughout divergence and reveal a pattern of heterogeneous genomic differentiation that is consistent with divergent selection on many unlinked loci. Third, using simulations of haplodiploid and diploid populations evolving according to the estimated divergence history of N. lecontei and N. pinetum, we found that divergent selection would lead to higher differentiation in haplodiploids. Based on these results, we hypothesize that haplodiploids undergo divergence-with-gene-flow and sympatric speciation more readily than diploids. Quantifying ERK-activity in response to inhibition of the BRAFV600E-MEK-ERK cascade using mathematical modelling BackgroundSimultaneous inhibition of multiple components of the BRAF-MEK-ERK cascade (vertical inhibition) has become a standard of care for treating BRAF-mutant melanoma. However, the molecular mechanism of how vertical inhibition synergistically suppresses intracellular ERK-activity, and consequently cell proliferation, are yet to be fully elucidated.

MethodsWe develop a mechanistic mathematical model that describes how the mutant BRAF-inhibitor, dabrafenib, and the MEK-inhibitor, trametinib, affect BRAFV600E-MEK-ERK signalling. The model is based upon a system of chemical reactions that describes cascade signalling dynamics. Using mass action kinetics, the chemical reactions are re-expressed as ordinary differential equations which are parameterised by in vitro data and solved numerically to obtain the temporal evolution of cascade component concentrations.

ResultsThe model provides a quantitative method to compute how dabrafenib and trametinib can be used in combination to synergistically inhibit ERK-activity in BRAFV600E-mutant melanoma cells. The model elucidates molecular mechanisms of vertical inhibition of the BRAFV600E-MEK-ERK cascade and delineates how elevated BRAF concentrations generate drug resistance to dabrafenib and trametinib. The computational simulations further suggest that elevated ATP levels could be a factor in drug resistance to dabrafenib.

ConclusionsThe model can be used to systematically motivate which dabrafenib-trametinib dose-combinations, for treating BRAFV600E-mutated melanoma, warrant experimental investigation. Mutations in coral soma and sperm imply lifelong stem cell renewal and cell lineage selection In many animals, the germline differentiates early in embryogenesis, so only mutations that accumulate in germ cells are inherited by offspring1. Exceptions to this developmental process may indicate that other mechanisms have evolved to limit the effects of deleterious mutation accumulation2. Stony corals are animals that can live for hundreds of years3 and have long been thought to produce gametes from somatic tissue4. To clarify conflicting evidence about germline-soma distinction in corals, we sequenced high coverage, full genomes with technical replicates for parent coral branches and their sperm pools. We identified post-embryonic single nucleotide variants (SNVs) unique to each parent branch, then checked if each SNV was shared by the respective sperm pool: 26% of post-embryonic SNVs were shared by the sperm but 74% were not. We also identified germline SNVs, those that were present in the sperm but not in the parent. These data suggest that self-renewing stem cells in corals differentiate into germ and soma throughout the adult life of the colony, with SNV rates and patterns differing markedly in stem, soma, and germ lineages. In addition to informing the important place in the evolutionary spectrum from non-Weismannian to Weismmanian animals that corals occupy, these insights inform how corals may generate adaptive diversity necessary in the face of global climate change. 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. Regulation of the PKD2 channel function by TACAN Autosomal dominant polycystic kidney disease (ADPKD) is caused by mutations in membrane receptor PKD1 or cation channel PKD2. TACAN (also named TMEM120A), recently reported as an ion channel in neuron cells for mechano and pain sensing, is also distributed in diverse non-neuronal tissues such as kidney, heart and intestine, suggesting its involvement in other functions. In this study, we found that TACAN is in complex with PKD2 in native renal cell lines. Using the two-electrode voltage clamp in Xenopus oocytes we found that TACAN inhibited the channel activity of PKD2 gain-of-function mutant F604P. The first and last transmembrane domains of TACAN were found to interact with the PKD2 C-and N-terminal portions, respectively. We showed that the TACAN N-terminus acted as a blocking peptide and that TACAN inhibits the PKD2 function through the PKD2/TACAN binding. By patch clamping in mammalian cells, we found that TACAN inhibits both the single channel conductance and open probability of PKD2 and mutant F604P. PKD2 co-expressed with TACAN, but not PKD2 alone, exhibited pressure sensitivity. Furthermore, we also found that TACAN aggravates PKD2-dependent tail curvature and pronephric cysts in larval zebrafish, in support of the in vitro inhibitory effects of TACAN. In summary, this study revealed that TACAN acts as a PKD2 inhibitor and mediates mechano sensitivity of the PKD2/TACAN channel complex. Balancing sensitivity and specificity in preclinical research The success of scientific discovery in preclinical drug development is based on the different roles of exploration and confirmation. Exploration requires sensitivity while confirmation needs to be specific. Based on published empirical data, we simulated two preclinical research trajectories to examine their effectiveness to detect true effects while keeping animal numbers low. We compare a conventional p-value based approach with a method based on an a priori determined smallest effect size of interest (SESOI). Using a SESOI increases transition rates from exploration to confirmation and results in higher detection rates across the trajectory. Specificity in the SESOI trajectory increases if the prior probability of true effects is low. We conclude that employing a SESOI is superior to a p-value based approach. Based on our findings, we propose a reconsideration of planning and conducting preclinical experiments, especially when the prior probability of true hypotheses is low. Hippocampal and medial prefrontal cortices encode structural task representations following progressive and interleaved training schedules Memory generalisations may be underpinned by either encoding- or retrieval-based generalisation mechanisms and different training schedules may bias some learners to favour one of these mechanisms over the other. We used a transitive inference task to investigate whether generalisation is influenced by progressive vs randomly interleaved training, and overnight consolidation. On consecutive days, participants learnt pairwise discriminations from two transitive hierarchies before being tested during fMRI. Inference performance was consistently better following progressive training, and for pairs further apart in the transitive hierarchy. BOLD pattern similarity correlated with hierarchical distances in the left hippocampus (HIP) and medial prefrontal cortex (MPFC) following both training schedules. These results are consistent with the use of structural representations that directly encode hierarchical relationships between task features. However, such effects were only observed in the MPFC for recently learnt relationships and were strongest in participants who were least proficient at the inference task. We conclude that humans preferentially employ encoding-based mechanisms to store map-like relational codes that can be used for memory generalisation. These codes are expressed in the HIP and MPFC following both progressive and interleaved training but are not sufficient for accurate inference. Critical re-evaluation of the slope factor of the operational model of agonism Background and PurposeAlthough being a relative term, agonist efficacy is a cornerstone in the proper assessment of agonist selectivity and signalling bias. The operational model of agonism (OMA) has become successful in the determination of agonist efficacies and ranking them. In 1985, Black et al. introduced the slope factor to the OMA to make it more flexible and allow for fitting steep as well as flat concentration-response curves.

Experimental ApproachFunctional analysis of OMA demonstrates that the slope factor implemented by Black et al. affects relationships among parameters of the OMA, leading to wrong estimates of operational efficacy and affinity of an agonist. In contrast, functional analysis of OMA with Hill coefficient as a slope factor has not this flaw.

Key ResultsFitting of the OMA with Black et al. slope factor to concentration-response curves of several models (homotropic allosteric modulation, substrate inhibition and non-competitive auto-inhibition) gives wrong estimates of operational efficacy and affinity. In contrast, fitting of the OMA modified by the Hill coefficient to concentration-response curves of these several models gives correct estimates of operational efficacy and affinity.

Conclusion and ImplicationsThe OMA modified by the Hill coefficient can be used for correct ranking of agonism and subsequent analysis, like quantification of signalling bias. Gradients in cellular metabolism regulate development of tonotopy along the cochlea. In vertebrates with elongated auditory organs, mechanosensory hair cells (HCs) are organised such that complex sounds are broken down into their component frequencies along the basal-to-apical long (tonotopic) axis. To generate the frequency-specific characteristics required at the appropriate positions, nascent HCs must determine their tonotopic properties during development. This relies on complex signalling between the developing HC and its local niche along the axis within the auditory organ. Here we apply NAD(P)H fluorescence lifetime imaging (FLIM) and live imaging of mitochondria to reveal metabolic gradients along the tonotopic axis of the developing cochlea. We further show that re-shaping these gradients during development, by inhibiting cytosolic glucose metabolism, alters normal Bmp7 and Chordin like-1 signalling leading to flattening of tonotopic patterning along the axis. Our work supports a causal link between morphogen signalling and metabolic reprogramming in specifying tonotopic morphologies in auditory HCs. Structurally and functionally distinct early antibody responses predict COVID-19 disease trajectory and mRNA vaccine response A damaging inflammatory response is strongly implicated in the pathogenesis of severe COVID-19 but mechanisms contributing to this response are unclear. In two prospective cohorts, early non-neutralizing, afucosylated, anti-SARS-CoV-2 IgG predicted progression from mild, to more severe COVID-19. In contrast to the antibody structures that predicted disease progression, antibodies that were elicited by mRNA SARS-CoV-2 vaccines were low in Fc afucosylation and enriched in sialylation, both modifications that reduce the inflammatory potential of IgG. To study the biology afucosylated IgG immune complexes, we developed an in vivo model which revealed that human IgG-Fc{gamma}R interactions can regulate inflammation in the lung. Afucosylated IgG immune complexes induced inflammatory cytokine production and robust infiltration of the lung by immune cells. By contrast, vaccine elicited IgG did not promote an inflammatory lung response. Here, we show that IgG-Fc{gamma}R interactions can regulate inflammation in the lung and define distinct lung activities associated with the IgG that predict severe COVID-19 and protection against SARS-CoV-2.

One Sentence SummaryDivergent early antibody responses predict COVID-19 disease trajectory and mRNA vaccine response and are functionally distinct in vivo. LPS-TLR4 pathway exaggerates alcohol-induced acute-on-chronic liver failure via provoking NETs formation Backgrounds: Intrahepatic infiltration of neutrophils is a character of alcoholic acute-on-chronic liver failure (AACLF) and neutrophil extracellular traps (NETs) are an important strategy for neutrophils to fix and kill invading microorganisms. Intestinal bacteria and the gut-liver axis have been thought to play a key role in many liver diseases also including AACLF. However, whether NETs appear in AACLF and play a role in AACLF is still unsure. Methods: WT, NE KO, and TLR4 KO mice were used to build the AACLF model, and the intestinal bacteria were eliminated at the same time and LPS was given. Then the formation of NETs and AACLF related markers were detected. Results: The serum MPO-DNA and LPS concentration was increased in AACLF patients and a correlation was revealed between these two indexes. More intrahepatic NETs formed in AACLF mice by testing MPO-DNA, Cit H3, and NE. These markers decreased with gut detergent and restored markers with gut detergent plus LPS supplement. While NETs formation failed to change with gut microbiome or combine LPS supplement in TLR4 KO mice. As we tested AACLF related characters, liver injury, intrahepatic fat deposition, inflammation, and fibrosis alleviated with depletion of NE. These related marks were also attenuated with gut sterilization by antibiotics and recovered with combined treatment with antibiotics plus LPS. But the liver injury, intrahepatic fat, fibro deposition, and liver inflammation-related markers did show a significant difference in TLR4 KO mice when they received the same treatment. Conclusion: Intestinal-derived LPS promotes NETs formation in AACLF through the TLR4 pathway and further accelerates the AACLF process by NETs. 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. Retraction ATPase motors from three orthologous type IVa pilus systems support promiscuous retraction of the Vibrio cholerae competence pilus Bacterial surface appendages called type IVa pili (T4aP) promote diverse activities including DNA uptake, twitching motility, and virulence. These activities rely on the ability of T4aP to dynamically extend and retract from the cell surface. Dynamic extension relies on a motor ATPase commonly called PilB. Most T4aP also rely on specific motor ATPases, commonly called PilT and PilU, to dynamically and forcefully retract. Here, we systematically assess whether motor ATPases from three orthologous T4aP can functionally complement Vibrio cholerae mutants that lack their endogenous motors. We found that the PilT and PilU retraction ATPases from the three T4aP systems tested are promiscuous and promote retraction of the V. cholerae competence T4aP despite a high degree of sequence divergence. In contrast, the orthologous extension ATPases from the same T4aP systems were not able to mediate extension of the V. cholerae competence T4aP despite exhibiting a similar degree of sequence divergence. Also, we show that one of the PilT orthologs characterized does not support PilU-dependent retraction and provide some data to indicate that the C-terminus of PilT is important for PilU-dependent retraction. Together, our data suggest that retraction ATPases may have maintained a high degree of promiscuity for promoting retraction of T4aP, while extension ATPases may have evolved to become specific for their cognate systems.

IMPORTANCEOne way that bacteria interact with their environments is via hair-like appendages called type IVa pili (T4aP). These appendages dynamically extend and retract from the cell surface via the action of distinct ATPase motors. T4aP are present in diverse bacterial species. Here, we demonstrate that retraction motors from three T4aP are promiscuous, and capable of promoting retraction of a heterologous T4aP system. By contrast, the extension ATPase motors from these same T4aP systems are specific and cannot promote extension of a heterologous T4aP. Thus, these results suggest that T4aP extension may be more tightly regulated compared to T4aP retraction. 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. A multi-state dynamic process confers mechano-adaptation to the bacterial flagellar motor Adaptation is a defining feature of living systems. The bacterial flagellar motor adapts to changes in the external mechanical environment by adding or removing torque-generating stator units. However, the molecular mechanism for mechanosensitive motor remodeling remains unclear. Here, we induced stator disassembly using electrorotation, followed by the time-dependent assembly of the individual stator units into the motor after electrorotation was terminated. From these experiments, we extracted detailed statistics of the dwell times that comprise the stochastic dynamics of the binding and unbinding of stator units. Dwell times reveal multiple timescales, indicating the existence of multiple binding states of the stator units. Based on these results, we propose a minimal model in which the stator unit can occupy four different states - two bound states with very different rates of unbinding, a diffusive unbound state, and a transiently detached state. Our minimal model quantitatively explains multiple features of the experimental data and allows us to determine the transition rates among all four states. Our experiments and modeling suggest a mechanism of mechano-adaptive remodeling of the bacterial flagellar motor in which torque generated by bound stator units controls their effective unbinding rate by modulating the transition between the two bound states. Furthermore, the binding rate of stator units with the motor has a non-monotonic dependence on the number of bound units, likely because of two counteracting effects of motor rotation on the binding process. A novel two-generation approach to the population dynamics of gregarious parasitoids Parasitoids, as important natural enemies, help maintain balance in natural ecosystems. Their population dynamics is generally predicted from the number of individuals. Here, using gregarious parasitoids as models, we show that this traditional approach omits one important parameter: the mothers manipulation of offspring fertility due to the clutch size-body size- fertility correlation. As a result of this correlation, when females deliberately adjust the number of offspring laid in a host, they determine not only the number but also the body sizes and reproductive potentials of those offspring. For the first time in the model species Anaphes flavipes, we determined the parasitoids offspring fertility from clutch size. Using this, we experimentally clarified the advantage of specific clutch size combinations and we show that identical fertility in the F1 generation can lead to distinctly different fertility values in the F2 generation. Even with the same number of hosts, lower fertility in the F1 generation can cause higher fertility in the F2 generation. Based on these results, we propose a novel two-generation approach which includes the clutch size-body size-fertility correlation. Our novel approach provides a new perspective for determining the individual fitness levels of gregarious parasitoids with new options for the modelling of parasitoid population dynamics. Data-Driven Network Models for Genetic Circuits From Time-Series Data with Incomplete Measurements Synthetic gene networks are frequently conceptualized and visualized as static graphs. This view of biological programming stands in stark contrast to the transient nature of biomolecular interaction, which is frequently enacted by labile molecules that are often unmeasured. Thus, the network topology and dynamics of synthetic gene networks can be difficult to verify in vivo or in vitro, due to the presence of unmeasured biological states. Here we introduce the dynamical structure function as a new mesoscopic, data-driven class of models to describe gene networks with incomplete measurements. We introduce a network reconstruction algorithm and a code base for reconstructing the dynamical structure function from data, to enable discovery and visualization of graphical relationships in a genetic circuit diagram as time-dependent functions rather than static, unknown weights. We prove a theorem, showing that dynamical structure functions can provide a data-driven estimate of the size of crosstalk fluctuations from an idealized model. We illustrate this idea with numerical examples. Finally, we show how data-driven estimation of dynamical structure functions can explain failure modes in two experimentally implemented genetic circuits, a historical genetic circuit and a new E. coli based transcriptional event detector.