Patterns of connectome variability in autism across five functional activation tasks. Findings from the LEAP project BackgroundAutism spectrum disorder (autism) is a complex neurodevelopmental condition with pronounced behavioural, cognitive, and neural heterogeneities across individuals. Here, our goal was to characterise heterogeneity in autism by identifying patterns of neural diversity as reflected in BOLD fMRI in the way individuals with autism engage with a varied array of cognitive tasks.

MethodsAll analyses were based on the EU-AIMS/AIMS-2-TRIALS multisite Longitudinal European Autism Project (LEAP) with participants with autism and typically developing controls (TD) between 6 and 30 years of age. We employed a novel task-potency approach which combines the unique aspects of both resting-state fMRI and task-fMRI to quantify task-induced variations in the functional connectome. Normative modelling was used to map atypicality of features on an individual basis with respect to their distribution in neurotypical control participants. We applied robust out-of-sample canonical correlation analysis (CCA) to relate connectome data to behavioural data.

ResultsDeviation from the normative ranges of global functional connectivity was greater for individuals with autism compared to TD in each fMRI task paradigm (all tasks p<0.001). The similarity across individuals of the deviation pattern was significantly increased in autistic relative to TD individuals (p<0.002). The CCA identified significant and robust brainbehavior covariation between functional connectivity atypicality and autism-related behavioral features.

ConclusionsIndividuals with autism engage with tasks in a globally atypical way, but the particular spatial pattern of this atypicality is nevertheless similar across tasks. Atypicalities in the tasks originate mostly from prefrontal cortex and default mode network regions, but also speech and auditory networks. We show, moving forward, sophisticated modeling methods such as task-potency and normative modeling will prove key to unravelling complex heterogeneous conditions like autism.
Motor planning brings human primary somatosensory cortex into action-specific preparatory states Motor planning plays a critical role in producing fast and accurate movement. Yet, the neural processes that occur in human primary motor and somatosensory cortex during planning, and how they relate to those during movement execution, remain poorly understood. Here we used 7T functional magnetic resonance imaging (fMRI) and a delayed movement paradigm to study single finger movement planning and execution. The inclusion of no-go trials and variable delays allowed us to separate what are typically overlapping planning and execution brain responses. Although our univariate results show widespread deactivation during finger planning, multivariate pattern analysis revealed finger-specific activity patterns in contralateral primary somatosensory cortex (S1), which predicted the planned finger action. Surprisingly, these activity patterns were as informative as those found in contralateral primary motor cortex (M1). Control analyses ruled out the possibility that the detected information was an artifact of subthreshold movements during the preparatory delay. Furthermore, we observed that finger-specific activity patterns during planning were highly correlated to those during execution. These findings reveal that motor planning activates the specific S1 and M1 circuits that are engaged during the execution of a finger press, while activity in both regions is overall suppressed. We propose that preparatory states in S1 may improve movement control through changes in sensory processing or via direct influence of spinal motor neurons.

Significance statementMotor planning is important for good behavioral performance, yet it is unclear which neural processes underlie the preparation of the nervous system for an upcoming movement. Using high-resolution functional neuroimaging, we investigated how motor planning for finger presses changes the activity state in primary motor and primary somatosensory cortex, and how brain responses during planning and execution relate to each other. We show that planning leads to finger-specific activation in both M1 and S1, which is highly similar to the finger-specific activity patterns elicited during execution. Our findings suggest that S1 is being specifically prepared for an upcoming action, either to actively contribute to the outflowing motor command or to enable action-specific sensory gating.
Ultrahigh-Throughput Ambient MS: Direct Analysis at 22 Samples per Second by Infrared Matrix-Assisted Laser Desorption Electrospray Ionization Mass Spectrometry Infrared Matrix-Assisted Laser Desorption Electrospray Ionization (IR-MALDESI) mass spectrometry is an ambient-direct sampling method being developed for high-throughput, label-free, biochemical screening of large-scale compound libraries. Here, we report the development of an ultrahigh-throughput continuous motion IR-MALDESI sampling approach capable of acquiring data at rates up to 22.7 samples per second in a 384-well microtiter plate. At top speed, less than 1% analyte carryover is observed from well-to-well and signal intensity relative standard deviations (RSD) of 11.5% and 20.9% for 3 M 1-hydroxymidazolam and 12 M dextrorphan, respectively, are achieved. The ability to perform parallel kinetics studies on 384 samples with ~30s time resolution using an isocitrate dehydrogenase 1 (IDH1) enzyme assay is shown. Finally, we demonstrate the repeatability and throughput of our approach by measuring 115,200 samples from 300 microtiter plate reads consecutively over 5.54 hours with RSDs under 8.14% for each freshly introduced plate. Taken together, these results demonstrate the use of IR-MALDESI at sample acquisition rates that surpass other currently reported direct sampling mass spectrometry approaches used for high throughput compound screening.

For Table of Contents Only

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=125 SRC="FIGDIR/small/465730v3_ufig1.gif" ALT="Figure 1">
View larger version (49K):
org.highwire.dtl.DTLVardef@5f916corg.highwire.dtl.DTLVardef@4473b4org.highwire.dtl.DTLVardef@82e8forg.highwire.dtl.DTLVardef@128b57b_HPS_FORMAT_FIGEXP  M_FIG C_FIG 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. GPI-anchored FGF directs cytoneme-mediated bidirectional signaling to self-regulate tissue-specific dispersion How signaling proteins generate a multitude of information to organize tissue patterns is critical to understanding morphogenesis. In Drosophila, FGF produced in wing-disc cells regulates the development of the disc-associated air-sac-primordium (ASP). Here, we show that FGF is Glycosylphosphatidylinositol-anchored to the producing cell surface and that this modification both inhibits free FGF secretion and activates target-specific bidirectional FGF-FGFR signaling through cytonemes. FGF-source and recipient ASP cells extend cytonemes that present FGF and FGFR on their surfaces and reciprocally recognize each other over distance by contacting through CAM-like FGF-FGFR binding. Contact-mediated FGF-FGFR binding induces bidirectional signaling, which, in turn, promotes ASP and source cells to polarize cytonemes toward each other and reinforce signaling contacts. Subsequent un-anchoring of FGFR-bound-FGF from the source cell membrane dissociates cytoneme contacts and delivers FGF target-specifically to ASP cytonemes for paracrine functions. Thus, GPI-anchored FGF organizes both source and recipient cells and self-regulates its cytoneme-mediated tissue-specific dispersion and signaling. Juvenile Gouldian Finches (Erythrura gouldiae) form sibling sub-groups during social integration The formation of social relationships in complex groups is critical in shaping patterns of social organization and behavioral development. In many birds, young individuals remain dependent on their parents for extended periods but must abruptly transition to navigating interactions in the wider group after independence. While lack of social relationships during this period is detrimental in the development of later social skills, little is known about the social relationships juveniles form after independence in many bird species. In this study we describe patterns of social interactions in juvenile Gouldian Finches (Erythrura gouldiae) during transitions from family groups into flocks of unfamiliar individuals. A group of 20 juveniles from 4 families were introduced into two flocks. After introductions juveniles showed a gradient of approach rates with most approaches directed towards siblings, followed by juvenile peers, adult females, and lastly adult males. Significant preferences for siblings resulted in the emergence of sibling subgroups within the larger social network. This active self-assortment of siblings suggests that sibling sub-groups are an important bridge linking social connections within the family and the wider group. Such findings suggests that sibling relationships have a critical role in the socialization after independence, as well as structuring the social organization of Gouldian finch flocks. Nej1XLF interacts with Sae2CTIP to inhibit Dna2 mediated resection at DNA doublestrand break The two major pathways of DNA double strand break (DSB) repair, non-homologous end-joining (NHEJ) and homologous recombination (HR), are highly conserved from yeast to mammals. The regulation of 5 DNA resection controls repair pathway choice and influences repair outcomes. Nej1 was first identified as a canonical NHEJ factor involved in stimulating the ligation of broken DNA ends and more recently, it was shown to be important for DNA end-bridging and inhibiting 5 resection mediated by Dna2-Sgs1. Nej1 interacts with Sae2 and this impacts DSB repair in three ways. First, Nej1 inhibits MRX-Sae2 interactions and Sae2 localization to a DSB. Second, Nej1 inhibits Sae2-dependent recruitment of Dna2 in the absence of Sgs1. Third, NEJ1 and SAE2 showed an epistatic relationship for DNA end-bridging, an event that restrains the broken ends and reduces the frequency of genomic deletions from developing at the DSB. Deletion of NEJ1 suppressed the synthetic lethality of sae2{Delta} sgs1{Delta} and was dependent on the nuclease activity of Dna2. These Nej1 functions promote end-joining DSB repair, but could also be relevant for controlling resection initiation during HR repair.

HighlightsO_LINej1 physically interacts with Sae2 and inhibits end-resection at a DSB.
C_LIO_LINej1 inhibits Sae2 interactions with the MRX complex.
C_LIO_LINej1 inhibits Sae2-dependent recruitment of Dna2 to a DSB.
C_LIO_LINEJ1 and SAE2 are epistatic for DNA end-bridging.
C_LIO_LIDeletion of NEJ1 suppresses the synthetic lethality of sae2{Delta} sgs1{Delta}, which is dependent on Dna2 nuclease activity.
C_LI 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. The Tabula Sapiens: a multiple organ single cell transcriptomic atlas of humans Molecular characterization of cell types using single cell transcriptome sequencing is revolutionizing cell biology and enabling new insights into the physiology of human organs. We created a human reference atlas comprising nearly 500,000 cells from 24 different tissues and organs, many from the same donor. This atlas enabled molecular characterization of more than 400 cell types, their distribution across tissues and tissue specific variation in gene expression. Using multiple tissues from a single donor enabled identification of the clonal distribution of T cells between tissues, the tissue specific mutation rate in B cells, and analysis of the cell cycle state and proliferative potential of shared cell types across tissues. Cell type specific RNA splicing was discovered and analyzed across tissues within an individual. Multiscale mechanistic modelling of the host defense in invasiveaspergillosis reveals leukocyte activation and iron acquisition as driversof infection outcome Aspergillus species are ubiquitous environmental molds, the spores of which are inhaled daily by most humans. Immunocompromised hosts can develop an invasive infection resulting in high mortality. There is, therefore, a pressing need for host-centric therapeutics for this infection. To address this need, we created a multi-scale computational model of the infection, focused on its interaction with the innate immune system and iron, a critical nutrient for the pathogen. The model, parameterized using published data, was found to recapitulate a wide range of biological features and was experimentally validated in vivo. Conidial swelling was identified as critical in fungal strains with high growth, whereas the siderophore secretion rate is an essential prerequisite for establishment of the infection in low-growth strains. In immunocompetent hosts, high growth, high swelling probability, and impaired leukocyte activation lead to a high conidial germination rate. Similarly, in neutropenic mice, high fungal growth was achieved through synergy between high growth rate, high swelling probability, slow leukocyte activation, and high siderophore secretion. In summary, the model reveals a small set of parameters related to fungal growth, iron acquisition, and leukocyte activation as key determinants of the fate of the infection. Computation of the mitochondrial age distribution along the axon length We describe a compartmental model of mitochondrial transport in axons, which we apply to compute mitochondrial age at different distances from the soma. The model predicts that the average mitochondrial age at the tip of an axon with a length of 1 cm is approximately 22 hours. The mitochondria are youngest closest to the soma and their age scales approximately linearly with distance from the soma. To the best of the authors knowledge, this is the first attempt to predict the spatial distribution of mitochondrial age within an axon. A sensitivity study of the mean age of mitochondria to various model parameters is also presented. Allopurinol Blocks The Formation And Progression Of Aortic Aneurysm In A Mouse Model Of Marfan Syndrome 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 TCF11 Has a Potent Tumor-Repressing Effect than Its Prototypic Nrf1α by Definition of both Similar yet Different Regulatory Profiles, with a Striking Disparity from Nrf2 Nrf1 and Nrf2, as two principal CNC-bZIP transcription factors, regulate similar but different targets involved in a variety of biological functions for maintaining cell homeostasis and organ integrity. Of note, the unique topobiological behavior of Nrf1 makes its functions more complicated than Nrf2, because it is allowed for alternatively transcribing and selectively splicing to yield multiple isoforms (e.g., TCF11, Nrf1). In order to gain a better understanding of their similarities and differences in distinct regulatory profiles, all four distinct cell models for stably expressing TCF11, TCF11{Delta}N, Nrf1 or Nrf2 have been herein established by an Flp-In T-REx-293 system and then identified by transcriptomic sequencing. Further analysis revealed that Nrf1 and TCF11 have similar yet different regulatory profiles, although both contribute basically to positive regulation of their co-targets, which are disparate from those regulated by Nrf2. Such disparity in those gene regulation by Nrf1 and Nrf2 was further corroborated by scrutinizing comprehensive functional annotation of their specific and/or common target genes. Conversely, the mutant TCF11{Delta}N, resulting from a deletion of the N-terminal amino acids 2-156 from TCF11, resembles Nrf2 with the largely consistent structure and function. Interestingly, our further experimental evidence demonstrates that TCF11 acts as a potent tumor-repressor relative to Nrf1, albeit both isoforms possess a congruous capability to prevent malignant growth of tumor and upregulate those genes critical for improving the survival of patients with hepatocellular carcinoma. metapredict: a fast, accurate, and easy-to-use predictor of consensus disorder and structure Intrinsically disordered proteins and protein regions make up a substantial fraction of many proteomes where they play a wide variety of essential roles. A critical first step in understanding the role of disordered protein regions in biological function is to identify those disordered regions correctly. Computational methods for disorder prediction have emerged as a core set of tools to guide experiments, interpret results, and develop hypotheses. Given the multiple different predictors available, consensus scores have emerged as a popular approach to mitigate biases or limitations of any single method. Consensus scores integrate the outcome of multiple independent disorder predictors and provide a per-residue value that reflects the number of tools that predict a residue to be disordered. Although consensus scores help mitigate the inherent problems of using any single disorder predictor, they are computationally expensive to generate. They also necessitate the installation of multiple different software tools, which can be prohibitively difficult. To address this challenge, we developed a deep-learning-based predictor of consensus disorder scores. Our predictor, metapredict, utilizes a bidirectional recurrent neural network trained on the consensus disorder scores from 12 proteomes. By benchmarking metapredict using two orthogonal approaches, we found that metapredict is among the most accurate disorder predictors currently available. Metapredict is also remarkably fast, enabling proteome-scale disorder prediction in minutes. Importantly, metapredict is fully open source and is distributed as a Python package, a collection of command-line tools, and a web server, maximizing the potential practical utility of the predictor. We believe metapredict offers a convenient, accessible, accurate, and high-performance predictor for single-proteins and proteomes alike.

Statement of SignificanceIntrinsically disordered regions are found across all kingdoms of life where they play a variety of essential roles. Being able to accurately and quickly identify disordered regions in proteins using just the amino acid sequence is critical for the appropriate design and interpretation of experiments. Despite this, performing large-scale disorder prediction on thousands of sequences is challenging using extant disorder predictors due to various difficulties including general installation and computational requirements. We have developed an accurate, high-performance and easy-to-use predictor of protein disorder and structure. Our predictor, metapredict, was designed for both proteome-scale analysis and individual sequence predictions alike. Metapredict is implemented as a collection of local tools and an online web server, and is appropriate for both seasoned computational biologists and novices alike. Contribution Evenness: A functional redundancy metric sensitive to trait stability in microbialcommunities The concept of functional redundancy has received considerable attention in both the macroecology and microbial ecology literature. As a result, multiple metrics of functional redundancy have been proposed. These vary in how they weight trait levels, species abundance, functional richness, and species richness. Here we present a new functional redundancy metric tailored for community-aggregated traits, which are traits that are quantified at the community level and can be quantitatively partitioned among species. We call this metric Contribution Evenness (CE) because it measures how evenly species contribute to a community-aggregated trait. As CE is an evenness measurement, it ranges from 0 and 1, where 0 corresponds to a single species contributing to a community-aggregated trait and 1 corresponds to all species contributing equally. Using in silico simulations of species extinctions, we demonstrate that CE reflects the stability of an ecosystem function to species extinction, a hypothesized ecological consequence of functional redundancy. As a positive control and to illustrate how CE can be used with sequence data, we analyzed the functional redundancy of eight nitrogen-transforming pathways using 2,631 metagenome-assembled genomes from 47 marine sites. CE for marine nitrogen cycle marker genes was consistent with our qualitative understanding of which nitrogen pathways are most functionally redundant in the ocean. We found that, on average, the NH4+ assimilation pathway was the most functionally redundant (0.44 {+/-} 0.08) while dissimilatory nitrate reduction was the least redundant (0.005 {+/-} 0.005). As demonstrated here, CE provides a promising framework for measuring trait stability in microbiomes. Antibacterial potency of Type VI amidase effector toxins is dependent on substrate topology and cellular context Members of the bacterial T6SS amidase effector (Tae) superfamily of toxins are delivered between competing bacteria to degrade cell wall peptidoglycan. Although Taes share a common substrate, they exhibit distinct antimicrobial potency across different competitor species. To investigate the molecular basis governing these differences, we quantitatively defined the functional determinants of Tae1 from Pseudomonas aeruginosa PAO1 using a combination of nuclear magnetic resonance (NMR) and a high-throughput in vivo genetic approach called deep mutational scanning (DMS). As expected, combined analyses confirmed the role of critical residues near the Tae1 catalytic center. Unexpectedly, DMS revealed substantial contributions to enzymatic activity from a much larger, ring-like functional hot spot extending around the entire circumference of the enzyme. Comparative DMS across distinct growth conditions highlighted how functional contribution of different surfaces is highly context-dependent, varying alongside composition of targeted cell walls. These observations suggest that Tae1 engages with the intact cell wall network through a more distributed three-dimensional interaction interface than previously appreciated, providing an explanation for observed differences in antimicrobial potency across divergent Gram-negative competitors. Further binding studies of several Tae1 variants with their cognate immunity protein demonstrate that requirements to maintain protection from Tae activity may be a significant constraint on the mutational landscape of tae1 toxicity in the wild. In total, our work reveals that Tae diversification has likely been shaped by multiple independent pressures to maintain interactions with binding partners that vary across bacterial species and conditions. A Novel Proposal for an Index for Regional Cerebral Perfusion Pressure - A Theoretical Approach Using Fluid Dynamics Cerebral blood flow (CBF) / cerebral blood volume (CBV) ratio derived by [15O] H2O/ CO2 and CO positron emission tomography (PET) examination has been used empirically as an index for cerebral perfusion pressure (CPP). However, this equation lacks a theoretical basis and furthermore, has not been confirmed to be proportionate to CPP, as direct measurement of local CPP is not practical.

We have developed a new index for CPP using the Poiseuille equation based on a simple model. Our model suggests that CBF/CBV2 is proportionate to CPP and that it is mathematically a more accurate index than CBF/CBV. Sablefish (Anoplopoma fimbria Pallas, 1814) plasma biochemistry and hematology reference intervals including blood cell morphology Plasma biochemistry and hematology reference intervals are integral health assessment tools in all medical fields, including aquatic animal health. As sablefish (Anoplopoma fimbria) are becoming aquaculturally and economically more important, this manuscript provides essential reference intervals (RI) for their plasma biochemistry and hematology along with reference photomicrographs of blood cells in healthy, fasted sablefish. Blood cell morphology can differ between fish species. In addition, blood cell counts and blood chemistry can vary between fish species, demographics, water conditions, seasons, diets, and culture systems, which precludes the use of RIs from other fish species. For this study, blood was collected for plasma biochemistry and hematology analysis between June 20 and July 18, 2019, from healthy, yearling sablefish, hatched and reared in captivity on a commercial diet. Overnight fast of 16-18 hours did not sufficiently reduce lipids in the blood, which led to visible lipemia and frequent rupture of blood cells during analysis. Therefore, sablefish should be fasted for 24 to 36 hours before blood is collected to reduce hematology artifacts or possible reagent interference in plasma biochemistry analysis. Lymphocytes were the most dominant leukocytes (98%), while eosinophils were rare, and basophils were not detected in sablefish. Neutrophils were very large cells with Dohle bodies. In mammals and avian species, Dohle bodies are usually signs of toxic change from inflammation, but no such association was found in these fish. In conclusion, lipemia can interfere with sablefish blood analysis, and available removal methods should be evaluated as fasting for up to 36 h might not always be feasible. Also, more studies are required to establish RI for different developmental stages and rearing conditions. Agonist efficiency from concentration-response curves: structural implications and applications Agonists are evaluated by a concentration-response curve (CRC), with a midpoint (EC50) that indicates potency, a high-concentration asymptote that indicates efficacy and a low-concentration asymptote that indicates constitutive activity. A third agonist attribute, efficiency ({eta}), is the fraction of binding energy that is applied to the conformational change that activates the receptor. We show that {eta} can be calculated from EC50 and the asymptotes of a CRC derived from either single-channel or whole-cell responses. For 20 agonists of skeletal muscle nicotinic receptors, the distribution of {eta} values is bimodal with population means at 51% (including ACh, nornicotine and DMPP) and 40% (including epibatidine, varenicline and cytisine). The value of {eta} is related inversely to the size of the agonists head-group, with high-versus low-efficiency ligands having an average volume of 70 [A]3 versus 102 [A]3. Most binding site mutations have only a small effect on ACh efficiency except for Y190A (35%), W149A (60%) and those at G153 (42%). If {eta} is known, the midpoint and high-concentration asymptote can be calculated from each other. Hence, an entire CRC can be estimated from the response to a single agonist concentration, and efficacy can be estimated from EC50 of a CRC that has been normalized to 1. Given {eta}, the level of constitutive activity can be estimated from a single CRC.

Statement of significanceReceptors are molecular machines that convert chemical energy from agonist binding into mechanical energy of a global conformational change that generates a cell response. Agonists are distinguished by their potency (proportional to affinity) and efficacy, but also by the efficiency at which their binding energy is applied to receptor activation. Here we show that agonist efficiency can be estimated from a single concentrationresponse curve, and estimate efficiencies of 20 nicotinic receptor agonists. These have a bimodal distribution with larger agonists belonging to the lower-efficiency population. We further show that mutations of some binding site residues alter efficiency, and that knowledge of agonist efficiency simplifies and extends dose-response curve analysis. CellKb Immune: a manually curated database of mammalian hematopoietic marker gene sets for rapid cell type identification Using single-cell RNA sequencing, it is now possible to get the detailed gene expression pattern in each cell for millions of cells in a single experiment under varying conditions. Annotating cells by assigning a type or state to them is important for correct interpretation of results. In spite of current advances in technology, assigning cell types in single-cell datasets remains a bottleneck due to the lack of a comprehensive reference database and an accurate search method in a single tool. CellKb Immune is a knowledgebase of manually collected, curated and annotated marker gene sets representing cell types in the immune response in mouse. It allows users to find matching cell types in literature given a list of genes. We evaluated the contents and search methodology of CellKb Immune using a leave-one-out approach. We further used CellKb Immune to annotate previously defined marker gene sets from a bulk RNA-seq dataset to confirm its accuracy and coverage. CellKb Immune provides an easy-to-use online tool with a reliable method to find matching cell types and annotate cells in single-cell experiments. It is available at https://www.cellkb.com/immune. The structurome of a Clostridium difficile phage and the remarkable accurate prediction of its novel phage receptor-binding protein As natural bacterial predators, bacteriophages have the potential to be developed to tackle antimicrobial resistance, but our exploitation of them is limited by understanding their vast uncharacterised genetic diversity1,2. Fascinatingly, this genetic diversity reflects many ways that phages can make proteins, performing similar functions that together form the familiar phage particle. Critical to infection are phage receptor-binding proteins (RBPs) that bind bacterial  receptors and initiate bacterial entry3. Here we identified and characterised Gp22, a novel RBP for phage CDHS-1 that infects pathogenic C. difficile, but that had no recognisable RBPs. We showed that Gp22 antibodies neutralised CDHS-1 infection and used immunogold-labelling and transmission electron microscopy to identify their location on the capsid. The Gp22 three-dimensional structure was resolved by X-ray crystallography revealing a new RBP class with an N-terminal L-shaped -helical superhelix domain and a C-terminal Mg2+-binding domain. The findings provide novel insights into C. difficile phage biology and phage-host interactions. This will facilitate optimal phage development and future engineering strategies4,5. Furthermore, the AlphaFold2-predicted Gp22 structure, which was strikingly accurate, paves the way for a structurome based transformation and guidance of future phage studies where many proteins lack sequence homology but have recognisable protein structures. Damaraland mole-rats are not obligate cooperative breeders Damaraland mole-rats (Fukomys damarensis) are usually viewed as a eusocial or obligate cooperative breeder in which successful reproduction is dependent on help from closely related group members. However, because longitudinal studies of mole-rats in their natural environment are uncommon, the extent to which successful reproduction by breeders relies on assistance from other group members is unclear, and for non-breeders, the immediate and delayed fitness consequences of staying and helping are poorly understood. Using data from a 7-year field study of marked individuals, we investigate whether increasing group size confers reproductive and survival benefits to breeders and non-breeders, and explore the different routes through which individuals acquire breeding positions. We show that solitary natal dispersal was the most common route to breeding for both sexes and that the inheritance of dominant breeding positions was uncommon in both sexes. After dispersing, females typically settled alone in new burrow systems where they enjoyed high survival rates and remained in good body condition - often for several years - before being joined by males. In contrast to most obligately cooperative species, pairs of potential breeders reproduced successfully without helpers and experimentally formed pairs had the same reproductive success as larger established groups. Though larger breeding groups recruited slightly more pups on average, our data suggest that neither survival nor reproduction depend on the presence of non-breeding helpers, indicating that Damaraland mole-rats are not obligate cooperative breeders. We suggest that extended philopatry and group living in Damaraland mole-rats have evolved because of the high costs and constraints of dispersal rather than because of strong indirect benefits accrued through cooperative behaviour and that similarities between their breeding systems and those of obligatorily eusocial insects have been over-emphasized.

Significance StatementThe social mole-rats are often considered eusocial mammals in which successful reproduction depends on assistance from non-breeding helpers. In this study we show that in wild Damaraland mole-rats, the presence of non-breeders is associated with both costs and benefits and that nascent breeding pairs show high reproductive success despite the lack of non-breeding helpers. These findings indicate that Damaraland mole-rats are not obligate cooperative breeders and suggest that similarities between their breeding systems and those of obligatorily eusocial insects have been over-emphasized. Multiplex knockout of trichome-regulating MYB duplicates in hybrid poplar using a single gRNA As the focus for CRISPR edited plants moves from proof-of-concept to real world applications, precise gene manipulation will increasingly require concurrent multiplex editing for polygenic traits. A common approach for editing across multiple sites is to design one gRNA per target; however, this complicates construct assembly and increases the possibility of off-target mutations. In this study, we utilized one gRNA to target MYB186, a known positive trichome regulator, as well as its paralogs MYB138 and MYB38 at a consensus site for mutagenesis in Populus tremula x P. alba INRA 717-1B4. Unexpected duplications of MYB186 and MYB138 resulted in a total of eight alleles for the three targeted genes in the hybrid poplar. Deep sequencing and PCR analyses confirmed editing across all eight targets in nearly all of the resultant glabrous mutants, ranging from small indels to large genomic dropouts, with no off-target activity detected at four potential sites. This highlights the effectiveness of a single gRNA targeting conserved exonic regions for multiplex editing. Additionally, cuticular wax and whole leaf analyses showed a complete absence of triterpenes in the trichomeless mutants, hinting at a previously undescribed role for the non-glandular trichomes of poplar.

ONE SENTENCE SUMMARYTargeting conserved sequences with a single gRNA allowed efficient mutagenesis of a multigene family and the recovery of trichomeless and triterpene-free poplar mutants. We get by with a little help from our friends: shared adaptive variation provides a bridge to novel ecological specialists during adaptive radiation Adaptive radiations involve astounding bursts of phenotypic, ecological, and species diversity. However, the microevolutionary processes that underlie the origins of these bursts are still poorly understood. We report the discovery of an intermediate  wide-mouth scale-eating ecomorph in a sympatric radiation of Cyprinodon pupfishes, illuminating the transition from a widespread algae-eating generalist to a novel microendemic scale-eating specialist. We first show that this ecomorph occurs in sympatry with generalist C. variegatus and scale-eating specialist C. desquamator on San Salvador Island, Bahamas, but is genetically differentiated, morphologically distinct, and often consumes scales. We then compared the timing of selective sweeps on shared and unique adaptive variants in trophic specialists to characterize their adaptive walk. Shared adaptive regions swept first in both the specialist desquamator and the intermediate  wide-mouth ecomorph, followed by unique sweeps of introgressed variation in  wide-mouth and de novo variation in desquamator. The two scale-eating populations additionally shared 9% of their hard selective sweeps with molluscivores C. brontotheroides, despite no single common ancestor among specialists. Our work provides a new microevolutionary framework for investigating how major ecological transitions occur and illustrates how both shared and unique genetic variation can provide a bridge for multiple species to access novel ecological niches. Unevolved proteins from modern and prebiotic amino acids manifest distinct structural profiles Natural proteins represent numerous but tiny structure/function islands in a vast ocean of possible protein sequences, most of which has not been explored by either biological evolution or research. Recent studies have suggested this uncharted sequence space possesses surprisingly high structural propensity, but development of an understanding of this phenomenon has been awaiting a systematic high-throughput approach.

Here, we designed, prepared, and characterized two combinatorial protein libraries consisting of randomized proteins, each 105 residues in length. The first library constructed proteins from the entire canonical alphabet of 20 amino acids. The second library used a subset of only 10 residues (A,S,D,G,L,I,P,T,E,V) that represent a consensus view of plausibly available amino acids through prebiotic chemistry. Our study shows that compact conformations resistant to proteolysis are (i) abundant (up to 40%) in random sequence space, (ii) independent of general Hsp70 chaperone system activity, and (iii) not granted solely by "late" and complex amino acid additions. The Hsp70 chaperone system effectively increases solubility and refoldability of the canonical alphabet but has only a minor impact on the "early" library. The early alphabet proteins are inherently more soluble and refoldable, possibly assisted by the cell-like environment in which these assays were performed.

Our work indicates that both early and modern amino acids are predisposed to supporting protein structure (either in forms of oligomers or globular/molten globule structures) and that protein structure may not be a unique outcome of evolution. Diversity of SARS-CoV-2 genome among various strain identified in Lucknow, Uttar Pradesh Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has emerged as a significant challenge worldwide. Rapid genome sequencing of SARS-CoV-2 is going on across the globe to detect mutations and genomic modifications in SARS-CoV-2. In this study, we have sequenced twenty-three SARS-CoV-2 positive samples collected during the first pandemic from the state of Uttar Pradesh, India. We observed a range of already reported mutations (2-22), including; D614G, L452R, Q613H, Q677H, T1027I in the S gene; S194L in the N gene; Q57H, L106F, T175I in the ORF3. Few unreported mutations such as P309S in the ORF1ab gene; T379I in the N gene; and L52F, V77I in the ORF3a gene were also detected. Phylogenetic genome analysis showed similarity with other SARS-CoV-2 viruses reported from Uttar Pradesh. The observed mutations may be associated with SARS-CoV-2 virus pathogenicity or disease severity.