Mean-field approximation of network of biophysical neurons driven by conductance-based ion exchange Numerous network and whole brain modeling approaches make use of mean-field models. Their relative simplicity allows studying network dynamics at a large scale. They correspond to lumped descriptions of neuronal assemblies connected via synapses. mean-field models do not consider the ionic composition of the extracellular space, which can change in physiological and pathological conditions, with strong effects on neuron activity. Here we derive a mean-field model of a population of Hodgkin-Huxley type neurons, which links the neuronal intra- and extra-cellular ion concentrations to the mean membrane potential and the mean synaptic input in terms of the synaptic conductance. The model can generate various physiological brain activities including multi-stability at resting states, as well as pathological spiking and bursting behaviors, and depolarization block. The results from the analytical solution of the mean-field model agree with the mean behavior of numerical simulations of large-scale networks of neurons. The mean-field model is analytically exact for non-autonomous ion concentration variables and provides a mean-field approximation in the thermodynamic limit, for locally homogeneous mesoscopic networks of biophysical neurons driven by an ion-exchange mechanism. These results may provide the missing link between high-level neural mass approaches which are used in the brain network modeling and physiological parameters that drive the neuronal dynamics.

Significance StatementIn this study, we applied mathematical formalism to estimate the mean-field behaviors of a large neuronal ensemble taking into account the ion-exchange between the intracellular and extracellular space. The model demonstrates different brain activities including resting state, spiking behavior, and seizure, as a function of the extracellular ion concentration. The relevant parameter regime for different brain activities is extracted mainly in terms of ion concentration and heterogeneity of individual neurons in the network. This neural mass model enables studying the influence of changes in extracellular ionic conditions on whole brain dynamics in health and disease. The effect of external stimulus current and conductance-based coupling of neural masses are also analyzed.
Reading Modality Shapes Reading Network in Proficient Blind Readers of Braille The neural basis of reading is highly consistent across many languages and scripts. Are there alternative neural routes to reading? How does the sensory modality of symbols (tactile vs. visual) influence their neural representations? We examined these questions by comparing reading of visual print (sighted group, n=19) and tactile Braille (congenitally blind group, n=19). Blind and sighted readers were presented with written (words, consonant strings, non-letter shapes) and spoken stimuli (words, backward speech) that varied in word-likeness. Consistent with prior work, the ventral occipitotemporal cortex (vOTC) was active during Braille and visual reading. A posterior/anterior vOTC word-form gradient was observed only in sighted readers with more anterior regions preferring larger orthographic units (words). No such gradient was observed in blind readers. Consistent with connectivity predictions, in blind compared to sighted readers, posterior parietal cortices were recruited to a greater degree and contained word-preferring patches. Lateralization of Braille in blind readers was predicted by laterality of spoken language and reading hand. The effect of spoken language increased along a cortical hierarchy, whereas effect of reading-hand waned. These results suggested that the neural basis of reading is influenced by symbol modality and spoken language and support connectivity-based views of cortical function.
Reduced glutathione decreases cell adhesion and increases cell volume. Glutathione is the most abundant thiol in animal cells. Reduced glutathione (GSH) is a major intracellular antioxidant neutralizing free radicals and detoxifying electrophiles. It plays important roles in many cellular processes, including cell differentiation, proliferation, and apoptosis. In the present study we demonstrate that extracellular concentration of reduced glutathione markedly increases cell volume within few hours, in a dose-response manner. Pre-incubation of cells with BSO, the inhibitor of {gamma}-glutamylcysteine synthetase, responsible for the first step in intracellular glutathione synthesis did not change the effect of reduced glutathione on cell volume suggesting a mechanism limited to the interaction of extracellular reduced glutathione on cell membrane. Similarly, inhibition of {gamma}-glutamylcyclotransferase involved in intracellular glutamate production had no effect on the action of reduced glutathione. Oxidized glutathione exerted no effect on cell volume. Results show that reduced GSH decreases cell adhesion resulting in an increased cell volume. Since many cell types are able to export GSH, the present results suggest that this could be a fundamental self-regulation of cell volume, giving the cells a self-control on their adhesion proteins. The Functional Nanopore Screen: A Versatile High-throughput Assay to Study and Engineer Protein Nanopores in Escherichia coli Nanopores comprise a versatile class of membrane proteins that carry out a range of key physiological functions and are increasingly developed for different biotechnological applications. Yet, a capacity to study and engineer protein nanopores by combinatorial means has so far been hampered by a lack of suitable assays that combine sufficient experimental resolution with throughput. Addressing this technological gap, the Functional Nanopore (FuN) screen now provides a quantitative and dynamic read-out of nanopore assembly and function in E. coli. The assay is based on genetically-encoded fluorescent protein sensors that resolve the nanopore-dependent influx of Ca2+ across the inner membrane of E. coli. Illustrating its versatile capacity, the FuN screen is first applied to dissect the molecular features that underlie the assembly and stability of nanopores formed by the S2168 holin. In a subsequent step, nanopores are engineered by recombining the transmembrane module of S2168 with different ring-shaped oligomeric protein structures that feature defined hexa-, hepta- and octameric geometries. Library screening highlights substantial plasticity in the ability of the S2168 transmembrane module to oligomerize in alternative geometries while the functional properties of the resultant nanopores can be fine-tuned through the identity of the connecting linkers. Overall, the FuN screen is anticipated to facilitate both fundamental studies and complex nanopore engineering endeavors with many potential applications in biomedicine, biotechnology and synthetic biology. Inhibitors of ERp44, PDIA1, and AGR2 induce disulfide-mediated oligomerization of Death Receptors 4 and 5 and cancer cell death Breast cancer mortality remains unacceptably high, indicating a need for safer and more effective therapeutic agents. Disulfide bond Disrupting Agents (DDAs) were previously identified as a novel class of anticancer compounds that selectively kill cancers that overexpress the Epidermal Growth Factor Receptor (EGFR) or its family member HER2. DDAs kill EGFR+ and HER2+ cancer cells via the parallel downregulation of EGFR, HER2, and HER3 and activation/oligomerization of Death Receptors 4 and 5 (DR4/5). However, the mechanisms by which DDAs mediate these effects are unknown. Affinity purification analyses employing biotinylated-DDAs reveal that the Protein Disulfide Isomerase (PDI) family members AGR2, PDIA1, and ERp44 are DDA target proteins. Further analyses demonstrate that shRNA-mediated knockdown of AGR2 and ERp44, or expression of ERp44 mutants, enhance basal and DDA-induced DR5 oligomerization. DDA treatment of breast cancer cells disrupts PDIA1 and ERp44 mixed disulfide bonds with their client proteins. Together, the results herein reveal DDAs as the first small molecule, active site inhibitors of AGR2 and ERp44, and demonstrate roles for AGR2 and ERp44 in regulating the activity, stability, and localization of DR4 and DR5, and activation of Caspase 8. The maize gene maternal derepression of r1 (mdr1) encodes a DNA glycosylase that demethylates DNA and reduces siRNA expression in endosperm Demethylation of transposons can activate expression of nearby genes and cause imprinted gene expression in endosperm, and it is hypothesized to lead to expression of transposon siRNAs that reinforce silencing in the next generation through transfer either into egg or embryo. Here we describe maternal derepression of r1 (mdr1), which encodes a DNA glycosylase with homology to Arabidopsis DEMETER and which is partially responsible for demethylation of thousands of regions in endosperm. Instead of promoting siRNA expression in endosperm, MDR1 activity inhibits it. Methylation of most repetitive DNA elements in endosperm is not significantly affected by MDR1, with an exception of Helitrons. While maternally-expressed imprinted genes preferentially overlap with MDR1 demethylated regions, the majority of genes that overlap demethylated regions are not imprinted. Double mutant megagametophytes lacking both MDR1 and its close homolog DNG102 result in early seed failure, and double mutant microgametophytes fail pre-fertilization. These data establish DNA demethylation by glycosylases as essential in maize endosperm and pollen and suggest that neither transposon repression nor genomic imprinting are its main function in endosperm. The genome of the cereal pest Sitophilus oryzae: a transposable element haven BackgroundThe rice weevil Sitophilus oryzae is one of the most important agricultural pests, causing extensive damage to cereal in fields and to stored grains. S. oryzae has an intracellular symbiotic relationship (endosymbiosis) with the Gram-negative bacterium Sodalis pierantonius and is a valuable model to decipher host-symbiont molecular interactions.

ResultsWe sequenced the Sitophilus oryzae genome using a combination of short and long reads to produce the best assembly for a Curculionidae species to date. We show that S. oryzae has undergone successive bursts of transposable element (TE) amplification, representing 72% of the genome. In addition, we show that many TE families are transcriptionally active, and changes in their expression are associated with insect endosymbiotic state. S. oryzae has undergone a high gene expansion rate, when compared to other beetles. Reconstruction of host-symbiont metabolic networks revealed that, despite its recent association with cereal weevils (30 Kyear), S. pierantonius relies on the host for several amino acids and nucleotides to survive and to produce vitamins and essential amino-acids required for insect development and cuticle biosynthesis.

ConclusionsHere we present the genome of an agricultural pest beetle, which may act as a foundation for pest control. In addition, S. oryzae may be a useful model for endosymbiosis, and studying TE evolution and regulation, along with the impact of TEs on eukaryotic genomes. Stress-dependent cell stiffening by tardigrade tolerance proteins through reversible formation of cytoskeleton-like filamentous network and gel-transition Tardigrades are able to tolerate almost complete dehydration by entering a reversible ametabolic state called anhydrobiosis and resume their animation upon rehydration. Dehydrated tardigrades are exceptionally stable and withstand various physical extremes. Although trehalose and late embryogenesis abundant (LEA) proteins have been extensively studied as potent protectants against dehydration in other anhydrobiotic organisms, tardigrades produce high amounts of tardigrade-unique protective proteins. Cytoplasmic-abundant heat-soluble (CAHS) proteins are uniquely invented in the lineage of eutardigrades, a major class of the phylum Tardigrada and are essential for their anhydrobiotic survival. However, the precise mechanisms of their action in this protective role are not fully understood. In the present study, we first postulated the presence of tolerance proteins that form protective condensates via phase separation in a stress-dependent manner and searched for tardigrade proteins that reversibly form condensates upon dehydration-like stress. Through comprehensive analysis, we identified 336 such proteins, collectively dubbed "Desolvataion-induced ReversiblY condensing Proteins (DRYPs)". Unexpectedly, we rediscovered CAHS proteins as highly enriched in DRYPs, 3 of which were major components of DRYPs. We revealed that these CAHS proteins reversibly polymerize into many cytoskeleton-like filaments depending on hyperosmotic stress in cultured cells and undergo reversible gel-transition in vitro. CAHS filamentation increases cell stiffness to resist deformation and improves resistance to dehydration-like stress. The conserved putative helical C-terminal region is necessary and sufficient for filament formation by CAHS proteins, and mutations disrupting the secondary structure of this region impaired both the filament formation and the gel transition. On the basis of these results, we propose that CAHS proteins are novel cytoskeletal proteins that form filamentous networks and undergo gel-transition in a stress-dependent manner to provide on-demand physical stabilization of cell integrity against deformative forces during dehydration and also contribute to the exceptional physical stability in a dehydrated state. Oxford Nanopore R10.4 long-read sequencing enables near-perfect bacterial genomes from pure cultures and metagenomes without short-read or reference polishing Long-read Oxford Nanopore sequencing has democratized microbial genome sequencing and enables the recovery of highly contiguous microbial genomes from isolates or metagenomes. However, to obtain near-perfect genomes it has been necessary to include short-read polishing to correct insertions and deletions derived from homopolymer regions. Here, we show that Oxford Nanopore R10.4 can be used to generate near-perfect microbial genomes from isolates or metagenomes without shortread or reference polishing. Extracellular Matrix Hydrogels Promote Expression of Muscle-Tendon Junction Proteins Muscle and tendon injuries are prevalent and range from minor sprains and strains to traumatic, debilitating injuries. However, the interactions between these tissues during injury and recovery remain unclear. Three-dimensional tissue models that incorporate both tissues and a physiologically relevant junction between muscle and tendon may help understand how the two tissues interact. Here, we use tissue specific extracellular matrix (ECM) derived from muscle and tendon to determine how cells of each tissue interact with the microenvironment of the opposite tissue resulting in junction specific features. ECM materials were derived from the Achilles tendon and gastrocnemius muscle, decellularized, and processed to form tissue specific pre-hydrogel digests. ECM materials were unique in respect to protein composition and included many types of ECM proteins, not just collagens. After digestion and gelation, ECM hydrogels had similar complex viscosities which were less than type I collagen hydrogels at the same concentration. C2C12 myoblasts and tendon fibroblasts were cultured in tissuespecific ECM conditioned media or encapsulated in tissue-specific ECM hydrogels to determine cell-matrix interactions and the effects on a muscle-tendon junction marker, paxillin. ECM conditioned media had only a minor effect on upregulation of paxillin in cells cultured in monolayer. However, cells cultured within ECM hydrogels had 50-70% higher paxillin expression than cells cultured in type I collagen hydrogels. Contraction of the ECM hydrogels varied by the type of ECM used. Subsequent experiments with varying density of type I collagen (and thus contraction) showed no correlation between paxillin expression and the amount of gel contraction, suggesting that a constituent of the ECM was the driver of paxillin expression in the ECM hydrogels. In addition, the extracellular matrix protein type XXII collagen had similar expression patterns as paxillin, with smaller effect sizes. Using tissue specific ECM allowed for the de-construction of the cell-matrix interactions similar to muscletendon junctions to study the expression of MTJ specific proteins.

Impact StatementThe muscle-tendon junction is an important feature of muscle-tendon units; however, despite crosstalk between the two tissue types, it is overlooked in current research. Deconstructing the cell-matrix interactions will provide the opportunity to study significant junction specific features and markers that should be included in tissue models of the muscletendon unit, while gaining a deeper understanding of the natural junction. This research aims to inform future methods to engineer a more relevant multi-tissue platform to study the muscletendon unit. Combining Plasma Extracellular Vesicle Let-7b-5p, miR-184 and Circulating miR-22-3p Levels for NSCLC Diagnosis and Drug Resistance Prediction Low-dose computed tomography (LDCT) Non-Small Cell Lung (NSCLC) screening is associated with high false-positive rates, leading to unnecessary expensive and invasive follow ups. There is a need for minimally invasive approaches to improve the accuracy of NSCLC diagnosis. Plasma extracellular vesicle (EV) and circulating microRNA (miRNA) have been proposed as cancer screening biomarkers. However, the identification of highly sensitive and broadly predictive core miRNA signatures remains a challenge. Also, how these systemic and diverse miRNAs impact cancer drug response is not well understood. Using an integrative approach, we examined plasma EV and circulating miRNA isolated from NSCLC patients versus screening controls with a similar risk profile. We found that combining EV (Hsa-miR-184, Let-7b-5p) and circulating (Hsa-miR-22-3p) miRNAs abundance robustly discriminates between NSCLC patients and high-risk cancer-free controls. Diagnosed NSCLC patients harboring sensitizing mutations in epidermal growth factor receptor EGFR (T790M, L578R) are treated with Osimertinib, a potent tyrosine kinase inhibitor (TKI). Nearly all patients develop TKI resistance via complex mechanisms and progress. We found that Hsa-miR-22-3p, Hsa-miR-184, and Let-7b-5p functionally converge on WNT/{beta}catenin and mTOR/AKT signaling axes, known cancer therapy resistance signals. Targeting Hsa-miR-22-3p and Hsa-miR-184 desensitized EGFR-mutated (T790M, L578R) NSCLC cells to Osimertinib. These findings suggest that the expression levels of circulating hsa-miR-22-3p combined with EV hsa-miR-184 and Let-7b-5p levels potentially define a core biomarker signature for improving the accuracy of NSCLC diagnosis. Importantly, these biomarkers have the potential to enable prospective identification of patients who are at risk of responding poorly to Osimertinib alone but likely to benefit from Osimertinib/AKT blockade combination treatments. SARS-CoV-2 infection induces soluble platelet activation markers and PAI-1 in the early moderate stage of COVID-19 IntroductionCoagulation dysfunction and thromboembolism emerge as strong comorbidity factors in severe COVID-19. However, it is unclear when particularly platelet activation markers and coagulation factors dysregulated during the pathogenesis of COVID-19. Here, we sought to assess the levels of coagulation and platelet activation markers at moderate and severe stages of COVID-19 to understand the pathogenesis.

MethodsTo understand this, hospitalized COVID-19 patients with (severe cases that required intensive care) or without pneumonia (moderate cases) were recruited. Phenotypic and molecular characterizations were performed employing basic coagulation tests including PT, APTT, D-Dimer and TFPI. The flow cytometry-based multiplex assays were performed to assess FXI, anti-thrombin, prothrombin, fibrinogen, FXIII, P-selectin, sCD40L, plasminogen, tissue-plasminogen activator (tPA), plasminogen activator inhibitor-1 (PAI-1) and D-Dimer.

ResultsThe investigations revealed induction of plasma P-selectin and CD40 ligand (sCD40L) in moderate COVID-19 cases which were significantly abolished with the progression of COVID-19 severity. Moreover, a profound reduction in plasma tissue factor pathway inhibitor (TFPI) and FXIII were identified particularly in the severe COVID-19. Further analysis revealed fibrinogen induction in both moderate and severe patients. Interestingly, an elevated PAI-1 more prominently in moderate, and tPA particularly in severe COVID-19 cases were observed. Particularly, the levels of fibrinogen and tPA directly correlated with the severity of the disease.

ConclusionsIn summary, induction of soluble P-selectin, sCD40L, fibrinogen and PAI-1 in moderate COVID-19 cases suggests the activation of platelets and coagulation system before patients require intensive care. These findings would help in designing better thromboprophylaxis to limit the COVID-19 severity. Newly developed ad hoc molecular assays show how eDNA can witness and anticipate the monk seal recolonization of central Mediterranean The monk seal Monachus monachus is the most endangered pinniped worldwide and the only one found in the Mediterranean, where its distribution and abundance have suffered a drastic decline in the last few decades. Data on its status are scattered due to both its rarity and evasiveness, and records are biased towards occasional, mostly coastal, encounters. Nowadays molecular techniques allow us to detect and quantify minute amounts of DNA traces released in the environment (eDNA) by any organism. We present three qPCR-assays targeting the monk seal mitogenome. The assays were soundly tested on an extensive and diversified sample set (n=73), including positive controls from Madeira breeding population collected during the peak of abundance, and two opportunistic Mediterranean eDNA-sample collections (offshore/coastal) from on-going projects. Monk seal DNA was detected in 47.2% and 66.7% of the samples collected in the Tyrrhenian from a ferry platform (2018-2019) and in the Pelagie archipelago -Strait of Sicily- (2020) respectively, anticipating (up to 2 year) visual observations occurred subsequently in proximity of the sampled areas. In the Tyrrhenian, detection occurrence increased between 2018 and 2019. Monk seal DNA recoveries were commoner in night-time ferry-samples, suggesting nocturnal predatory activity in pelagic waters. The proposed technique provides a non-invasive and yet highly-sensitive tool for defining the monk seal actual distribution and home range, its recovery rate and pinpoint coastal/offshore localities where prioritizing conservation, research, citizen science and education initiatives. TeamTree analysis: a new approach to evaluate scientific production Advances in science and technology depend on the work of research teams and the publication of results through peer-reviewed articles representing a growing socio-economic resource. Current methods to mine the scientific literature regarding a field of interest focus on content, but the workforce credited by authorship remains largely unexplored, and appropriate measures of scientific production are debated. Here, a new bibliometric approach named TeamTree analysis is introduced that visualizes the development and composition of the workforce driving a field. A new citation-independent measure that scales with the H index estimates impact based on publication record, genealogical ties and collaborative connections. This author-centered approach complements existing tools to mine the scientific literature and to evaluate research across disciplines. An NLR Integrated Domain toolkit to identify plant pathogen effector targets Plant resistance genes (or NLR "Nod-like Receptors") are known to contain atypical domains procuring them with a decoy capacity. Some of these integrated domains (or ID) allow the plant to lure the virulence determinants ("effectors") of pathogens and triggering a specific NLR immune reaction.

In this work, our goal was to generate a library of known IDs that could be screened with plant pathogen effectors in order to identify putative new effector virulence targets and NLR-effector pairs.

We curated the IDs contained in NLRs from seven model and crop plant species. We cloned 52 IDs representing 31 distinct Pfam domains. This library was screened for interaction by yeast-two-hybrid with a set of 31 conserved Ralstonia solanacearum type III effectors. This screening and the further in planta interaction assay allowed us to identify three interactions, involving different IDs (kinase, DUF3542, WRKY) and two type III effectors (RipAE and PopP2).

PopP2 was found to physically interact with ID#85, an atypical WRKY domain integrated in the GmNLR-ID85 NLR protein from Soybean. Using a imaging method in living plant cells, we showed that PopP2 associates with ID#85 in the nucleus. But unlike the known WRKY-containing Arabidopsis RRS1-R NLR receptor, this newly identified soybean WRKY domain could not be acetylated by PopP2 and its atypical sequence (WRKYGKR) also probably renders it inefficient in plant immunity triggering.

This ID toolkit is available for screening with other plant pathogen effectors and should prove useful to discover new effectors targets and potentially engineer new plant resistance genes. Morphology of immatures of the thelytokous ant, Monomorium triviale Wheeler (Formicidae: Myrmicinae: Solenopsidini) with descriptions of the extraordinary last-instar queen larvae The ant genus Monomorium is one of the most species-rich but taxonomically problematic groups in the hyperdiverse subfamily Myrmicinae. An East Asian species, M. triviale Wheeler, produces both reproductive queens and sterile workers via obligate thelytokous parthenogenesis. Here, we describe the immature forms of M. triviale based on light and scanning electron microscopy observations, with a note on the striking caste dimorphism in the last larval instar. The last-instar queen larvae were easily recognized by their large size, "aphaenogastroid" body shape, and rows of doorknob-like tubercles on the lateral and dorsal body surface. This type of queen-specific structure has not been found in ants in general, let alone congeneric species found in Japan. In stark contrast to the queen larvae, worker larvae showed a "pheidoloid" body shape and a body surface similar to other ants. The worker larvae were estimated to have three instars, consistent with previously described congeners. The pupae of both castes had no cocoon, a characteristic commonly described in other Myrmicinae species. In total, the developmental period from egg to adult worker averaged 59 days under 25{degrees}C. We discuss possible functions of the tubercles of queen larvae based on previous studies. A quartz crystal microbalance method to quantify the size of hyaluronan and other glycosaminoglycans on surfaces Hyaluronan (HA) is a major component of peri- and extra-cellular matrices and plays important roles in many biological processes such as cell adhesion, proliferation and migration. The abundance, size distribution and presentation of HA dictate its biological effects and are also useful indicators of pathologies and disease progression. Methods to assess the molecular mass of free-floating HA and other glycosaminoglycans (GAGs) are well established. In many biological and technological settings, however, GAGs are displayed on surfaces, and methods to obtain the size of surface-attached GAGs are lacking. Here, we present a method to size HA that is end-attached to surfaces. The method is based on the quartz crystal microbalance with dissipation monitoring (QCM-D) and exploits that the softness and thickness of films of grafted HA increase with HA size. These two quantities are sensitively reflected by the ratio of the dissipation shift ({Delta}D) and the negative frequency shift (-{Delta}f) measured by QCM-D upon the formation of HA films. Using a series of size-defined HA preparations, ranging in size from [~]2 kDa tetrasaccharides to [~]1 MDa polysaccharides, we establish a monotonic yet non-linear standard curve of the {Delta}D/-{Delta}f ratio as a function of HA size, which reflects the distinct conformations adopted by grafted HA chains depending on their size and surface coverage. We demonstrate that the standard curve can be used to determine the mean size of HA, as well as other GAGs, such as chondroitin sulfate and heparan sulfate, of preparations of previously unknown size in the range from 1 to 500 kDa, with a resolution of better than 10%. For polydisperse samples, our analysis shows that the process of surface-grafting preferentially selects smaller GAG chains, and thus reduces the average size of GAGs that are immobilised on surfaces comparative to the original solution sample. Our results establish a quantitative method to size HA and other GAGs grafted on surfaces, and also highlight the importance of sizing GAGs directly on surfaces. The method should be useful for the development and quality control of GAG-based surface coatings in a wide range of research areas, from molecular interaction analysis to biomaterials coatings. Wolbachia load variation in Drosophila is more likely caused by drift than by host genetic factors Symbiosis is a continuum of long-term interactions ranging from mutualism to parasitism, according to the balance between costs and benefits for the protagonists. The density of endosymbionts is, in both cases, a key factor that determines both the transmission of symbionts and the host extended phenotype, and is thus tightly regulated within hosts. However, the evolutionary and molecular mechanisms underlying bacterial density regulation are currently poorly understood. In this context, the symbiosis between the fruit fly and its intracellular bacteria Wolbachia (wMelPop strain) is particularly interesting to study. Although vertically transmitted, the symbiont is pathogenic, and a positive correlation between virulence and wMelPop density is observed. In addition, the number of repeats of a bacterial genomic region -Octomom-is positively correlated with Wolbachia density, underlying a potential genetic mechanism that controls bacterial density. Interestingly, the number of repeats varies between host individuals, but most likely also within them. Such genetic heterogeneity within the host could promote conflicts between bacteria themselves and with the host, notably by increasing within-host competition between symbiont genotypes through a process analogous to the tragedy of the commons. To characterize the determinisms at play in the regulation of bacterial density, we first introgressed wMelPop in different genetic backgrounds of D. melanogaster, and found different density levels and Octomom copy numbers in each host lineage. To determine whether such variations reflect a host genetic determinism on density regulation through Octomom copy number selection, we replicated the introgressions and performed reciprocal crosses on the two Drosophila populations with the most extreme density levels. In both experiments, we detected an absence of directionality in the patterns of infection, associated with a strong instability of these patterns across generations. Given that bacterial density was highly correlated with Octomom copy numbers in all experiments, these results rather suggest a strong influence of drift and a random increase in the frequency of certain bacterial variants. We then discuss how drift, both on the symbiont population during transmission and on the host population, could limit the efficiency of selection in such a symbiotic system, and the consequences of drift on the regulation of density and composition of bacterial populations. 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. Intestine-enriched apolipoprotein b orthologs are required for stem cell differentiation and regeneration in planarians Lipid metabolism plays an instructive role in regulating stem cell state and differentiation. However, the roles of lipid mobilization and utilization in stem cell-driven regeneration are unclear. Planarian flatworms readily restore missing tissue due to injury-induced activation of pluripotent somatic stem cells called neoblasts. Here, we identify two intestine-enriched orthologs of apolipoprotein b, apob-1 and apob-2, which mediate transport of neutral lipid stores from the intestine to target tissues including neoblasts, and are required for tissue homeostasis and regeneration. Inhibition of apob function by RNAi causes head regression and lysis in uninjured animals, and delays body axis re-establishment and regeneration of multiple organs in amputated fragments. Furthermore, apob RNAi causes expansion of the population of differentiating neoblast progeny and dysregulates expression of genes enriched in differentiating and mature cells in eight major cell type lineages. We conclude that intestine-derived lipids serve as a source of metabolites required for neoblast differentiation. 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. Muscle Activation Profile During Perturbed Walking is Modulated According to Body State During an unexpected loss of balance, avoiding a fall requires people to readjust their footing rapidly and effectively. We investigated the activation patterns of the ankle and knee muscles, and muscle fiber type recruitment resulting from unannounced, medio-lateral (i.e., right/left) horizontal-surface walking perturbations in twenty healthy adults (27.00{+/-}2.79 years, 10 females). Surface electromyography (sEMG) total spectral power for specific frequency bands (40-60Hz, 60-150Hz, 150-250Hz, 250-400Hz and 400-1000Hz), from tibialis anterior (TA) and vastus lateralis (VL) muscles were analyzed. Compared to non-perturbed walking, we found a significant increase in the total spectral power of lower-extremity muscles during the first 3 seconds after perturbation. When two feet were on the ground in time of perturbation we found a different muscle fiber type recruitment pattern between VL and TA muscles. This was not significant for perturbations implemented when one foot was on the ground. Our findings suggest that muscle operating frequency is modulated in real time to fit body state and somatosensory input in the context of functional goal requirements such as a rapid change of footing in response to unexpected loss of balance in single and double-support phases of gait.

New & NoteworthyTo study muscle spectral profiles in response to loss of balance, we investigated the dynamics of muscle spectral power changes, across different frequency bands after unannounced mechanical perturbations during walking. We showed increased activation of high-frequency motor units of the lower-limb muscles, subside 3 seconds after perturbation. Differences in power increase of specific frequency bands suggest that muscle activation is modulated in real time to fit body state in the context of functional goal requirements. Towards a unified model of naive T cell dynamics across the lifespan Naive CD4 and CD8 T cells are cornerstones of adaptive immunity, but the dynamics of their establishment early in life and how their kinetics change as they mature following release from the thymus are poorly understood. Further, due to the diverse signals implicated in naive T cell survival, it has been a long-held and conceptually attractive view that they are sustained by active homeostatic control as thymic activity wanes. Here we employ multiple experimental systems to identify a unified model of naive CD4 and CD8 T cell population dynamics across the mouse lifespan. We infer that both subsets divide rarely and progressively increase their survival capacity with cell age. Strikingly, this simple model captures naive CD4 T cell dynamics throughout life. In contrast, we find that newly generated naive CD8 T cells are lost more rapidly during the first 3-4 weeks of life, likely due to increased recruitment into memory. We find no evidence for elevated division rates in neonates, or for feedback regulation of naive T cell numbers at any age. We show how confronting mathematical models with diverse datasets can reveal a quantitative and remarkably simple picture of naive T cell dynamics from birth into old age. Fishing for biodiversity by balanced harvesting Fisheries are damaging, and seemingly incompatible with the conservation of marine ecosystems. Yet fish are an important source of food, and support the lives of many people in coastal communities. This paper considers an idea that a moderate intensity of fishing, appropriately scaled across species, could help in maintaining biodiversity, rather than reducing it. The scaling comes from an intuition that rates of fishing mortality of species should be kept in line with production rates of the species, a notion known as balanced harvesting. This places species conservation and exploitation on an equal footing in a single equation, showing quantitatively the relative levels of fishing mortality that species of different abundance can support. Using a dynamic model of a marine ecosystem, we give numerical evidence showing for the first time that fishing, if scaled in this way, can protect rarer species, while allowing some exploitation of species with greater production. This works because fishing mortality rates, when scaled by production, are density-dependent. Such fishing, operating adaptively to follow species production rates over time, contains a feedback that would help to protect species from overfishing in the presence of uncertainty about how marine ecosystems work. Bone marrow haematopoietic stem cells influence liver homeostatic networks and cancer development In fatty liver disease, systemic homeostasis is perturbed. While pre-clinical models are used to understand its pathogenesis, translating this knowledge to patients is difficult. However, by focusing on the most preserved homeostasis systems between species and models, novel disease dimensions can be unearthed. We interrogated core liver gene co-expression networks in a mouse model of liver cancer following dietary challenge. Immunometabolic modules showed temporal changes under the influence of diet duration and aging. The behaviour of immune network in tumours mirrored their counterparts in non-tumour liver. A high immune response network was associated with a lower tumour burden in mice and humans. This module in mice was enriched for genes related to haematopoietic cell differentiation. Consistently, the bone marrow haematopoietic stem and progenitor cells response was reflective of the liver immune response. Linking haematopoiesis to hepatic homeostasis uncovers a hitherto unexplored dimension of tissue crosstalk that can inform pathogenesis. Decoding silence in free recall In experiments on free recall from lists of items, not all memory retrievals are necessarily reported. Previous studies investigated unreported retrievals by attempting to induce their externalization. We show that, without any intervention, their statistics may be directly estimated through a model-free analysis of inter-response times - the silent intervals between recalls. A delay attributable to unreported recalls emerges in three situations: if the final item was already recalled ("silent recency effect"); if the item that, within the list, follows the latest recalled item was already recalled ("silent contiguity effect"); and in sequential recalls within highly performing trials ("sequential slowdown"). We then turn to reproducing all these effects by a minimal model where the discarding of memories ("bouncing") occurs either if they are repetitious or, in strategically organized trials, if they are not sequential. Based on our findings, we propose various approaches to further probing the submerged dynamics of memory retrieval.