Higher-order epistasis and phenotypic prediction Contemporary high-throughput mutagenesis experiments are providing an increasingly detailed view of the complex patterns of genetic interaction that occur between multiple mutations within a single protein or regulatory element. By simultaneously measuring the effects of thousands of combinations of mutations, these experiments have revealed that the genotype-phenotype relationship typically reflects genetic interactions not only between pairs of sites, but also higher-order interactions among larger numbers of sites. However, modeling and understanding these higher-order interactions remains challenging. Here, we present a method for reconstructing sequence-to-function mappings from partially observed data that can accommodate all orders of genetic interaction. The main idea is to make predictions for unobserved genotypes that match the type and extent of epistasis found in the observed data. This information on the type and extent of epistasis can be extracted by considering how phenotypic correlations change as a function of mutational distance, which is equivalent to estimating the fraction of phenotypic variance due to each order of genetic interaction (additive, pairwise, three-way, etc.). Using these estimated variance components, we then define an empirical Bayes prior that in expectation matches the observed pattern of epistasis, and reconstruct the genotype-phenotype mapping by conducting Gaussian process regression under this prior. To demonstrate the power of this approach, we present an application to the antibody-binding domain GB1 and also provide a detailed exploration of a dataset consisting of high-throughput measurements for the splicing efficiency of human pre-mRNA 5' splice sites, for which we also validate our model predictions via additional low-throughput experiments.
Microbial genomic trait evolution is dominated by frequent and rare pulsed evolution On the macroevolutionary timescale, does trait evolution proceed gradually or by rapid bursts (pulses) separated by prolonged periods of stasis or slow evolution? Although studies have shown pulsed evolution is prevalent in animals, our knowledge about the tempo and mode of evolution across the tree of life is very limited. This long-standing debate calls for a test in bacteria and archaea, the most ancient and diverse forms of life with unique population genetic properties (asexual reproduction, large population sizes, short generation times, high dispersal rates and extensive lateral gene transfers). Using a likelihood-based framework, we analyzed evolutionary patterns of four microbial genomic traits (genome size, genome GC%, 16S rRNA GC% and the nitrogen use in proteins) on a broad macroevolutionary timescale. Our model fitting of phylogenetic comparative data shows that pulsed evolution is not only present, but also prevalent and predominant in microbial genomic trait evolution. Interestingly, for the first time, we detected two distinct types of pulsed evolution (small frequent and large rare jumps) that are predicted by the punctuated equilibrium and quantum evolution theories. Our findings suggest that major bacterial lineages could have originated in quick bursts and pulsed evolution is a common theme across the tree of life despite the drastically different population genetic properties of bacteria, archaea and eukaryotes.
Gradient-mixing LEGO robots for purifying DNA origami nanostructures of multiple components by rate-zonal centrifugation DNA origami purification is critical in emerging applications of functionalized DNA nanostructures from basic fundamental biophysics, nanorobots to therapeutics. Advances in DNA origami purification have led to the establishment of rate-zonal centrifugation (RZC) as a scalable, high-yield, and contamination-free approach to purifying DNA origami nanostructures. In RZC purification, a linear density gradient is created using viscous agents, such as glycerol and sucrose, to separate molecules based on their mass and shape during high-rpm centrifugation. However, current methods for creating density gradients are typically time-consuming because of their reliance on slow passive diffusion. Here, we built a LEGO gradient mixer to rapidly create a quasi-continuous density gradient with minimal layering of concentrations using simple rotational motion. We found that rotating two layers of different concentrations at an angle can reduce the diffusion time from a few hours to mere minutes. The instrument needed to perform the movement can be constructed from low-cost components, such as Arduino and LEGO Mindstorms pieces, and has comparable efficacy to commercial gradient mixers currently available. Our results demonstrate that the creation of a linear density gradient can be achieved with minimal labor, time, and cost with this machine. With the recent advances in DNA origami production, we anticipate our findings to further improve the viability of scaling up DNA origami purification in grams quantities. Our simple process enables automated large-scale purification of functionalized DNA origami more feasible in resource-constrained settings. Crowdsourced identification of multi-target kinase inhibitors for RET- and TAU- based disease: the Multi-Targeting Drug DREAM Challenge A continuing challenge in modern medicine is the identification of safer and more efficacious drugs. Precision therapeutics, which have one molecular target, have been long promised to be safer and more effective than traditional therapies. This approach has proven to be challenging for multiple reasons including lack of efficacy, rapidly acquired drug resistance, and narrow patient eligibility criteria. An alternative approach is the development of drugs that address the overall disease network by targeting multiple biological targets ( polypharmacology). Rational development of these molecules will require improved methods for predicting single chemical structures that target multiple drug targets. To address this need, we developed the Multi-Targeting Drug DREAM Challenge, in which we challenged participants to predict single chemical entities that target pro-targets but avoid anti-targets for two unrelated diseases: RET-based tumors and a common form of inherited Tauopathy. Here, we report the results of this DREAM Challenge and the development of two neural network-based machine learning approaches that were applied to the challenge of rational polypharmacology. Together, these platforms provide a potentially useful first step towards developing lead therapeutic compounds that address disease complexity through rational polypharmacology.

Author SummaryMany modern drugs are developed with the goal of modulating a single cellular pathway or target. However, many drugs are, in fact,  dirty; they target multiple cellular pathways or targets. This phenomenon is known as multi-targeting or polypharmacology. While some strive to develop  cleaner therapeutics that eliminate secondary targets, recent work has shown that multi-targeting therapeutics have key advantages for a variety of diseases. However, while multi-targeting drugs that affect a precisely-defined profile of targets may be more effective, it is difficult to computationally predict which molecules have desirable target profiles. Here, we report the results of a competitive crowdsourcing project (the Multi-Targeting Drug DREAM Challenge), where we challenged participants to predict chemicals that have desired target profiles for cancer and neurodegenerative disease. Population differences in the length and early-life dynamics of telomeres among European pied flycatchers Telomere length and shortening rate are increasingly used as biomarkers for long-term costs in ecological and evolutionary studies because of their relationships with survival and fitness. Both early-life conditions and growth, and later-life stressors can create variation in telomere shortening rate. Studies on between- population telomere length and dynamics are scarce, despite the expectation that populations exposed to varying environmental constraints would present divergent telomere length patterns. The pied flycatcher (Ficedula hypoleuca) is a passerine bird breeding across Eurasia (from Spain to western Siberia) and migrating through the Iberian Peninsula to spend the non-breeding period in sub-Saharan Africa. Thus, different populations show marked differences in migration distance. We studied the large-scale variation of telomere length and early-life dynamics in the pied flycatcher by comparing six European populations across a north- south gradient (Finland, Estonia, England, and Spain) predicting negative effect of migration distance on adult telomere length, and of nestling growth on nestling telomere dynamics. There were clear population differences in telomere length, with English birds from mid-latitudes having the longest telomeres. Telomere length did not thus show consistent latitudinal variation and was not linearly linked to differences in migration distance. Early-life telomere shortening rate tended to vary between populations. Fast growth was associated with shorter telomeres in the early life, but faster nestling growth affected telomeres more negatively in northern than southern populations. While the sources of between-population differences in telomere-related biology remain to be more intensively studied, our study illustrates the need to expand telomere studies at the between-population level. Experimental evolution in TP53 deficient human gastric organoids recapitulates tumorigenesis The earliest events during human tumor initiation are poorly characterized but may hold clues as to how to detect and prevent malignancy. Here we model occult pre-neoplasia by engineering TP53 deficiency, a common and early event in gastric cancer, into human gastric organoids. By performing experimental evolution in multiple clonally derived cultures over two years we define causal relationships between this common initiating genetic lesion and resulting phenotypes. TP53 loss elicited progressive aneuploidy, including copy number alterations and complex structural variants that are common in gastric cancers and which follow preferred temporal orders. Longitudinal single cell sequencing of TP53 deficient gastric organoids similarly indicates progression towards malignant transcriptional programs. Moreover, lineage tracing with expressed cellular barcodes demonstrates reproducible dynamics whereby initially rare subclones with shared transcriptional programs repeatedly attain clonal dominance. This powerful platform for experimental evolution exposes stringent selection, clonal interference and a striking degree of phenotypic convergence in pre-malignant epithelial organoids, implying that the earliest stages of tumorigenesis may be predictable while illuminating evolutionary constraints and barriers to malignant transformation.

Graphical abstract

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=135 SRC="FIGDIR/small/487529v2_ufig1.gif" ALT="Figure 1">
View larger version (63K):
org.highwire.dtl.DTLVardef@122debborg.highwire.dtl.DTLVardef@4cd2daorg.highwire.dtl.DTLVardef@1a720d6org.highwire.dtl.DTLVardef@15ccaa7_HPS_FORMAT_FIGEXP  M_FIG C_FIG Modeling microRNA-driven post-transcriptional regulation by using exon-intron split analysis (EISA) in pigs The contribution of microRNAs (miRNAs) to mRNA regulation has often been explored by post hoc selection of downregulated genes and determining whether they harbor binding sites for miRNAs of interest. This approach, however, does not discriminate whether these mRNAs are also downregulated at the transcriptional level. Here, we have characterized the transcriptional and post-transcriptional changes of mRNA expression in two porcine tissues: gluteus medius muscle of fasted and fed Duroc gilts and adipose tissue of lean and obese Duroc-Gottingen minipigs. Exon-intron split analysis (EISA) of RNA-seq data allowed us to identify downregulated mRNAs with high post-transcriptional signals in fed or obese states, and we assessed whether they harbor binding sites for upregulated miRNAs in any of these two physiological states. We found 26 downregulated mRNAs with high post-transcriptional signals in the muscle of fed gilts and 21 of these were predicted targets of upregulated miRNAs also in the fed state. For adipose tissue, 44 downregulated mRNAs in obese minipigs displayed high post-transcriptional signals, and 25 of these were predicted targets of miRNAs upregulated in the obese state. These results suggest that the contribution of miRNAs to mRNA repression is more prominent in the skeletal muscle system. Finally, we identified several genes that may play relevant roles in the energy homeostasis of the pig skeletal muscle (DKK2 and PDK4) and adipose (SESN3 and ESRRG) tissues. By differentiating transcriptional from post-transcriptional changes in mRNA expression, EISA provides a valuable view about the regulation of gene expression, complementary to canonical differential expression analyses. Unbiased integration of single cell transcriptome replicates Single cell transcriptomic approaches are becoming mainstream, with replicate experiments commonly performed with the same single cell technology. Methods that enable integration of these datasets by removing batch effects while preserving biological information are required for unbiased data interpretation. Here we introduce Canek for this purpose. Canek leverages information from mutual nearest neighbor to combine local linear corrections with cell-specific non-linear corrections within a fuzzy logic framework. Using a combination of real and synthetic datasets, we show that Canek corrects batch effects while introducing the least amount of bias compared with competing methods. Canek is computationally efficient and can easily integrate thousands of single-cell transcriptomes from replicated experiments. Creatine transporter deficient rat model show motor dysfunction, cerebellar alterations and muscle creatine deficiency without muscle atrophy Creatine (Cr) is a nitrogenous organic acid and plays roles as fast phosphate energy buffer to replenish ATP, osmolyte, antioxidant, neuromodulator, and as a compound with anabolic and ergogenic properties in muscle. Cr is taken from the diet or endogenously synthetized by the enzymes AGAT and GAMT, and specifically taken up by the transporter SLC6A8. Loss-of-function mutations in the genes encoding for the enzymes or the transporter cause Cerebral Creatine Deficiency Syndromes (CCDS). CCDS are characterized by brain Cr deficiency, intellectual disability with severe speech delay, behavioral troubles, epilepsy and motor dysfunction. Among CCDS, the X-linked Cr transporter deficiency (CTD) is the most prevalent with no efficient treatment so far. Different animal models of CTD show reduced brain Cr levels, cognitive deficiencies and together they cover other traits similar to those of patients. However, motor function was poorly explored in CTD models and some controversies in the phenotype exist in comparison with CTD patients. Our recently described Slc6a8Y389C knock-in (KI) rat model of CTD showed mild impaired motor function, morphological alterations in cerebellum, reduced muscular mass, Cr deficiency and increased guanidinoacetate content in muscle, although no consistent signs of muscle atrophy. Our results indicate that such motor dysfunction co-occurred with both nervous and muscle dysfunction, suggesting that muscle strength and performance as well as neuronal connectivity might be affected by this Cr deficiency in muscle and brain. Face dissimilarity judgements are predicted by representational distance in morphable and image-computable models Human vision is attuned to the subtle differences between individual faces. Yet we lack a quantitative way of predicting how similar two face images look, or whether they appear to show the same person. Principal-components-based 3D morphable models are widely used to generate stimuli in face perception research. These models capture the distribution of real human faces in terms of dimensions of physical shape and texture. How well does a "face space" defined to model the distribution of faces as an isotropic Gaussian explain human face perception? We designed a behavioural task to collect dissimilarity and same/different identity judgements for 232 pairs of realistic faces. The stimuli densely sampled geometric relationships in a face space derived from principal components of 3D shape and texture (Basel Face Model, BFM). We then compared a wide range of models in their ability to predict the data, including the BFM from which faces were generated, a 2D morphable model derived from face photographs, and image-computable models of visual perception. Euclidean distance in the BFM explained both similarity and identity judgements surprisingly well. In a comparison against 14 alternative models, we found that BFM distance was competitive with representational distances in state-of-the-art image-computable deep neural networks (DNNs), including a novel DNN trained on BFM identities. Models describing the distribution of facial features across individuals are not only useful tools for stimulus generation. They also capture important information about how faces are perceived, suggesting that human face representations are tuned to the statistical distribution of faces. Prediction of survival after partial hepatectomy using a physiologically based pharmacokinetic model of indocyanine green liver function tests The evaluation of hepatic function and functional capacity of the liver are essential tasks in hepatology as well as in hepatobiliary surgery. Indocyanine green (ICG) is a widely applied test compound that is used in clinical routine to evaluate hepatic function. Important questions for the functional evaluation with ICG in the context of hepatectomy are how liver disease such as cirrhosis alters ICG elimination, and if postoperative survival can be predicted from preoperative ICG measurements. Within this work a physiologically based pharmacokinetic (PBPK) model of ICG was developed and applied to the prediction of the effects of a liver resection under various degrees of cirrhosis. For the parametrization of the computational model and validation of model predictions a database of ICG pharmacokinetic data was established. The model was applied (i) to study the effect of liver cirrhosis and liver resection on ICG pharmacokinetics; and (ii) to evaluate the model-based prediction of postoperative ICG-R15 as a measure for postoperative outcome. Key results are the accurate prediction of changes in ICG pharmacokinetics caused by liver cirrhosis and postoperative changes of ICG-elimination after liver resection, as validated with a wide range of data sets. Based on the PBPK model, individual survival after liver resection could be classified, demonstrating its potential value as a clinical tool. 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. RORγt serine 182 tightly regulates T cell heterogeneity to maintain mucosal homeostasis and restrict tissue inflammation Intestine homeostasis is maintained by the delicate balance of Th17 effector cells and Treg cells. Dysregulation of these cell populations contributes to inflammation, tissue damage, and chronic conditions. ROR{gamma}t is essential for the differentiation of Th17 and a subset of Treg (ROR{gamma}t+ Treg) cells involved in intestinal inflammation. ROR{gamma}t belongs to the nuclear receptor family of transcription factors with hinge regions that are highly flexible for co-activator/co-repressor interactions. Serine 182 at the hinge region of ROR{gamma}t is phosphorylated. This study aims to uncover how S182 on ROR{gamma}t contributes to mucosal homeostasis and diseases. We used CRISRP technology to generate a phosphor-null knock-in mutant mouse line (ROR{gamma}tS182A) to assess its role in intestine physiology. scRNA-seq was performed on WT and ROR{gamma}tS182A cohoused littermates to evaluate colonic T cell heterogeneity under steady state and colitis settings. Single-cell transcriptomics revealed that ROR{gamma}tS182 maintains colonic T cell heterogeneity under steady state, without interfering T cell development and differentiation. In inflamed tissues, ROR{gamma}tS182 simultaneously restricts IL-1{beta}-mediated Th17 activities and promotes anti-inflammatory cytokine IL-10 production in LT-like Treg cells. Phospho-null ROR{gamma}tS182A knock-in mice challenged with DSS induced colitis and EAE experienced delayed recovery and exacerbated pathology. The double switch role of ROR{gamma}tS182 is critical in resolving T cell-mediated inflammation and provides a potential therapeutic target to combat autoimmune diseases. Natural variations in reproductive timing of brood size and X-chromosome nondisjunction in Caenorhabditis elegans Reproductive timing underlies self-sperm production and X-chromosome nondisjunction rate in Caenorhabditis elegans. These traits can be measured as brood size and male-production rate, and brood size decreases and male-production rate increases as worms age. These phenotypic changes occur simultaneously; however, whether these changes share common underlying genetic architectures still remains unclear. To enhance our understanding of reproductive timing in C. elegans, we measured and compared these two reproductive timing-related traits of 96 wild strains during early, late and total reproductive periods. We found that the two traits exhibited natural phenotypic variation with few outliers and a similar reproductive timing pattern as previous reports. The brood size and the male-production rate were not genetically correlated, implying that they might have different genetic architectures and that meiosis quality can be maintained despite increased progeny-production in the late reproductive period caused by more self-sperm. We also identified loci and candidate genetic variants significantly associated with male-production rate in the late and total reproductive periods. Our results provide an insight into life history traits in wild C. elegans strains. 2,3-cAMP treatment mimics stress molecular response in Arabidopsis thaliana The role of the RNA degradation product 2,3-cyclic adenosine monophosphate (2,3-cAMP) is poorly understood. Recent studies have identified 2,3-cAMP in plant material and determined its role in stress signaling. The level of 2,3-cAMP increases upon wounding, dark, and heat, and 2,3-cAMP by binding to an RNA-binding protein, Rbp47b, promotes stress granule (SG) assembly. To gain further mechanistic insight into 2,3-cAMP function, we used a multi-omics approach combining transcriptomics, metabolomics, and proteomics to dissect Arabidopsis response to 2,3-cAMP treatment. We demonstrated that 2,3-cAMP is metabolized into adenosine, suggesting that the well-known cyclic nucleotide-adenosine pathway from human cells might also exist in plants. Transcriptomic analysis revealed only minor overlap between 2,3-cAMP-and adenosine-treated plants, suggesting that these molecules act through independent mechanisms. Treatment with 2,3-cAMP changed the levels of hundreds of transcripts, proteins, and metabolites, many previously associated with plant stress responses including protein and RNA degradation products, glucosinolates, chaperones and SG components. Finally, we demonstrated that 2,3-cAMP treatment influences the movement of processing bodies, supporting the role of 2,3-cAMP in the formation and motility of membraneless organelles. Preprints in motion: tracking changes between preprint posting and journal publication during a pandemic Amidst the COVID-19 pandemic, preprints in the biomedical sciences are being posted and accessed at unprecedented rates, drawing widespread attention from the general public, press and policymakers for the first time. This phenomenon has sharpened longstanding questions about the reliability of information shared prior to journal peer review. Does the information shared in preprints typically withstand the scrutiny of peer review, or are conclusions likely to change in the version of record? We assessed preprints from bioRxiv and medRxiv that had been posted and subsequently published in a journal through 30th April 2020, representing the initial phase of the pandemic response. We utilised a combination of automatic and manual annotations to quantify how an article changed between the preprinted and published version. We found that the total number of figure panels and tables changed little between preprint and published articles. Moreover, the conclusions of 7.2% of non-COVID-19-related and 17.2% of COVID-19-related abstracts undergo a discrete change by the time of publication, but the majority of these changes do not qualitatively change the conclusions of the paper. Improved radiation expression profiling in blood by sequential application of sensitive and specific gene signatures PurposeCombinations of expressed genes can discriminate radiation-exposed from normal control blood samples by machine learning based signatures (with 8 to 20% misclassification rates). These signatures can quantify therapeutically-relevant as well as accidental radiation exposures. The prodromal symptoms of Acute Radiation Syndrome (ARS) overlap those present in Influenza and Dengue Fever infections. Surprisingly, these human radiation signatures misclassified gene expression profiles of virally infected samples as false positive exposures. The present study investigates these and other confounders, and then mitigates their impact on signature accuracy.

MethodsThis study investigated recall by previous and novel radiation signatures independently derived from multiple Gene Expression Omnibus datasets on common and rare non-malignant blood disorders and blood-borne infections (thromboembolism, S. aureus bacteremia, malaria, sickle cell disease, polycythemia vera, and aplastic anemia). Normalized expression levels of signature genes are used as input to machine learning-based classifiers to predict radiation exposure in other hematological conditions.

ResultsExcept for aplastic anemia, these blood-borne disorders modify the normal baseline expression values of genes present in radiation signatures, leading to false-positive misclassification of radiation exposures in 8 to 54% of individuals. Shared changes, predominantly in DNA damage response and apoptosis-related gene transcripts in radiation and confounding hematological conditions, compromise the utility of these signatures for radiation assessment. These confounding conditions (sickle cell disease, thromboembolism, S. aureus bacteremia, malaria) induce neutrophil extracellular traps, initiated by chromatin decondensation, DNA damage response and fragmentation followed by programmed cell death. Riboviral infections (for example, Influenza or Dengue fever) have been proposed to bind and deplete host RNA binding proteins, inducing R-loops in chromatin. R-loops that collide with incoming replication forks can result in incompletely repaired DNA damage, inducing apoptosis and releasing mature virus. To mitigate the effects of confounders, we evaluated predicted radiation-positive samples with novel gene expression signatures derived from radiation-responsive transcripts encoding secreted blood plasma proteins whose expression levels are unperturbed by these conditions.

ConclusionsThis approach identifies and eliminates misclassified samples with underlying hematological or infectious conditions, leaving only samples with true radiation exposures. Diagnostic accuracy is significantly improved by selecting genes that maximize both sensitivity and specificity in the appropriate tissue using combinations of the best signatures for each of these classes of signatures. Estrogen-related receptor alpha and Rplp1-dependent translation coordinately regulate starvation response and decrease NASH progression Although general translation declines during fasting, maintaining the translation of a subset or proteins is necessary for metabolic homeostasis and cell viability. Using unbiased proteome analysis of hepatic cells during starvation, we identified a novel pathway in which Esrra-mediated transcription of Rplp1-dependent translation of lysosomal proteins declined during early starvation and recovered after prolonged starvation to restore autophagy-lysosome function. Interestingly, hepatic Esrra-Rplp1-dependent translation rate of lysosomal proteins also was impaired in patients and mice with non-alcoholic steatohepatitis (NASH), and translational response to starvation was dysregulated in mice with NASH. Remarkably, activation of Esrra pharmacologically, genetically, or by alternate day fasting restored protein translation, increased expression of lysosomal proteins, induced autophagy, and reduced lipotoxicity, inflammation, and fibrosis in cell culture and in vivo models of NASH. Thus, hepatic Esrra is essential for ribosome-dependent translation of lysosomal proteins during starvation, and prevention of lipotoxicity and progression in NASH.

Lay summaryFasting for weight loss improves non-alcoholic fatty liver disease (NAFLD) and non-alcoholic steatohepatitis (NASH); however, the mechanism is not well understood. Our study shows that a nuclear protein, estrogen related receptor alpha (Esrra), increases ribosome-mediated translation of autophagy and lysosome proteins during chronic starvation to maintain essential metabolic pathways for cell survival. Surprisingly, this translational pathway is impaired during NASH with reduced lysosome-autophagy activity accompanied by increased inflammation and fibrosis gene expression in the liver. Pharmacologic, genetic and dietary activation of Esrra decreases lipid-mediated toxicity in liver cells as well as inflammation and fibrosis in livers from mice with NASH. These findings suggest that the Esrra-ribosome-lysosome pathway is important for liver response to fasting and NASH and thus may be a good therapeutic target for the treatment of NASH. 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. Intertissue mechanical interactions shape the olfactory circuit in zebrafish While the chemical signals guiding neuronal migration and axon elongation have been extensively studied, the influence of mechanical cues on these processes remains poorly studied in vivo. Here, we investigate how mechanical forces exerted by surrounding tissues steer neuronal movements and axon extension during the morphogenesis of the olfactory placode in zebrafish. We mainly focus on the mechanical contribution of the adjacent eye tissue, which develops underneath the placode through extensive evagination and invagination movements. Using quantitative analysis of cell movements and biomechanical manipulations, we show that the developing eye exerts lateral traction forces on the olfactory placode through extracellular matrix, mediating proper morphogenetic movements and axon extension within the placode. Our data shed new light on the key participation of intertissue mechanical interactions in the sculpting of neuronal circuits. Hierarchical clustering optimizes the tradeoff between compositionality and expressivity of task structures for flexible reinforcement learning A hallmark of human intelligence, but challenging for reinforcement learning (RL) agents, is the ability to compositionally generalise, that is, to recompose familiar knowledge components in novel ways to solve new problems. For instance, when navigating in a city, one needs to know the location of the destination and how to operate a vehicle to get there, whether it be pedalling a bike or operating a car. In RL, these correspond to the reward function and transition function, respectively. To compositionally generalize, these two components need to be transferable independently of each other: multiple modes of transport can reach the same goal, and any given mode can be used to reach multiple destinations. Yet there are also instances where it can be helpful to learn and transfer entire structures, jointly representing goals and transitions, particularly whenever these recur in natural tasks (e.g., given a suggestion to get ice cream, one might prefer to bike, even in new towns). Prior theoretical work has explored how, in model-based RL, agents can learn and generalize task components (transition and reward functions). But a satisfactory account for how a single agent can simultaneously satisfy the two competing demands is still lacking. Here, we propose a hierarchical RL agent that learns and transfers individual task components as well as entire structures (particular compositions of components) by inferring both through a non-parametric Bayesian model of the task. It maintains a factorised representation of task components through a hierarchical Dirichlet process, but it also represents different possible covariances between these components through a standard Dirichlet process. We validate our approach on a variety of navigation tasks covering a wide range of statistical correlations between task components and show that it can also improve generalisation and transfer in more complex, hierarchical tasks with goal/subgoal structures. Finally, we end with a discussion of our work including how this clustering algorithm could conceivably be implemented by cortico-striatal gating circuits in the brain. A novel form of macropinocytosis mediates ultra-rapid transfer of pathological alpha-synuclein to lysosomes The nervous system spread of alpha-synuclein fibrils leads to Parkinsons disease (PD) and other synucleinopathies, yet the mechanisms underlying internalization and cell-to-cell transfer are enigmatic. Here we use confocal and superresolution microscopy, subcellular fractionation and electron microscopy (EM) of immunogold labelled alpha-synuclein pre-formed fibrils (PFF) to demonstrate that this toxic protein species enters cells using a novel form of ultra-rapid macropinocytosis with transfer to lysosomes in as little as 2 minutes, an unprecedented cell biological kinetic for lysosomal targeting. PFF uptake circumvents classical endosomal pathways and is independent of clathrin. Immunogold-labelled PFF are seen at the highly curved inward edge of membrane ruffles, in newly formed macropinosomes, and in lysosomes. While many of the fibrils remain in lysosomes that continue to take up PFF for hours, a portion are transferred to neighboring naive cells on the external face of vesicles, likely exosomes. These data indicate that PFF uses a novel internalization mechanism as a component of cell-to-cell propagation. Response of soil and vegetation in a warm-temperate Pine forest to intensive biomass harvests, phosphorus fertilisation, and wood ash application. Background and AimsConcerns about climate change and carbon economy have prompted the promotion of alternative energy sources, including forest-based bioenergy. An evaluation of the environmental consequences of intensive harvests (stumps and roots, and also branches and foliage) for energy wood supply, and use of wood-ash recycling as a compensatory practice, helps in the evaluation of the use of forest biomass for energy production.

MethodsWe made use of records from a split-plot experimental site crossing four different intensities of biomass harvesting (Stem-Only Harvest [SOH], Aboveground Additional Harvest [AAH], Belowground Additional Harvest [BAH], and Whole-Tree Harvest [WTH]) and three compensation methods (control [C], wood ash application [A] and phosphorus fertilisation [P]) to evaluate, in the 11-years-old stand (maritime pine: Pinus pinaster) that followed the biomass exports of the former stand, their effects on nutrient budgets, tree growth, soil fertility, chemical properties and soil carbon. This site is located in a forest on a poor soil, under a warm temperate climate (SW France).

Key resultsDespite their low additional biomass exports (+10% for AAH to +34% for WTH), the non-conventional harvest practices exported much higher quantities of nutrients than the conventional SOH technique (for example +145% for N and K in WTH). Consequently, these treatments had negative effects on the soil nutritive status. Additional biomass harvests impacted the soil organic matter content, with negative effects on P-organic, soil cation exchange capacity, exchangeable Ca, and most extractible nutrients. However, data suggested that tree growth and foliage nutrient content had not yet been significantly impacted by harvest treatments, whereas tree nutritional status was improved by P-fertiliser or wood ash. As expected, we observed a positive effect of wood ash application on soil pH and nutrient content but, like additional harvests, wood ash application decreased the pool of soil organic carbon (~10% of the initial stock with ~7% of N-total losses).

ConclusionsOverall, this factorial experiment showed that exporting more forest biomass due to the additional harvesting of tree canopies, stumps and roots had negative consequences on the ecosystem biogeochemistry. Additional harvests have aggravated the poverty of the already oligotrophic soil, and decreased the soil organic carbon content. Importantly, applying nutrients as fertiliser or wood ash did not compensate for all the negative impacts of biomass exports and the method of wood ash recycling in forests could even decrease the soil organic carbon. Antibacterial potency of Type VI amidase effector toxins is dependent on substrate topology and cellular context Members of the bacterial T6SS amidase effector (Tae) superfamily of toxins are delivered between competing bacteria to degrade cell wall peptidoglycan. Although Taes share a common substrate, they exhibit distinct antimicrobial potency across different competitor species. To investigate the molecular basis governing these differences, we quantitatively defined the functional determinants of Tae1 from Pseudomonas aeruginosa PAO1 using a combination of nuclear magnetic resonance (NMR) and a high-throughput in vivo genetic approach called deep mutational scanning (DMS). As expected, combined analyses confirmed the role of critical residues near the Tae1 catalytic center. Unexpectedly, DMS revealed substantial contributions to enzymatic activity from a much larger, ring-like functional hot spot extending around the entire circumference of the enzyme. Comparative DMS across distinct growth conditions highlighted how functional contribution of different surfaces is highly context-dependent, varying alongside composition of targeted cell walls. These observations suggest that Tae1 engages with the intact cell wall network through a more distributed three-dimensional interaction interface than previously appreciated, providing an explanation for observed differences in antimicrobial potency across divergent Gram-negative competitors. Further binding studies of several Tae1 variants with their cognate immunity protein demonstrate that requirements to maintain protection from Tae activity may be a significant constraint on the mutational landscape of tae1 toxicity in the wild. In total, our work reveals that Tae diversification has likely been shaped by multiple independent pressures to maintain interactions with binding partners that vary across bacterial species and conditions. Hundreds of viral families in the healthy infant gut The gut microbiome (GM) is shaped through infancy and plays a major role in determining susceptibility to chronic inflammatory diseases later in life. Bacteriophages (phages) are known to modulate bacterial populations in numerous ecosystems, including the gut. However, virome data is difficult to analyse because it mostly consists of unknown viruses, i.e. viral dark matter. Here, we manually resolved the viral dark matter in the largest human virome study published to date. Fecal viromes from a cohort of 647 infants at 1 year of age were deeply sequenced and analysed through successive rounds of clustering and curation. We uncovered more than ten thousand viral species distributed over 248 viral families falling within 17 viral order-level clades. Most of the defined viral families and orders were novel and belonged to the Caudoviricetes viral class. Bacterial hosts were predicted for 79% of the viral species using CRISPR spacers, including those in metagenomes from the same fecal samples. While Bacteroides-infecting Crassphages were present, novel viral families were more predominant, including phages infecting Clostridiales and Bifidobacterium. Phage lifestyles were determined for more than three thousand caudoviral species. Lifestyles were homogeneous at the family level for 149 Caudoviricetes families, including 32 families that were found to be virulent, while 117 were temperate. Virulent phage families were more abundant but temperate ones were more diverse and widespread. Together, the viral families found in this study represent a major expansion of existing bacteriophage taxonomy. Chemomechanically voxelated niches for programmable histogenesis Stem-cell derived tissue models are commonly cultured under globally-delivered stimuli that trigger histogenesis via self-organizing activity. However, the culture of such tissue models is prone to stochastic behavior, limiting the reproducibility of cellular composition and resulting in non-physiological architectures. To overcome these shortcomings, we developed a method for printing cell niche microenvironments with microstructured cues that mediate local histogenic processes, including mechanosensing and differentiation of selected cell types. Microstructured cues include independently tunable mechano-chemical properties, with conjugated peptides, proteins, and morphogens across a range of Youngs moduli. By rationally designing niches, we mediate the structure of tissues derived from stem-cell-progenitor sources, including a bone-fat assembly from stromal mesenchyme, and embryonic tissues derived from hiPSC. We show that microstructured cues can recapitulate mechano-chemical signals resembling early embryonic histogenesis. This outcome includes a role for niche mechanics in human embryonic organization, where soft niche mechanics bias markers of mesendodermal differentiation and epithelial-to-mesenchymal-transition (EMT), as well as a demonstration of a material-mediated morphogen signaling centers able to induce foci of mesenchymal and EMT differentiation. Thus, microstructured materials can mediate local histogenic processes to enhance the structure and composition of tissue models.