Hierarchical clustering optimizes the tradeoff between compositionality and expressivity of task structures for flexible reinforcement learning A hallmark of human intelligence, but challenging for reinforcement learning (RL) agents, is the ability to compositionally generalise, that is, to recompose familiar knowledge components in novel ways to solve new problems. For instance, when navigating in a city, one needs to know the location of the destination and how to operate a vehicle to get there, whether it be pedalling a bike or operating a car. In RL, these correspond to the reward function and transition function, respectively. To compositionally generalize, these two components need to be transferable independently of each other: multiple modes of transport can reach the same goal, and any given mode can be used to reach multiple destinations. Yet there are also instances where it can be helpful to learn and transfer entire structures, jointly representing goals and transitions, particularly whenever these recur in natural tasks (e.g., given a suggestion to get ice cream, one might prefer to bike, even in new towns). Prior theoretical work has explored how, in model-based RL, agents can learn and generalize task components (transition and reward functions). But a satisfactory account for how a single agent can simultaneously satisfy the two competing demands is still lacking. Here, we propose a hierarchical RL agent that learns and transfers individual task components as well as entire structures (particular compositions of components) by inferring both through a non-parametric Bayesian model of the task. It maintains a factorised representation of task components through a hierarchical Dirichlet process, but it also represents different possible covariances between these components through a standard Dirichlet process. We validate our approach on a variety of navigation tasks covering a wide range of statistical correlations between task components and show that it can also improve generalisation and transfer in more complex, hierarchical tasks with goal/subgoal structures. Finally, we end with a discussion of our work including how this clustering algorithm could conceivably be implemented by cortico-striatal gating circuits in the brain.
Juvenile Gouldian Finches (Erythrura gouldiae) form sibling sub-groups during social integration The formation of social relationships in complex groups is critical in shaping patterns of social organization and behavioral development. In many birds, young individuals remain dependent on their parents for extended periods but must abruptly transition to navigating interactions in the wider group after independence. While lack of social relationships during this period is detrimental in the development of later social skills, little is known about the social relationships juveniles form after independence in many bird species. In this study we describe patterns of social interactions in juvenile Gouldian Finches (Erythrura gouldiae) during transitions from family groups into flocks of unfamiliar individuals. A group of 20 juveniles from 4 families were introduced into two flocks. After introductions juveniles showed a gradient of approach rates with most approaches directed towards siblings, followed by juvenile peers, adult females, and lastly adult males. Significant preferences for siblings resulted in the emergence of sibling subgroups within the larger social network. This active self-assortment of siblings suggests that sibling sub-groups are an important bridge linking social connections within the family and the wider group. Such findings suggests that sibling relationships have a critical role in the socialization after independence, as well as structuring the social organization of Gouldian finch flocks.
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. Beta-cell specific Insr deletion promotes insulin hypersecretion and improves glucose tolerance prior to global insulin resistance Abstract

Insulin receptor (Insr) protein can be found at higher levels in pancreatic {beta}-cells than in most other tissues, but the consequences of {beta}-cell insulin resistance remain enigmatic. Ins1cre allele was used to delete Insr specifically in {beta}-cells of both female and male mice. Experimental mice were compared to Ins1cre-containing littermate controls at multiple ages and on multiple diets. RNA-seq of purified recombined {beta}-cells revealed transcriptomic consequences of Insr loss, which differed between female and male mice. Action potential and calcium oscillation frequencies were increased in Insr knockout {beta}- cells from female, but not male mice, whereas only male {beta}InsrKO mice had reduced ATP-coupled oxygen consumption rate and reduced expression of genes involved in ATP synthesis. Female {beta}InsrKO and {beta}InsrHET mice exhibited elevated insulin release in perifusion experiments, during hyperglycemic clamps, and following i.p. glucose challenge. Deletion of Insr did not alter {beta}-cell area up to 9 months of age, nor did it impair hyperglycemia-induced proliferation. Based on our data, we adapted a mathematical model to include {beta}-cell insulin resistance, which predicted that {beta}-cell Insr knockout would improve glucose tolerance depending on the degree of whole-body insulin resistance. Indeed, glucose tolerance was significantly improved in female {beta}InsrKO and {beta}InsrHET mice when compared to controls at 9, 21 and 39 weeks, and also in insulin-sensitive 4-week old males. We did not observe improved glucose tolerance in older male mice or in high fat diet-fed mice, corroborating the prediction that global insulin resistance obscures the effects of {beta}-cell specific insulin resistance. The propensity for hyperinsulinemia was associated with mildly reduced fasting glucose and increased body weight. We further validated our main in vivo findings using the Ins1-CreERT transgenic line and found that male mice had improved glucose tolerance 4 weeks after tamoxifen-mediated Insr deletion. Collectively, our data show that loss of {beta}-cell Insr contributes to glucose-induced hyperinsulinemia, thereby improving glucose homeostasis in otherwise insulin sensitive dietary and age contexts. New multiscale characterisation methodology for effective determination of isolation-structure-function relationship of extracellular vesicles. Extracellular vesicles (EVs) have been lauded as next generation medicines, but very few EV-based therapeutics have progressed to clinical use. Limited clinical translation is largely due to technical barriers that hamper our ability to mass-produce EVs, i.e. to isolate, purify and characterise them effectively. Technical limitations in comprehensive characterisation of EVs leads to unpredicted biological effects of EVs. Here, using a range of optical and non-optical techniques, we showed that the differences in molecular composition of EVs isolated using two isolation methods correlated with the differences in their biological function. Our results demonstrated that the isolation method determines the composition of isolated EVs at single and sub-population levels. Besides the composition, we measured for the first time the dry mass and predicted sedimentation of EVs. These parameters were shown to correlate well with the biological and functional effects of EVs on single cell and cell cultures. We anticipate that our new multiscale characterisation approach, which goes beyond traditional experimental methodology, will support fundamental understanding of EVs as well as elucidate the functional effects of EVs in in vitro and in vivo studies. Our findings and methodology will be pivotal for developing optimal isolation methods and establishing EVs as mainstream therapeutics and diagnostics. This innovative approach is applicable to a wide range of sectors including biopharma and biotechnology as well as to regulatory agencies. Morphology of immatures of the thelytokous ant, Monomorium triviale Wheeler (Formicidae: Myrmicinae: Solenopsidini) with descriptions of the extraordinary last-instar queen larvae The ant genus Monomorium is one of the most species-rich but taxonomically problematic groups in the hyperdiverse subfamily Myrmicinae. An East Asian species, M. triviale Wheeler, produces both reproductive queens and sterile workers via obligate thelytokous parthenogenesis. Here, we describe the immature forms of M. triviale based on light and scanning electron microscopy observations, with a note on the striking caste dimorphism in the last larval instar. The last-instar queen larvae were easily recognized by their large size, "aphaenogastroid" body shape, and rows of doorknob-like tubercles on the lateral and dorsal body surface. This type of queen-specific structure has not been found in ants in general, let alone congeneric species found in Japan. In stark contrast to the queen larvae, worker larvae showed a "pheidoloid" body shape and a body surface similar to other ants. The worker larvae were estimated to have three instars, consistent with previously described congeners. The pupae of both castes had no cocoon, a characteristic commonly described in other Myrmicinae species. In total, the developmental period from egg to adult worker averaged 59 days under 25{degrees}C. We discuss possible functions of the tubercles of queen larvae based on previous studies. Inhibitory proteins block substrate access by occupying the active site cleft of Bacillus subtilis intramembrane metalloprotease SpoIVFB Intramembrane proteases function in numerous signaling pathways that impact health, but knowledge about regulation of intramembrane proteolysis is limited. We examined inhibition of intramembrane metalloprotease SpoIVFB by proteins BofA and SpoIVFA. We found that BofA residues in and near a predicted transmembrane segment are required for SpoIVFB inhibition, and cross-linking experiments indicated that this transmembrane segment occupies the SpoIVFB active site cleft. SpoIVFA is also required for SpoIVFB inhibition. The inhibitory proteins block access of the substrate N-terminal Proregion to the membrane-embedded SpoIVFB active site, based on additional cross-linking experiments; however, the inhibitory proteins did not prevent interaction between the substrate C-terminal region and the SpoIVFB soluble domain. A structural model was built of SpoIVFB in complex with BofA and parts of SpoIVFA and substrate, using partial homology and constraints from cross-linking and co-evolutionary analyses. The model predicts that conserved BofA residues interact to stabilize a transmembrane segment and a membrane-embedded C-terminal region. SpoIVFA is predicted to bridge the BofA C-terminal region and SpoIVFB, forming a membrane-embedded inhibition complex. Our results reveal a novel mechanism of intramembrane protease inhibition with clear implications for relief from inhibition in vivo and design of inhibitors as potential therapeutics. Preprints in motion: tracking changes between preprint posting and journal publication during a pandemic Amidst the COVID-19 pandemic, preprints in the biomedical sciences are being posted and accessed at unprecedented rates, drawing widespread attention from the general public, press and policymakers for the first time. This phenomenon has sharpened longstanding questions about the reliability of information shared prior to journal peer review. Does the information shared in preprints typically withstand the scrutiny of peer review, or are conclusions likely to change in the version of record? We assessed preprints from bioRxiv and medRxiv that had been posted and subsequently published in a journal through 30th April 2020, representing the initial phase of the pandemic response. We utilised a combination of automatic and manual annotations to quantify how an article changed between the preprinted and published version. We found that the total number of figure panels and tables changed little between preprint and published articles. Moreover, the conclusions of 7.2% of non-COVID-19-related and 17.2% of COVID-19-related abstracts undergo a discrete change by the time of publication, but the majority of these changes do not qualitatively change the conclusions of the paper. Plant inositol-phosphate-glycans and a fucosylated xyloglucan oligosaccharide are accumulated upon Arabidopsis thaliana/ Botrytis cinerea infection In mammals, insulin is involved in controlling blood glucose levels and its role in modulating immunity is being more and more documented. This hormone promotes the release of inositolphosphate glycans (IPG) which act as mediators. In plants, one IG has already been identified in plant culture cells (Smith and Fry, 1999; Smith et al., 1999) but, to our knowledge, no IPG have been yet identified. Here, we discovered 7 IPG that are accumulated upon Arabidopsis thaliana-Botrytis cinerea interaction, concomitantly with oligogalacturonides and a fucosylated xyloglucan oligosaccharide. Further structural characterization showed that they come from the hydrolysis of polar heads of Serie A to H glycosylinositol phosphorylceramides presumably via a phospholipase C activity. Taken together with the emerging role of insulin as immune regulator, these results question the role of IPG as damage associated molecular pattern both in animal and plant kingdoms. Allopatric divergence of cooperators limits cheating range and effects of a defector mutation Social defectors may meet diverse cooperators. Genotype-by-genotype interactions may constrain the ranges of cooperators upon which particular defectors can cheat, limiting cheater spread. The bacterium Myxococcus xanthus undergoes cooperative multicellular development, but some developmental defectors cheat on cooperators during this process. In this study, interactions between a cheater disrupted at the signaling gene csgA and allopatrically diversified cooperators reveal a very small cheating range. Expectedly, the cheater failed to cheat on all natural-isolate cooperators owing to non-cheater-specific antagonisms. Surprisingly, lab-evolved cooperators that diverged from their cheating-susceptible ancestor by fewer than 20 mutations and without experiencing cheating had already exited the csgA mutants cheating range. Cooperators might also diversify in the potential for a mutation to reduce expression of cooperative trait or generate a cheating phenotype. A new csgA mutation constructed in several highly diverged cooperators generated diverse sporulation phenotypes, ranging from a complete defect to no defect, indicating that genetic backgrounds can limit the set of genomes in which a mutation creates a defector. Our results suggest that natural populations feature geographic mosaics of cooperators that have diversified in their susceptibility to particular cheaters and in the phenotypes generated by any given cooperation-gene mutation.

Significance statementSelection on cooperators exploited by obligate cheaters can induce evolution of resistance to cheating. Here we show that cooperators can also rapidly evolve immunity to cheating simply as a byproduct of evolutionary divergence in environments in which cooperation and cheating at the focal trait do not occur because the trait is not expressed. We also find that differences in the genomic context in which a cooperation-gene mutation arises can profoundly alter its phenotypic effect and determine whether the mutation generates a social defect at all - a pre-requisite for obligate cheating. These findings suggest that general divergence of social populations under a broad range of environmental conditions can restrict both the set of mutations that might generate social defectors in the first place and the host range of such defectors once they arise. Cargo-mediated mechanisms reduce inter-motor mechanical interference, promote load-sharing and enhance processivity in teams of molecular motors In cells, multiple molecular motors work together as teams to carry cargoes such as vesicles and organelles over long distances to their destinations by stepping along a network of cytoskeletal filaments. How motors that typically mechanically interfere with each other, work as teams is unclear. Here we explored the possibility that purely physical mechanisms may potentially enhance teamwork, both at the single motor and cargo level. To explore these mechanisms, we developed a 3D dynamical simulation of cargo transport along microtubules by teams of canonically non-cooperative kinesin-1 motors. We accounted for cargo membrane fluidity by explicitly simulating the Brownian dynamics of motors on the cargo surface and considered both the load and ATP dependence of single motor functioning. We showed explicitly, for the first time, that surface fluidity leads to the reduction of negative mechanical interference between kinesins, enhancing load sharing thereby decreasing single motor off-rates and increasing processivity. Remarkably, we also showed that, independent of fluidity, increasing numbers of bound motors pull the cargo closer to the microtubule, increasing the on-rates of individual unbound motors, resulting in a cooperative increase in bound motor numbers that depends on 3D cargo geometry. At the cargo level, surface fluidity makes more motors available for binding, though this effect is significant only at low ATP or high motor density. Interestingly, we find that the fluidity induced reduction in mechanical interference dominates over the increased availability of motors for typical physiological and in vitro conditions, this allowing us to reconcile different experimental results in different regimes. Finally, we show that these effects can altogether result in enhanced mechanical efficiency and tunable cargo run-lengths for teams of molecular motors over physiological ranges of fluidity with implications for new experimental validation efforts, transport in vivo, artificial cargo design and multi-motor-driven mechanical processes in general.

Author summaryIn cells, multiple molecular motors work together as teams to carry cargoes such as vesicles and organelles over long distances to their destinations by stepping along a network of protein filaments. How do intrinsically non-cooperative motors function as teams within cells? In this paper, we show, using computer simulations, that the fluid surfaces of cellular cargo reduce the mechanical interference between motors allowing better load sharing and decreasing their detachment rates, thereby increasing the distance over which they can carry cargo. We additionally found that motors pull cargo closer to the filament increasing the attachment rates of other unbound motors. These effects act synergistically with an increased availability of motors due to fluidity to further increase travel distances. Our simulation model accounted for the dynamics of the motors on the cargo surface and their load and ATP dependent kinetics, allowing us to connect single motor properties to overall transport. Our work on understanding how teamwork arises in mechanically coupled motors sheds new light on cellular processes, reconciles existing observations, encourages new experimental validation efforts and can also suggest new ways of improving mechanical efficiency and transport of artificial cargo powered by motor teams. Mathematical Modeling and Analysis of Mitochondrial Retrograde Signaling Dynamics Mitochondria, semi-autonomous eukaryotic organelles, participate in energy production and metabolism, making mitochondrial quality control crucial. As most mitochondrial proteins are encoded by nuclear genes, quality control depends on proper mitochondria-nucleus communication, designated mitochondrial retrograde signaling. Early studies focused on retrograde signaling participants and specific gene knockouts. However, mitochondrial signal modulation remains elusive. Using yeast, we simulated signal propagation following mitochondrial damage and proposed a mathematical model based ordinary differential equations. Mitochondrial retrograde signaling decisions were described using a Boolean model. Dynamics were analyzed through an ordinary differential equation-based model and extended to evaluate the model response to noisy damage signals. Simulation revealed localized protein concentration dynamics, including waveforms, frequency response, and robustness under noise. Retrograde signaling is bistable with three localized steady states, and increased damage compromises robustness. We elucidated mitochondrial retrograde signaling, thus providing a basis for drug design against yeast and fungi.

Author SummaryThe yeast RTG pathway regulates mitochondrial metabolism and mitochondrial quality through passing mitochondrial signal to the nucleus to modulate gene expressions. Using microscopic data of RTG proteins translocation from cytosol to nucleus, the parameters were found by fitting 16 knockout conditions simulated from the proposed differential equation-based model. Further, through dose response, ultrasensitivity, frequency response and noise interference, we demonstrated a switch-like property of RTG pathway activation, a capacity charging property of RTG protein translocation, a low pass filter property of signals in different frequencies, and noise amplification due to mitochondrial signal. The proposed mathematical model enables us to understand the dynamics and mechanisms in mitochondrial retrograde signaling in yeast and provides potential antifungal treatment strategies. Slow Cortical Waves through Cyclicity Analysis Fine-grained understanding of dynamics in cortical networks is crucial in unpacking brain function. Here, we introduce a novel analytical method to characterize the dynamic interaction between distant brain regions, and apply it to data from the Human Connectome Project.

Resting-state fMRI results in time series recordings of the activity of different brain regions, which are aperiodic and lacking a base frequency. Cyclicity Analysis, a novel technique robust with respect to time-reparametrizations, is effective in recovering temporal ordering of such time series along a circular trajectory without assuming any time-scale. Our analysis detected slow cortical waves of activity propagating across the brain with consistent lead-lag relationships between specific brain regions. We also observed short bursts of task-modulated strong temporal ordering that dominate overall lead-lag relationships between pairs of regions in the brain. Our results suggest the possible role played by slow waves of information transmission between brain regions that underlie emergent cognitive function. Systems-wide analysis of the GATC-binding nucleoid-associated protein Gbn and its impact on Streptomyces development Bacterial chromosome structure is organized by a diverse group of proteins collectively referred to as nucleoid-associated proteins (NAPs). Many NAPs have been well studied in Streptomyces, including Lsr2, HupA, HupS, and sIHF. Here, we show that SCO1839 represents a novel family of Actinobacteria NAPs and recognizes a consensus sequence consisting of GATC followed by (A/T)T. The protein was designated Gbn for GATC-binding NAP. Deletion of gbn led to alterations in development and antibiotic production in Streptomyces coelicolor. Chromatin immunoprecipitation sequencing (ChIP-Seq) detected more than 2800 binding regions, encompassing some 3600 GATCWT motifs, which comprise 55% of all such motifs in the S. coelicolor genome. DNA binding of Gbn in vitro increased DNA stiffness but not compaction, suggesting a role in regulation rather than chromosome organization. Transcriptomics analysis showed that Gbn binding generally leads to reduced gene expression. The DNA binding profiles were nearly identical between vegetative and aerial growth. Exceptions are SCO1311 and SCOt32, for a tRNA editing enzyme and a tRNA that recognises the rare leucine codon CUA, respectively, which nearly exclusively bound during vegetative growth. Taken together, our data show that Gbn is a highly pleiotropic NAP that impacts growth and development in streptomycetes.

IMPORTANCEA large part of the chemical space of bioactive natural products is derived from Actinobacteria. Many of the biosynthetic gene clusters for these compounds are cryptic, in others words, they are expressed in nature but not in the laboratory. Understanding the global regulatory networks that control gene expression is key to the development of approaches to activate this biosynthetic potential. Chromosome structure has a major impact on the control of gene expression. In bacteria, the organization of chromosome structure is mediated by a diverse group of proteins referred to collectively as nucleoid-associated proteins (NAPs), which play an important role in the control of gene expression, nucleoid structure and DNA repair. We here present the discovery of a novel and extremely pleiotropic NAP, which we refer to as Gbn. Gbn is a sporulation-specific protein that occurs only in the Actinobacteria and binds to GATC sequences, with a subtle but broad effect on global gene expression. The discovery of Gbn is a new step towards better understanding of how gene expression and chromosome structure is governed in antibiotic-producing streptomycetes. Changes to the mtDNA copy number during yeast culture growth We show that the mtDNA copy number in growing cultures of the yeast Saccharomyces cerevisiae increases by a factor of up to 4, being lowest ([~]10 per haploid genome) and stable during rapid fermentative growth, and highest at the end of the respiratory phase. When yeast are grown on glucose, the onset of the mtDNA copy number increase coincides with the early stages of the diauxic shift, and the increase continues through respiration. A lesser yet still substantial copy number increase occurs when yeast are grown on a nonfermentable carbon source, i.e. when there is no diauxic shift. The mtDNA copy number increase during and for some time after the diauxic shift is not driven by an increase in cell size. The copy number increase occurs in both haploid and diploid strains, but is markedly attenuated in a diploid wild isolate that is a ready sporulator. Strain-to-strain differences in mtDNA copy number are least apparent in fermentation and most apparent in late respiration or stationary phase. While changes in mitochondrial morphology and function were previously known to accompany changes in physiological state, it had not been previously shown that the mtDNA copy number changes substantially over time in a clonal growing culture. The mtDNA copy number in yeast is therefore a highly dynamic phenotype. Fiber orientation downsampling compromises the computation of white matter tract-related deformation Incorporating neuroimaging-revealed structural details into finite element (FE) head models opens vast new opportunities to better understand brain injury mechanisms. Recently, growing efforts have been made to integrate fiber orientation from diffusion tensor imaging (DTI) into FE models to predict white matter (WM) tract-related deformation that is biomechanically characterized by tract-related strains. Commonly used approaches often downsample the spatially enriched fiber orientation to match the FE resolution with one orientation per element (i.e., element-wise orientation implementation). However, the validity of such downsampling operation and corresponding influences on the computed tract-related strains remain elusive. To address this, the current study proposed a new approach to integrate voxel-wise fiber orientation from one DTI atlas (isotropic resolution of 1 mm3) into FE models by embedding orientations from multiple voxels within one element (i.e., voxel-wise orientation implementation). By setting the responses revealed by the newly proposed voxel-wise orientation implementation as the reference, we evaluated the reliability of two previous downsampling approaches by examining the downsampled fiber orientation and the computationally predicted tract-related strains secondary to one concussive impact. Two FE models with varying element sizes (i.e., 6.37 {+/-} 1.60 mm and 1.28 {+/-} 0.55 mm, respectively) were incorporated. The results showed that, for the model with a large voxel-mesh resolution mismatch, the downsampled element-wise fiber orientation, with respect to its voxel-wise counterpart, exhibited an absolute deviation over 30{degrees} across the WM/gray matter interface and the pons regions. Accordingly, this orientation deviation compromised the computation of tract-related strains with normalized root-mean-square errors up to 30% and underestimated the peak tract-related strains up to 10%. For the other FE model with finer meshes, the downsampling-induced effects were lower, both on the fiber orientation and tract-related strains. Taken together, the voxel-wise orientation implementation is recommended in future studies as it leverages the DTI-delineated fiber orientation to a larger extent than the element-wise orientation implementation. Thus, this study yields novel insights on integrating neuroimaging-revealed fiber orientation into FE models and may better inform the computation of WM tract-related deformation, which are crucial for advancing the etiological understanding and computational predictability of brain injury.

Graphic abstract

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=142 SRC="FIGDIR/small/471622v2_ufig1.gif" ALT="Figure 1">
View larger version (55K):
org.highwire.dtl.DTLVardef@12ebda9org.highwire.dtl.DTLVardef@1a5d624org.highwire.dtl.DTLVardef@8e068borg.highwire.dtl.DTLVardef@4f8ff5_HPS_FORMAT_FIGEXP  M_FIG C_FIG Loci for prediction of penicillin and tetracycline susceptibility in Neisseria gonorrhoeae: a genome wide association study BackgroundWhile Neisseria gonorrhoeae poses an urgent public health threat because of increasing antimicrobial resistance, much of the circulating population remains susceptible to historical treatment regimens. Point-of-care diagnostics that report susceptibility could allow for reintroduction of these regimens, but development of such diagnostics has been limited to ciprofloxacin, for which susceptibility can be predicted from a single locus.

MethodsWe assembled a dataset of 12,045 N. gonorrhoeae genomes with phenotypic resistance data for tetracycline (n = 3,611) and penicillin (n = 6,935). Using conditional genome wide association studies (GWAS), we sought to define genetic variants associated with susceptibility to penicillin and tetracycline. We evaluated the sensitivity and specificity of these variants for predicting susceptibility and non-resistance in our collection of gonococcal genomes.

FindingsIn our conditional penicillin GWAS, the presence of a genetic variant defined by a non-mosaic penA allele without an insertion at codon 345 was significantly associated with penicillin susceptibility and had the highest negative effect size of significant variants (p = 5.0 x 10-14, {beta} = -2.5). In combination with the absence of blaTEM, this variant predicted penicillin susceptibility with high specificity (99.8%) and modest sensitivity (36.7%). For tetracycline, the wild type allele at rpsJ codon 57, encoding valine, was significantly associated with tetracycline susceptibility (p = 5.6 x 10-16, {beta} = -1.6) after conditioning on the presence of tetM. The combination of rpsJ codon 57 allele and tetM absence predicted tetracycline susceptibility with high specificity (97.2%) and sensitivity (88.7%).

InterpretationAs few as two genetic loci can predict susceptibility to penicillin and tetracycline in N. gonorrhoeae with high specificity. Molecular point-of-care diagnostics targeting these loci have the potential to increase available treatments for gonorrhea.

FundingNational Institute of Allergy and Infectious Diseases, the National Science Foundation, and the Smith Family Foundation

Research in ContextO_ST_ABSEvidence before this studyC_ST_ABSWe searched PubMed with the terms "Neisseria gonorrhoeae" and "diagnostic" or "assay" plus "penicillin" or "tetracycline" for reports in any language published up to July 1, 2021. We additionally searched for "Neisseria gonorrhoeae" and "genome wide association study". We found that previously proposed molecular diagnostics for penicillin and tetracycline susceptibility either exclusively focused on plasmid-mediated resistance (i.e., targeting blaTEM or tetM) or did not include variants in genes encoding antibiotic targets (e.g., did not include penA or rpsJ). Targets for molecular surveillance have focused on resistance-associated alleles rather than susceptibility-associated alleles. We did not find any previous penicillin or tetracycline conditional genome wide association studies (GWAS) in N. gonorrhoeae.

Added value of this studyTo identify targets for molecular diagnostics that predict penicillin and tetracycline susceptibility, we conducted GWAS conditioning on the presence of plasmid-mediated resistance determinants to detect chromosomal loci with the highest association with susceptibility. We discovered a sequence (penA_01) that differentiates susceptible isolates from those with a resistance-associated insertion at codon 345 and from those with mosaic penA alleles, which is associated with penicillin susceptibility. We also found that rpsJ codon 57 was the chromosomal locus contributing the most to tetracycline susceptibility. The combination of these chromosomal loci and the absence of plasmid encoded determinants predicts penicillin and tetracycline susceptibility with high specificity in both a large global collection of N. gonorrhoeae and a validation dataset consisting of recently published genomes from CDCs Gonococcal Isolate Surveillance Program (GISP) surveillance collected in 2018.

Implications of all the available evidenceThe chromosomal loci penA_01 and rpsJ codon 57 in combination with plasmid loci blaTEM and tetM are candidates for the development of point-of-care molecular diagnostics for penicillin and tetracycline susceptibility. The loci may be combined with the currently available ciprofloxacin susceptibility diagnostics to predict susceptibility to multiple antibiotics. Additionally, our study suggests that conditional GWAS focused on variants associated with susceptibility may be a promising approach to identify minimal sets of loci for molecular diagnostics and surveillance. Host Variation in Type I Interferon Signaling Genes (MX1), CCR5{triangleup}32, and MHC Class I Alleles in Treated HIV+ Non-Controllers Predict Viral Reservoir Size Prior genomewide association studies have identified variation in MHC Class I alleles and CCR5{Delta}32 as genetic predictors of viral control, especially in "elite" controllers, individuals who remain virally suppressed in the absence of therapy. We analyzed custom host whole exome sequencing and direct HLA typing from 202 ART-suppressed HIV+ non-controllers in relation to four measures of the peripheral CD4+ T cell reservoir: HIV intact DNA, total (t)DNA, unspliced (us)RNA, and RNA/DNA. Linear mixed models were adjusted for potential covariates including age, sex, nadir CD4+ T cell count, pre-ART HIV RNA, timing of ART initiation from date of detected HIV infection, and duration of ART suppression. Previously reported HLA "protective" B*57:01 (in untreated elite controllers) was significantly associated with lower usRNA (p=4.2x10-4, Benjamini-Hochberg false discovery rate [FDR] q=3.3x10-3). The favorable CCR5{Delta}32 mutation was also associated with lower tDNA and usRNA, primarily in Europeans (p=4.3x10-3 and 8.7x10-3). Several SNPs in interferon stimulated gene, MX1, were associated with HIV total DNA (p=1.3x10-7; q=0.02), and the direction of these associations paralleled predicted eQTL (expression quantitative trait) MX1 gene expression. Among Europeans, gene set enrichment analysis identified genes involved in "negative regulation of host cell viral entry" to be significantly associated with usRNA (q=0.03). These findings suggest a potential role for interferon signaling, MHC class I alleles, and CCR5{Delta}32 in maintaining the reservoir size in HIV non-controllers. These data need to be further validated in functional studies and replication cohorts. 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. A network approach to genetic circuit designs As genetic circuits become more sophisticated, the size and complexity of data about their designs increases. This data captured goes beyond monolithic genetic sequences and towards circuit modularity and functional details, which are beneficial for analyzing circuit performance and establishing design automation techniques. However, the accessibility, visualisation and usability of design data (and metadata) have received relatively little attention to date. Here, we present a method to turn circuit designs into networks and showcase its potential to enhance the utility of design data. Since networks are dynamic structures, initial graphs can be interactively shaped into sub-networks of relevant information based on requirements such as abstraction, hierarchy and protein interactions. Additionally, several visual changes can be applied, such as colouring or clustering nodes based on types (e.g., genes or promoters), resulting in easier comprehension from a user perspective. This approach allows circuit designs to be coupled to other networks, such as metabolic pathways or implementation protocols captured in graph-like formats. Therefore, we advocate using networks to structure, access and improve synthetic biology information. LY-CoV1404 (bebtelovimab) potently neutralizes SARS-CoV-2 variants SARS-CoV-2 neutralizing monoclonal antibodies (mAbs) can reduce the risk of hospitalization when administered early during COVID-19 disease. However, the emergence of variants of concern has negatively impacted the therapeutic use of some authorized mAbs. Using a high throughput B-cell screening pipeline, we isolated a highly potent SARS-CoV-2 spike glycoprotein receptor binding domain (RBD)-specific antibody called LY-CoV1404 (also known as bebtelovimab). LY-CoV1404 potently neutralizes authentic SARS-CoV-2 virus, including the prototype, B.1.1.7, B.1.351 and B.1.617.2). In pseudovirus neutralization studies, LY-CoV1404 retains potent neutralizing activity against numerous variants including B.1.1.7, B.1.351, B.1.617.2, B.1.427/B.1.429, P.1, B.1.526, B.1.1.529, and the BA.2 subvariant and retains binding to spike proteins with a variety of underlying RBD mutations including K417N, L452R, E484K, and N501Y. Structural analysis reveals that the contact residues of the LY-CoV1404 epitope are highly conserved with the exception of N439 and N501. Notably, the binding and neutralizing activity of LY-CoV1404 is unaffected by the most common mutations at these positions (N439K and N501Y). The breadth of reactivity to amino acid substitutions present among current VOC together with broad and potent neutralizing activity and the relatively conserved epitope suggest that LY-CoV1404 has the potential to be an effective therapeutic agent to treat all known variants causing COVID-19.

In BriefLY-CoV1404 is a potent SARS-CoV-2-binding antibody that neutralizes all known variants of concern and whose epitope is rarely mutated.

HighlightsO_LILY-CoV1404 potently neutralizes SARS-CoV-2 authentic virus and known variants of concern including the B.1.1.529 (Omicron), the BA.2 Omicron subvariant, and B.1.617.2 (Delta) variants
C_LIO_LINo loss of potency against currently circulating variants
C_LIO_LIBinding epitope on RBD of SARS-CoV-2 is rarely mutated in GISAID database
C_LIO_LIBreadth of neutralizing activity and potency supports clinical development
C_LI Co-infection of Candidatus Piscichlamydia trichopodus (order Chlamydiales) and Henneguya sp (Myxosporea, Myxobolidae) in snakeskin gourami Trichopodus pectoralis (Regan 1910) The present study describes a simultaneous infection of a novel Chlamydia-like organism (CLO) with a Myxozoa parasite, Henneguya sp. in snakeskin gourami Trichopodus pectoralis in Thailand. A new CLO is proposed "Candidatus Piscichlamydia trichopodus" (CPT) based on 16S rRNA phylogenetic analysis. Systemic intracellular CPT infection was confirmed by histological examination, in situ hybridization, PCR assay, and sequencing of 16S rRNA. This novel pathogen belongs to the order Chlamydiales but differs in certain aspects from other species. The histopathological changes associated with CPT infection were different from the typical pathological lesions of epitheliocystis caused by previously known CLO. Unlike other CLO, CPT localized in the connective tissue rather than in the epithelial cells and formed smaller clumps of intracellular bacteria that stained dark blue with hematoxylin. On the other hand, typical myxospores of the genus Henneguya with tails were observed in the gill sections. Infection with Henneguya sp. resulted in extensive destruction of the gill filaments, most likely leading to respiratory distress. Due to the frequency of co-infections and the unavailability of culture methods for CLO and Henneguya sp, it was difficult to determine which pathogens were directly responsible for the associated mortality. However, co-infections may increase the negative impact on the host and the severity of the disease. Given the commercial importance of the snakeskin gourami and its significant aquaculture potential, the findings of this study are important for further studies on disease prevention. Selection towards lighter seeds in Danish grasslands over an eight-year period Time-series analysis of cover data from 236 Danish grassland sites demonstrated community selection on seed mass. Across four grassland habitat types and during an eight-year period there was significant community selection towards lighter seeds. In the analysis, the continuous seed size variable was used to group plant species into functional types. This method was chosen in order to account for the sampling process of the abundance data, and is in contrast to most other analyses of trait selection, where the community weighted mean of the traits is used as the dependent variable. Generally, little is known on the underlying causes of possible plant trait community selection in semi-natural grasslands that are not subjected to successional processes. A non-canonical convergence of carbohydrate and glutamine metabolism is required after metabolic rewiring in a solid environment The fate of pyruvate, which is modulated by the activity of the mitochondrial pyruvate carrier (MPC), is a defining metabolic feature in many cancers. Diffuse large B-cell lymphomas (DLBCLs) are a genetically and metabolically heterogeneous cancer. Although MPC expression and activity differed between DLBCL subgroups, mitochondrial pyruvate was uniformly consumed by glutamate pyruvate transaminase 2 (GPT2) to support -ketoglutarate production as part of glutaminolysis. This led us to discover that glutamine exceeds pyruvate as a carbon source for the tricarboxylic acid (TCA) cycle in DLBCLs. Furthermore, we found that MPC inhibition unexpectedly leads to decreased glutaminolysis, which is contrary to previous observations in other cell types. We also discovered that MPC inhibition and depletion only decreased DLBCL proliferation in an extracellular matrix (ECM) environment and in vivo xenografts, but not in the typical DLBCL suspension environment. We also have found that the metabolic profile of DLBCL cells in ECM is markedly different from cells in suspension environment. Thus, we report that besides the canonical glutamate dehydrogenase (GDH)-mediated glutaminolysis, the non-canonical GPT2 mediated consumption and assimilation of glutamine and pyruvate in DLBCLs enables their proliferation in an extracellular environment-dependent manner.

HIGHLIGHTSO_LIGlutamine, but not glucose, is a major carbon source for the tricarboxylic acid cycle in diffuse large B-cell lymphomas.
C_LIO_LIMitochondrial pyruvate supports glutaminolysis in diffuse large B-cell lymphomas by supplying pyruvate for glutamate pyruvate transaminase 2-mediated  -ketoglutarate production.
C_LIO_LIMitochondrial pyruvate carrier inhibition leads to decreased glutaminolysis in diffuse large B-cell lymphomas.
C_LIO_LI -ketoglutarate production is important for diffuse large B-cell lymphoma proliferation in a solid extracellular matrix environment.
C_LIO_LIMitochondrial pyruvate carrier activity supports diffuse large B-cell lymphoma proliferation in a solid extracellular matrix environment and in mouse xenografts.
C_LI β-catenin and canonical Wnts in Hydra pattern formation: Insights from mathematical modelling into two distinct systems Formation of the body axes and the subsequent formation of the apical termini are fundamental steps during animal development. Nuclear {beta}-catenin and canonical Wnt3 have been identified as major players in both these processes in Hydra. Based on molecular knowledge of canonical Wnt signalling directly linking nuclear {beta}-catenin and Wnt3 activity, de novo axis formation and head formation have frequently been assumed to be part of the same pattern formation system. In this study, we investigated this relationship by combining new model simulations with available experimental results. We demonstrated that nuclear {beta}-catenin and Wnt3 likely contribute to two distinct de novo pattern formation systems in Hydra, organizing development and differentiation on two different spatial scales. In particular, the results suggested that nuclear {beta}-catenin acts at the scale of the whole body to control axis formation, whereas canonical Wnt3 signalling is involved in a downstream pathway responsible for small-scale patterning of the head. These results suggest that the underlying mechanisms may also be more complex than previously assumed in other animals where axis formation has been ascribed to canonical Wnt signalling. ASURAT: functional annotation-driven unsupervised clustering of single-cell transcriptomes MotivationSingle-cell RNA sequencing (scRNA-seq) analysis reveals heterogeneity and dynamic cell transitions. However, conventional gene-based analyses require intensive manual curation to interpret the biological implications of computational results. Hence, a theory for efficiently annotating individual cells is necessary.

ResultsWe present ASURAT, a computational pipeline for simultaneously performing unsupervised clustering and functional annotation of disease, cell type, biological process, and signaling pathway activity for single-cell transcriptomic data, using correlation graph-based decomposition of genes based on database-derived functional terms. We validated the usability and clustering performance of ASURAT using scRNA-seq datasets for human peripheral blood mononuclear cells, which required fewer manual curations than existing methods. Moreover, we applied ASURAT to scRNA-seq and spatial transcriptome datasets for small cell lung cancer and pancreatic ductal adenocarcinoma, identifying previously overlooked subpopulations and differentially expressed genes. ASURAT is a powerful tool for dissecting cell subpopulations and improving biological interpretability of complex and noisy transcriptomic data.

AvailabilityA GPLv3-licensed implementation of ASURAT is on GitHub (https://github.com/keita-iida/ASURAT).