Excitatory and inhibitory motor cortical clusters account for balance, variability, and task performance Behaving animals have to process partial or incomplete information to plan movements. Neural activity and behavior of highly trained animals are strikingly variable across behavioral trials. In the present study, we show that incomplete information is reflected in both neural and behavioral variability. We study a mechanism for such variability in spiking neural network models with cluster topologies that enables multistability and attractor dynamics. Multistable attractors have been studied in spiking neural networks through clusters of strongly interconnected excitatory neurons. However, we show that existing models suffer from weak robustness and fail to reproduce the experimental findings. To overcome these problems, we introduce a novel architecture with excitatory and inhibitory clusters of spiking neurons. In contrast to earlier prominent models of excitatory clusters, this model ensures the local balance of excitation and inhibition and shows robust multistability across a wide range of network parameters. Furthermore, with the appropriate stimulation of network clusters, this network topology enables qualitatively and quantitatively reproducing in vivo firing rates, variability dynamics and behavioral reaction times for different task conditions as observed in recordings from the motor cortex of behaving monkeys.
Pax3 repairs a neural circuit through a program of directed axon outgrowth Repairing damaged or dysfunctional human brain circuits remains an ongoing challenge for biomedical science. While surviving neuronal networks can be reorganised after lesion, for example by neurotrophins, these new connections are disorganised and rarely produce clinical improvement. Here we investigate how to promote axonal growth while retaining correct cellular targeting. We show that, in response to brain-derived neurotrophic factor (BDNF) in target-tissue, potential reinnervating neurons upregulate Pax3. Pax3 in turn increases polysialic acid-neural cell adhesion molecule (PSA-NCAM) on their axon terminals, facilitating their outgrowth and pathfinding, and resulting in correctly-targeted neural circuit repair in the mature nervous system. This is a novel role for Pax3, which we confirmed by showing its expression in afferent neurons is essential for spontaneous and BDNF-induced reinnervation in the developing and mature brains, respectively. Together these results suggest that Pax3 contributes to a repair program, in which axon growth is promoted and direction signaling maintained. These data advance our understanding for accurately rebuilding neural circuits: restricting growth-promotion to potential afferent neurons, as opposed to stimulating the whole circuit, allows axon growth without impairing its guidance.
Anti- SARS-CoV-2 Receptor Binding Domain Antibody Evolution after mRNA Vaccination Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection produces B-cell responses that continue to evolve for at least one year. During that time, memory B cells express increasingly broad and potent antibodies that are resistant to mutations found in variants of concern1. As a result, vaccination of coronavirus disease 2019 (COVID-19) convalescent individuals with currently available mRNA vaccines produces high levels of plasma neutralizing activity against all variants tested1, 2. Here, we examine memory B cell evolution 5 months after vaccination with either Moderna (mRNA-1273) or Pfizer- BioNTech (BNT162b2) mRNA vaccines in a cohort of SARS-CoV-2 naive individuals. Between prime and boost, memory B cells produce antibodies that evolve increased neutralizing activity, but there is no further increase in potency or breadth thereafter. Instead, memory B cells that emerge 5 months after vaccination of naive individuals express antibodies that are similar to those that dominate the initial response. While individual memory antibodies selected over time by natural infection have greater potency and breadth than antibodies elicited by vaccination, the overall neutralizing potency of plasma is greater following vaccination. These results suggest that boosting vaccinated individuals with currently available mRNA vaccines will increase plasma neutralizing activity but may not produce antibodies with breadth equivalent to those obtained by vaccinating convalescent individuals. Characterization of methylation profiles in spontaneous preterm birth placental villous tissue Preterm birth is a global public health crisis which results in significant neonatal and maternal mortality. Yet little is known regarding the molecular mechanisms of idiopathic spontaneous preterm birth, and we have few diagnostic markers for adequate assessment of placental development and function. Previous studies of placental pathology and our transcriptomics studies suggest a role for placental maturity in idiopathic spontaneous preterm birth. It is known that placental DNA methylation changes over gestation. We hypothesized that if placental hypermaturity is present in our samples, we would observe a unique idiopathic spontaneous preterm birth DNA methylation profile potentially driving the gene expression differences we previously identified in our placental samples. Our results indicate the idiopathic spontaneous preterm birth DNA methylation pattern mimics the term birth methylation pattern suggesting hypermaturity. Only seven significant differentially methylated regions fitting the idiopathic spontaneous preterm birth specific (relative to the controls) profile were identified, indicating unusually high similarity in DNA methylation between idiopathic spontaneous preterm birth and term birth samples. We identified an additional 1,718 significantly methylated regions in our gestational age matched controls were the idiopathic spontaneous preterm birth DNA methylation pattern mimics the term birth methylation pattern, again indicating a striking level of similarity between the idiopathic spontaneous preterm birth and term birth samples. Pathway analysis of these regions revealed differences in genes within the WNT and Cadherin signaling pathways, both of which are essential in placental development and maturation. Taken together, these data demonstrate that the idiopathic spontaneous preterm birth samples are molecularly more mature than expected given their respective gestational age which likely impacts birth timing. Reconstitution of minimal motility system based on Spiroplasma swimming by expressing two bacterial actins in synthetic minimal bacterium Motility is one of the most important features of life, but its evolutionary origin remains unknown. In this study, we focused on Spiroplasma, commensal, or parasitic bacteria. They swim by switching the helicity of a ribbon-like cytoskeleton that comprises six proteins, each of which evolved from a nucleosidase and bacterial actin called MreB. We expressed these proteins in a synthetic, non-motile minimal bacterium, JCVI-syn3.0B, whose reduced genome was computer-designed and chemically synthesized. The synthetic bacterium exhibited swimming motility with features characteristic of Spiroplasma swimming. Moreover, some combinations of the two proteins produced a helical cell shape and swimming, suggesting that the swimming originated from the differentiation and coupling of bacterial actins, and we also obtained a minimal system for motility of the synthetic bacterium.

One-Sentence SummaryThe minimal system comprised two bacterial actins that provided cell helicity and swimming to the synthetic minimal bacterium. Spatial organization of chromosomes leads to heterogeneous chromatin motion and drives the liquid- or gel-like dynamical behavior of chromatin Chromosome organization and dynamics are involved in regulating many fundamental processes such as gene transcription and DNA repair. Experiments unveiled that chromatin motion is highly heterogeneous inside cell nuclei, ranging from a liquid-like, mobile state to a gel-like, rigid regime. Using polymer modeling, we investigate how these different physical states and dynamical heterogeneities may emerge from the same structural mechanisms. We found that the formation of topologically-associating domains (TADs) is a key driver of chromatin motion heterogeneity. In particular, we demonstrated that the local degree of compaction of the TAD regulates the transition from a weakly compact, fluid state of chromatin to a more compact, gel state exhibiting anomalous diffusion and coherent motion. Our work provides a comprehensive study of chromosome dynamics and a unified view of chromatin motion enabling to interpret the wide variety of dynamical behaviors observed experimentally across different biological conditions, suggesting that the  liquid or  solid behaviour of chromatin are in fact two sides of the same coin. Engineered disulfide reveals structural dynamics of locked SARS-CoV-2 spike The spike (S) protein of SARS-CoV-2 has been observed in three distinct pre-fusion conformations: locked, closed and open. Of these, the function of the locked conformation remains poorly understood. Here we engineered a SARS-CoV-2 S protein construct "S-R/x3" to arrest SARS-CoV-2 spikes in the locked conformation by a disulfide bond. Using this construct we determined high-resolution structures confirming that the x3 disulfide bond has the ability to stabilize the otherwise transient locked conformations. Structural analyses reveal that wild-type SARS-CoV-2 spike can adopt two distinct locked-1 and locked-2 conformations. For the D614G spike, based on which all variants of concern were evolved, only the locked-2 conformation was observed. Analysis of the structures suggests that rigidified domain D in the locked conformations interacts with the hinge to domain C and thereby restrains RBD movement. Structural change in domain D correlates with spike conformational change. We propose that the locked-1 and locked-2 conformations of S are present in the acidic high-lipid cellular compartments during virus assembly and egress. In this model, release of the virion into the neutral pH extracellular space would favour transition to the closed or open conformations. The dynamics of this transition can be altered by mutations that modulate domain D structure, as is the case for the D614G mutation, leading to changes in viral fitness. The S-R/x3 construct provides a tool for the further structural and functional characterization of the locked conformations of S, as well as how sequence changes might alter S assembly and regulation of receptor binding domain dynamics. Differences in gait stability and acceleration characteristics between healthy young and older females Our aim was to evaluate differences in gait acceleration intensity, variability and stability of feet and trunk between older females and young females using inertial sensors. Twenty older females (OF; mean age 68.4, SD 4.1 years) and eighteen young females (YF; mean age 22.3, SD 1.7 years) were asked to walk straight for 100 meters at their preferred speed, while wearing inertial sensors on heels and lower back. We calculated spatiotemporal measures, foot and trunk acceleration characteristics and their variability, as well as trunk stability using the local divergence exponent (LDE). Two-way analysis of variance (including the factors foot and age), Students t-test, and Mann-Whitney U test were used to compare statistical differences of measures between groups. Cohens d effects were calculated for each variable. Foot maximum vertical acceleration and amplitude, trunk-foot vertical acceleration attenuation, as well as their variability were significantly smaller in OF than in YF. In contrast, trunk mediolateral acceleration amplitude, maximum vertical acceleration, and amplitude, as well as their variability were significantly larger in OF than in YF. Moreover, OF showed lower stability (i.e., higher LDE values) in mediolateral acceleration, mediolateral and vertical angular velocity of the trunk. Even though we measured healthy older females, these participants showed lower vertical foot accelerations with higher vertical trunk acceleration, lower trunk-foot vertical acceleration attenuation, less gait stability, and more variability of the trunk, and hence, were more likely to fall. These findings suggest that instrumented gait measurements may help for early detection of changes or impairments in gait performance, even before this can be observed by clinical eye or gait speed. The initial response of females towards congeneric males matches the propensity to hybridise in Ophthalmotilapia. AO_SCPLOWBSTRACTC_SCPLOWCichlid radiations often harbour closely related species with overlapping niches and distribution ranges. Such species sometimes hybridise in nature, which raises the question how can they coexist. This also holds for the Tanganyika mouthbrooders Ophthalmotilapia ventralis and O. nasuta. Earlier studies found indications of asymmetrical hybridisation with females of O. ventralis accepting males of O. nasuta, but not the other way around. We hypothesised that this was due to differences in the capacity for species recognition. Given the higher propensity of O. ventralis females towards hybridisation, we expect a reduced ability for species recognition in O. ventralis females, compared to O. nasuta females. We staged two experiments, one focusing on 22 female O. nasuta and one on 21 female O. ventralis. These fish were placed in one half of a tank and briefly exposed to a conspecific or a heterospecific male, a conspecific female, or nothing (control). Female response was evaluated by scoring six tracking parameters and by noting the occurrence of ten discrete behaviours before and during the encounter. Females always responded to the presence of another fish by approaching it. Remarkably, for both O. nasuta and O. ventralis, we did not find a different response between encounters with conspecific males and females. However, in agreement with our hypothesis, O. nasuta females behaved differently towards conspecific or heterospecific males, whereas O. ventralis females did not. When presented with a heterospecific male, O. nasuta females performed a lower number of  ram behaviours. Additionally, they never displayed the  flee behaviour, a component of the species mating repertoire that was seen in all but one of the presentations with a conspecific male. Our findings show that differences in species recognition at first encounter predict to a large degree the outcome of the mating process, even in the absence of mating behaviour. Resolving when (and where) the Thylacine went extinct Like the Dodo and Passenger Pigeon before it, the predatory marsupial Thylacine (Thylacinus cynocephalus), or  Tasmanian tiger, has become an iconic symbol of human-caused extinction. The last captive animal died in 1936, but even today reports of the Thylacines possible ongoing survival in remote regions of Tasmania are newsworthy and capture the publics imagination. Extirpated from mainland Australia in the mid-Holocene, the large island of Tasmania became the species final stronghold. Following European settlement in the 1800s, the Thylacine was heavily persecuted and pushed to the margins of its range, although many sightings were reported thereafter--even well beyond the 1930s. To gain a new depth of insight into the extinction of the Thylacine, we assembled an exhaustive database of 1,237 observational records from Tasmania (from 1910 onwards), quantified their uncertainty, and charted the patterns these revealed. We also developed a new method to visualize the species 20th-century spatio-temporal dynamics, to map potential post-bounty refugia and pinpoint the most-likely location of the final persisting subpopulation. A direct reading of the high-quality records (confirmed kills and captures, in combination with sightings by past Thylacine hunters and trappers, wildlife professionals and experienced bushmen) implies a most-likely extinction date within four decades following the last capture (i.e., 1940s to 1970s). However, uncertainty modelling of the entire sighting record, where each observation is assigned a probability and the whole dataset is then subject to a sensitivity analysis, suggests that extinction might have been as recent as the late 1980s to early 2000s, with a small chance of persistence in the remote south-western wilderness areas. Beyond the intrinsically fascinating problem of reconstructing the final fate of the Thylacine, the new spatio-temporal mapping of extirpation developed herein would also be useful for conservation prioritization and search efforts for other rare taxa of uncertain status. A parametric logistic equation with light flux and medium concentration for cultivation planning of microalgae Microalgae are considered to be promising producers of bioactive chemicals, feeds, and fuels from carbon dioxide by photosynthesis. Thus, the prediction of microalgal growth profiles is important for the planning of cost-effective and sustainable cultivation-harvest cycles. This paper proposes a mathematical model capable of predicting the effect of light flux into culture and medium concentration on the growth profiles of microalgae by incorporating these growth-limiting factors into the logistic equation. The specific form of the equation is derived based on the experimentally measured growth profiles of Monoraphidium sp., a microalgal strain isolated by the authors, under 16 conditions consisting of the combinations of incident light fluxes into culture and initial medium concentrations. Using a cross-validation method, it is shown that the proposed model has the ability to predict necessary incident light flux into culture and initial medium concentration for harvesting target biomass at target time. Finally, the model-guided cultivation planning is performed and is evaluated by comparing the result with the experimental data. Host species identity shapes the diversity and structure of insect microbiota As for most of the life that inhabits our planet, microorganisms play an essential role in the fitness of insects, including nutrition, reproduction, defence, and many other functions. More recently, we assisted in an exponential growth of studies describing the taxonomical composition of bacterial communities across insects phylogeny. However, there is still an outstanding question that needs to be answered: which factors contribute most to shape insects microbiomes? This study tries to find an answer to this question by taking advantage of publicly available sequencing data and reanalysing over 4,000 samples of insect-associated bacterial communities under a common framework. Results suggest that insect taxonomy has a wider impact on the structure and diversity of their associated microbial communities than the other factors considered (diet, sex, life stage, sample origin, and treatment). However, when specifically testing for signatures of co-diversification of insect species and their microbiota, analyses found weak support for this, suggesting that while insect species strongly drive the structure and diversity of insect microbiota, the diversification of those microbial communities did not follow their hosts phylogeny. Furthermore, a parallel survey of the literature highlights several methodological limitations that need to be considered in future research endeavours. Decoding the porcine developing spatial processing system and production of human entorhinal stellate cell-like cells by a direct programming approach. Classic studies investigating how and when the entorhinal cortex (component of the memory processing system of the brain) develops have been based on traditional thymidine autoradiography and histological techniques. In this study, we take advantage of modern technologies to trace at a high resolution, the cellular complexity of the developing porcine medial entorhinal cortex by using single-cell profiling. The postnatal medial entorhinal cortex comprises 4 interneuron, 3 pyramidal neuron and 2 stellate cell populations which emerge from intermediate progenitor and immature neuron populations. We discover four MGE-derived interneurons and one CGE-derived interneuron population as well as several IN progenitors. We also identify two oligodendrocyte progenitor populations and three populations of oligodendrocytes. We perform a proof-of-concept experiment demonstrating that porcine scRNA-seq data can be used to develop novel protocols for producing human entorhinal cells in-vitro. We identified six transcription factors (RUNX1A1, SOX5, FOXP1, MEF2C, TCF4, EYA2) important in neurodevelopment and differentiation from one RELN+ stellate cell population. Using a lentiviral vector approach, we reprogrammed human induced pluripotent stem cells into stellate cell-like cells which expressed RELN, SATB2, LEF1 and BCL11B. Our findings contribute to the understanding of the formation of the brains cognitive memory and spatial processing system and provides proof-of-concept for the production of entorhinal cells from human pluripotent stem cells in-vitro. 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. Molecular diversity maintained by long-term balancing selection in mating loci defines multiple mating types in fungi Balancing selection, an evolutionary force that retains genetic diversity, has been detected in multiple genes and organisms, such as the sexual mating loci in fungi. However, to quantify the strength of balancing selection and define the mating-related genes require a large number of specimens. In tetrapolar basidiomycete fungi, sexual type is determined by two unlinked loci, MATA and MATB. Genes in both loci defines mating type identity, control successful mating and completion of the life cycle. These loci are usually highly diverse. Previous studies have speculated, based on culture crosses, that species of the non-model genus Trichaptum (Hymenochaetales, Basidiomycota) possess a tetrapolar mating system, with multiple alleles. Here, we sequenced a hundred and eighty specimens of three Trichaptum species. We characterized the chromosomal location of MATA and MATB, the molecular structure of MAT regions and their allelic richness. Our sequencing effort was sufficient to molecularly characterize multiple MAT alleles segregating before the speciation event of Trichaptum species. Our analyses suggested that long-term balancing selection has generated trans-species polymorphisms. Mating sequences were classified in different allelic classes based on an amino acid identity (AAI) threshold supported by phylogenetics. The inferred allelic information mirrored the outcome of in vitro crosses, thus allowing us to support the degree of allelic divergence needed for successful mating. Even with the high amount of divergence, key amino acids in functional domains are conserved. The observed allelic classes could potentially generate 14,560 different mating types. We conclude that the genetic diversity of mating in Trichaptum loci is due to long-term balancing selection, with limited recombination and duplication activity. Our large number of sequenced specimens highlighted the importance of sequencing multiple individuals from different species to detect the mating-related genes, the mechanisms generating diversity and the evolutionary forces maintaining them.

Author summaryFungi have complex mating systems, and basidiomycete fungi can encode thousands of mating types. Individuals with the same mating type cannot mate. This sexual system has evolved to facilitate sexual mating, increasing the chances to recombine into advantageous allelic combination and prune deleterious alleles. We explored the genomes of hundred and eighty specimens, combined with experimental mating studies of selected specimens, from a non-model organism (Trichaptum). We characterized the genomic regions controlling sex. The mating ability of the specimens confirmed the role of the mating alleles observed in the genomic data. The detailed analyses of many specimens allowed us to observe gene duplication and rearrangements within the mating loci, increasing the diversity within these loci. We supported previous suggestions of balancing selection in this region, an evolutionary force that maintains genomic diversity. These results supports that our fungal specimens are prone to outcross, which might facilitate the adaptation to new conditions. Targeted Gene Editing in Plants using CRISPR-Cas9, Single-Stranded DNA Oligonucleotides, and Protoplast Regeneration Versatile genome editing can be facilitated by the insertion of DNA sequences into specific locations. Current protocols involving clustered regularly interspaced short palindromic repeats (CRISPR) and CRISPR-associated (Cas) proteins rely on low efficiency homology-directed repair or non-homologous end joining with modified double-stranded DNA oligonucleotides as donors. Our simple protocol eliminates the need for expensive equipment, chemical and enzymatic donor DNA modification or plasmid construction by using polyethylene glycol-calcium to deliver non-modified single-stranded DNA oligonucleotides and CRISPR-Cas9 ribonucleoprotein into protoplasts. Plants regenerated via edited protoplasts achieved targeted insertion frequencies of up to 50.0% in Nicotiana benthamiana and 13.6% in rapid cycling Brassica oleracea without antibiotic selection. Using a 60-nt donor containing 27 nt in each homologous arm, 6 of 22 regenerated N. benthamiana plants showed targeted insertions, and one contained a precise insertion of a 6-bp HindIII site. The inserted sequences were transmitted to the next generation and invite the possibility of future exploration of versatile genome editing by targeted DNA insertion in plants. Erythropoietin supplementation ameliorates the immunological and hematological deficiencies of lysinuric protein intolerance in mice SLC7A7 encodes for y+LAT1, a cationic amino acid transporter expressed on the basolateral membrane of epithelial cells and macrophages. Mutations in human SLC7A7 cause lysinuric protein intolerance (LPI), a rare autosomal disorder with a variety of clinical symptoms. Citrulline and a low protein diet are widely used to treat the main symptoms of LPI; however, the hematological and immunological abnormalities in LPI remain poorly understood and difficult to manage by established treatments. Here, we show that erythropoietin levels are significantly diminished in mice lacking Slc7a7, which reduces the number of erythroblast precursors, resulting in mature red blood cells with a smaller mean corpuscular volume. Ex vivo analysis demonstrated that erythrocytes from Slc7a7-deficient mice were phagocytosed more by bone marrow-derived macrophages. Gene expression analysis of sorted red pulp macrophages (RPMs) from Slc7a7-deficient mice indicated that this cell population had higher erythrophagocytic activity, but preserved heme catabolic activity. The increase in erythrophagocytosis resulted in iron accumulation in macrophages and tissues, and ultimately led to RPM apoptosis. Surprisingly, ablation of Slc7a7 in myeloid cells failed to trigger tissue iron accumulation or changes in the RPM population, indicating that the expression of this gene in macrophages is not key for the development of iron accumulation in tissue. Finally, restoration of serum erythropoietin levels in Slc7a7-deficient mice re-established normal erythropoiesis and erythrophagocytosis and prevented the accumulation of iron in tissues. Thus, erythropoietin has a fundamental role in the immune-hematological alterations of LPI in mice, revealing a new crucial role for this hormone in this disease. 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. 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. 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. Phenotypic plasticity is aligned with phenological adaptation on both micro- and macroevolutionary timescales In seasonally-variable environments, phenotypic plasticity in phenology may be critical for adaptation to fluctuating environmental (temperature) conditions. Using an 18-generation longitudinal dataset from natural insect (damselfly) populations, we show that phenology has strongly advanced. Individual fitness data suggest this is likely an adaptive response towards a thermally-dependent fitness optimum. A laboratory experiment revealed that developmental plasticity qualitatively matches environmental-dependence of selection, partially explaining observed phenological advance. Expanding our analysis to the macroevolutionary level, we use a database of over 1-million occurrence records and spatiotemporally-matched temperature data from 49 Swedish Odonate species to infer macroevolutionary dynamics of phenology. Phenological plasticity was more closely aligned with adaptation for species that have recently colonized northern latitudes, with more phenological mismatch at lower latitudes. Our results show that phenological plasticity plays a key role in microevolutionary dynamics within a single species, and such plasticity may have facilitated post-Pleistocene range expansion in this insect clade. Quantifying ERK-activity in response to inhibition of the BRAFV600E-MEK-ERK cascade using mathematical modelling BackgroundSimultaneous inhibition of multiple components of the BRAF-MEK-ERK cascade (vertical inhibition) has become a standard of care for treating BRAF-mutant melanoma. However, the molecular mechanism of how vertical inhibition synergistically suppresses intracellular ERK-activity, and consequently cell proliferation, are yet to be fully elucidated.

MethodsWe develop a mechanistic mathematical model that describes how the mutant BRAF-inhibitor, dabrafenib, and the MEK-inhibitor, trametinib, affect BRAFV600E-MEK-ERK signalling. The model is based upon a system of chemical reactions that describes cascade signalling dynamics. Using mass action kinetics, the chemical reactions are re-expressed as ordinary differential equations which are parameterised by in vitro data and solved numerically to obtain the temporal evolution of cascade component concentrations.

ResultsThe model provides a quantitative method to compute how dabrafenib and trametinib can be used in combination to synergistically inhibit ERK-activity in BRAFV600E-mutant melanoma cells. The model elucidates molecular mechanisms of vertical inhibition of the BRAFV600E-MEK-ERK cascade and delineates how elevated BRAF concentrations generate drug resistance to dabrafenib and trametinib. The computational simulations further suggest that elevated ATP levels could be a factor in drug resistance to dabrafenib.

ConclusionsThe model can be used to systematically motivate which dabrafenib-trametinib dose-combinations, for treating BRAFV600E-mutated melanoma, warrant experimental investigation. Genomic regions associate with major axes of variation driven by gas exchange and leaf construction traits in cultivated sunflower (Helianthus annuus L.) Stomata and leaf veins play an essential role in transpiration and the movement of water throughout leaves. These traits are thus thought to play a key role in the adaptation of plants to drought and a better understanding of the genetic basis of their variation and coordination could inform efforts to improve drought tolerance. Here, we explore patterns of variation and covariation in leaf anatomical traits and analyze their genetic architecture via genome-wide association (GWA) analyses in cultivated sunflower (Helianthus annuus L.). Traits related to stomatal density and morphology as well as lower order veins were manually measured from digital images while the density of minor veins was estimated using a novel deep learning approach. Leaf, stomatal, and vein traits exhibited numerous significant correlations that generally followed expectations based on functional relationships. Correlated suites of traits could further be separated along three major principal component (PC) axes that were heavily influenced by variation in traits related to gas exchange, leaf hydraulics, and leaf construction. While there was limited evidence of colocalization when individual traits were subjected to GWA analyses, major multivariate PC axes that were most strongly influenced by several traits related to gas exchange or leaf construction did exhibit significant genomic associations. These results provide insight into the genetic basis of leaf trait covariation and showcase potential targets for future efforts aimed at modifying leaf anatomical traits in sunflower.

Significance StatementUsing traditional and automated/high-throughput (using a novel deep learning approach) phenotyping methods we studied leaf anatomical variation in sunflower. Genome-wide association (GWA) analyses identified numerous genomic regions underlying individual trait variation and regions underlying major multivariate axes of phenotypic variation. These results illustrate the value of employing a multivariate approach to GWA analyses and shed light on the extent to which leaf trait (co-)variation can be genetically decoupled to explore novel phenotypic space. Cell cycle-independent integration of stress signals promotes Non-G1/G0 quiescence entry Cellular quiescence is a non-proliferative state required for cell survival under stress and during development. In most quiescent cells, proliferation is stopped in a reversible state of low Cdk1 kinase activity; in many organisms, however, quiescent states with high Cdk1 activity can also be established through still uncharacterized stress or developmental mechanisms. Here, we used a microfluidics approach coupled to phenotypic classification by machine learning to identify stress pathways associated with starvation-triggered high-Cdk1 quiescent states in Saccharomyces cerevisiae. We found that low- and high-Cdk1 quiescent states shared a core of stress-associated processes, such as autophagy, protein aggregation, and mitochondrial upregulation, but differed in the nuclear accumulation of the stress transcription factors Xbp1, Gln3, and Sfp1. The decision between low- or high-Cdk1 quiescence was controlled by cell cycle-independent accumulation of Xbp1, which acted as a time-delayed integrator of the duration of stress stimuli. Our results show how cell cycle-independent stress-activated factors promote cellular quiescence outside of G1/G0. Group-walk, a rigorous approach to group-wise false discovery rate analysis by target-decoy competition Target-decoy competition (TDC) is a commonly used method for false discovery rate (FDR) control in the analysis of tandem mass spectrometry data. This type of competitionbased FDR control has recently gained significant popularity in other fields after Barber and Candes laid its theoretical foundation in a more general setting that included the feature selection problem. In both cases, the competition is based on a head-to-head comparison between an (observed) target score and a corresponding decoy (knockoff) score. However, the effectiveness of TDC depends on whether the data is homogeneous, which is often not the case: in many settings, the data consists of groups with different score profiles or different proportions of true nulls. In such cases, applying TDC while ignoring the group structure often yields imbalanced lists of discoveries, where some groups might include relatively many false discoveries and other groups include relatively very few. On the other hand, as we show, the alternative approach of applying TDC separately to each group does not rigorously control the FDR.

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

Group-walk is available at https://cran.r-project.org/web/packages/groupwalk/index.html Native proline-rich motifs exploit sequence context to target actin-remodeling Ena/VASP proteins The human proteome is replete with short linear motifs (SLiMs) of 4-6 residues that are critical for protein-protein interactions, yet the importance of the sequence surrounding such motifs is underexplored. We devised a proteomic screen to systematically examine the influence of SLiM sequence context on protein-protein interactions. Focusing on the EVH1 domain of ENAH, an actin regulator that is upregulated in invasive cancers, we screened 36-residue proteome-derived peptides for binding. We discovered a pocket on the ENAH EVH1 domain that diverged from its orthologs to recognize extended SLiMs, and we found that proteins with two EVH1-binding SLiMs can wrap around a single domain. We also found that the ciliary protein PCARE uses an extended 23-residue region to obtain higher affinity than any known ENAH EVH1-binding motif. Our screen provides a way to uncover the effects of broader proteomic context on motif-mediated interactions, revealing diverse mechanisms of contextual control over EVH1 interactions and establishing that SLiMs cant be fully understood outside of their native context.