Bias correction for inverse variance weighting Mendelian randomization BackgroundInverse-variance weighted two-sample Mendelian randomization (IVW-MR) is the most widely used approach that utilizes genome-wide association studies (GWAS) summary statistics to infer the existence and the strength of the causal effect between an exposure and an outcome. Estimates from this approach can be subject to different biases due to the use of weak instruments and winners curse, which can change as a function of the overlap between the exposure and outcome samples.

MethodsWe developed a method (MRlap) that simultaneously considers weak instrument bias and winners curse, while accounting for potential sample overlap. Assuming spike-and-slab genomic architecture and leveraging LD-score regression and other techniques, we could analytically derive, reliably estimate, and hence correct for the bias of IVW-MR using association summary statistics only.

ResultsWe tested our approach using simulated data for a wide range of realistic settings. In all the explored scenarios, our correction reduced the bias, in some situations by as much as 30 folds. When applied to real data on obesity-related exposures, we observed significant differences between IVW-based and corrected effects, both for non-overlapping and fully overlapping samples. Additionally, our results are consistent with the fact that the strength of the biases will decrease as the sample size increases. Finally, we also showed that the overall bias is also dependent on the genetic architecture of the exposure, and traits with low heritability and/or high polygenicity are more strongly affected.

ConclusionsOur method not only reduces bias in causal effect estimation but also enables the use of much larger GWAS sample sizes, by allowing for potentially overlapping samples.
Dosage sensitivity in Pumilio1-associated diseases involves two distinct mechanisms Mutations in the RNA-binding protein (RBP) Pumilio1 (PUM1) can cause dramatically different phenotypes. Mild mutations that reduce PUM1 levels by 25% lead to a mild, adult-onset ataxia, whereas more severe mutations that reduce PUM1 levels 40-60% produce an infantile syndrome involving multiple developmental delays and seizures. Why this difference in expression should cause such different phenotypes has been unclear; known PUM1 targets are de-repressed to equal degrees in both diseases. We therefore sought to identify PUM1s protein partners in the murine brain. We identified a number of putative interactors involved in different aspects of RNA metabolism such as silencing, alternative splicing, and polyadenylation. We find that PUM1 haploinsufficiency alters the stability of several interactors and disrupts the regulation of targets of those interactors, whereas mild PUM1 loss only de-represses PUM1-specific targets. We validated these phenomena in patient-derived cell lines and show that normalizing PUM1 levels rescues the levels of interactors and their targets. We therefore propose that dosage sensitivity does not necessarily reflect a linear change in levels but can involve distinct mechanisms. Studying the interactors of RBPs in vivo will be necessary to understand their functions in neurological diseases.
Ribosome actuator via the SpoT/ppGpp pathway to mitigate gene overexpression burden Gene overexpression causes a burden on ribosomes, which decreases the expression of other genes and growth rate. Here, we engineer a genetic ribosome actuator that increases gene expression rate and growth rate by over 150% and 80%, respectively. The actuator expresses the hydrolysis domain of the SpoT enzyme (SpoTH) in a bacterial strain with high basal levels of ppGpp, thus increasing ribosome availability. We use this actuator in a feedforward control arrangement, wherein SpoTH is co-expressed with a gene of interest (GOI). Without SpoTH co-expression, activation of the GOI decreases the expression of a constitutive gene and growth rate by over 80% and 40%, respectively. By contrast, with SpoTH co-expression, the ribosome binding site (RBS) of SpoTH can be tuned to achieve inappreciable change of either the expression of the constitutive gene or growth rate when the GOI is activated to the same level. We finally employ the actuator to keep growth rate constant while expressing dCas9 to a level that would otherwise decrease growth rate by more than 40% due to toxicity. This actuator will be relevant in applications where high cellular protein concentrations are desired without growth defects or global reduction of gene expression. 3CAC: improving the classification of phages and plasmids in metagenomic assemblies using assembly graphs MotivationBacteriophages and plasmids usually coexist with their host bacteria in microbial communities and play important roles in microbial evolution. Accurately identifying sequence contigs as phages, plasmids, and bacterial chromosomes in mixed metagenomic assemblies is critical for further unravelling their functions. Many classification tools have been developed for identifying either phages or plasmids in metagenomic assemblies. However, only two classifiers, PPR-Meta and viralVerify, were proposed to simultaneously identify phages and plasmids in mixed metagenomic assemblies. Due to the very high fraction of chromosome contigs in the assemblies, both tools achieve high precision in the classification of chromosomes but perform poorly in classifying phages and plasmids. Short contigs in these assemblies are often wrongly classified or classified as uncertain.

ResultsHere we present 3CAC, a new three-class classifier that improves the precision of phage and plasmid classification. 3CAC starts with an initial three-class classification generated by existing classifiers and improves the classification of short contigs and contigs with low confidence classification by using proximity in the assembly graph. Evaluation on simulated metagenomes and on real human gut microbiome samples showed that 3CAC outperformed PPR-Meta and viralVerify in both precision and recall, and increased F1-score by 10-60 percentage points.

AvailabilityThe 3CAC software is available on https://github.com/Shamir-Lab/3CAC.

Contactrshamir@tau.ac.il

Supplementary informationSupplementary data are available at Bioinformatics online. ATP synthase is essential for high conductance but not low conductance Mitochondrial Permeability Transition A sudden increase in permeability of the mitochondrial inner membrane, known as mitochondrial permeability transition (PT), is the central event responsible for cell death and tissue damage in conditions such as stroke and heart attack. PT is caused by the opening of the Cyclosporin A (CSA) dependent calcium-induced large pore, the Permeability Transition Pore (PTP). The molecular details of PTP opening remain incompletely understood. We utilized holographic microscopy imaging to directly assess PTP induction in HAP1 living cells. We found that calcium stress caused CSA-sensitive depolarization of the mitochondrial inner membrane. PTP did not occur until after nearly complete mitochondrial membrane depolarization. In cells lacking c and {delta} subunits of the ATP synthase, we observed CSA-sensitive membrane depolarization but not high-conductance PTP. We propose that PTP activation occurs as a two-phase process, with the initial membrane depolarization followed by the second phase of large pore opening that requires ATP synthase. SARS-CoV-2 spike-glycoprotein processing at S1/S2 and S2' and shedding of the ACE2 viral receptor: roles of Furin and TMPRSS2 and implications for viral infectivity and cell-to-cell fusion The spike (S)-protein of SARS-CoV-2 binds ACE2 and requires proteolytic "priming" at PRRAR685{downarrow} into S1 and S2 (cleavage at S1/S2), and "fusion-activation" at a S2 site for viral entry. In vitro, Furin cleaved peptides mimicking the S1/S2 cleavage site more efficiently than at the putative S2, whereas TMPRSS2 inefficiently cleaved both sites. In HeLa cells Furin-like enzymes mainly cleaved at S1/S2 during intracellular protein trafficking, and S2 processing by Furin at KPSKR815{downarrow} was strongly enhanced by ACE2, but not for the optimized S2 KRRKR815{downarrow} mutant (S2), whereas individual/double KR815AA mutants were retained in the endoplasmic reticulum. Pharmacological Furin-inhibitors (Boston Pharmaceuticals, BOS-inhibitors) effectively blocked endogenous S-protein processing in HeLa cells. Furthermore, we show using pseudotyped viruses that while entry by a "pH-dependent" endocytosis pathway in HEK293 cells did not require Furin processing at S1/S2, a "pH-independent" viral entry in lung-derived Calu-3 cells was sensitive to inhibitors of Furin (BOS) and TMPRSS2 (Camostat). Consistently, these inhibitors potently reduce infectious viral titer and cytopathic effects, an outcome enhanced when both compounds were combined. Quantitative analyses of cell-to-cell fusion and spike processing revealed the key importance of the Furin sites for syncytia formation. Our assays showed that TMPRSS2 enhances fusion and proteolysis at S2 in the absence of cleavage at S1/S2, an effect that is linked to ACE2 shedding by TMPRSS2. Overall, our results indicate that Furin and TMPRSS2 play synergistic roles in generating fusion-competent S-protein, and in promoting viral entry, supporting the combination of Furin and TMPRSS2 inhibitors as potent antivirals against SARS-CoV-2.

IMPORTANCESARS-CoV-2 is the etiological agent of COVID-19 that resulted in >5 million deaths. The spike protein (S) of the virus directs infection of the lungs and other tissues by binding the angiotensin-converting enzyme 2 (ACE2) receptor. For effective infection, the S-protein is cleaved at two sites: S1/S2 and S2. Cleavage at S1/S2, induces a conformational change favoring the recognition of ACE2. The S2 cleavage is critical for cell-to-cell fusion and virus entry into host cells. Our study contributes to a better understanding of the dynamics of interaction between Furin and TMPRSS2 during SARS-CoV-2 entry and suggests that the combination of a non-toxic Furin inhibitor with a TMPRSS2 inhibitor could significantly reduce viral entry in lung cells, as evidenced by an average synergistic [~]95% reduction of viral infection. This represents a powerful novel antiviral approach to reduce viral spread in individuals infected by SARS-CoV-2 or future related coronaviruses. P5A-ATPase controls the ER translocation of Wnt in neuronal migration Wnt family are conserved secretory proteins required for developmental patterning and tissue homeostasis. The mechanisms underlying intracellular maturation and intercellular signal transduction of Wnt proteins have been extensively studied. However, how Wnt is targeted to the endoplasmic reticulum (ER) for processing and secretion remains elusive. Here we report that CATP-8/P5A-ATPase directs neuronal migration non-cell autonomously in C. elegans by regulating EGL-20/Wnt biogenesis. CATP-8 functions as a translocase to translocate EGL-20/Wnt nascent polypeptide into the ER by interacting with the hydrophobic core region of EGL-20 signal sequence. Such regulation of Wnt biogenesis by P5A-ATPase is conserved in human cells. These findings reveal physiological roles of P5A-ATPase in neural development and identify Wnt proteins as direct substrates of P5A-ATPase for ER translocation. Protein scaffold-based multimerization of soluble ACE2 efficiently blocks SARS-CoV-2 infection in vitro and in vivo Soluble ACE2 (sACE2) decoy receptors are promising agents to inhibit SARS-CoV-2, as their efficiency is less likely to be affected by common escape mutations in viral proteins. However, their success may be limited by their relatively poor potency. To address this challenge, we developed a large decoy library of sACE2 fusion proteins, generated with several protease inhibitors or multimerization tags. Among these decoys, multimeric sACE2 consisting of SunTag or MoonTag systems, which were originally utilized for signal amplification or gene activation systems, were extremely effective in neutralizing SARS-CoV-2 in pseudoviral systems and in clinical isolates. These novel sACE2 fusion proteins exhibited greater than 100-fold SARS-CoV-2 neutralization efficiency, compared to monomeric sACE2. SunTag or MoonTag in combination with a more potent version of sACE2, which has multiple point mutations for greater binding (v1), achieved near complete neutralization at a sub-nanomolar range, comparable with clinical monoclonal antibodies. Pseudoviruses bearing mutant versions of Spike, alpha, beta, gamma or delta variants, were also neutralized efficiently with SunTag or MoonTag fused sACE2(v1). Finally, therapeutic treatment of sACE2(v1)-MoonTag provided protection against SARS-CoV-2 infection in an in vivo mouse model. Overall, we suggest that the superior activity of the sACE2-SunTag or sACE2-MoonTag fusions is due to the greater occupancy of the multimeric sACE2 receptors on Spike protein as compared to monomeric sACE2. Therefore, these highly potent multimeric sACE2 decoy receptors may offer a promising treatment approach against SARS-CoV-2 infections.

One Sentence SummaryMultimerization of sACE2 markedly enhanced the neutralization of SARS-CoV-2 by blocking multiple viral spike proteins simultaneously. Genome editing excisase origins illuminated by somatic genome of Blepharisma Massive DNA excision occurs regularly in ciliates, ubiquitous microbial eukaryotes with somatic and germline nuclei in the same cell. Tens of thousands of internally eliminated sequences (IESs) scattered throughout a copy of the ciliate germline genome are deleted during development of the streamlined somatic genome. Blepharisma represents one of the two earliest diverging ciliate classes, and, unusually, has dual pathways of somatic nuclear development, making it ideal for investigating the functioning and evolution of these processes. Here, we report the somatic genome assembly of Blepharisma stoltei strain ATCC 30299 (41 Mb), arranged as numerous alternative telomere-capped minichromosomes. This genome encodes eight PiggyBac transposase homologs liberated from transposons. All are subject to purifying selection, but just one, the putative IES excisase, has a complete catalytic triad. We propose PiggyBac homologs were ancestral excisases that enabled evolution of extensive, natural genome editing. Coordination of two opposite flagella allows high-speed swimming and active turning of individual zoospores Phytophthora species cause diseases in a large variety of plants and represent a serious agricultural threat, leading, every year, to multibillion dollar losses. Infection occurs when these biflagellated zoospores move across the soil at their characteristic high speed and reach the roots of a host plant. Despite the relevance of zoospore spreading in the epidemics of plant diseases, characteristics of individual swimming of zoospores have not been fully investigated. It remains unknown about the characteristics of two opposite beating flagella during translation and turning, and the roles of each flagellum on zoospore swimming. Here, combining experiments and modeling, we show how these two flagella contribute to generate thrust when beating together, and identify the mastigonemes-attached anterior flagellum as the main source of thrust. Furthermore, we find that turning involves a complex active process, in which the posterior flagellum temporarily stops, while the anterior flagellum keeps on beating and changes its pattern from sinusoidal waves to power and recovery strokes, similar to Chlamydomonass breaststroke, to reorient its body to a new direction. Our study is a fundamental step towards a better understanding of the spreading of plant pathogens motile forms, and shows that the motility pattern of these biflagellated zoospores represents a distinct eukaryotic version of the celebrated "run-and-tumble" motility class exhibited by peritrichous bacteria. Social integration predicts survival in female white-faced capuchin monkeys Across multiple species of social mammals, a growing number of studies have found that that individual sociality is associated with survival. In long-lived species, like primates, lifespan is one of the main components of fitness. We used 18 years of data from the Lomas Barbudal Monkey project to quantify social integration in 11 capuchin (Cebus capucinus) groups and tested whether female survivorship was associated with females tendencies to interact with three types of partners: (1) all group members, (2) adult females, (3) adult males. We found strong evidence that females who engaged more with other females in affiliative interactions and foraged in close proximity experienced increased survivorship. We found some weak evidence that females might also benefit from engaging in more support in agonistic contexts with other females. These benefits were evident in models that account for the females rank and group size. Female interactions with all group members also increased survival, but the estimates of the effects were more uncertain. In interactions with adult males, only females who provided more grooming to males survived longer. The results presented here suggest that social integration may result in survival-related benefits. Females might enjoy these benefits through exchanging grooming for other currencies, such as coalitionary support or tolerance. NPEPPS regulates intracellular import and sensitivity to cisplatin by interaction with volume regulated anion channels Despite routine use of platinum-based chemotherapeutics in cancer treatment, there remains a need to improve efficacy and patient selection. Multi-omic assessment of human bladder cancer cell lines and their cisplatin resistant derivatives and whole-genome CRISPR screens identified NPEPPS, the puromycin-sensitive aminopeptidase as a novel driver of cisplatin resistance. NPEPPS depletion increased cisplatin import and sensitization of resistant cells in vitro and in vivo. Pharmacologic inhibition of NPEPPS in cells and chemoresistant, patient-derived tumor organoids improved response to cisplatin. NPEPPS was found in complex with volume regulated anion channel (VRAC) subunits LRRC8A and LRRC8D, whose loss is known to enhance resistance to cisplatin. Depletion of LRRC8A the only obligate subunit for normal VRAC function abrogated the effect of NPEPPS-mediated cisplatin import. Our findings describe the first mechanism by which VRACs can be targeted for therapeutic benefit.

HIGHLIGHTSO_LIMulti-omic screening found NPEPPS is a driver of cisplatin resistance in cancer
C_LIO_LINPEPPS is in protein complex with LRRC8A and D, volume regulated anion channel (VRAC) subunits
C_LIO_LILRRC8A and D loss increases cisplatin resistance and reduce its intracellular levels
C_LIO_LINPEPPS regulates cisplatin import and sensitivity through VRACs
C_LIO_LINPEPPS function is inhibited by small molecules offering potential for clinical translation
C_LI The spatial configuration of biotic interactions shapes coexistence-area relation-ships in an annual plant community The increase of species richness with area is a universal phenomenon on Earth. How-ever, this observation contrasts with our poor understanding of how these species-area relationships (SARs) emerge from the collective effects of area, spatial heterogeneity, and local interactions. By combining a structuralist approach with five years of empirical observations in a highly-diverse Mediterranean grassland, we show that spatial hetero-geneity plays a little role in the accumulation of species richness with area in our system. Instead, as we increase the sampled area more species combinations are realized, and they coexist mainly due to direct pairwise interactions rather than by changes in single-species dominance or by indirect interactions. We also identify a small set of transient species with small population sizes that are consistently found across spatial scales. These findings empirically support the importance of the architecture of species interactions together with stochastic events for driving coexistence- and species-area relationships. A Novel Organoid Model of In Vitro Spermatogenesis Using Human Induced Pluripotent Stem Cells Infertility is thought to be caused by genetic mutations and dysfunction in the cellular niche where spermatogenesis takes place. An understanding of the specialized cellular processes which drive spermatogenesis is needed to develop treatments; however, the development of in vitro systems to study these cells has been hindered by our reliance on rarely available human testicular tissues for research. Human induced pluripotent stem cells (hiPSCs) can be used to derive human testicular-like cells, and thus provide an avenue for the development of in vitro testicular model systems. Therefore, this study set out to engineer a human testicular tissue model using hiPSCs for the first time. We demonstrate the ability of hiPSC-derived testicular cells to self-organize and mature into testicular-like tissues using organoid culture. Moreover, we show that hiPSC-derived testicular organoids promote testicular somatic cell maturation and spermatogenesis up to the post-meiotic spermatid stage. These hiPSC-derived testicular organoids have the potential to replace rarely available primary testicular tissues to further infertility research in an in vitro setting. A potent human monoclonal antibody with pan-neutralizing activities directly dislocates S trimer of SARS-CoV-2 through binding both up and down forms of RBD The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has caused a global pandemic of novel corona virus disease (COVID-19). The neutralizing monoclonal antibodies (mAbs) targeting the receptor binding domain (RBD) of SARS-CoV-2 are among the most promising strategies to prevent and treat COVID-19. However, SARS-CoV-2 variants of concern (VOCs) profoundly reduced the efficacies of most of mAbs and vaccines approved for clinical use. Herein, we demonstrated mAb 35B5 efficiently neutralizes both wild-type (WT) SARS-CoV-2 and VOCs, including B.1.617.2 (delta) variant, in vitro and in vivo. Cryo-electron microscopy (cryo-EM) revealed that 35B5 neutralizes SARS-CoV-2 by targeting a unique epitope that avoids the prevailing mutation sites on RBD identified in circulating VOCs, providing the molecular basis for its pan-neutralizing efficacy. The 35B5-binding epitope could also be exploited for the rational design of a universal SARS-CoV-2 vaccine. 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. Untargeted metabolomic analyses reveal the diversity and plasticity of the specialized metabolome in seeds of different Camelina sativa genotypes Despite the essential role of Specialized Metabolites (SMs) in the interaction of plants with the environment, studying the ability of crop seeds to produce these protective compounds has been neglected. Furthermore, seeds produce a myriad of SMs providing an interesting model to investigate their diversity and plasticity. Camelina sativa gained a lot of interest in the past few years as rustic oil seed crop. A characterization of seed SM landscapes in six camelina genotypes grown in the field and harvested during five growing seasons has been undertaken in this work. This allowed a comprehensive annotation of seed SMs combining analyses that cluster SMs based on their chemical structures and co-accumulation patterns. These data showed broad effects of the environment on the stimulation of the seed-specialized metabolome. Among well annotated compounds, flavonols were identified as the metabolic class characterised by high plasticity, revealing significant variable accumulation according to the year and/or the genotype. Lastly, a deeper characterisation of primary metabolites and lipids in two selected genotypes has been performed. We showed that, in addition to flavonols, alkaloids and glucosinolates displayed a higher phenotypic plasticity with respect to most of the primary metabolites, including some sugars and major storage compounds such as fatty acids, proteins and most lipid classes (e.g. DAG, TAG), but similar plasticity compared to free aminoacids and carboxylic acids. This work highlighted major and unexplored effects of the environment on the seed specialized metabolome demonstrating that seeds exhibit a dynamic and plastic metabolism, with an impact on seed quality.

Significance statementSeeds produce a myriad of Specialized Metabolites (SMs) with an essential role in the interaction of plants with the environment. We characterized SM landscapes, primary metabolites and lipid composition in the seeds of camelina genotypes grown in the open field in five consecutive growing seasons. Our results showed the predominant effect of the environment on the regulation of the seed - specialized metabolome, with a potential impact on seed quality of camelina that may also occur in other oilseed crops. Universality of evolutionary trajectories under arbitrary competition dynamics The assumption of constant population size is central in population genetics. It led to a large body of results, that are robust to modeling choices and that have proven successful to understand evolutionary dynamics. In reality, allele frequencies and population size are both determined by the interaction between a population and the environment. Relaxing the constant-population assumption have two big drawbacks. It increases the technical difficulty of the analysis, and it requires specifying a mechanism for the saturation of the population size, possibly making the results contingent on model details. Here, we develop a framework that encompasses a great variety of systems with an arbitrary mechanism for population growth limitation. By using techniques based on scale separation for stochastic processes, we are able to calculate analytically properties of evolutionary trajectories, such as the fixation probability. Remarkably, these properties assume a universal form with respect to our framework, which depends on only three parameters related to the inter-generation timescale, the invasion fitness, and the carrying capacity of the strains. In other words, different systems, such as Lotka-Volterra or a chemostat model (contained in our framework), share the same evolutionary outcomes after a proper re-mapping of their parameters. An important and surprising consequence of our results is that the direction of selection can be inverted, with a population evolving to reach lower values of invasion fitness. Loss of SUMO-specific protease 2 causes isolated glucocorticoid deficiency by blocking adrenal cortex zonal transdifferentiation SUMOylation is a dynamic posttranslational modification, that provides fine-tuning of protein function involved in the cellular response to stress, differentiation, and tissue development. In the adrenal cortex, an emblematic endocrine organ that mediates adaptation to physiological demands, the SUMOylation gradient is inversely correlated with the gradient of cellular differentiation raising important questions about its role in functional zonation and the response to stress. Considering that SUMO-specific protease 2 (SENP2), a deSUMOylating enzyme, is upregulated by ACTH/PKA signalling within the zona Fasciculata (zF), we generated mice with adrenal-specific Senp2 loss to address these questions. Disruption of SENP2 activity in steroidogenic cells leads to specific hypoplasia of the zF, a blunted reponses to ACTH and isolated glucocorticoid deficiency. Mechanistically, overSUMOylation resulting from SENP2 loss shifts the balance between ACTH/PKA and WNT/{beta}-catenin signalling leading to repression of PKA activity and ectopic activation of {beta}-catenin. At the cellular level, this blocks transdifferentiation of {beta}-catenin-positive zona Glomerulosa cells into zF cells and sensitises them to premature apoptosis. Our findings indicate that the SUMO pathway is critical for adrenal homeostasis and stress responsiveness. 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. Malaria parasites utilize pyrophosphate to fuel an essential proton pump in the ring stage and the transition to trophozoite stage The malaria parasite relies on anaerobic glycolysis for energy supply when growing inside RBCs as its mitochondrion does not produce ATP. The ring stage lasts [~] 22 hours and is traditionally thought to be metabolically quiescent. However, recent studies show that the ring stage is active for several energy-costly processes including gene transcription/translation, protein export, and movement inside the RBC. It has remained unclear if a low glycolytic flux can meet the energy demand of the ring stage. Here we show that the metabolic by-product, pyrophosphate, is a critical energy source for the development of the ring stage and its transition to the trophozoite stage. During early phases of the asexual development, the parasite utilizes Plasmodium falciparum vacuolar pyrophosphatase 1 (PfVP1), an ancient pyrophosphate-driven proton pump, to pump protons across the parasite plasma membrane to maintain the membrane potential and cytosolic pH. Conditional deletion of PfVP1 leads to delayed ring stage development and a complete blockage of the ring-to-trophozoite transition, which can be partially rescued by Arabidopsis thaliana vacuolar pyrophosphatase 1, but not by the soluble pyrophosphatase from Saccharomyces cerevisiae. Proton-pumping pyrophosphatases are absent in humans, which highlights the possibility of developing highly selective VP1 inhibitors against the malaria parasite. Glucose inhibits haemostasis and accelerates diet-induced hyperlipidaemia in zebrafish larvae Hyperglycaemia damages the microvasculature in part through the reduced recruitment of immune cells and interference with platelet signalling, leading to poor wound healing and accelerated lipid deposition in mammals. We investigated the utility of zebrafish larvae to model the effect of glucose on neutrophil and macrophage recruitment to a tail wound, wound-induced haemostasis, and chicken egg yolk feed challenge-induced hyperlipidaemia by supplementing larvae with exogenous glucose by immersion or injection. Neither method of glucose supplementation affected the recruitment of neutrophils and macrophages following tail transection. Glucose injection reduced thrombocyte retention and fibrin plug formation while only thrombocyte retention was reduced by glucose immersion following tail transection. We observed accelerated lipid accumulation in glucose-injected larvae challenged with high fat chicken egg yolk feeding. Our study identifies conserved and divergent effects of high glucose on inflammation, haemostasis, and hyperlipidaemia in zebrafish larvae compared to mammals. "Adopt a Microorganism" methodology enhances the discourse richness of high school students in both hybrid learning and remote education The Internet has changed the way teachers and students access information and build knowledge. The recent COVID-19 pandemic has created challenges for both teachers and students and a demand for new methodologies of remote learning. In the life sciences, mixing online content with practical activities represents an even greater challenge. In microbiology, the implementation of an active teaching methodology, the #Adopt project, based on the social network Facebook(R), represents an excellent option for connecting remote education with classroom activities. In 2020, the version applied in high school, "Adopt a Microorganism", was adapted to meet the demands of emergency remote education owing to the suppression of face-to-face activities caused by the pandemic. In the present study, we assessed how the change in methodology impacted the discourse richness of students from high school integrated with technical education in the Business Administration program of the Federal Institute of Sao Paulo, Sorocaba Campus. Three questionnaires related to the groups of microorganisms (Archaea, Bacteria, Virus, Fungi, and Protozoan) were applied. The students responses in the 2019 and 2020 classes were compared concerning content richness and multiplicity of concepts through the application of the Shannon diversity index, an approach that is generally used to assess biodiversity in different environments. The observed results suggest that remote learning provided students with a conceptual basis and richness of content equivalent to that achieved by students subjected to the hybrid teaching model. In conclusion, this study suggests that the #Adopt project methodology increases students discourse richness in microbiology even without face-to-face traditional classes. Neuromodulators generate multiple context-relevant behaviors in a recurrent neural network by shifting activity hypertubes Mood, arousal, and other internal states can drastically alter behavior, even in identical external circumstances -- a cold glass of water when you are thirsty is much more desirable than when you are sated. Neuromodulators are critical controllers of such neural states, with dysfunctions linked to various neuropsychiatric disorders. Although biological aspects of neuromodulation have been well studied, the computational principles underlying how large-scale neuromodulation of distributed neural populations shifts brain states remain unclear. We use recurrent neural networks to model how synaptic weight modulation -- an important function of neuromodulators -- can achieve nuanced alterations in neural computation, even in a highly simplified form. We find that under structural constraints like those in brains, this provides a fundamental mechanism that can increase the computational capability and flexibility of a neural network by enabling overlapping storage of synaptic memories able to generate diverse, even diametrically opposed, behaviors. Our findings help explain how neuromodulators "unlock" specific behaviors by creating task-specific hypertubes in the space of neural activities and motivate more flexible, compact and capable machine learning architectures. Modular Dynamic Biomolecular Modelling with Bond Graphs: The Unification of Stoichiometry, Thermodynamics, Kinetics and Data. Renewed interest in dynamic simulation models of biomolecular systems has arisen from advances in genome-wide measurement and applications of such models in biotechnology and synthetic biology. In particular, genome-scale models of cellular metabolism beyond the steady state are required in order to represent transient and dynamic regulatory properties of the system. Development of such whole-cell models requires new modelling approaches. Here we propose the energy-based bond graph methodology, which integrates stoichiometric models with thermo-dynamic principles and kinetic modelling. We demonstrate how the bond graph approach intrinsically enforces thermodynamic constraints, provides a modular approach to modelling, and gives a basis for estimation of model parameters leading to dynamic models of biomolecular systems. The approach is illustrated using a well-established stoichiometric model of Escherichia coli (E. coli) and published experimental data.