DeepVelo: Data-driven Single-cell Transcriptomic Deep Velocity Field Learning with Neural Ordinary Differential Equations Recent advances in single-cell RNA sequencing technology have provided unprecedented opportunities to simultaneously measure the gene expression profile and transcriptional velocity of individual cells, enabling us to sample gene regulatory network dynamics along developmental trajectories. However, traditional methods have faced challenges in modeling gene expression dynamics within individual cells due to sparse, non-linear (e.g., obligate heterodimer transcription factors), and high-dimensional measurements. Here, we present DeepVelo, a neural-network-based ordinary differential equation model that can learn non-linear, high-dimensional single-cell transcriptome dynamics and describe continuous gene expression changes within individual cells across time. We applied DeepVelo to multiple published datasets from different technical platforms and demonstrated its utility to 1) formulate transcriptome dynamics on different timescales, 2) measure the instability of cell states, and 3) identify developmental driver genes upstream of a signaling cascade. Benchmarking against state-of-the-art methods shows that DeepVelo can improve velocity field representation accuracy by at least 50% in out-of-sample cells. Further, perturbation studies revealed that single-cell dynamical systems may exhibit properties similar to those of chaotic systems. In summary, DeepVelo allows for the data-driven discovery of differential equations that delineate single-cell transcriptome dynamics.

TeaserEmbedding neural networks into ordinary differential equations to model gene expression changes within single cells across time.
Discovering Rare Cell Types through Information-based Dimensionality Reduction Dimensionality reduction is crucial to summarizing the complex transcriptomic landscape of single cell datasets for downstream analyses. However, current dimensionality reduction approaches favor large cellular populations defined by many genes, at the expense of smaller and more subtly-defined populations. Here, we present surprisal component analysis (SCA), a technique that leverages the information-theoretic notion of surprisal for dimensionality reduction, and demonstrate its ability to improve the representation of clinically important populations that are indistinguishable using existing pipelines. For example, in cytotoxic T-cell data, SCA cleanly separates the gamma-delta and MAIT cell subpopulations, which are not detectable via PCA, ICA, scVI, or a wide array of specialized rare cell recovery tools. We also show that, when used instead of PCA, SCA improves downstream imputation to more accurately restore mRNA dropouts and recover important gene-gene relationships. SCAs information-theoretic paradigm opens the door to more meaningful signal extraction, with broad applications to the study of complex biological tissues in health and disease.
Clonal Dynamics of Cancer-Initiating Cells in a Mouse Genetic Mosaic Model of Serous Ovarian Cancer High-grade serous ovarian cancer (HGSOC) is the most lethal gynecologic malignancy. Early detection remains the most promising strategy to decrease mortality and requires a deep understanding of the cancer initiation process. Recent breakthroughs identified the fallopian tube (FT) Pax8+ cells as a cell-of-origin for HGSOC. However, the progressive trajectory from an initiated Pax8+ cell toward malignancy remains elusive. Here, we use a mouse genetic mosaic system called MADM (Mosaic Analysis with Double Markers), which generates sporadic mutant cells labeled with GFP, to trace the clonal dynamics of HGSOC development. Surprisingly, we observed dichotomous clonal expansion kinetics of individual cancer-initiating Pax8+ cells; while the vast majority stalled immediately, a few significantly expanded. Clonal analysis of wildtype cells showed concordant dichotomy, suggesting that a small fraction of Pax8+ cells intrinsically carry higher expansion potential. Furthermore, although all wildtype clones stalled shortly after development, some mutant clones continued to progress, indicating that oncogenic mutations prolong clonal expansion beyond the cell-intrinsic potential. Finally, intra-clonal fate mapping revealed that progressing mutant clones have a higher propensity to produce Pax8+ cells, linking imbalanced cell differentiation to disease progression. In summary, our studies reveal that the earliest clonal expansion during HGSOC initiation emerges when oncogenic mutations exaggerate the intrinsic expansion potential of a rare Pax8+ subset and bias its differentiation potential, warranting future investigation as a readout for early detection and cancer prevention.

SIGNIFICANCE STATEMENTEffective early detection is the most promising approach to manage cancer, but it requires a thorough understanding of cancer initiation and early progression. Serous ovarian cancer has been found to originate from Pax8+ cells in the fallopian tube. To assess the clonal dynamics of cancer-initiating cells, we established a mouse genetic mosaic model that enables clonal analysis of mutant Pax8+ cells in the fallopian tube. We uncovered that the earliest clonal event in ovarian cancer development emerges from a rare subset of stem/progenitor-like cells with intrinsic expansion potential, and that oncogenic mutations exaggerate and prolong this potential to promote cancer initiation. Our findings pinpoint a precise target for future early-detection and cancer prevention studies.

HIGHLIGHTSO_LIOvarian cancer initiation revealed by a mouse genetic mosaic model.
C_LIO_LIDichotomous expansion capacity of initiated cells at the clonal level.
C_LIO_LISynergies between cell-intrinsic properties and onco-mutations fuels clonal expansion.
C_LIO_LIImbalanced differentiation within highly expansive clones leads to Pax8 clusters.
C_LI 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. Trypanosomal variant surface glycoprotein expression in human African trypanosomiasis patients Trypanosoma brucei gambiense, an extracellular protozoan parasite, is the primary causative agent of human African trypanosomiasis. T. b. gambiense is endemic to West and Central Africa, where it is transmitted by the bite of infected tsetse flies. In the bloodstream of an infected host, the parasite evades antibody recognition by altering the Variant Surface Glycoprotein (VSG) that forms a dense coat on its cell surface through a process known as antigenic variation. Each VSG has a variable N-terminal domain exposed to the host and a less variable C-terminal domain that is partially hidden from host antibodies. Our lab developed VSG-seq, a targeted RNA-seq method, to study VSG gene expression in T. brucei. Studies using VSG-seq to characterize antigenic variation in a mouse model have revealed marked diversity in VSG expression within parasite populations, but this finding has not yet been validated in a natural human infection. Here, we used VSG-seq to analyze VSGs expressed in the blood of 12 patients infected with T. b. gambiense. The number of VSGs identified per patient ranged from one to 14; notably, two VSGs were shared by more than one patient. Analysis of expressed VSG N-terminal domain types revealed that 82% of expressed VSGs encoded a type B N-terminus, a bias not seen in datasets from other T. brucei subspecies. C-terminal types in T. b. gambiense infection were also restricted. These results demonstrate a bias either in the underlying VSG repertoire of T. b. gambiense or in this subspecies expression of VSGs during infection. This work demonstrates the feasibility of using VSG-seq to study antigenic variation in human infections and highlights the importance of understanding VSG repertoires in the field.

Author SummaryHuman African Trypanosomiasis is a neglected tropical disease primarily caused by the extracellular parasite Trypanosoma brucei gambiense. To avoid elimination by the host, these parasites repeatedly replace their dense surface coat of Variant Surface Glycoprotein (VSG). Despite the important role of VSGs in prolonging infection, VSG expression during natural human infections is poorly understood. A better understanding of natural VSG gene expression dynamics can clarify the mechanisms that T. brucei uses to alter its VSG coat and improve how trypanosomiasis is diagnosed in humans. We analyzed the expressed VSGs detected in the blood of patients with trypanosomiasis. Our findings indicate that parasites in both natural and experimental infections express diverse sets of VSGs. "Adopt a Microorganism" methodology enhances the discourse richness of high school students in both hybrid learning and remote education 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. Markov Chain Abstractions of Electrochemical Reaction-Diffusion in Synaptic Transmission for Neuromorphic Computing Progress in computational neuroscience towards understanding brain function is challenged both by the complexity of molecular-scale electrochemical interactions at the level of individual neurons and synapses and the dimensionality of network dynamics across the brain covering a vast range of spatial and temporal scales. Our work abstracts an existing highly detailed, biophysically realistic 3D reaction-diffusion model of a chemical synapse to a compact internal state space representation that maps onto parallel neuromorphic hardware for efficient emulation at a very large scale and offers near-equivalence in input-output dynamics while preserving biologically interpretable tunable parameters. 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. Population bottlenecks constrain host microbiome diversity and genetic variation impeding fitness It is becoming increasingly clear that microbial symbionts influence key aspects of their hosts fitness, and vice versa. This may fundamentally change our thinking about how microbes and hosts interact in influencing fitness and adaptation to changing environments. Here we explore how reductions in population size commonly experienced by threatened species influence microbiome diversity. Consequences of such reductions are normally interpreted in terms of a loss of genetic variation, increased inbreeding and associated inbreeding depression. However, fitness effects of population bottlenecks might also be mediated through microbiome diversity, such as through loss of functionally important microbes. Here we utilise 50 Drosophila melanogaster lines with different histories of population bottlenecks to explore these questions. The lines were phenotyped for egg-to-adult viability and their genomes sequenced to estimate genetic variation. The bacterial 16S rRNA gene was amplified in these lines to investigate microbial diversity. We found that 1) host population bottlenecks constrained microbiome richness and diversity, 2) core microbiomes of hosts with low genetic variation were constituted from subsets of microbiomes found in flies with higher genetic variation, 3) both microbiome diversity and host genetic variation contributed to host population fitness, 4) connectivity and robustness of bacterial networks was low in the inbred lines regardless of host genetic variation, 5) reduced microbial diversity was associated with weaker evolutionary responses of hosts in stressful environments, and 6) these effects were unrelated to Wolbachia density. These findings suggest that population bottlenecks reduce hologenomic variation (combined host and microbial genetic variation). Thus, while the current biodiversity crisis focuses on population sizes and genetic variation of eukaryotes, an additional focal point should be the microbial diversity carried by the eukaryotes, which in turn may influence host fitness and adaptability with consequences for the persistence of populations.

Author summaryIt is becoming increasingly clear that organisms and the microbes that live on or in them - their microbiome - affect each other in profound ways that we are just beginning to understand. For instance, a diverse microbiome can help maintain metabolic functions or fight pathogens causing diseases. A disrupted microbiome may be especially critical for animals and plants that occur in low numbers because of threats from e.g. human exploitation or climate change, as they may already suffer from genetic challenges such as inbreeding and reduced evolutionary potential. The importance of such a reduction in population size, called a bottleneck, on the microbial diversity and the potential interactive effects on host health remains unexplored. Here we experimentally test these associations by investigating the microbiomes of 50 inbred or non-inbred populations of vinegar flies. We found that restricting the population size constrain the hosts genetic variation and simultaneously decreases the diversity of the microbiome that they harbor, and that both effects were detrimental to host fitness. The microbial communities in inbred host populations were less robust than in their non-inbred counterparts, suggesting that we should increasingly consider the microbiome diversity, which may ultimately influence the health and persistence of threatened species. Modulation of P2X4 pore closure by magnesium,potassium, and ATP The P2X4 receptor plays a prominent role in cellular responses to extracellular ATP. Through classical all-atom molecular dynamics (MD) simulations totaling 24 {micro}s we have investigated how metal-complexed ATP stabilizes the channels open state and prevents its closing. We have identified two metal-binding sites, magnesium (Mg2+) and potassium (K+), one at the intersection of the three subunits in the ectodomain (MBS1) and the second one near the ATP binding site (MBS2), similar to those characterized in Gulf coast P2X. Our data indicate that when Mg2+ and K+ ions are complexed with ATP, the channel is locked into an open state. Interestingly, irrespective of the number of bound ATP molecules, Mg2+ ions bound to the MBS2 resisted collapsing of the open state protein to a closed state by stabilizing the ATP-protein interactions. However, when Mg2+ in the MBS2 was replaced with K+ ions, as might be expected when in equilibrium with an extracellular solution, the interactions between the subunits were weakened and we found evidence of pore collapse. This collapse was apparent when fewer than two ATP were bound to MBS2 in the presence of K+. Therefore, the different capacities of common cations to stabilize the channel may underlie a mechanism governing P2X4 channel gating in physiological systems. This study provides structural insights into the differential modulation of ATP activation of P2X4 by Mg2+ and K+. A northern range extension of a Canadian species of Special Concern, Dielis pilipes (Hymenoptera: Scoliidae), in the Okanagan Valley of British Columbia The only known Canadian records of the yellow scarab hunter wasp, Dielis pilipes (Saussure) (Hymenoptera: Scoliidae), are from the southern Okanagan and Similkameen valleys of British Columbia. We report a 25-km northern range extension of the species into the Ponderosa Pine Biogeoclimatic Ecosystem Classification zone, collected in an unmanaged agricultural field in Summerland, British Columbia. This finding is of conservation importance and has implications for natural biological control of ten-lined June beetles, Polyphylla decemlineata (Say) and P. crinita LeConte (Coleoptera: Scarabaeidae), incidental agricultural pests in the Okanagan. Plucking a string or playing a G? Choice type impacts human reinforcement learning In reinforcement learning (RL) experiments, participants learn to make rewarding choices in response to different stimuli; RL models use outcomes to estimate stimulus-response values which change incrementally. RL models consider any response type indiscriminately, ranging from more concretely defined motor choices (e.g. pressing a key with the index finger), to more general choices that can be executed in a number of ways (e.g. getting dinner at the restaurant). But does the learning process vary as a function of the choice type? In Experiment 1, we show that it does: participants were slower and less accurate in learning correct choices of a general format compared to learning more concrete, motor actions. Using computational modeling, we show that two mechanisms contribute to this. First, there was evidence of irrelevant credit assignment: the values of motor actions interfered with the values of other choice dimensions, resulting in more incorrect choices when the correct response is not defined by a single motor action; second, information integration for relevant general choices was slower. In Experiment 2, we replicated and further extended the findings from Experiment 1, by showing that slowed learning was attributable to weaker working memory use, rather than slowed RL learning. In both experiments we ruled out the explanation that the difference in performance between two condition types was driven by difficulty/different levels of complexity. We conclude that defining a more abstract choice space used by multiple learning systems for credit assignment recruits executive resources, limiting how much such processes then contribute to fast learning. Atrazine induced In vivo immunotoxicity in Bivalves Atrazine is ubiquitously used broad-spectrum herbicide to control the weeds in agriculture. The present study aimed to evaluate the acute toxicity and immunotoxicity of Atrazine in two ecologically and economically important bivalves. Acute toxicity of atrazine evaluated in triplicates by taking control and six experimental groups each comprising of 30 animals and treated with a range of atrazine from 2 PPM to 12 PPM for 96 hours. Mortalities were recorded for every 24 hours until 96 hours and data analyzed by one-way ANOVA and Dunnett T-test. The results indicated a significant increase in mortalities with increase in dose and time of exposure in both species. The values of LC50 were determined as 6.10 PPM and 4.90 PPM respectively for Perna viridis and Paphia malabarica. Furthermore, the immunotoxic potential of atrazine assessed by treating mussels and clams with the five sub-lethal doses of atrazine for 14 days and quantifying the viability of hemocytes by using simple yet reliable Tryphan blue exclusion assay. The results of the present study suggest atrazine-induced immunotoxicity by decreasing the number of viable hemocytes in bivalves. Hemocytes with phagocytic function are indispensable to confer innate immunity in bivalves, decreased viability of these cells leads to compromised immunity. This study is first of its kind to implicate atrazine with the immunotoxicity in bivalves. Unifying Transcriptome Regulation and Cellular Operational Efficiency: Application of the Locality and Caching Principles via a Cell-to-Computer Analogy Gene expression is time-consuming, and the delay from transcription activation to produced proteins is sequentially longer from bacteria to yeast and to humans. How human cells overcome the delay and attain operational efficiency, i.e., quick proteomic response to signals, is not well understood. The computer has endured the same system latency issue due to much slower information retrieval (hard drive (HD) to memory and to central processing unit (CPU)) than CPU execution, and mitigated it via efficient memory management, namely, the spatiotemporal locality principles that control specialized user programs and the permanent caching of core system functions - the operating system (OS) kernel. In this study, we unified gene expression and HD-memory-CPU information flow as instances of the Shannon information theory, that support the respective system operations and both consist of three components: information storage, the execution/decoding step, and the channel for the dynamic storage-to-execution information flow; the gene expression machinery and their regulators, and the OS kernel, were deemed as the respective channels. This abstraction prompted a multi-omic comparative analysis, leading to three observations. First, the temporal locality principle explains the mRNA stabilization-by-translation regulatory mechanism and controls specialized cellular functions. Second, the caching principle explains cytoplasmic mRNA sequestration and the defiance of the locality principle by highly sequestered mRNAs. Third, strikingly, in both systems, the caching principle controls the information channels; similar to permanent caching of OS kernel, basic gene expression machinery and their regulators are the top most sequestered mRNAs. In short, the locality and the caching principles act in the cells to control specialized and core cellular functions, respectively. The results close the knowledge gap between, and thus unify, transcriptome regulation and cellular operational latency mitigation. Genome-wide association study of musical beat synchronization demonstrates high polygenicity Moving in synchrony to the beat is a fundamental component of musicality. Here, we conducted a genome-wide association study (GWAS) to identify common genetic variants associated with beat synchronization in 606,825 individuals. Beat synchronization exhibited a highly polygenic architecture, with sixty-nine loci reaching genome-wide significance (p<5x10-8) and SNP-based heritability (on the liability scale) of 13%-16%. Heritability was enriched for genes expressed in brain tissues, and for fetal and adult brain-specific gene regulatory elements, underscoring the role of central nervous system-expressed genes linked to the genetic basis of the trait. We performed validations of the self-report phenotype (through internet-based experiments) and of the GWAS (polygenic scores for beat synchronization were associated with patients algorithmically classified as musicians in medical records of a separate biobank). Genetic correlations with breathing function, motor function, processing speed, and chronotype suggest shared genetic architecture with beat synchronization and provide avenues for new phenotypic and genetic explorations. First Report of Lactoria cornuta from Karnataka, West coast of India, Role of tropical cyclones in Range Extension India is known for its rich biodiversity and is fortunate to have several endemic species from different classes of vertebrates. India is home to 7.5% of the global fish diversity, with 91 endemic species of ray finned fishes, the actinopterygians. Some fish species have never been reported until recently, one such example is the long-horned Cowfish (Ostracidae), best known for inhabiting only certain areas of the Indo-Pacific oceans. However, this fish has been reported from Bay of Bengal, Costs of Goa and Lakshadweep. This is the first time this fish has been found on the Coasts of Karnataka, west coast of India. Given its morphology, the migration seems highly improbable since it can only endure lethargic swimming. In this article, we have explored the role of other external forces that could have contributed to its range extension in the east and west coast of India. In this paper we are exploring the role of Super Cyclones as a vital force in determining the expansion of its range. The fact of collection longhorn cowfish for the first time from the coast of Karwar coincides with the passage of Cyclone Ockhi on the west coast of India this suggests the possible role of super cyclone in its range extension. A year later, another specimen was recovered from a location very close to the first, which indicates the successful establishment of Lacturia cornuta in its new environment. Related morphometric and meristic studies of our specimens are consistent with previous studies reported in the Bay of Bengal. Getting more out of FLAG-Tag co-immunoprecipitation mass spectrometry experiments using FAIMS. Co-immunoprecipitation of proteins coupled to mass spectrometry is critical for the understanding of protein interaction networks. In instances where a suitable antibody is not available, it is common to graft synthetic tags onto target protein sequences and allowing the use of commercially available antibodies for affinity purification. A common approach is through FLAG-Tag co-immunoprecipitation. To allow the selective elution of protein complexes, competitive displacement using a large molar excess of the tag peptides is often carried out. Yet, this creates downstream challenges for the mass spectrometry analysis due to the presence of large quantities of these peptides. Here, we demonstrate that Field Asymmetric Ion Mobility Spectrometry (FAIMS), a gas phase ion separation device prior to mass spectrometry analysis can be applied to FLAG-Tag co-immunoprecipitation experiment to increase the depth of protein coverage. By excluding these abundant tag peptides, we were able to observe deeper coverage of interacting proteins and as a result, deeper biological insights, without the need for additional sample handling or altering sample preparation protocols. Reticular Dysgenesis-associated Adenylate Kinase 2 deficiency causes failure of myelopoiesis through disordered purine metabolism Reticular Dysgenesis is a particularly grave form of severe combined immunodeficiency that affects the adaptive and innate immune system. Patients suffer from congenital neutropenia, lymphopenia, and deafness. The disease is caused by biallelic loss of function in mitochondrial Adenylate Kinase 2 (AK2). AK2 mediates the phosphorylation of AMP to ADP, as substrate for ATP synthesis. Accordingly, declining oxidative phosphorylation has been postulated as the driver of disease pathology. The mechanistic basis, however, remains incompletely understood. Single cell RNA-sequencing of patient bone marrow cells implicated altered RNA catabolism and ribonucleoprotein synthesis in the pathogenesis of Reticular Dysgenesis. To investigate these findings, we developed a disease model based on CRISPR-mediated disruption of the AK2 gene in primary human hematopoietic stem cells. We found that AK2-deficient myeloid progenitor cells not only have compromised mitochondrial energy metabolism and increased AMP levels, but also NAD+ and aspartate depletion, metabolites that rely on TCA-cycle activity for regeneration and synthesis. Furthermore, AK2-deficient cells exhibited strikingly increased levels of the purine nucleotide precursor IMP, decreased cellular RNA content, ribosome subunit expression, protein synthesis and a profoundly hypo-proliferative phenotype. The rise in IMP levels stemmed from increased AMP deamination. Pharmacologic inhibition of AMP deaminase normalized IMP levels in AK2-deficient cells, but further aggravated the disease phenotype, pointing to AMP catabolism as a possible metabolic adaptation to mitigate AMP-mediated toxicity. Inducing an adenosine disequilibrium in control cells produced a similar myeloid maturation defect.

This study shows that AK2 deficiency globally curtailed mitochondrial metabolism resulting in NAD+ and aspartate deficiency and disordered purine metabolism. AMP accumulation and its detrimental effects on ribonucleotide synthesis capacity may contribute to the failure of myelopoiesis in Reticular Dysgenesis. 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. Designer Scaffolds for Interfacial Bioengineering In regenerative medicine, the healing of the interfacial zone between tissues is a major challenge, yet approaches for studying the complex microenvironment of this interface remain lacking. Here, we study these complex living interfaces by manufacturing modular "blocks" of naturally porous decellularized plant-derived scaffolds of varying shapes and sizes with a computer numerical controlled mill. We demonstrate how each scaffold can be seeded with different cell types and easily assembled in a manner akin to LEGO bricks to create an engineered tissue interface (ETI). Cells were observed to migrate across the interface formed between an empty scaffold and a scaffold pre-seeded with cells. However, when both scaffolds contained cells, only a shallow cross-over zone of cell infiltration formed at the interface. As a proof-of-concept study, we utilized ETIs to investigate the interaction between lab grown bone and connective tissues. Consistent with the above, a cross-over zone of the two distinct cell types formed at the interface between scaffolds, otherwise the populations remained distinct. Finally, we demonstrate how ETIs are biocompatible in vivo, becoming vascularized and integrated into surrounding tissue after implantation. This work creates new tissue design avenues for understanding biological processes or the development of synthetic artificial tissues. Optogenetic dissection of actomyosin-dependent mechanics underlying tissue fluidity Rapid epithelial tissue flows are essential to building and shaping developing embryos. However, it is not well understood how the mechanical properties of tissues and the forces driving them to flow are jointly regulated to accommodate rapid tissue remodeling. To dissect the roles of actomyosin in the mechanics of epithelial tissue flows, here we use two optogenetic tools, optoGEF and optoGAP, to manipulate Rho/Rho-kinase signaling and actomyosin contractility in the germband epithelium, which flows via convergent extension movements during Drosophila body axis elongation. The ability to perturb actomyosin in the tissue allows us to analyze the effects of actomyosin on cell rearrangements, tissue tensions, and tissue mechanical properties. We find that either optogenetic activation or deactivation of Rho1 signaling and actomyosin contractility at the apical surface of the germband disrupts cell rearrangements and tissue-level flows. Rho1 activation leads to poorly oriented rearrangements that are associated with a redistribution of myosin II from the junctional to the medial-apical domain, whereas Rho1 deactivation leads to fewer, slower cell rearrangements that are associated with decreased junctional and medial myosin. By probing mechanical tensions in the tissue using laser ablation and inferring tissue mechanical properties from cell packings, we find that actomyosin influences both the anisotropic forces that drive tissue flow and the mechanical properties of the tissue that resist flow, leading to complex relationships between actomyosin activity and tissue fluidity. Moreover, our results link the subcellular distribution of myosin II to tissue tension and cell packings, revealing how junctional and medial myosin have differential roles in promoting and orienting cell rearrangements during tissue flow. An IPI based immune prognostic model for diffuse large B-cell lymphoma BackgroundInternational Prognostic Index (IPI) was widely used to better discriminate prognosis of patients with diffuse large B-cell lymphoma (DLBCL). However, there is a significant need to identify novel valuable biomarkers in the context of targeted therapies, such as immune checkpoint blockade (ICB) therapy.

MethodsGene expression data and clinical information of DLBCL were obtained from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) datasets. 371 immune-related hub genes in DLBCL patients with different IPI levels were identified by weighted gene co-expression network analysis (WGCNA), and 8 genes were selected to construct an IPI-based immune prognostic model (IPI-IPM). Afterward, the genetic, somatic mutational and molecular profiles of IPI-IPM subgroups were analyzed, as well as the potential clinical response of ICB in different IPI-IPM subgroups.

ResultsThe IPI-IPM was constructed based on the expression of CMBL, TLCD3B, SYNDIG1, ESM1, EPHA3, HUNK, PTX3 and IL12A, where high-risk patients had shorter overall survival (OS) than low-risk patients, consistent with the results in the GEO cohorts. The comprehensive results showed that high IPI-IPM risk scores were correlated with immune-related signaling pathways, high KMT2D and CD79B mutation rates, as well as up-regulation of inhibitory immune checkpoints including PD-L1, BTLA and SIGLEC7, indicating more potential response to ICB therapy.

ConclusionThe IPI-IPM has independent prognostic significance for DLBCL patients, which provides an immunological perspective to elucidate the mechanisms on tumor progression, also sheds a light on developing immunotherapy for DLBCL. Actin isovariant ACT7 controls root meristem development in Arabidopsis through modulating auxin and ethylene responses Meristem is the most functionally dynamic part in a plant body. The shaping of the meristem requires constant cell division and cell elongation, which are influenced by hormones and cell cytoskeletal component, actin. Although the roles of hormones in modulating meristem development have been extensively studied, the role of actin in this process is still elusive. Using the single and double mutants of the vegetative class actin, we demonstrate that actin isovariant ACT7 plays an important role in root meristem development. In absence of ACT7, but not ACT8 and ACT2, depolymerization of actin was observed. Consistently, act7 mutant showed reduced cell division, cell elongation, and meristem length. Intracellular distribution and trafficking of auxin transport proteins in the actin mutants revealed that ACT7 specifically functions in root meristem to facilitate the trafficking of auxin efflux carriers PIN1 and PIN2, and consequently the transport of auxin. Compared with act7, act7act8 double mutant shows slightly enhanced phenotypic response and altered intracellular trafficking. The altered distribution of auxin in act7 and act7act8 affects the response of the roots to ethylene but not to cytokinin. Collectively, our results suggest that ACT7 dependent auxin-ethylene response plays a key role in controlling Arabidopsis root meristem development. Drug screening at single-organoid resolution via bioprinting and interferometry High-throughput drug screening is an established approach to investigate tumor biology and identify therapeutic leads. Traditional platforms for high-throughput screening use two-dimensional cultures of immortalized cell lines which do not accurately reflect the biology of human tumors. More clinically relevant model systems, such as three-dimensional tumor organoids, can be difficult to screen and scale. For example, manually seeded organoids coupled to destructive endpoint assays allow for the characterization of response to treatment, but do not capture the transitory changes and intra-sample heterogeneity underlying clinically observed resistance to therapy. We therefore developed a pipeline to generate bioprinted tumor organoids linked to label-free, real-time imaging via high-speed live cell interferometry (HSLCI) and machine learning-based quantitation of individual organoids. Bioprinting cells gives rise to 3D organoid structures that preserve tumor histology and gene expression. HSLCI imaging in tandem with machine learning-based image segmentation and organoid classification tools enables accurate, label-free parallel mass measurements for thousands of bioprinted organoids. We demonstrate that our method quantitatively identifies individual organoids as insensitive, transiently sensitive, or persistently sensitive to specific treatments. This opens new avenues for rapid, actionable therapeutic selection using automated tumor organoid screening. N-VEGF, the forgotten isoform of VEGF-A, activate cellular stress survival circuits Vascular endothelial growth factor A (VEGF-A) is a secreted protein that stimulates angiogenesis in response to hypoxia. Under hypoxic conditions, a non-canonical long isoform called L-VEGF is concomitantly expressed with VEGF-A. Once translated, L-VEGF it is proteolytically cleaved to generate N-VEGF and VEGF-A. Interestingly, while VEGF-A is secreted and affects the surrounding cells, N-VEGF is mobilized to the nucleus. This suggests that N-VEGF participates in transcriptional response to hypoxia. In this study, we performed a series of complementary experiments to examine the functional role of N-VEGF. Strikingly, we found that the mere expression of N-VEGF followed by its hypoxia-independent mobilization to the nucleus was sufficient to induce key genes associated with angiogenesis, such as Hif1, VEGF-A isoforms, as well as genes associated with cell survival under hypoxia. Complementarily, when N-VEGF was genetically depleted, key hypoxia-induced genes were downregulated and cells were significantly susceptible to hypoxia-mediated apoptosis. This is the first reports of N-VEGF serving as an autoregulatory arm of VEGF-A. Further experiments will be needed to determine the role of N-VEGF in cancer and embryogenesis. Rho Kinase regulates neutrophil NET formation that is involved in UVB-induced skin inflammation ObjectiveUltraviolet B (UVB) is an important trigger of skin inflammation and lupus with leukocyte recruitment to inflamed skin. We recently reported the involvement of neutrophil NETosis in UVB-induced skin inflammation, and that NETotic nuclear envelope rupture is driven by PKC-mediated nuclear lamin B disassembly. To address the role of Actin cytoskeleton in NETosis, we investigated the effects of Rho kinase (ROCK) and its downstream actomyosin cytoskeletal networks on PKC nuclear translocation and NET formation, as well as their involvement in UVB-induced skin inflammation.

MethodsWe studied the dynamic changes of ROCK and actomyosin cytoskeletal networks during NETosis induction and their involvement in PKC nuclear translocation. Using mice with hematopoietic-specific ROCK1 deficiency, we investigated the effects of ROCK1 deficiency on NETosis, and its involvement in UVB-induced skin inflammation.

ResultsOur time course studies demonstrated the dynamic changes of actin polymerization and ROCK activation, support the role of actin cytoskeleton in nuclear translocation of cytosolic PKC in early stage of NETosis induction. Inhibition of actin polymerization or key molecules of the ROCK/MLCK/myosin pathway decreased PKC nuclear translocation and NET formation. Genetic deficiency of ROCK1, inhibited NETosis ex vivo and in vivo, decreased extracellular display of NET-associated IL-17A, TNF, IFN{gamma}, and IFN in inflamed skin, which were correlated with the ameliorated skin inflammation in UVB-irradiated mice with hematopoietic-specific ROCK1 deficiency.

ConclusionsROCK regulated NETosis through modulation of PKC nuclear translocation via actomyosin cytoskeletal networks in neutrophils. ROCK1 deficiency ameliorated UVB- induced skin inflammation by attenuation of NETosis and NET-associated cytokines.