Human cells contain myriad excised linear introns with potential functions in gene regulation and as RNA biomarkers By using TGIRT-seq, we identified >8,500 short full-length excised linear intron (FLEXI) RNAs in human cells. Subsets of FLEXIs accumulated in a cell-type specific manner, and [~]200 corresponded to agotrons or mirtrons or encoded snoRNAs. Analysis of CLIP-seq datasets identified potential interactions between FLEXIs and >100 different RNA-binding proteins (RBPs), 53 of which had binding sites in [&ge;]30 different FLEXIs. In addition to proteins that function in RNA splicing, these 53 RBPs included transcription factors, chromatin remodeling proteins, and cellular growth regulators that impacted FLEXI host gene alternative splicing and/or mRNA levels in knockdown datasets. We computationally identified six groups of RBPs whose binding sites were enriched in different subsets of FLEXIs: AGO1-4 and DICER associated with agotrons and mirtrons; AATF, DKC1, NOLCI, and SMNDC1 associated with snoRNA-encoding FLEXIs; two different combinations of alternative splicing factors found in stress granules; and two novel RBP-intron combinations, one including LARP4 and PABC4, which function together in the cytoplasm to regulate ribosomal protein translation. Our results suggest a model in which proteins involved in transcriptional regulation, alternative splicing, or post-splicing secondary functions bind and stabilize cell-type specific subsets of FLEXIs that perform different biological functions and have potential utility as biomarkers.
Full-length Genome of a Ogataea polymorpha strain CBS4732 ura3Δ reveals large duplicated segments in subtelomeric regions BackgroundCurrently, methylotrophic yeasts (e.g., Pichia pastoris, Ogataea polymorpha, and Candida boindii) are subjects of intense genomics studies in basic research and industrial applications. In the genus Ogataea, most research is focused on three basic O. polymorpha strains--CBS4732, NCYC495, and DL-1. However, the relationship between CBS4732, NCYC495, and DL-1 remains unclear, as the genomic differences between them have not be exactly determined without their high-quality complete genomes. As a nutritionally deficient mutant derived from CBS4732, the O. polymorpha strain CBS4732 ura3{Delta} (named HU-11) is being used for high-yield production of several important proteins or peptides. HU-11 has the same reference genome as CBS4732 (noted as HU-11/CBS4732), because the only genomic difference between them is a 5-bp insertion.

ResultsIn the present study, we have assembled the full-length genome of O. polymorpha HU-11/CBS4732 using high-depth PacBio and Illumina data. Long terminal repeat (LTR) retrotransposons, rDNA, 5 and 3 telomeric, subtelomeric, low complexity and other repeat regions were curated to improve the genome quality. Particularly, we detected large duplicated segments (LDSs) in the subtelomeric regions and exactly determined all the structural variations (SVs) between CBS4732 and NCYC495.

New findings mainly include(1) the genomic differences between HU-11/CBS4732 and NCYC495 include single nucleotide polymorphisms, small insertions and deletions, and only three SVs; (2) six genes were incorporated into CBS4732 from Cyberlindnera jadinii by horizontal gene transfer and may bring HU-11/CBS4732 new biological functions or physiological properties; (3) many recombination events may have occurred on chromosome 4 and 5 of CBS4732 and NCYC495 ancestors and two large segments were acquired by CBS4732 and NCYC495 from chromosome 6 and C. jadinii during recombination, respectively; and (4) the genome expansion in methylotrophic yeasts is mainly driven by large segment duplication in subtelomeric regions.

ConclusionsThe present study preliminarily revealed the complex relationship between CBS4732, NCYC495, and DL-1. The new findings provide new opportunities for in-depth understanding of genome evolution in methylotrophic yeasts and lay the foundations for the industrial applications of O. polymorpha CBS4732, NCYC495, DL-1, and their derivative strains. The full-length genome of the O. polymorpha strain HU-11/CBS4732 should be included into the NCBI RefSeq database for future studies of Ogataea spp..
Mechanisms Underlying WNT-mediated Priming of Human Embryonic Stem Cells Embryogenesis is guided by a limited set of signaling pathways dynamically expressed in different places. How a context dependent signaling response is generated has been a central question of developmental biology, which can now be addressed with in vitro models of human embryos that are derived from embryonic stem cells (hESCs) called gastruloids. Our previous work demonstrated that during early self-organization of gastruloids, cells chronicle signaling hierarchy. Only cells that have been exposed (primed) by WNT signaling can respond to subsequent Activin exposure and differentiate to mesendodermal (ME) fates. Here, we show that WNT priming does not alter SMAD2 binding nor its chromatin opening, but rather, acts by inducing the expression of the SMAD2 co-factor, EOMES. Expression of EOMES is sufficient to replace WNT upstream of Activin-mediated ME differentiation, thus unveiling the mechanistic basis for priming and cellular memory in early development. Temporal Scaling of Human Scalp-Recorded Potentials Much of human behaviour is governed by common processes that unfold over varying timescales. Standard event-related potential analysis assumes fixed-latency responses relative to experimental events. However, recent single unit recordings in animals have revealed neural activity scales to span different durations during behaviours demanding flexible timing. Here, we employed a general linear modelling approach using a novel combination of fixed-duration and variable-duration regressors to unmix fixed-time and scaled-time components in human magneto/electroencephalography (M/EEG) data. We use this to reveal consistent temporal scaling of human scalp-recorded potentials across four independent EEG datasets, including interval perception, production, prediction and value-based decision making. Between-trial variation in the temporally scaled response predicts between-trial variation in subject reaction times, demonstrating the relevance of this temporally scaled signal for temporal variation in behaviour. Our results provide a general approach for studying flexibly timed behaviour in the human brain. The compound of baicalein, wogonin, and oroxylin-A inhibiting EMT in A549 cell line via PI3K/AKT-TWIST1-glycolysis pathway Non-small cell lung cancer (NSCLC) is a worldwide disease with high morbidity and mortality, which is most derived from its metastasis. Some researches show that the epithelial-mesenchymal transition (EMT) process promotes lung cancer cells migration and invasion, leading to NSCLC metastasis. Total Flavonoid Aglycones Extract (TFAE) isolated from Scutellaria baicalensis was reported to inhibit tumor growth and induce apoptosis. In this study, we found that baicalein, wogonin, and oroxylin-A were the active compounds of TFAE. After reconstructing with these three compounds (baicalein (65.8%), wogonin (21.2%), and oroxylin-A (13.0%)), the reconstructed TFAE (reTFAE) present inhibitory effect on the EMT process of A549 cells. Then, bioinformatic technology was employed to elucidate the potential pharmacodynamic mechanism network of reTFAE. We identified the relationship between reTFAE and PI3K/Akt signaling pathways, with TWIST1 as the key protein. LY294002, the inhibitor of PI3K/Akt signaling pathway, and knock-down TWIST1 could significantly enhance the efficacy of reTFAE, with increasing expression of epithelial markers and decreasing expression of mesenchymal markers in A549 cells at the same time. Furthermore, stable isotope dimethyl-labeled proteomics technology was conducted to complement the follow-up mechanism that the EMT-inhibition process may be realized through glycolysis pathway. In conclusion, we claim that TWIST1-targeted flavonoids could provide a new strategy to inhibit EMT progress for the treatment of NSCLC. P2X2 receptor subunit interfaces are missense variant hotspots where mutations tend to increase apparent ATP affinity Background and PurposeP2X receptors (P2XRs) are trimeric ligand-gated ion channels (LGICs) that open a cation-selective pore in response to ATP binding to their large extracellular domain (ECD). The seven known P2XR subtypes can assemble as homo- or heterotrimeric complexes and contribute to numerous physiological functions, including nociception, inflammation and hearing. The overall structure of P2XRs is well established, but little is known about the spectrum and prevalence of human genetic variations and the functional implications in specific domains.

Experimental ApproachHere we examine the impact of P2X2 receptor (P2X2R) inter-subunit interface missense variants identified in the human population or through structural predictions. We test both single and double mutants through electrophysiological and biochemical approaches.

Key resultsWe demonstrate that predicted ECD inter-subunit interfaces display a higher-than-expected density of missense variations and that the majority of mutations that disrupt putative inter-subunit interactions result in channels with higher apparent ATP affinity. Lastly, we show that double mutants at the subunit interface show significant energetic coupling, especially if located in close proximity.

Conclusions and ImplicationsWe provide the first structural mapping of the mutational burden across the human population in a LGIC and show that the density of missense mutations is constrained between protein domains, indicating evolutionary selection at the domain level. Our data may indicate that, unlike other LGICs, P2X2Rs have evolved an intrinsically high threshold for activation, possibly to allow for additional modulation or as a cellular protection mechanism against overstimulation.

Bullet point summary What is already known: O_LIP2X2 receptors are ATP-activated ion channels implicated in hearing and nociceptice pathways
C_LI
 What this study adds: O_LIA structural mapping of missense variants observed in the human population
C_LIO_LIWe identify the intersubunit-interface as a variant hotspot and decipher functional impact of mutations
C_LI
 Clinical significance: O_LIThe development of both inhibitors and activators of P2X2 receptor function may be required
C_LI Ribosome actuator via the SpoT/ppGpp pathway to mitigate gene overexpression burden Gene overexpression causes a burden on ribosomes, which decreases the expression of other genes and growth rate. Here, we engineer a genetic ribosome actuator that increases gene expression rate and growth rate by over 150% and 80%, respectively. The actuator expresses the hydrolysis domain of the SpoT enzyme (SpoTH) in a bacterial strain with high basal levels of ppGpp, thus increasing ribosome availability. We use this actuator in a feedforward control arrangement, wherein SpoTH is co-expressed with a gene of interest (GOI). Without SpoTH co-expression, activation of the GOI decreases the expression of a constitutive gene and growth rate by over 80% and 40%, respectively. By contrast, with SpoTH co-expression, the ribosome binding site (RBS) of SpoTH can be tuned to achieve inappreciable change of either the expression of the constitutive gene or growth rate when the GOI is activated to the same level. We finally employ the actuator to keep growth rate constant while expressing dCas9 to a level that would otherwise decrease growth rate by more than 40% due to toxicity. This actuator will be relevant in applications where high cellular protein concentrations are desired without growth defects or global reduction of gene expression. Dynamic modulation of enhancer responsiveness by core promoter elements in living Drosophila embryo Regulatory interactions between enhancers and core promoters are fundamental for the temporal and spatial specificity of gene expression in development. The central role of core promoters is to initiate productive transcription in response to enhancers activation cues. However, it has not been systematically assessed how individual core promoter elements affect the induction of transcriptional bursting by enhancers. Here, we provide evidence that each core promoter element dynamically alters "enhancer responsiveness" by changing functional parameters of transcriptional bursting in developing Drosophila embryos. Quantitative live imaging analysis revealed that the timing and the continuity of burst induction are common regulatory steps on which core promoter elements impact. We also show that core promoters use the upstream TATA box as an amplifier of transcriptional bursting upon interaction with enhancers. Genome editing analysis of the pair-rule gene fushi tarazu further revealed that the endogenous TATA box and the downstream core promoter element are both essential for the correct expression and the function of key developmental genes in early embryos. We suggest that core promoter elements serve as a key regulatory module in converting enhancer activity into transcription dynamics during animal development. Heterogeneous recruitment abilities to RNA polymerases generate nonlinear scaling of gene expression level with cell volume While most genes expression levels are proportional to cell volumes, some genes exhibit nonlinear scaling between their expression levels and cell volume. Therefore, their mRNA and protein concentrations change as the cell volume increases, which often have crucial biological functions such as cell-cycle regulation. However, the biophysical mechanism underlying the nonlinear scaling between gene expression and cell volume is still unclear. In this work, we show that the nonlinear scaling is a direct consequence of the heterogeneous recruitment abilities of promoters to RNA polymerases based on a gene expression model at the whole-cell level. Those genes with weaker (stronger) recruitment abilities than the average ability spontaneously exhibit superlinear (sublinear) scaling with cell volume. Analysis of the promoter sequences and the nonlinear scaling of Saccharomyces cerevisiaes mRNA levels shows that motifs associated with transcription regulation are indeed enriched in genes exhibiting nonlinear scaling, in concert with our model. 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. Temperature Increase Drives Critical Slowing Down of Fish Ecosystems Fish ecosystems perform ecological functions that are critically important for the sustainability of marine ecosystems, such as global food security and carbon stock. During the 21st century, significant global warming caused by climate change has created pressing challenges for fish ecosystems that threaten species existence and global ecosystem health. Here, we study a coastal fish community in Maizuru Bay, Japan, and investigate the relationships between fluctuations of ST, abundance-based species interactions and salient fish biodiversity. Observations show that a local 20% increase in temperature from 2002 to 2014 underpins a long-term reduction in fish diversity ([~]25%) played out by some native and invasive species (e.g. Chinese wrasse) becoming exceedingly abundant; this causes a large decay in commercially valuable species (e.g. Japanese anchovy) coupled to an increase in ecological productivity. The fish community is analyzed considering five temperature ranges to understand its atemporal seasonal sensitivity to ST changes, and long-term trends. An optimal information flow model is used to reconstruct species interaction networks that emerge as topologically different for distinct temperature ranges and species dynamics. Networks for low temperatures are more scale-free compared to ones for intermediate (15-20{degrees}C) temperatures in which the fish ecosystem experiences a first-order phase transition in interactions from locally stable to metastable and globally unstable for high temperatures states as suggested by abundance-spectrum transitions. The dynamic dominant eigenvalue of species interactions shows increasing instability for competitive species (spiking in summer due to intermediate-season critical transitions) leading to enhanced community variability and critical slowing down despite higher time-point resilience. Native competitive species whose abundance is distributed more exponentially have the highest total directed interactions and are keystone species (e.g. Wrasse and Horse mackerel) for the most salient links with cooperative decaying species. Competitive species, with higher eco-climatic memory and synchronization, are the most affected by temperature and play an important role in maintaining fish ecosystem stability via multitrophic cascades (via cooperative-competitive species imbalance), and as bioindicators of change. More climate-fitted species follow temperature increase causing larger divergence divergence between competitive and cooperative species. Decreasing dominant eigenvalues and lower relative network optimality for warmer oceans indicate fishery more attracted toward persistent oscillatory states, yet unpredictable, with lower cooperation, diversity and fish stock despite the increase in community abundance due to non-commercial and venomous species. We emphasize how changes in species interaction organization, primarily affected by temperature fluctuations, are the backbone of biodiversity dynamics and yet for functional diversity in contrast to taxonomic richness. Abundance and richness manifest gradual shifts while interactions show sudden shift. The work provides data-driven tools for analyzing and monitoring fish ecosystems under the pressure of global warming or other stressors. Abundance and interaction patterns derived by network-based analyses proved useful to assess ecosystem susceptibility and effective change, and formulate predictive dynamic information for science-based fishery policy aimed to maintain marine ecosystems stable and sustainable. Naming Human Diseases: Ethical Principles of Curating Exclusive Substitute for Inopportune Nosology BackgroundIn the medical sphere, understanding naming conventions strengthen the integrity of naming human diseases remains nominal rather than substantial yet. Since the current nosology-based standard for human diseases could not offer a one-size-fits-all corrective mechanism, many idiomatic but flawed names frequently appear in scientific literature and news outlets at the cost of sociocultural impacts.

ObjectiveWe attempt to examine the ethical oversights of current naming practices and propose heuristic rationales and approaches to determine a pithy name instead of an inopportune nosology.

MethodsFirst, we examined the compiled global online news volumes and emotional tones on some inopportune nosology like German measles, Middle Eastern Respiratory Syndrome, Spanish flu, Hong Kong flu, and Huntingtons disease in the wake of COVID-19. Second, we prototypically scrutinize the lexical dynamics and pathological differentials of German measles and common synonyms by leveraging the capacity of the Google Books Ngram Corpus. Third, we demonstrated the empirical approaches to curate an exclusive substitute for an anachronistic nosology German measles based on deep learning models and post-hoc explanations.

ResultsThe infodemiological study shows that the public informed the offensive names with extremely negative tones in textual and visual narratives. The findings of the historiographical study indicate that many synonyms of German measles did not survive, while German measles became an anachronistic usage, and rubella has taken the dominant place since 1994. The PubMedBERT model could identify rubella as a potential substitution for German measles with the highest semantic similarity. The results of the semantic drift experiments further indicate that rubella tends to survive during the ebb and flow of semantic drift.

ConclusionsOur findings indicate that the nosological evolution of anachronistic names could result in sociocultural impacts without a corrective mechanism. To mitigate such impacts, we introduce some ethical principles for formulating an improved naming scheme. Based on deep learning models and post-hoc explanations, our illustrated experiments could provide hallmark references to the remedial mechanism of naming practices and pertinent credit allocations. An epigenomic shift in amygdala marks the transition to maternal behaviors in alloparenting virgin female mice In many species, adults care for young offspring that are not their own, a phenomenon called alloparenting. However, most nonparental adults must be sensitized by repeated or extended exposures to newborns before they will robustly display parental-like behaviors. To capture neurogenomic events underlying the transition to active parental caring behaviors, we analyzed brain gene expression and chromatin profiles of virgin female mice co-housed with mothers during pregnancy and after birth. After an initial display of antagonistic behaviors and a surge of defense-related gene expression, we observed a dramatic shift in the chromatin landscape specifically in amygdala of the pup-exposed virgin females, accompanied by a dampening of anxiety-related gene expression. This epigenetic shift coincided with hypothalamic expression of the oxytocin gene and the emergence of behaviors and gene expression patterns classically associated with maternal care. The results outline a neurogenomic program associated with dramatic behavioral changes and suggest molecular networks relevant to human postpartum mental health. 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. Regulation of the PKD2 channel function by TACAN Autosomal dominant polycystic kidney disease (ADPKD) is caused by mutations in membrane receptor PKD1 or cation channel PKD2. TACAN (also named TMEM120A), recently reported as an ion channel in neuron cells for mechano and pain sensing, is also distributed in diverse non-neuronal tissues such as kidney, heart and intestine, suggesting its involvement in other functions. In this study, we found that TACAN is in complex with PKD2 in native renal cell lines. Using the two-electrode voltage clamp in Xenopus oocytes we found that TACAN inhibited the channel activity of PKD2 gain-of-function mutant F604P. The first and last transmembrane domains of TACAN were found to interact with the PKD2 C-and N-terminal portions, respectively. We showed that the TACAN N-terminus acted as a blocking peptide and that TACAN inhibits the PKD2 function through the PKD2/TACAN binding. By patch clamping in mammalian cells, we found that TACAN inhibits both the single channel conductance and open probability of PKD2 and mutant F604P. PKD2 co-expressed with TACAN, but not PKD2 alone, exhibited pressure sensitivity. Furthermore, we also found that TACAN aggravates PKD2-dependent tail curvature and pronephric cysts in larval zebrafish, in support of the in vitro inhibitory effects of TACAN. In summary, this study revealed that TACAN acts as a PKD2 inhibitor and mediates mechano sensitivity of the PKD2/TACAN channel complex. 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. The isoforms of pyruvate kinase act as nutrient sensors for the β-cell KATP channel Insulin release from pancreatic {beta}-cells is triggered by the metabolism-dependent closure of ATP-sensitive K+ (KATP) channels. Here, we used isoform-specific deletion of pyruvate kinase (PK) to identify the unique roles of PKm1 and PKm2 in the {beta}-cell nutrient response. In the absence of PK activity, cytosolic ATP/ADP was insufficient to close KATP channels. Both PK isoforms were found on the plasma membrane, where they endow KATP channels with nutrient sensitivity and specificity. Although PKm1 accounts for >90% of cellular activity, PKm2 controls glucosedependent KATP closure. Amino acids, which generate mitochondrial phosphoenolpyruvate (PEP) without allosterically activating PKm2, depend on PKm1 to raise ATP/ADP and close KATP channels. Shifting the {beta}-cell from PKm2 to PKm1 correspondingly increased amino acid sensitivity and boosted secretory output. In addition to these membrane depolarizing "on-switches," {beta}-cell deletion of PCK2 exposed a mitochondrial PEP-dependent "off-switch" that facilitates Ca2+ extrusion. Restricting PCK2 prolonged Ca2+ influx and increased insulin secretion. In summary, it is not the KATP channel alone, but the PK-KATP microcompartment that confers nutrient sensitivity to the {beta}-cell plasma membrane. The differential response of the PK isoforms to the glycolytic and mitochondrial sources of PEP influences {beta}-cell nutrient preference, and controls the oscillatory cycle regulating insulin secretion. Substrate specialization in microbes is driven by biochemical constraints of dynamic flux sensing Central carbon metabolism is highly conserved across microbial species, but can catalyze very different pathways depending on the organism and their ecological niche. Here, we study the dynamic re-organization of central metabolism after switches between the two major opposing pathway configurations of central carbon metabolism, glycolysis and gluconeogenesis in Escherichia coli, Pseudomonas aeruginosa and Pseudomonas putida. We combined growth dynamics and dynamic changes of intracellular metabolite levels with a coarse-grained model that integrates fluxes, regulation, protein synthesis and growth and uncovered fundamental limitations of the regulatory network: after nutrient shifts, metabolite concentrations collapse to their equilibrium, rendering the cell unable to sense which direction the flux is supposed to flow through the metabolic network. The cell can partially alleviate this by picking a preferred direction of regulation at the expense of increasing lag times in the opposite direction. Moreover, decreasing both lag times simultaneously comes at the cost of reduced growth rate or higher futile cycling between metabolic enzymes. These three trade-offs can explain why microorganisms specialize for either glycolytic or gluconeogenic substrates and can help elucidate the complex growth patterns exhibited by different microbial species.

Graphical synopsis O_FIG_DISPLAY_L [Figure 1] M_FIG_DISPLAY C_FIG_DISPLAY

Standfirst textMicrobes face a series of fundamental trade-offs that limit their ability to optimize simultaneously for both glycolytic and gluconeogenic growth.

Bullet pointsO_LILag times between glycolysis and gluconeogenesis show asymmetry in many microbes: A long lag in one direction, but a short lag in the other.
C_LIO_LILong lag times are caused by an inability to sense fluxes after nutrient shifts.
C_LIO_LIWith existing regulation, lag time asymmetry can only be overcome by reducing either growth rate or increasing futile cycling in metabolism.
C_LI Water levels primarily drive variation in photosynthesis and nutrient use of scrub Red Mangroves in the southeastern Florida Everglades We investigated how mangrove-island micro-elevation (i.e., habitat: center vs. edge) affects tree physiology in a scrub mangrove forest of the southeastern Everglades. We measured leaf gas exchange rates of scrub Rhizophora mangle L. trees monthly during 2019, hypothesizing that CO2 assimilation (Anet) and stomatal conductance (gsw) would decline with increasing water levels and salinity, expecting more-considerable differences at mangrove-island edges than centers, where physiological stress is greatest. Water levels varied between 0 and 60 cm from the soil surface, rising during the wet season (May-October) relative to the dry season (November-April). Porewater salinity ranged from 15 to 30 ppt, being higher at mangrove-island edges than centers. Anet maximized at 15.1 {micro}mol m-2 s-1, and gsw was typically <0.2 mol m-2 s-1, both of which were greater in the dry than the wet season and greater at island centers than edges, with seasonal variability being roughly equal to variation between habitats. After accounting for season and habitat, water level positively affected Anet in both seasons but did not affect gsw. Our findings suggest that inundation stress (i.e., water level) is the primary driver of variation in leaf gas exchange rates of scrub mangroves in the Florida Everglades, while also constraining Anet more than gsw. The interaction between inundation stress due to permanent flooding and habitat varies with season as physiological stress is alleviated at higher-elevation mangrove-island center habitats during the dry season. Freshwater inflows during the wet season, increase water levels and inundation stress at higher-elevation mangrove-island centers, but also potentially alleviate salt and sulfide stress in soils. Thus, habitat heterogeneity leads to differences in nutrient and water acquisition and use between trees growing in island centers versus edges, creating distinct physiological controls on photosynthesis, which likely affect carbon flux dynamics of scrub mangroves in the Everglades. Fast and Accurate Distance-based Phylogenetic Placement using Divide and Conquer Phylogenetic placement of query samples on an existing phylogeny is increasingly used in molecular ecology, including sample identification and microbiome environmental sampling. As the size of available reference trees used in these analyses continues to grow, there is a growing need for methods that place sequences on ultra-large trees with high accuracy. Distance-based placement methods have recently emerged as a path to provide such scalability while allowing flexibility to analyze both assembled and unassembled environmental samples. In this paper, we introduce a distance-based phylogenetic placement method, APPLES-2, that is more accurate and scalable than existing distance-based methods and even some of the leading maximum likelihood methods. This scalability is owed to a divide-and-conquer technique that limits distance calculation and phylogenetic placement to parts of the tree most relevant to each query. The increased scalability and accuracy enables us to study the effectiveness of APPLES-2 for placing microbial genomes on a data set of 10,575 microbial species using subsets of 381 marker genes. APPLES-2 has very high accuracy in this setting, placing 97% of query genomes within three branches of the optimal position in the species tree using 50 marker genes. Our proof of concept results show that APPLES-2 can quickly place metagenomic scaffolds on ultra-large backbone trees with high accuracy as long as a scaffold includes tens of marker genes. These results pave the path for a more scalable and widespread use of distance-based placement in various areas of molecular ecology. Liver environment-imposed constraints diversify movement strategies of liver-localized CD8 T cells Pathogen-specific CD8 T cells face the problem of finding rare cells that present their cognate antigen either in the lymph node or infected tissue. While quantitative details of T cell movement strategies in some tissues such lymph nodes or skin have been relatively well characterized, we still lack quantitative understanding of T cell movement in many other important tissues such as the spleen, lung, liver, and gut. Furthermore, how environmental constraints influence movement patterns of T cells in tissues remains incompletely characterized. Liver is one of the major organs in mice and humans and is the target site for many pathogens such as hepatitis B and C viruses in humans and Plasmodium parasites in multiple mammalian species. We developed a protocol to generate stable numbers of liver-located CD8 T cells, used intravital microscopy to record movement patterns of CD8 T cells in livers of live mice, and analyzed these and previously published data using well-established statistical and computational methods. We show that in most of our experiments Plasmodium-specific liver-localized CD8 T cells perform correlated random walks characterized by transiently superdiffusive displacement with persistence times of 10-15 min that exceed those observed for T cells in lymph nodes. Liver-localized CD8 T cells typically crawl on the lumenal side of liver sinusoids (i.e., are in the blood); simulating T cell movement in digital structures derived from the liver sinusoids illustrates that liver structure alone is sufficient to explain the relatively long superdiffusive displacement of T cells. In experiments when CD8 T cells in the liver poorly attach to the sinusoids (e.g., one week post immunization with radiation-attenuated Plasmodium sporozoites) T cells also undergo Levy flights - large displacements occurring due to cell deattaching from the endothelium, floating with the blood flow, and re-attaching at another location. Our analysis, thus, provides quantitative details of movement patterns of liver-localized CD8 T cells and illustrates how structural and physiological details of the tissue may impact T cell movement patterns. Impairment of SARS-CoV-2 spike-glycoprotein maturation and fusion-activity by nitazoxanide: an effect independent of spike variants emergence SARS-CoV-2, the causative agent of COVID-19, has caused an unprecedented global health crisis. The SARS-CoV-2 spike, a surface-anchored trimeric class-I fusion-glycoprotein essential for viral entry, represents a key target for developing vaccines and therapeutics capable of blocking virus invasion. The emergence of SARS-CoV-2 spike-variants that facilitate virus spread and may affect vaccine efficacy highlights the need to identify novel antiviral strategies for COVID-19 therapy. Here we demonstrate that nitazoxanide, an antiprotozoal agent with recognized broad-spectrum antiviral activity, interferes with SARS-CoV-2 spike biogenesis, hampering its maturation at an endoglycosidase H-sensitive stage. Engineering multiple SARS-CoV-2 variant-pseudoviruses and utilizing quantitative cell-cell fusion assays, we show that nitazoxanide-induced spike modifications hinder progeny virion infectivity as well as spike-driven pulmonary cell-cell fusion, a critical feature of COVID-19 pathology. Nitazoxanide, being equally effective against the ancestral SARS-CoV-2 Wuhan-spike and different emerging variants, including the Delta variant of concern, may represent a useful tool in the fight against COVID-19 infections. Material-mediated histogenesis using mechano-chemically microstructured cell niches Stem-cell derived tissue models are commonly cultured under globally-delivered stimuli that trigger histogenesis via self-organizing activity. However, the culture of such tissue models is prone to stochastic behavior, limiting the reproducibility of cellular composition and resulting in non-physiological architectures. To overcome these shortcomings, we developed a method for printing cell niche microenvironments with microstructured cues that mediate local histogenic processes, including mechanosensing and differentiation of selected cell types. Microstructured cues include independently tunable mechano-chemical properties, with conjugated peptides, proteins, and morphogens across a range of Youngs moduli. By rationally designing niches, we mediate the structure of tissues derived from stem-cell-progenitor sources, including a bone-fat assembly from stromal mesenchyme, and embryonic tissues derived from hiPSC. We show that microstructured cues can recapitulate mechano-chemical signals resembling early embryonic histogenesis. This outcome includes a role for niche mechanics in human embryonic organization, where soft niche mechanics bias markers of mesendodermal differentiation and epithelial-to-mesenchymal-transition (EMT), as well as a demonstration of a material-mediated morphogen signaling centers able to induce foci of mesenchymal and EMT differentiation. Thus, microstructured materials can mediate local histogenic processes to enhance the structure and composition of tissue models. 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. Bioinformatic and mechanistic analysis of the palmerolide PKS-NRPS biosynthetic pathway from the microbiome of an Antarctic ascidian Complex interactions exist between microbiomes and their hosts. Increasingly, defensive metabolites that have been attributed to host biosynthetic capability are now being recognized as products of host-associated microbes. These unique metabolites often have bioactivity targets in human disease and can be purposed as pharmaceuticals. Polyketides are a complex family of natural products that often serve as defensive metabolites for competitive or pro-survival purposes for the producing organism, while demonstrating bioactivity in human diseases as cholesterol lowering agents, anti-infectives, and anti-tumor agents. Marine invertebrates and microbes are a rich source of polyketides. Palmerolide A, a polyketide isolated from the Antarctic ascidian Synoicum adareanum, is a vacuolar-ATPase inhibitor with potent bioactivity against melanoma cell lines. The biosynthetic gene clusters (BGCs) responsible for production of secondary metabolites are encoded in the genomes of the producers as discrete genomic elements. A candidate palmerolide BGC was identified from a S. adareanum microbiome-metagenome based on a high degree of congruence with a chemical structure-based retrobiosynthetic prediction. Protein family homology analysis, conserved domain searches, active site and motif identification were used to identify and propose the function of the [~]75 kb trans-acyltransferase (AT) polyketide synthase-non-ribosomal synthase (PKS-NRPS) domains responsible for the stepwise synthesis of palmerolide A. Though PKS systems often act in a predictable co-linear sequence, this BGC includes multiple trans-acting enzymatic domains, a non-canonical condensation termination domain, a bacterial luciferase-like monooxygenase (LLM), and is found in multiple copies within the metagenome-assembled genome (MAG). Detailed inspection of the five highly similar pal BGC copies suggests the potential for biosynthesis of other members of the palmerolide chemical family. This is the first delineation of a biosynthetic gene cluster from an Antarctic microbial species, recently proposed as Candidatus Synoicohabitans palmerolidicus. These findings have relevance for fundamental knowledge of PKS combinatorial biosynthesis and could enhance drug development efforts of palmerolide A through heterologous gene expression. From fitting the average to fitting the individual: A cautionary tale for mathematicalmodelers An outstanding challenge in the clinical care of cancer is moving from a one-size-fits-all approach that relies on population-level statistics towards personalized therapeutic design. Mathematical modeling is a powerful tool in treatment personalization, as it allows for the incorporation of patient-specific data so that treatment can be tailor-designed to the individual. Herein, we work with a mathematical model of murine cancer immunotherapy that has been previously-validated against the average of an experimental dataset. We ask the question: what happens if we try to use this same model to perform personalized fits, and therefore make individualized treatment recommendations? Typically, this would be done by choosing a single fitting methodology, and a single cost function, identifying the individualized best-fit parameters, and extrapolating from there to make personalized treatment recommendations. Our analyses show the potentially problematic nature of this approach, as predicted personalized treatment response proved to be sensitive to the fitting methodology utilized. We also demonstrate how a small amount of the right additional experimental measurements could go a long way to improve consistency in personalized fits. Finally, we show how quantifying the robustness of the average response could also help improve confidence in personalized treatment recommendations.

Author summaryAs we enter the era of healthcare where personalized medicine becomes a more common approach to treating cancer patients, harnessing the power of mathematical models will only become more essential. Using a preclinical dataset on cancer immunotherapy, we explore the challenges and limitations that arise when trying to move from fitting and making predictions for the population-level average, to fitting and making predictions for an individual. We find that the standard of approach of picking a single fitting methodology and a single cost function may end up having limited predictive value when applied to individual data. We also show how having a small amount of the right additional experimental data, and establishing the robustness of average treatment response, can help improve confidence in personalized model predictions.