Nanoscale resolution imaging of the whole mouse embryos and larval zebrafish using expansion microscopy Nanoscale resolution imaging of whole vertebrates is required for a systematic understanding of human diseases, but this has yet to be realized. Expansion microscopy (ExM) is an attractive option for achieving this goal, but the expansion of whole vertebrates has not been demonstrated due to the difficulty of expanding hard body components. Here, we demonstrate whole-body ExM, which enables nanoscale resolution imaging of anatomical structures, proteins, and endogenous fluorescent proteins (FPs) of whole zebrafish larvae and mouse embryos by expanding them fourfold. We first show that post-digestion decalcification and digestion kinetics matching are critical steps in the expansion of whole vertebrates. Then, whole-body ExM is combined with the improved pan-protein labeling approach to demonstrate the three-dimensional super-resolution imaging of antibody- or FP-labeled structures and all major anatomical structures surrounding them. We also show that whole-body ExM enables visualization of the nanoscale details of neuronal structures across the entire body.
Neural Flip-Flops II: The Role of Cascaded Oscillators in Short-Term Memory, EEGs, and Epilepsy By enabling many brain structures state changes, the explicit cascaded oscillators proposed here can generate the rhythmic neural activity found in EEGs. The function of such synchronization in information processing systems is timing error avoidance. The narrow requirement for the oscillator input pulse duration suggests a possible relationship to the abnormal electrical activity characteristic of epileptic seizures. Together, flip-flops and synchronization by oscillators suggest a resolution to the longstanding controversy of whether short-term memory depends on neurons firing persistently or in brief, coordinated bursts.

The proposed cascade of oscillators consists of a ring oscillator and four toggle flip-flops connected in sequence. The novel oscillator and toggle are composed of three and six neurons, respectively. Their operation depends only on minimal properties of excitatory and inhibitory inputs.

The hypothesis that cascaded oscillators produce EEG phenomena implies that the distribution of EEG frequencies is determined by just two parameters, the mean (d) and standard deviation ({sigma}d) of the delay times of neurons that make up the initial oscillators in the cascades. For example, if d and {sigma}d are measured in milliseconds, the boundary separating the alpha and beta frequency bands is O_FD O_INLINEFIG[Formula]C_INLINEFIGC_FD

With 4 and 1.5 ms being the best available estimates for d and {sigma}d, respectively, this predicted boundary value is 14.9 Hz, which is within the range of commonly cited estimates obtained empirically from EEGs. Four and 1.5 ms also accurately predict the peaks and other boundaries of the five major EEG frequency bands.

Significance statementThe neuronal model proposed here implies several major aspects of electroencephalography. The matched periods of neural activity found in EEGs, and their wide distribution across the brain and across the frequency spectrum, follow from selective pressure for a biologically useful function: timing error avoidance for diverse brain functions in the trade-off between speed and accuracy. This activity can be achieved with a simple organization of synaptic connections and minimal neuron capabilities of excitation and inhibition. The multimodal distribution of EEG frequencies is an explicit function of the mean and variance of neuron delay times. The model suggests a relationship to epileptic seizures and a resolution to a short-term memory controversy. Two EEG characteristics make other models implausible.
Revisiting regulatory decoherence and phenotypic integration: accounting for temporal bias in co-expression analyses Environment can alter the degree of phenotypic variation and covariation, potentially influencing evolutionary trajectories. However, environment-driven changes in phenotypic variation remain understudied. In an effort to exploit the abundance of RNASequencing data now available, an increasing number of ecological studies rely on population-level correlation to characterize the plastic response of the entire transcriptome and to identify environmentally responsive molecular pathways. These studies are fundamentally interested in identifying groups of genes that respond in concert to environmental shifts. We show that population-level differential co-expression exhibits biases when capturing changes of regulatory activity and strength in rice plants responding to elevated temperature. One possible cause of this bias is regulatory saturation, the observation that detectable co-variance between a regulator and its target may be low as their transcript abundances are induced. This phenomenon appears to be particularly acute for rapid-onset environmental stressors. However, our results suggest that temporal correlations may be a reliable means to detect transient regulatory activity following rapid onset environmental perturbations such as temperature stress. Such temporal bias is likely to confound the studies of phenotypic integration, where high-order organismal traits are hypothesized to be more integrated with strong correlation under stressful conditions, while recent transcriptome studies exhibited weaker coexpression between genes under stressful conditions. Collectively, our results point to the need to account for the nuances of molecular interactions and the possibly confounding effects that these can introduce into conventional approaches to study transcriptome datasets. Correlational networking guides the discovery of cryptic natural product biosynthetic enzymes Bacterial natural product biosynthetic genes, canonically clustered, have been increasingly found to rely on cryptic enzymes encoded elsewhere in the genome for completion of biosynthesis. The study and application of lanthipeptides are frequently hindered by cryptic proteases required for final maturation. Here, we establish a global correlation network bridging the gap between lanthipeptide precursors and cryptic proteases. Applying our analysis to 161,954 bacterial genomes, we establish 5,209 correlations between precursors and cryptic proteases, with 91 prioritized. We use network predictions and co-expression analysis to reveal a previously missing protease for the maturation of class I lanthipeptide paenilan. We further discover widely distributed bacterial M16B metallopeptidases of previously unclear biological function as a new family of lanthipeptide proteases. We show the involvement of a pair of bifunctional M16B proteases in the production of novel class III lanthipeptides with high substrate specificity. Together, these results demonstrate the strength of our correlational networking approach to the discovery of cryptic lanthipeptide proteases and potentially other missing enzymes for natural products biosynthesis.

Graphical abstract O_FIG_DISPLAY_L [Figure 1] M_FIG_DISPLAY C_FIG_DISPLAY Imbalanced segregation of recombinant haplotypes in hybrid populations reveals inter- and intrachromosomal Dobzhansky-Muller incompatibilities Dobzhansky-Muller incompatibilities (DMIs) are a major component of reproductive isolation between species. DMIs imply negative epistasis, and are exposed when two diverged populations hybridize. Mapping the locations of DMIs has largely relied on classical genetic mapping. Approaches to date are hampered by low power and the challenge of identifying DMI loci on the same chromosome, because strong initial linkage of parental haplotypes weakens statistical tests. Here, we propose new statistics to infer negative epistasis from haplotype frequencies in hybrid populations. When two divergent populations hybridize, the variance in heterozygosity at two loci decreases faster with time at DMI loci than at random pairs of loci. When two populations hybridize at near-even admixture proportions, the deviation of the observed variance from its expectation becomes negative for the DMI pair. This negative deviation enables us to detect intermediate to strong negative epistasis both within and between chromosomes. In practice, the detection window in hybrid populations depends on the demographic scenario, the recombination rate, and the strength of epistasis. When the initial proportion of the two parental populations is uneven, only strong DMIs can be detected with our method unless migration prevents parental haplotypes from being lost. We use the new statistics to infer candidate DMIs from three hybrid populations of swordtail fish. We identify numerous new DMI candidates, some of which are inferred to interact with several loci within and between chromosomes. Moreover, we discuss our results in the context of an expected enrichment in intrachromosomal over interchromosomal DMIs.

Author SummaryGenetic incompatibility in the form of (Bateson-)Dobzhansky-Muller incompatibilities (DMIs) is an important component of reproductive isolation between species. However, the evolutionary role of DMIs during the process of speciation is contentious. DMIs occur when two or more genetic variants interact to reduce their carriers fitness. Once recombination combines incompatible variants in hybrids, selection acts to remove these variants from the population. One step towards addressing the evolutionary role of DMIs is to quantify the prevalence of DMIs in incipient and hybridizing species. Here, we present statistics that are sensitive to the resulting recombinant imbalance and that can indicate the location of DMIs in hybrid genomes in various demographic scenarios. We use simulations to show that the time window during which a DMI is detectable depends on its genomic location, the severity of the DMI, and the populations demography. Importantly, our statistic distinguishes genetic associations arising due to physical linkage from those arising due to gene interactions, which allows for the inference of both inter- and intrachromosomal DMIs. Applying our statistics to three hybrid populations of swordtail fish, we confirm previously known DMIs and identify new candidate incompatibilities. Medium-throughput zebrafish optogenetic platform identifies deficits in subsequent neural activity following brief early exposure to cannabidiol and Δ-9-tetrahydrocannabinol In the light of legislative changes and the widespread use of cannabis as a recreational and medicinal drug, delayed effects of cannabis upon brief exposure during embryonic development are of high interest as early pregnancies often go undetected. Here, zebrafish embryos were exposed to cannabidiol (CBD) and {Delta}-9-tetrahydrocannabinol l (THC) until the end of gastrulation (1-10 hours post-fertilization) and analyzed later in development (4-5 days post-fertilization). In order to measure neural activity, we implemented CaMPARI (Calcium-Modulated Photoactivatable Ratiometric Integrator) and optimized the protocol for a 96-well format complemented by locomotor analysis. Our results revealed that neural activity was decreased by CBD more than THC. At higher doses, both cannabinoids could dramatically reduce neural activity and locomotor activity. Interestingly, the decrease was more pronounced when CBD and THC were combined. At the receptor level, CBD-mediated reduction of locomotor activity was partially prevented using cannabinoid type 1 and 2 receptor inhibitors. Overall, we report that CBD toxicity occurs via two cannabinoid receptors and is synergistically enhanced by THC exposure to negatively impact neural activity late in larval development. Future studies are warranted to reveal other cannabinoids and receptors involved in this pathway to understand the subsequent health implications of cannabis consumption on fetal development. Identification and analysis of novel variants in SARS-COV-2 genomes isolated from the Kingdom of Bahrain The challenges imposed by the ongoing outbreak of severe acute respiratory syndrome coronavirus-2 affects every aspect of our modern world, ranging from our health to our socio-economic needs. Our existence highly depends on the vaccines availability, which demands in-depth research of the available strains and their mutations. In this work, we have analyzed all the available SERS-CoV2 genomes isolated from the Kingdom of Bahrain in terms of their variance and origin analysis. We have predicted various known and unique mutations in the SERS-CoV2 isolated from Bahrain. The complexity of the phylogenetic tree and dot plot representation of the strains mentioned above with other isolates of Asia indicates the versatility and multiple origins of Bahrains SERS-CoV2 isolates. We have also identified two high impact spike mutations from these strains which increase the virulence of SARS-CoV2. Our research could have a high impact on vaccine development and distinguishes the source of SERS-CoV2 in the Kingdom of Bahrain. SPIRO - the automated Petri plate imaging platform designed by biologists, for biologists The imaging of plant seedlings, fungal mycelia and bacterial colonies grown on Petri plates is common in phenotyping assays, and is typically done manually despite the procedures being time-consuming and laborious. The main reason for this is the still limited availability of existing automated phenotyping tools and facilities. Additionally, constructing a custom-made automated solution is a daunting task for most research groups specializing in biology.

Here, we describe SPIRO, the Smart Plate Imaging Robot, an automated platform that acquires time-lapse photos of up to four vertically oriented Petri plates in a single experiment. SPIRO was designed for biologists by biologists; thus, its assembly does not require experience in engineering or programming and its operation is sufficiently intuitive to be carried out without training. SPIRO has a small footprint optimal for fitting into standard incubators for plants and microbes, and is equipped with an LED light source for imaging in the dark, thus allowing acquisition of photos under optimal growth conditions. SPIROs web-based user interface allows setting up experiments and downloading data remotely, without interfering with samples growth. SPIROs excellent image quality is optimal for automated image processing, which we demonstrate with two semi-automated assays for analysis of commonly used phenotypic traits: seed germination and root growth.

Moreover, the robot can be easily customized for a specific use, as all information about SPIRO, including the models for 3D-printed structural components, control software, and scripts for image analysis, are released under permissive open-source licenses. The interplay of chromatin phase separation and lamina interactions in nuclear organization The genetic material of eukaryotes is segregated into transcriptionally active euchromatin and silent heterochromatin compartments. The spatial arrangement of chromatin compartments evolves over the course of cellular life in a process that remains poorly understood. The latest nuclear imaging experiments reveal a number of dynamical signatures of chromatin which are reminiscent of active multi-phase liquids. This includes the observations of viscoelastic response, coherent motions, Ostwald ripening and coalescence of chromatin compartments. There is also growing evidence that liquid-liquid phase separation of protein and nucleic acid components is the underlying mechanism for the dynamical behavior of chromatin. In order to dissect the organizational and dynamical implications of chromatins liquid behavior, we have devised a phenomenological field-theoretic model of nucleus as a multi-phase condensate of liquid chromatin types. Employing the liquid chromatin model of Drosophila nucleus, we have carried out an extensive set of simulations with an objective to shed light on the dynamics and chromatin patterning observed in the latest nuclear imaging experiments. Our simulations reveal the emergence of experimentally detected mesoscale chromatin channels and spheroidal droplets which arise from the dynamic interplay of chromatin type to type interactions and intermingling of chromosomal territories. We also quantitatively reproduce coherent motions of chromatin domains observed in displacement correlation spectroscopy measurements which are explained within the framework of our model by phase separation of chromatin types operating within constrained intra-chromosomal and inter-chromosomal boundaries. Finally, we illuminate the role of heterochromatin-lamina interactions in the nuclear organization by showing that these interactions enhance the mobility of euchromatin and indirectly introduce correlated motions of heterochromatin droplets.

SIGNIFICANCE STATEMENTThe latest super-resolution imaging experiments have revealed a surprisingly dynamic and stochastic nature of chromatin in eukaryotic nuclei which is reminiscent of multi-phase fluid behavior. As a result, the concepts from the theory of complex fluids such as phase separation, viscoelasticity, and droplet nucleation have found widespread utility in understanding salient features of nuclear organization. In order to understand and disentangle the complex interplay of forces that contribute to the emergent patterns of organization and dynamics, we have devised a phenomenological field-theoretic model of nucleus as a multi-phase condensate of liquid chromatin types. Armed with a mesoscopic model of nuclear chromatin, we have shed light on the distinct dynamical and structural contributions of chromatin type interactions intermingling of chromosomal territories and lamina binding. We also shed light on the dynamical heterogeneity and coherent motions of chromatin domains which are fully captured by an interplay of micro-phase separation of chromatin types and lamina binding. Fungicide resistance characterised across seven modes of action in a Botrytis cinerea isolated from Australian vineyards BACKGROUNDBotrytis bunch rot, caused by Botrytis cinerea, is an economically important disease of grapes in Australia and across grape growing regions worldwide. Control of this disease relies on canopy management and the application of fungicides. Fungicide application can lead to the selection of resistant B. cinerea populations, which has an adverse effect on management of the disease. Characterising the distribution and severity of resistant B. cinerea populations is needed to inform resistance management strategies.

RESULTSIn this study, 725 isolates were sampled from 75 Australian vineyards during 2013 - 2016 and were screened against seven fungicides with different modes of action (MOAs). The resistance frequencies for azoxystrobin, boscalid, fenhexamid, fludioxonil, iprodione, pyrimethanil and tebuconazole were 5, 2.8, 2.1, 6.2, 11.6, 7.7 and 2.9% respectively. Nearly half of the resistant isolates (43.8%) were resistant to more than one of the fungicides tested. The frequency of vineyards with at least one isolate simultaneously resistant to 1, 2, 3, 4 or 5 fungicides was 19.5, 7.8, 6.5, 10.4 and 2.6%. Resistance was associated with previously published genotypes in CytB (G143A), SdhB (H272R/Y), Erg27 (F412S), Mrr1 (D354Y), Bos1 (I365S, N373S + Q369P, I365S + D757N) and Pos5 (V273I, P319A, L412F/V). Novel genotypes were also described in Mrr1 (S611N, D616G) Pos5 (V273L) and Cyp51 (P347S). Expression analysis was used to characterise fludioxonil resistant isolates exhibiting overexpression (6.3-9.6-fold) of the ABC transporter gene AtrB (MDR1 phenotype).

CONCLUSIONResistance frequencies were lower when compared to most previously published surveys of B. cinerea resistance in grape and other crops. Nevertheless, continued monitoring of critical MOAs used in Australian vineyards is recommended. Population differences in the length and early-life dynamics of telomeres among European pied flycatchers Telomere length and shortening rate are increasingly used as biomarkers for long-term costs in ecological and evolutionary studies because of their relationships with survival and fitness. Both early-life conditions and growth, and later-life stressors can create variation in telomere shortening rate. Studies on between- population telomere length and dynamics are scarce, despite the expectation that populations exposed to varying environmental constraints would present divergent telomere length patterns. The pied flycatcher (Ficedula hypoleuca) is a passerine bird breeding across Eurasia (from Spain to western Siberia) and migrating through the Iberian Peninsula to spend the non-breeding period in sub-Saharan Africa. Thus, different populations show marked differences in migration distance. We studied the large-scale variation of telomere length and early-life dynamics in the pied flycatcher by comparing six European populations across a north- south gradient (Finland, Estonia, England, and Spain) predicting negative effect of migration distance on adult telomere length, and of nestling growth on nestling telomere dynamics. There were clear population differences in telomere length, with English birds from mid-latitudes having the longest telomeres. Telomere length did not thus show consistent latitudinal variation and was not linearly linked to differences in migration distance. Early-life telomere shortening rate tended to vary between populations. Fast growth was associated with shorter telomeres in the early life, but faster nestling growth affected telomeres more negatively in northern than southern populations. While the sources of between-population differences in telomere-related biology remain to be more intensively studied, our study illustrates the need to expand telomere studies at the between-population level. Parallel Characterization of cis-Regulatory Elements for Multiple Genes Using CRISPRpath Current pooled CRISPR screens for cis-regulatory elements (CREs) can only accommodate one gene based on its expression level. Here, we describe CRISPRpath, a scalable screening strategy for parallelly characterizing CREs of genes linked to the same biological pathway and converging phenotypes. We demonstrate the ability of CRISPRpath for simultaneously identifying functional enhancers of six genes in the 6-thioguanine-induced DNA mismatch repair pathway using both CRISPR interference (CRISPRi) and CRISPR nuclease (CRISPRn) approaches. 60% of the identified enhancers are known promoters with distinct epigenomic features compared to other active promoters, including increased chromatin accessibility and interactivity. Furthermore, by imposing different levels of selection pressure, CRISPRpath can distinguish enhancers exerting strong impact on gene expression from those exerting weak impact. Our results offer a nuanced view of cis-regulation and demonstrate that CRISPRpath can be leveraged for understanding the complex gene regulatory program beyond transcriptional output at scale. 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. HYBRID REACTION-DIFFUSION AND CLOCK-AND-WAVEFRONT MODEL FOR THE ARREST OF OSCILLATIONS IN THE SOMITOGENESIS SEGMENTATION CLOCK AO_SCPLOWBSTRACTC_SCPLOWThe clock and wavefront paradigm is arguably the most widely accepted model for explaining the embryonic process of somitogenesis. According to this model, somitogenesis is based upon the interaction between a genetic oscillator, known as segmentation clock, and a differentiation wavefront, which provides the positional information indicating where each pair of somites is formed. Shortly after the clock and wavefront paradigm was introduced, Meinhardt presented a conceptually different mathematical model for morphogenesis in general, and somitogenesis in particular. Recently, Cotterell et al. rediscovered an equivalent model by systematically enumerating and studying small networks performing segmentation. Cotterell et al. called it a progressive oscillatory reaction-diffusion (PORD) model. In the Meinhardt-PORD model, somitogenesis is driven by short-range interactions and the posterior movement of the front is a local, emergent phenomenon, which is not controlled by global positional information. With this model, it is possible to explain some experimental observations that are incompatible with the clock and wavefront model. However the Meinhardt-PORD model has some important disadvantages of its own. Namely, it is quite sensitive to fluctuations and depends on very specific initial conditions (which are not biologically realistic). In this work, we propose an equivalent Meinhardt-PORD model, and then amend it to couple it with a wavefront consisting of a receding morphogen gradient. By doing so, we get a hybrid model between the Meinhardt-PORD and the clock-and-wavefront ones, which overcomes most of the deficiencies of the two originating models.

Somitogenesis, the process by which somites are formed, is an essential developmental stage in many vertebrates. This process occurs with a strikingly regular periodicity, that is preserved among embryos of a single species. The clock and wavefront paradigm is arguably the most widely accepted model for explaining somitogenesis. However, it is incapable of explaining some experimental facts, like the appearance of somites in the absence of an external wavefront (i.e. a receding morphogen gradient). Shortly after the clock and wavefront paradigm was introduced, Meinhardt presented a conceptually different mathematical model for morphogenesis in general, and somitogenesis in particular. Recently, Cotterell et al. rediscovered an equivalent model by systematically enumerating and studying small networks performing segmentation, and called it a progressive oscillatory reaction-diffusion (PORD) model. The Meinhardt-PORD model tackles some of the deficiencies of the clock and wave-front models, but it has some serious issues of its own. In the present work, we introduce an equivalent Meinhardt-PORD model, and then amend it to couple it with a receding morphogen gradient. By doing so, we get a hybrid model that incorporates characteristics of the Meinhardt-PORD and clock-and-wavefront models. We show that this hybrid model under-goes a bifurcation, from a stable to an unstable limit cycle, as the value of the parameter accounting for a background regulatory input (associated to the receding morphogen gradient) decreases. This bifurcation allows the model to explain why somites can form in the absence of an external wavefront, reassesses the role of the receding morphogen gradient as a conductor for somitogenesis, and makes the model behavior robust to random fluctuations, as well as independent from specific initial conditions (the latter, are two of the weak points of the Meinhardt-PORD model). We argue that this findings provide convincing evidence that reaction-diffusion and positional information (receding morphogen gradient) mechanisms could work together in somitogenesis. Propelling and perturbing appendages together facilitate strenuous ground self-righting Terrestrial animals must self-right when overturned on the ground, but this locomotor task is strenuous. To do so, the discoid cockroach often pushes its wings against the ground to begin a somersault which rarely succeeds. As it repeatedly attempts this, the animal probabilistically rolls to the side to self- right. During winged self-righting, the animal flails its legs vigorously. Here, we studied whether wing opening and leg flailing together facilitate strenuous ground self-righting. Adding mass to increase hind leg flailing kinetic energy increased the animals self-righting probability. We then developed a robot with similar strenuous self-righting behavior and used it as a physical model for systematic experiments. The robots self-righting probability increased with wing opening and leg flailing amplitudes. A potential energy landscape model revealed that, although wing opening did not generate sufficient kinetic energy to overcome the high pitch potential energy barrier to somersault, it reduced the barrier for rolling, facilitating the small kinetic energy from leg flailing to probabilistically overcome it to self-right. The model also revealed that the stereotyped body motion during self-righting emerged from physical interaction of the body and appendages with the ground. Our work demonstrated the usefulness of potential energy landscape for modeling self-righting transitions.

Impact statementAnimals can, and robots should, use different types of appendages together to propel and perturb themselves to self-right when overturned, a strenuous yet crucial locomotor task. An integrative analysis to predict the Active Ingredients and explore Polypharmacological Mechanisms of Orthosiphon stamineus Benth. for human health BackgroundOrthosiphon stamineus Benth. was a dietary supplement or a traditional folk herb with widespread clinical application, but it was lack of comprehensive understanding of its active ingredients and polypharmacological mechanisms. The aim of this work was to systematically study its natural compounds and the molecular mechanisms of the O. stamineus via network pharmacology.

MethodsCompounds from O. stamineus were collected by literature retrieval and evaluated by SwissADME with the physicochemical property (ADMET model) and the likelihood for a natural medicine. The connection of active ingredients and target genes was built and confirmed by Cytoscape and AutoDock vina. Then, Gene Ontology (GO) enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were separately performed to obtain a more in-depth understanding of O. stamineus. Finally, the relationship among active ingredients, targets, and diseases was built to clarify the polypharmacological mechanisms and found relevant active substances for further drug discovery.

ResultsA total of 159 compounds of O. stamineus were collected, and 22 potentially active ingredients were screened out. The Ingredient-Target Interaction Network was built with 12 flavonoids, 3 diterpenes, 3 phenols, and 4 volatile oils, and 65 targets. The Docking analysis indicated that the ingredient-target interaction network was reliable; most ligand-receptor had a strong binding affinity (lowest binding energy: -6.9 kcal/mol). After pathway analysis, 185 significant biological processes and 36 signal pathways were found, and the ingredient-target-disease network of O. stamineus was constructed for polypharmacological mechanisms.

ConclusionOur study clarified the polypharmacological mechanisms via the relationship among active ingredients, targets, and diseases and provided better guidance for subsequent experiments and potential active ingredients for drug discovery or health promotion. A MALDI-TOF assay identifies nilotinib as an inhibitor of inflammation in acute myeloid leukaemia Inflammatory responses are important in cancer, particularly in the context of monocyte-rich aggressive myeloid neoplasm. We developed a label-free cellular phenotypic drug discovery assay to identify anti-inflammatory drugs in human monocytes derived from acute myeloid leukaemia (AML), by tracking several biological features ionizing from only 2,500 cells using matrix-assisted laser desorption/ionization-time of flight (MALDI-TOF) mass spectrometry. A proof-of-concept screen showed that the BCR-ABL inhibitor nilotinib, but not the structurally similar imatinib, blocks inflammatory responses. In order to identify the cellular (off-)targets of nilotinib, we performed thermal proteome profiling (TPP). Unlike imatinib, nilotinib and other later generation BCR-ABL inhibitors inhibit the p38-MK2/3 signalling axis which suppressed the expression of inflammatory cytokines, cell adhesion and innate immunity markers in activated human monocytes derived from AML. Thus, our study provides a tool for the discovery of new anti-inflammatory drugs, which could contribute to the treatment of inflammation in myeloid neoplasms and other diseases.

Key PointsLabel-free cell-based assay identifies new anti-inflammatory drugs using MALDI-TOF MS. Nilotinib reduces inflammation by inhibition of MAPK14-MK2/3 signalling axis in AML. Data-Driven Network Models for Genetic Circuits From Time-Series Data with Incomplete Measurements Synthetic gene networks are frequently conceptualized and visualized as static graphs. This view of biological programming stands in stark contrast to the transient nature of biomolecular interaction, which is frequently enacted by labile molecules that are often unmeasured. Thus, the network topology and dynamics of synthetic gene networks can be difficult to verify in vivo or in vitro, due to the presence of unmeasured biological states. Here we introduce the dynamical structure function as a new mesoscopic, data-driven class of models to describe gene networks with incomplete measurements. We introduce a network reconstruction algorithm and a code base for reconstructing the dynamical structure function from data, to enable discovery and visualization of graphical relationships in a genetic circuit diagram as time-dependent functions rather than static, unknown weights. We prove a theorem, showing that dynamical structure functions can provide a data-driven estimate of the size of crosstalk fluctuations from an idealized model. We illustrate this idea with numerical examples. Finally, we show how data-driven estimation of dynamical structure functions can explain failure modes in two experimentally implemented genetic circuits, a historical genetic circuit and a new E. coli based transcriptional event detector. Heteologous saRNA-Prime, DNA Dual-Antigen-Boost SARS-CoV-2 Vaccination Elicits Robust Cellular Immunogenicity and Cross-Variant Neutralizing Antibodies We assessed if immune responses are enhanced in CD-1 mice by heterologous vaccination with two different nucleic acid-based COVID-19 vaccines: a next-generation human adenovirus serotype 5 (hAd5)-vectored dual-antigen spike (S) and nucleocapsid (N) vaccine (AdS+N) and a self-amplifying and -adjuvanted S RNA vaccine (SASA S) delivered by a nano-lipid carrier. The AdS+N vaccine encodes S modified with a fusion motif to increase cell-surface expression. The N antigen is modified with an Enhanced T-cell Stimulation Domain (N-ETSD) to direct N to the endosomal/lysosomal compartment and increase MHC class I and II stimulation potential. The S sequence in the SASA S vaccine comprises the D614G mutation, two prolines to stabilize S in the prefusion conformation, and 3 glutamines in the furin cleavage region to increase cross-reactivity across variants. CD-1 mice received vaccination by homologous and heterologous prime > boost combinations. Humoral responses to S were the highest with any regimen including the SASA S vaccine, and IgG bound to wild type and Delta (B.1.617.2) variant S1 at similar levels. An AdS+N boost of an SASA S prime particularly enhanced both CD4+ and CD8+ T-cell responses to both wild type and Delta S peptides relative to all other vaccine regimens. Sera from mice receiving SASA S homologous or heterologous vaccination were found to be highly neutralizing of all pseudovirus strains tested: Wuhan, Beta, Delta, and Omicron strain. The findings here support the clinical testing of heterologous vaccination by an SASA S > AdS+N regimen to provide increased protection against emerging SARS-CoV-2 variants. A Novel Proposal for an Index for Regional Cerebral Perfusion Pressure - A Theoretical Approach Using Fluid Dynamics Cerebral blood flow (CBF) / cerebral blood volume (CBV) ratio derived by [15O] H2O/ CO2 and CO positron emission tomography (PET) examination has been used empirically as an index for cerebral perfusion pressure (CPP). However, this equation lacks a theoretical basis and furthermore, has not been confirmed to be proportionate to CPP, as direct measurement of local CPP is not practical.

We have developed a new index for CPP using the Poiseuille equation based on a simple model. Our model suggests that CBF/CBV2 is proportionate to CPP and that it is mathematically a more accurate index than CBF/CBV. iCn3D: From Web-based 3D Viewer to Structure Analysis Tool in Batch Mode iCn3D was initially developed as a web-based 3D molecular viewer. It then evolved from visualization into a full-featured interactive structural analysis software. It became a collaborative research instrument through the sharing of permanent, shortened URLs that encapsulate not only annotated visual molecular scenes, but also all underlying data and analysis scripts in a FAIR manner. More recently, with the growth of structural databases, the need to analyze large structural datasets systematically led us to use Python scripts and convert the code to be used in Node.js scripts. We showed a few examples of Python scripts at https://github.com/ncbi/icn3d/tree/master/icn3dpython to export secondary structures or PNG images from iCn3D. Users just need to replace the URL in the Python scripts to export other annotations from iCn3D. Furthermore, any interactive iCn3D feature can be converted into a Node.js script to be run in batch mode, enabling an interactive analysis performed on one or a handful of protein complexes to be scaled up to analysis features of large ensembles of structures. Currently available Node.js analysis scripts examples are available at https://github.com/ncbi/icn3d/tree/master/icn3dnode. This development will enable ensemble analyses on growing structural databases such as AlphaFold or RoseTTAFold on one hand and Electron Microscopy on the other. In this paper, we also review new features such as DelPhi electrostatic potential, 3D view of mutations, alignment of multiple chains, assembly of multiple structures by realignment, dynamic symmetry calculation, 2D cartoons at different levels, interactive contact maps, and use of iCn3D in Jupyter Notebook as described at https://pypi.org/project/icn3dpy. Creatine transporter deficient rat model show motor dysfunction, cerebellar alterations and muscle creatine deficiency without muscle atrophy Creatine (Cr) is a nitrogenous organic acid and plays roles as fast phosphate energy buffer to replenish ATP, osmolyte, antioxidant, neuromodulator, and as a compound with anabolic and ergogenic properties in muscle. Cr is taken from the diet or endogenously synthetized by the enzymes AGAT and GAMT, and specifically taken up by the transporter SLC6A8. Loss-of-function mutations in the genes encoding for the enzymes or the transporter cause Cerebral Creatine Deficiency Syndromes (CCDS). CCDS are characterized by brain Cr deficiency, intellectual disability with severe speech delay, behavioral troubles, epilepsy and motor dysfunction. Among CCDS, the X-linked Cr transporter deficiency (CTD) is the most prevalent with no efficient treatment so far. Different animal models of CTD show reduced brain Cr levels, cognitive deficiencies and together they cover other traits similar to those of patients. However, motor function was poorly explored in CTD models and some controversies in the phenotype exist in comparison with CTD patients. Our recently described Slc6a8Y389C knock-in (KI) rat model of CTD showed mild impaired motor function, morphological alterations in cerebellum, reduced muscular mass, Cr deficiency and increased guanidinoacetate content in muscle, although no consistent signs of muscle atrophy. Our results indicate that such motor dysfunction co-occurred with both nervous and muscle dysfunction, suggesting that muscle strength and performance as well as neuronal connectivity might be affected by this Cr deficiency in muscle and brain. Itaconate and its derivatives repress C2C12 myogenesis A Krebs cycle intermediate metabolite, itaconate, has gained attention as a potential antimicrobial and autoimmune disease treatment due to its anti-inflammatory effects. While itaconate and its derivatives pose an attractive therapeutic option for the treatment of inflammatory diseases, the effects outside the immune system still remain limited, particularly in the muscle. Therefore, we endeavored to determine if itaconate signaling impacts muscle differentiation. Utilizing the well-established C2C12 model of in vitro myogenesis, we evaluated the effects of itaconate and its derivatives on transcriptional and protein markers of muscle differentiation as well as mitochondrial function. We found itaconate and the derivatives dimethyl itaconate and 4-octyl itaconate disrupt differentiation media-induced myogenesis. A primary biological effect of itaconate is a succinate dehydrogenase (SDH) inhibitor. We find the SDH inhibitors dimethyl malonate and harzianopyridone phenocopie the anti-myogenic effects of itaconate. Furthermore, we find treatment with exogenous succinate results in blunted myogenesis. Together our data indicate itaconate and its derivatives interfere with in vitro myogenesis, potentially through inhibition of SDH and subsequent succinate accumulation. More importantly, our findings suggest the therapeutic potential of itaconate and its derivatives could be limited due to deleterious effects on myogenesis. Liver ductal organoids reconstruct intrahepatic biliary trees in decellularized liver grafts Three-dimensional scaffolds decellularized from native organs are a promising technique to establish engineered liver grafts and overcome the current shortage of donor organs. However, limited sources of bile duct cells and inappropriate cell distribution in bioengineered liver grafts have hindered their practical application. Organoid technology is anticipated to be an excellent tool for the advancement of regenerative medicine. In the present study, we reconstructed intrahepatic bile ducts in a rat decellularized liver graft by recellularization with liver ductal organoids. Using an ex vivo perfusion culture system, we demonstrated the biliary characteristics of repopulated mouse liver organoids, which maintained bile duct markers and reconstructed biliary tree-like networks with luminal structures. We also established a method for the co-recellularization with engineered bile ducts and primary hepatocytes, revealing the appropriate cell distribution to mimic the native liver. We then utilized this model in human organoids to demonstrate the reconstructed bile ducts. Our results show that liver ductal organoids are a potential cell source for bile ducts from bioengineered liver grafts using three-dimensional scaffolds. 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. Lowland plant arrival in alpine ecosystems facilitates a decrease in soil carbon content under experimental climate warming Climate warming is releasing carbon from soils around the world1-3, constituting a positive climate feedback. Warming is also causing species to expand their ranges into new ecosystems4-9. Yet, in most ecosystems, whether range expanding species will amplify or buffer expected soil carbon loss is unknown10. Here we used two whole-community transplant experiments and a follow-up glasshouse experiment to determine whether the establishment of herbaceous lowland plants in alpine ecosystems influences soil carbon content under warming. We found that warming (transplantation to low elevation) led to a negligible decrease in alpine soil carbon content, but its effects became significant and 52% {+/-} 31% (mean {+/-} 95% CIs) larger after lowland plants were introduced at low density into the ecosystem. We present evidence that decreases in soil carbon content likely occurred via lowland plants increasing rates of root exudation, soil microbial respiration and CO2 release under warming. Our findings suggest that warming-induced range expansions of herbaceous plants have the potential to alter climate feedbacks from this system, and that plant range expansions among herbaceous communities may be an overlooked mediator of warming effects on carbon dynamics.