Diverse classes of constraints enable broader applicability of a linear programming-based dynamic metabolic modeling framework Current metabolic modeling tools suffer from a variety of limitations, from scalability to simplifying assumptions, that preclude their use in many applications. We recently created a modeling framework, LK-DFBA, that addresses a key gap: capturing metabolite dynamics and regulation while retaining a potentially scalable linear programming structure. Key to this frameworks success are the linear kinetics and regulatory constraints imposed on the system. However, while the linearity of these constraints reduces computational complexity, it may not accurately capture the behavior of many biochemical systems. Here, we developed three new classes of LK-DFBA constraints to better model interactions between metabolites and the reactions they regulate. We tested these new approaches on several synthetic and biological systems, and also performed the first-ever comparison of LK-DFBA predictions to experimental data. We found that no single constraint approach was optimal across all systems examined, and systems with the same topological structure but different parameters were often best modeled by different types of constraints. However, we did find that the optimal constraint approach was generally robust to local perturbations of the system, indicating that just a single wild-type dataset could allow identification of the ideal constraint for a given system. These results suggest that the availability of multiple constraint approaches will allow LK-DFBA to model a wider range of metabolic systems.
E156G and Arg158, Phe-157/del mutation in NTD of spike protein in B.1.617.2 lineage of SARS-CoV-2 leads to immune evasion through antibody escape Emerging variants of SARS-CoV-2 with better immune escape mechanisms and higher transmissibility remains a persistent threat across the globe. B.1.617.2 (Delta) variant was first emerged from Maharashtra, India in December, 2020. This variant is classified to be a major cause and concern of the second wave of COVID-19 in India. In the present study, we explored the genomic and structural basis of this variant through computational analysis, protein modelling and molecular dynamics (MD) simulations approach. B.1.617.2 variant carried E156G and Arg158, Phe-157/del mutations in NTD of spike protein. These mutations in N-terminal domain (NTD) of spike protein of B.1.617.2 variant revealed more rigidity and reduced flexibility compared to spike protein of Wuhan isolate. Further, docking and MD simulation study with 4A8 monoclonal antibody which was reported to bind NTD of spike protein suggested reduced binding of B.1.617.2 spike protein compared to that of spike protein of Wuhan isolate. The results of the present study demonstrate the possible case of immune escape and thereby fitness advantage of the new variant and further warrants demonstration through experimental evidence. Our study identified the probable mechanism through which B.1.617.2 variant is more pathogenically evolved with higher transmissibility as compared to the wild-type.



O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=93 SRC="FIGDIR/small/447321v3_ufig1.gif" ALT="Figure 1">
View larger version (23K):
org.highwire.dtl.DTLVardef@b92cborg.highwire.dtl.DTLVardef@1d261f7org.highwire.dtl.DTLVardef@11da73eorg.highwire.dtl.DTLVardef@1cef6ca_HPS_FORMAT_FIGEXP  M_FIG C_FIG
Remarked Suppression of Aβ42 Protomer-protomer Dissociation Reaction via Pentamer Dimerization Elucidated by Molecular Dynamics Simulation Multimeric protein complexes are molecular apparatuses to regulate biological systems and often determine their fate. Among proteins forming such molecular assemblies, amyloid proteins have drawn attention over a half-century since amyloid fibril formation of these proteins is supposed to be a common pathogenic cause for neurodegenerative diseases. This process is triggered by the accumulation of fibril-like aggregates, while the microscopic mechanisms are mostly elusive due to technical limitation of experimental methodologies in individually observing each of diverse aggregate species in the aqueous solution. We then addressed this problem by employing atomistic molecular dynamics simulations for the paradigmatic amyloid protein, amyloid-{beta} (1-42) (A{beta}42). Seven different dimeric forms of oligomeric A{beta}42 fibril-like aggregate in aqueous solution, ranging from tetramer to decamer, were considered. We found additive effects of the size of these fibril-like aggregates on their thermodynamic stability and have clarified kinetic suppression of protomer-protomer dissociation reactions at and beyond the point of pentamer dimer formation. This observation was obtained from the specific combination of the A{beta}42 protomer structure and the physicochemical condition that we here examined, while it is worthwhile to recall that several amyloid fibrils take dimeric forms of their protomers. We could thus conclude that the stable formation of fibril-like protomer dimer should be involved in a turning point where rapid growth of amyloid fibrils is triggered. Induction of exaggerated cytokine production in human peripheral blood mononuclear cells by a recombinant SARS-CoV-2 spike glycoprotein S1 and its inhibition by dexamethasone An understanding of the pathological inflammatory mechanisms involved in SARS-CoV-2 virus infection is necessary in order to discover new molecular pharmacological targets for SARS-CoV-2 cytokine storm. In this study, the effects of a recombinant SARS-CoV-2 spike glycoprotein S1 was investigated in human peripheral blood mononuclear cells (PBMCs). Stimulation of PBMCs with spike glycoprotein S1 (100 ng/mL) resulted in significant elevation in the production of TNF, IL-6, IL-1{beta} and IL-8. However, pre-treatment with dexamethasone (100 nM) caused significant reduction in the release of these cytokines. Further experiments revealed that S1 stimulation of PBMCs increased phosphorylation of NF-{kappa}B p65 and I{kappa}B, and I{kappa}B degradation. DNA binding of NF-{kappa}B p65 was also significantly increased following stimulation with spike glycoprotein S1. Treatment of PBMCs with dexamethasone (100 nM) or BAY11-7082 (1 {micro}M) resulted in inhibition of spike glycoprotein S1-induced NF-{kappa}B activation. Activation of p38 MAPK by S1 was blocked in the presence of dexamethasone and SKF 86002. CRID3, but not dexamethasone pre-treatment produced significant inhibition of S1-induced activation of NLRP3/caspase-1. Further experiments revealed that S1-induced increase in the production of TNF, IL-6, IL-1{beta} and IL-8 was reduced in the presence of BAY11-7082 and SKF 86002, while CRID3 pre-treatment resulted in the reduction of IL-1{beta} production. These results suggest that SARS-CoV-2 spike glycoprotein S1 stimulated PBMCs to release pro-inflammatory cytokines through mechanisms involving activation of NF-{kappa}B, p38 MAPK and NLRP3 inflammasome. It is proposed that the clinical benefits of dexamethasone in COVID-19 is possibly due to its anti-inflammatory activity in reducing SARS-CoV-2 cytokine storm. Lower limb in sprinters with larger relative mass but not larger normalized moment of inertia Inhomogeneous muscle hypertrophy in response to daily motor executions is a morphological adaptation in humans. Inhomogeneous morphological adaptation would change the mass distribution, which in turn may affect the mechanical difficulty in moving from an inertia perspective; however, the morphological adaptation from the inertia perspective has not been examined. Here we show no corresponding differences in the normalized mass and normalized moment of inertia between the sprinters and the untrained non-sprinters. We analyzed the fat- and water-separated magnetic resonance images from the lower limbs of 11 male sprinters (100 m best time of 10.44-10.83 s) and 12 untrained non-sprinters. The lower-limb relative mass was significantly larger in sprinters (18.7{+/-}0.7% body mass) than in non-sprinters (17.6{+/-}0.6% body mass), while the normalized moment of inertia of the lower limb around the hip in the anatomical position was similar between them (0.044{+/-}0.002 vs. 0.042{+/-}0.002 [a. u.]). The thigh relative mass in sprinters (12.9{+/-}0.4% body mass) was significantly larger than that in non-sprinters (11.9{+/-}0.4% body mass), whereas the shank and foot relative masses were similar. We revealed that the adaptation to sprinting in humans does not relatively increase the mechanical difficulty to swing the lower limb in terms of inertia, even though the lower-limb mass increases, reflecting muscular hypertrophy. This is suggested to be one of the significant advantages of human morphology for adaptability to explosive motor tasks.

NEW & NOTEWORTHYWe compared the lower-limb inertia properties between male sprinters and untrained non-sprinters by using the fat- and water-separated magnetic resonance images. Sprinters had lower limbs with a larger relative mass but not a larger normalized moment of inertia than non-sprinters. We suggest one of the practical advantages of human morphological characteristics for their adaptation to explosive motor tasks, which provides a novel perspective to the literature for understanding human morphological adaptability. Unevolved proteins from modern and prebiotic amino acids manifest distinct structural profiles Natural proteins represent numerous but tiny structure/function islands in a vast ocean of possible protein sequences, most of which has not been explored by either biological evolution or research. Recent studies have suggested this uncharted sequence space possesses surprisingly high structural propensity, but development of an understanding of this phenomenon has been awaiting a systematic high-throughput approach.

Here, we designed, prepared, and characterized two combinatorial protein libraries consisting of randomized proteins, each 105 residues in length. The first library constructed proteins from the entire canonical alphabet of 20 amino acids. The second library used a subset of only 10 residues (A,S,D,G,L,I,P,T,E,V) that represent a consensus view of plausibly available amino acids through prebiotic chemistry. Our study shows that compact conformations resistant to proteolysis are (i) abundant (up to 40%) in random sequence space, (ii) independent of general Hsp70 chaperone system activity, and (iii) not granted solely by "late" and complex amino acid additions. The Hsp70 chaperone system effectively increases solubility and refoldability of the canonical alphabet but has only a minor impact on the "early" library. The early alphabet proteins are inherently more soluble and refoldable, possibly assisted by the cell-like environment in which these assays were performed.

Our work indicates that both early and modern amino acids are predisposed to supporting protein structure (either in forms of oligomers or globular/molten globule structures) and that protein structure may not be a unique outcome of evolution. Recombination marks the evolutionary dynamics of a recently endogenized retrovirus All vertebrate genomes have been colonized by retroviruses along their evolutionary trajectory. While endogenous retroviruses (ERVs) can contribute important physiological functions to contemporary hosts, such benefits are attributed to long-term co-evolution of ERV and host because germline infections are rare and expansion is slow, because the host effectively silences them. The genomes of several outbred species including mule deer (Odocoileus hemionus) are currently being colonized by ERVs, which provides an opportunity to study ERV dynamics at a time when few are fixed. Because we have locus-specific data on the distribution of cervid endogenous retrovirus (CrERV) in populations of mule deer, in this study we determine the molecular evolutionary processes acting on CrERV at each locus in the context of phylogenetic origin, genome location, and population prevalence. A mule deer genome was de novo assembled from short and long insert mate pair reads and CrERV sequence generated at each locus. CrERV composition and diversity have recently measurably increased by horizontal acquisition of a new retrovirus lineage. This new lineage has further expanded CrERV burden and CrERV genomic diversity by activating and recombining with existing CrERV. Resulting inter-lineage recombinants endogenized and subsequently retrotransposed. CrERV loci are significantly closer to genes than expected if integration were random and gene proximity might explain the recent expansion by retrotransposition of one recombinant CrERV lineage. Thus, in mule deer, retroviral colonization is a dynamic period in the molecular evolution of CrERV that also provides a burst of genomic diversity to the host population. Single amino-acid mutation in a Drosophila melanogaster ribosomal protein: an insight in uL11 transcriptional activity The ribosomal protein uL11 is located at the basis of the ribosome P-stalk and plays a paramount role in translational efficiency. In addition, no mutant for uL11 is available suggesting that this gene is haplo-insufficient as many other Ribosomal Protein Genes (RPGs). We have previously shown that overexpression of Drosophila melanogaster uL11 enhances the transcription of many RPGs and Ribosomal Biogenesis genes (RiBis) suggesting that uL11 might globally regulate the level of translation through its transcriptional activity. Moreover, uL11 trimethylated on lysine 3 (uL11K3me3) interacts with the chromodomain of the Enhancer of Polycomb and Trithorax Corto, and both proteins co- localize with RNA Polymerase II at many sites on polytene chromosomes. These data have led to the hypothesis that the N-terminal end of uL11, and more particularly the trimethylation of lysine 3, supports the extra-ribosomal activity of uL11 in transcription. To address this question, we mutated the lysine 3 codon using a CRISPR/Cas9 strategy and obtained several lysine 3 mutants. We describe here the first mutants of D. melanogaster uL11. Unexpectedly, the uL11K3A mutant, in which the lysine 3 codon is replaced by an alanine, displays a genuine Minute phenotype known to be characteristic of RPG deletions (longer development, low fertility, high lethality, thin and short bristles) whereas the uL11K3Y mutant, in which the lysine 3 codon is replaced by a tyrosine, is unaffected. In agreement, the rate of translation decreases in uL11K3A but not in uL11K3Y. Co-immunoprecipitation experiments show that the interaction between uL11 and the Corto chromodomain is impaired by both mutations. RNA- seq analysis from wing imaginal discs shows enrichment in the GO categories "glutathione metabolism" for up-regulated genes in both uL11K3A and uL11K3Y mutants and "regulation of transcription" for down-regulated genes in uL11K3A only. Analysis of the cis-regulatory sequences of these genes suggests that uL11 might regulate transcription of target genes in concert with the couple of transcription factors Mad/Med that mediate response to the Bone Morphogenetic Protein (BMP) signaling pathway. Citation needed? Wikipedia and the COVID-19 pandemic BackgroundWith the COVID-19 pandemics outbreak, millions flocked to Wikipedia for updated information. Amid growing concerns regarding an "infodemic", ensuring the quality of information is a crucial vector of public health. Investigating if and how Wikipedia remained up to date and in line with science is key to formulating strategies to counter misinformation. Using citation analyses, we asked: which sources informed Wikipedias COVID-19-related articles before and during the pandemics first wave (January-May 2020).

ResultsWe found that coronavirus-related articles referenced trusted media sources and high-quality academic research. Moreover, despite a surge in COVID-19 preprints, Wikipedia had a clear preference for open-access studies published in respected journals and made little use of preprints. Building a timeline of English COVID-19 articles from 2001-2020 revealed a nuanced trade-off between quality and timeliness. It further showed how preexisting articles on key topics related to the virus created a framework for integrating new knowledge. Supported by a rigid sourcing policy, this "scientific infrastructure" facilitated contextualization and regulated the influx of new information. Lastly, we constructed a network of DOI-Wikipedia articles, which showed the shifting landscape of pandemic-related knowledge on Wikipedia and how academic citations create a web of shared knowledge supporting topics like COVID-19 vaccine development.

ConclusionsUnderstanding how scientific research interacts with the digital knowledge-sphere during the pandemic provides insight into how Wikipedia can facilitate access to science. It also reveals how, aided by what we term its "citizen encyclopedists", it successfully fended off COVID-19 disinformation and how this unique model may be deployed in other contexts. A new method for the sampling and preservation of placental specimens in low-resource settings for the identification of P. falciparum and analysis of nucleic acids Collection, preservation, and shipment of histological specimens in low-resource settings is challenging. We present a novel method that achieved excellent preservation of placental specimens from rural Mali by using formalin fixation, ethanol dehydration, and long-term storage in a solar-powered freezer. Sample preservation success was 92%, permitting evaluation of current and past malaria infection, anemia, placental maturity, and inflammation. Using RNAscope(R) hybridization we were able to visualize cell-specific gene expression patterns in the formalin-fixed paraffin-embedded (FFPE) specimens. Additionally, our method entailed mirrored sampling from the two cut faces of a cotyledon, one for the FFPE workflows and the other for storage in RNAlater and RNA-seq. The isoforms of pyruvate kinase act as nutrient sensors for the β-cell KATP channel Insulin release from pancreatic {beta}-cells is triggered by the metabolism-dependent closure of ATP-sensitive K+ (KATP) channels. Here, we used isoform-specific deletion of pyruvate kinase (PK) to identify the unique roles of PKm1 and PKm2 in the {beta}-cell nutrient response. In the absence of PK activity, cytosolic ATP/ADP was insufficient to close KATP channels. Both PK isoforms were found on the plasma membrane, where they endow KATP channels with nutrient sensitivity and specificity. Although PKm1 accounts for >90% of cellular activity, PKm2 controls glucosedependent KATP closure. Amino acids, which generate mitochondrial phosphoenolpyruvate (PEP) without allosterically activating PKm2, depend on PKm1 to raise ATP/ADP and close KATP channels. Shifting the {beta}-cell from PKm2 to PKm1 correspondingly increased amino acid sensitivity and boosted secretory output. In addition to these membrane depolarizing "on-switches," {beta}-cell deletion of PCK2 exposed a mitochondrial PEP-dependent "off-switch" that facilitates Ca2+ extrusion. Restricting PCK2 prolonged Ca2+ influx and increased insulin secretion. In summary, it is not the KATP channel alone, but the PK-KATP microcompartment that confers nutrient sensitivity to the {beta}-cell plasma membrane. The differential response of the PK isoforms to the glycolytic and mitochondrial sources of PEP influences {beta}-cell nutrient preference, and controls the oscillatory cycle regulating insulin secretion. 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. Variation in female songbird state determines signal strength needed to evoke copulation It is the female response to male signals that determines courtship success. In most songbirds, females control reproduction via the copulation solicitation display (CSD), an innate, stereotyped posture produced in direct response to male displays. Because CSD can be elicited in the absence of males by the presentation of recorded song, CSD production enables investigations into the effects of underlying signal features and behavioral state on female mating preferences. Using computer vision to quantify CSD trajectory in female brown-headed cowbirds (Molothrus ater), we show that both song quality and a females internal state predict CSD production, as well as the onset latency and duration of the display. We also show that CSD can be produced in a graded fashion based on both signal strength and internal state. These results emphasize the importance of underlying receiver state in determining behavioral responses and suggest that female responsiveness acts in conjunction with male signal strength to determine the efficacy of male courtship. MS2AI: Automated repurposing of public peptide LC-MS data for machine learning applications MotivationLiquid-chromatography mass-spectrometry (LC-MS) is the established standard for analyzing the proteome in biological samples by identification and quantification of thousands of proteins. Machine learning (ML) promises to considerably improve the analysis of the resulting data, however, there is yet to be any tool that mediates the path from raw data to modern ML applications. More specifically, ML applications are currently hampered by three major limitations: (1) absence of balanced training data with large sample size; (2) unclear definition of sufficiently information-rich data representations for e.g. peptide identification; (3) lack of benchmarking of ML methods on specific LC-MS problems.

ResultsWe created the MS2AI pipeline that automates the process of gathering vast quantities of mass spectrometry (MS) data for large scale ML applications. The software retrieves raw data from either in-house sources or from the proteomics identifications database, PRIDE. Subsequently, the raw data is stored in a standardized format amenable for ML encompassing MS1/MS2 spectra and peptide identifications. This tool bridges the gap between MS and AI, and to this effect we also present an ML application in the form of a convolutional neural network for the identification of oxidized peptides.

AvailabilityAn open source implementation of the software can be found freely available for non-commercial use at https://gitlab.com/roettgerlab/ms2ai.

Contactveits@bmb.sdu.dk

Supplementary informationSupplementary data are available at Bioinformatics online. Frameshift and wild-type proteins are highly similar because the genetic code and genomes were optimized for frameshift tolerance Frameshift protein sequences encoded by alternative reading frames of coding genes have been considered meaningless, and frameshift mutations have been considered of little importance for the molecular evolution of coding genes and proteins. However, functional frameshifts have been found widely existing. It was puzzling how a frameshift protein kept its structure and functionality while its amino-acid sequence was changed substantially. Here we show that frame similarities between frameshifts and wild types are higher than random similarities and are defined at the genetic code, gene, and genome levels. In the standard genetic code, frameshift codon substitutions are more conservative than random substitutions. The frameshift tolerability of the standard genetic code ranks in the top 2.0-3.5% of alternative genetic codes, showing that the genetic code is nearly optimal for frameshift tolerance. Furthermore, frameshift-resistant codons (codon pairs) appear more frequently than expected in many genes and certain genomes, showing that the frameshift optimality is reflected not only in the genetic code but more importantly, in its allowance of further optimizing the frameshift tolerance of a particular gene or genome, which shed light on the role of frameshift mutations in molecular and genomic evolution. Fitness-related traits are maximized in recently introduced, slow-growing populations of a global invasive clam Many species are shifting their ranges being forced to rapidly respond to novel stressful environmental conditions. Colonizing individuals experience strong selective forces that favor the expression of life history traits notably affecting dispersal and reproductive rates in newly invaded habitats. Limited information is currently available on trait variation within the invasive range despite being critical for understanding ecological and evolutionary factors that drive the process of range expansion of invasive species. Here we evaluated life history shifts of the widely introduced Asian clam Corbicula, within its invaded range. Through an exhaustive literature search, we obtained data for 17 invasive Corbicula populations from different ecosystems worldwide. We tested the relationship between population and individual parameters relevant to the process of range expansion. Our main results are that recently introduced Corbicula populations were characterized by (i) low density and low rate of population increase, (ii) clams reproduced earlier in slow-growing populations, and (iii) density had no effect on population increase. All Corbicula populations analyzed in this study, which are fixed for one genotype (lineage Form A/R), experienced different selective environments in the introduced range. These findings support the perspective that adaptive phenotypic plasticity favored the expression of traits that maximize fitness in recently established populations, which faced stronger r-selective forces relative to long-established ones. We discuss the role of plasticity in facilitating rapid adaptation and increasing the likelihood of populations to overcome difficulties associated with low densities and low population increase in newly invaded areas. Herbivore life histories are altered by drought stress in their hosts plants Drought associated with climate change can stress plants, altering their interactions with phytophagous arthropods. Dryness not only impacts cultivated plants but also their parasites, which in some cases are favoured by drought. Herbivorous arthropods feeding on drought-stressed plants typically produce bigger offspring and develop faster. However, how much responses to abiotic factors differ among populations of a species remains poorly documented. Here, we evaluate variability among populations of a major agricultural pest, the two spotted spider mite, Tetranychus urticae, in response to drought-stress. We compare key life-history parameters of twelve populations that originate from climates ranging from wet and cool Atlantic locations to medium to dry hot Mediterranean locations. We evaluated how plant drought stress affects four life history traits: development time, fecundity, sex-ratio and emigration rate in an experiment comparing well-watered and drought-stressed bean plants. Mites feeding on drought-stressed plants developed faster and attempted to leave leaves less often, and young females were more fecund. The mites from wet temperate climates exhibited greater plasticity between the two water regimes than mites originating from dryer and hot climates. Thus, climate in the area of origin influences mite response to drought. Decoding personalized motor cortical excitability states from human electroencephalography Brain state-dependent transcranial magnetic stimulation (TMS) requires real-time identification of cortical excitability states. Current approaches deliver TMS during brain states that correlate with motor cortex (M1) excitability at the group level. Here, we hypothesized that machine learning classifiers could successfully discriminate between high and low M1 excitability states in individual participants using information obtained from low-density electroencephalography (EEG) signals. To test this, we analyzed a publicly available dataset that delivered 600 single TMS pulses to the right M1 during EEG and electromyography (EMG) recordings in 20 healthy adults. Multivariate pattern classification was used to discriminate between brain states during which TMS evoked small and large motor-evoked potentials (MEPs). Results show that personalized classifiers successfully discriminated between low and high M1 excitability states in 80% of tested participants. MEPs elicited during classifier-predicted high excitability states were significantly larger than those elicited during classifier-predicted low excitability states in 90% of tested participants. Personalized classifiers did not generalize across participants. Overall, results show that individual participants exhibit unique brain activity patterns which predict low and high M1 excitability states and that these patterns can be efficiently captured using low-density EEG signals. Our findings suggest that deploying individualized classifiers during brain state-dependent TMS may enable fully personalized neuromodulation in the future. The cyclic dinucleotide 2'3'-cGAMP induces a broad anti-bacterial and anti-viral response in the sea anemone Nematostella vectensis In mammals, cyclic dinucleotides (CDNs) bind and activate STING to initiate an anti-viral type I interferon response. CDNs and STING originated in bacteria and are present in most animals. By contrast, interferons are believed to have emerged in vertebrates; thus, the function of CDN signaling in invertebrates is unclear. Here, we use a CDN, 23-cGAMP, to activate immune responses in a model cnidarian invertebrate, the starlet sea anemone Nematostella vectensis. Using RNA-Seq, we found that 23-cGAMP induces robust transcription of both anti-viral and anti-bacterial genes in N. vectensis. Many of the anti-viral genes induced by 23-cGAMP are homologs of vertebrate interferon stimulated genes, implying that the interferon response predates the evolution of interferons. Knockdown experiments identified a role for NF-{kappa}B in specifically inducing anti-bacterial genes downstream of 23-cGAMP. Some of these putative anti-bacterial genes were also found to be induced during Pseudomonas aeruginosa infection. We characterized the protein product of one of the putative anti-bacterial genes, the N. vectensis homolog of Dae4, and found that it has conserved anti-bacterial activity. This work suggests that a broad anti-bacterial and anti-viral transcriptional response is an evolutionarily ancestral output of 23-cGAMP signaling in animals.

Significance statementCyclic dinucleotides are signaling molecules that originated in bacteria and were subsequently acquired and co-opted by animals for immune signaling. The major cyclic dinucleotide signaling pathway in mammals results in the production of anti-viral molecules called interferons. Invertebrates such as sea anemones lack interferons, and thus it was unclear whether cyclic dinucleotide signaling would play a role in immunity in these animals. Here we report that in the anemone Nematostella vectensis, cyclic dinucleotides activate both anti-viral and anti-bacterial immune responses, and do so through a conserved pathway. These results provide insights into the evolutionary origins of innate immunity, and suggest a broader ancestral role for cyclic dinucleotide signaling that evolved toward more specialized anti-viral functions in mammals. A high throughput lipidomics method using scheduled multiple reaction monitoring Lipid compositions of cells, tissues and bio-fluids are complex, with varying concentrations and structural diversity, which makes their identification challenging. Newer methods for comprehensive analysis of lipids are thus necessary. Herein, we propose a targeted-mass spectrometry based method for large-scale lipidomics using a combination of variable retention time window and relative dwell time weightage. Using this, we detected more than 1000 lipid species, including structural isomers. The limit of detection varied from femtomolar to nanomolar range and the coefficient of variance <30% for 849 lipid species. We used this method to identify lipids altered due to Vitamin B12 deficiency and found that the levels of lipids with {omega}-3 fatty acid chains decreased while those with {omega}-6 increased. This method enables identification of by far the largest number of lipid species with structural isomers in a single experiment and would significantly advance our understanding of the role of lipids in biological processes. OncoLoop: A network-based precision cancer medicine framework At present, prioritizing cancer treatments at the individual patient level remains challenging, and performing co-clinical studies using patient-derived models in real-time is often not feasible. To circumvent these challenges, we introduce OncoLoop, a precision medicine framework to predict and validate drug sensitivity in human tumors and their pre-existing high-fidelity (cognate) model(s) by leveraging perturbational profiles of clinically-relevant oncology drugs. As proof-of-concept, we applied OncoLoop to prostate cancer (PCa) using a series of genetically-engineered mouse models (GEMMs) that recapitulate a broad spectrum of disease states, including castration-resistant, metastatic, and neuroendocrine prostate cancer. Interrogation of published cohorts using Master Regulator (MR) conservation analysis revealed that most patients were represented by at least one cognate GEMM-derived tumor (GEMM-DT). Drugs predicted to invert MR activity in patients and their cognate GEMM-DTs were successfully validated, including in two cognate allografts and one cognate patient-derived xenograft (PDX). OncoLoop is a highly generalizable framework that can be extended to other cancers and potentially other diseases.

Significance StatementOncoLoop is a transcriptomic-based experimental and computational framework that can support rapid-turnaround co-clinical studies to identify and validate drugs for individual patients, which can then be readily adapted to clinical practice. This framework should be applicable in many cancer contexts for which appropriate models and drug perturbation data are available. Poly ADP-ribosylation of SET8 leads to aberrant H4K20 methylation domains in mammalian cells In mammalian cells, SET8 mediated Histone H4 Lys 20 monomethylation (H4K20me1) has been implicated in regulating mitotic condensation, DNA replication, DNA damage response, and gene expression. Here we show SET8, the only known enzyme for H4K20me1 is post-translationally poly ADP-ribosylated by PARP1 on lysine residues. PARP1 interacts with SET8 in a cell cycle-dependent manner. Poly ADP-ribosylation on SET8 renders it catalytically compromised, and degradation via ubiquitylation pathway. Knockdown of PARP1 led to an increase of SET8 protein levels, leading to aberrant H4K20me1 and H4K20me3 domains genome-wide. H4K20me1 is associated with higher gene transcription levels while the increase of H4K20me3 levels was predominant in DNA repeat elements correlating with modulation of chromatin states. Hence, SET8 mediated chromatin remodeling in mammalian cells are modulated by poly ADP-ribosylation by PARP1. Non-Newtonian Blood Rheology Impacts Left Atrial Stasis in Patient-Specific Simulations The lack of mechanically effective contraction of the left atrium (LA) during atrial fibrillation (AF) disturbs blood flow, increasing the risk of thrombosis and ischemic stroke. Thrombosis is most likely in the left atrial appendage (LAA), a small narrow sac where blood is prone to stagnate. Slow flow promotes the formation of erythrocyte aggregates in the LAA, also known as rouleaux, causing viscosity gradients that are usually disregarded in patient-specific simulations. To evaluate these non-Newtonian effects, we built atrial models derived from 4D computed tomography scans of patients and carried out computational fluid dynamics simulations using the Carreau-Yasuda constitutive relation. We examined six patients, three of whom had AF and LAA thrombosis or a history of transient ischemic attacks (TIAs). We modeled the effects of hematocrit and rouleaux formation kinetics by varying the parameterization of the Carreau-Yasuda relation and modulating non-Newtonian viscosity changes based on residence time. Comparing non-Newtonian and Newtonian simulations indicates that slow flow in the LAA increases blood viscosity, altering secondary swirling flows and intensifying blood stasis. While some of these effects can be subtle when examined using instantaneous metrics like shear rate or kinetic energy, they are manifested in the blood residence time, which accumulates over multiple heartbeats. Our data also reveal that LAA blood stasis worsens when hematocrit increases, offering a potential new mechanism for the clinically reported correlation between hematocrit and stroke incidence. In summary, we submit that hematocrit-dependent non-Newtonian blood rheology should be considered in calculating patient-specific blood stasis indices by computational fluid dynamics. A genomic assessment of the marine-speciation paradox within the toothed whale superfamily Delphinoidea The importance of post-divergence gene flow in speciation has been documented across a range of taxa in recent years, and may have been especially widespread in highly mobile, wide-ranging marine species, such as cetaceans. Here, we studied individual genomes from nine species across the three families of the toothed whale superfamily Delphinoidea (Delphinidae, Phocoenidae, Monodontidae). To investigate the role of post-divergence gene flow in the speciation process, we used a multifaceted approach, including: (i) phylogenomics, (ii) the distribution of shared derived alleles, and (iii) demographic inference. We found the divergence of lineages within Delphinoidea did not follow a process of pure bifurcation, but was much more complex. Sliding-window phylogenomics reveal a high prevalence of discordant topologies within the superfamily, with further analyses indicating these discordances arose due to both incomplete lineage sorting and gene flow. D-statistics, D-foil, and f-branch analyses supported gene flow between members of Delphinoidea, with the vast majority of gene flow occurring as ancient interfamilial events. Demographic analyses provided evidence that introgressive gene flow has likely ceased between all species pairs tested, despite reports of contemporary interspecific hybrids. Our study provides the first steps towards resolving the large complexity of speciation within Delphinoidea; we reveal the prevalence of ancient interfamilial gene flow events prior to the diversification of each family, and suggests that contemporary hybridisation events may be disadvantageous, as hybrid individuals do not appear to contribute to the parental species gene pools. Multiple brace root phenotypes promote anchorage and limit root lodging in maize Plant mechanical failure (lodging) causes global yield losses of 7-66% in cereal crops. We have previously shown that the above-ground nodal roots (brace roots) in maize are critical for anchorage. However, it is unknown how brace root phenotypes vary across genotypes and the functional consequence of this variation. This study quantifies the contribution of brace roots to anchorage, brace root traits, plant height, and root lodging susceptibility in 52 maize inbred lines. We show that the contribution of brace roots to anchorage and root lodging susceptibility varies among genotypes and this contribution can be explained by plant architectural variation. Additionally, supervised machine learning models were developed and show that multiple plant architectural phenotypes can predict the contribution of brace roots to anchorage and root lodging susceptibility. Together these data define the plant architectures that are important in lodging resistance and show that the contribution of brace roots to anchorage is a good proxy for root lodging susceptibility. Microbial biomass, composition, and functions are responsible for the differential removal of trace organic chemicals in biofiltration systems Biofiltration processes help to remove trace organic chemicals (TOrCs) both in wastewater and drinking water treatment systems. However, the detailed TOrCs biotransformation mechanisms as well as the underlying drivers behind the variability of site specific transformation processes remain elusive. In this study, we used laboratory batch incubations to investigate the biotransformation of 51 TOrCs in eight bioactive filter materials of different origins treating a range of waters, from wastewater effluents to drinking water. Microscopy, 16S rRNA amplicon and whole metagenome sequencing for assessing associations between the biotransformation rate constants, microbial composition and genetic potential complemented chemical analysis. We observed strong differences in the mean global removal of TOrCs between the individual sand filters (-1.4% to 58%), which were mirrored in overall biomass, microbial community composition, and enzyme encoding genes. From the six investigated biomass markers, ATP turned out to be a major predictor of the mean global biotransformation rate, while compound specific biotransformations were correlated with the microbial community composition. High biomass ecosystems were indicated in our systems by a dominance of Nitrospirae, but individual TOrC biotransformation was statistically connected to rare taxa (< 2%) such as Hydrogenophaga, or indiviudal functions such as the enoyl-CoA hydratase/3-hydroxyacyl-CoA dehydrogenase encoding genes. In general, this study provides new insights into so far rarely addressed variability of TOrCs biotransformation. We propose novel biological indicators for the removal performance of TOrCs in biofiltration systems, highlighting the role of living biomass in predicting and normalizing the global transformation, and the role of the microbial community for the individual transformation of TOrCs in engineered and natural systems.

Contribution to the Field StatementTrace organic chemicals (TOrCs) are an emerging problem in the aquatic environment that has attracted global attention over the last decade. Recent research efforts on this topic have increased our knowledge on the transformation of TOrCs and various technologies have been developed to improve their removal. In this study, we investigated a wide range of biotransformation of TOrCs by eight sand filter materials from wastewater and water treatment plants. Biotransformation rate constants were calculated using first-order kinetics to evaluate TOrC removal performance. We reevaluated the role of biomass and could thus explain a greater part of the global TOrC removal performance. The remaining variation in removal rates of individual compounds correlated with the microbiome of the biofilter. Rare biosphere lineages and specific enzyme categories genes were correlated with the removal of certain compounds. In summary, our research identified future indicators for successful biotransformation of TOrCs in biofilter systems.