Metacells untangle large and complex single-cell transcriptome networks The exponential scaling of scRNA-seq data represents an important hurdle for downstream analyses. Here we develop a coarse-graining framework where highly similar cells are merged into metacells. We demonstrate that metacells not only preserve but often improve the results of downstream analyses including visualization, clustering, differential expression, cell type annotation, gene correlation, imputation, RNA velocity and data integration. By capitalizing on the redundancy inherent to scRNA-seq data, metacells significantly facilitate and accelerate the construction and interpretation of single-cell atlases, as demonstrated by the integration of 1.46 million cells from COVID-19 patients in less than two hours on a standard desktop.
3CAC: improving the classification of phages and plasmids in metagenomic assemblies using assembly graphs MotivationBacteriophages and plasmids usually coexist with their host bacteria in microbial communities and play important roles in microbial evolution. Accurately identifying sequence contigs as phages, plasmids, and bacterial chromosomes in mixed metagenomic assemblies is critical for further unravelling their functions. Many classification tools have been developed for identifying either phages or plasmids in metagenomic assemblies. However, only two classifiers, PPR-Meta and viralVerify, were proposed to simultaneously identify phages and plasmids in mixed metagenomic assemblies. Due to the very high fraction of chromosome contigs in the assemblies, both tools achieve high precision in the classification of chromosomes but perform poorly in classifying phages and plasmids. Short contigs in these assemblies are often wrongly classified or classified as uncertain.

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

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

Contactrshamir@tau.ac.il

Supplementary informationSupplementary data are available at Bioinformatics online.
A Hierarchy of Biomolecular Proportional-Integral-Derivative Feedback Controllers for Robust Perfect Adaptation and Dynamic Performance Proportional-Integral-Derivative (PID) feedback controllers have been the most widely used controllers in industry for almost a century due to their good performance, simplicity, and ease of tuning. Motivated by their success in various engineering disciplines, PID controllers recently found their way into synthetic biology, where the design of feedback molecular control systems has been identified as an important goal. In this paper, we consider the mathematical realization of PID controllers via biomolecular interactions. We propose an array of topologies that offer a compromise between simplicity and high performance. We first demonstrate that different Proportional-Integral (PI) controllers exhibit different capabilities for enhancing the dynamics and reducing variance (cell-to-cell variability). Next, we introduce several derivative controllers that are realized based on incoherent feedforward loops acting in a feedback configuration. Alternatively, we show that differentiators can be realized by placing molecular integrators in a negative feedback loop--an arrangement that can then be augmented by PI components to yield PID feedback controllers. We demonstrate that the derivative component can be exploited for enhancing system stability, dramatically increasing the molecular control systems dynamic performance, and for reducing the noise effect on the output of interest. The PID controller features are established through various deterministic and stochastic analyses as well as numerical simulations. Finally, we provide an experimental demonstration using a recently developed hybrid setup, the cyberloop, where the controller is implemented in silico to control a biological genetic circuit in vivo. The large array of novel biomolecular PID controllers introduced here forms a basis for the design and construction of advanced high-performance biomolecular control systems that robustly regulate the dynamics of living systems. 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. Comparison of hybrid learning and remote education in the implementation of the 'Adopt a Microorganism' methodology The Internet has changed the way teachers and students access information and build knowledge. The recent COVID-19 pandemic has created challenges for both teachers and students and a demand for new methodologies of remote learning. In the life sciences, mixing online content with practical activities represents an even greater challenge. In microbiology, the implementation of an active teaching methodology, the #Adopt project, based on the social network Facebook(R), represents an excellent option for connecting remote education with classroom activities. In 2020, the version applied in high school, "Adopt a Microorganism", was adapted to meet the demands of emergency remote education owing to the suppression of face-to-face activities caused by the pandemic. In the present study, we assessed how the change in methodology impacted the discourse richness of students from high school integrated with technical education in the Business Administration program of the Federal Institute of Sao Paulo, Sorocaba Campus. Three questionnaires related to the groups of microorganisms (Archaea, Bacteria, Virus, Fungi, and Protozoan) were applied. The students responses in the 2019 and 2020 classes were compared concerning content richness and multiplicity of concepts through the application of the Shannon diversity index, an approach that is generally used to assess biodiversity in different environments. The observed results suggest that remote learning provided students with a conceptual basis and richness of content equivalent to that achieved by students subjected to the hybrid teaching model. In conclusion, this study suggests that the #Adopt project methodology increases students discourse richness in microbiology even without face-to-face traditional classes. Rostral Anterior Cingulate Activations inversely relate to Reward Payoff Maximation & predict Depressed Mood Choice selection strategies and decision making are typically investigated using multiple-choice gambling paradigms that require participants to maximize reward payoff. However, research shows that performance in such paradigms suffers from individual biases towards the frequency of gains to choose smaller local gains over larger longer term gain, also referred to as melioration. Here, we developed a simple two-choice reward task, implemented in 186 healthy human adult subjects across the adult lifespan to understand the behavioral, computational, and neural bases of payoff maximization versus melioration. The observed reward choice behavior on this task was best explained by a reinforcement learning model of differential future reward prediction. Simultaneously recorded and source-localized electroencephalography (EEG) showed that diminished theta-band activations in the right rostral anterior cingulate cortex (rACC) correspond to greater reward payoff maximization, specifically during the presentation of cumulative reward information at the end of each task trial. Notably, these activations (greater rACC theta) predicted depressed mood symptoms, thereby showcasing a reward processing marker of potential clinical utility.

Significance StatementThis study presents cognitive, computational and neural (EEG-based) analyses of a rapid reward-based decision-making task. The research has the following three highlights. 1) It teases apart two core aspects of reward processing, i.e. long term expected value maximization versus immediate gain frequency melioration based choice behavior. 2) It models reinforcement learning based behavioral differences between individuals showing that observed performance is best explained by differential extents of reward prediction. 3) It investigates neural correlates in 186 healthy human subjects across the adult lifespan, revealing specific theta band cortical source activations in right rostral anterior cingulate as correlates for maximization that further predict depressed mood across subjects. A Wnt/β-catenin signaling driven adipocyte population required for beiging and glucose homeostasis Wnt/{beta}-catenin signaling has been well established as a potent inhibitor of adipogenesis. Here, we identified a population of adipocytes that exhibit persistent activity of Wnt/{beta}-catenin signaling, as revealed by the Tcf/Lef-Gfp reporter allele, in embryonic and adult mouse fat depots, named as Wnt+ adipocytes. We showed that the {beta}-catenin-mediated signaling activation in these cells is Wnt ligand- and receptor-independent but relies on AKT/mTOR pathway and is essential for cell survival. Such adipocytes are distinct from classical ones in transcriptomic and genomic signatures and can be induced from various sources of mesenchymal stromal cells including human cells. Genetic lineage-tracing and targeted cell ablation studies revealed that these adipocytes convert into beige adipocytes directly and are also required for beige fat recruitment under thermal challenge, demonstrating both cell autonomous and non-cell autonomous roles in adaptive thermogenesis. Furthermore, mice bearing targeted ablation of these adipocytes exhibited glucose intolerance, while mice receiving exogenously supplied such cells manifested enhanced glucose utilization. Our studies uncover a unique adipocyte population in regulating beiging in adipose tissues and systemic glucose homeostasis. Allosteric activation of SARS-CoV-2 RdRp by remdesivir triphosphate and other phosphorylated nucleotides The catalytic subunit of SARS-CoV-2 RNA-dependent RNA polymerase (RdRp), Nsp12, has a unique NiRAN domain that transfers nucleoside monophosphates to the Nsp9 protein. The NiRAN and RdRp modules form a dynamic interface distant from their catalytic sites and both activities are essential for viral replication. We report that codon-optimized (for the pause-free translation) Nsp12 exists in inactive state in which NiRAN/RdRp interactions are broken, whereas translation by slow ribosomes and incubation with accessory Nsp7/8 subunits or NTPs partially rescue RdRp activity. Our data show that adenosine and remdesivir triphosphates promote synthesis of A-less RNAs, as does ppGpp, while amino acid substitutions at the NiRAN/RdRp interface augment activation, suggesting that ligand binding to the NiRAN catalytic site modulates RdRp activity. The existence of allosterically-linked nucleotidyl transferase sites that utilize the same substrates has important implications for understanding the mechanism of SARS-CoV-2 replication and design of its inhibitors.

HighlightsO_LICodon-optimization of Nsp12 triggers misfolding and activity loss
C_LIO_LISlow translation, accessory Nsp7 and Nsp8 subunits, and NTPs rescue Nsp12
C_LIO_LINon-substrate nucleotides activate RNA chain synthesis, likely via NiRAN domain
C_LIO_LICrosstalk between two Nsp12 active sites that bind the same ligands
C_LI Liberation from fibrogenesis or tumorigenesis via cellular senescence by a novel small molecule BackgroundUncontrolled proliferative diseases such as fibrosis or cancer can be fatal. We previously found that a compound containing a chromone scaffold, ONG41008, had potent anti-fibrotic effects in diseased human lung myofibroblasts but not normal human lung fibroblasts.

MethodsWe investigated the effects of ONG41008 on tumor cells, and compared these effects with those in pathogenic myofibrotic cells and normal fibroblasts cells.

FindingsStimulation of A549 lung carcinoma epithelial cells with ONG41008 resulted in cellular senescence, indicating that dysregulated cell proliferation is common to fibrotic cells and tumor cells. Replicative senescence of A549 cells resulted in multinucleation, which was followed by oncogene-induced senescence. There was significant upregulation of expression and nuclear translocation of p-TP53 and p16 in ONG41008-treated A549 cells, and all cells died after 72 hr. Similar effects occurred after ONG41008 treatment in several human aggressive cancer cell lines such as PANC1, MCF7, PC3, or primary non-small cell lung carcinoma cells. Unlike cisplatin, ONG41008 was not toxic to normal human lung fibroblasts or primary prostate epithelial cells, suggesting ONG41008 can distinguish the intracellular microenvironment between normal cells and aged or diseased cells. This effect might occur as a result of the increased NAD/NADH ratio or increased lactate dehydrogenase levels in aged or diseased cells.

InterpretationTo our best knowledge, this is the first study to show that a small molecule can arrest uncontrolled proliferation during fibrogenesis or tumorigenesis in vitro. ONG41008 could be a potential drug for a broad range of fibrotic or tumorigenic diseases.

Research in ContextO_ST_ABSEvidence before this studyC_ST_ABSThe notion that aging is a disease and that diseases occur as a consequence of aging was first put forward by David Sinclair and colleagues (Aging Cell;14(4):497-510). There are biological systems that provide evidence for this notion; for example, somatic cells can revert to embryonic cells, producing younger somatic cells. This phenomenon underlies induced pluripotent stem cells. Another example is that some types of jellyfish can live forever. These examples suggest that a counter-aging program exists in animals. Human diseases are the manifestations of cell aging generated by the accumulation of somatic mutations. Aged and pathogenic cells are senescent, so a drug that specifically targeted senescent cells might initiate a cellular program that could ameliorate age-associated disease. Indeed, the kinase inhibitor dasatinib induces cellular senescence (Clin Ther 2007 29:2289-2308). In 2017, two drugs that target senescent cells were identified: quercetin and fisetin. These drugs selectively kill senescent cells, and are referred to as senolytic drugs or senotherapeutics (Aging 2017 8;9(3):955-963). Although it is established that senescent cells accumulate in cancer and idiopathic pulmonary fibrosis (Nat Commun 2017 23;8 :14532), the effect of senolytic drugs in these diseases is largely unknown.

Added value of this studyThis study characterized a novel drug, termed ONG41008, which was found to have both senogenic and senolytic effects in cell-based assays. ONG41008 induced senescence in myofibroblasts and several cancer cell lines representative of aggressive human cancers, which was followed by cell death. Importantly, ONG41008 exhibited essentially no toxicity on normal human lung fibroblasts or primary prostate epithelial cells.

Implication of all the evidenceBased on our results, we believe that ONG41008 is a potent inducer of cellular senescence (replicative senescence and oncogene-induced senescence) and causes arrest of uncontrolled, pathogenic proliferation of myofibroblasts or cancer cells. Fatal neuroinvasion and SARS-CoV-2 tropism in K18-hACE2 mice is partially independent on hACE2 expression Animal models recapitulating distinctive features of severe COVID-19 are critical to enhance our understanding of SARS-CoV-2 pathogenesis. Transgenic mice expressing human angiotensin-converting enzyme 2 (hACE2) under the cytokeratin 18 promoter (K18-hACE2) represent a lethal model of SARS-CoV-2 infection. The precise mechanisms of lethality in this mouse model remain unclear. Here, we evaluated the spatiotemporal dynamics of SARS-CoV-2 infection for up to 14 days post-infection. Despite infection and moderate pneumonia, rapid clinical decline or death of mice was invariably associated with viral neuroinvasion and direct neuronal injury (including brain and spinal neurons). Neuroinvasion was observed as early as 4 dpi, with virus initially restricted to the olfactory bulb supporting axonal transport via the olfactory neuroepithelium as the earliest portal of entry. No evidence of viremia was detected suggesting neuroinvasion occurs independently of entry across the blood brain barrier. SARS-CoV-2 tropism was not restricted to ACE2-expressing cells (e.g., AT1 pneumocytes), and some ACE2-positive lineages were not associated with the presence of viral antigen (e.g., bronchiolar epithelium and brain capillaries). Detectable ACE2 expression was not observed in neurons, supporting overexpression of ACE2 in the nasal passages and neuroepithelium as more likely determinants of neuroinvasion in the K18-hACE2 model. Although our work incites caution in the utility of the K18-hACE2 model to study global aspects of SARS-CoV-2 pathogenesis, it underscores this model as a unique platform for exploring the mechanisms of SARS-CoV-2 neuropathogenesis that may have clinical relevance acknowledging the growing body of evidence that suggests COVID-19 may result in long-standing neurologic consequences.

IMPORTANCECOVID-19 is predominantly a respiratory disease caused by SARS-CoV-2 that has infected more than 191 million people with over 4 million fatalities (2021-07-20). The development of animal models recapitulating distinctive features of severe COVID-19 is critical to enhancing our understanding of SARS-CoV-2 pathogenesis and in the evaluation of vaccine and therapeutic efficacy. Transgenic mice expressing human angiotensin-converting enzyme 2 (hACE2) under the cytokeratin 18 promoter (K18-hACE2) represent a lethal model of SARS-CoV-2 infection. Here, we show lethality of this model is invariably associated with viral neuroinvasion linked with viral replication and assembly. Importantly, pneumonia albeit invariably present was generally moderate with the absence of culturable infectious virus at peak neuroinvasion. The dynamics of viral neuroinvasion and pneumonia were only partially dependent on hACE2. Overall, this study provides an in-depth sequential characterization of the K18-hACE2 model following SARS-CoV-2 infection, highlighting its significance to further study the mechanisms of SARS-CoV-2 neuropathogenesis. An early signaling GARP transcription factor regulates the entry point for differentiation in Giardia Transcriptional regulation of differentiation is critical for parasitic pathogens to adapt to environmental changes and regulate transmission. In response to encystation stimuli, Giardia lamblia cells shift from G1+S phase of the cell cycle to G2+M. By 2-4 hours, cyst wall proteins are upregulated, indicating that key regulatory steps occur within the first four hours of encystation. However, all characterized transcription factors (TFs) in Giardia have only been investigated at later time points of encystation. How TFs initiate encystation and link it to the cell cycle remains enigmatic. Here, we systematically screened six putative early upregulated TFs for nuclear localization, established their dynamic expression profiles, and determined their functional role in regulating encystation. We found a critical repressor, GLP4 that increases rapidly after 30 min of encystation stimuli and downregulates encystation specific markers, including Cyst Wall Proteins (CWPs) and enzymes in the cyst N-acetyl galactosamine (GalNAc) pathway. Depletion of GLP4 increases cyst production. Importantly, we observe that G2+M cells exhibit higher levels of CWP1 resulting from the activation of MYB2, a TF previously linked to encystation in Giardia. GLP4 upregulation occurs in G1+S cells, suggesting a role in repressing MYB2 and encystation specific genes in the G1+S phase of the cell cycle. Furthermore, we demonstrate that depletion of GLP4 upregulates MYB2 and promotes encystation while overexpression of GLP4 downregulates MYB2 and represses encystation. Together, these results suggest that Giardia employs a dose-dependent transcriptional response that involves the cell cycle regulated repressor GLP4 to orchestrate MYB2 and entry into the encystation pathway.

ImportanceTransition between life cycle stages is a common feature among parasitic pathogens and its regulation must be optimized to balance persistence of infection with transmission. The early transcription factors (TFs) regulating commitment to differentiate are totally unknown in Giardia. In this work, we identified GLP4, a previously uncharacterized GARP-like TF, as an early-acting transcriptional repressor that inhibits G1+S cells from entering the encystation pathway. GLP4 is therefore a key regulator controlling the balance between proliferative growth and terminal differentiation into infective cysts. Global diversity and distribution of prophages are lineage-specific within the Ralstonia solanacearum plant pathogenic bacterium species complex Ralstonia solanacearum is a destructive plant pathogenic bacterium and the causative agent of bacterial wilt disease, infecting over 200 plant species worldwide. In addition to chromosomal genes, its virulence is mediated by mobile genetic elements including integrated DNA of bacteriophages, i.e., prophages, which may carry fitness-associated auxiliary genes or modulate host gene expression. Although experimental studies have characterised several prophages that shape R. solanacearum virulence, the global diversity, distribution, and wider functional gene content of R. solanacearum prophages is unknown. In this study, prophages were identified in a diverse collection of 192 R. solanacearum draft genome assemblies originating from six continents. Prophages were identified bioinformatically and their diversity investigated using genetic distance measures, gene content, GC, and total length. Prophage distribution was characterised using metadata on R. solanacearum geographic origin and lineage classification (phylotypes), and their functional gene content was assessed by identifying putative prophage-encoded auxiliary genes. In total, 343 intact prophages were identified, forming ten genetically distinct clusters. These included five prophage clusters belonging to the Inoviridae, Myoviridae, and Siphoviridae phage families, and five uncharacterised clusters, possibly representing novel, previously undescribed phages. The prophages had broad geographical distribution being present across multiple continents. However, they were generally host phylogenetic lineage-specific, and overall, prophage diversity was proportional to the genetic diversity of their hosts. The prophages contained a myriad of auxiliary genes involved in metabolism and virulence of both phage and bacteria. Our results show that while R. solanacearum prophages are highly diverse globally, they make lineage-specific contributions to the R. solanacearum accessory genome, which could have resulted from shared coevolutionary history. An environmentally responsive transcriptional state modulates cell identities during root development Roots are fundamental organs for plant development and response to their environment: they anchor the plant to its growth substrate, uptake nutrients and water vital to plant growth, and can sense and respond to a variety of biotic and abiotic stresses. The architecture of root systems and their growth are known to be strongly affected by the environmental conditions found in the soil. However, the acquisition of cell identities at the root meristem is still mainly viewed as ontogenetically driven, where a small number of stem cells generate all the cell types through stereotyped divisions followed by differentiation, along a simple developmental trajectory. The extent to which environmental cues precisely shape and affect these developmental trajectories remains an open question. We used single-cell RNA-seq, combined with spatial mapping, to deeply explore the trajectories of cell states at the tip of Arabidopsis roots, known to contain multiple developing lineages. Surprisingly, we found that most lineage trajectories exhibit a stereotyped bifid topology with two developmental trajectories rather than one. The formation of one of the trajectories is driven by a strong and specific activation of genes involved in the responses to various environmental stimuli, that affects only of a subset of the cells in multiple cell types simultaneously, revealing another layer of patterning of cell identities in the root that is independent of cell ontogeny. We demonstrate the robustness of this environmentally responsive transcriptional state by showing that it is present in a mutant where cell type identities are greatly perturbed, as well as in different Arabidopsis ecotypes. We also show that the root can adapt the proportion of cells that acquire this particular state in response to environmental signals such as nutrient availability. The discovery of this cell state reveals new layers of cell identity that may underpin the adaptive potential of plant development. A novel nongenetic murine model of hyperglycemia and hyperlipidemia-associated aggravated atherosclerosis ObjectiveAtherosclerosis, the main pathology underlying cardiovascular diseases is accelerated in diabetic patients. Genetic mouse models require breeding efforts which are time-consuming and costly. Our aim was to establish a new nongenetic model of inducible metabolic risk factors that mimics hyperlipidemia, hyperglycemia, or both and allows the detection of phenotypic differences dependent on the metabolic stressor(s).

Methods and ResultsWild-type mice were injected with gain-of-function PCSK9D377Y (proprotein convertase subtilisin/kexin type 9) mutant adeno-associated viral particles (AAV) and streptozotocin and fed either a high-fat diet (HFD) for 12 or 20 weeks or a high-cholesterol/high-fat diet (Paigen diet, PD) for 8 weeks. To evaluate atherosclerosis, two different vascular sites (aortic sinus and the truncus of the brachiocephalic artery) were examined in the mice. Combined hyperlipidemic and hyperglycemic (HGHCi) mice fed a HFD or PD displayed characteristic features of aggravated atherosclerosis when compared to hyperlipidemia (HCi HFD or PD) mice alone. Atherosclerotic plaques of HGHCi HFD animals were larger, showed a less stable phenotype (measured by the increased necrotic core area, reduced fibrous cap thickness, and less -SMA-positive area) and had more inflammation (increased plasma IL-1{beta} level, aortic pro-inflammatory gene expression, and MOMA-2-positive cells in the BCA) after 20 weeks of HFD. Differences between the HGHCi and HCi HFD models were confirmed using RNA-seq analysis of aortic tissue, revealing that significantly more genes were dysregulated in mice with combined hyperlipidemia and hyperglycemia than in the hyperlipidemia-only group. The HGHCi-associated genes were related to pathways regulating inflammation (increased Cd68, iNos, and Tnfa expression) and extracellular matrix degradation (Adamts4 and Mmp14). When comparing HFD with PD, the PD aggravated atherosclerosis to a greater extent in mice and showed plaque formation after 8 weeks. Hyperlipidemic and hyperglycemic mice fed a PD (HGHCi PD) showed less collagen (Sirius red) and increased inflammation (CD68-positive cells) within aortic plaques than hyperlipidemic mice (HCi PD). HGHCi-PD mice represent a directly inducible hyperglycemic atherosclerosis model compared with HFD-fed mice, in which atherosclerosis is severe by 8 weeks.

ConclusionWe established a nongenetically inducible mouse model allowing comparative analyses of atherosclerosis in HCi and HGHCi conditions and its modification by diet, allowing analyses of multiple metabolic hits in mice. Effects of the social environment on movement-integrated habitat selection O_LIMovement links the distribution of habitats with the social environment of animals using those habitats; yet integrating movement, habitat selection, and socioecology remains an opportunity for further study.
C_LIO_LIHere, our objective was to disentangle the roles of habitat selection and social association as drivers of collective movement in a gregarious ungulate. To accomplish this objective, we (1) assessed whether socially familiar individuals form discrete social communities and whether social communities have high spatial, but not necessarily temporal, overlap; and (2) we modelled the relationship between collective movement and selection of foraging habitats using socially informed integrated step selection analysis.
C_LIO_LIWe used social network analysis to assign individuals to social communities and determine short and long-term social preference among individuals. Using integrated step selection functions (iSSF), we then modelled the effect of social processes, i.e., nearest neighbour distance and social preference, and movement behaviour on patterns of habitat selection.
C_LIO_LIBased on assignment of individuals to social communities and home range overlap analyses, individuals assorted into discrete social communities, and these communities had high spatial overlap. By unifying social network analysis with iSSF, we identified movement-dependent social association, where individuals foraged with more familiar individuals, but moved collectively with any between foraging patches.
C_LIO_LIOur study demonstrates that social behaviour and space use are inter-related based on spatial overlap of social communities and movement-dependent habitat selection. Movement, habitat selection, and social behaviour are linked in theory. Here, we put these concepts into practice to demonstrate that movement is the glue connecting individual habitat selection to the social environment.
C_LI Somatic uniparental disomy mitigates the most damaging EFL1 allele combination in Shwachman-Diamond syndrome Shwachman-Diamond syndrome (SDS; OMIM: #260400) is caused by variants in SBDS (Shwachman-Bodian-Diamond syndrome gene), which encodes a protein that plays an important role in ribosome assembly. Recent reports suggest that recessive variants in EFL1 are also responsible for SDS. However, the precise genetic mechanism that leads to EFL1-induced SDS remains incompletely understood. Here we present three unrelated Korean SDS patients that carry biallelic pathogenic variants in EFL1 with biased allele frequencies, resulting from a bone marrow-specific somatic uniparental disomy (UPD) in chromosome 15. The recombination events generated cells that were homozygous for the relatively milder variant, allowing for the evasion of catastrophic physiological consequences. Still, the milder EFL1 variant was solely able to impair 80S ribosome assembly and induce SDS features in cell line, zebrafish, and mouse models. The loss of EFL1 resulted in a pronounced inhibition of terminal oligo-pyrimidine element-containing ribosomal protein transcript 80S assembly. Therefore, we propose a more accurate pathogenesis mechanism of EFL1 dysfunction that eventually leads to aberrant translational control and ribosomopathy. Conserved organ-specific microbial assemblages in different populations of a terrestrial crab Brachyuran crabs originated in the oceans and evolved specific morphological and physiological adaptations to live in freshwater, intertidal and even terrestrial habitats but the role of a selection mechanism involving symbiotic microorganisms long these colonization processes are not known. In this work we investigated the associated microbiota of three populations of a terrestrial brachyuran crab, Chiromantes haematocheir, to find evidence of a conserved crab-specific microbiome unrelated to the population of origin and dissimilar from environmental microbial assemblages. Bacterial 16S rRNA gene and fungal ITS sequences were obtained from selected crab organs and environmental matrices to profile microbial communities. In spite of the presence of truly marine larval stages and no gregarious behaviour, favouring microbiota exchanges, we found common, organ-specific microbiota, associated to the gut and the gills of the crabs (with more than 15% of the genera detected specifically enriched only in one organ). Our results suggest an early establishment of a new common, stable microbiota in the transition from water to land. Dimerization of kringle 1 domain from hepatocyte growth factor/scatter factor provides a potent minimal MET receptor agonist. Hepatocyte Growth Factor/Scatter Factor (HGF/SF) and its cognate receptor MET play several essential roles in embryogenesis and regeneration in post-natal life of epithelial organs such as liver, kidney, lung, and pancreas, prompting a strong interest in harnessing HGF/SF-MET signalling for regeneration of epithelial organs after acute or chronic damage. The limited stability and tissue diffusion of native HGF/SF, however, which reflect the tightly controlled, local mechanism of action of the morphogen, have led to a major search for HGF/SF mimics for therapy. In this work, we describe the rational design, production and characterisation of K1K1, a novel minimal MET agonist consisting of two copies of the kringle 1 domain of HGF/SF placed in tandem. K1K1 is highly stable and displays biological activities equivalent or superior to native HGF/SF in a variety of in vitro assay systems and in a mouse model of liver disease. These data suggest that this engineered ligand may find wide applications in acute and chronic diseases of the liver and other epithelial organs dependent on MET activation. A Eukaryotic Circuit for Secrete-and-Sense Autonomy Cancers represent complex autonomous robust systems, displaying interconnectivity with feedback control. Autonomous growth is fueled by a cancer cells ability to  secrete-and-sense growth factors: a poorly understood phenomenon. Using an integrated systems and experimental approach, here we dissect the impact of a feedback-coupled GTPase circuit within the secretory pathway that imparts secretion-coupled autonomy. The circuit is assembled when the Ras-superfamily monomeric GTPase Arf1, and the heterotrimeric GTPase Giabg and their corresponding GAPs and GEFs are coupled by GIV/Girdin, a protein that is known to fuel aggressive traits in diverse cancers. One forward and two key negative feedback loops within the circuit create closed-loop control (CLC), allow the two GTPases to coregulate each other, and convert the expected switch-like behavior of Arf1-dependent secretion into an unexpected dose response alignment behavior of sensing and secretion. Such behavior translates into cell survival that is self-sustained by stimuli-proportionate secretion. Findings highlight how enhanced coupling of two biological switches in cancer cells is critical for multiscale feedback control to achieve secretion-coupled autonomy.

GRAPHIC ABSTRACT

O_FIG O_LINKSMALLFIG WIDTH=163 HEIGHT=200 SRC="FIGDIR/small/436048v3_ufig1.gif" ALT="Figure 1">
View larger version (51K):
org.highwire.dtl.DTLVardef@183ff5corg.highwire.dtl.DTLVardef@f21de7org.highwire.dtl.DTLVardef@1893683org.highwire.dtl.DTLVardef@159d966_HPS_FORMAT_FIGEXP  M_FIG C_FIG In BriefThis work delivers an experimentally validated dynamical systems model of the cooperativity between two distinct classes of biological switches in eukaryotic cells and reveals the basis of secretion-coupled cellular autonomy in cancer.

HIGHLIGHTSO_LIModeling and experimental approaches were used to dissect a coupled GTPase circuit
C_LIO_LICoupling enables closed loop feedback and mutual control of GTPases
C_LIO_LICoupling generates dose response alignment behavior of sensing and secretion.
C_LIO_LICoupling is critical for multiscale feedback control to achieve secretion-coupled autonomy.
C_LI Response of soil and vegetation in a warm-temperate Pine forest to intensive biomass harvests, phosphorus fertilisation, and wood ash application. Background and AimsConcerns about climate change and carbon economy have prompted the promotion of alternative energy sources, including forest-based bioenergy. An evaluation of the environmental consequences of intensive harvests (stumps and roots, and also branches and foliage) for energy wood supply, and use of wood-ash recycling as a compensatory practice, helps in the evaluation of the use of forest biomass for energy production.

MethodsWe made use of records from a split-plot experimental site crossing four different intensities of biomass harvesting (Stem-Only Harvest [SOH], Aboveground Additional Harvest [AAH], Belowground Additional Harvest [BAH], and Whole-Tree Harvest [WTH]) and three compensation methods (control [C], wood ash application [A] and phosphorus fertilisation [P]) to evaluate, in the 11-years-old stand (maritime pine: Pinus pinaster) that followed the biomass exports of the former stand, their effects on nutrient budgets, tree growth, soil fertility, chemical properties and soil carbon. This site is located in a forest on a poor soil, under a warm temperate climate (SW France).

Key resultsDespite their low additional biomass exports (+10% for AAH to +34% for WTH), the non-conventional harvest practices exported much higher quantities of nutrients than the conventional SOH technique (for example +145% for N and K in WTH). Consequently, these treatments had negative effects on the soil nutritive status. Additional biomass harvests impacted the soil organic matter content, with negative effects on P-organic, soil cation exchange capacity, exchangeable Ca, and most extractible nutrients. However, data suggested that tree growth and foliage nutrient content had not yet been significantly impacted by harvest treatments, whereas tree nutritional status was improved by P-fertiliser or wood ash. As expected, we observed a positive effect of wood ash application on soil pH and nutrient content but, like additional harvests, wood ash application decreased the pool of soil organic carbon (~10% of the initial stock with ~7% of N-total losses).

ConclusionsOverall, this factorial experiment showed that exporting more forest biomass due to the additional harvesting of tree canopies, stumps and roots had negative consequences on the ecosystem biogeochemistry. Additional harvests have aggravated the poverty of the already oligotrophic soil, and decreased the soil organic carbon content. Importantly, applying nutrients as fertiliser or wood ash did not compensate for all the negative impacts of biomass exports and the method of wood ash recycling in forests could even decrease the soil organic carbon. VIP plasma levels associate with survival in severe COVID-19 patients, correlating with protective effects in SARS-CoV-2-infected cells. Infection by SARS-CoV-2 may elicit uncontrolled and damaging inflammatory responses. Thus, it is critical to identify compounds able to inhibit virus replication and thwart the inflammatory reaction. Here, we show that the plasma levels of the immunoregulatory neuropeptide VIP are elevated in patients with severe COVID-19, correlating with reduced inflammatory mediators and with survival on those patients. In vitro, VIP and PACAP, highly similar neuropeptides, decreased the SARS-CoV-2 genome replication in human monocytes and viral production in lung epithelial cells, also reducing cell death. Both neuropeptides inhibited the production of proinflammatory mediators in lung epithelial cells and in monocytes. VIP and PACAP prevented in monocytes the SARS-CoV-2-induced activation of NF-kB and SREBP1 and SREBP2, transcriptions factors involved in proinflammatory reactions and lipid metabolism, respectively. They also promoted CREB activation, a transcription factor with antiapoptotic activity and negative regulator of NF-kB. Specific inhibition of NF-kB and SREBP1/2 reproduced the anti-inflammatory, antiviral and cell death protection effects of VIP and PACAP. Our results support further clinical investigations of these neuropeptides against COVID-19. Elucidation of SARS-CoV-2 budding mechanisms through molecular dynamics simulations of M and E protein complexes SARS-CoV-2 and other coronaviruses pose major threats to global health, yet computational efforts to understand them have largely overlooked the process of budding, a key part of the coronavirus life cycle. When expressed together, coronavirus M and E proteins are sufficient to facilitate budding into the ER-Golgi intermediate compartment (ERGIC). To help elucidate budding, we ran atomistic molecular dynamics (MD) simulations using the Feig laboratorys refined structural models of the SARS-CoV-2 M protein dimer and E protein pentamer. Our MD simulations consisted of M protein dimers and E protein pentamers in patches of membrane. By examining where these proteins induced membrane curvature in silico, we obtained insights around how the budding process may occur. Multiple M protein dimers acted together to induce global membrane curvature through protein-lipid interactions while E protein pentamers kept the membrane planar. These results could eventually help guide development of antiviral therapeutics which inhibit coronavirus budding. An Extended Motif in the SARS-CoV-2 Spike Modulates Binding and Release of Host Coatomer in Retrograde Trafficking {beta}-Coronaviruses such as SARS-CoV-2 hijack coatomer protein-I (COPI) for spike protein retrograde trafficking to the progeny assembly site in endoplasmic reticulum-Golgi intermediate compartment (ERGIC). However, limited residue-level details are available into how the spike interacts with COPI. Here we identify an extended COPI binding motif in the spike that encompasses the canonical K-x-H dibasic sequence. This motif demonstrates selectivity for COPI subunit. Guided by an in silico analysis of dibasic motifs in the human proteome, we employ mutagenesis and binding assays to show that the spike motif terminal residues are critical modulators of complex dissociation, which is essential for spike release in ERGIC. COPI residues critical for spike motif binding are elucidated by mutagenesis and crystallography and found to be conserved in the zoonotic reservoirs, bats, pangolins, camels, and in humans. Collectively, our investigation on the spike motif identifies key COPI binding determinants with implications for retrograde trafficking. Genome annotation of Caenorhabditis briggsae by TEC-RED identifies new exons, paralogs, and conserved and novel operons The nematode Caenorhabditis briggsae is routinely used in comparative and evolutionary studies involving its well-known cousin C. elegans. The C. briggsae genome sequence has accelerated research by facilitating the generation of new resources, tools, and functional studies of genes. While substantial progress has been made in predicting genes and start sites, experimental evidence is still lacking in many cases. Here, we report an improved annotation of the C. briggsae genome using the Trans-spliced Exon Coupled RNA End Determination (TEC-RED) technique. In addition to identifying the 5 ends of expressed genes, we have discovered operons and paralogs. In summary, our analysis yielded 10,243 unique 5 end sequence tags with matches in the C. briggsae genome. Of these, 6,395 were found to represent 4,252 unique genes along with 362 paralogs and 52 previously unknown exons. These genes included 14 that are exclusively trans-spliced in C. briggsae when compared with C. elegans orthologs. A major contribution of this study is the identification of 493 operons, of which two-thirds are fully supported by tags. In addition, two SL1-type operons were discovered. Interestingly, comparisons with C. elegans showed that only 40% of operons are conserved. Of the remaining operons, 73 are novel, including 12 that entirely lack orthologs in C. elegans. Further analysis revealed that four of the 12 novel operons are conserved in C. nigoni. Altogether, the work described here has significantly advanced our understanding of the C. briggsae system and serves as a rich resource to aid biological studies involving this species. The first sheep graph-based pan-genome reveals the spectrum of structural variations and their effects on tail phenotypes Structural variations (SVs) are a major contributor to genetic diversity and phenotypic variations, but their prevalence and functions in domestic animals are largely unexplored. Here, we assembled 26 haplotype-resolved genome assemblies from 13 genetically diverse sheep using PacBio HiFi sequencing. We constructed a graph-based ovine pan-genome and discovered 142,422 biallelic insertions and deletions, 7,028 divergent alleles and 13,419 multiallelic variations. We then used a graph-based approach to genotype the biallelic SVs in 684 individuals from 45 domestic breeds and two wild species. Integration with RNA-seq data allows to identify candidate expression-associated SVs. We demonstrate a direct link of SVs and phenotypes by localizing the putative causative insertion in HOXB13 gene responsible for the long-tail trait and identifying multiple large SVs associated with the fat-tail. Beyond generating a benchmark resource for ovine structural variants, our study highlights that animal genetic research will greatly benefit from using a pan-genome graph rather than a single reference genome. Mitochondrial pyruvate supports lymphoma proliferation by fueling a non-canonical glutamine metabolism pathway The fate of pyruvate, which is modulated by the activity of the mitochondrial pyruvate carrier (MPC), is a defining metabolic feature in many cancers. Diffuse large B-cell lymphomas (DLBCLs) are a genetically and metabolically heterogeneous cancer. Although MPC expression and activity differed between DLBCL subgroups, mitochondrial pyruvate was uniformly consumed by glutamate pyruvate transaminase 2 (GPT2) to support -ketoglutarate production as part of glutaminolysis. This led us to discover that glutamine exceeds pyruvate as a carbon source for the tricarboxylic acid (TCA) cycle in DLBCLs. Furthermore, we found that MPC inhibition unexpectedly leads to decreased glutaminolysis, which is contrary to previous observations in other cell types. We also discovered that MPC inhibition and depletion only decreased DLBCL proliferation in an extracellular matrix (ECM) environment and in vivo xenografts, but not in the typical DLBCL suspension environment. We also have found that the metabolic profile of DLBCL cells in ECM is markedly different from cells in suspension environment. Thus, we report that besides the canonical glutamate dehydrogenase (GDH)-mediated glutaminolysis, the non-canonical GPT2 mediated consumption and assimilation of glutamine and pyruvate in DLBCLs enables their proliferation in an extracellular environment-dependent manner.

HIGHLIGHTSO_LIGlutamine, but not glucose, is a major carbon source for the tricarboxylic acid cycle in diffuse large B-cell lymphomas.
C_LIO_LIMitochondrial pyruvate supports glutaminolysis in diffuse large B-cell lymphomas by supplying pyruvate for glutamate pyruvate transaminase 2-mediated  -ketoglutarate production.
C_LIO_LIMitochondrial pyruvate carrier inhibition leads to decreased glutaminolysis in diffuse large B-cell lymphomas.
C_LIO_LI -ketoglutarate production is important for diffuse large B-cell lymphoma proliferation in a solid extracellular matrix environment.
C_LIO_LIMitochondrial pyruvate carrier activity supports diffuse large B-cell lymphoma proliferation in a solid extracellular matrix environment and in mouse xenografts.
C_LI