TemplateFlow: FAIR-sharing of multi-scale, multi-species brain models Reference anatomies of the brain and corresponding atlases play a central role in experimental neuroimaging workflows and are the foundation for reporting standardized results. The choice of such references --i.e., templates-- and atlases is one relevant source of methodological variability across studies, which has recently been brought to attention as an important challenge to reproducibility in neuroscience. TemplateFlow is a publicly available framework for human and nonhuman brain models. The framework combines an open database with software for access, management, and vetting, allowing scientists to distribute their resources under FAIR --findable, accessible, interoperable, reusable-- principles. TemplateFlow supports a multifaceted insight into brains across species, and enables multiverse analyses testing whether results generalize across standard references, scales, and in the long term, species, thereby contributing to increasing the reliability of neuroimaging results.
Rostral Anterior Cingulate Activations inversely relate to Reward Payoff Maximation & predict Depressed Mood Choice selection strategies and decision making are typically investigated using multiple-choice gambling paradigms that require participants to maximize reward payoff. However, research shows that performance in such paradigms suffers from individual biases towards the frequency of gains to choose smaller local gains over larger longer term gain, also referred to as melioration. Here, we developed a simple two-choice reward task, implemented in 186 healthy human adult subjects across the adult lifespan to understand the behavioral, computational, and neural bases of payoff maximization versus melioration. The observed reward choice behavior on this task was best explained by a reinforcement learning model of differential future reward prediction. Simultaneously recorded and source-localized electroencephalography (EEG) showed that diminished theta-band activations in the right rostral anterior cingulate cortex (rACC) correspond to greater reward payoff maximization, specifically during the presentation of cumulative reward information at the end of each task trial. Notably, these activations (greater rACC theta) predicted depressed mood symptoms, thereby showcasing a reward processing marker of potential clinical utility.

Significance StatementThis study presents cognitive, computational and neural (EEG-based) analyses of a rapid reward-based decision-making task. The research has the following three highlights. 1) It teases apart two core aspects of reward processing, i.e. long term expected value maximization versus immediate gain frequency melioration based choice behavior. 2) It models reinforcement learning based behavioral differences between individuals showing that observed performance is best explained by differential extents of reward prediction. 3) It investigates neural correlates in 186 healthy human subjects across the adult lifespan, revealing specific theta band cortical source activations in right rostral anterior cingulate as correlates for maximization that further predict depressed mood across subjects.
The cardiac TRPA1 channel drives calcium-mediated mechano-arrhythmogenesis SUMMARY PARAGRAPHMaintenance of cardiac function involves a regulatory loop in which electrical excitation causes the heart to contract through excitation-contraction coupling (ECC),1 and the mechanical state of the heart directly affects its electrical activity through mechano-electric coupling (MEC).2 However, in pathological states such as acute ischaemia that alter early or late electro-mechanical coordination (i.e., disturbances in ECC or repolarisation-relaxation coupling, RRC), MEC may contribute to the initiation and / or sustenance of arrhythmias (mechano-arrhythmogenesis).3 The molecular identity of specific factor(s) underlying mechano-arrhythmogenesis in acute ischaemia, however, remain undefined.4 By rapid stretch of rabbit single left ventricular cardiomyocytes, we show that upon ATP-sensitive potassium channel-induced alterations of RRC, overall vulnerability to mechano-arrhythmogenesis is increased, with mechano-sensitive5-11 transient receptor potential kinase ankyrin 1 (TRPA1) channels12 acting as the molecular driver through a Ca2+-mediated mechanism. Specifically, TRPA1 activation drives stretch-induced excitation and creates a substrate for self-sustained arrhythmias, which are maintained by increased cytosolic free Ca2+ concentration ([Ca2+]i) and spontaneous [Ca2+]i oscillations. This TRPA1-dependent mechano-arrhythmogenesis involves microtubules, and can be prevented by block of TRPA1 or buffering of [Ca2+]i. Thus, in cardiac pathologies with disturbed RRC dynamics and / or augmented TRPA1 activity, TRPA1 may represent an anti-arrhythmic target with untapped therapeutic potential.13-17 SPOT-Contact-Single: Improving Single-Sequence-Based Prediction of Protein Contact Map using a Transformer Language Model MotivationAccurate prediction of protein contact-map is essential for accurate protein structure and function prediction. As a result, many methods have been developed for protein contact map prediction. However, most methods rely on protein-sequence-evolutionary information, which may not exist for many proteins due to lack of naturally occurring homologous sequences. Moreover, generating evolutionary profiles is computationally intensive. Here, we developed a contact-map predictor utilizing the output of a pre-trained language model ESM-1b as an input along with a large training set and an ensemble of residual neural networks.

ResultsWe showed that the proposed method makes a significant improvement over a single-sequence-based predictor SSCpred with 15% improvement in the F1-score for the independent CASP14-FM test set. It also outperforms evolutionary-profile-based methods TrRosetta and SPOT-Contact with 48.7% and 48.5% respective improvement in the F1-score on the proteins without homologs (Neff=1) in the independent SPOT-2018 set. The new method provides a much faster and reasonably accurate alternative to evolution-based methods, useful for large-scale prediction.

AvailabilityStand-alone-version of SPOT-Contact-Single is available at https://github.com/jas-preet/SPOT-Contact-Single. Direct prediction can also be made at https://sparks-lab.org/server/spot-contact-single. The datasets used in this research can also be downloaded from the GitHub.

Contactjaspreetsingh2@griffithuni.edu.au, k.paliwal@griffith.edu.au, and zhouyq@szbl.ac.cn

Supplementary informationSupplementary data are available at Bioinformatics online. Mutagenic distinction between the receptor-binding and fusion subunits of the SARS-CoV-2 spike glycoprotein and its upshot We observe that a residue R of the spike glycoprotein of SARS-CoV-2 which has mutated in one or more of the current Variants of Concern or Interest or under Monitoring rarely participates in a backbone hydrogen bond if R lies in the S1 subunit and usually participates in one if R lies in the S2 subunit. A partial explanation for this based upon free energy is explored as a potentially general principle in the mutagenesis of viral glycoproteins. This observation could help target future vaccine cargos for the evolving coronavirus as well as more generally. A study of the Delta and Omicron variants suggests that Delta was an energetically necessary intermediary in the evolution from Wuhan-Hu-1 to Omicron. Loss of histone methyltransferase SETD1B in oogenesis results in the redistribution of genomic histone 3 lysine 4 trimethylation Histone 3 lysine 4 trimethylation (H3K4me3) is an epigenetic mark found at gene promoters and CpG islands. H3K4me3 is essential for mammalian development, yet mechanisms underlying its genomic targeting are poorly understood. H3K4me3 methyltransferases SETD1B and MLL2 are essential for oogenesis. We investigated changes in H3K4me3 in Setd1b conditional knockout (cKO) oocytes using ultra-low input ChIP-seq, with comparisons to DNA methylation and gene expression analyses. H3K4me3 was redistributed in Setd1b cKO oocytes showing losses at active gene promoters associated with downregulated gene expression. Remarkably, many regions also gained H3K4me3, in particular those that were DNA hypomethylated, transcriptionally inactive and CpG-rich, which are hallmarks of MLL2 targets. Consequently, loss of SETD1B disrupts the balance between MLL2 and de novo DNA methyltransferases in determining the epigenetic landscape during oogenesis. Our work reveals two distinct, complementary mechanisms of genomic targeting of H3K4me3 in oogenesis, with SETD1B linked to gene expression and MLL2 to CpG content.

Graphical AbstractIn oogenesis, SETD1B and CXXC1 target H3K4me3 to actively transcribed gene promoters, while MLL2 targets transcriptionally inactive regions based on underlying CpG composition (upper panel). When SETD1B is ablated, H3K4me3 is lost at a subset of active promoters, resulting in downregulation of transcription (lower panel). Loss of SETD1B alters the activity of MLL2, permitting MLL2 to deposit H3K4me3 at CpG-rich regions, many of which should otherwise be DNA methylated. Thus, it is evident that MLL2 and de novo DNMTs compete for genomic occupancy late in oogenesis, and loss of SETD1B disrupts the balance of these mechanisms.

O_FIG_DISPLAY_L [Figure 1] M_FIG_DISPLAY C_FIG_DISPLAY A Novel Organoid Model of In Vitro Spermatogenesis Using Human Induced Pluripotent Stem Cells Infertility is thought to be caused by genetic mutations and dysfunction in the cellular niche where spermatogenesis takes place. An understanding of the specialized cellular processes which drive spermatogenesis is needed to develop treatments; however, the development of in vitro systems to study these cells has been hindered by our reliance on rarely available human testicular tissues for research. Human induced pluripotent stem cells (hiPSCs) can be used to derive human testicular-like cells, and thus provide an avenue for the development of in vitro testicular model systems. Therefore, this study set out to engineer a human testicular tissue model using hiPSCs for the first time. We demonstrate the ability of hiPSC-derived testicular cells to self-organize and mature into testicular-like tissues using organoid culture. Moreover, we show that hiPSC-derived testicular organoids promote testicular somatic cell maturation and spermatogenesis up to the post-meiotic spermatid stage. These hiPSC-derived testicular organoids have the potential to replace rarely available primary testicular tissues to further infertility research in an in vitro setting. The multi-level phenotypic impact of synonymous substitutions: heterologous gene expression in human cells. Redundancy in the genetic code allows for differences in transcription and/or translation efficiency between sequences carrying synonymous polymorphisms, potentially leading to phenotypic changes. It is commonly admitted that the evolution codon usage bias ("CUB", the over-representation of certain codons in a genome, a gene or in positions along a gene) are driven by a combination of neutral and selective processes, but their relative contribution is a matter of debate, especially in mammals. Particularly, integrative studies quantifying the phenotypic impact of CUB at different molecular and cellular levels are lacking. Here we report a multiscale analysis of the effects of synonymous codon recoding during heterologous gene expression in human cells.

Six synonymous versions of the shble antibiotic resistance gene were generated, fused to a fluorescent reporter, and independently expressed in HEK293 cells. Multiscale phenotype was assessed by means of: i) mRNA-to-DNA and protein-to-mRNA ratios for each shble version; ii) cellular fluorescence, using flow cytometry, as a proxy for single-cell level expression; and iii) real-time cell proliferation in absence or presence of antibiotic, as a proxy for the cell fitness.

We show that differences in CUB strongly impact the molecular and the cellular phenotype: i) they result in large differences in mRNA and in protein levels, as well in mRNA-to-protein ratio; ii) they introduce splicing events not predicted by current algorithms; iii) they lead to reproducible phenotypic heterogeneity; iv) they lead to a trade-off between the benefit of antibiotic resistance and the burden of heterologous expression.

We interpret that CUB modulate mRNA availability and suitability for translation in human cells, leading to differences in protein levels and eventually eliciting phenotypic differences.

AUTHOR SUMMARYThe genetic code is redundant, with several codons encoding for the same amino acid. These synonymous codons are not used with equal frequencies. Instead, codon usage bias (CUB) varies between species, genes, and even positions along a gene. At each of these levels, CUB are shaped by the overall balance between mutational biases and selection forces. To shed light on the molecular mechanisms underlying molecular and organism phenotypes, integrative studies quantifying the phenotypic impact of CUB at different levels of biological integration are necessary. Here, we monitored the multiscale changes induced by six synonymous versions of an antibiotic resistance gene independently expressed in a human cell line. We show that: 1. both mRNA levels, protein levels, and their ratios are affected by CUB; 2. potential effects on mRNA (splicing and availability for translation) seem to constitute the main level of action; 3. cell fitness is severely impacted by a trade-off between the burden of heterologous expression and the benefits of antibiotic resistance. This integrative study provides new insights on translation regulation and the associated phenotypic impact in human cells, associated to CU. 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. Phylogenetic reconstruction of ancestral ecological networks through time for pierid butterflies and their host plants The study of herbivorous insects underpins much of the theory that concerns the evolution of species interactions. In particular, Pieridae butterflies and their host plants have served as a model system for studying evolutionary arms-races. To learn more about the coevolution of these two clades, we reconstructed ancestral ecological networks using stochastic mappings that were generated by a phylogenetic model of host-repertoire evolution. We then measured if, when, and how two ecologically important structural features of the ancestral networks (modularity and nestedness) evolved over time. Our study shows that as pierids gained new hosts and formed new modules, a subset of them retained or recolonized the ancestral host(s), preserving connectivity to the original modules. Together, host-range expansions and recolonizations promoted a phase transition in network structure. Our results demonstrate the power of combining network analysis with Bayesian inference of host-repertoire evolution to understand changes in complex species interactions over time. Construction of a new chromosome-scale, long-read reference genome assembly for the Syrian hamster, Mesocricetus auratus BackgroundThe Syrian hamster (Mesocricetus auratus) has been suggested as a useful mammalian model for a variety of diseases and infections, including infection with respiratory viruses such as SARS-CoV-2. The MesAur1.0 genome assembly was generated in 2013 using whole-genome shotgun sequencing with short-read sequence data. Current more advanced sequencing technologies and assembly methods now permit the generation of near-complete genome assemblies with higher quality and greater continuity.

FindingsHere, we report an improved assembly of the M. auratus genome (BCM_Maur_2.0) using Oxford Nanopore Technologies long-read sequencing to produce a chromosome-scale assembly. The total length of the new assembly is 2.46 Gbp, similar to the 2.50 Gbp length of a previous assembly of this genome, MesAur1.0. BCM_Maur_2.0 exhibits significantly improved continuity with a scaffold N50 that is 6.7 times greater than MesAur1.0. Furthermore, 21,616 protein coding genes and 10,459 noncoding genes are annotated in BCM_Maur_2.0 compared to 20,495 protein coding genes and 4,168 noncoding genes in MesAur1.0. This new assembly also improves the unresolved regions as measured by nucleotide ambiguities, where approximately 17.11% of bases in MesAur1.0 were unresolved compared to BCM_Maur_2.0 in which the number of unresolved bases is reduced to 3.00%.

ConclusionsAccess to a more complete reference genome with improved accuracy and continuity will facilitate more detailed, comprehensive, and meaningful research results for a wide variety of future studies using Syrian hamsters as models. The hexose transporter SWEET5 confers galactose sensitivity to Arabidopsis pollen germination via the galactokinase GALK Galactose is an abundant and essential sugar used for the biosynthesis of many macromolecules in different organisms, including plants. Galactose metabolism is tightly and finely controlled since excess galactose and derivatives are inhibitory. In Arabidopsis, root growth and pollen germination were strongly inhibited by excess galactose. However, the mechanism of galactose induced inhibition during pollen germination remains obscure. In this study, we characterized a plasma-membrane localized transporter, AtSWEET5, that transports glucose and galactose. SWEET5 protein level started to accumulate at the tricellular stage of pollen development and peaked in mature pollen before rapidly declining after pollen was germinated. SWEET5 levels are responsible for the dosage-dependent sensitivity of galactose and galactokinase (GALK) is essential for the inhibitory effects of galactose during pollen germination. Overall, SWEET5 and GALK contribute to the maintenance of galactose metabolic homeostasis during pollen germination. SWEET5 serves as a major low-affinity hexose transporter at the early stage of pollen germination.

One-sentence summarySWEET5 mediates pollen galactose sensitivity via GALK that is required for efficient galactose uptake in pollen during pollen germination. Plant host shift: The viability and persistence of a polyphagous, non-native insect Assessing the effects of a plant-host shift is important for monitoring insect populations over long time periods and for interventions in a conservation or pest management framework. In a heterogeneous environment, individuals may disperse between sources and sinks in order to persist. Here we propose a single-species two-patch model that aims to capture the generational movement of an insect that exhibits density-dependent dispersal, to see how shifting between hosts could alter its viability and asymptotic dynamics. We then analyse the stability and persistence properties of the model and further validate it using parameter estimates derived from laboratory experiments. In order to evaluate the potential of this model, we applied it to Drosophila suzukii (Diptera: Drosophilidae), which has become a harmful pest in several countries around the world. Although many studies have investigated the preference and attractiveness of potential hosts on this invasive drosophilid, no studies thus far have investigated whether a shift of fruit host could affect such a species ecological viability or spatiotemporal persistence. The model results show that a shift in host choice can significantly affect the growth potential and fecundity of a species such as D. suzukii, which ultimately could aid such invasive populations in their ability to persist within a changing environment. Clonal evolution in serially passaged Cryptococcus neoformans x deneoformans hybrids reveals a heterogenous landscape of genomic change Cryptococcus neoformans x deneoformans hybrids (also known as serotype AD hybrids) are basidiomycete yeasts that are common in a clinical setting. Like many hybrids, the AD hybrids are largely locked at the F1 stage and are mostly unable to undergo normal meiotic reproduction. However, these F1 hybrids, which display a high ([~]10%) sequence divergence are known to genetically diversify through mitotic recombination and aneuploidy, and this diversification may be adaptive. In this study, we evolved a single AD hybrid genotype in six diverse environments by serial passaging and then used genome resequencing of evolved clones to determine evolutionary mechanisms of adaptation. The evolved clones generally increased fitness after passaging, accompanied by an average of 3.3 point mutations, 2.9 loss of heterozygosity (LOH) events, and 0.7 trisomic chromosomes per clone. LOH occurred through nondisjunction of chromosomes, crossing over consistent with break-induced replication, and gene conversion, in that order of prevalence. The breakpoints of these recombination events were significantly associated with regions of the genome with lower sequence divergence between the parents and clustered in subtelomeric regions, notably in regions that had undergone introgression between the two parental species. Parallel evolution was observed, particularly through repeated homozygosity via nondisjunction, yet there was little evidence of environment-specific parallel change for either LOH, aneuploidy, or mutations. These data show that AD hybrids have both a remarkable genomic plasticity and yet are challenged in the ability to recombine through sequence divergence and chromosomal rearrangements, a scenario likely limiting the precision of adaptive evolution to novel environments. Multiple spillovers and onward transmission of SARS-Cov-2 in free-living and captive White-tailed deer Many animal species are susceptible to SARS-CoV-2 and could potentially act as reservoirs, yet transmission of the virus in non-human free-living animals has not been documented. White-tailed deer (Odocoileus virginianus), the predominant cervid in North America, are susceptible to SARS-CoV-2 infection, and experimentally infected fawns can transmit the virus. To test the hypothesis that SARS-CoV-2 may be circulating in deer, we tested 283 retropharyngeal lymph node (RPLN) samples collected from 151 free-living and 132 captive deer in Iowa from April 2020 through December of 2020 for the presence of SARS-CoV-2 RNA. Ninety-four of the 283 deer (33.2%; 95% CI: 28, 38.9) samples were positive for SARS-CoV-2 RNA as assessed by RT-PCR. Notably, between November 23, 2020 and January 10, 2021, 80 of 97 (82.5%; 95% CI 73.7, 88.8) RPLN samples had detectable SARS-CoV-2 RNA by RT-PCR. Whole genome sequencing of the 94 positive RPLN samples identified 12 SARS-CoV-2 lineages, with B.1.2 (n = 51; 54.5%), and B.1.311 (n = 19; 20%) accounting for ~75% of all samples. The geographic distribution and nesting of clusters of deer and human lineages strongly suggest multiple zooanthroponotic spillover events and deer-to-deer transmission. The discovery of sylvatic and enzootic SARS-CoV-2 transmission in deer has important implications for the ecology and long-term persistence, as well as the potential for spillover to other animals and spillback into humans. These findings highlight an urgent need for a robust and proactive "One Health" approach to obtaining a better understanding of the ecology and evolution of SARS-CoV-2.

One-Sentence SummarySARS-CoV-2 was detected in one-third of sampled white-tailed deer in Iowa between September 2020 and January of 2021 that likely resulted from multiple human-to-deer spillover and deer-to-deer transmission events. 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. Glucose dynamics during ozone exposure measured using radiotelemetry: Stress drivers BackgroundStress-related neurobehavioral and metabolic disorders are associated with altered circulating adrenal-derived hormones and hyperglycemia. Temporal assessment of glucose and these hormones is critical for insights on an individuals health.

ObjectivesHere we used implantable-telemetry in rats to assess real-time changes in circulating glucose during and after exposure to the air pollutant ozone, and link responses to circulating neuroendocrine stress and metabolic hormones. We also compared rodent glucose and corticosterone (cortisol in humans) responses to humans exposed to ozone.

MethodsFirst, using a cross-over design, we monitored glucose levels during single or repeated ozone exposures (0.0, 0.2, 0.4 and 0.8-ppm) and non-exposure periods in male Wistar-Kyoto-rats implanted with glucose-telemeters. A second cohort of un-implanted rats was exposed to ozone (0.0, 0.4 or 0.8-ppm) for 30-min, 1-hour, 2-hour, or 4-hour with hormones measured immediately after exposure. Then we assessed glucose metabolism in sham and adrenalectomized rats with or without pharmacological interventions of adrenergic and glucocorticoid receptors. Finally, we assessed glucose and cortisol in serum samples from a clinical study involving exposure of human volunteers to air or 0.3 ppm ozone.

ResultsOzone (0.8-ppm) caused hyperglycemia and hypothermia beginning 90-min into exposure, with reversal of effects 4-6 hours post-exposure. Glucose monitoring during four daily 4-hour ozone exposures revealed duration of hyperglycemia, adaptation, and diurnal variations. Ozone-induced hyperglycemia was preceded by increased adrenocorticotropic hormone, corticosterone, and epinephrine, but depletion of thyroid-stimulating, prolactin, and luteinizing hormones. Hyperglycemia was inhibited in rats that were adrenalectomized and/or treated with glucocorticoid inhibitor. The depletion of cortisol was dampened in humans exposed to ozone during intermittent exercise.

DiscussionWe demonstrate for the first time the temporality of neuroendocrine-stress-mediated biological sequalae responsible for ozone-induced metabolic dysfunction as exposure occurs. Real-time glucose monitoring with stress hormones assessment may be useful in identifying interactions among pollutants and stress-related illnesses. Testing for phylogenetic signal in single-cell RNA-seq data AO_SCPLOWBSTRACTC_SCPLOWPhylogenetic methods are emerging as a useful tool to understand cancer evolutionary dynamics, including tumor structure, heterogeneity, and progression. Most currently used approaches utilize either bulk whole genome sequencing (WGS) or single-cell DNA sequencing (scDNA-seq) and are based on calling copy number alterations and single nucleotide variants (SNVs). Here we explore the potential of single-cell RNA sequencing (scRNA-seq) to reconstruct cancer evolutionary dynamics. scRNA-seq is commonly applied to explore differential gene expression of cancer cells throughout tumor progression. The method exacerbates the single-cell sequencing problem of low yield per cell with uneven expression levels. This accounts for low and uneven sequencing coverage and makes SNV detection and phylogenetic analysis challenging. In this paper, we demonstrate for the first time that scRNA-seq data contains sufficient evolutionary signal and can be utilized in phylogenetic analyses. We explore and compare results of such analyses based on both expression levels and SNVs called from scRNA-seq data. Both techniques are shown to be useful for reconstructing phylogenetic relationships between cells, reflecting the clonal composition of a tumor. Both standardized expression values and SNVs appear to be equally capable of reconstructing a similar pattern of phylogenetic relationship. This pattern is stable even when phylogenetic uncertainty is taken in account. Our results open up a new direction of somatic phylogenetics based on scRNA-seq data. Further research is required to refine and improve these approaches to capture the full picture of somatic evolutionary dynamics in cancer. LPS-TLR4 pathway exaggerates alcohol-induced acute-on-chronic liver failure via provoking NETs formation Backgrounds: Intrahepatic infiltration of neutrophils is a character of alcoholic acute-on-chronic liver failure (AACLF) and neutrophil extracellular traps (NETs) are an important strategy for neutrophils to fix and kill invading microorganisms. Intestinal bacteria and the gut-liver axis have been thought to play a key role in many liver diseases also including AACLF. However, whether NETs appear in AACLF and play a role in AACLF is still unsure. Methods: WT, NE KO, and TLR4 KO mice were used to build the AACLF model, and the intestinal bacteria were eliminated at the same time and LPS was given. Then the formation of NETs and AACLF related markers were detected. Results: The serum MPO-DNA and LPS concentration was increased in AACLF patients and a correlation was revealed between these two indexes. More intrahepatic NETs formed in AACLF mice by testing MPO-DNA, Cit H3, and NE. These markers decreased with gut detergent and restored markers with gut detergent plus LPS supplement. While NETs formation failed to change with gut microbiome or combine LPS supplement in TLR4 KO mice. As we tested AACLF related characters, liver injury, intrahepatic fat deposition, inflammation, and fibrosis alleviated with depletion of NE. These related marks were also attenuated with gut sterilization by antibiotics and recovered with combined treatment with antibiotics plus LPS. But the liver injury, intrahepatic fat, fibro deposition, and liver inflammation-related markers did show a significant difference in TLR4 KO mice when they received the same treatment. Conclusion: Intestinal-derived LPS promotes NETs formation in AACLF through the TLR4 pathway and further accelerates the AACLF process by NETs. Heterogeneous interactions and polymer entropy decide organization and dynamics of chromatin domains Chromatin is known to be organized into multiple domains of varying sizes and compaction. While these domains are often imagined as static structures, they are highly dynamic and show cell-to-cell variability. Since processes such as gene regulation and DNA replication occur in the context of these domains, it is important to understand their organization, fluctuation and dynamics. To simulate chromatin domains, one requires knowledge of interaction strengths among chromatin segments. Here, we derive interaction strength parameters from experimentally known contact maps and use them to predict chromatin organization and dynamics. Taking two domains on the human chromosome as examples, we investigate its 3D organization, size/shape fluctuations, and dynamics of different segments within a domain, accounting for hydrodynamic effects. Considering different cell types, we quantify changes in interaction strengths and chromatin shape fluctuations in different epigenetic states. Perturbing the interaction strengths systematically, we further investigate how epigenetic-like changes can alter the spatio-temporal nature of the domains. Our results show that heterogeneous weak interactions are crucial in determining the organization of the domains. Computing effective stiffness and relaxation times, we investigate how perturbations in interactions affect the solid-like and liquid-like nature of chromatin domains. Quantifying dynamics of chromatin segments within a domain, we show how the competition between polymer entropy and interaction energy influence the timescales of loop formation and maintenance of stable loops. A nontoxic tetrodotoxin analog, 5,6,11-trideoxytetrodotoxin, as an olfactory chemoattractant for the grass puffer Toxic puffers contain the potent neurotoxin, tetrodotoxin (TTX). Although TTX is considered to serve as a defense substance, previous behavioral studies have demonstrated that TTX acts as an attractive pheromone for some toxic puffers. To elucidate the physiological mechanism of putative pheromonal action of TTX, we examined whether grass puffers Takifugu alboplumbeus can detect TTX. Electroolfactogram (EOG) results suggest that the olfactory epithelium (OE) of grass puffers responded to a type of TTX analog (5,6,11-trideoxyTTX), although it did not respond to TTX. We also examined the attractive action of 5,6,11-trideoxyTTX on grass puffers by recording their swimming behavior under dark conditions. Grass puffers preferred to stay on the side of the aquarium where 5,6,11-trideoxyTTX was administered, and their swimming speed decreased. Additionally, odorant-induced labeling of olfactory sensory neurons by immunohistochemistry against neural activity marker (phosphorylated extracellular signal regulated kinase; pERK) revealed that labeled olfactory sensory neurons were localized in the region surrounding "islets" where there was considered as nonsensory epithelium. 5,6,11-trideoxyTTX has been known to accumulate in grass puffers, but its toxicity is much lower (almost nontoxic) than TTX. Our results suggest that toxic puffers may positively use this unknown function TTX analog as an olfactory chemoattractant. Active antithrombin glycoforms are selectively physiosorbed on plasma extracellular vesicles Antithrombin (AT) is a glycoprotein produced by the liver and a principal antagonist of active clotting proteases. A deficit in AT function leads to AT qualitative deficiency, challenging to diagnose. Here we report that active AT may travel physiosorbed on the surface of plasma extracellular vesicles (EVs), contributing to form the "EV-protein corona". The corona is enriched in specific AT glycoforms, thus suggesting glycosylation to play a key role in AT partitioning between EVs and plasma. Differences in AT glycoform composition of the corona of EVs separated from plasma of healthy and AT qualitative deficiency-affected subjects were also noticed. This suggests deconstructing the plasma into its nanostructured components, as extracellular vesicles, could help to unravel pathophysiological mechanisms otherwise undiscovered. Once-daily feeding is associated with better health in companion dogs: Results from the Dog Aging Project A variety of diets have been studied for possible anti-aging effects. In particular, studies of intermittent fasting and time-restricted feeding in laboratory rodents have found evidence of beneficial health outcomes. Companion dogs represent a unique opportunity to study diet in a large mammal that shares human environments. The Dog Aging Project has been collecting data on thousands of companion dogs of all different ages, sizes, and breeds since 2019. We leveraged this diverse cross-sectional dataset to investigate associations between feeding frequency and cognitive function (n = 10,474) as well as nine broad categories of health conditions (n = 24,238). Controlling for sex, age, breed, and other potential confounders, we found that dogs fed once daily rather than more frequently had lower mean scores on a cognitive dysfunction scale, and lower odds of having gastrointestinal, dental, orthopedic, kidney/urinary, and liver/pancreas disorders. Therefore, we find that once-daily feeding is associated with better health in multiple domains. Future research with longitudinal data can provide stronger evidence for a possible causal effect of feeding frequency on health in companion dogs. Genome-wide analysis of the apple family 1 glycosyltransferases identified a flavonoid-modifying UGT, MdUGT83L3, which is targeted by MdMYB88 and contributes to stress adaptation Nowadays, the plant family 1 UDP-glycosyltransferases (UGTs) are more and more investigated owing to their contribution to the plant secondary metabolism and diverse biological roles. Apple (Malus domestica) is one of the most widely cultivated fruit trees with great economical importance. However, poor knowledge is known about the apple UGTs. In this study, we identified 229 members of family 1 through a genome-wide analysis of the apple UGTs, which were clustered into 18 groups, from A to R. We also performed detailed analysis to 34 apple UGTs with quantitavive RT-PCR, and discovered a number of stress-regulated UGTs. Among them, we characterized the role of MD09G1064900, also named as MdUGT83L3, which is significantly induced by salt and cold. In vivo analysis showed that it has high activity towards cyanidin, and moderate activity towards quercetin and keampferol. The transgenic callus and regenerated apple plants overexpressing MdUGT83L3 showed enhanced tolerance to salt and cold treatments. Overexpression of MdUGT83L3 also increased anthocyanin accumulation in the callus tissues and enhanced ROS clearing upon exposed to salt and cold stresses. Futhermore, via yeast-two-hybrid assay, EMSA and CHIP analysis, we also found that MdUGT83L3 could be directly regualted by MdMYB88. Our study indicated that MdUGT83L3, under the regulation of MdMYB88, plays important roles in salt and cold stress adaptation via modulating flavonoid metabolism in apple. A tool for assessing alignment of biomedical data repositories with open, FAIR, citation and trustworthy principles Increasing attention is being paid to the operation of biomedical data repositories in light of efforts to improve how scientific data is handled and made available for the long term. Multiple groups have produced recommendations for functions that biomedical repositories should support, with many using requirements of the FAIR data principles as guidelines. However, FAIR is but one set of principles that has arisen out of the open science community. They are joined by principles governing open science, data citation and trustworthiness, all of which are important aspects for biomedical data repositories to support. Together, these define a framework for data repositories that we call OFCT: Open, FAIR, Citable and Trustworthy. Here we developed an instrument using the open source PolicyModels toolkit that attempts to operationalize key aspects of OFCT principles and piloted the instrument by evaluating eight biomedical community repositories listed by the NIDDK Information Network (dkNET.org). Repositories included both specialist repositories that focused on a particular data type or domain, in this case diabetes and metabolomics, and generalist repositories that accept all data types and domains. The goal of this work was both to obtain a sense of how much the design of current biomedical data repositories align with these principles and to augment the dkNET listing with additional information that may be important to investigators trying to choose a repository, e.g., does the repository fully support data citation? The evaluation was performed from March to November 2020 through inspection of documentation and interaction with the sites by the authors. Overall, although there was little explicit acknowledgement of any of the OFCT principles in our sample, the majority of repositories provided at least some support for their tenets. Implicit Laplacian of Enhanced Edge: An Unguided Algorithm for Accurate and Automated Quantitative Analysis of Cytoskeletal Images The eukaryotic cytoskeleton plays essential roles in cell signaling and trafficking, which is broadly associated with immunity and diseases of human and plants. To date, most analyses aiming at defining the temporal and spatial dynamics of the cytoskeleton have relied on qualitative and quantitative analysis of fluorescence images to describe cytoskeletal function. While state-of-the-art, these approaches have limitations: the diverse shape and brightness of the cytoskeleton cause considerable inaccuracy in both human-driven and automated approaches, and the widely adopted image projection process (3D to 2D) leads to bias and information loss. Here, we describe the development and application of Implicit Laplacian of Enhanced Edge (ILEE), an unguided approach that uses a 2D/3D-compatible local thresholding algorithm for the quantitative evaluation of cytoskeletal status and organization at high performance. Using ILEE, we constructed a Python library to enable automated cytoskeletal image analysis, providing numerous biologically-interpretable indices measuring the density, bundling, severing, branching, and directionality of the cytoskeleton. The data presented herein demonstrate that ILEE resolves the defects of classic cytoskeleton analysis approaches, enables the measurement of novel cytoskeletal features, and yields quantitatively descriptive data with superior accuracy, stability, and robustness. We released the ILEE algorithm as an open-source library and further developed a Google Colab interface as a community resource.