Nucleic Transformer: Deep Learning on Nucleic Acids with Self-attention and Convolutions AO_SCPLOWBSTRACTC_SCPLOWMuch work has been done to apply machine learning and deep learning to genomics tasks, but these applications usually require extensive domain knowledge and the resulting models provide very limited interpretability. Here we present the Nucleic Transformer, a conceptually simple but effective and interpretable model architecture that excels in a variety of DNA/RNA tasks. The Nucleic Transformer processes nucleic acid sequences with self-attention and convolutions, two deep learning techniques that have proved dominant in the fields of computer vision and natural language processing. We demonstrate that the Nucleic Transformer can be trained in both supervised and unsupervised fashion without much domain knowledge to achieve high performance with limited amounts of data in Escherichia coli promoter classification, viral genome identification, and degradation properties of COVID-19 mRNA vaccine candidates. Additionally, we showcase extraction of promoter motifs from learned attention and how direct visualization of self-attention maps assists informed decision making using deep learning models.
Group-walk, a rigorous approach to group-wise false discovery rate analysis by target-decoy competition Target-decoy competition (TDC) is a commonly used method for false discovery rate (FDR) control in the analysis of tandem mass spectrometry data. This type of competitionbased FDR control has recently gained significant popularity in other fields after Barber and Candes laid its theoretical foundation in a more general setting that included the feature selection problem. In both cases, the competition is based on a head-to-head comparison between an (observed) target score and a corresponding decoy (knockoff) score. However, the effectiveness of TDC depends on whether the data is homogeneous, which is often not the case: in many settings, the data consists of groups with different score profiles or different proportions of true nulls. In such cases, applying TDC while ignoring the group structure often yields imbalanced lists of discoveries, where some groups might include relatively many false discoveries and other groups include relatively very few. On the other hand, as we show, the alternative approach of applying TDC separately to each group does not rigorously control the FDR.

We developed Group-walk, a procedure that controls the FDR in the target-decoy / knockoff setting while taking into account a given group structure. Group-walk is derived from the recently developed AdaPT -- a general framework for controlling the FDR with sideinformation. We show using simulated and real datasets that when the data naturally divides into groups with different characteristics Group-walk can deliver consistent power gains that in some cases are substantial. These groupings include the precursor charge state (4% more discovered peptides at 1% FDR threshold), the peptide length (3.6% increase) and the mass difference due to modifications (26% increase).

Group-walk is available at https://cran.r-project.org/web/packages/groupwalk/index.html
Genome-wide association study of musical beat synchronization demonstrates high polygenicity Moving in synchrony to the beat is a fundamental component of musicality. Here, we conducted a genome-wide association study (GWAS) to identify common genetic variants associated with beat synchronization in 606,825 individuals. Beat synchronization exhibited a highly polygenic architecture, with sixty-nine loci reaching genome-wide significance (p<5x10-8) and SNP-based heritability (on the liability scale) of 13%-16%. Heritability was enriched for genes expressed in brain tissues, and for fetal and adult brain-specific gene regulatory elements, underscoring the role of central nervous system-expressed genes linked to the genetic basis of the trait. We performed validations of the self-report phenotype (through internet-based experiments) and of the GWAS (polygenic scores for beat synchronization were associated with patients algorithmically classified as musicians in medical records of a separate biobank). Genetic correlations with breathing function, motor function, processing speed, and chronotype suggest shared genetic architecture with beat synchronization and provide avenues for new phenotypic and genetic explorations. Simulation of foraging behavior using a decision-making agent with Bayesian and inverse Bayesian inference: Temporal correlations and power laws in displacement patterns It has been stated that in human migratory behavior, the step length series may have temporal correlation and that there is some relationship between this time dependency and the fact that the frequency distribution of step length follows a power-law distribution. Furthermore, the frequency of occurrence of the step length in some large marine organisms has been found to switch between power-law and exponential distributions, depending on the difficulty of prey acquisition. However, to date it has not been clarified how the aforementioned three phenomena arise: the positive correlation created in the step length series, the relation between the positive correlation of the step length series and the form of an individuals step length distribution, and the switching between power-law and exponential distributions depending on the abundance of prey. This study simulated foraging behavior using the Bayesian decision-making agent simultaneously performing both knowledge learning and knowledge-based inference to analyze how the aforementioned three phenomena arise. In the agent with learning and inference, past experiences were stored as hypotheses (knowledge) and they were used in current foraging behavior; at the same time, the hypothesis continued to be updated based on new experiences. The simulation results show that the agent with both learning and inference has a mechanism that simultaneously causes all the phenomena. Repeated truncation of a modular antimicrobial peptide gene for neural context Antimicrobial peptides (AMPs) are host-encoded antibiotics that combat invading pathogens. These genes commonly encode multiple products as post-translationally cleaved polypeptides. Recent studies have highlighted roles for AMPs in neurological contexts suggesting functions for these defence molecules beyond infection. During our immune study characterizing the antimicrobial peptide gene Baramicin, we recovered multiple Baramicin paralogs in Drosophila melanogaster and other species, united by their N-terminal IM24 domain. Not all paralogs were immune-induced. Here, through careful dissection of the Baramicin familys evolutionary history, we find that these non-immune paralogs result from repeated events of duplication and subsequent truncation of the coding sequence from an immune-inducible ancestor. These truncations leave only the IM24 domain as the prominent gene product. Surprisingly, using mutation and targeted gene silencing we demonstrate that two such genes are adapted for function in neural contexts in D. melanogaster. We also show enrichment in the head for independent Baramicin genes in other species. The Baramicin evolutionary history reveals that the IM24 Baramicin domain is not strictly useful in an immune context. We thus provide a case study for how an AMP-encoding gene might play dual roles in both immune and non-immune processes via its multiple peptide products. We reflect on these findings to highlight a blind spot in the way researchers approach AMP research in in vivo contexts.

Significance statementAntimicrobial peptides are immune proteins recently implicated in neurological roles. To date little attention has been paid to the contributions of different gene products in this function. Here we show that an antimicrobial peptide gene encodes multiple products with either immune-specific or neurological roles. Noradrenergic control of non-REM sleep substates One promising approach towards understanding what makes sleep vulnerable in disease is to look at how wake-promoting mechanisms affect healthy sleep. Wake-promoting neuronal activity is inhibited during non-REM sleep (NREMS). However, many mammalian species, including humans, show recurrent moments of sleep fragility during which sensory reactivity is elevated. Wake-related neuronal activity could thus remain active in NREMS, but its roles in dynamic variations of sensory reactivity remain unknown. Here, we demonstrate that mouse NREMS is a brain state with recurrent fluctuations of the wake-promoting neurotransmitter noradrenaline on the [~]50-seconds time-scale. These fluctuations occurred around mean noradrenaline levels greater than the ones of quiet wakefulness, while they declined steeply in REMS. They coincided with a clustering of sleep spindle rhythms in the forebrain and with heart rate variations. We addressed the origins of these fluctuations by using closed-loop optogenetic locus coeruleus (LC) activation or inhibition timed to moments of low and high spindle activity during NREMS. We could suppress, lock or entrain sleep spindle clustering or heart rate variations, demonstrating that both fore- and hindbrain-projecting LC neurons show synchronized infraslow activity variations in natural NREMS. Noradrenergic modulation of thalamic but not cortical circuits was required for sleep spindle clustering and involved noradrenaline release into primary sensory and reticular thalamic nuclei that activated both 1- and {beta}-adrenergic receptors to cause slowly decaying membrane depolarizations. Noradrenergic signaling by LC, primarily known for attention promotion in wakefulness, renders mammalian NREMS more  wake-like on the close-to-minute-time scale through sustaining thalamocortical and autonomic sensory arousability. Naturally occurring combinations of receptors from single cell transcriptomics in endothelial cells VEGF inhibitor drugs have been successful, especially in ophthalmology, but not all patients respond to them. Combinations of drugs are likely to be needed for more effective therapies of angiogenesis-related diseases. In this paper we describe naturally occurring combinations of receptors in endothelial cells that might help to understand how cells communicate and to identify targets for drug combinations. We also develop and share a new software tool called DECNEO to identify them.

Single-cell gene expression data are used to identify a set of co-expressed endothelial cell receptors, conserved among species (mice and humans) and enriched, within a network, of connections to up-regulated genes. This set includes several receptors previously shown to play a role in angiogenesis. Multiple statistical tests from large datasets, including an independent validation set, support the reproducibility, evolutionary conservation and role in angiogenesis of these naturally occurring combinations of receptors. We also show tissue-specific combinations and, in the case of choroid endothelial cells, consistency with both well-established and recent experimental findings, presented in a separate paper.

The results and methods presented here advance the understanding of signaling to endothelial cells. The methods are generally applicable to the decoding of intercellular combinations of signals. Nascent alt-protein chemoproteomics reveals a repressor of ribosome biogenesis Many unannotated microproteins and alternative proteins (alt-proteins) have recently been found to be co-encoded with canonical proteins, but few of their functions are known. Motivated by the hypothesis that alt-proteins undergoing active or stress-induced synthesis could play important cellular roles, here, we developed a chemoproteomic pipeline to identify nascent alt-proteins in human cells. We identified 22 actively translated unannotated alt-proteins, one of which is upregulated after DNA damage stress. We further defined MINAS-60 (MIcroprotein that Negatively regulates ASsembly of the pre-60S ribosomal subunit), a nucleolar localized alt-protein co-encoded with human RBM10.Depletion of MINAS-60 increases the amount of the mature 60S ribosomal subunit, consequently upregulating global protein synthesis and cell proliferation by repressing late-stage pre-60S assembly and export of the 60S ribosome subunit to the cytoplasm. Together, these results implicate MINAS-60 as a repressor of ribosome biogenesis, and demonstrate that chemoproteomics can enable generation of functional hypotheses for uncharacterized alt-proteins. Priority conservation areas and a global population estimate for the Critically Endangered Philippine Eagle derived from modelled range metrics using remote sensing habitat characteristics Many range-restricted taxa are currently experiencing population declines yet lack fundamental information regarding distribution and population size. Establishing baseline estimates for both these key biological parameters is however critical for directing conservation planning for at-risk range-restricted species. The International Union for the Conservation of Nature (IUCN) Red List uses three range metrics that define species distributions and inform extinction risk assessments: extent of occurrence (EOO), area of occupancy (AOO) and area of habitat (AOH). However, calculating all three metrics using standard IUCN approaches relies on a geographically representative sample of locations, which for rare species is often spatially biased. Here, we apply model-based interpolation using Species Distribution Models (SDMs), correlating occurrences with remote-sensing covariates, to calculate IUCN range metrics, protected area coverage and a global population estimate for the Critically Endangered Philippine Eagle (Pithecophaga jefferyi). Our final range wide continuous SDM had high predictive accuracy (Continuous Boyce Index = 0.927) and when converted to a binary model estimated an AOH = 23,185 km2, a maximum EOO = 605,759 km2, a minimum EOO = 272,272 km2, with an AOO = 53,867 km2. Based on inferred habitat from the AOH metric, we estimate a global population of 318 breeding pairs (range: 258-362 pairs), or 636 mature individuals, across the Philippine Eagle global range. Protected areas covered 34 % of AOH, 15 % less than the target representation, with the continuous model identifying key habitat as priority conservation areas. We demonstrate that even when occurrences are geographically biased, robust habitat models can be built that enable quantification of baseline IUCN range metrics, protected area coverage, and a population size estimate. In the absence of adequate location data for many rare and threatened taxa, our method is a promising spatial modelling tool with widespread applications, in particular for island endemics facing high extinction risk. A novel stem cell type at the basal side of the subventricular zone maintains adult neurogenesis According to the current consensus, neural stem cells (NSCs) apically contacting the lateral ventricle generate differentiated progenitors by rare asymmetric divisions or by relocating to the basal side of the ventricular-subventricular zone V-SVZ. Both processes will then ultimately lead to the generation of adult-born olfactory bulb (OB) interneurons. In contrast to this view, we here found that adult-born OB interneurons largely derive from an additional NSC type resident in the basal V-SVZ. Despite being both capable of self-renewal and long-term quiescence, apical and basal NSCs differ in Nestin expression, primary cilia extension and frequency of cell division. The expression of Notch-related genes also differed between the two NSC groups and Notch-activation was greatest in apical NSCs. Apical downregulation of Notch-effector Hes1 decreased Notch activation while increasing proliferation across the niche and neurogenesis from apical NSCs. Underscoring their different roles in neurogenesis, lactation-dependent increase in neurogenesis was paralleled by extra activation of basal but not apical NSCs. Thus, basal NSCs support OB neurogenesis whereas apical NSCs impart Notch-mediated lateral inhibition across the V-SVZ. 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. Lower limb in sprinters with larger relative mass but not larger normalized moment of inertia Inhomogeneous muscle hypertrophy in response to daily motor executions is a morphological adaptation in humans. Inhomogeneous morphological adaptation would change the mass distribution, which in turn may affect the mechanical difficulty in moving from an inertia perspective; however, the morphological adaptation from the inertia perspective has not been examined. Here we show no corresponding differences in the normalized mass and normalized moment of inertia between the sprinters and the untrained non-sprinters. We analyzed the fat- and water-separated magnetic resonance images from the lower limbs of 11 male sprinters (100 m best time of 10.44-10.83 s) and 12 untrained non-sprinters. The lower-limb relative mass was significantly larger in sprinters (18.7{+/-}0.7% body mass) than in non-sprinters (17.6{+/-}0.6% body mass), while the normalized moment of inertia of the lower limb around the hip in the anatomical position was similar between them (0.044{+/-}0.002 vs. 0.042{+/-}0.002 [a. u.]). The thigh relative mass in sprinters (12.9{+/-}0.4% body mass) was significantly larger than that in non-sprinters (11.9{+/-}0.4% body mass), whereas the shank and foot relative masses were similar. We revealed that the adaptation to sprinting in humans does not relatively increase the mechanical difficulty to swing the lower limb in terms of inertia, even though the lower-limb mass increases, reflecting muscular hypertrophy. This is suggested to be one of the significant advantages of human morphology for adaptability to explosive motor tasks.

NEW & NOTEWORTHYWe compared the lower-limb inertia properties between male sprinters and untrained non-sprinters by using the fat- and water-separated magnetic resonance images. Sprinters had lower limbs with a larger relative mass but not a larger normalized moment of inertia than non-sprinters. We suggest one of the practical advantages of human morphological characteristics for their adaptation to explosive motor tasks, which provides a novel perspective to the literature for understanding human morphological adaptability. Recombination marks the evolutionary dynamics of a recently endogenized retrovirus All vertebrate genomes have been colonized by retroviruses along their evolutionary trajectory. While endogenous retroviruses (ERVs) can contribute important physiological functions to contemporary hosts, such benefits are attributed to long-term co-evolution of ERV and host because germline infections are rare and expansion is slow, because the host effectively silences them. The genomes of several outbred species including mule deer (Odocoileus hemionus) are currently being colonized by ERVs, which provides an opportunity to study ERV dynamics at a time when few are fixed. Because we have locus-specific data on the distribution of cervid endogenous retrovirus (CrERV) in populations of mule deer, in this study we determine the molecular evolutionary processes acting on CrERV at each locus in the context of phylogenetic origin, genome location, and population prevalence. A mule deer genome was de novo assembled from short and long insert mate pair reads and CrERV sequence generated at each locus. CrERV composition and diversity have recently measurably increased by horizontal acquisition of a new retrovirus lineage. This new lineage has further expanded CrERV burden and CrERV genomic diversity by activating and recombining with existing CrERV. Resulting inter-lineage recombinants endogenized and subsequently retrotransposed. CrERV loci are significantly closer to genes than expected if integration were random and gene proximity might explain the recent expansion by retrotransposition of one recombinant CrERV lineage. Thus, in mule deer, retroviral colonization is a dynamic period in the molecular evolution of CrERV that also provides a burst of genomic diversity to the host population. Quercetin handles cellular oxidant / antioxidant systems and mitigates immunosenescence hallmarks in human PBMCs: an in vitro study Immunosenescence, oxidative stress, and low vaccine efficacy are important symptoms of aging. The goal of our study was to see if quercetin had anti-aging and stimulating effects on PBMC immune cells in vitro. In the presence of concanavalin a, PBMCs were isolated from healthy elderly and young people and cultured in a complete RPMI-1640 medium supplemented with quercetin. Cell proliferation was assessed using the MTT colorimetric assay after a 48-hour incubation period. Spectrophophotometric assays were used to assess oxidative biomarkers (PC, MDA, and GSH). The ELISA method was used to determine the amount of IL-2 released. A Griess reagent was used to investigate iNOS activity. When compared to young control cells, aged PBMCs had lower proliferation potency, lower IL-2 and NO release, and higher MDA and PC levels. Importantly, quercetin-treated aged PBMCs have a high proliferative response comparable to young cells, restored iNOS activity, and increased levels of GSH antioxidant defences. In comparison to untreated aged PBMCs, treated PBMCs have lower lipo-oxidative damage but higher PC levels. Quercetin may be used as a promising dietary vaccinal adjuvant in the elderly, it has significant effects in reducing immunosenescence hallmarks, as well as mitigating the lipo-oxidative stress in PBMCs cells.

Practical abstractQuercetin is a common dietary polyphenol found in fruits and vegetables that is well known for its antioxidant activity in radical scavenging and anticancer properties characterized by immune system stimulation. In vitro, quercetin boosts immune responses in the elderly and reduces immunosenescence symptoms, according to our findings. As a result, this molecule could be a safe and promising vaccine adjuvant for boosting immunity and reducing aging complications in immune cells as well as reducing oxidative stress in PBMC cells. A tool for assessing alignment of biomedical data repositories with open, FAIR, citation and trustworthy principles Increasing attention is being paid to the operation of biomedical data repositories in light of efforts to improve how scientific data is handled and made available for the long term. Multiple groups have produced recommendations for functions that biomedical repositories should support, with many using requirements of the FAIR data principles as guidelines. However, FAIR is but one set of principles that has arisen out of the open science community. They are joined by principles governing open science, data citation and trustworthiness, all of which are important aspects for biomedical data repositories to support. Together, these define a framework for data repositories that we call OFCT: Open, FAIR, Citable and Trustworthy. Here we developed an instrument using the open source PolicyModels toolkit that attempts to operationalize key aspects of OFCT principles and piloted the instrument by evaluating eight biomedical community repositories listed by the NIDDK Information Network (dkNET.org). Repositories included both specialist repositories that focused on a particular data type or domain, in this case diabetes and metabolomics, and generalist repositories that accept all data types and domains. The goal of this work was both to obtain a sense of how much the design of current biomedical data repositories align with these principles and to augment the dkNET listing with additional information that may be important to investigators trying to choose a repository, e.g., does the repository fully support data citation? The evaluation was performed from March to November 2020 through inspection of documentation and interaction with the sites by the authors. Overall, although there was little explicit acknowledgement of any of the OFCT principles in our sample, the majority of repositories provided at least some support for their tenets. Neutralization potential of Covishield vaccinated individuals sera against B.1.617.1 Covishield comprises the larger proportion in the vaccination program in India. Hence, it is of utmost importance to understand neutralizing capability of vaccine against the B.1.617.1 variant which is considered to responsible for surge of the cases in India. The neutralizing-antibody (NAb) titer against B.1.167.1 and prototype B.1 variant (D614G) was determined of the vaccine sera (4 weeks after second dose) of COVID-19 naive subjects (n=43) and COVID-19 recovered subjects (n=18). The results demonstrated that sera of COVID-19 recovered subjects (n=18) who received two doses of Covishield have higher NAb response compared to the COVID-19 naive with a significant difference (p<0.0001) in NAb titer against B.1 and B.1.617.1 In-spite of reduction in the neutralizing titer against B.1.617.1 variant; Covishield vaccine-induced antibodies are likely to be protective to limit the severity and mortality of the disease in the vaccinated individuals. FGF8-mediated signaling regulates tooth developmental pace during odontogenesis The developing human and mouse teeth constitute an ideal model system to study the regulatory mechanism underlying organ growth control due to the fact that their teeth share highly conserved and well-characterized developmental processes and their developmental tempo vary notably. In the current study, we manipulated heterogenous recombination between human and mouse dental tissues and demonstrate that the dental mesenchyme dominates the tooth developmental tempo and FGF8 could be a critical player during this developmental process. Forced activation of FGF8 signaling in the dental mesenchyme of mice promoted cell proliferation, prevented cell apoptosis via p38 and perhaps PI3K-Akt intracellular signaling, and impelled the transition of the cell cycle from G1- to S-phase in the tooth germ, resulting in the slowdown of the tooth developmental pace. Our results provide compelling evidence that extrinsic signals can profoundly affect tooth developmental tempo and the dental mesenchymal FGF8 could be a pivotal factor in controlling developmental pace in a non-cell-autonomous manner during mammalian odontogenesis. Multiplex knockout of trichome-regulating MYB duplicates in hybrid poplar using a single gRNA As the focus for CRISPR edited plants moves from proof-of-concept to real world applications, precise gene manipulation will increasingly require concurrent multiplex editing for polygenic traits. A common approach for editing across multiple sites is to design one gRNA per target; however, this complicates construct assembly and increases the possibility of off-target mutations. In this study, we utilized one gRNA to target MYB186, a known positive trichome regulator, as well as its paralogs MYB138 and MYB38 at a consensus site for mutagenesis in Populus tremula x P. alba INRA 717-1B4. Unexpected duplications of MYB186 and MYB138 resulted in a total of eight alleles for the three targeted genes in the hybrid poplar. Deep sequencing and PCR analyses confirmed editing across all eight targets in nearly all of the resultant glabrous mutants, ranging from small indels to large genomic dropouts, with no off-target activity detected at four potential sites. This highlights the effectiveness of a single gRNA targeting conserved exonic regions for multiplex editing. Additionally, cuticular wax and whole leaf analyses showed a complete absence of triterpenes in the trichomeless mutants, hinting at a previously undescribed role for the non-glandular trichomes of poplar.

ONE SENTENCE SUMMARYTargeting conserved sequences with a single gRNA allowed efficient mutagenesis of a multigene family and the recovery of trichomeless and triterpene-free poplar mutants. Re-distribution of antibiotic within a community is behind the unpredictable change in efficacy observed in classic laboratory assays Microbes associate in nature forming complex communities, but they are often studied in purified form. Here I show that differences in antibiotic sensitivity between co-existing bacterial species can modify the diffusion of antibiotic molecules, and alter the drug efficacy measured in conventional in vitro assays. I developed a model, validated experimentally using an artificial two-species community, showing antibiotic molecules re-distribute evenly across all species if they have similar drug sensitivity. The result is lower drug efficacy with respect to conventional in vitro assays. However, when species have different sensitivities, the re-distribution of antibiotic is uneven and most sensitive species accumulate more drug. The result is increased drug efficacy against these species. Despite this discrepancy with conventional laboratory assays, my data suggests the changes in drug efficacy are predictable and provides insight into how communities can affect the diffusion of other molecules in the environment. Chemotherapy coupled to macrophage inhibition leads to enhanced T and B cell infiltration and durable triple negative breast cancer regression Immunosuppressive elements within the tumor microenvironment such as Tumor Associated Macrophages (TAMs) can present a barrier to successful anti-tumor responses by cytolytic T cells. We employed preclinical syngeneic p53 null mouse models of triple negative breast cancer (TNBC) to develop a treatment regimen that harnessed the immunostimulatory effects of low-dose cyclophosphamide coupled with the pharmacologic inhibition of TAMs using either a small molecule CSF1R inhibitor or an anti-CSF1R antibody. This therapeutic combination was used to successfully treat several highly aggressive TNBC murine mammary tumors and lung metastasis. Using this regimen and single cell RNA sequencing we characterized tumor infiltrating lymphocytes (TILs) including helper T cells and antigen-presenting B cells that were highly enriched in good responders to combination therapy. Using high dimensional imaging techniques, we identified the close spatial localization of B220+ CD86+ activated B cells and CD4+ T cells in tertiary lymphoid structures that were present up to 6 weeks post-treatment in one model that also exhibited long-term tumor regression post-treatment. We also characterized the transcriptional and metabolic heterogeneity of TAMS in these two closely related claudin-low/mesenchymal subtype tumor models with differential treatment responses. A murine TAM signature derived from the T12 model is highly expressed and conserved in human claudin-low breast cancers, and high expression of the T12 signature correlated with reduced overall survival. This T12 tumor TAM signature may help identify human claudin-low breast cancer patients that will benefit from the combination of cyclophosphamide and anti-CSF1R therapy. These studies illustrate the complexity of the tumor immune microenvironment and highlight different immune responses that result from rationale combinations of immunotherapy.

SignificanceA treatment regimen that harnessed the immunostimulatory effects of cyclophosphamide coupled with the inhibition of CSF1R was used to successfully treat several highly aggressive claudin-low TNBC murine mammary tumors and lung metastasis. Spatial organization of chromosomes leads to heterogeneous chromatin motion and drives the liquid- or gel-like dynamical behavior of chromatin Chromosome organization and dynamics are involved in regulating many fundamental processes such as gene transcription and DNA repair. Experiments unveiled that chromatin motion is highly heterogeneous inside cell nuclei, ranging from a liquid-like, mobile state to a gel-like, rigid regime. Using polymer modeling, we investigate how these different physical states and dynamical heterogeneities may emerge from the same structural mechanisms. We found that the formation of topologically-associating domains (TADs) is a key driver of chromatin motion heterogeneity. In particular, we demonstrated that the local degree of compaction of the TAD regulates the transition from a weakly compact, fluid state of chromatin to a more compact, gel state exhibiting anomalous diffusion and coherent motion. Our work provides a comprehensive study of chromosome dynamics and a unified view of chromatin motion enabling to interpret the wide variety of dynamical behaviors observed experimentally across different biological conditions, suggesting that the  liquid or  solid behaviour of chromatin are in fact two sides of the same coin. Conserved organ-specific microbial assemblages in different populations of a terrestrial crab Brachyuran crabs originated in the oceans and evolved specific morphological and physiological adaptations to live in freshwater, intertidal and even terrestrial habitats but the role of a selection mechanism involving symbiotic microorganisms long these colonization processes are not known. In this work we investigated the associated microbiota of three populations of a terrestrial brachyuran crab, Chiromantes haematocheir, to find evidence of a conserved crab-specific microbiome unrelated to the population of origin and dissimilar from environmental microbial assemblages. Bacterial 16S rRNA gene and fungal ITS sequences were obtained from selected crab organs and environmental matrices to profile microbial communities. In spite of the presence of truly marine larval stages and no gregarious behaviour, favouring microbiota exchanges, we found common, organ-specific microbiota, associated to the gut and the gills of the crabs (with more than 15% of the genera detected specifically enriched only in one organ). Our results suggest an early establishment of a new common, stable microbiota in the transition from water to land. Lipid membranes modulate the activity of RNA through sequence-dependent interactions RNA is a ubiquitous biomolecule that can serve as both catalyst and information carrier. Understanding how RNA bioactivity is controlled is crucial for elucidating its physiological roles and potential applications in synthetic biology. Here we show that lipid membranes can act as RNA organization platforms, introducing a novel mechanism for ribo-regulation. The activity of R3C ribozyme can be modified by the presence of lipid membranes, with direct RNA-lipid interactions dependent on RNA nucleotide content, base pairing and length. In particular, the presence of guanine in short RNAs is crucial for RNA-lipid interactions, and G-quadruplex formation further promotes lipid binding. Lastly, by artificially modifying the R3C substrate sequence to enhance membrane binding we generated a lipid-sensitive ribozyme reaction with riboswitch-like behavior. These findings introduce RNA-lipid interactions as a tool for developing synthetic riboswitches and novel RNA-based lipid biosensors, and bear significant implications for RNA World scenarios for the origin of life. 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. Engineered bacteria detect tumor DNA Advances in bacterial engineering have catalysed the development of living cell diagnostics and therapeutics1-3, including microbes that respond to diseases such as gut inflammation4, intestinal bleeding5, pathogens6 and hypoxic tumors7. Bacteria can easily access the entire gastrointestinal tract via oral administration8, and they can produce outputs that can be noninvasively measured in stool4 or urine7. Cellular memory, such as bistable switches4,9,10 or genomic rearrangement11, has been used to allow bacteria to store information over time. However, living biosensors have not yet been engineered to detect specific DNA sequences or mutations from outside the cell. Here, we engineer naturally competent Acinetobacter baylyi to detect donor DNA from the genomes of colorectal cancer (CRC) cells and organoids. We characterize the functionality of the biosensors in vitro with co-culture assays and then validate in vivo with sensor bacteria delivered orally or rectally into mice injected with orthotopic donor CRC organoids. We observe horizontal gene transfer from the tumor to the sensor bacteria in vivo, allowing their detection in stool. The sensor bacteria achieved 100% discrimination between mice with and without CRC using both delivery methods. Our findings establish a framework for biosensing applications that require the detection of mutations or organisms within environments that are difficult to sample. In addition, the platform can be readily expanded to include in situ production and delivery of therapeutic payloads at the detection site. Non-Newtonian Blood Rheology Impacts Left Atrial Stasis in Patient-Specific Simulations The lack of mechanically effective contraction of the left atrium (LA) during atrial fibrillation (AF) disturbs blood flow, increasing the risk of thrombosis and ischemic stroke. Thrombosis is most likely in the left atrial appendage (LAA), a small narrow sac where blood is prone to stagnate. Slow flow promotes the formation of erythrocyte aggregates in the LAA, also known as rouleaux, causing viscosity gradients that are usually disregarded in patient-specific simulations. To evaluate these non-Newtonian effects, we built atrial models derived from 4D computed tomography scans of patients and carried out computational fluid dynamics simulations using the Carreau-Yasuda constitutive relation. We examined six patients, three of whom had AF and LAA thrombosis or a history of transient ischemic attacks (TIAs). We modeled the effects of hematocrit and rouleaux formation kinetics by varying the parameterization of the Carreau-Yasuda relation and modulating non-Newtonian viscosity changes based on residence time. Comparing non-Newtonian and Newtonian simulations indicates that slow flow in the LAA increases blood viscosity, altering secondary swirling flows and intensifying blood stasis. While some of these effects can be subtle when examined using instantaneous metrics like shear rate or kinetic energy, they are manifested in the blood residence time, which accumulates over multiple heartbeats. Our data also reveal that LAA blood stasis worsens when hematocrit increases, offering a potential new mechanism for the clinically reported correlation between hematocrit and stroke incidence. In summary, we submit that hematocrit-dependent non-Newtonian blood rheology should be considered in calculating patient-specific blood stasis indices by computational fluid dynamics.