G2 stem cells orchestrate time-directed, long-range coordination of calcium signaling during skin epidermal regeneration Skin epidermal homeostasis is maintained via constant regeneration by stem cells, which must communicate to balance their self-renewal and differentiation. A key molecular pathway, Ca2+ signaling has been implicated as a signal integrator in developing and wounded epithelial tissues[1, 2, 3, 4]. Yet how stem cells carry out this signaling across a regenerative tissue remains unknown due to significant challenges in studying signaling dynamics in live mice, limiting our understanding of the mechanisms of stem cell communication during homeostasis. To interpret high dimensional signals that have complex spatial and temporal patterns, we combined optimized imaging of Ca2+ signaling in thousands of epidermal stem cells in living mice with a new machine learning tool, Geometric Scattering Trajectory Homology (GSTH). Using a combination of signal processing, data geometry, and topology, GSTH captures patterns of signaling at multiple scales, either between direct or distant stem cell neighbors. Here we show that epidermal stem cells display dynamic intercellular Ca2+ signaling among neighborhoods of up to 10 cells that is surprisingly coordinated and directed through time across a pool of thousands of stem cells. We find that this collective coordination is an emergent property of the stem cell compartment, distinct from excitatory quiescent neuronal tissues. We demonstrate that cycling stem cells, specifically G2 cells, govern homeostatic patterns of Ca2+ signaling. Stem cells in different cell cycle stages dynamically regulate localization of the gap junction component Connexin43 (Cx43). Lastly, we uncouple global from local communication and identify Cx43 as the molecular mediator necessary for connectivity between local signaling neighborhoods. This work provides resolution in how stem cells at different stages of the cell cycle communicate and how that diversity of phases is essential for tissue wide communication and signaling flow during epidermal regeneration. Our approach provides a framework to investigate stem cell populations and their signaling dynamics, previously not possible.
Molecular mechanisms of stress-induced reactivation in mumps virus condensates Many viruses can establish long-term persistent infections in the human host and cause chronic diseases. Here we combine cell biology, whole-cell proteomics and cryo-electron tomography to uncover how cellular stress disrupts the host-virus equilibrium in persistent infection and induces viral replication, using a model of negative-stranded RNA viruses, the mumps virus. We show that phosphorylation of the largely disordered viral phosphoprotein coincides with increased partitioning of viral polymerase into pre-formed liquid-like condensates and the formation of a stable replication machinery. By obtaining the first atomic models for the authentic mumps virus nucleocapsid, we elucidate a concomitant conformational change that exposes the viral genome to its replication machinery. These events that occur within viral condensates upon stress, together with concerted down-regulation of the host antiviral defense, provide an environment that supports up-regulation of viral replication and constitute a stress-mediated switch that disrupts the host-virus equilibrium in persistent infection.

In BriefA multi-scale approach uncovers molecular and structural basis of how cellular stress provokes activation of persistent viral infection mediated by biomolecular condensates.
Balancing sensitivity and specificity in preclinical research The success of scientific discovery in preclinical drug development is based on the different roles of exploration and confirmation. Exploration requires sensitivity while confirmation needs to be specific. Based on published empirical data, we simulated two preclinical research trajectories to examine their effectiveness to detect true effects while keeping animal numbers low. We compare a conventional p-value based approach with a method based on an a priori determined smallest effect size of interest (SESOI). Using a SESOI increases transition rates from exploration to confirmation and results in higher detection rates across the trajectory. Specificity in the SESOI trajectory increases if the prior probability of true effects is low. We conclude that employing a SESOI is superior to a p-value based approach. Based on our findings, we propose a reconsideration of planning and conducting preclinical experiments, especially when the prior probability of true hypotheses is low. 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. SERCA function is impaired in skeletal and cardiac muscles from young DBA/2J mdx mice The DBA/2J (D2) mdx mouse has emerged as a more severe model of Duchenne muscular dystrophy when compared to the traditional C57BL/10 (C57) mdx mouse. Here, we questioned whether sarco(endo)plasmic reticulum Ca2+-ATPase (SERCA) function would differ in muscles from young D2 and C57 mdx mice. In gastrocnemius muscles, both D2- and C57 mdx mice exhibited signs of impaired Ca2+ uptake, however, this was more severe in D2 mdx mice. Maximal SERCA activity was lowered only in D2 mdx gastrocnemius muscles and not C57 mdx muscles. Furthermore, in the left ventricle and diaphragm, Ca2+ uptake was impaired in C57 mdx muscles with lowered rates of Ca2+ uptake compared with C57 WT mice, whereas in muscles from D2 mdx mice, rates of Ca2+ uptake were unattainable due to the severe impairments in their ability to transport Ca2+. Overall, our study demonstrates that SERCA function is drastically impaired in young D2 mdx mice. Autotaxin is a Potential Link between Genetic Risk Factors and Pathogenesis of Systematic Lupus Erythematosus in Plasmacytoid Dendritic Cells BackgroundThe importance of autotaxin, an enzyme that catalyzes lysophospholipid production, has recently been recognized in various diseases, including cancer and autoimmune diseases. Herein, we examined the role of autotaxin in systemic lupus erythematosus (SLE), utilizing data from ImmuNexUT, a comprehensive database consisting of transcriptome data and expression quantitative trait locus (eQTL) data of immune cells from patients with immune-mediated disorders.

MethodsSerum autotaxin concentrations in patients with SLE and healthy controls (HCs) were compared. The transcriptome data of patients with SLE and age-and sex-matched HCs were obtained from ImmuNexUT. The expression of ENPP2, the gene encoding autotaxin, was examined in peripheral blood immune cells. Next, weighted gene correlation network analysis (WGCNA) was performed to identify genes with expression patterns similar to ENPP2. The ImmuNexUT eQTL database and public epigenomic databases were used to infer the relationship between autotaxin and pathogenesis of SLE.

ResultsAutotaxin levels were elevated in the serum of patients with SLE compared to HCs. Furthermore, the expression of ENPP2 was higher in plasmacytoid dendritic cells (pDCs) than in other immune cell subsets, and its expression was elevated in pDCs of patients with SLE compared to HCs. In WGCNA, ENPP2 belonged to a module that correlated with disease activity. This module was enriched in interferon-associated genes and included genes whose expression was influenced by single-nucleotide polymorphisms associated with SLE, suggesting that it is a key module connecting genetic risk factors of SLE with disease pathogenesis. Analysis utilizing the ImmuNexUT eQTL database and public epigenomic databases suggested that the increased expression of ENPP2 in pDCs from patients with SLE may be caused by increased expression of interferon-associated genes and increased binding of STAT3 complexes to the regulatory region of ENPP2.

ConclusionsAutotaxin may play a critical role in connecting genetic risk factors of SLE to disease pathogenesis in pDCs. Neuronal Activation of the Gastrointestinal Tract Shapes the Gut Environment in Mice The enteric nervous system (ENS) coordinates responses that shape the intestinal milieu, and participates in bi-directional communication with various organs including the brain. Tools to study, and knowledge of, ENS development and function lag behind brain research. Herein, we deploy recombinant adeno-associated viral (rAAV) vectors with enhanced tropism for the gut to map and activate gut associated neurons in mice with spatial and temporal resolution. rAAV-mediated fluorescent labelling coupled with whole-tissue clearing methods in the small intestine and colon reveal novel neuronal architectures in the gut of mice. Viral delivery of genetically-encoded calcium indicators enable in vivo measurement of ENS activity. Further, we employed chemogentics to specifically activate gut neurons that express choline acetyltransferase (ChAT+) or tyrosine hydroxylase (TH+). Targeted activation of ChAT+ or TH+ neuronal populations associated with the gastrointestinal (GI) tract alters the intestinal transcriptome, mouse and microbial proteomes, metagenome, and fecal metabolome. We reveal previously unknown and broad roles for enteric neurons in modulating GI physiology, mucosal immunity, and microbiome structure, and provide evidence for novel functions by the ENS such as regulating fungal colonization and shaping of bile acid profiles in the gut. These advanced experimental tools and rich datasets with multi-parameter characterization of the gut ecosystem may facilitate expanded understanding and study of the ENS. Stimulus manipulations permit activation of fiber subpopulations in the mouse and rat vagus Afferent and efferent vagal fibers mediate bidirectional communication between the brain and visceral organs. Small, unmyelinated C-afferents constitute the majority of vagal fibers, play critical roles in numerous interoceptive circuits and autonomic reflexes in health and disease and may contribute to the efficacy and safety of vagus nerve stimulation (VNS). Selective engagement of C-afferents with electrical stimuli has not been feasible, due to the default fiber recruitment order: larger fibers first, smaller fibers last. Here, we determine and optimize an electrical stimulus that selectively engages vagal C-afferents. Intermittent KHz-frequency electrical stimulation (KES) activates motor and, preferentially, sensory vagal neurons in the brainstem. During KES, asynchronous activity of C-afferents increases, while that of larger fibers remains largely unchanged. In parallel, KES effectively blocks excitability of larger fibers while moderately suppressing excitability of C-afferents. By compiling selectivity indices in individual animals, we find that optimal KES parameters for C-afferents are >5KHz frequency and 7-10 times engagement threshold (xT) intensity in rats, 15-25xT in mice. These effects can be explained in computational models by how sodium channel responses to KES are shaped by axonal size and myelin. Our results indicate that selective engagement of vagal C-afferents is attainable by intermittent KES. Juvenile Gouldian Finches (Erythrura gouldiae) form sibling sub-groups during social integration The formation of social relationships in complex groups is critical in shaping patterns of social organization and behavioral development. In many birds, young individuals remain dependent on their parents for extended periods but must abruptly transition to navigating interactions in the wider group after independence. While lack of social relationships during this period is detrimental in the development of later social skills, little is known about the social relationships juveniles form after independence in many bird species. In this study we describe patterns of social interactions in juvenile Gouldian Finches (Erythrura gouldiae) during transitions from family groups into flocks of unfamiliar individuals. A group of 20 juveniles from 4 families were introduced into two flocks. After introductions juveniles showed a gradient of approach rates with most approaches directed towards siblings, followed by juvenile peers, adult females, and lastly adult males. Significant preferences for siblings resulted in the emergence of sibling subgroups within the larger social network. This active self-assortment of siblings suggests that sibling sub-groups are an important bridge linking social connections within the family and the wider group. Such findings suggests that sibling relationships have a critical role in the socialization after independence, as well as structuring the social organization of Gouldian finch flocks. Influence of DNA extraction kits on freshwater fungal DNA metabarcoding BackgroundEnvironmental DNA (eDNA) metabarcoding is a common technique for efficient biodiversity monitoring, especially of microbes. Recently, the usefulness of aquatic eDNA in monitoring the diversity of both terrestrial and aquatic fungi has been suggested. In eDNA studies, different experimental factors, such as DNA extraction kits or methods, can affect the subsequent analyses and the results of DNA metabarcoding. However, few methodological studies have been carried out on eDNA of fungi, and little is known about how experimental procedures can affect the results of biodiversity analysis. In this study, we focused on the effect of DNA extraction method on fungal DNA metabarcoding using freshwater samples obtained from rivers and lakes.

MethodsDNA was extracted from freshwater samples using the DNeasy PowerSoil kit, which is mainly used to extract microbial DNA from soil, and the DNeasy Blood & Tissue kit, which is commonly used for eDNA studies on animals. We then compared PCR inhibition and fungal DNA metabarcoding results; i.e., operational taxonomic unit (OTU) number and composition of the extracted samples.

ResultsNo PCR inhibition was detected in any of the samples, and no significant differences in the number of OTUs and OTU compositions were detected between the samples processed using different kits. These results indicate that both DNA extraction kits may provide similar diversity results for the river and lake samples evaluated in this study. Therefore, it may be possible to evaluate the diversity of fungi using a unified experimental method, even with samples obtained for diversity studies on other taxa such as those of animals. Morphology of immatures of the thelytokous ant, Monomorium triviale Wheeler (Formicidae: Myrmicinae: Solenopsidini) with descriptions of the extraordinary last-instar queen larvae The ant genus Monomorium is one of the most species-rich but taxonomically problematic groups in the hyperdiverse subfamily Myrmicinae. An East Asian species, M. triviale Wheeler, produces both reproductive queens and sterile workers via obligate thelytokous parthenogenesis. Here, we describe the immature forms of M. triviale based on light and scanning electron microscopy observations, with a note on the striking caste dimorphism in the last larval instar. The last-instar queen larvae were easily recognized by their large size, "aphaenogastroid" body shape, and rows of doorknob-like tubercles on the lateral and dorsal body surface. This type of queen-specific structure has not been found in ants in general, let alone congeneric species found in Japan. In stark contrast to the queen larvae, worker larvae showed a "pheidoloid" body shape and a body surface similar to other ants. The worker larvae were estimated to have three instars, consistent with previously described congeners. The pupae of both castes had no cocoon, a characteristic commonly described in other Myrmicinae species. In total, the developmental period from egg to adult worker averaged 59 days under 25{degrees}C. We discuss possible functions of the tubercles of queen larvae based on previous studies. Annotation of snoRNA abundance across human tissues reveals complex snoRNA-host gene relationships BackgroundSmall nucleolar RNAs (snoRNAs) are mid-size non-coding RNAs required for ribosomal RNA modification, implying a ubiquitous tissue distribution linked to ribosome synthesis. However, increasing numbers of studies identify extra-ribosomal roles of snoRNAs in modulating gene expression, suggesting more complex snoRNA abundance patterns. Therefore, there is a great need for mapping the snoRNome in different human tissues as the blueprint for snoRNA functions.

ResultsWe used a low structure bias RNA-Seq approach to accurately quantify snoRNAs and compare them to the entire transcriptome in seven healthy human tissues (breast, ovary, prostate, testis, skeletal muscle, liver and brain). We identified 475 expressed snoRNAs categorized in two abundance classes that differ significantly in their function, conservation level and correlation with their host gene: 390 snoRNAs are uniformly expressed and 85 are enriched in the brain or reproductive tissues. Most tissue-enriched snoRNAs are embedded in lncRNAs and display strong correlation of abundance with them, whereas uniformly expressed snoRNAs are mostly embedded in protein-coding host genes and are mainly non- or anticorrelated with them. 59% of the non-correlated or anticorrelated protein-coding host gene/snoRNA pairs feature dual-initiation promoters, as opposed to only 16% of the correlated non-coding host gene/snoRNA pairs.

ConclusionsOur results demonstrate that snoRNAs are not a single homogeneous group of housekeeping genes but include highly regulated tissue-enriched RNAs. Indeed, our work indicates that the architecture of snoRNA host genes varies to uncouple the host and snoRNA expressions in order to meet the different snoRNA abundance levels and functional needs of human tissues. Data-Driven Network Models for Genetic Circuits From Time-Series Data with Incomplete Measurements Synthetic gene networks are frequently conceptualized and visualized as static graphs. This view of biological programming stands in stark contrast to the transient nature of biomolecular interaction, which is frequently enacted by labile molecules that are often unmeasured. Thus, the network topology and dynamics of synthetic gene networks can be difficult to verify in vivo or in vitro, due to the presence of unmeasured biological states. Here we introduce the dynamical structure function as a new mesoscopic, data-driven class of models to describe gene networks with incomplete measurements. We introduce a network reconstruction algorithm and a code base for reconstructing the dynamical structure function from data, to enable discovery and visualization of graphical relationships in a genetic circuit diagram as time-dependent functions rather than static, unknown weights. We prove a theorem, showing that dynamical structure functions can provide a data-driven estimate of the size of crosstalk fluctuations from an idealized model. We illustrate this idea with numerical examples. Finally, we show how data-driven estimation of dynamical structure functions can explain failure modes in two experimentally implemented genetic circuits, a historical genetic circuit and a new E. coli based transcriptional event detector. Joint interactions with humans may pose a higher risk of zoonotic outbreaks than interactions with conspecifics among wildlife populations at human-wildlife interfaces O_LIPandemics caused by wildlife-origin pathogens, like COVID-19, highlight the importance of understanding the ecology of zoonotic transmission and outbreaks among wildlife populations at human-wildlife interfaces. To-date, the relative effects of human-wildlife and wildlife-wildlife interactions on the likelihood of such outbreaks remain unclear.
C_LIO_LIIn this study, we used social network analysis and epidemiological Susceptible Infected Recovered (SIR) models, to track zoonotic outbreaks through wild animals joint propensities to engage in social-ecological co-interactions with humans, and their social grooming interactions with conspecifics.
C_LIO_LIWe collected behavioral and demographic data on 10 groups of macaques (Macaca spp.) living in (peri)urban environments across Asia. Outbreak sizes predicted by the SIR models were related to structural features of the social networks, and particular properties of individual animals connectivity within those networks.
C_LIO_LIOutbreak sizes were larger when the first-infected animal was highly central, in both types of networks. Across host-species, particularly for rhesus and bonnet macaques, the effects of network centrality on outbreak sizes were stronger through macaques human co-interaction networks compared to grooming networks.
C_LIO_LIOur findings, independent of pathogen-transmissibility, suggest that wildlife populations in the Anthropocene are vulnerable to zoonosis more so due to their propensities to aggregate around anthropogenic factors than their gregariousness with conspecifics. Thus, the costs of zoonotic outbreaks may outweigh the potential/perceived benefits of jointly interacting with humans to procure anthropogenic food. From One Health perspectives, animals that consistently interact with both humans and conspecifics across time and space are useful targets for disease spillover assessments and control.
C_LI In Vivo Multi-Day Calcium Imaging of CA1 Hippocampus in Freely Moving Rats Reveals a High Preponderance of Place Cells with Consistent Place Fields Calcium imaging using GCaMP indicators and miniature microscopes has been used to image cellular populations during long timescales and in different task phases, as well as to determine neuronal circuit topology and organization. Because the hippocampus (HPC) is essential for tasks of memory, spatial navigation, and learning, calcium imaging of large populations of HPC neurons can provide new insight on cell changes over time during these tasks. All reported HPC in vivo calcium imaging experiments have been done in mouse. However, rats have many behavioral and physiological experimental advantages over mice. In this paper, we present the first (to our knowledge) in vivo calcium imaging from CA1 hippocampus in freely moving male rats. Using the UCLA Miniscope, we demonstrate that, in rat, hundreds of cells can be visualized and held across weeks. We show that calcium events in these cells are highly correlated with periods of movement, with few calcium events occurring during periods without movement. We additionally show that an extremely large percent of cells recorded during a navigational task are place cells (77.3{+/-}5.0%, surpassing the percent seen during mouse calcium imaging), and that these cells enable accurate decoding of animal position and can be held over days with consistent place fields in a consistent spatial map. A detailed protocol is included, and implications of these advancements on in vivo imaging and place field literature are discussed.

Significance statementIn vivo calcium imaging in freely moving animals allows the visualization of cellular activity across days. In this paper, we present the first in vivo Ca2+ recording from CA1 hippocampus in freely moving rats. We demonstrate that hundreds of cells can be visualized and held across weeks, and that calcium activity corresponds to periods of movement. We show that a high percentage (77.3{+/-}5.0%) of imaged cells are place cells, and that these place cells enable accurate decoding and can be held stably over days with little change in field location. Because the hippocampus is essential for many tasks involving memory, navigation, and learning, imaging of large populations of HPC neurons can shed new insight on cellular activity changes and organization. Shade-induced transcription of PIF-Direct-Target Genes precedes H3K4-trimethylation chromatin modification rises The phytochrome (phy)-PIF (Phytochrome Interacting Factor) sensory module perceives and transduces light signals to Direct-Target Genes (DTGs), which then drive the adaptational responses in plant growth and development, appropriate to the prevailing environment. These signals include the first exposure of etiolated seedlings to sunlight upon emergence from subterranean darkness, and the change in color of the light that is filtered through, or reflected from, neighboring vegetation ( shade). Previously, we identified three broad categories of rapidly signal-responsive genes: those repressed by light and conversely induced by shade; those repressed by light, but subsequently unresponsive to shade; and those responsive to shade only. Here, we investigate the potential role of epigenetic chromatin modifications in regulating these contrasting patterns of phy-PIF module-induced expression of DTGs. Using RNA-seq and ChlP-seq, time-resolved profiling of transcript and histone 3 lysine 4 trimethylation (H3K4me3) levels, respectively, we show that, whereas the initial dark-to-light transition triggers a rapid, apparently temporally-coincident decline of both parameters, the light-to-shade transition induces similarly rapid increases in transcript levels that precede increases in H3K4me3 levels. Together with other recent findings, these data raise the possibility that, rather than being causal in the shade-induced expression changes, H3K4me3 may function to buffer the rapidly fluctuating shade/light switching that is intrinsic to vegetational canopies under natural sunlight conditions. An integrative approach leads to the discovery of a novel anti-leukemic peptide from human milk Chemotherapy in childhood leukemia is associated with late morbidity in leukemic survivors, while certain patient subsets are relatively resistant to standard chemotherapy. It is therefore important to identify new agents with sensitivity and selectivity towards leukemic cells, while having less systemic toxicity. Peptide-based therapeutics has gained much attention during the last few years. Here, we used an integrative workflow combining mass spectrometric peptide library construction, in silico anticancer peptide screening, and in vitro leukemic cell studies to discover a novel anti-leukemic peptide having 3+charges and alpha-helical structure, namely HMP-S7, from human breast milk. HMP-S7 showed cytotoxic activity against four distinct leukemic cell lines in a dose-dependent manner but had no effect on solid malignancies or representative normal cells. HMP-S7 induced leukemic cell death by penetrating the plasma membrane to enter the cytoplasm and cause leakage of lactate dehydrogenase, thus acting in a membranolytic manner. Importantly, HMP-S7 exhibited anti-leukemic effect against patient-derived leukemic cells ex vivo. In conclusion, HMP-S7 is a selective anti-leukemic peptide with promise which requires further validation in preclinical and clinical studies.

TeaserIn silico screening of naturally occurring human milk peptides discovers a new anticancer peptide that kills leukemic cells in vitro and ex vivo. Nonideal polar residues in the core of the coiled coil of FtsLB are critical for modulating its stability and activation in cell division The FtsLB complex is a key regulator of bacterial cell division, existing in either an off or on state which supports the activation of septal peptidoglycan synthesis. In Escherichia coli, residues known to be critical for this activation are located in a region near the C-terminal end of the periplasmic coiled-coil domain of FtsLB, raising questions about the precise role of this conserved domain in the activation mechanism. Here, we investigate an unusual cluster of polar amino acids found within the core of the FtsLB coiled coil. We hypothesized that these amino acids likely reduce the structural stability of the domain and thus may be important for governing conformational changes. We found that mutating these positions to hydrophobic residues increased the thermal stability of FtsLB but caused cell division defects, suggesting that the coiled-coil domain is a "detuned" structural element. In addition, we identified suppressor mutations within the polar cluster, indicating that the precise identity of the polar amino acids is important for fine-tuning the structural balance between the off and on states. We propose a revised structural model of the tetrameric FtsLB (named the "Y-model") in which the periplasmic domain splits into a pair of coiled-coil branches. In this configuration, the hydrophilic terminal moieties of the polar amino acids remain more favorably exposed to water than in the original four-helix bundle model ("I-model"). We propose that a shift in this architecture, dependent on its marginal stability, is involved in activating the FtsLB complex and triggering septal cell wall reconstruction. A distinct brain expression profile for genes within large introgression deserts and under positive selection in Homo sapiens Analyses of ancient DNA from extinct hominins have provided unique insights into the complex evolutionary history of Homo sapiens, intricately related to that of the Neanderthals and the Denisovans as revealed by several instances of admixture events. These analyses have also allowed the identification of introgression deserts: genomic regions in our species that are depleted of  archaic haplotypes. The presence of genes like FOXP2 in these deserts has been taken to be suggestive of brain-related functional differences between Homo species. Here, we seek a deeper characterization of these regions, taking into account signals of positive selection in our lineage. Analyzing publicly available transcriptomic data from the human brain at different developmental stages, we found that structures outside the cerebral neocortex, in particular the cerebellum, the striatum and the mediodorsal nucleus of the thalamus show the most divergent transcriptomic profiles when considering genes within large introgression deserts and under positive selection. Elucidating Microscopic Events Driven by GTP Hydrolysis Reaction in Ras-GAP System with Semi-reactive Molecular Dynamics Simulation: Alternative Role of Phosphate Binding Loop as Mechanical Energy Storage ATPase and GTPase have been widely found as chemical energy-mechanical work transducers, whereas the physicochemical mechanisms are not satisfactorily understood. We addressed the problem by examining John Ross conjecture that repulsive Coulomb interaction between ADP/GDP and inorganic phosphate (Pi) does the mechanical work upon the system. We effectively simulated the consequence of GTP hydrolysis reaction in a complex system of Rat sarcoma (Ras) and GTPase activation protein (GAP) in the framework of classical molecular dynamics by switching force field parameters between the reactant and product systems. We then observed ca. 5 kcal/mol raise of potential energy about the phosphate-binding loop (P-loop) in Ras protein, indicating that the mechanical work generated via the GTP hydrolysis is converted into the local interaction energy and stored in the P-loop. Interestingly, this local energy storage in the P-loop depends on neither impulsive nor consecutive collisions of GDP and Pi with P-loop. Instead, GTP-GDP conversion itself does work on the Ras system, elevating the potential energy. These observations encourage us to challenge a conjecture previously given by Ross. We assert that triphosphate nucleotide hydrolyses do mechanical work by producing emergent steric interaction accompanied with relaxation, namely, a shift of biomolecular system to non-equilibrium state on the reshaped potential energy landscape. Recalling the universality of the P-loop motif among GTPases and ATPases, the observations that we obtained through this study would progress physicochemical understanding for the operating principles of GTP/ATP hydrolysis-driven biological nano-machines. A chemical screen based on an interruption of zebrafish gastrulation identifies the HTR2C inhibitor Pizotifen as a suppressor of EMT-mediated metastasis Metastasis is responsible for approximately 90% of cancer-associated mortality but few models exist that allow for rapid and effective screening of anti-metastasis drugs. Current mouse models of metastasis are too expensive and time consuming to use for rapid and high-throughput screening. Therefore, we created a unique screening concept utilizing conserved mechanisms between zebrafish gastrulation and cancer metastasis for identification of potential anti-metastatic drugs. We hypothesized that small chemicals that interrupt zebrafish gastrulation might also suppress metastatic progression of cancer cells and developed a phenotype-based chemical screen to test the hypothesis. The screen used epiboly, the first morphogenetic movement in gastrulation, as a marker and enabled 100 chemicals to be tested in five hours. The screen tested 1280 FDA-approved drugs and identified Pizotifen, an antagonist for serotonin receptor 2C (HTR2C) as an epiboly-interrupting drug. Pharmacologic and genetic inhibition of HTR2C suppressed metastatic progression in a mouse model. Blocking HTR2C with Pizotifen restored epithelial properties to metastatic cells through inhibition of Wnt-signaling. In contrast, HTR2C induced epithelial to mesenchymal transition (EMT) through activation of Wnt-signaling and promoted metastatic dissemination of human cancer cells in a zebrafish xenotransplantation model. Taken together, our concept offers a novel platform for discovery of anti-metastasis drugs. Transcriptional Control of Brain Tumour Stem Cell by a Carbohydrate Binding Protein Brain tumour stem cells (BTSCs) and intratumoural heterogeneity represent major challenges in glioblastoma therapy. Here, we report that the LGALS1 gene, encoding the carbohydrate binding protein, galectin1, is a key regulator of BTSCs and glioblastoma resistance to therapy. Genetic deletion of LGALS1 alters BTSC gene expression profiles and results in downregulation of gene sets associated with mesenchymal subtype of glioblastoma. Using a combination of pharmacological and genetic approaches, we establish that inhibition of LGALS1 signalling in BTSCs impairs self-renewal, suppresses tumourigenesis, prolongs lifespan, and improves glioblastoma response to ionizing radiation in preclinical animal models. Mechanistically, we show that LGALS1 is a direct transcriptional target of STAT3 with its expression robustly regulated by the ligand OSM. Importantly, we establish that galectin1 forms a complex with the transcription factor HOXA5 to reprogram BTSC transcriptional landscape. Our data unravel an oncogenic signalling pathway by which galectin1/HOXA5 complex maintains BTSCs and promotes glioblastoma. A novel Boolean network inference strategy to model early hematopoiesis aging Hematopoietic stem cell (HSC) aging is a multifactorial event that leads to changes in HSC properties and function. These changes are intrinsically coordinated and affect the early hematopoiesis, involving hematopoietic stem and progenitor cells (HSPCs). The objective of this work is to better understand the mechanisms and factors controlling these changes. We have therefore developed an original strategy to construct a Boolean network of genes explaining the priming and homeostasis of HSCs (graphical abstract). Based on our previous scRNA-seq data, we performed an exhaustive analysis of the transcriptional network and identified active transcription modules or regulons along the differentiation trajectory of selected HSPC states. This global view of transcriptional regulation led us to focus on 15 components, 13 selected TFs (Tal1, Fli1, Gata2, Gata1, Zfpm1, Egr1, Junb, Ikzf1, Myc, Cebpa, Bclaf1, Klf1, Spi1) and 2 complexes regulating the ability of HSC to cycle (CDK4/6 - Cyclin D and CIP/KIP). We then defined the connections controlling the differentiation dynamics of HSC states and constructed an influence graph between the TFs involved in the dynamics by mixing observations from our scRNA-seq data and knowledge from the literature. Then, using answer set programming (ASP) and in silico perturbation analysis, we obtained a Boolean model which is the solution of a Boolean satisfiability problem. Finally, perturbation of the model based on age-related changes revealed important regulations, such as the overactivation of Egr1 and Junb or the loss of Cebpa activation by Gata2, which were found to be relevant for the myeloid bias of aged HSC. Our work shows the efficiency of the combination of manual and systematic methods to elaborate a Boolean model. The developed strategy led to the proposal of new regulatory mechanisms underlying the differentiation bias of aged HSCs, explaining the decreased transcriptional priming of HSCs to all mature cell types except megakaryocytes. CDK9 and PP2A regulate RNA polymerase II transcription termination and coupled RNA maturation. CDK9 is a critical kinase required for the productive transcription of protein-coding genes by RNA polymerase II (pol II). As part of P-TEFb, CDK9 phosphorylates the carboxyl-terminal domain (CTD) of pol II and elongation factors, including SPT5, which allows pol II to elongate past the early elongation checkpoint (EEC) encountered soon after initiation. We show that, in addition to halting pol II at the EEC, loss of CDK9 activity causes premature termination of transcription across the last exon, loss of polyadenylation factors from chromatin, and loss of polyadenylation of nascent transcripts. Inhibition of the phosphatase PP2A abrogates the premature termination and loss of polyadenylation caused by CDK9 inhibition, indicating that this kinase/phosphatase pair regulates transcription elongation and RNA processing at the end of protein-coding genes. Our phosphoproteomic analyses after CDK9 inhibition, using either DRB or an analog-sensitive CDK9 cell line, confirm the splicing factor SF3B1 as an additional key target of this kinase. These results emphasize the important roles that CDK9 plays in coupling transcription elongation and termination to RNA maturation downstream of the EEC. hECA: the cell-centric assembly of a cell atlas Single-cell omics data can characterize multifaceted features of massive cells and bring significant insights to biomedical researches. The accumulation of single-cell data provides growing resources for constructing atlases for all cells of a human organ or the whole body. The true assembly of a cell atlas should be cell-centric rather than file-centric. We proposed a unified information framework enabling seamless cell-centric data assembly and developed a human Ensemble Cell Atlas (hECA) as an instance. hECA version 1.0 assembled scRNA-seq data across multiple studies into one orchestrated data repository. It contains 1,093,299 labeled cells and metadata from 116 published human single-cell studies, covering 38 human organs and 11 systems. We invented three methods of applications based on the cell-centric assembly: "In data" cell sorting enables targeted data retrieval in the full atlas with customizable logic expressions; The "quantitative portraiture" system provides a multi-view presentation of biological entities (organs, cell types, and genes) of multiple granularities; The customizable reference creation allows users to use the cell-centric assembly to generate references for their own cell type annotations. Case studies on agile construction of user-defined sub-atlases and "in data" investigation of CAR-T off-targets in multiple organs showed the great potential of cell-centric atlas assembly. Compression for population genetic data through finite-state entropy We improve the efficiency of population genetic file formats and GWAS computation by leveraging the distribution of sample ordering in population-level genetic data. We identify conditional exchangeability of these data, recommending finite state entropy algorithms as an arithmetic code naturally suited to population genetic data. We show between 10% and 40% speed and size improvements over dictionary compression methods for population genetic data such as Zstd and Zlib in computation and and decompression tasks. We provide a prototype for genome-wide association study with finite state entropy compression demonstrating significant space saving and speed comparable to the state-of-the-art.