State-dependent network interactions differentially gate sensory input at the motor and command neuron level in Caenorhabditis elegans Neural responses are influenced by both external stimuli and internal network states. While network states have been linked to behavioral and stimulus states, little is known about how sensory inputs are filtered by whole-brain activity to affect motor and command neurons. Here, we recorded whole-brain activity of Caenorhabditis elegans experiencing bacterial food stimuli, and modeled how sensory inputs affect motor and command neurons in a network state-dependent manner. First, we classified active neurons into six functional clusters: two sensory neuron clusters (ON, OFF), and four motor/command neuron clusters (AVA, RME, SMDD, SMDV). Using encoding models, we found that ON and OFF sensory neurons that respond to onset and removal of bacteria, respectively, employ different adaptation strategies. Next, we used decoding models to show that bacterial onset and removal differentially drive AVA and RME cluster activity. To explore state-dependent effects on AVA and RME clusters, we developed a model that identified network states and fitted submodels for each state to predict how each of the six functional clusters drive AVA and RME cluster activity. We also identified network states in which AVA and RME clusters were either largely unperturbed by or receptive to bacterial sensory input. Furthermore, this model allowed us to disentangle the state-dependent contributions of stimulus timescales and bacterial content to neural activity. Collectively, we present an interpretable approach for modeling network dynamics that goes beyond implication of neurons in particular states, and moves toward explicitly dissecting how neural populations work together to produce state dependence.

Significance StatementA major function of the brain is to transform sensory information into behavior. As the first receiver of sensory input, sensory neuron activity is often most correlated with stimulus features. However, this high-fidelity representation of sensory input becomes diluted as it travels to downstream neurons, where sensory information is integrated with network activity. By the time sensory information reaches motor neurons, it is often difficult to dissociate the influence of sensory input from the influence of network activity. Here, we describe a method that is fully interpretable such that we can show how neural populations on a whole-brain scale interact to produce network states. From there, we can attribute motor neuron activity to network history and sensory input.
Human cortical areas are sparsely connected: Combining histology with diffusion MRI to estimate the absolute number of axons The tracts between cortical areas are conceived as playing a central role in cortical information processing, but their actual numbers have never been determined in humans. Here we estimate the absolute number of axons linking cortical areas from a whole-cortex diffusion-MRI (dMRI) connectome, calibrated using the histologically-measured callosal fiber density. Median connectivity is estimated as ~6200 axons between cortical areas within-hemisphere and ~1300 axons inter-hemispherically, with axons connecting functionally-related areas surprisingly sparse. For example, we estimate that <5% of the axons in the trunk of the arcuate and superior longitudinal fasciculi connect Wernickes and Brocas areas. These results suggest that detailed information is transmitted between cortical areas either via linkage of the dense local connections or via rare, extraordinarily privileged long-range connections.
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. Mate sharing in male stump-tailed macaques as a possible case of coalition-like behavior to modify the group-wise fitness distribution In multi-male multi-female groups of animals, male reproductive success is often skewed toward top-ranking males. Coalition formation by middle- to low-rankers can be seen as a collaborative effort to modify the distribution of reproductive success within the group, so that they gain more than they would do otherwise. It has been pointed out, on the other hand, that a coalition of top-ranking males could hardly be profitable in the sense that they would gain little additional benefit from making collaborative effort. Here we report our novel observation of facultative sharing of mating opportunities among males in a wild population of stump-tailed macaques (Macaca arctoides) as a possible case of coalition-like behavior in which dominant males jointly guard females from mating with subordinate males and actively share mating opportunities within the allies. First, we report our novel observation of facultative sharing of mating opportunities in male stump-tailed macaques, where two or more males remain in close proximity to and copulate with a female in turn without contesting or sneaking. Second, considering the kind of coalition formation in which dominant males collaboratively exclude subordinates from mating competition and thereby strengthen the reproductive skew that already exists, we specify, by means of mathematical modeling, the condition for this kind of coalition formation to be selectively favored. Finally, we derive predictions about the occurrence of the coalition-like behavior depending on ecological and demographic factors, and test them empirically using data from the five groups of stump-tailed macaques in our study population. 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. Glucose inhibits haemostasis and accelerates diet-induced hyperlipidaemia in zebrafish larvae Hyperglycaemia damages the microvasculature in part through the reduced recruitment of immune cells and interference with platelet signalling, leading to poor wound healing and accelerated lipid deposition in mammals. We investigated the utility of zebrafish larvae to model the effect of glucose on neutrophil and macrophage recruitment to a tail wound, wound-induced haemostasis, and chicken egg yolk feed challenge-induced hyperlipidaemia by supplementing larvae with exogenous glucose by immersion or injection. Neither method of glucose supplementation affected the recruitment of neutrophils and macrophages following tail transection. Glucose injection reduced thrombocyte retention and fibrin plug formation while only thrombocyte retention was reduced by glucose immersion following tail transection. We observed accelerated lipid accumulation in glucose-injected larvae challenged with high fat chicken egg yolk feeding. Our study identifies conserved and divergent effects of high glucose on inflammation, haemostasis, and hyperlipidaemia in zebrafish larvae compared to mammals. Trehalose increases tomato drought tolerance, induces defenses, and increases resistance to bacterial wilt disease Ralstonia solanacearum causes bacterial wilt disease, leading to severe crop losses. Xylem sap from R. solanacearum-infected tomato is enriched in the disaccharide trehalose. Water-stressed plants also accumulate trehalose, which increases drought tolerance via abscisic acid (ABA) signaling. Because R. solanacearum-infected plants suffer reduced water flow, we hypothesized that bacterial wilt physiologically mimics drought stress, which trehalose could mitigate. We found that R. solanacearum-infected plants differentially expressed drought-associated genes, including those involved in ABA and trehalose metabolism, and had more ABA in xylem sap. Consistent with this, treating tomato roots with ABA reduced both stomatal conductance and stem colonization by R. solanacearum. Treating roots with trehalose increased xylem sap ABA and reduced plant water use by lowering stomatal conductance and temporarily improving water use efficiency. Trehalose treatment also upregulated expression of salicylic acid (SA)-dependent tomato defense genes; increased xylem sap levels of SA and other antimicrobial compounds; and increased bacterial wilt resistance of SA-insensitive NahG tomato plants. Additionally, trehalose treatment increased xylem concentrations of jasmonic acid and related oxylipins. Finally, trehalose-treated plants were substantially more resistant to bacterial wilt disease. Together, these data show that exogenous trehalose reduced both water stress and bacterial wilt disease and triggered systemic disease resistance, possibly through a Damage Associated Molecular Pattern (DAMP) response pathway. This suite of responses revealed unexpected linkages between plant responses to biotic and abiotic stress and suggested that R. solanacearum-infected plants increase trehalose to improve water use efficiency and increase wilt disease resistance. The pathogen may degrade trehalose to counter these efforts. Together, these results suggest that treating tomatoes with exogenous trehalose could be a practical strategy for bacterial wilt management. Genomic databanks and targeted assays help characterise domestic mosquito incursions Biosecurity strategies aiming to restrict the spread of invasive pests can benefit from knowing where new incursions have come from. This knowledge can be acquired using genomic databanks, where genetic variation in incursion samples is compared against variation in reference samples. Here, we use genomic databanks to investigate domestic incursions of two mosquito species in Australia. We used a deep learning method to trace a 2021 invasion of Aedes aegypti in Tennant Creek, Northern Territory, to Townsville, Queensland, and to trace Ae. albopictus incursions between several invaded and uninvaded islands in the Torres Strait. We observed high precision of tracing despite 30-70 generations separating incursions and source samples. Targeted assays also improved precision, in this case by comparing Wolbachia infection data to mitochondrial DNA variation in Ae. aegypti. For Ae. albopictus, incursions detected in 2021 were traced to different islands than incursions detected in 2019. Patterns of relatedness and inbreeding indicated that Tennant Creek was likely invaded by Ae. aegypti from one family, whereas Torres Strait incursions involved distinct kinship groups. Our results highlight the value of repeated tracing with genomic databanks that remain informative over years, and demonstrate how additional targeted assays (e.g. Wolbachia) can improve inferences. Herbivore life histories are altered by drought stress in their hosts plants Drought associated with climate change can stress plants, altering their interactions with phytophagous arthropods. Dryness not only impacts cultivated plants but also their parasites, which in some cases are favoured by drought. Herbivorous arthropods feeding on drought-stressed plants typically produce bigger offspring and develop faster. However, how much responses to abiotic factors differ among populations of a species remains poorly documented. Here, we evaluate variability among populations of a major agricultural pest, the two spotted spider mite, Tetranychus urticae, in response to drought-stress. We compare key life-history parameters of twelve populations that originate from climates ranging from wet and cool Atlantic locations to medium to dry hot Mediterranean locations. We evaluated how plant drought stress affects four life history traits: development time, fecundity, sex-ratio and emigration rate in an experiment comparing well-watered and drought-stressed bean plants. Mites feeding on drought-stressed plants developed faster and attempted to leave leaves less often, and young females were more fecund. The mites from wet temperate climates exhibited greater plasticity between the two water regimes than mites originating from dryer and hot climates. Thus, climate in the area of origin influences mite response to drought. Basal Gp78-dependent mitophagy promotes mitochondrial health and limits mitochondrial ROS production Mitochondria are major sources of cytotoxic reactive oxygen species (ROS) that contribute to cancer progression. Mitophagy, the selective elimination of mitochondria by autophagy, monitors and maintains mitochondrial health and integrity, eliminating ROS-producing mitochondria. However, mechanisms underlying mitophagic control of mitochondrial homeostasis under basal conditions remain poorly understood. Gp78 E3 ubiquitin ligase is an endoplasmic reticulum membrane protein that induces mitochondrial fission and mitophagy of depolarized mitochondria. Here, we report that CRISPR/Cas9 knockout of Gp78 in HT-1080 fibrosarcoma cells increased mitochondrial volume and rendered cells resistant to carbonyl cyanide m-chlorophenyl hydrazone (CCCP)-induced mitophagy. These effects were phenocopied by knockdown of the essential autophagy protein ATG5 in wild-type HT-1080 cells. Use of the mito-Keima mitophagy probe confirmed that Gp78 promoted both basal and damage-induced mitophagy. Application of a spot detection algorithm (SPECHT) to GFP-mRFP tandem fluorescent-tagged LC3 (tfLC3)-positive autophagosomes reported elevated autophagosomal maturation in wild-type HT-1080 cells relative to Gp78 knockout cells, predominantly in proximity to mitochondria. Mitophagy inhibition by either Gp78 knockout or ATG5 knockdown reduced mitochondrial potential and increased mitochondrial ROS. Live cell analysis of tfLC3 in HT-1080 cells showed the preferential association of autophagosomes with mitochondria of reduced potential. Basal Gp78-dependent mitophagic flux is therefore selectively associated with reduced potential mitochondria promoting maintenance of a healthy mitochondrial population and limiting ROS production. 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. Spatial proteomics reveals subcellular reorganization in human keratinocytes exposed to UVA light The effects of UV light on the skin have been extensively investigated. However, systematic information about how exposure to UVA light, the least energetic but the most abundant UV radiation reaching the Earth, shapes the subcellular organization of proteins is lacking. Using subcellular fractionation, mass spectrometry-based proteomics, machine learning algorithms, immunofluorescence, and functional assays, we mapped the subcellular reorganization of the proteome of human keratinocytes in response to UVA light. Our workflow quantified and assigned subcellular localization for over 1600 proteins, of which about 200 were found to redistribute upon UVA exposure. Reorganization of the proteome affected modulators of signaling pathways, cellular metabolism, and DNA damage response. Strikingly, mitochondria were identified as one of the main targets of UVA-induced stress. Further investigation demonstrated that UVA induces mitochondrial fragmentation, up-regulates redox-responsive proteins and attenuates respiratory rates. These observations emphasize the role of this radiation as a potent metabolic stressor in the skin. The Consequences of Differential Origin Licensing Dynamics in Distinct Chromatin Environments MCM complexes are loaded onto chromosomes to license DNA replication origins in G1 phase of the cell cycle, but it is not yet known how mammalian MCM complexes are adequately distributed to both euchromatin and heterochromatin. To address this question, we combined time-lapse live-cell imaging with fixed cell immunofluorescence imaging of single human cells to quantify the relative rates of MCM loading in heterochromatin and euchromatin at different times within G1. We report here that MCM loading in euchromatin is faster than in heterochromatin in very early G1, but surprisingly, heterochromatin loading accelerates relative to euchromatin loading in middle and late G1. These different loading dynamics require ORCA-dependent differences in ORC distribution during G1. A consequence of heterochromatin origin licensing dynamics is that cells experiencing a truncated G1 phase from premature cyclin E expression enter S phase with under-licensed heterochromatin, and DNA damage accumulates preferentially in heterochromatin in the subsequent S/G2 phase. Thus G1 length is critical for sufficient MCM loading, particularly in heterochromatin, to ensure complete genome duplication and to maintain genome stability. Micro topographical instruction of bacterial attachment, biofilm formation and in vivo host response Bio-instructive materials that prevent bacterial biofilm formation and drive an appropriate host immune response have the potential to significantly reduce the burden of medical device-associated infections. As bacterial surface attachment is sensitive to surface topography, we experimentally surveyed 2,176 combinatorially generated shapes using an unbiased high-throughput micro topographical polystyrene polymer chip screen. This identified topographies that reduced colonization in vitro by up to 15-fold compared with a flat surface for both motile and non-motile bacterial pathogens. Equivalent reductions were achieved on polyurethane, a polymer commonly used in medical devices. Using machine learning methods, a set of design rules based on generalisable topo-descriptors was established for predicting bacteria-resistant micro topographies. In a murine foreign body infection model, anti-attachment topographies were shown to be refractory to P. aeruginosa colonization and serendipitously, the fibrotic response to the implant was reduced, offering control of foreign body response using simple topographical patterning of non-eluting implants. 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. Simultaneous suppression of ribosome biogenesis and Tor activation by TRIM-NHL proteins promotes terminal differentiation Proper tissue development and homeostasis depends on the balance between growth and terminal differentiation, but the mechanisms coordinating these processes remain elusive. Accumulating evidence indicates that ribosome biogenesis (RiBi) and protein synthesis, two of the most energy-consuming cellular processes sustaining growth, are tightly regulated and yet can be uncoupled during stem cell fate transitions. Here, using the Drosophila adult female germline stem cell (GSC) and larval neuroblast (NB) systems, we show that Mei-P26 and Brat, two Drosophila TRIM-NHL paralogs of the mammalian TRIM32 protein family, are responsible for uncoupling RiBi and protein synthesis during GSC and NB differentiation, respectively. Mei-P26 and Brat modify the metabolism of differentiating cells by activating the Target of rapamycin (Tor) kinase to promote translation, while concomitantly repressing RiBi. Depletion of Mei-P26 or Brat results in excessive cellular growth and defective terminal differentiation, which can be rescued by ectopic activation of Tor together with suppression of RiBi. Our results indicate that the anabolic reprogramming established by TRIM-NHL activity by uncoupling RiBi and translation activities creates the conditions required for terminal differentiation. Virus-host interactions and genetic diversity of Antarctic sea ice bacteriophages Despite generally appreciated significant roles of microbes in sea ice and polar waters, detailed studies of virus-host systems from such environments have been so far limited by only a few isolates. Here, we investigated infectivity under various conditions, infection cycles, and genetic diversity of Antarctic sea ice bacteriophages: PANV1, PANV2, OANV1, and OANV2. The phages infect common sea ice bacteria belonging to the genera Paraglaciecola or Octadecabacter. Although the phages are marine and cold-active, replicating at 0-5 {degrees}C, they all survived temporal incubations at [&ge;]30 {degrees}C and remained infectious without any salts or supplemented only with magnesium, suggesting a robust virion assembly maintaining integrity under a wide range of conditions. Host recognition in the cold proved to be effective, and the release of progeny viruses occurred as a result of cell lysis. The analysis of viral genome sequences showed that nearly half of the gene products of each virus are unique, highlighting that sea ice harbors unexplored virus diversity. Based on predicted genes typical for tailed double-stranded DNA phages, we suggest placing the four studied viruses in the class Caudoviricetes. Searching against viral sequences from metagenomic assemblies revealed that related viruses are not restricted to Antarctica, but also found in distant marine environments.

ImportanceVery little is known about sea ice microbes despite the significant role of sea ice in the global oceans as well as microbial input into biogeochemical cycling. Studies on the sea ice viruses have been typically limited to -omics-based approaches and microscopic examinations of sea ice samples. Up to date, only four cultivable viruses have been isolated from Antarctic sea ice. Our study of these unique isolates advances the understanding of the genetic diversity of viruses in sea ice environments, their interactions with host microbes, and possible links to other biomes. Such information contributes to more accurate future sea ice biogeochemical models. The addition of FAIMS Increases Targeted Proteomics Sensitivity from FFPE Tumor Biopsies Mass spectrometry-based targeted proteomics allows objective protein quantitation of clinical biomarkers from a single section of formalin-fixed, paraffin-embedded (FFPE) tumor tissue biopsies. We combined high-field asymmetric waveform ion mobility spectrometry (FAIMS) and parallel reaction monitoring (PRM) to increase assay sensitivity. The modular nature of the FAIMS source allowed direct comparison of the performance of FAIMS-PRM to PRM. Limits of quantitation were determined by spiking synthetic peptides into a human spleen matrix. In addition, 20 clinical samples were analyzed using FAIMS-PRM and the quantitation of HER2 was compared with that obtained with the Ventana immunohistochemistry assay. FAIMS-PRM improved the overall signal-to-noise ratio over that from PRM and increased assay sensitivity in FFPE tissue analysis for four (HER2, EGFR, cMET, and KRAS) of five proteins of clinical interest. FAIMS-PRM enabled sensitive quantitation of basal HER2 expression in breast cancer samples classified as HER2 negative by immunohistochemistry. Furthermore, we determined the degree of FAIMS-dependent background reduction and showed that this correlated with an improved lower limit of quantitation with FAIMS. FAIMS-PRM is anticipated to benefit clinical trials in which multiple biomarker questions must be addressed and the availability of tumor biopsy samples is limited. RIG-I immunotherapy overcomes radioresistance in p53-positive malignant melanoma Radiation therapy induces cytotoxic DNA damage, which results in cell-cycle arrest and activation of cell-intrinsic death pathways, but its application has been limited by the radioresistance of tumors, such as in malignant melanoma. RIG-I is a cytosolic immune receptor expressed in all somatic cells, including tumor cells, with a key role in sensing viral RNA. RIG-I specific oligonucleotide ligands elicit a robust cell-intrinsic antiviral response and immunogenic cell death in tumor cells and are being tested in clinical trials. Nonetheless, their potential to overcome radioresistance has not yet been explored. Here, we demonstrate that activation of RIG-I enhances the extent and immunogenicity of irradiation-induced tumor cell death in human and murine melanoma cell lines in vitro and improved survival in the murine B16 melanoma model. Pathway analysis of transcriptomic data revealed a central role for p53 downstream of the combination treatment, which was corroborated using p53-/- B16 cells. In vivo, the effect of irradiation on immune-cell activation and inhibition of tumor growth was absent in mice carrying p53-/- B16 tumors, while the response to RIG-I stimulation in those mice was maintained. Our results identify p53 as pivotal for the synergistic antitumoral effect of RIG-I and irradiation, resulting in potent induction of immunogenic tumor-cell death. Thus, the administration of RIG-I ligands in combination with radiotherapy is a promising therapeutic approach to treating radioresistant tumors with a functional p53 pathway, such as malignant melanoma. The ecological consequences and evolution of retron-mediated suicide as a way to protect Escherichia coli from being killed by phage In 1984, retrons were described as DNA sequences that code for a reverse transcriptase and a unique single-stranded DNA/RNA hybrid called multicopy single-stranded DNA (msDNA). It would not be until 2020 that a function was shown, when compelling evidence was presented that retrons can protect bacteria from infections by lytic phage via an abortive infection mechanism. When infected with the virulent mutant of the phage lambda, {lambda}VIR, and to a lesser extent other phages, a retron designated Ec48 is activated, the Escherichia coli bearing this element dies and the infecting phage is lost. With the aid of a mathematical model, we explored the a priori conditions under which retrons will protect bacterial populations from invasion by phage and when retron-bearing bacteria will evolve in populations without this element. Using E. coli with and without Ec48, and {lambda}VIR we estimated the parameters of this model and tested the hypotheses generated from our analysis of its properties. Our models and experiments suggest cells expressing retron-mediated abortive infection have no advantage over cells lacking this system as a defense against phage, and we question whether the protection against phage infection is the unique selective force responsible for their evolution in bacteria. Synthetic Cooling Agents in US-marketed E-cigarette Refill Liquids and Popular Disposable E-cigarettes: Chemical Analysis and Risk Assessment. BackgroundMenthol, through its cooling and pleasant sensory effects, facilitates smoking and tobacco product initiation, resulting in the high popularity of mint/menthol-flavored E-cigarettes. More recently, E-cigarette vendors started marketing synthetic cooling agents as additives that impart a cooling effect but lack a characteristic minty odor. Knowledge about content of synthetic coolants in US-marketed E-cigarette products and associated health risks is limited.

MethodsE-liquid vendor sites were searched with the terms "koolada", "kool/cool", "ice", or WS-3/WS-23, denoting individual cooling agents, and relevant refill E-liquids were purchased. "Ice" flavor varieties of Puffbar, the most popular disposable E-cigarette brand, were compared with non-"Ice" varieties. E-liquids were characterized, and synthetic coolants quantified using GC/MS. Margin of exposure (MOE), a risk assessment parameter, was calculated to assess the risk associated with synthetic coolant exposure from E-cigarette use.

ResultsWS-3 was detected in 24/25 refill E-liquids analyzed. All Puffbar flavor varieties contained either WS-23 (13/14) or WS-3 (5/14), in both "Ice"- and non-"Ice" flavors. Modeling consumption of WS-3 from vaped E-liquids, resulted in MOEs below the safe margin of 100 for most daily use scenarios. MOEs for WS-23 were <100 for 10/13 Puffbar flavors in all use scenarios. Puffbar power specifications are identical to Juul devices.

ConclusionsSynthetic cooling agents (WS-3/WS-23) were present in US-marketed E-cigarettes, at levels that may result in consumer exposures exceeding safety thresholds set by regulatory agencies. Synthetic coolants are not only found in mint-or menthol-flavored products, but also in fruit- and candy-flavored products, including popular disposable E-cigarette products such as Puffbar.

ImplicationsSynthetic cooling agents are widely used in "kool/cool"- and "ice"-flavored E-liquids and in E-liquids without these labels, both as a potential replacement for menthol or to add cooling  notes to non-menthol flavors. These agents may be used to bypass current and future regulatory limits on menthol content in tobacco products, and not just E-cigarettes. Since synthetic cooling agents are odorless, they may not fall under the category of "characterizing flavor", potentially circumventing regulatory measures based on this concept. Regulators need to consider the additional health risks associated with exposure to synthetic cooling agents. 'rbioacc': an R-package to analyse toxicokinetic data O_LI rbioacc is an  R-package dedicated to the analysis of experimental data collected from bioaccumulation tests during which organisms are exposed to a chemical (exposure phase) and then put into a clean media (depuration phase). Internal concentrations are regularly measured over time all along the experiment.
C_LIO_LI rbioacc provides ready-to-use functions to visualize and fully analyze such data. Under a Bayesian framework, this package fits a generic one-compartment toxicokinetic (TK) model automatically built from the data. It provides TK parameter estimates (appropriate uptake and elimination rates) and bioaccumulation metrics (e.g., BCF, BSAF, BMF). All parameter estimates, bioaccumulation metrics as well as predictions of internal concentrations into organisms are delivered with their uncertainty.
C_LIO_LIBioaccumulation metrics are primarily provided in support of environmental risk assessment, in full compliance with regulatory requirements required to approve marketing applications of chemical substances.
C_LIO_LIThis paper gives brief worked examples of the use of  rbioacc from data collected through standard bioaccumulation tests, and publicly available within the scientific literature. These examples constitute step-by-step user-guides to analyze any new data set, uploaded in the right format.
C_LI Structures of the Omicron spike trimer with ACE2 and an anti-Omicron antibody: mechanisms for the high infectivity, immune evasion and antibody drug discovery The Omicron variant of SARS-CoV-2 has rapidly become the dominant infective strain and the focus efforts against the ongoing COVID-19 pandemic. Here we report an extensive set of structures of the Omicron spike trimer by its own or in complex with ACE2 and an anti-Omicron antibody. These structures reveal that most Omicron mutations are located on the surface of the spike protein, which confer stronger ACE2 binding by nearly 10 folds but become inactive epitopes resistant to many therapeutic antibodies. Importantly, both RBD and the closed conformation of the Omicron spike trimer are thermodynamically unstable, with the melting temperature of the Omicron RBD decreased by as much as 7{degrees}C, making the spiker trimer prone to random open conformations. An unusual RBD-RBD interaction in the ACE2-spike complex unique to Omicron is observed to support the open conformation and ACE2 binding, serving the basis for the higher infectivity of Omicron. A broad-spectrum therapeutic antibody JMB2002, which has completed Phase 1 clinical trial, is found to interact with the same two RBDs to inhibit ACE2 binding, in a mode that is distinguished from all previous antibodies, thus providing the structural basis for the potent inhibition of Omicron by this antibody. Together with biochemical data, our structures provide crucial insights into higher infectivity, antibody evasion and inhibition of Omicron. G protein-coupled estrogen receptor activation by Bisphenol-A disrupts protection from apoptosis conferred by estrogen receptors ERα and ERβ in pancreatic beta cells 17{beta}-estradiol protects pancreatic {beta}-cells from apoptosis via the estrogen receptors ER, ER{beta} and GPER. Conversely, the endocrine disruptor Bisphenol-A (BPA), which exerts multiple effects in this cell type via the same estrogen receptors, increased basal apoptosis. The molecular initiated events that trigger these opposite actions have yet to be identified. We demonstrated that combined genetic downregulation and pharmacological blockade of each estrogen receptor increased apoptosis to a different extent. The increase in apoptosis induced by BPA was diminished by the pharmacological blockade or the genetic silencing of GPER, and it was partially reproduced by the GPER agonist G1. BPA and G1-induced apoptosis were abolished upon pharmacological inhibition, silencing of ER and ER{beta}, or in dispersed islet cells from ER{beta} knockout (BERKO) mice. Yet, the ER and ER{beta} agonists, PPT and DPN, respectively, had no effect on beta cell viability. To exert their biological actions, ER and ER{beta} form homodimers and heterodimers. Molecular dynamic simulations together with proximity ligand assay and coimmunoprecipitation experiments indicated that the interaction of BPA with ER and ER{beta} as well as the GPER activation by G1 decreased ER{beta} heterodimers. We propose that ER{beta} heterodimers play an antiapoptotic role in beta cells and that BPA- and G1-induced decrease in ER{beta} heterodimers leads to beta cell apoptosis. Unveiling how different estrogenic chemicals affect the crosstalk among estrogen receptors should help to identify diabetogenic endocrine disruptors.

HighlightsO_LIPharmacological blockade and gene silencing of estrogen receptors ER, ER{beta} and GPER indicate that they are antiapoptotic in basal conditions.
C_LIO_LIGPER activation by G1 and BPA triggered apoptosis via a crosstalk with ER and ER{beta}.
C_LIO_LIBPA interaction with ER and ER{beta} as well as GPER activation decreased ER{beta} heterodimers, which was associated to increased apoptosis.
C_LIO_LIThis pathway represents a novel molecular initiating event underlying the pro-apoptotic effect of BPA
C_LIO_LIThe EndoC-{beta}H1 cell line may be a valid model of human {beta}-cells for identifying diabetogenic pollutants.
C_LI Widespread redundancy in -omics profiles of cancer mutation states In studies of cellular function in cancer, researchers are increasingly able to choose from many -omics assays as functional readouts. Choosing the correct readout for a given study can be difficult, and which layer of cellular function is most suitable to capture the relevant signal may be unclear. In this study, we consider prediction of cancer mutation status (presence or absence) from functional -omics data as a representative problem. Since functional signatures of cancer mutation have been identified across many data types, this problem presents an opportunity to quantify and compare the ability of different -omics readouts to capture signals of dysregulation in cancer. The TCGA Pan-Cancer Atlas contains genetic alteration data including somatic mutations and copy number variants (CNVs), as well as several -omics data types. From TCGA, we focus on RNA sequencing, DNA methylation arrays, reverse phase protein arrays (RPPA), microRNA, and somatic mutational signatures as -omics readouts.

Across a collection of genes recurrently mutated in cancer, RNA sequencing tends to be the most effective predictor of mutation state. Surprisingly, we found that for many of the genes we considered, one or more other data types are approximately equally effective predictors. Performance was more variable between mutations than it was between data types for the same mutation, and there was often little difference between the top data types. We also found that combining data types into a single multi-omics model provided little or no improvement in predictive ability over the best individual data type. Based on our results, for the design of studies focused on the functional outcomes of cancer mutations, there are often multiple -omics types that can serve as effective readouts, although gene expression seems to be a reasonable default option. Extended Graphical Lasso for Multiple Interaction Networks for High Dimensional Omics Data There has been a spate of interest in association networks in biological and medical research, for example, genetic interaction networks. In this paper, we propose a novel method, the extended joint hub graphical lasso (EDOHA), to estimate multiple related interaction networks for high dimensional omics data across multiple distinct classes. To be specific, we construct a convex penalized log likelihood optimization problem and solve it with an alternating direction method of multipliers (ADMM) algorithm. The proposed method can also be adapted to estimate interaction networks for high dimensional compositional data such as microbial interaction networks. The performance of the proposed method in the simulated studies shows that EDOHA has remarkable advantages in recognizing class-specific hubs than the existing comparable methods. We also present three applications of real datasets. Biological interpretations of our results confirm those of previous studies and offer a more comprehensive understanding of the underlying mechanism in disease.

Author summaryReconstruction of multiple association networks from high dimensional omics data is an important topic, especially in biology. Previous studies focused on estimating different networks and detecting common hubs among all classes. Integration of information over different classes of data while allowing difference in the hub nodes is also biologically plausible. Therefore, we propose a method, EDOHA, to jointly construct multiple interaction networks with capacity in finding different hub networks for each class of data. Simulation studies show the better performance over conventional methods. The method has been demonstrated in three real world data.