Evaluation of US state pollinator plans using 3 evidence-based policymaking frameworks After the US federal government created a national pollinator protection plan in 2015, many states followed with their own. Since their goal is to promote pollinating insect conservation, we wanted to know whether the state plans are using best practices for evidence-based science policy. In early 2019 we found and downloaded every existing, publicly available US state pollinator protection plan. We then used content analysis to assess the goals, scope, and implementation of state-level pollinator protection plans across the US. This analysis was conducted using three distinct frameworks for evidence-based policymaking: US Department of Interior Adaptive Resources Management (ARM), US Environmental Protection Agency management pollinator protection plan (MP3) guidance, and Pew Trusts Pew-MacAthur Results First Project elements of evidence-based state policymaking (PEW) framework. Then we scored them using the framework criteria, to assess whether the plans were using known best practices for evidence based policymaking. Of the 31 states with a state pollinator plan, Connecticut was the state with the lowest total score across the three evaluation frameworks. The state with the highest overall scores, across the three frameworks, was Missouri. Most states did not score highly on the majority of the frameworks. Overall, many state plans were lacking policy elements that address monitoring, evaluation, and adjustment. These missing elements impact states ability to achieve their conservation goals. Our results indicate that states can improve their pollinator conservation policies to better match evidence-based science policy guidance, regardless of which framework is used.
Effects of phenotypic variation on consumer coexistence and prey community structure A popular idea in ecology is that trait variation among individuals from the same species may promote the coexistence of competing species. However, theoretical and empirical tests of this idea have yielded inconsistent findings. We manipulated intraspecific trait diversity in a ciliate competing with a nematode for bacterial prey in experimental microcosms. We found that intraspecific trait variation inverted the original competitive hierarchy to favor the consumer with variable traits, ultimately resulting in competitive exclusion. This competitive outcome was driven by foraging traits (size, speed, and directionality) that increased the ciliates fitness ratio and niche overlap with the nematode. The interplay between consumer trait variation and competition resulted in non-additive cascading effects - mediated through prey defense traits - on prey community assembly. Our results suggest that predicting consumer competitive population dynamics and the assembly of prey communities will require understanding the complexities of trait variation within consumer species.
Distinguishing excess mutations and increased cell death based on variant allele frequencies Tumors often harbor orders of magnitude more mutations than healthy tissues. The increased number of mutations may be due to an elevated mutation rate or frequent cell death and correspondingly rapid cell turnover, or a combination of the two. It is difficult to disentangle these two mechanisms based on widely available bulk sequencing data, where sequences from individual cells are intermixed and, thus, the cell lineage tree of the tumor cannot be resolved. Here we present a method that can simultaneously estimate the cell turnover rate and the rate of mutations from bulk sequencing data. Our method works by simulating tumor growth and finding the parameters with which the observed data can be reproduced with maximum likelihood. Applying this method to a real tumor sample, we find that both the mutation rate and the frequency of death may be high.

Author SummaryTumors frequently harbor an elevated number of mutations, compared to healthy tissue. These extra mutations may be generated either by an increased mutation rate or the presence of cell death resulting in increased cellular turnover and additional cell divisions for tumor growth. Separating the effects of these two factors is a nontrivial problem. Here we present a method which can simultaneously estimate cell turnover rate and genomic mutation rate from bulk sequencing data. Our method is based on the estimation of the parameters of a generative model of tumor growth and mutations. Applying our method to a human hepatocellular carcinoma sample reveals an elevated per cell division mutation rate and high cell turnover. Using total RNA quality metrics for time since deposition estimates in degrading bloodstains Determining the time since deposition (TSD) of bloodstains would provide forensic scientists with critical information regarding the timeline of the events involving bloodshed. The physicochemical changes occurring to biomolecules as a bloodstain degrades can be used to approximate the TSD of bloodstains. Our study aims to quantify the timewise degradation trends and temperature dependence found in total RNA from bloodstains without the use of amplification, expanding the scope of the RNA TSD research which has previously targeted mRNA molecules. Whole bovine blood was stored in plastic microcentrifuge tubes at 21{degrees}C or 4{degrees}C and tested over different timepoints spanning one week. Nine RNA metrics were visually assessed and quantified using linear and mixed models; the RNA Integrity Number equivalent (RINe) and the DV200 demonstrated strong negative trends over time and statistical independence. The RINe model fit was high (R2 = 0.60), and while including the biological replicate as a random effect increased the fit for all RNA metrics, no significant differences were found between biological replicates stored at the same temperature for the RINe and DV200 metrics. Importantly, this suggests that these standardized metrics can likely be directly compared between scenarios and individuals, with DV200 having an inflection point at [~]28 hrs. This study provides a novel approach for blood TSD estimates, producing metrics that are not affected by inter-individual variation and improving our understanding of the rapid degradation occurring in bloodstains.

HIGHLIGHTSO_LIAmplification-free analysis of total RNA in degrading bloodstains.
C_LIO_LIShort-term RNA degradation assessment using high-resolution size measurements.
C_LIO_LITotal RNA quality and quantity metrics were assessed across a one-week.
C_LIO_LITotal RNA quality metrics demonstrated the strongest timewise trends.
C_LIO_LIBiological replicates produced similar results for RNA quality metrics.
C_LI 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. Affordable Caenorhabditis elegans tracking system for classroom use For decades the nematode C. elegans has served as an outstanding research organism owing to its unsurpassed experimental amenability. This advantage has also made this tiny worm an attractive vehicle for science instruction across higher learning institutions. However, the prohibitive cost associated with the automated behavioral assessment of these animals remains an obstacle preventing their full adoption in undergraduate and high school settings. To improve this situation, we developed an inexpensive worm tracking system for use by high school interns and undergraduate students. Over the past two years this tracker has been successfully used by undergraduate students in our introductory Cell and Molecular lab (BSC220) at Illinois State University. Here we describe and demonstrate the use of our inexpensive worm tracking system. Phylogenomic analyses of Alismatales shed light into adaptations to aquatic environments Land plants first evolved from freshwater algae, and flowering plants returned to water as early as the Cretaceous and multiple times subsequently. Alismatales is the largest clade of aquatic angiosperms including all marine angiosperms, as well as terrestrial plants. We used Alismatales to explore plant adaptation to aquatic environments by analyzing a data set that included 95 samples (89 Alismatales species) covering four genomes and 91 transcriptomes (59 generated in this study). To provide a basis for investigating adaptation, we assessed phylogenetic conflict and whole-genome duplication (WGD) events in Alismatales. We recovered a relationship for the three main clades in Alismatales as (Tofieldiaceae, Araceae) + core Alismatids. We also found phylogenetic conflict among the backbone of the three main clades that could be explained by incomplete lineage sorting and introgression. Overall, we identified 18 putative WGD events across Alismatales. One of them occurred at the most recent common ancestor of core Alismatids, and three occurred at seagrass lineages. We also found that lineage and life-form were both important for different evolutionary patterns for the genes related to freshwater and marine adaptation. For example, several light or ethylene-related genes were lost in the seagrass Zosteraceae, but are present in other seagrasses and freshwater species. Stomata-related genes were lost in both submersed freshwater species and seagrasses. Nicotianamine synthase genes, which are important in iron intake, expanded in both submersed freshwater species and seagrasses. Our results advance the understanding of the adaptation to aquatic environments and whole-genome duplications using phylogenomics. The Skeletal Muscle Circadian Clock Regulates Titin Splicing Through RBM20 Circadian rhythms are maintained by a cell autonomous, transcriptional-translational feedback loop known as the molecular clock. While previous research suggests a role of the molecular clock in regulating skeletal muscle structure and function, no mechanisms have connected the molecular clock to sarcomere filaments. Utilizing inducible, skeletal muscle specific, Bmal1 knockout (iMSBmal1-/-) mice, we showed that knocking out skeletal muscle clock function alters titin isoform expression through both RNAseq and LC-MS. This alteration in titins spring length resulted in sarcomere length heterogeneity. We demonstrate the direct link between altered titin splicing and sarcomere length in vitro using U7 snRNPs that truncate the region of titin altered in iMSBmal1-/- muscle. Finally, we identified a mechanism whereby the skeletal muscle clock regulates titin isoform expression through transcriptional regulation of Rbm20, a potent splicing regulator of titin. Our findings demonstrate the importance of the skeletal muscle circadian clock in maintaining titin isoform through regulation of RBM20 expression. Because circadian rhythm disruption is a feature of many chronic diseases, our results highlight a novel pathway that could be targeted to maintain skeletal muscle structure and function in a range of pathologies. First evidence of in vitro cytotoxic effects of marine microlitter on Merluccius merluccius and Mullus barbatus, two Mediterranean commercial fish species Marine litter is composed mainly of plastics and is recognized as a serious threats to marine ecosystems. Ecotoxicological approaches have started elucidating the potential severity of microplastics (MPs) in controlled laboratory studies with pristine materials but no information exist on marine environmental microlitter as a whole. Here, we characterized the litter in the coastal Northern Tyrrhenian sea and in the stomach of two fish species of socio-economic importance, and exposed primary cell cultures of mucosal and lymphoid organs to marine microlitter for evaluating possible cytotoxic effects. An average of 0.30 {+/-} 0.02 microlitter items m-3 was found in water samples. FT-IR analysis revealed that plastic particles, namely HDPE, polyamide and polypropylene were present in 100% and 83.3% of Merluccius merluccius and Mullus barbatus analyzed, which overall ingested 14.67 {+/-} 4.10 and 5.50 {+/-} 1.97 items/individual, respectively. Moreover, microlitter was confirmed as a vector of microorganisms. Lastly, the apical end-point of viability was found to be significantly reduced in splenic cells exposed in vitro to two microlitter conditions. Considering the role of the spleen in the mounting of adaptive immune responses, our results warrant more in-depth investigations for clarifying the actual susceptibility of these two species to anthropogenic microlitter.

HighlightsO_LI0.30 {+/-} 0.02 microlitter items m-3 were found at the surface of coastal Northern Tyrrhenian sea
C_LIO_LI14.67 {+/-} 4.10 and 5.50 {+/-} 1.97 items/individual were retrieved from the stomach of hakes and mullets
C_LIO_LIThe ingested microlitter contained plastic items
C_LIO_LIMicrolitter was validated as a carrier of bacteria, fungi and flagellates
C_LIO_LISplenic cells exposed to two microlitter conditions for 72 hours suffered cytotoxicity
C_LI Different contribution of the monkey prefrontal and premotor dorsal cortex in decision making during a transitive inference task Several studies have reported similar neural modulations between brain areas of the frontal cortex, such as the dorsolateral prefrontal (DLPFC) and the premotor dorsal (PMd) cortex, in tasks requiring encoding of the abstract rules for selecting the proper action. Here we compared the neuronal modulation of the DLPFC and PMd of monkeys trained to choose the higher rank from a pair of abstract images (target item), selected from an arbitrarily rank-ordered set (A>B>C>D>E>F) in the context of a transitive inference task. Once acquired by trial-and-error, the ordinal relationship between pairs of adjacent images (i.e., A>B; B>C; C>D; D>E; E>F), monkeys were tested in indicating the ordinal relation between items of the list not paired during learning. During these decisions, we observed that the choice accuracy increased and the reaction time decreased as the rank difference between the compared items enhanced. This result is in line with the hypothesis that after learning, the monkeys built an abstract mental representation of the ranked items, where rank comparisons correspond to the items position comparison on this representation. In both brain areas, we observed higher neuronal activity when the target item appeared in a specific location on the screen with respect to the opposite position and that this difference was particularly enhanced at lower degrees of difficulty. By comparing the time evolution of the activity of the two areas, we observed that the neural encoding of target item spatial position occurred earlier in the DLPFC than in the PMd. 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. Shelterin components modulate nucleic acids condensation and phase separation in the context of telomeric DNA Telomeres are nucleoprotein complexes that protect the ends of chromosomes and are essential for chromosome stability in Eukaryotes. In cells, individual telomeres form distinct globules of finite size that appear to be smaller than expected for bare DNA. Moreover, upon changes in their protein composition, telomeres can cluster to form telomere-induced-foci (TIFs) or co-localize with promyelocytic leukemia (PML) nuclear bodies. The physical basis for collapse of individual telomeres and coalescence of multiple ones remains unclear, as does the relationship between these two phenomena. By combining single-molecule measurements, optical microscopy, turbidity assays, and simulations, we show that the telomere scaffolding protein TRF2 can condense individual DNA chains and drives coalescence of multiple DNA molecules, leading to phase separation and the formation of liquid-like droplets. Addition of the TRF2 binding protein hRap1 modulates phase boundaries and tunes the specificity of solution demixing while simultaneously altering the degree of DNA compaction. Our results suggest that the condensation of single telomeres and formation of biomolecular condensates containing multiple telomeres are two different outcomes driven by the same set of molecular interactions. Moreover, binding partners, such as other telomere components, can alter those interactions to promote single-chain DNA compaction over multiple-chain phase separation. Estimation of cell cycle kinetics in higher plant root meristem links organ position with cellular fate and chromatin structure Plant root development is a complex spatial-temporal process that originates in the root apical meristem (RAM). To shape the organs structure signaling between the different cells and cell files must be highly coordinated. Thereby, diverging kinetics of chromatin remodeling and cell growth in these files need to be integrated and balanced by differential cell growth and local differences in cell proliferation frequency. Understanding the local differences in cell cycle duration in the RAM and its correlation with chromatin organization is crucial to build a holistic view on the different regulatory processes and requires a quantitative estimation of the chromatin geometry and underlying mitotic cell cycle phases timing at every cell file and every position. Unfortunately, so far precise methods for such analysis are missing.

This study presents a robust and straightforward pipeline to determine in parallel the duration of cell cycles key stages in all cell layers of a plants root and their nuclei organization. The methods combine marker-free techniques based on the detection of the nucleus, deep analysis of the chromatin phase transition, incorporation of 5-ethynyl-2'-deoxyuridine (EdU), and mitosis with a deep-resolution plant phenotyping platform to analyze all key cell cycle events kinetics.

In the Arabidopsis thaliana L. RAM S-phase duration was found to be as short as 20-30 minutes in all cell files. The subsequent G2-phase duration however depends on the cell type/position and varies from 3.5 hours in the pericycle to more than 4.5 hours in the epidermis. Overall, S+G2+M duration in Arabidopsis under our condition is 4 hours in the pericycle and up to 5.5 hours in the epidermis.

Endocycle duration was determined as the time required to achieve 100% EdU index in the transition zone and estimated to be in the range of 3-4 hours.

Besides Arabidopsis, we show that the presented technique is applicable also to root tips of other dicot and monocot plants (tobacco (Nicotiana tabacum L.), tomato (Lycopersicon esculentum L.) and wheat (Triticum aestivum L.). Paradoxical Roles of Peritumoral Myofibroblasts in Intrahepatic Cholangiocarcinoma Growth and Metastasis Intrahepatic cholangiocarcinoma (iCCA) is an aggressive type of primary liver cancer characterized by its highly desmoplastic stroma. Compared to the ample knowledge on cancer-associated fibroblast residing within iCCA tumor mass, little is known about the fibrotic response from the tumor-surrounding liver and its role in iCCA development. In this study, we find a significant accumulation of peritumoral myofibroblasts (pMFs) in both patient iCCA and tumors from an iCCA orthotopic allograft mouse model. Using tumor and liver spheroid coculture we show that iCCA-liver interaction induces rapid pMF accumulation at the interface. We find pMFs placed around iCCA spheroids exert a strong suppressive effect on tumor cells growth, in contrast to the pro-proliferative effect of MFs mixed within tumor spheroids. However, prolonged iCCA-pMF interaction elicits tumor cell dissemination in vitro. We find an upregulation of Vcam1 in tumor cells in the early phase of iCCA-pMF interaction both in vitro and in vivo which is downregulated when tumor cells disseminate. Blocking Vcam1 activity in iCCA allograft mouse models slows primary tumor growth but lead to increased tumor metastasis. Our data suggest that pMFs are beyond simple pro- or anti-tumorigenic in iCCA, with the ability to suppress tumor growth but elicit tumor cell dissemination. Paternal aging affects developmental convergence of vocal behavior in individual mice: machine learning-driven analyses of individuality Infant crying is an innate communicative behavior that is frequently impaired in certain neurodevelopmental disorders (NDDs). Since advanced paternal age is a reported risk factor for NDDs in offspring, we evaluated the impact of a fathers age on early vocal development in C57BL/6J mice. We recorded and applied a unique combination of computational analyses to ultrasonic vocalizations (USVs) emitted by mouse pups sired by young and aged fathers. Our data showed that advanced paternal age reduced the number and duration of USVs, and altered the syllable composition in pups. Moreover, pups born to young fathers showed convergent vocal characteristics with a rich repertoire during postnatal development, while those born to aged fathers exhibited more divergent vocal patterns with limited repertoire. Principal component analysis in conjunction with clustering analysis demonstrated that pups from aged fathers deviated from typical trajectories of vocal development, which were considered as atypical individuals. Thus, our study indicates that advanced paternal age has a significant effect on offsprings early vocal development. It is suggested that the trajectories of vocal development could be a useful marker of the NDD-like phenotype associated with the advanced paternal age. In addition, our comprehensive computational analysis described here is an effective approach to characterize the altered individual diversity relevant to neurodevelopmental disorders.

One Sentence SummaryAdvanced paternal age affects vocal development in early postnatal mice, with more pups showing atypical developmental trajectories. In vitro pectate lyase activity and carbon uptake assays and whole genome sequencing of Bacillus amyloliquefaciens subsp. plantarum strains for a pectin defective pathway The pectin lyase activity of 59 Bacillus amyloliquefaciens subsp. plantarum (Bap) strains was tested in vitro on Pectate Agar (PA) and Tris-Spizizen Salts (TSS) medium. Bap strains were cultured on TSA medium and washed three times with sterile water before the inoculation on PA media. Higher and lower pectate lyase activity were observed in six (AP193, AP203, AP299, AP80, AP102, and AP52) and four (AP 194, AP214, AP215, and AP305) Bap strains compared to other Bap strains. A total of 12 Bap strains (AP67, AP71, AP77, AP78, AP85, AP102, AP108, AP135, AP143, AP189, AP192, and AP193) grew vigorously on TSS medium. A total of six Bap strains (AP194, AP204, AP214, AP216, AP219, and HD73) had lower growth compared to other Bap strains. Pectin (1%) were used for in vitro PA and TSS medium. Pectate lyase and utilization activity were not found in Bacillus thuringiensis subsp. kurstaki strain HD73 compared to Bap strains. A draft genome sequence for strains AP194 and AP214 that were negative for pectin utilization were generated using an Illumina MiSeq. RAST analysis revealed that the pectin-associated gene altronate hydrolase (uxaA) absent in AP 214 strain. Multiple amino acid alignments of exuT and uxuB gene sequence showed dissimilarities among AP194, AP214, and reference Bap strains. 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. Transcriptomic analysis of temporal shifts in berry development between two grapevine cultivars of the Pinot family reveals potential genes controlling ripening time BackgroundGrapevine cultivars of the Pinot family represent clonally propagated mutants with major phenotypic and physiological differences, such as different colour or shifted ripening time, as well as changes in important viticultural traits. Specifically, the cultivars  Pinot Noir (PN) and  Pinot Noir Precoce (PNP, early ripening) flower at the same time, but vary in the beginning of berry ripening (veraison) and, consequently, harvest time. In addition to genotype, seasonal climatic conditions (i.e. high temperatures) also affect ripening times. To reveal possible regulatory genes that affect the timing of veraison onset, we investigated differences in gene expression profiles between PN and PNP throughout berry development with a closely meshed time series and over two separate years.

ResultsThe difference in the duration of berry formation between PN and PNP was quantified to be approximately two weeks under the growth conditions applied, using plant material with a proven PN and PNP clonal relationship. Clusters of co-expressed genes and differentially expressed genes (DEGs) were detected which reflect the shift in the timing of veraison onset. Functional annotation of these DEGs fit to observed phenotypic and physiological changes during berry development. In total, we observed 3,342 DEGs in 2014 and 2,745 DEGs in 2017 between PN and PNP, with 1,923 DEGs across both years. Among these, 388 DEGs were identified as veraison-specific and 12 were considered as berry ripening time regulatory candidates. The expression profiles revealed two candidate genes for ripening time control which we designated VviRTIC1 and VviRTIC2 (VIT_210s0071g01145 and VIT_200s0366g00020, respectively). These genes likely contribute the phenotypic differences observed between PN and PNP.

ConclusionsMany of the 1,923 DEGs show highly similar expression profiles in both cultivars if the patterns are aligned according to developmental stage. In our work, putative genes differentially expressed between PNP and PN which could control ripening time as well as veraison-specific genes were identified. We point out connections of these genes to molecular events during berry development and discuss potential candidate genes which may control ripening time. Two of these candidates were observed to be differentially expressed in the early berry development phase. Several down-regulated genes during berry ripening are annotated as auxin response factors / ARFs. Conceivably, general changes in auxin signaling may cause the earlier ripening phenotype of PNP. Exploring functional protein covariation across single cells using nPOP Many biological processes, such as the cell division cycle, are reflected in protein covariation across single cells. This covariation can be quantified and interpreted by single-cell proteomics with sufficiently high throughput and accuracy. Toward this goal, we developed the nano-ProteOmic sample Preparation (nPOP) method for single-cell proteomics. nPOP uses piezo acoustic dispensing to isolate individual cells in 300 picoliter volumes and performs all subsequent preparation steps in small droplets on a fluorocarbon-coated slide. This design enables simultaneous sample preparation of thousands of single cells, including lysing, digesting, and labeling individual cells in volumes below 20 nl. We used nPOP to prepare 1,888 single cells and 128 negative controls in a single batch. Their analysis enabled quantifying the covariation between thousands of proteins and cell-cycle protein markers. Many protein sets covaried with the cell cycle similarly across all cell types and states, reflecting cell-type independent cell cycle functions. However, the cell cycle covariation of other protein sets differed markedly between cell types, even within subpopulation of melanoma cells expressing markers for drug-resistance priming. The cells expressing these markers accumulated in the G1 phase of the cell cycle and exhibited different covariation of enzymes catabolizing glucose. These results demonstrate that protein covariation across single cells may reveal functionally concerted biological differences between closely related cell states.

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=59 SRC="FIGDIR/small/441211v3_ufig1.gif" ALT="Figure 1">
View larger version (53K):
org.highwire.dtl.DTLVardef@1e23b1org.highwire.dtl.DTLVardef@b93faeorg.highwire.dtl.DTLVardef@3c00forg.highwire.dtl.DTLVardef@362b41_HPS_FORMAT_FIGEXP  M_FIG C_FIG 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. SARS-CoV-2 Variants Reveal Features Critical for Replication in Primary Human Cells Since entering the human population, SARS-CoV-2 (the causative agent of COVID-19) has spread worldwide, causing >100 million infections and >2 million deaths. While large-scale sequencing efforts have identified numerous genetic variants in SARS-CoV-2 during its circulation, it remains largely unclear whether many of these changes impact adaptation, replication or transmission of the virus. Here, we characterized 14 different low-passage replication-competent human SARS-CoV-2 isolates representing all major European clades observed during the first pandemic wave in early 2020. By integrating viral sequencing data from patient material, virus stocks, and passaging experiments, together with kinetic virus replication data from non-human Vero-CCL81 cells and primary differentiated human bronchial epithelial cells (BEpCs), we observed several SARS-CoV-2 features that associate with distinct phenotypes. Notably, naturally-occurring variants in Orf3a (Q57H) and nsp2 (T85I) were associated with poor replication in Vero-CCL81 cells but not in BEpCs, while SARS-CoV-2 isolates expressing the Spike D614G variant generally exhibited enhanced replication abilities in BEpCs. Strikingly, low-passage Vero-derived stock preparation of 3 SARS-CoV-2 isolates selected for substitutions at positions 5/6 of E, and were highly attenuated in BEpCs, revealing a key cell-specific function to this region. Rare isolate-specific deletions were also observed in the Spike furin-cleavage site during Vero-CCL81 passage, but these were rapidly selected against in BEpCs, underscoring the importance of this site for SARS-CoV-2 replication in primary human cells. Overall, our study uncovers sequence features in SARS-CoV-2 variants that determine cell-specific virus replication, and highlights the need to monitor SARS-CoV-2 stocks carefully when phenotyping newly emerging variants or potential variants-of-concern. Cell type-specific inference of differential expression in spatial transcriptomics Spatial transcriptomics enables spatially resolved gene expression measurements at near single-cell resolution. There is a pressing need for computational tools to enable the detection of genes that are differentially expressed (DE) within specific cell types across tissue context. We show that current approaches cannot learn cell type-specific DE due to changes in cell type composition across space and the fact that measurement units often detect transcripts from more than one cell type. Here, we introduce a statistical method, Cell type-Specific Inference of Differential Expression (C-SIDE), that identifies cell type-specific patterns of differential gene expression while accounting for localization of other cell types. We model spatial transcriptomics gene expression as an additive mixture across cell types of general log-linear cell type-specific expression functions. This approach provides a unified framework for defining and identifying gene expression changes in a wide-range of relevant contexts: changes due to pathology, anatomical regions, physical proximity to specific cell types, and cellular microenvironment. Furthermore, our approach enables statistical inference across multiple samples and replicates when such data is available. We demonstrate, through simulations and validation experiments on Slide-seq and MER-FISH datasets, that our approach accurately identifies cell type-specific differential gene expression and provides valid uncertainty quantification. Lastly, we apply our method to characterize spatially-localized tissue changes in the context of disease. In an Alzheimers mouse model Slide-seq dataset, we identify plaque-dependent patterns of cellular immune activity. We also find a putative interaction between tumor cells and myeloid immune cells in a Slide-seq tumor dataset. We make our C-SIDE method publicly available as part of the open source R package https://github.com/dmcable/spacexr. Depletion of Sun1/2 induces heterochromatin accrual in mesenchymal stem cells during adipogenesis Critical to the mechano-regulation of mesenchymal stem cells (MSC), Linker of the Nucleoskeleton and Cytoskeleton (LINC) complex transduces cytoskeletal forces to the nuclei. The LINC complex contains outer nuclear membrane Nesprin proteins that associate with the cytoskeleton and their inner nuclear membrane couplers, Sun proteins. In addition to coupling Nesprin-associated cytoskeletal elements to inner nuclear membrane, Sun proteins also function in regulating nuclear mechanics and chromatin tethering to inner nuclear membrane. This suggests that release of LINC-mediated cytoskeletal connections from cell nuclei may have different effects on chromatin organization and MSC differentiation than those due to ablation of intra-nuclear Sun proteins. To test this hypothesis we compared Sun1/2 depletion with expression of the dominant-negative KASH domain (dnKASH) that inhibits Nesprin-Sun association. In cells cultured under adipogenic conditions, disconnecting LINC from the cytoskeleton through dnKASH expression, increased adipogenic gene expression and fat droplet formation; heterochromatin H3K27 and H3K9 methylation was unaltered. In contrast, Sun1/2 depletion inhibited adipogenic gene expression and fat droplet formation; as well the anti-adipogenic effect of Sun1/2 depletion was accompanied by increased H3K9me3, which was enriched on Adipoq, silencing this fat locus. We conclude that releasing the nucleus from cytoskeletal constraints via dnKASH accelerates adipogenesis while depletion of Sun1/2 increases heterochromatin accrual on adipogenic genes in a fashion independent of LINC complex function. Therefore, while these two approaches both disable LINC functions, their divergent effects on the epigenetic landscape indicate they cannot be used interchangeably to study mechanical regulation of cell differentiation. Structure-Functional-Selectivity Relationship Studies on A-86929 Analogs and Small biaryl Fragments toward Discovery of biased D1 agonists Dopamine regulates normal functions such as movement, reinforcement learning, and cognition, and its dysfunction has been implicated in multiple psychiatric and neurological disorders. Dopamine acts through the D1- (D1R and D5R) and D2-class (D2R, D3R and D4R) of seven transmembrane receptors, and activates both G-protein- and {beta}-arrestin-dependent signaling pathways, to mediate its physiological effects. Current dopamine receptor-based therapies are used to ameliorate motor deficits in Parkinsons disease, or as antipsychotic medications for schizophrenia. These drugs show efficacy for ameliorating only some symptoms caused by dopamine dysfunction and are plagued by debilitating side-effects. Studies in primates and rodents have shown that shifting the balance of dopamine receptor signaling towards the arrestin pathway can be beneficial for inducing normal movement, while reducing motor side-effects such as dyskinesias, and can be efficacious at enhancing cognitive function compared to balanced agonists. Several structure-activity-relationship (SAR) studies have embarked on discovering {beta}-arrestin-biased dopamine agonists, focused on D2 partial agonists, non-catechol D1 agonists, and mixed D1/D2R dopamine receptor agonists. Here, we describe an SAR study to identify novel D1R {beta}-arrestin biased ligands using A-86929, a high-affinity D1R catechol agonist, as a core scaffold. Previously described and novel analogs of A-86929 were synthesized and screened in vitro for structure-functional-selectivity relationships (SFSR) studies to identify chemical motifs responsible for {beta}-arrestin biased activity at both D1 and D2Rs. Most of the A-86929 analogs screened were G protein biased but none of them were exclusively arrestin-biased. Additionally, various catechol aryl fragments were designed and synthesized. Other compounds surveyed included hydroxyl and chloro analogs of dopamine to test for hydrogen bonding and ionic interactions. Some of these small molecular probes displayed weak bias towards the {beta}-arrestin pathway. Continued in-depth SFSR studies informed by structure determination, molecular modeling, and mutagenesis studies will facilitate discovery of potent and efficacious arrestin-biased dopamine receptor ligands. Ensovibep, a novel trispecific DARPin candidate that protects against SARS-CoV-2 variants SARS-CoV-2 has infected millions of people globally and continues to undergo evolution. Emerging variants can be partially resistant to vaccine induced immunity and therapeutic antibodies, emphasizing the urgent need for accessible, broad-spectrum therapeutics. Here, we report a comprehensive study of ensovibep, the first trispecific clinical DARPin candidate, that can simultaneously engage all three units of the spike protein trimer to potently inhibit ACE2 interaction, as revealed by structural analyses. The cooperative binding of the individual modules enables ensovibep to retain inhibitory potency against all frequent SARS-CoV-2 variants, including Omicron BA.1 and BA.2, as of February 2022. Moreover, viral passaging experiments show that ensovibep, when used as a single agent, can prevent development of escape mutations comparably to a cocktail of monoclonal antibodies (mAb). Finally, we demonstrate that the very high in vitro antiviral potency also translates into significant therapeutic protection and reduction of pathogenesis in Roborovski dwarf hamsters infected with either the SARS-CoV-2 wild-type or the Alpha variant. In this model, ensovibep prevents fatality and provides substantial protection equivalent to the standard of care mAb cocktail. These results support further clinical evaluation and indicate that ensovibep could be a valuable alternative to mAb cocktails and other treatments for COVID-19. Structural basis for Nirmatrelvir in vitro efficacy against SARS-CoV-2 variants The COVID-19 pandemic continues to be a public health threat with emerging variants of SARS-CoV-2. Nirmatrelvir (PF-07321332) is a reversible, covalent inhibitor targeting the main protease (Mpro) of SARS-CoV-2 and the active protease inhibitor in PAXLOVID (nirmatrelvir tablets and ritonavir tablets). We evaluated the in vitro catalytic activity and in vitro potency of nirmatrelvir against the main protease (Mpro) of prevalent variants of concern (VOC) or variants of interest (VOI): Alpha (, B.1.1.7), Beta ({beta}, B.1.351), Delta ({delta}, B1.617.2), Gamma ({gamma}, P.1), Lambda ({lambda}, B.1.1.1.37/C37), Omicron (o, B.1.1.529) as well as the original Washington or wildtype strain. These VOC/VOI carry prevalent mutations at varying frequencies in the Mpro specifically for: , {beta},{gamma} (K90R), {lambda} (G15S) and o (P132H). In vitro biochemical enzymatic assay characterization of the enzyme kinetics of the mutant Mpros demonstrate that they are catalytically comparable to wildtype. Nirmatrelvir has similar potency against each mutant Mpro including P132H that is observed in the Omicron variant with a Ki of 0.635 nM as compared to a Ki of 0.933nM for wildtype. The molecular basis for these observations were provided by solution-phase structural dynamics and structural determination of nirmatrelvir bound to the o, {lambda} and {beta} Mpro at 1.63 - 2.09 [A] resolution. These in vitro data suggest that PAXLOVID has the potential to maintain plasma concentrations of nirmatrelvir many-fold times higher than the amount required to stop the SARS-CoV-2 VOC/VOI, including Omicron, from replicating in cells (1).