Evolutionary analysis of THAP9 transposase: conserved regions, novel motifs THAP9 is a transposable element-derived gene that encodes the THAP9 protein, which is homologous to the Drosophila P-element transposase (DmTNP) and can cut and paste DNA. However, the exact functional role of THAP9 is unknown. Here, we perform evolutionary analysis and extensive in silico characterization of THAP9, including predicting domains and putative post-translational modification sites. We predict previously unreported mammalian-specific post-translational modification sites that may play a role in the subcellular localization of THAP9. We also observe that although THAP9 has evolved under a strong pervasive purifying selection, yielding high conservation of THAP9, there are distinct class-specific conservation patterns of key functional residues in certain domains. Furthermore, investigation of THAP9 expression profiles in various cancer and matched normal datasets demonstrated underexpression and overexpression in testicular cancers and thymic epithelial tumors, respectively, thus suggesting a possible role of THAP9 in cancer.
DeLUCS: Deep Learning for (Unsupervised) Clustering of DNA Sequences We present a novel Deep Learning method for the Unsupervised Clustering of DNA Sequences (DeLUCS) that does not require sequence alignment, sequence homology, or (taxonomic) identifiers. DeLUCS uses Frequency Chaos Game Representations (FCGR) of primary DNA sequences, and generates "mimic" sequence FCGRs to self-learn data patterns (genomic signatures) through the optimization of multiple neural networks. A majority voting scheme is then used to determine the final cluster assignment for each sequence. The clusters learned by DeLUCS match true taxonomic groups for large and diverse datasets, with accuracies ranging from 77% to 100%: 2,500 complete vertebrate mitochondrial genomes, at taxonomic levels from sub-phylum to genera; 3,200 randomly selected 400 kbp-long bacterial genome segments, into clusters corresponding to bacterial families; three viral genome and gene datasets, averaging 1,300 sequences each, into clusters corresponding to virus subtypes. DeLUCS significantly outperforms two classic clustering methods (K-means++ and Gaussian Mixture Models) for unlabelled data, by as much as 47%. DeLUCS is highly effective, it is able to cluster datasets of unlabelled primary DNA sequences totalling over 1 billion bp of data, and it bypasses common limitations to classification resulting from the lack of sequence homology, variation in sequence length, and the absence or instability of sequence annotations and taxonomic identifiers. Thus, DeLUCS offers fast and accurate DNA sequence clustering for previously intractable datasets.
The mechanism of full activation of tumor suppressor PTEN at the phosphoinositide-enriched membrane Tumor suppressor PTEN, the second most highly mutated protein in cancer, dephosphorylates signaling lipid PIP3 produced by PI3Ks. Excess PIP3 promotes cell proliferation. The mechanism at the membrane of this pivotal phosphatase is unknown hindering drug discovery. Exploiting explicit solvent simulations, we tracked full-length PTEN trafficking from the cytosol to the membrane. We observed its interaction with membranes composed of zwitterionic phosphatidylcholine, anionic phosphatidylserine, and phosphoinositides, including signaling lipids PIP2 and PIP3. We tracked its moving away from the zwitterionic and getting absorbed onto anionic membrane that harbors PIP3. We followed it localizing on microdomains enriched in signaling lipids, as PI3K does, and observed PIP3 allosterically unfolding the N-terminal PIP2 binding domain, positioning it favorably for the polybasic motif interaction with PIP2. Finally, we determined PTEN catalytic action at the membrane, all in line with experimental observations, deciphering the mechanisms of how PTEN anchors to the membrane and restrains cancer. 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. Linking null models to hypotheses to improve nestedness analysis O_LINestedness is a common pattern in metacommunities and interaction networks, whose causes are still subject to investigation. Nestedness inference is challenging because, beyond calculating an index, we need to compare observed and expected values using null models. There are many different null models and the choice between them strongly affects test outcomes. Furthermore, there is no theoretical basis to guide this choice.
C_LIO_LIHere, we propose a different look at the meaning of nestedness that improves our understanding of its causes and directly links null models to hypotheses. The nestedness of a matrix is the combination of marginal sum inequality and high overlap. The higher the overlap, the more predictable the cell values by marginal sums. Therefore, the only way a process can increase nestedness is by promoting marginal sum inequalities, without concomitantly introducing preferences. We propose a combination of null models to test for different topological hypotheses. In these null models, interactions are randomized following a probability matrix. The probability matrix of the equiprobable model excludes marginal sum inequalities and consequently all nestedness-generating mechanisms, providing a distribution of expected values for nestedness significance tests. The probability matrix of the proportional model, by conserving marginal sums and maximizing overlap, conserves nestedness-generating mechanisms and excludes nestedness-disrupting mechanisms, thereby yielding the expected nestedness values for fully nested matrices.
C_LIO_LIWe advance a new protocol to test for nestedness in binary and weighted matrices, in which those two null models represent different topological hypotheses. Our protocol helps to tell apart cases in which nestedness is different from expected by chance and cases in which nestedness represents the best topological archetype. We also evaluate the efficiency of applying this protocol using some common nestedness indices. Analyses based on weighted data proved to be highly efficient and informative, while binary analyses induced severe distortions. Finally, we empirically illustrate our approach using a plant-pollinator network.
C_LIO_LIThrough a shift of perspective, our approach reconciliates contradictions in null model analysis and provides a guideline for nestedness interpretation.
C_LI Nuclear alpha-synuclein is present in the human brain and is modified in dementia with Lewy bodies Dementia with Lewy bodies is pathologically defined by the cytoplasmic accumulation of alpha-synuclein within neuronal cells in the brain. Alpha-synuclein is predominately pre-synaptic, but has been reported present in various subcellular compartments in cell and animal models. In particular, nuclear alpha-synuclein is evident in-vitro and in disease models and has been associated with altered DNA integrity, gene transcription, nuclear homeostasis. However, owing to various factors, the presence of alpha-synuclein in the nuclei of human brain cells remains controversial, as does its role in synucleinopathies.

Here, we close this gap and provide a unique demonstration confirming the presence of nuclear alpha-synuclein in post-mortem brain tissue obtained from cases of dementia with Lewy bodies as well as from controls via immunohistochemistry, immunoblot, and label-free mass-spectrometry.

Discrete intra-nuclear alpha-synuclein puncta reactive against phosphorylated serine 129-alpha-synuclein and pan-alpha-synuclein antibodies were observed in cortical neurons and non-neuronal cells in fixed brain sections and in isolated nuclear preparations from Dementia with Lewy bodies cases and matched controls. Subsequent biochemical analysis of subcellular fractionated tissue confirmed alpha-synuclein as present in a nuclear fraction at levels ~ 10-fold lower than in the cytoplasm. Critically, however, an increase in monomeric nuclear alpha-synuclein phosphorylated as serine 129 was observed in cases of dementia with Lewy bodies alongside higher molecular weight pan- and phosphorylation reactive alpha-synuclein species, consistent with the formation of intranuclear phosphorylated alpha-synuclein oligomers. Furthermore, the presence of nuclear alpha-synuclein was confirmed via label free mass spectrometry, as 6 unique alpha-synuclein derived peptide sequences were identified in nuclear fractions (71.4% sequence coverage).

Collectively, our data confirm the presence of nuclear alpha-synuclein in human brain tissue and describe nuclear pathology associated with dementia with Lewy bodies. These findings address a major controversy in the synucleinopathy field by confirming the presence of nuclear alpha-synuclein in autoptic human brain tissue and, for the first time, identify that alpha-synuclein is aggregated into novel and potentially pathological assemblies in the nucleus as part of the disease process associated with dementia with Lewy bodies and thus may contribute to the disease phenotype. PTEN regulates adipocyte progenitor growth, differentiation and replicative aging The tumor suppressor phosphatase and tensin homolog (PTEN) negatively regulates the insulin signaling pathway. Germline PTEN pathogenic variants cause PTEN Hamartoma Tumor Syndrome (PHTS), associated with lipoma development in children. It remains unclear which mechanisms trigger this aberrant adipose tissue growth. Adipocyte progenitor cells (APCs) lose their capacity to differentiate into adipocytes during continuous culture, while APCs from PHTS patients lipomas retain their adipogenic potential over a prolonged period. To investigate the role of PTEN in adipose tissue development we performed functional assays and RNA sequencing of control and PTEN knockdown APCs. Reduction of PTEN levels using siRNA or CRISPR lead to an enhanced proliferation and differentiation of APCs. FOXO1 was downregulated on the mRNA level while inactivation through phosphorylation increased. FOXO1 phosphorylation initiates the expression of the lipogenesis activating transcription factor SREBP1. SREBP1 levels were higher after PTEN knockdown and may account for the enhanced adipogenesis. To validate this we overexpressed constitutively active FOXO1 in PTEN CRISPR cells and found reduced adipogenesis, accompanied by a SREBP1 downregulation. We observed that PTEN levels were upregulated during long term culture of wild type APCs. PTEN CRISPR cells showed less senescence compared to controls and the senescence marker CDKN1A (p21) was downregulated in PTEN knockdown cells. Cellular senescence was the most significantly enriched pathway found in RNA sequencing of PTEN knockdown vs. control cells. These results provide evidence that PTEN is involved in the regulation of APCs proliferation, differentiation and senescence, thereby contributing to aberrant adipose tissue growth in PHTS patients. Artificial seawater based long-term culture of colonial ascidians Tunicates are highly diverse marine invertebrate filter-feeders that are vertebrates closest relatives. These organisms, despite a drastically different body plan during their adulthood, have a tissue complexity related to that of vertebrates. Ascidians, which compose most of the Tunicata, are benthic sessile hermaphrodites that reproduce sexually through a motile tadpole larval stage. Over half of the known ascidians species are able to reproduce asexually by budding, typically leading to the formation of colonies where animals, called zooids, are interconnected through an external vascular system. In addition, colonial ascidians are established models for important biological processes including allorecognition, immunobiology, aging, angiogenesis and whole-body regeneration. However, the current paucity in breeding infrastructures limits the study of these animals to coastal regions.

To promote a wider scientific spreading and popularity of colonial ascidians, we have developed a flexible recirculating husbandry setup for their long-term in-lab culture. Our system is inspired both by the flow-through aquariums used by coastal ascidian labs, as well as by the recirculating in-lab systems used for zebrafish research. Our hybrid system thus combines colony breeding, water filtering and food culturing in a semi-automated system where specimens develop on hanging microscopy glass slides. Temperature, light/dark cycles, flow speed and feeding rates can be controlled independently in four different breeding environments to provide room for species-specific optimization as well as for running experiments. This setup is complemented with a quarantine for the acclimatization of wild isolates.

Herein we present our success in breeding Botrylloides diegensis, a species of colonial ascidians, for more than 3 years in recirculating artificial seawater over 600 km away from their natural habitat. We show that colonies adapt well to in-lab culturing provided that a suitable marine microbiome is present, and that a specific strain can be isolated, propagated and efficiently used for research over prolonged periods of time. The flexible and modular structure of our system can be scaled and adapted to the needs of specific species, such as Botryllus schlosseri, as well as of particular laboratory spaces. Overall, we show that Botrylloides diegensis can be proficiently bred in-land and suggest that our results can be extended to other species of colonial ascidians to promote research on these fascinating animals.

HighlightsO_LIFirst artificial seawater based recirculating aquaculture for colonial ascidians
C_LIO_LIOver 3 years of continuous breeding
C_LIO_LISemi-automated setup with minimized maintenance
C_LIO_LIGood biomass production for strain propagation
C_LIO_LI4 different culture conditions for optimized breeding for species of interest
C_LI The compound of baicalein, wogonin, and oroxylin-A inhibits EMT in A549 cell line via PI3K/AKT-TWIST1-glycolysis pathway Non-small cell lung cancer (NSCLC) is a worldwide disease with high morbidity and mortality, which is most derived from its metastasis. Some researches show that the epithelial-mesenchymal transition (EMT) process promotes lung cancer cells migration and invasion, leading to NSCLC metastasis. Total Flavonoid Aglycones Extract (TFAE) isolated from Scutellaria baicalensis was reported to inhibit tumor growth and induce apoptosis. In this study, we found that baicalein, wogonin, and oroxylin-A were the active compounds of TFAE. After reconstructing with these three compounds (baicalein (65.8%), wogonin (21.2%), and oroxylin-A (13.0%)), the reconstructed TFAE (reTFAE) present inhibitory effect on the EMT process of A549 cells. Then, bioinformatic technology was employed to elucidate the potential pharmacodynamic mechanism network of reTFAE. We identified the relationship between reTFAE and PI3K/Akt signaling pathways, with TWIST1 as the key protein. LY294002, the inhibitor of PI3K/Akt signaling pathway, and knock-down TWIST1 could significantly enhance the efficacy of reTFAE, with increasing expression of epithelial markers and decreasing expression of mesenchymal markers in A549 cells at the same time. Furthermore, stable isotope dimethyl-labeled proteomics technology was conducted to complement the follow-up mechanism that the EMT-inhibition process may be realized through glycolysis pathway. In conclusion, we claim that TWIST1-targeted flavonoids could provide a new strategy to inhibit EMT progress for the treatment of NSCLC. 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. Boosting L-type Ca2+ channel activity in tuberous sclerosis mitigates excess glutamate receptors Tuberous sclerosis complex (TS) is a dominant, multisystem disorder with devastating neurological symptoms. Approximately 85% of TS patients suffer from epilepsy over their lifespan and roughly 25-50% of those patients develop Autism Spectrum Disorder (1, 2). Current seizure therapies are effective in some, but not all, and often have significant risk factors associated with their use (1, 3). Thus, there is a critical need for new medication development or drug repositioning. Herein, we leveraged proteomic signatures of epilepsy and ASD, often comorbid in TS, to utilize an in silco approach to identify new drug therapies for TS-related seizures. We have discovered that activation of L-type voltage dependent calcium channels (L-VDCC) by Bay-K8644 (BayK) in a preclinical mouse model of TS rescues the excess expression of ionotropic, AMPA-subtype glutamate (GluA) receptors. As added proof of BayK working through L-VDCC to regulate GluA levels, we found that increasing expression of alpha2delta2 (2{delta}2), an auxiliary calcium channel subunit that boosts L-VDCC surface expression, similarly lowers the surface expression of dendritic GluA in TS. These BayK-induced molecular alterations may potentially improve the quality of life of humans suffering from TS.

Significance StatementCausal mechanisms of Tuberous Sclerosis (TS)-associated neurological disorders are under-characterized and treatment options are lacking. Using a computational approach of mTOR/DJ-1 target mRNAs to predict new medications, we report that boosting L-type voltage-dependent Ca2+ channel (L-VDCC) activity in a preclinical TS mouse model that exhibits a deficit in dendritic L-VDCC activity ameliorates key molecular pathologies that are predicted to underlie seizures. Restoring the mTOR/DJ-1 pathway upstream of L-VDCC, therefore, may serve as a new therapeutic avenue to mitigate seizures and mortality in TS. Cellular-level phenotyping of tumor-immune microenvironment (TiME) in patients in vivo reveals distinct inflammation and endothelial anergy signatures Immunotherapies have shown unprecedented clinical benefits in several malignancies1-3. However, clinical responses remain variable and unpredictable, indicating the need to develop predictive platforms that can improve patient stratification4. Phenotyping of tumors into hot, altered, or cold5 based on assessment of only T-cell infiltration in static tumor biopsies provides suboptimal prediction of immunotherapy response6,7. In vivo dynamic mechanisms within the tumor microenvironment such as tumor angiogenesis and leukocyte trafficking5,8,9 also play a central role in modulating anti-tumor immunity and therefore immunotherapy response. Here, we report novel tumor immune microenvironment (TiME) phenotyping in vivo in patients with non-invasive spatially-resolved cellular-level imaging based on endogenous contrast. Investigating skin cancers as a model, with reflectance confocal microscopy (RCM) imaging10, we determined four major phenotypes with variable prevalence of vasculature (Vasc) and inflammation (Inf) features: VaschiInfhi, VaschiInflo, VascloInfhi and Vascmed/hiInflo. The VaschiInfhi phenotype correlates with high immune activation, exhaustion, and vascular signatures while VaschiInflo with endothelial anergy and immune exclusion. Automated quantification of TiME features demonstrates moderate-high accuracy and correlation with corresponding gene expression. Prospectively analyzed response to topical immunotherapy show highest response in VascloInfhi, and reveals the added value of vascular features in predicting treatment response. Our novel in vivo cellular-level imaging and phenotyping approach can potentially advance our fundamental understanding of TiME, develop robust predictors for immunotherapy outcomes and identify novel targetable pathways in future. Modulation of P2X4 pore closure by magnesium,potassium, and ATP The P2X4 receptor plays a prominent role in cellular responses to extracellular ATP. Through classical all-atom molecular dynamics (MD) simulations totaling 24 {micro}s we have investigated how metal-complexed ATP stabilizes the channels open state and prevents its closing. We have identified two metal-binding sites, magnesium (Mg2+) and potassium (K+), one at the intersection of the three subunits in the ectodomain (MBS1) and the second one near the ATP binding site (MBS2), similar to those characterized in Gulf coast P2X. Our data indicate that when Mg2+ and K+ ions are complexed with ATP, the channel is locked into an open state. Interestingly, irrespective of the number of bound ATP molecules, Mg2+ ions bound to the MBS2 resisted collapsing of the open state protein to a closed state by stabilizing the ATP-protein interactions. However, when Mg2+ in the MBS2 was replaced with K+ ions, as might be expected when in equilibrium with an extracellular solution, the interactions between the subunits were weakened and we found evidence of pore collapse. This collapse was apparent when fewer than two ATP were bound to MBS2 in the presence of K+. Therefore, the different capacities of common cations to stabilize the channel may underlie a mechanism governing P2X4 channel gating in physiological systems. This study provides structural insights into the differential modulation of ATP activation of P2X4 by Mg2+ and K+. Actomyosin activity-dependent apical targeting of Rab11 vesicles reinforces apical constriction During tissue morphogenesis, cell shape changes resulting from cell-generated forces often require active regulation of intracellular trafficking. How mechanical stimuli influence intracellular trafficking and how such regulation impacts tissue mechanics are not fully understood. In this study, we identify an actomyosin dependent mechanism involving Rab11- mediated trafficking in regulating apical constriction in the Drosophila embryo. During Drosophila mesoderm invagination, apical actin and Myosin II (actomyosin) contractility induces apical accumulation of Rab11-marked vesicle-like structures ("Rab11 vesicles") by promoting a directional bias in dynein mediated vesicle transport. At the apical domain, Rab11 vesicles are enriched near the adherens junctions (AJs). The apical accumulation of Rab11 vesicles is essential to prevent fragmented apical AJs, breaks in the supracellular actomyosin network and a reduction in the apical constriction rate. This Rab11 function is separate from its role in promoting apical Myosin II accumulation. These findings suggest a feedback mechanism between actomyosin activity and Rab11-mediated intracellular trafficking that regulates the force generation machinery during tissue folding. Increased susceptibility to Chrysanthemum Yellows phytoplasma infection in Atcals7 ko plants is accompanied by enhanced expression of carbohydrate transporters. Callose deposition around sieve pores, under control of callose synthase 7 (AtCALS7), has been interpreted as a mechanical response to limit pathogen spread in phytoplasma-infected plants. Wild-type and Atcals7ko mutants were therefore employed to unveil the mode of involvement of CALS7 in the plants response to phytoplasma infection. The fresh weights of healthy and CY-(Chrysanthemum Yellows) phytoplasma-infected Arabidopsis wild-type and mutant plants indicated two superimposed effects of the absence of CALS7: a partial impairment of photo-assimilate transport and a stimulated phytoplasma proliferation as illustrated by a significantly increased phytoplasma titre in Atcal7ko mutants. Further studies solely dealt with the effects of CALS7 absence on phytoplasma growth. Phytoplasma infection affected sieve-element substructure to a larger extent in mutants than in wild-type plants, which was also true for the levels of some free carbohydrates. Moreover, infection induced a similar upregulation of gene expression of enzymes involved in sucrose cleavage (AtSUS5, AtSUS6) and transmembrane transport (AtSWEET11) in mutants and wild-type plants, but an increased gene expression of carbohydrate transmembrane transporters (AtSWEET12, AtSTP13, AtSUC3) in infected mutants only. It remains still unclear how the absence of AtCALS7 leads to gene upregulation and how an increased intercellular mobility of carbohydrates and possibly effectors contributes to a higher susceptibility. It is also unclear if modified sieve-pore structures in mutants allow a better spread of phytoplasmas giving rise to higher titre.

Author SummaryPhytoplasma infections are one of the most limiting factors for production of important crops all over the world. Phytoplasma disease epidemics can be handled mainly by insect-vector control using insecticides. Basic information about plant-phytoplasma interactions are still limited, nevertheless it is necessary to design new management and breeding strategies aimed to obtain more tolerant or resistant cultivars. Phytoplasmas are obligate intracellular parasites restricted to the phloem sieve tubes. Callose deposition at the sieve plates has been described since the  70s as a mechanical defence process to limit pathogen spread by occluding sieve pores. Studies reported 40 years later demonstrated that callose at the sieve pores is also involved in sieve-pore development and function and, hence, in mass-flow regulation, carbohydrate metabolism and distribution, and plant growth. Here, we reported on the role(s) of sieve-element callose in phytoplasma-infected Arabidopsis, using a mutant lacking AtCALS7, the enzyme responsible for callose synthesis in the sieve elements. The results indicate that loss of AtCAL7 appears to confer increased susceptibility to phytoplasma infection, due to alterations in expression of genes involved in sugar metabolism and membrane transport. In the long run, the identification of plant resistance or susceptibility traits against phytoplasmas will allow a complete re-organization of chemical control strategies, with obvious opportunities of reducing insecticide burden. A unified mechanism for innate and learned visual landmark guidance in the insect central complex Insects can navigate efficiently in both novel and familiar environments, and this requires flexiblity in how they are guided by sensory cues. A prominent landmark, for example, can ellicit strong innate behaviours (attraction or menotaxis) but can also be used, after learning, as a specific directional cue as part of a navigation memory. However, the mechanisms that allow both pathways to co-exist, interact or override each other are largely unknown. Here we propose a model for the behavioural integration of innate and learned guidance based on the neuroanatomy of the central complex (CX), adapted to control landmark guided behaviours. We consider a reward signal provided either by an innate attraction to landmarks or a long-term visual memory in the mushroom bodies (MB) that modulates the formation of a local vector memory in the CX. Using an operant strategy for a simulated agent exploring a simple world containing a single visual cue, we show how the generated short-term memory can support both innate and learned steering behaviour. In addition, we show how this architecture is consistent with the observed effects of unilateral MB lesions in ants that cause a reversion to innate behaviour. We suggest the formation of a directional memory in the CX can be interpreted as transforming rewarding (positive or negative) sensory signals into a mapping of the environment that describes the geometrical attractiveness (or repulsion). We discuss how this scheme might represent an ideal way to combine multisensory information gathered during the exploration of an environment and support optimized cue integration.

1 Author summaryIn this paper, we modeled the neural pathway allowing insects to perform landmark guided behaviours using their internal compass. First, we observed the intrinsic property of the connectome, extracted from drosophila online database, between the internal compass neurons and the steering neurons to support an oriented behaviour towards a single landmark. Then, we proposed and evaluated an adaptation of the bees path integration neural circuit, to sustain flexible landmark guidance behaviours such as attraction or menotaxis. We showed the model ability to form a memory during the exploration of the local environment to support both innate or learned navigation behaviour using a single landmark in the environment. In addtion, we demonstrated the transformation of a simple goodness/badness signal, from innate or long-term memory pathways, into an oriented steering signal that could be applied to other sensory pathways. Furthermore, by reproducing lesion experiments in the mushroom bodies of wood ants we highlight the consistency of the model with biological observation. We then discuss the different emergent properties and the potential outcome that this local, and operant, memory supports. A burst of regulatory and protein innovation at the origin of placental mammals drove the emergence of placenta and underpins divergent early pregnancy strategies in modern mammals. Summary paragraphThe origin of live birth in mammals [~]148 million years ago was a dramatic shift in reproductive strategy, yet the molecular changes that established mammal viviparity are largely unknown. Although progesterone receptor signalling predates the origin of mammals andis highly conserved in, and critical for, successful mammal pregnancy, it alone cannot explain the origin and subsequent diversity of implantation strategies throughout the placental mammal radiation. MiRNAs are known to be flexible and dynamic regulators with a well established role in the pathophysiology of mammal placenta. We propose that a dynamic core microRNA (miRNA) network originated early in placental mammal evolution, responds to conserved mammal pregnancy cues (e.g. progesterone), and facilitates species-specific responses. Here we identify 13 miRNAs that arose at the origin of placental mammals and were subsequently retained in all descendent lineages. The expression of these 13 miRNAs in response to early pregnancy molecules is regulated in a species-specific manner in endometrial epithelia of species with extremes of implantation strategies. Furthermore, these 13 miRNAs preferentially target 84 proteins under positive selective pressure on ancestral eutherian. Discovery of this core "live-birth" toolkit and specifically adapted proteins helps explain the origin and evolution of the placenta in mammals. RNA degradation of blood reveals a time diagnostic window in a controlled environment Determining the time since deposition (TSD) of bloodstains would provide forensic scientists with critical information regarding the timeline of the events involving bloodshed. The physicochemical changes occurring to biomolecules as a bloodstain degrades can be used to approximate the TSD of bloodstains. Our study aims to quantify the timewise degradation trends and temperature dependence found in total RNA from bloodstains without the use of amplification, expanding the scope of the RNA TSD research which has previously targeted mRNA molecules. Whole bovine blood was stored in plastic microcentrifuge tubes at 21{degrees}C or 4{degrees}C and tested over different timepoints spanning one week. Nine RNA metrics were visually assessed and quantified using linear and mixed models; the RNA Integrity Number equivalent (RINe) and the DV200 demonstrated strong negative trends over time and statistical independence. The RINe model fit was high (R2 = 0.60), and while including the biological replicate as a random effect increased the fit for all RNA metrics, no significant differences were found between biological replicates stored at the same temperature for the RINe and DV200 metrics. Importantly, this suggests that these standardized metrics can likely be directly compared between scenarios and individuals, with DV200 having an inflection point at [~]28 hrs. This study provides a novel approach for blood TSD estimates, producing metrics that are not affected by inter-individual variation and improving our understanding of the rapid degradation occurring in bloodstains.

HIGHLIGHTSO_LIAmplification-free analysis of total RNA in degrading bloodstains.
C_LIO_LIShort-term RNA degradation assessment using high-resolution size measurements.
C_LIO_LITotal RNA quality and quantity metrics were assessed across a one-week.
C_LIO_LITotal RNA quality metrics demonstrated the strongest timewise trends.
C_LIO_LIBiological replicates produced similar results for RNA quality metrics.
C_LI Yersinia remodels epigenetic histone modifications in human macrophages Various pathogens systematically reprogram gene expression in macrophages, but the underlying mechanisms are largely unknown. We investigated whether the enteropathogen Yersinia enterocolitica alters chromatin states to reprogram gene expression in primary human macrophages. Genome-wide chromatin immunoprecipitation (ChIP) seq analyses showed that pathogen-associated molecular patterns (PAMPs) induced up- or down-regulation of histone modifications (HMod) at approximately 14500 loci in promoters and enhancers. Effectors of Y. enterocolitica reorganized about half of these dynamic HMod, with the effector YopP being responsible for about half of these modulatory activities. The reorganized HMod were associated with genes involved in immune response and metabolism. Remarkably, the altered HMod also associated with 61 % of all 534 known Rho GTPase pathway genes, revealing a new level in Rho GTPase regulation and a new aspect of bacterial pathogenicity. Changes in HMod were associated to varying degrees with corresponding gene expression, e. g. depending on chromatin localization and cooperation of the HMod. In summary, infection with Y. enterocolitica remodels HMod in human macrophages to modulate key gene expression programs of the innate immune response.

Author SummaryHuman pathogenic bacteria can affect epigenetic histone modifications to modulate gene expression in host cells. However, a systems biology analysis of this bacterial virulence mechanism in immune cells has not been performed. Here we analyzed genome-wide epigenetic histone modifications and associated gene expression changes in primary human macrophages infected with enteropathogenic Yersinia enterocolitica. We demonstrate that Yersinia virulence factors extensively modulate histone modifications and associated gene expression triggered by the pathogen-associated molecular patterns (PAMPs) of the bacteria. The epigenetically modulated genes are involved in several key pathways of the macrophage immune response, including the Rho GTPase pathway, revealing a novel level of Rho GTPase regulation by a bacterial pathogen. Overall, our findings provide an in-depth view of epigenetic and gene expression changes during host-pathogen interaction and might have further implications for understanding of the innate immune memory in macrophages. Dipeptides are minimalistic but sufficient for liquid-liquid phase separation Liquid-liquid phase separation (LLPS) of proteins mediates the assembly of biomolecular condensates involved in physiological and pathological processes. Identifying the minimalistic building blocks and the sequence determinant of protein phase separation is of urgent importance but remains challenging due to the enormous sequence space and difficulties of existing methodologies in characterizing the phase behavior of ultrashort peptides. Here we demonstrate computational tools to efficiently quantify the microscopic fluidity and density of liquid-condensates/solid-aggregates and the temperature-dependent phase diagram of peptides. Utilizing our approaches, we comprehensively predict the LLPS abilities of all 400 dipeptide combinations of coded amino acids based on 492 micro-second molecular dynamics simulations, and observe the occurrences of spontaneous LLPS. We identify 54 dipeptides that form solid-like aggregates and three categories of dipeptides with high LLPS propensity. Our predictions are validated by turbidity assays and differential interference contrast (DIC) microscopy on four representative dipeptides (WW, QW, GF, and VI). Phase coexistence diagrams are constructed to explore the temperature dependence of LLPS. Our results reveal that aromatic moieties are crucial for a dipeptide to undergo LLPS, and hydrophobic and polar components are indispensable. We demonstrate for the first time that dipeptides, minimal but complete, possess multivalent interactions sufficient for LLPS, suggesting that LLPS is a general property of peptides/proteins, independent of their sequence length. This study provides a computational and experimental approach to the prediction and characterization of the phase behavior of minimalistic peptides, and will be helpful for understanding the sequence-dependence and molecular mechanism of protein phase separation.

SignificanceProtein liquid-liquid phase separation (LLPS) is associated with human health and diseases. Identifying the minimalistic building blocks and sequence determinants of LLPS is of urgent importance but remains computationally challenging partially due to the lack of methodologies characterizing the liquid condensates. Herein we provide approaches to evaluate LLPS ability of dipeptides, and screen all 400 dipeptides by MD simulations combined with multi-bead-per-residue models which capture key interactions driving LLPS that are missing in one-bead-per-residue models. Three categories of LLPS dipeptides are identified and the experimentally-verified QW dipeptide is by far the smallest LLPS system. Our results suggest that dipeptides, minimal but complete, possess multivalent interactions sufficient for LLPS, and LLPS is a general property of peptides/proteins, independent of their length. Evidence for a hydrogen sulfide-sensing E3 ligase in yeast In yeast, control of sulfur amino acid metabolism relies upon Met4, a transcription factor which activates the expression of a network of enzymes responsible for the biosynthesis of cysteine and methionine. In times of sulfur abundance, the activity of Met4 is repressed via ubiquitination by the SCFMet30 E3 ubiquitin ligase, but the mechanism by which the F-box protein Met30 senses sulfur status to tune its E3 ligase activity remains unresolved. Herein, we show that Met30 responds to flux through the transsulfuration pathway to regulate the MET gene transcriptional program. In particular, Met30 is responsive to the biological gas hydrogen sulfide, which is sufficient to induce ubiquitination of Met4 in vivo. Additionally, we identify important cysteine residues in Met30s WD-40 repeat region that sense the availability of sulfur in the cell. Our findings reveal how SCFMet30 dynamically senses the flow of sulfur metabolites through the transsulfuration pathway to regulate synthesis of these special amino acids. Network Inference from Perturbation Time Course Data Networks underlie much of biology from subcellular to ecological scales. Yet, understanding what experimental data are needed and how to use them for unambiguously identifying the structure of even small networks remains a broad challenge. Here, we integrate a dynamic least squares framework into established modular response analysis (DL-MRA), that specifies sufficient experimental perturbation time course data to robustly infer arbitrary two and three node networks. DL-MRA considers important network properties that current methods often struggle to capture: (i) edge sign and directionality; (ii) cycles with feedback or feedforward loops including self-regulation; (iii) dynamic network behavior; (iv) edges external to the network; and (v) robust performance with experimental noise. We evaluate the performance of and the extent to which the approach applies to cell state transition networks, intracellular signaling networks, and gene regulatory networks. Although signaling networks are often an application of network reconstruction methods, the results suggest that only under quite restricted conditions can they be robustly inferred. For gene regulatory networks, the results suggest that incomplete knockdown is often more informative than full knockout perturbation, which may change experimental strategies for gene regulatory network reconstruction. Overall, the results give a rational basis to experimental data requirements for network reconstruction and can be applied to any such problem where perturbation time course experiments are possible. 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. Modular evolution of secretion systems and virulence plasmids in a bacterial species complex BackgroundMany bacterial taxa are species complexes and uncertainties regarding the organization of their genetic diversity challenge research efforts. We utilized Agrobacterium tumefaciens, a taxon known for its phytopathogenicity and applications in transformation, as a study system and devised strategies for investigating genome diversity and evolution of species complexes.

ResultsWe utilized 35 genome assemblies to achieve a comprehensive and balanced sampling of A. tumefaciens. Our confident inference of gene content and core-genome phylogeny supported a quantitative guideline for delineating 12 species and allowed for robust investigations of genes critical in fitness and ecology. For the type VI secretion system (T6SS) involved in interbacterial competition and thought to be conserved, we detected multiple losses and one horizontal gene transfer. For the tumor-inducing plasmids (pTi) and pTi-encoded type IV secretion system (T4SS) that are essential for agrobacterial phytopathogenicity, we uncovered novel diversity and hypothesized their involvement in shaping this species complex. Intriguingly, for both T6SS and T4SS, genes encoding structural components are highly conserved, whereas extensive diversity exists for genes encoding effectors and other proteins.

ConclusionsWe demonstrated that the combination of a phylogeny-guided sampling scheme and an emphasis on high-quality assemblies provides a cost-effective approach for robust analysis in evolutionary genomics. Our strategies for multi-level investigations at scales that range from whole-genomes to intragenic domains and phylogenetic depths of between- and within-species are applicable to other bacteria. Finally, modularity observed in the molecular evolution of genes and domains is useful for inferring functional constraints and informing experimental works. 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. ScRNA-seq Expression of IFI27 and APOC2 Identifies Four Alveolar Macrophage Superclusters in Healthy BALF Alveolar macrophages (AMs) reside on the luminal surface of the airways and alveoli, ensuring proper gas exchange by ingesting cellular debris and pathogens, and regulating inflammatory responses. Therefore, understanding the heterogeneity and diverse roles played by AMs and recruited monocytes is critical for treating airway diseases. We performed single-cell RNA sequencing on 113,213 bronchoalveolar lavage cells from four healthy and three uninflamed cystic fibrosis subjects and identified four AM superclusters (families) based on the expression of IFI27 and APOC2 genes. These 4 AM families have at least eight distinct functional members (subclusters) named after their differentially expressed gene(s): IGF1, CCL18, CXCL5, Cholesterol, Chemokine, Metallothionein, Interferon and small-cluster AMs. Interestingly, the Chemokine cluster further divides with each subcluster selectively expressing a unique combination of chemokines. One of the most striking observations, besides the heterogeneity, is the conservation of AM family members in relatively equal ratio across all AM superclusters and individuals. Finally, transcriptional data and TotalSeq technology were used to investigate cell surface markers that distinguish resident AMs from recruited monocytes. Overall, functional specializations of the AM subclusters suggest that there are AM niches with defined programming states, highlighting a clear division of labor.