Evolution of thermal physiology alters the projected range of threespine stickleback under climate change Species distribution models (SDMs) are widely used to predict range shifts but could be unreliable under climate change scenarios because they do not account for evolution. The thermal physiology of a species is a key determinant of range and thus incorporating thermal trait evolution into SDMs might be expected to alter projected ranges. We identified a genetic basis for physiological and behavioural traits that evolve in response to temperature change in natural populations of threespine stickleback (Gasterosteus aculeatus). Using these data, we created geographic range projections using a mechanistic niche area approach under two climate change scenarios. Under both scenarios, trait data was either static ( no evolution models), allowed to evolve at observed evolutionary rates ( evolution models), or allowed to evolve at a rate of evolution scaled by the trait variance that is explained by quantitative trait loci (QTL;  scaled evolution models). We show that incorporating these traits and their evolution substantially altered the projected ranges for a widespread panmictic marine population, with over 7-fold increases in area under climate change projections when traits are allowed to evolve. Evolution-informed SDMs should improve the precision of forecasting range dynamics under climate change, and aid in their application to management and the protection of biodiversity.
How much does it take to be old? Modelling the time since the last harvesting to infer the distribution of overmature forests in France AimThe distribution of overmature forests in metropolitan France is poorly known, with only a few well-studied prominent sites, and has never been evaluated countrywide. Here, we modelled French forest reserves time since the last harvesting operation - a proxy for forest maturity - then inferred the current statistical distribution of overmature forests (i.e. forests over 50 years without harvesting) in France.

LocationMetropolitan France

MethodsWe used inventories from forest reserves and managed forests to calibrate a generalised linear mixed model explaining the time since the last harvesting with selected structural attributes and environmental variables. We then projected this model on the independent National Forest Inventory dataset. We thus obtained an updated estimation of the proportion and a rough distribution of overmature forest stands in metropolitan France.

ResultsWe found that high basal area of very large trees, high volumes of standing and downed deadwood, high diversity of tree-related microhabitats and more marginally diversity of decay stages best characterized the time since the last harvesting. Volumes of stumps and high density of coppices translating legacy of past forest management also distinguished more overmature plots. Our projection yielded an estimated 3% of French forests over 50 years without harvesting mostly located in more inaccessible areas (i.e. mountainous areas) and a promising proportion of future overmature forests if left unharvested.

Main conclusionsOur study showed that the time since the last harvesting is a good proxy for a combination of stand structure attributes key in characterising overmature temperate forests. It gives the first robust statistical estimate of the proportion of overmature forests and may serve to report on their status in metropolitan France. Our method could be implemented at a larger spatial scale, notably in countries with accessible National Forest Inventory and calibration data, to produce indicators at international level.
swCAM: estimation of subtype-specific expressions in individual samples with unsupervised sample-wise deconvolution MotivationComplex biological tissues are often a heterogeneous mixture of several molecularly distinct cell or tissue subtypes. Both subtype compositions and expressions in individual samples can vary across different biological states or conditions. Computational deconvolution aims to dissect patterns of bulk gene expression data into subtype compositions and subtype-specific expressions. Typically, existing deconvolution methods can only estimate averaged subtype-specific expressions in a population, while detecting differential expressions or co-expression networks in particular subtypes requires unique subtype expression estimates in individual samples. Different from population-level deconvolution, however, individual-level deconvolution is mathematically an underdetermined problem because there are more variables than observations.

ResultsWe report a sample-wise Convex Analysis of Mixtures (swCAM) method that can estimate subtype proportions and subtype-specific expressions in individual samples from bulk tissue transcriptomes. We extend our previous CAM framework to include a new term accounting for between-sample variations and formulate swCAM as a nuclear-norm and{ell} 2,1-norm regularized matrix factorization problem. We determine hyperparameter values using a cross-validation scheme with random entry exclusion and obtain a swCAM solution using an efficient alternating direction method of multipliers. The swCAM is implemented in open-source R scripts. Experimental results on realistic simulation data show that swCAM can accurately estimate subtype-specific expressions in individual samples and successfully extract co-expression networks in particular subtypes that are otherwise unobtainable using bulk expression data. Application of swCAM to bulk-tissue data of 320 samples from bipolar disorder patients and controls identified changes in cell proportions, expression and coexpression modules in patient neurons. Mitochondria related genes showed significant changes suggesting an important role of energy dysregulation in bipolar disorder.

Availability and implementationThe R Scripts of swCAM is freely available at https://github.com/Lululuella/swCAM. A users guide and a vignette are provided.

Contactyuewang@vt.edu

Supplementary informationSupplementary data are available at Bioinformatics online. 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. Diet effects on mouse meiotic recombination: a warning for recombination studies Meiotic recombination is a critical process for sexually reproducing organisms. This exchange of genetic information between homologous chromosomes during meiosis is important not only because it generates genetic diversity, but also because it is often required for proper chromosome segregation. Consequently, the frequency and distribution of crossovers are tightly controlled to ensure fertility and offspring viability. However, in many systems it has been shown that environmental factors can alter the frequency of crossover events. Two studies in flies and yeast point to nutritional status affecting the frequency of crossing over. However, this question remains unexplored in mammals. Here we test how crossover frequency varies in response to diet in Mus musculus males. We use immunohistochemistry to estimate crossover frequency in multiple genotypes under two diet treatments. Our results indicate that while crossover frequency was unaffected by diet in some strains, other strains were sensitive even to small composition changes between two common laboratory chows. Therefore, recombination is both resistant and sensitive to certain dietary changes in a strain-dependent manner and, hence, this response is genetically determined. Our study is the first to report a nutrition effect on genome-wide levels of recombination. Moreover, our work highlights the importance of controlling diet in recombination studies and may point to diet as a potential source of variability among studies, which is relevant for reproducibility. Conformational dynamics and allosteric modulation of the SARS-CoV-2 spike Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infects cells through binding to angiotensin-converting enzyme 2 (ACE2). This interaction is mediated by the receptor-binding domain (RBD) of the viral spike (S) glycoprotein. Structural and dynamic data have shown that S can adopt multiple conformations, which controls the exposure of the ACE2-binding site in the RBD. Here, using single-molecule Forster resonance energy transfer (smFRET) imaging we report the effects of ACE2 and antibody binding on the conformational dynamics of S from the Wuhan-1 strain and the B.1 variant (D614G). We find that D614G modulates the energetics of the RBD position in a manner similar to ACE2 binding. We also find that antibodies that target diverse epitopes, including those distal to the RBD, stabilize the RBD in a position competent for ACE2 binding. Parallel solution-based binding experiments using fluorescence correlation spectroscopy (FCS) indicate antibody-mediated enhancement of ACE2 binding. These findings inform on novel strategies for therapeutic antibody cocktails. Robust and Remote Measurement of Heart Rate Based on A Surveillance Camera Arrhythmia is a marked symptom of many cardiovascular diseases (CVDs). The accurate and timely detection of heart rate can greatly reduce the harm of arrhythmia to people. However, it is still a challenge to robustly and remotely measure heart rate in daily life due to the changing environmental conditions during measurement, such as the varying light intensity, the movement of people, and the uncertain distance between the sensor to people. In this study, we propose a method to accurately measure human heart rate within a distance of 4.5 meters under different light intensities by simply using a surveillance camera. After a 20-second color video of a persons hand is captured by the camera, a method based on the Fast Fourier Transform (FFT) algorithm is designed to extract the blood volume pulse wave to calculate the heart rate. According to the comparison between the real heart rate and results measured by electrocardiography (ECG), the proposed method achieves an accuracy of 95.8% when the measurement is performed within a distance of 4.5 meters and 90% when within 5.0 meters. Our experiments show that when the illuminance varies between 100-1000 lux (lighting level indoor), we still get the correct results. Our experiments also demonstrate that the proposed method accurately obtain heart rate even when the light intensity is below 32 lux (300-500 lux in a workplace environment). The methods strong adaptability to changing environmental conditions makes it applicable to many scenarios, such as homes of the elderly, classrooms, and other public spaces. Neural Mechanism of Hunger Gated Food-seeking and Evaluating The physiological need for energy evokes motivated feeding behaviours that help to ensure survival. However, the neural mechanisms underlying the generation of food motivation remain poorly understood. We investigated these mechanisms by subdividing feeding-related motivated behaviours into food-seeking, evaluating, and swallowing. Micro-endoscopic results indicated that neurons containing leptin receptors (LepRs) in the lateral hypothalamus (LH) are the major food-specific subpopulation of LHGABA neurons. Optogenetic manipulation of LHLepR neurons bidirectionally regulated both food-seeking and evaluating. Furthermore, micro-endoscope data revealed that distinct LHLepR neurons encode seeking and evaluating. Computational modelling analysis demonstrated that LHLepR neurons encode motivation, whereas neurons containing agouti-related peptide and neuropeptide Y (AgRP/NPY) encode the need for food. Additionally, slice studies revealed that NPY decreases inhibitory input to LHLepR neurons via LHGABA interneurons. This mechanism explains the permissive gate role of hunger (food need) in seeking/evaluating motivation. Together, the present study provides a comprehensive neural mechanism of how physiological needs drive distinct motivated behaviours. 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. Extensive duplication and convergent evolution of an antimicrobial peptide gene 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. Diverse stem-chondrichthyan oral structures and evidence for an independently acquired acanthodid dentition The teeth of sharks famously form a series of transversely-organised files with conveyor-belt replacement that are borne directly on the jaw cartilages, in contrast to the dermal plate-borne dentition of bony fishes that undergoes site-specific replacement. A major obstacle in understanding how this system evolved is the poorly understood relationships of the earliest chondrichthyans and the profusion of morphologically and terminologically diverse bones, cartilages, splints and whorls that they possess. Here we use tomographic methods to investigate mandibular structures in several early branching  acanthodian-grade stem-chondrichthyans. We show that the dentigerous jaw bones of disparate genera of ischnacanthids are united by a common construction, being growing bones with non-shedding dentition. Mandibular splints, which support the ventro-lateral edge of the Meckels cartilage in some taxa, are formed from dermal bone and may be an acanthodid synapomorphy. We demonstrate that the teeth of Acanthodopsis are borne directly on the mandibular cartilage and that this taxon is deeply nested within an edentulous radiation, representing an unexpected independent origin of teeth. Many or even all of the range of unusual oral structures may be apomorphic, but they should nonetheless be considered when building hypotheses of tooth and jaw evolution, both in chondrichthyans and more broadly. The "LINC" in between Delta40p53-miRNA axis in the regulation of cellular processes Recently, we showed that {Delta}40p53, the translational isoform of p53, can inhibit cell growth independently of p53 by regulating microRNAs. Here, we explored the role of {Delta}40p53 in regulating the lncRNA-miRNA axis. Based on early research from our laboratory, we selected LINC00176 for further exploration. LINC00176 levels were affected by the overexpression and knockdown of {Delta}40p53. Further, under DNA damage, ER stress, and glucose deprivation, LINC00176 was upregulated in HCT116 p53-/- cells (harboring only {Delta}40p53) compared to HCT116 p53+/+ cells. Additionally, ChIP and RNA stability assays revealed that {Delta}40p53 transactivates LINC00176 transcriptionally and stabilizes it post-transcriptionally. We ectopically overexpressed and knocked down LINC00176 in HCT116 p53-/- cells, which affected proliferation, cell viability, and the expression of epithelial markers. Finally, RNA pulldown revealed that LINC00176 sequesters several putative miRNA targets. These results provide important insights into the pivotal role of {Delta}40p53 in regulating the lncRNA-miRNA axis.

Graphical Abstract

O_FIG O_LINKSMALLFIG WIDTH=171 HEIGHT=200 SRC="FIGDIR/small/489891v2_ufig1.gif" ALT="Figure 1">
View larger version (31K):
org.highwire.dtl.DTLVardef@1035f85org.highwire.dtl.DTLVardef@74185borg.highwire.dtl.DTLVardef@1eb3704org.highwire.dtl.DTLVardef@168ebab_HPS_FORMAT_FIGEXP  M_FIG O_FLOATNOGraphical Abstract:C_FLOATNO Model depicting the regulation of {Delta}40p53-LINC00176-miRNA axis

(A) In absence of {Delta}40p53, LINC00176 levels do not increase which leads to increased pool of free miRNAs which probably binds to several mRNA targets preventing their function, ultimately resulting in increased cell proliferation. (B) However, when {Delta}40p53 is present it upregulates LINC00176 levels both at the level of transcription and post-transcription, which sequesters oncogenic miRNAs preventing it to bind to their respective mRNA targets, which ultimately results in decreased cell proliferation.

C_FIG FACT regulates pluripotency through distal regulation of gene expression in murine embryonic stem cells The FACT complex is a conserved histone chaperone with essential roles in transcription and histone deposition. FACT is essential in pluripotent and cancer cells, but otherwise dispensable for most mammalian cell types. FACT deletion or inhibition can block induction of pluripotent stem cells, yet the mechanism through which FACT regulates cell fate decisions remains unclear. To determine this mechanism, we used inducible depletion of FACT subunit SPT16 in murine embryonic stem cells paired with genomic factor localization, nascent transcription, and chromatin accessibility analyses. Over a timecourse of SPT16 depletion, nucleosomes invade loci bound by master pluripotency factors and gene-distal DNaseI hypersensitive sites. Simultaneously, transcription of Pou5f1 (OCT4), Sox2, Nanog, and enhancer RNAs produced at the genes' associated enhancers are downregulated, suggesting that FACT regulates expression of the pluripotency factors themselves. We find that FACT maintains cellular pluripotency through a precise nucleosome-based regulatory mechanism for appropriate expression of both coding and non-coding transcripts associated with pluripotency. 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. ALLOPURINOL BLOCKS THE AORTIC ANEURYSM IN A MOUSE MODEL OF MARFAN SYNDROME VIA REDUCING AORTIC OXIDATIVE STRESS BackgroundThe pathogenesis and progression of aortic aneurysm in Marfan syndrome (MFS) involves dysregulated TGF-{beta} and nitric oxide signaling, altered hemodynamics, and biomechanical forces. Increasing evidence indicates that redox stress participates in MFS aortopathy, though its contribution is not well established. We reported elevated reactive oxygen species (ROS) formation and NADPH oxidase NOX4 upregulation in MFS patients and mice aortic samples. Here we address the contribution of xanthine oxidoreductase (XOR), which is transcribed as dehydrogenase (XDH) and then post-translationally modified to oxidase (XO), catabolizing mainly purines into uric acid and ROS.

Methods and ResultsWe report that in MFS patients, XOR protein levels increased in aortic samples. In MFS mice (Fbn1C1041G/+), XOR mRNA transcripts augmented in the aorta of young but not older MFS mice. Moreover, the enzymatic activity of XO form increased with respect to XDH. The administration of the XOR inhibitor allopurinol (ALO) blocked the progression of aortic root aneurysm in MFS mice. ALO was also protective when administrated before the appearance of aneurysm. MFS-induced elastic fiber fragmentation, fibrotic remodeling, nuclear translocation of pNRF2, and increased 3-nitrotyrosine levels in the aortic tunica media as well as endothelial dysfunction were all prevented by allopurinol. The large overproduction of H2O2 in MFS aorta did not occur in allopurinol-treated mice whose uric acid plasma levels remained unaltered.

ConclusionsThis study strengthens the concept that redox stress is an important determinant of aortic aneurysm formation and progression in MFS and supports a clinical trial for ALO in the pharmacological treatment of MFS aortopathy.

GRAPHICAL ABSTRACT

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=146 SRC="FIGDIR/small/464182v3_ufig1.gif" ALT="Figure 1">
View larger version (31K):
org.highwire.dtl.DTLVardef@eab02aorg.highwire.dtl.DTLVardef@1d3b39eorg.highwire.dtl.DTLVardef@1b54d0forg.highwire.dtl.DTLVardef@1d61ad0_HPS_FORMAT_FIGEXP  M_FIG C_FIG Spatial and Texture Analysis of Root System Distribution with Earth Mover's Distance (STARSEED) PurposeRoot system architectures are complex, multidimensional, and challenging to characterize effectively for agronomic and ecological discovery.

MethodsWe propose a new method, Spatial and Texture Analysis of Root SystEm distribution with Earth movers Distance (STARSEED), for comparing root system distributions that incorporates spatial information through a novel application of the Earth Movers Distance (EMD).

ResultsWe illustrate that the approach captures the response of sesame root systems for different genotypes and soil moisture levels. STARSEED provides quantitative and visual insights into changes that occur in root architectures across experimental treatments.

ConclusionSTARSEED can be generalized to other plants and provides insight into root system architecture development and response to varying growth conditions not captured by existing root architecture metrics and models. The code and data for our experiments are publicly available: https://github.com/GatorSense/STARSEED. Artificial seawater based long-term culture of colonial ascidians Tunicates are highly diverse marine invertebrate filter-feeders that are vertebrates closest relatives. These organisms, despite a drastically different body plan during their adulthood, have a tissue complexity related to that of vertebrates. Ascidians, which compose most of the Tunicata, are benthic sessile hermaphrodites that reproduce sexually through a motile tadpole larval stage. Over half of the known ascidians species are able to reproduce asexually by budding, typically leading to the formation of colonies where animals, called zooids, are interconnected through an external vascular system. In addition, colonial ascidians are established models for important biological processes including allorecognition, immunobiology, aging, angiogenesis and whole-body regeneration. However, the current paucity in breeding infrastructures limits the study of these animals to coastal regions.

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

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

HighlightsO_LIFirst artificial seawater based recirculating aquaculture for colonial ascidians
C_LIO_LIOver 3 years of continuous breeding
C_LIO_LISemi-automated setup with minimized maintenance
C_LIO_LIGood biomass production for strain propagation
C_LIO_LI4 different culture conditions for optimized breeding for species of interest
C_LI Macroscopic control of synchronous electrical signaling with chemically-excited gene expression Excitable cells can convert electrical signals into chemical outputs to facilitate the active transport of information across larger distances. This electrical-to-chemical conversion requires a tightly regulated expression of ion channels. Alterations of ion channel expression provide landmarks of numerous pathological diseases, such as cardiac arrhythmia, epilepsy, or cancer. Although the activity of ion channels can be locally regulated by external light or chemical stimulus, it remains challenging to coordinate the expression of ion channels on extended spatial-temporal scales in a non-invasive manner. Here, we have engineered yeast S. cerevisiae to read and convert local chemical concentrations into a dynamic electrical field distributed across cell populations. The core mechanism encodes a chemically-excitable dual-feedback gene circuit that precisely tunes the expression domain of potassium channels, globally coordinating cyclic firing of the plasma membrane potential (PMP). We demonstrate that this mechanism leverages an engineered constitutively open bacterial potassium channel KcsA to directly couple chemical stimuli with ion flux through gene expression and it can interface with the host ion channels through the pulsatile production of toxins. Our study provides a robust synthetic transcriptional toolbox underlying the conversion of local chemical environments into spatiotemporally organized electrical impulses for various cellular engineering, synthetic biology, and potential therapeutic applications. Protein cost minimization promotes the emergence of coenzyme redundancy Coenzymes distribute a variety of chemical moieties throughout cellular metabolism, participating in group (e.g., phosphate, acyl) and electron transfer. For a variety of reactions requiring acceptors or donors of specific resources, there often exist degenerate sets of molecules (e.g., NAD(H) and NADP(H)) that carry out similar functions. Although the physiological roles of various coenzyme systems are well established, it is unclear what selective pressures may have driven the emergence of coenzyme redundancy. Here we use genome-wide metabolic modeling approaches to decompose the selective pressures driving enzymatic specificity for either NAD(H) or NADP(H) in the metabolic network of Escherichia coli. We found that few enzymes are thermodynamically constrained to using a single coenzyme, and in principle, a metabolic network relying on only NAD(H) is feasible. However, structural and sequence analyses revealed widespread conservation of residues that retain selectivity for either NAD(H) or NADP(H), suggesting that additional forces may shape specificity. Using a model accounting for the cost of oxidoreductase enzyme expression, we found that coenzyme redundancy universally reduces the minimal amount of protein required to catalyze coenzyme-coupled reactions, inducing individual reactions to strongly prefer one coenzyme over another when reactions are near thermodynamic equilibrium. We propose that protein minimization generically promotes coenzyme redundancy, and that coenzymes typically thought to exist in a single pool (e.g., CoA) may exist in more than one form (e.g., dephospho-CoA).

Significance statementMetabolism relies on a small class of molecules (coenzymes) that serve as universal donors and acceptors of key chemical groups and electrons. Although metabolic networks crucially depend on structurally redundant coenzymes (e.g., NAD(H) and NADP(H)) associated with different enzymes, the criteria that led to the emergence of this redundancy remain poorly understood. Our combination of modeling, and structural and sequence analysis indicates that coenzyme redundancy is not essential for metabolism, but rather an evolved strategy promoting efficient usage of enzymes when biochemical reactions are near equilibrium. Our work suggests that early metabolism may have operated with fewer coenzymes, and that adaptation for metabolic efficiency may have driven the rise of coenzyme diversity in living systems. Extensive characterization of HIV-1 reservoirs reveals links to plasma viremia before and during analytical treatment interruption The HIV-1 reservoir is composed of cells harboring latent proviruses that are capable of contributing to viremia upon antiretroviral treatment (ART) interruption. Although this reservoir is known to be maintained by clonal expansion, the contribution of large, infected cell clones to residual viremia and viral rebound remains underexplored. Here, we conducted an extensive analysis on four ART-treated individuals who underwent an analytical treatment interruption (ATI). We performed subgenomic (V1-V3 env), near full-length proviral and integration site sequencing, and used multiple displacement amplification to sequence both the integration site and provirus from single HIV-infected cells. We found eight proviruses that could phylogenetically be linked to plasma virus obtained before or during the ATI. This study highlights a role for HIV-infected cell clones in the maintenance of the replication-competent reservoir and suggests that infected cell clones can directly contribute to rebound viremia upon ATI. Plucking a string or playing a G? Choice type impacts human reinforcement learning In reinforcement learning (RL) experiments, participants learn to make rewarding choices in response to different stimuli; RL models use outcomes to estimate stimulus-response values which change incrementally. RL models consider any response type indiscriminately, ranging from more concretely defined motor choices (e.g. pressing a key with the index finger), to more general choices that can be executed in a number of ways (e.g. getting dinner at the restaurant). But does the learning process vary as a function of the choice type? In Experiment 1, we show that it does: participants were slower and less accurate in learning correct choices of a general format compared to learning more concrete, motor actions. Using computational modeling, we show that two mechanisms contribute to this. First, there was evidence of irrelevant credit assignment: the values of motor actions interfered with the values of other choice dimensions, resulting in more incorrect choices when the correct response is not defined by a single motor action; second, information integration for relevant general choices was slower. In Experiment 2, we replicated and further extended the findings from Experiment 1, by showing that slowed learning was attributable to weaker working memory use, rather than slowed RL learning. In both experiments we ruled out the explanation that the difference in performance between two condition types was driven by difficulty/different levels of complexity. We conclude that defining a more abstract choice space used by multiple learning systems for credit assignment recruits executive resources, limiting how much such processes then contribute to fast learning. Physiological Substrates and Ontogeny-Specific Expression of the Ubiquitin Ligases MARCH1 and MARCH8 MARCH1 and MARCH8 are ubiquitin ligases that control the expression and trafficking of critical immunoreceptors. Understanding of their function is hampered by three major knowledge gaps: (i) it is unclear which cell types utilize these ligases; (ii) their level of redundancy is unknown; and (iii) most of their putative substrates have been described in cell lines, often overexpressing MARCH1 or MARCH8, and it is unclear which substrates are regulated by either ligase in vivo. Here we address these questions by systematically analyzing the immune cell repertoire of MARCH1- or MARCH8-deficient mice, and applying unbiased proteomic profiling of the plasma membrane of primary cells to identify MARCH1 and MARCH8 substrates. Only CD86 and MHC II were unequivocally identified as immunoreceptors regulated by MARCH1 and MARCH8, but each ligase carried out its function in different tissues. MARCH1 regulated MHC II and CD86 in professional and "atypical" antigen presenting cells of hematopoietic origin, whereas MARCH8 only operated in non-hematopoietic cells. Our results reveal that the range of cells constitutively endowed with antigen-presentation capacity is wider than generally appreciated. They also establish MARCH1 and MARCH8 as specialized regulators of CD4+ T cell immunity in two ontogenically distinct cellular compartments. An integrative approach leads to the discovery of a novel anti-leukemic peptide from human milk Chemotherapy in childhood leukemia is associated with late morbidity in leukemic survivors, while certain patient subsets are relatively resistant to standard chemotherapy. It is therefore important to identify new agents with sensitivity and selectivity towards leukemic cells, while having less systemic toxicity. Peptide-based therapeutics has gained much attention during the last few years. Here, we used an integrative workflow combining mass spectrometric peptide library construction, in silico anticancer peptide screening, and in vitro leukemic cell studies to discover a novel anti-leukemic peptide having 3+charges and alpha-helical structure, namely HMP-S7, from human breast milk. HMP-S7 showed cytotoxic activity against four distinct leukemic cell lines in a dose-dependent manner but had no effect on solid malignancies or representative normal cells. HMP-S7 induced leukemic cell death by penetrating the plasma membrane to enter the cytoplasm and cause leakage of lactate dehydrogenase, thus acting in a membranolytic manner. Importantly, HMP-S7 exhibited anti-leukemic effect against patient-derived leukemic cells ex vivo. In conclusion, HMP-S7 is a selective anti-leukemic peptide with promise which requires further validation in preclinical and clinical studies.

TeaserIn silico screening of naturally occurring human milk peptides discovers a new anticancer peptide that kills leukemic cells in vitro and ex vivo. Dynamics of SARS-CoV-2 mutations reveals regional-specificity and similar trends of N501 and High-Frequency mutation N501Y in different levels of control measures Coronavirus disease 2019 (COVID-19) is a contagious disease caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). This disease has spread globally, causing more than 161.5 million cases and 3.3 million deaths to date. Surveillance and monitoring of new mutations in the virus genome are crucial to our understanding of the adaptation of SARS-CoV-2. Moreover, how the temporal dynamics of these mutations is influenced by control measures and non-pharmaceutical interventions (NPIs) is poorly understood. Using 1 058 020 SARS-CoV-2 from sequenced COVID-19 cases from 98 countries (totaling 714 country-month combinations), we perform a normalization by COVID-19 cases to calculate the relative frequency of SARS-CoV-2 mutations and explore their dynamics over time. We found 115 mutations estimated to be present in more than 3 % of global COVID-19 cases and determined three types of mutation dynamics: High-Frequency, Medium-Frequency, and Low-Frequency. Classification of mutations based on temporal dynamics enable us to examine viral adaptation and evaluate the effects of implemented control measures in virus evolution during the pandemic. We showed that Medium-Frequency mutations are characterized by high prevalence in specific regions and/or in constant competition with other mutations in several regions. Finally, taking N501Y mutation as representative of High-Frequency mutations, we showed that level of control measure stringency negatively correlates with the effective reproduction number of SARS-CoV-2 with High-Frequency or not-High-Frequency and both follows similar trends in different levels of stringency. Optogenetic dissection of actomyosin-dependent mechanics underlying tissue fluidity Rapid epithelial tissue flows are essential to building and shaping developing embryos. However, it is not well understood how the mechanical properties of tissues and the forces driving them to flow are jointly regulated to accommodate rapid tissue remodeling. To dissect the roles of actomyosin in the mechanics of epithelial tissue flows, here we use two optogenetic tools, optoGEF and optoGAP, to manipulate Rho/Rho-kinase signaling and actomyosin contractility in the germband epithelium, which flows via convergent extension movements during Drosophila body axis elongation. The ability to perturb actomyosin in the tissue allows us to analyze the effects of actomyosin on cell rearrangements, tissue tensions, and tissue mechanical properties. We find that either optogenetic activation or deactivation of Rho1 signaling and actomyosin contractility at the apical surface of the germband disrupts cell rearrangements and tissue-level flows. Rho1 activation leads to poorly oriented rearrangements that are associated with a redistribution of myosin II from the junctional to the medial-apical domain, whereas Rho1 deactivation leads to fewer, slower cell rearrangements that are associated with decreased junctional and medial myosin. By probing mechanical tensions in the tissue using laser ablation and inferring tissue mechanical properties from cell packings, we find that actomyosin influences both the anisotropic forces that drive tissue flow and the mechanical properties of the tissue that resist flow, leading to complex relationships between actomyosin activity and tissue fluidity. Moreover, our results link the subcellular distribution of myosin II to tissue tension and cell packings, revealing how junctional and medial myosin have differential roles in promoting and orienting cell rearrangements during tissue flow. A novel Fiji/ImageJ plugin for the rapid analysis of blebbing cells When confined, cells have recently been shown to undergo a phenotypic switch to what has been termed, fast amoeboid (leader bleb-based) migration. However, as this is a nascent area of research, few tools are available for the rapid analysis of cell behavior. Here, we demonstrate that a novel Fiji/ImageJ-based plugin, Analyze_Blebs, can be used to quickly obtain cell migration parameters and morphometrics from time lapse images. As validation, we show that Analyze_Blebs can detect significant differences in cell migration and morphometrics, such as the largest bleb size, upon introducing different live markers of F-actin, including F-tractin and LifeAct tagged with green and red fluorescent proteins. We also demonstrate, using flow cytometry, that live markers increase total levels of F-actin. Furthermore, that F-tractin increases cell stiffness, which was found to correlate with a decrease in migration, thus reaffirming the importance of cell mechanics as a determinant of Leader Bleb-Based Migration (LBBM).

Highlight summaryA new plugin, Analyze_Blebs, enables the rapid analysis of cell migration and morphometrics of fast amoeboid cells. Morphometrics combined with cell stiffness are found to be predictive of confined, fast amoeboid (leader bleb-based), migration.