Combinatorial drug-microenvironment interaction mapping reveals cell-extrinsic drug resistance mechanisms and clinically relevant patient subgroups in CLL The tumour microenvironment and genetic alterations collectively influence drug efficacy in cancer, but current evidence is limited to small scale studies and systematic analyses are lacking. We chose Chronic Lymphocytic Leukaemia (CLL), the most common leukaemia in adults, as a model disease to study this complex interplay systematically. We performed a combinatorial assay using 12 drugs individually co-applied with each of 17 microenvironmental stimuli in 192 primary CLL samples, generating a comprehensive map of drug-microenvironment interactions in CLL. This data was combined with whole-exome sequencing, DNA-methylation, RNA-sequencing and copy number variant annotation. Our assay identified four distinct CLL subgroups that differed in their responses to the panel of microenvironmental stimuli. These subgroups were characterized by distinct clinical outcomes independently of known prognostic markers. We investigated the effect of CLL- specific recurrent genetic alterations on microenvironmental responses and identified trisomy 12 as an amplifier of multiple microenvironmental stimuli. We further quantified the impact of microenvironmental stimuli on drug response, confirmed known interactions such as Interleukin (IL) 4 mediated resistance to B cell receptor (BCR) inhibitors, and identified new interactions such as Interferon-{gamma} induced resistance to BCR inhibitors. Finally, we identified interactions which were limited to genetic subgroups. Resistance to chemotherapeutics, such as Fludarabine, induced by Toll-Like Receptor (TLR) agonists could be observed in IGHV unmutated patient samples and IGHV mutated samples with trisomy 12. In-vivo relevance was investigated in CLL-infiltrated lymph nodes, which showed increased IL4 and TLR signalling activity compared to healthy samples (p<0.001). High IL4 activity in lymph nodes correlated with faster disease progression (p=0.038).

We provide a publicly available resource (www.dietrichlab.de/CLL_Microenvironment/) which uncovers tumour cell extrinsic influences on drug response and disease progression in CLL, and how these interactions are modulated by cell intrinsic molecular features.
Striatin 3 and MAP4K4 cooperate towards oncogenic growth and tissue invasion in medulloblastoma MAP4K4 has been correlated with increased cell motility and reduced proliferation in mammalian cells. Consequences and regulation of this dichotomous functionality of MAP4K4 in tumor cells remained elusive. We find that MAP4K4 interacts with the STRN3 protein in medulloblastoma and that STRN3 and MAP4K4 exert opposing functions in Hippo tumor suppressor signaling and clonal growth. However, depletion of either STRN3 or MAP4K4 reduces migration and tissue invasion, and concomitant loss of both proteins halts tumor cell growth in the cerebellar tissue. Mechanistically, STRN3 couples MAP4K4 to the protein phosphatase 2A, which inactivates growth repressing activities of MAP4K4. In parallel, STRN3 enables growth factor-induced PKC{theta} activation and the direct phosphorylation of VASPS157 by MAP4K4, which are necessary for efficient migration and tissue invasion. Thus, STRN3 controls MAP4K4 functionality in growth factor-activated tumor cells, and the cooperation of the two proteins promotes invasiveness through nPKCs and VASP and increases proliferation via inactivation of Hippo signaling. Targeting of the STRN3-MAP4K4 cooperation could reactivate Hippo signaling and repress invasiveness in growth factor-driven tumors.

Graphical abstract

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=139 SRC="FIGDIR/small/442906v2_ufig1.gif" ALT="Figure 1">
View larger version (24K):
org.highwire.dtl.DTLVardef@102b3c1org.highwire.dtl.DTLVardef@6753ceorg.highwire.dtl.DTLVardef@1af2ff5org.highwire.dtl.DTLVardef@759127_HPS_FORMAT_FIGEXP  M_FIG C_FIG
R7 photoreceptor axon targeting depends on the relative levels of lost and found expression in R7 and its synaptic partners As neural circuits form, growing processes select the correct synaptic partners through interactions between cell surface proteins. The presence of such proteins on two neuronal processes may lead to either adhesion or repulsion; however, the consequences of mismatched expression have rarely been explored. Here we show that the Drosophila CUB-LDL protein Lost and found (Loaf) is required in the UV-sensitive R7 photoreceptor for normal axon targeting only when Loaf is also present in its synaptic partners. Although targeting occurs normally in loaf mutant animals, removing loaf from photoreceptors or expressing it in their postsynaptic neurons Tm5a/b or Dm9 in a loaf mutant causes mistargeting of R7 axons. Loaf localizes primarily to intracellular vesicles including endosomes. We propose that Loaf regulates the trafficking or function of one or more cell surface proteins, and an excess of these proteins on the synaptic partners of R7 prevents the formation of stable connections. 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. The interplay of chromatin phase separation and lamina interactions in nuclear organization The genetic material of eukaryotes is segregated into transcriptionally active euchromatin and silent heterochromatin compartments. The spatial arrangement of chromatin compartments evolves over the course of cellular life in a process that remains poorly understood. The latest nuclear imaging experiments reveal a number of dynamical signatures of chromatin which are reminiscent of active multi-phase liquids. This includes the observations of viscoelastic response, coherent motions, Ostwald ripening and coalescence of chromatin compartments. There is also growing evidence that liquid-liquid phase separation of protein and nucleic acid components is the underlying mechanism for the dynamical behavior of chromatin. In order to dissect the organizational and dynamical implications of chromatins liquid behavior, we have devised a phenomenological field-theoretic model of nucleus as a multi-phase condensate of liquid chromatin types. Employing the liquid chromatin model of Drosophila nucleus, we have carried out an extensive set of simulations with an objective to shed light on the dynamics and chromatin patterning observed in the latest nuclear imaging experiments. Our simulations reveal the emergence of experimentally detected mesoscale chromatin channels and spheroidal droplets which arise from the dynamic interplay of chromatin type to type interactions and intermingling of chromosomal territories. We also quantitatively reproduce coherent motions of chromatin domains observed in displacement correlation spectroscopy measurements which are explained within the framework of our model by phase separation of chromatin types operating within constrained intra-chromosomal and inter-chromosomal boundaries. Finally, we illuminate the role of heterochromatin-lamina interactions in the nuclear organization by showing that these interactions enhance the mobility of euchromatin and indirectly introduce correlated motions of heterochromatin droplets.

SIGNIFICANCE STATEMENTThe latest super-resolution imaging experiments have revealed a surprisingly dynamic and stochastic nature of chromatin in eukaryotic nuclei which is reminiscent of multi-phase fluid behavior. As a result, the concepts from the theory of complex fluids such as phase separation, viscoelasticity, and droplet nucleation have found widespread utility in understanding salient features of nuclear organization. In order to understand and disentangle the complex interplay of forces that contribute to the emergent patterns of organization and dynamics, we have devised a phenomenological field-theoretic model of nucleus as a multi-phase condensate of liquid chromatin types. Armed with a mesoscopic model of nuclear chromatin, we have shed light on the distinct dynamical and structural contributions of chromatin type interactions intermingling of chromosomal territories and lamina binding. We also shed light on the dynamical heterogeneity and coherent motions of chromatin domains which are fully captured by an interplay of micro-phase separation of chromatin types and lamina binding. Deletion of Rap-Phr systems in Bacillus subtilis influences in vitro biofilm formation and plant root colonization Natural isolates of the soil-dwelling bacterium Bacillus subtilis form robust biofilms under laboratory conditions and colonize plant roots. B. subtilis biofilm gene expression displays phenotypic heterogeneity that is influenced by a family of Rap-Phr regulatory systems. Most Rap-Phr systems in B. subtilis have been studied independently, in different genetic backgrounds and under distinct conditions, hampering true comparison of the Rap-Phr systems impact on bacterial differentiation. Here, we investigated each of the 12 Rap-Phr systems of B. subtilis NCIB 3610 for their effect on biofilm formation. By studying single {Delta}rap-phr mutants, we show that despite redundancy between the cell-cell communication systems, deletion of each of the 12 Rap-Phr systems influences matrix gene expression, which could possibly enable fine-tuning of the timing and level of matrix production in response to specific conditions. Furthermore, some of the {Delta}rap-phr mutants demonstrated altered biofilm formation in vitro and colonization of Arabidopsis thaliana roots, but not necessarily similarly in both processes, indicating that the pathways regulating matrix gene expression and other factors important for biofilm formation may be differently regulated under these distinct conditions.

Significance StatementNatural isolates of Bacillus subtilis form robust biofilms in vitro and on plant roots. The formation of these heterogeneous populations is regulated by diverse Rap-Phr systems. However, most Rap-Phr systems in B. subtilis have been studied independently and in different genetic backgrounds. Here, we report that all 12 Rap-Phr systems affect matrix gene expression, while some of them affect development of in vitro biofilms and plant root colonization. Our study highlights the importance of the Rap-Phr systems in environmental adaptation of B. subtilis, specifically during biofilm formation in the rhizosphere. Maternal microbiota Bifidobacterium promotes placental morphogenesis, nutrient transport and fetal growth in mice The gut microbiota plays a central role in regulating host metabolism. However, while substantial progress has been made in discerning how the microbiota influences host functions post birth and beyond, little is known about how key members of the maternal gut microbiota can influence feto-placental growth. Here, using germ-free and specific-pathogen-free mice, we demonstrate that the bacterium Bifidobacterium breve UCC2003 modulates maternal body adaptations, placental structure and nutrient transporter capacity, with implications for fetal metabolism and growth. The effects of B. breve UCC2003 on feto-placental growth are mediated, in part, by changes in the maternal and placental metabolome (i.e. acetate and carnitine). Histological analysis of the placenta confirmed that Bifidobacterium modifies placental structure via changes in Igf2P0, Dlk1, Mapk1 and Mapk14 expression. Additionally, B. breve UCC2003, acting through Slc2a1 and Fatp1-4 transporters, was shown to restore fetal glycaemia and fetal growth in association with changes in the fetal hepatic transcriptome. This study provides knowledge towards a novel and safe therapeutic strategy for treating pregnancy disorders via modulation of the maternal gut microbiota. Heterogeneous mutation rates and spectra in yeast hybrids Mutation rates and spectra vary between species and among populations. Hybridization can contribute to this variation, but its role remains poorly understood. Estimating mutation rates requires controlled conditions where the effect of natural selection can be minimized. One way to achieve this is through mutation accumulation experiments coupled with genome sequencing. Here we investigate 400 mutation accumulation lines initiated from 11 genotypes spanning intra-lineage, inter-lineage and interspecific crosses of the yeasts Saccharomyces paradoxus and S. cerevisiae and propagated for 770 generations. We find significant differences in mutation rates and spectra among crosses, which are not related to the level of divergence of parental strains but are specific to some genotype combinations. Differences in number of generations and departures from neutrality play a minor role, whereas polyploidy and loss of heterozygosity impact mutation rates in some of the hybrid crosses in an opposite way. HOW CONCERNING IS A SARS-COV-2 VARIANT OF CONCERN? COMPUTATIONAL PREDICTIONS AND THE VARIANTS LABELING SYSTEM We herein report a study to evaluate the use of computational prediction of SARS-CoV-2 genetic variations in improving the current variants labeling system. First, we reviewed the basis of the system developed by the World Health Organization (WHO) for the labeling of SARS-CoV-2 genetic variants and the adaptations made to it by the United States Center of Diseases Control (CDC). We observed that the labeling system is based upon the virus major attributes. However, we found that the labeling criteria of the SARS-CoV-2 variants derived from these attributes are not accurately defined and are used differently by the two health management agencies. Consequently, discrepancies exist between the labels given by WHO and CDC to same variants. Our observations suggest that giving the VOC label to a new variant is premature and might not be appropriate. Therefore, we carried out a comparative computational study to predict the effects of the mutations on the virus structure and functions of five VOCs. By linking these data to the criteria used by WHO and the CDC for variant labeling, we ascertained that comparative computational predictions of the impact of genetic variations are a better ground for rapid and more accurate labelling of SARS-CoV-2 variants. We propose to label all emergent variants VUM or VBM and to carry out computational predictive studies and thorough variants comparison, upon which more appropriate and informative labels can be attributed. Furthermore, harmonization of the variants labeling system would be globally beneficial to communicate about and fight COVID-19 pandemic. Adipose triglyceride lipase promotes prostaglandin-dependent actin remodeling by regulating substrate release from lipid droplets A key factor controlling oocyte quality and fertility is lipids. Even though lipid droplets (LDs) are crucial regulators of lipid metabolism, their roles in fertility are poorly understood. During Drosophila oogenesis, LD accumulation in nurse cells coincides with dynamic actin remodeling necessary for late-stage follicle morphogenesis and fertility. Loss of the LD-associated Adipose Triglyceride Lipase (ATGL) disrupts both actin bundle formation and cortical actin integrity, an unusual phenotype also seen when Pxt, the enzyme responsible for prostaglandin (PG) synthesis, is missing. Dominant genetic interactions and PG treatment of follicles in vitro reveal that ATGL and Pxt act in the same pathway to regulate actin remodeling, with ATGL upstream of Pxt. Further, lipidomic analysis detects arachidonic acid (AA) containing triglycerides in ovaries. Because AA is the substrate for Pxt, we propose that ATGL releases AA from LDs to drive PG synthesis necessary for follicle development. We also find that exogenous AA is toxic to follicles in vitro, and LDs modulate this toxicity. This leads to the model that LDs both sequester AA to limit toxicity, and release AA via ATGL to drive PG production. We speculate that the same pathways are conserved across organisms to regulate oocyte development and promote fertility. Improved detection of evolutionary selection highlights potential bias from different sequencing strategies in complex genomic-regions Balancing selection occurs when multiple alleles are kept at elevated frequencies in equilibrium due to opposing evolutionary pressures. A new statistical method was developed to test for selection using efficient Bayesian techniques. Selection signals in three different data sets, generated with variable sequencing technologies, were compared: clinical trios, HLA NGS typed samples, and whole-genome long-read samples. Genome-wide, selection was observed across multiple gene families whose biological functions favor diversification, revealing established targets as well as 45 novel genes under selection. Using high-resolution HLA typing and long-read sequencing data, for the characterization of the MHC, revealed strong selection in expected peptide-binding domains as well as previously understudied intronic and intergenic regions of the MHC. Surprisingly, SIRPA, demonstrated dramatic selection signal, second only to the MHC in most settings. In conclusion, employing novel statistical approaches and improved sequencing technologies is critical to properly analyze complex genomic regions. Autotaxin is a Potential Link between Genetic Risk Factors and Pathogenesis of Systematic Lupus Erythematosus in Plasmacytoid Dendritic Cells BackgroundThe importance of autotaxin, an enzyme that catalyzes lysophospholipid production, has recently been recognized in various diseases, including cancer and autoimmune diseases. Herein, we examined the role of autotaxin in systemic lupus erythematosus (SLE), utilizing data from ImmuNexUT, a comprehensive database consisting of transcriptome data and expression quantitative trait locus (eQTL) data of immune cells from patients with immune-mediated disorders.

MethodsSerum autotaxin concentrations in patients with SLE and healthy controls (HCs) were compared. The transcriptome data of patients with SLE and age-and sex-matched HCs were obtained from ImmuNexUT. The expression of ENPP2, the gene encoding autotaxin, was examined in peripheral blood immune cells. Next, weighted gene correlation network analysis (WGCNA) was performed to identify genes with expression patterns similar to ENPP2. The ImmuNexUT eQTL database and public epigenomic databases were used to infer the relationship between autotaxin and pathogenesis of SLE.

ResultsAutotaxin levels were elevated in the serum of patients with SLE compared to HCs. Furthermore, the expression of ENPP2 was higher in plasmacytoid dendritic cells (pDCs) than in other immune cell subsets, and its expression was elevated in pDCs of patients with SLE compared to HCs. In WGCNA, ENPP2 belonged to a module that correlated with disease activity. This module was enriched in interferon-associated genes and included genes whose expression was influenced by single-nucleotide polymorphisms associated with SLE, suggesting that it is a key module connecting genetic risk factors of SLE with disease pathogenesis. Analysis utilizing the ImmuNexUT eQTL database and public epigenomic databases suggested that the increased expression of ENPP2 in pDCs from patients with SLE may be caused by increased expression of interferon-associated genes and increased binding of STAT3 complexes to the regulatory region of ENPP2.

ConclusionsAutotaxin may play a critical role in connecting genetic risk factors of SLE to disease pathogenesis in pDCs. Intron-mediated induction of phenotypic heterogeneity Introns are universally present in the nuclear genomes of eukaryotes1. The budding yeast, an otherwise intron-poor species, preserves two sets of ribosomal protein (RP) genes differing primarily in their introns2-4. Despite recent findings on the role of RP introns under stress and starvation5-7, understanding the contribution of introns to ribosome regulation remains challenging. Here, combining isogrowth profiling8 with single-cell protein measurements9, we found that introns can mediate inducible phenotypic heterogeneity conferring a clear fitness advantage. Osmotic stress leads to bimodal expression of the small ribosomal subunit protein Rps22B, mediated by an intron in the 5 untranslated region of its transcript. The two resulting yeast subpopulations differ in their ability to cope with starvation. Low Rps22B protein levels resulted in prolonged survival under sustained starvation, while high Rps22B levels enabled cells to grow faster after transient starvation. Further, yeast growing at high sugar concentrations - similar to those in ripe grapes - exhibit bimodal Rps22B expression when approaching stationary phase. Differential intron-mediated regulation of RP genes thus provides a way to diversify the population when starvation looms in natural environments. Our findings reveal a new role for introns in inducing phenotypic heterogeneity in changing environments and suggest that duplicated RP genes in yeast contribute to resolving the evolutionary conflict between precise expression control and environmental responsiveness10. Long-term conservation effects of protected areas in stochastic population dynamics Terrestrial and marine protected areas are essential tools in mitigating anthropogenic impacts and promoting population persistence and resource sustainability. Adequately implemented protected areas (PAs) promote long-term conservation effects. Stochasticity causes fluctuations in the conservation effects of PAs, and so it is important to investigate the variabilities of these conservation effects to inform their long-term conservation effects. To investigate long-term conservation effects, I develop and analyze new models of stochastic processes that encompass the fluctuations generated by demographic or environmental stochasticity in PAs management. The stochastic model is built upon individual processes. In the model, density-independent mortality, migration between PAs and non-PAs, and site preferences characterize the features of the PA. The effect of PAs size is also examined. The long-term conservation effects are quantified using the coefficient of variation (CV) of population size in PAs, where a lower CV indicates higher robustness in stochastic variations. Typically, the results from this study demonstrate that sufficiently reduced density-independent mortality in PAs and high site preference and immigration rate of PA are likely to decrease the CV. However, different types of stochasticity induce rather different consequences: under demographic stochasticity, the CV is always reduced when PAs increase the population size therein, but an increased population size by PAs does not always decrease the CV under environmental stochasticity. The deterministic dynamics of the model are investigated, facilitating effective management decisions. Second Heart Field-derived Cells Contribute to Angiotensin II-mediated Ascending Aortopathies BackgroundThe ascending aorta is a common location for aneurysm and dissection. This aortic region is populated by a mosaic of medial and adventitial cells that are embryonically derived from either the second heart field (SHF) or the cardiac neural crest. SHF-derived cells populate areas that coincide with the spatial specificity of thoracic aortopathies. The purpose of this study was to determine whether and how SHF-derived cells contribute to ascending aortopathies.

MethodsAscending aortic pathologies were examined in patients with sporadic thoracic aortopathies and angiotensin II (AngII)-infused mice. Ascending aortas without overt pathology from AngII-infused mice were subjected to mass spectrometry assisted proteomics, and molecular features of SHF-derived cells were determined by single cell transcriptomic analyses. Genetic deletion of either low-density lipoprotein receptor-related protein 1 (Lrp1) or transforming growth factor-{beta} receptor 2 (Tgfbr2) in SHF- derived cells was conducted to examine the impact of SHF-derived cells on vascular integrity.

ResultsPathologies in human ascending aortic aneurysmal tissues were predominant in outer medial layers and adventitia. This gradient was mimicked in mouse aortas following AngII infusion that was coincident with the distribution of SHF-derived cells. Proteomics indicated that brief AngII infusion, prior to overt pathology, evoked downregulation of SMC proteins and differential expression of extracellular matrix proteins, including several LRP1 ligands. LRP1 deletion in SHF-derived cells augmented AngII-induced ascending aortic aneurysm and rupture. Single cell transcriptomic analysis revealed that brief AngII infusion decreased Lrp1 and Tgfbr2 mRNA abundance in SHF-derived cells and induced a unique fibroblast population with low abundance of Tgfbr2 mRNA. SHF-specific Tgfbr2 deletion led to embryonic lethality at E12.5 with dilatation of the outflow tract and retroperitoneal hemorrhage. Integration of proteomic and single cell transcriptomics results identified plasminogen activator inhibitor 1 (PAI1) as the most increased protein in SHF-derived SMCs and fibroblasts during AngII infusion. Immunostaining revealed a transmural gradient of PAI1 in both ascending aortas of AngII-infused mice and human ascending aneurysmal aortas that mimicked the gradient of medial and adventitial pathologies.

ConclusionSHF-derived cells exert a critical role in maintaining vascular integrity through LRP1 and TGF-{beta} signaling associated with increases of aortic PAI1.



O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=89 SRC="FIGDIR/small/930917v6_ufig1.gif" ALT="Figure 1">
View larger version (34K):
org.highwire.dtl.DTLVardef@185e3d9org.highwire.dtl.DTLVardef@131d702org.highwire.dtl.DTLVardef@11afa42org.highwire.dtl.DTLVardef@242786_HPS_FORMAT_FIGEXP  M_FIG C_FIG Clinical PerspectiveO_ST_ABSWhat is new?C_ST_ABS- SHF-derived SMCs and fibroblasts associate with AngII-induced aortic pathologies.
- AngII induces a distinct fibroblast sub-cluster that is less abundant for mRNAs related to major extracellular components and TGF{beta} ligands and receptors, but more abundant for proliferative genes.
- TGFBR2 deletion in SHF-derived cells are embryonic lethal with significant dilatation of the outflow tract in mice.
- SHF-specific deletion of LRP1 leads to aortic pathologies in mice, supporting the importance of SHF-derived cells in maintaining ascending aortic wall integrity.


What are the clinical implications?- Heterogeneity of the embryonic origins of SMCs and fibroblasts contributes to complex mechanisms of vasculopathy formation, which should be considered when investigating the pathogenesis of thoracic aortopathies. A pilot study on the home range and movement patterns of the Andean Fox Lycalopex culpaeus (Molina 1782) in Cotopaxi National Park, Ecuador This study reports movement patterns and home range estimates of an Andean fox (Lycalopex culpaeus) in Cotopaxi National Park in Ecuador, representing the first GPS-tagging of the species. The GPS functioned well during the 197-day tracking period. Home range sizes ranged between 4.9 - 8.1 km2, depending on the estimation method. Movement speeds averaged 0.17 km/hr at day vs. 0.87 km/hr. at night, and distance traveled averaged 0.23 km at day vs. 0.89 km at night. These preliminary results highlight the importance of collecting unbiased, high-quality data which enables an enhanced understanding on mammal behavior and human/animal interaction. Preprints in motion: tracking changes between preprint posting and journal publication during a pandemic Amidst the COVID-19 pandemic, preprints in the biomedical sciences are being posted and accessed at unprecedented rates, drawing widespread attention from the general public, press and policymakers for the first time. This phenomenon has sharpened longstanding questions about the reliability of information shared prior to journal peer review. Does the information shared in preprints typically withstand the scrutiny of peer review, or are conclusions likely to change in the version of record? We assessed preprints from bioRxiv and medRxiv that had been posted and subsequently published in a journal through 30th April 2020, representing the initial phase of the pandemic response. We utilised a combination of automatic and manual annotations to quantify how an article changed between the preprinted and published version. We found that the total number of figure panels and tables changed little between preprint and published articles. Moreover, the conclusions of 7.2% of non-COVID-19-related and 17.2% of COVID-19-related abstracts undergo a discrete change by the time of publication, but the majority of these changes do not qualitatively change the conclusions of the paper. A Novel Paradigm for Deep Reinforcement Learning of Biomimetic Systems Mechanisms behind neural control of movement have been an active area of research. Goal-directed movement is a common experimental setup used to understand these mechanisms and relevant neural pathways. On the one hand, optimal feedback control theory is used to model and make quantitative predictions of the coordinated activations of the effectors, such as muscles, joints or limbs. While on the other hand, evidence shows that higher centres such as Basal Ganglia and Cerebellum are involved in activities such as reinforcement learning and error correction. In this paper, we provide a framework to build a digital twin of relevant sections of the human spinal cord using our NEUROiD platform. The digital twin is anatomically and physiologically realistic model of the spinal cord at cellular, spinal networks and system level. We then build a framework to learn the supraspinal activations necessary to perform a simple goal directed movement of the upper limb. The NEUROiD model is interfaced to an Opensim model for all the musculoskeletal simulations. We use Deep Reinforcement Learning to obtain the supraspinal activations necessary to perform the goal directed movement. As per our knowledge, this is the first time an attempt is made to learn the stimulation pattern at the spinal cord level, especially by limiting the observation space to only the afferent feedback received on the Ia, II and Ib fibers. Such a setup results in a biologically realistic constrained environment for learning. Our results show that (1) Reinforcement Learning algorithm converges naturally to the triphasic response observed during goal directed movement (2) Increasing the complexity of the goal gradually helped to accelerate learning (3) Modulation of the afferent inputs were sufficient to execute tasks which were not explicitly learned, but were closely related to the learnt task. Out-of-Distribution Generalization from Labeled and Unlabeled Gene Expression Data for Drug Response Prediction Data discrepancy between preclinical and clinical datasets poses a major challenge for accurate drug response prediction based on gene expression data. Different methods of transfer learning have been proposed to address this data discrepancy. These methods generally use cell lines as source domains and patients, patient-derived xenografts, or other cell lines as target domains. However, they assume that they have access to the target domain during training or fine-tuning and they can only take labeled source domains as input. The former is a strong assumption that is not satisfied during deployment of these models in the clinic. The latter means these methods rely on labeled source domains which are of limited size. To avoid these assumptions, we formulate drug response prediction as an out-of-distribution generalization problem which does not assume that the target domain is accessible during training. Moreover, to exploit unlabeled source domain data, which tends to be much more plentiful than labeled data, we adopt a semi-supervised approach. We propose Velodrome, a semi-supervised method of out-of-distribution generalization that takes labeled and unlabeled data from different resources as input and makes generalizable predictions. Velodrome achieves this goal by introducing an objective function that combines a supervised loss for accurate prediction, an alignment loss for generalization, and a consistency loss to incorporate unlabeled samples. Our experimental results demonstrate that Velodrome outperforms state-of-the-art pharmacogenomics and transfer learning baselines on cell lines, patient-derived xenografts, and patients. Finally, we showed that Velodrome models generalize to different tissue types that were well-represented, under-represented, or completely absent in the training data. Overall, our results suggest that Velodrome may guide precision oncology more accurately. Design, Mutate, Screen: High-throughput creation of genetic clocks with different period-amplitude characteristics Directed evolution has become an invaluable tool in protein engineering and has also greatly influenced the construction of synthetic gene circuits. The ability to generate diversity at precise targets for directed evolution approaches has improved vastly, allowing researchers to create large, specific mutant libraries with relative ease. Screening approaches for large mutant libraries have similarly come a long way, especially when the desired behavior can easily be tested for with static, single time-point assays. For more complex gene circuits with dynamic phenotypes that change over time, directed evolution approaches to controlling and tuning circuit behavior have been hindered by the lack of sufficiently high-throughput screening methods to isolate variants with desired characteristics. Here we utilize directed mutagenesis and multiplexed microfluidics to develop a workflow for creating, screening and tuning dynamic gene circuits that operate at the population level. Specifically, we create a mutant library of an existing oscillator, the synchronized lysis circuit, and tune its dynamics while uncovering principles regarding its behavior. Lastly, we utilize this directed evolution workflow to construct a new synchronized genetic oscillator that exhibits robust dynamics over long time scales. Inter-trial Variations in EEG predict the individual differences in social tasks People experience events and form an impression of others in a way that is affected by social influence every day. In the present study, we designed a series of tasks centered on social influence to investigate peoples bias in following others opinions and its underlying neural predictors. Our results showed social conformity and proved that social influence-induced change can be predicted by the amount of inter-trial electroencephalogram (EEG) variations when people view others faces. This prediction effect is robust in the alpha-band over the right frontal and left occipital electrodes for negative influence. Inter-trial EEG variations can also predict the subsequent trust difference between negatively and positively influenced faces. Our findings suggest that higher Inter-trial EEG variations in the pre-influence task may serve as a predictor of high vulnerability to social influence. The present study provides a novel approach that considers both the stability of ones endogenous EEG and the variations in external task components to predict human social behaviors. 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. Complex population history of admixed sticklebacklineages Introgressive hybridization is an important process contributing to adaptation but challenging to identify, undermining the efforts to understand its role in evolution. On the other hand, many analytical methods assume direct descent from a single common ancestor, and multiple waves of colonisation and admixture can violate their assumptions and lead to seriously biased results. A detailed analysis of 888 whole-genome-sequenced nine-spined sticklebacks revealed a complex pattern of population ancestry involving multiple waves of colonization and introgression across northern Europe. The secondary contact zone between the two recognized lineages was found to be wider than anticipated, displaying a smooth gradient of foreign ancestry with some curious exceptions attributable to freshwater isolates preserving footprints of prehistoric colonisation events. The results provide a back-drop against which introgression and its role in local adaptation can be studied and highlight methodological challenges associated with hidden admixture events. Phytoene synthase 2 in tomato fruits remains functional and contributes to abscisic acid formation In ripening tomato fruits, the leaf-specific carotenoids biosynthesis mediated by phytoene synthase 2 (PSY2) is replaced by a fruit-specific pathway by the expression of two chromoplast-specific genes: phytoene synthase 1 (PSY1) and lycopene-{beta}-cyclase (CYCB). Though both PSY1 and PSY2 genes express in tomato fruits, the functional role of PSY2 is not known. To decipher whether PSY2-mediated carotenogenesis operates in ripening fruits, we blocked the in vivo activity of lycopene-{beta}-cyclases in fruits of several carotenoids and ripening mutants by CPTA (2-(4-Chlorophenylthio)triethylamine hydrochloride), an inhibitor of lycopene-{beta}-cyclases. The CPTA-treatment induced accumulation of lycopene in leaves, immature-green and ripening fruits. Even, in psy1 mutants V7 and r that are deficient in fruit-specific carotenoid biosynthesis, CPTA triggered lycopene accumulation but lowered the abscisic acid level. Differing from fruit-specific carotenogenesis, CPTA-treated V7 and r mutant fruits accumulated lycopene but not phytoene and phytofluene. The lack of phytoene and phytofluene accumulation was reminiscent of PSY2-mediated leaf-like carotenogenesis, where phytoene and phytofluene accumulation is never seen. The lycopene accumulation was associated with the partial transformation of chloroplasts to chromoplasts bearing thread-like crystalline structures. Our study uncovers the operation of a parallel carotenogenesis pathway mediated by PSY2 that provides precursors for abscisic acid biosynthesis in ripening tomato fruits.

Significance statementIt is believed that in ripening tomato fruits phytoene synthase 2 that drives carotenoid biosynthesis in leaves is redundant. Contrary to this, we show that in phytoene synthase 1 mutant fruit that is bereft of lycopene, the chemical inhibition of lycopene {beta}-cyclases triggers lycopene accumulation. Our results uncover that phytoene synthase 2 remains functional in ripening fruits and provides precursors for abscisic acid formation in fruits. The pH Dependence of Niclosamide Solubility, Dissolution, and Morphology Motivates Potentially Universal Mucin-Penetrating Nasal and Throat Sprays for COVID19, its Contagious Variants, and Other Respiratory Viral Infections MotivationWith the coronavirus pandemic still raging, prophylactic nasal and early treatment throat sprays could help prevent infection and reduce viral load. Niclosamide has the potential to treat a broad range of viral infections if local bioavailability is optimized as mucin-penetrating solutions as opposed to more traditional microparticle-based sprays that cannot penetrate the mucin.

ExperimentalpH-dependence of supernatant concentrations and dissolution rates of niclosamide were measured in buffered solutions by Nanodrop-UV/Vis-spectroscopy for niclosamide from different suppliers, as precipitated material, and as cosolvates. Data was compared to predictions from Henderson Hasselbalch and precipitation pH models. Optimal microscopy was used to observe the morphologies of precipitated and converted niclosamide.

ResultsSupernatant-concentrations of niclosamide increased with increasing pH, from 1.77uM at pH 3.66 to 30uM at pH 8, and more rapidly from 90uM at pH8.5 to 300uM at pH 9.1, reaching 641uM at pH 9.5. Logarithmic rates for dissolution increased by [~]3x for pHs 8.62 to 9.44. However, when precipitated from supersaturated solution, niclosamide equilibrated to much lower final supernatant concentrations, reflective of more stable polymorphs at each pH that were also apparent for niclosamide from other suppliers and cosolvates.

ConclusionsGiven niclosamides activity against COVID19, its more contagious variants, and other respiratory viral infections, these niclosamide solutions, that put the virus in lockdown, could represent universal prophylactic nasal and early treatment throat sprays. As solutions they would be the simplest and potentially most effective formulations from both an efficacy standpoint as well as manufacturing and distribution, with no cold chain. They now just need testing. Correlational networking guides the discovery of cryptic natural product biosynthetic enzymes Bacterial natural product biosynthetic genes, canonically clustered, have been increasingly found to rely on cryptic enzymes encoded elsewhere in the genome for completion of biosynthesis. The study and application of lanthipeptides are frequently hindered by cryptic proteases required for final maturation. Here, we establish a global correlation network bridging the gap between lanthipeptide precursors and cryptic proteases. Applying our analysis to 161,954 bacterial genomes, we establish 5,209 correlations between precursors and cryptic proteases, with 91 prioritized. We use network predictions and co-expression analysis to reveal a previously missing protease for the maturation of class I lanthipeptide paenilan. We further discover widely distributed bacterial M16B metallopeptidases of previously unclear biological function as a new family of lanthipeptide proteases. We show the involvement of a pair of bifunctional M16B proteases in the production of novel class III lanthipeptides with high substrate specificity. Together, these results demonstrate the strength of our correlational networking approach to the discovery of cryptic lanthipeptide proteases and potentially other missing enzymes for natural products biosynthesis.

Graphical abstract O_FIG_DISPLAY_L [Figure 1] M_FIG_DISPLAY C_FIG_DISPLAY