Increased biting rate and decreased Wolbachia density in irradiated mosquitoes BackgroundReleasing considerable numbers of radiation-sterilized males is a promising strategy to suppress mosquito vectors. However, releases may also include small percentages of biting females which translate to large numbers when releases are large. Currently, the effects of irradiation on the host-seeking and host-biting behaviors have not been exhaustively investigated. Information is also lacking regarding the effects of sterilizing treatment on the endosymbiotic bacterium Wolbachia, which is known to affect the vector competence of infected mosquitos.

MethodsTo ascertain the effects of irradiation on females, the pupae of two Aedes albopictus strains, differing in their natural or artificial Wolbachia infection type, and Ae. aegypti--which is not infected by Wolbachia--were treated with various doses of X-rays and monitored for key fitness parameters and biting behavior over a period of two weeks. The effect of radiation on Wolbachia was investigated by qPCR and FISH analysis.

ResultsPartial Ae. albopictus female sterility was achieved at 28 Gy but the number of weekly bites more than doubled compared to that of the controls. Radiation doses of 35 and 45 Gy completely inhibited progeny production but did not significantly affect the survival or flight ability of Ae. albopictus females and caused a tripling of the number of bites per female per week (compared to untreated controls). These results were also confirmed in Ae. aegypti after treatment at 50 Gy. Wolbachia density decreased significantly in 45-Gray-irradiated females, with the greatest decreases in the early irradiation group (26{+/-}2-hour-old pupae). Wolbachia density also decreased as adults aged. This trend was confirmed in ovaries but not in extra-ovarian tissues. FISH analysis showed a strongly reduced Wolbachia-specific fluorescence in the ovaries of 13{+/-}1-day-old females.

ConclusionsThese results suggest that, under SIT programs, the vector capacity of a target population could increase with the frequency of the irradiated females co-released with the sterile males due to an increased biting rate. In the context of a successful suppression, the related safety issues could be generally negligible, but they should be conservatively evaluated when large scale programs relying on imperfect sexing and high overflooding release ratios are run for long time in areas endemic for arboviral diseases. Also, the effects of irradiation on the vector competence deserve further investigation.
Metabolic reconstitution by a gnotobiotic microbiota varies over the circadian cycle The capacity of the intestinal microbiota to degrade otherwise indigestible diet components is known to greatly improve the recovery of energy from food. This has led to the hypothesis that increased digestive efficiency may underlie the contribution of the microbiome to obesity. OligoMM12-colonized gnotobiotic mice have a consistently higher fat-mass than germ-free or fully colonized counterparts. We therefore investigated their food intake, digestion efficiency, energy expenditure and respiratory quotient using a novel isolator-housed metabolic cage system which allows long-term measurements without contamination risk. This demonstrated that microbiota- released calories are perfectly balanced by decreased food intake in fully colonized versus gnotobiotic OligoMM12 and germ-free mice fed a standard chow diet, i.e., microbiota-released calories can in fact be well-integrated into appetite control. We also observed no significant difference in energy expenditure per gram lean mass between the different microbiota groups, suggesting that cumulative very small differences in energy balance, or altered energy storage must underlie fat accumulation in OligoMM12 mice. Consistent with altered energy storage, major differences were observed in the type of respiratory substrates used in metabolism over the circadian cycle: in germ-free mice the respiratory exchange ratio was consistently lower than that of fully colonized mice at all times of day, indicative of more reliance on fat and less on glucose metabolism. Intriguingly the RER of OligoMM12-colonized gnotobiotic mice phenocopied fully colonized mice during the dark (active/eating) phase but phenocopied germ-free mice during the light (fasting/resting) phase. Further, OligoMM12-colonized mice showed a germ-free-like drop in liver glycogen storage during the light cycle and both liver and plasma metabolomes of OligoMM12 mice clustered closely with germ-free mice. This implies the existence of microbiota functions that are required to maintain normal host metabolism during the resting/fasting phase of circadian cycle, and which are absent in the OligoMM12 consortium.
Bioactive constituents of Verbena officinalis regulate killing efficiency of primary human natural killer cells by accelerating killing processes BackgroundNatural killer (NK) cells play a key role in eliminating tumorigenic and pathogen-infected cells. Verbena officinalis (V. officinalis) has been used as a medical plant in traditional and modern medicine, exhibiting anti-tumor and anti-inflammation activity, but its role in NK cells remains unclear.

PurposeThe impact of bioactive constituents of V. officinalis on immune responses remains largely elusive. In this work we investigated the potential targets of V. officinalis and focused on killing efficiency and related functions of NK cells regulated by bioactive constituents of V. officinalis.

Study design/methodsWe used primary human NK cells from peripheral blood mononuclear cells. Potential regulatory roles of selected compounds in NK killing function-related genes were first analyzed by network pharmacology approaches. Potential engagement between compounds and NK receptors was examined by molecular docking assay. Killing efficiency of NK cells was determined with real-time killing assay and live-cell imaging. Proliferation was examined by CFSE staining. Expression of cytotoxic proteins was analyzed using flow cytometry. Lytic granule release was quantified by CD107a degranulation assay. Contact time required for killing and determination of serial killers were analyzed using live cell imaging results.

ResultsResults from network pharmacology approaches and molecular docking analyses predicted that the five bioactive compounds (acteoside, apigenin, kaempferol, verbenalin and hastatoside) from V. officinalis can regulated various key genes essential for NK effector functions and can bind to NK inhibitory receptors as well as activating receptors, indicating their potential role in regulating NK killing capacity. Analyses of NK killing efficiency revealed that out of these five compounds verbenalin exhibited an encouraging potency to enhance NK killing efficiency. Further investigation showed that verbenalin did not affect proliferation, expression of cytotoxic proteins, or lytic granule degranulation, but rather reduced contact time required for killing, enhanced total killing events per NK cell, and elevated the fraction of killing competent NK cells.

ConclusionsIn this work, we have identified verbenalin as the most potent bioactive constituent from V. officinalis to enhance NK cell killing efficiency via accelerating killing processes and converting non-killers to killing competent NK cells without affecting the lytic granule pathway. Enhancing CDK4/6 inhibitor therapy for medulloblastoma using nanoparticle delivery and scRNA-seq-guided combination with sapanisertib CDK4/6 inhibitors hold promise for brain tumor treatment, but efficacy has been limited by recurrence in both preclinical models and clinical trials. To address recurrence, we tested a nanoparticle formulation of the CDK4/6 inhibitor palbociclib (POx-palbo) in mice genetically-engineered to develop SHH-driven medulloblastoma. We then analyzed medulloblastomas in mice receiving palbociclib treatment, and compared the efficacy of combining palbociclib with specific inhibitors suggested by our analysis. POx-Palbo showed reduced toxicity compared to conventional palbociclib, was tolerable in parenteral administration, improved CNS pharmacokinetics, and extended survival of mice with medulloblastoma. Recurrence, however, remained problematic as fractions of tumor cells proliferated during therapy. ScRNA-seq identified a gene expression pattern unique to proliferating medulloblastoma cells in POx-Palbo-treated mice, marked by up-regulation of the glutamate transporter Slc1a2 and down-regulation of diverse ribosomal genes. Reduced mTORC1 signaling, suggested by ribosomal suppression in POx-Palbo-treated tumors was confirmed by decreased 4EBP1 phosphorylation (p4EBP1). Further reducing mTORC1 activity by combining POx-Palbo with the mTORC1 inhibitor sapanisertib produced mutually enhancing effects, with increased suppression of both pRB and p4EBP1, and prolonged mouse survival compared to either agent alone. In contrast, targeting cell cycle progression by combining POx-Palbo with the SHH-pathway inhibitor vismodegib, or with the replication-targeting agents gemcitabine or etoposide, failed to enhance efficacy. Our data show the potential of nanoparticle formulation and scRNA-seq analysis of resistance to improve brain tumor treatment, and identify POx-palbo plus sapanisertib as effective combinatorial therapy for SHH medulloblastoma. This combination may be appropriate for testing in patients with recurrence, who need new options. Balancing sensitivity and specificity in preclinical research The success of scientific discovery in preclinical drug development is based on the different roles of exploration and confirmation. Exploration requires sensitivity while confirmation needs to be specific. Based on published empirical data, we simulated two preclinical research trajectories to examine their effectiveness to detect true effects while keeping animal numbers low. We compare a conventional p-value based approach with a method based on an a priori determined smallest effect size of interest (SESOI). Using a SESOI increases transition rates from exploration to confirmation and results in higher detection rates across the trajectory. Specificity in the SESOI trajectory increases if the prior probability of true effects is low. We conclude that employing a SESOI is superior to a p-value based approach. Based on our findings, we propose a reconsideration of planning and conducting preclinical experiments, especially when the prior probability of true hypotheses is low. Inter-microscope comparability of dental microwear texture data obtained from different optical profilometers. Dental microwear texture analysis has become a well-established approach for dietary inference and reconstruction of mammals and other tetrapods, both extant and extinct. As the amount of available data grows immensely, researchers could benefit from combining data gathered by others to perform meta-analyses. However, different devices used to capture three-dimensional surface scans for DMTA are known to produce variation even when measuring the same surface. Here we compare DMTA data of 36 guinea pigs that received different diets in a controlled feeding experiment, measured on a confocal disc-scanning and a confocal laser-scanning microscope. We are testing different pre-analysis filtering protocols to mitigate differences. We find inter-microscopes and filter-related differences for the majority of analysed 41 DMTA parameters. Certain microscope-specific filter routines resulted in less differences than other pre-analysis protocols. We further identify DMTA parameters which were stable regardless of microscope or pre-analysis treatment. Overall, the results obtained on both microscopes show the same dietary differentiation between guinea pig feeding groups, which supports that DMTA is a suitable method to obtain repeatable, objective dietary inferences. We finally propose a roadmap to enhance data exchange and inter-lab comparability and collaboration in the future. A humanized nanobody phage display library yields potent binders of SARS CoV-2 spike Neutralizing antibodies targeting the SARS-CoV-2 spike protein have shown a great preventative/therapeutic potential. Here, we report a rapid and efficient strategy for the development and design of SARS-CoV-2 neutralizing humanized nanobody constructs with sub-nanomolar affinities and nanomolar potencies. CryoEM-based structural analysis of the nanobodies in complex with spike revealed two distinct binding modes. The most potent nanobody, RBD-1-2G(NCATS-BL8125), tolerates the N501Y RBD mutation and remains capable of neutralizing the B.1.1.7 (Alpha) variant. Molecular dynamics simulations provide a structural basis for understanding the neutralization process of nanobodies exclusively focused on the spike-ACE2 interface with and without the N501Y mutation on RBD. A primary human airway air-lung interface (ALI) ex vivo model showed that RBD-1-2G-Fc antibody treatment was effective at reducing viral burden following WA1 and B.1.1.7 SARS-CoV-2 infections. Therefore, this presented strategy will serve as a tool to mitigate the threat of emerging SARS-CoV-2 variants. Climate change and chill accumulation: implications for tree fruit production in cold-winter regions Winter chill accumulation is critical for the productivity and profitability of perennial tree fruit systems. Several studies have quantified the impacts of global warming on chill accumulation in the warmer production regions of the world, where insufficient chill events occur and their frequency is increasing. In contrast, we focus on a region with relatively cold winters-the Pacific Northwest United States (PNW)-where insufficient chill events are currently absent, and quantify the potential for introduction of these risks under climate change. Results identified spatial variation within the PNW, with chill accumulation projected to increase in some areas but decrease in others. There was also spatiotemporal variation in the driving factors of changes to chill accumulation. Even with decreases in chill accumulations, there are likely minimal issues with insufficient chill accumulation. However, delayed chill accumulation in combination with advances in the onset of heat accumulation can potentially shift the region from one where spring phenology is primarily forcing-driven to one where the dynamic interplay between chilling and forcing processes become important. These interactions might create production risks for varieties with high chill requirements, post mid-21st-century under high emissions scenarios. Future work should focus on understanding, modeling, and projecting responses across these overlapping chilling and forcing processes. Additionally, given significant spatial differences across a relatively small geographic range, it is also critical to understand and model these dynamics at a local landscape resolution for regions such as the PNW. Accuracy and Power Analysis of Social Networks Built From Count Data O_LIPower analysis is used to estimate the probability of correctly rejecting a null hypothesis for a given statistical model and dataset. Conventional power analyses assume complete information, but the stochastic nature of behavioural sampling can mean that true and estimated networks are poorly correlated. Power analyses do not currently take the effect of sampling into account. This could lead to inaccurate estimates of statistical power, potentially yielding misleading results.
C_LIO_LIHere we develop a method for computing network correlation: the correlation between an estimated social network and its true network, using a Gamma-Poisson model of social event rates for networks constructed from count data. We use simulations to assess how the level of network correlation affects the power of nodal regression analyses. We also develop a generic method of power analysis applicable to any statistical test, based on the concept of diminishing returns.
C_LIO_LIWe demonstrate that our network correlation estimator is both accurate and moderately robust to its assumptions being broken. We show that social differentiation, mean social event rate, and the harmonic mean of sampling times positively impacts the strength of network correlation. We also show that the required level of network correlation to achieve a given power level depends on many factors, but that 0.80 network correlation usually corresponds to around 0.80 power for nodal regression in ideal circumstances.
C_LIO_LIWe provide guidelines for using our network correlation estimator to verify the accuracy of networks built from count data, and to conduct power analysis. This can be used prior to data collection, in post hoc analyses, or even for subsetting networks in dynamic network analysis. The network correlation estimator and custom power analysis methods have been made available as an R package.
C_LI Distributed brain co-processor for tracking electrophysiology and behavior during electrical brain stimulation Early implantable epilepsy therapy devices provided open-loop electrical stimulation without brain sensing, computing, or an interface for synchronized behavioral inputs from patients. Recent epilepsy stimulation devices provide brain sensing but have not yet developed analytics for accurately tracking and quantifying behavior and seizures. Here we describe a distributed brain co-processor providing an intuitive bi-directional interface between patient, implanted neural stimulation and sensing device, and local and distributed computing resources. Automated analysis of continuous streaming electrophysiology is synchronized with patient reports using a hand-held device and integrated with distributed cloud computing resources for quantifying seizures, interictal epileptiform spikes, and patient symptoms during therapeutic electrical brain stimulation. The classification algorithms for interictal epileptiform spikes and seizures were developed and parameterized using long-term ambulatory data from 9 humans and 8 canines with epilepsy, and then implemented prospectively in out-of-sample testing in 2 pet canines and 4 humans with drug resistant epilepsy living in their natural environments. Live imaging YAP signaling in mouse embryo development YAP protein is a critical regulator of mammalian embryonic development. By generating a near-infrared fusion YAP reporter mouse line, we have achieved high-resolution live imaging of YAP localization during mouse embryonic development. We have validated the reporter by demonstrating its predicted responses to blocking LATs kinase activity or blocking cell polarity. By time lapse imaging preimplantation embryos, we revealed a mitotic reset behavior of YAP nuclear localization. We also demonstrated deep tissue live imaging in post-implantation embryos and revealed an intriguing nuclear YAP pattern in migrating cells. The YAP fusion reporter mice and imaging methods will open new opportunities for understanding dynamic YAP signaling in vivo in many different situations. Factors behind intraspecific competition during post-glacial recolonization: a Norway spruce case study O_LITo a large extent, the current distribution and population structure of many species were shaped by the cycles of isolation in glacial refugia and subsequent population expansions. Isolation in, and post-glacial expansion through, heterogeneous environments led to either neutral or adaptive divergence. Norway spruce is no exception, and its current distribution is the consequence of a complex interplay between evolutionary forces.
C_LIO_LIWe investigated population differentiation and adaptations of Norway spruce for juvenile growth, diameter of the trunk at breast height and its components: wood density and tracheid traits. Data from 4,461 phenotyped and genotyped Norway spruce from 396 half-sib families in two progeny tests were used to test for divergent selection in the framework of QST vs FST, and to carry out multiple Genome-Wide Association Studies.
C_LIO_LIWe show that (i) the macroscopic resultant trait (trunk diameter), unlike its microscopic components (tracheid dimensions) and juvenile growth, was under divergent selection, and (ii) adaptation predated the Last Glaciation Maximum.
C_LIO_LIAltogether, the current variation in these phenotypic traits in Norway spruce is better explained by local adaptation to ancestral environments than to current ones, where populations were partly pre-adapted mainly through growth-related traits.
C_LI An exon-intron split framework to prioritize miRNA-driven post-transcriptional regulatory signals and its application to study energy homeostasis in pigs The contribution of microRNAs (miRNAs) to mRNA regulation has often been explored by post hoc selection of downregulated genes and determining whether they harbor binding sites for miRNAs of interest. This approach, however, does not discriminate whether these mRNAs are also downregulated at the transcriptional level. Here, we have characterized the transcriptional and post-transcriptional changes of mRNA expression in two porcine tissues: gluteus medius muscle of fasted and fed Duroc gilts and adipose tissue of lean and obese Duroc-Gottingen minipigs. Exon-intron split analysis (EISA) of RNA-seq data allowed us to identify downregulated mRNAs with high post-transcriptional signals in fed or obese states, and we assessed whether they harbor binding sites for upregulated miRNAs in any of these two physiological states. We found 26 downregulated mRNAs with high post-transcriptional signals in the muscle of fed gilts and 21 of these were predicted targets of upregulated miRNAs also in the fed state. For adipose tissue, 44 downregulated mRNAs in obese minipigs displayed high post-transcriptional signals, and 25 of these were predicted targets of miRNAs upregulated in the obese state. These results suggest that the contribution of miRNAs to mRNA repression is more prominent in the skeletal muscle system. Finally, we identified several genes that may play relevant roles in the energy homeostasis of the pig skeletal muscle (DKK2 and PDK4) and adipose (SESN3 and ESRRG) tissues. By differentiating transcriptional from post-transcriptional changes in mRNA expression, EISA provides a valuable view about the regulation of gene expression, complementary to canonical differential expression analyses. Critical re-evaluation of the slope factor of the operational model of agonism Background and PurposeAlthough being a relative term, agonist efficacy is a cornerstone in the proper assessment of agonist selectivity and signalling bias. The operational model of agonism (OMA) has become successful in the determination of agonist efficacies and ranking them. In 1985, Black et al. introduced the slope factor to the OMA to make it more flexible and allow for fitting steep as well as flat concentration-response curves.

Experimental ApproachFunctional analysis of OMA demonstrates that the slope factor implemented by Black et al. affects relationships among parameters of the OMA, leading to wrong estimates of operational efficacy and affinity of an agonist. In contrast, functional analysis of OMA with Hill coefficient as a slope factor has not this flaw.

Key ResultsFitting of the OMA with Black et al. slope factor to concentration-response curves of several models (homotropic allosteric modulation, substrate inhibition and non-competitive auto-inhibition) gives wrong estimates of operational efficacy and affinity. In contrast, fitting of the OMA modified by the Hill coefficient to concentration-response curves of these several models gives correct estimates of operational efficacy and affinity.

Conclusion and ImplicationsThe OMA modified by the Hill coefficient can be used for correct ranking of agonism and subsequent analysis, like quantification of signalling bias. Stress alters the activity of Leydig cells dependently on the diurnal time BackgroundThe increasing amount of data points to the circadian timing system as an essential part of processes regulating androgen homeostasis. However, the relationship between stress response, timekeeping-, and steroidogenesis-related systems is unexplored.

ObjectiveThe purpose of the study was to analyze the stress-response of the testosterone-producing Leydig cells depending on the stressful events time.

Materials and methodsThe study was designed to follow the effects of 3-hour immobilization (IMO) applied at different periods during the day. The IMO performed once (1xIMO) or repeated in 10 consecutive days (10xIMO). Principal-component-analysis (PCA) followed the expression study of the clock and steroidogenic-related genes in Leydig cells.

ResultsBoth types of IMO in all investigated periods increased corticosterone and decreased testosterone blood level. Transcriptional analysis revealed different sensitivity to IMO events depending on the circadian time. The majority of steroidogenesis-related genes (Lhcgr, Cyp11a1, Cyp17a1, Hsd3b1/2) were down-regulated in the inactive but unchanged or even up-regulated in the active phase of the day. Both types of IMO potentiated the expression of clock elements Bmal1/BMAL1, Per1/PER1 regardless of the days stage and reduced Reverba in the inactive phase. The PCA confirmed a major shift, for both IMO-types, in the transcription of steroidogenesis and clock genes across the inactive/active phase. Further, the diurnal pattern of the glucocorticoid receptor (Nr3c1/GR) expression in Leydig cells was increased in the inactive phase due to 10xIMO. The observed time-dependent IMO-response of the Leydig cells correlated with different GR engagements.

DiscussionStress- and the circadian-system coordinatively shape Leydig cells physiology, assuming diverse GR engagement as a possible factor in mediating the diurnal effect of stress.

ConclusionThe Leydig cells stress-response depends on the time of the stressful situation, emphasizing the importance of circadian activity in supporting androgen homeostasis and male fertility. wnt16 exerts pleiotropic effects on bone and lean mass in zebrafish Osteoporosis, a disease of bone fragility, and sarcopenia, a condition of reduced muscle mass and strength, commonly occur in the same individual--a condition termed osteosarcopenia. Genetic variants at the CPED1-WNT16 locus have pleiotropic effects on human bone mineral density (BMD) and lean tissue mass. However, our current understanding of Wnt16 and its influence on bone fails to explain effects on lean mass at this locus. In this study, we sought to better understand the biological mechanism underlying pleiotropy at the CPED1-WNT16 locus on BMD and lean mass. Using single cell analysis, microCT imaging, and genetic approaches in zebrafish, we show that wnt16 exerts dual influence on muscle and bone, and provide evidence that this occurs through a function in somite development, a period in which precursors for both of these tissues are specified. We show that mutations in wnt16 alter vertebral size and lean mass during axial skeletogenesis. We show that wnt16+ cells transcriptionally resemble cells within the amniote dermomyotome, and provide evidence they function as myogenic precursors during embryonic development. We also provide evidence that WNT16 is a gene of major effect on lean mass at the CPED-WNT16 locus. Our studies support the potential for wnt16 to exert pleiotropic effects on bone and lean mass through a function in myogenic precursors. Nuclear NAD+-biosynthetic enzyme NMNAT1 facilitates survival of developing retinal neurons Despite mounting evidence that the mammalian retina is exceptionally reliant on proper NAD+ homeostasis for health and function, the specific roles of subcellular NAD+ pools in retinal development, maintenance, and disease remain obscure. Here, we show that deletion of the nuclear-localized NAD+ synthase nicotinamide mononucleotide adenylyltransferase-1 (NMNAT1) in the developing murine retina causes early and severe degeneration of photoreceptors and select inner retinal neurons via multiple distinct cell death pathways. This severe phenotype is associated with disruptions to retinal central carbon metabolism, purine nucleotide synthesis, and amino acid pathways. Furthermore, large-scale transcriptomics reveals dysregulation of a collection of photoreceptor and synapse-specific genes in NMNAT1 knockout retinas prior to detectable morphological or metabolic alterations. Collectively, our study reveals previously unrecognized complexity in NMNAT1-associated retinal degeneration and suggests a yet-undescribed role for NMNAT1 in gene regulation during photoreceptor terminal differentiation. Antithrombin glycoforms are selectively physiosorbed on plasma extracellular vesicles Antithrombin (AT) is a glycoprotein produced by the liver and a principal antagonist of active clotting proteases. A deficit in AT function leads to AT qualitative deficiency, challenging to diagnose. Here we report that active AT may travel physiosorbed on the surface of plasma extracellular vesicles (EVs), contributing to form the "EV-protein corona". The corona is enriched in specific AT glycoforms, thus suggesting glycosylation to play a key role in AT partitioning between EVs and plasma. Differences in AT glycoform composition of the corona of EVs separated from plasma of healthy and AT qualitative deficiency-affected subjects were also noticed. This suggests deconstructing the plasma into its nanostructured components, as extracellular vesicles, could help to unravel pathophysiological mechanisms otherwise undiscovered. Re-convolving the compositional landscape of primary and recurrent glioblastoma reveals prognostic andtargetable tissue states Glioblastoma is an aggressive diffusely infiltrating neoplasm that spreads beyond surgical resection margins, where it intermingles with non-neoplastic brain cells. This complex microenvironment harbouring infiltrating glioma and non-neoplastic brain cells is the origin of tumor recurrence. Thus, understanding the cellular and molecular features of the glioma microenvironment is therapeutically and prognostically important. Here, we used single-nucleus RNA sequencing (snRNAseq) of primary and recurrent glioma to define three compositional  tissue-states rooted in cohabitation of cell-types and transcriptional states. These comprise a state featuring A) abundant normal brain microenvironment cells, B) reactive/inflamed brain microenvironment, and C) cellular/proliferative tumor. All these states exhibited variable degrees of infiltration by glioma cells. Spatial transcriptomics confirmed that the cell-types and transcriptomics states which compositionally cohabitate tissue states also colocalize in space. Tissue states are clinically significant because they correlate with radiographic, histopathologic, and prognostic features. We found that enrichment of tissue state B signature correlated with shorter survival. Importantly, we found that our compositionally defined tissue states are enriched in distinct metabolic pathways. One such pathway is fatty acid biosynthesis, which was enriched in tissue state B - a state enriched in recurrent glioblastoma and composed of astrocyte-like/mesenchymal glioma cells, reactive astrocytes resembling those seen in neurodegeneration, and monocyte-like myeloid cells. We showed that treating acute slices of GBM with a fatty acid synthesis inhibitor is sufficient to deplete the transcriptional signature of tissue state B. Our findings define a novel compositional approach to the glioma-infiltrated tissue which allows us to discover prognostic and targetable features, paving the way to new mechanistic and therapeutic discoveries. P2X2 receptor subunit interfaces are missense variant hotspots where mutations tend to increase apparent ATP affinity Background and PurposeP2X receptors (P2XRs) are trimeric ligand-gated ion channels (LGICs) that open a cation-selective pore in response to ATP binding to their large extracellular domain (ECD). The seven known P2XR subtypes can assemble as homo- or heterotrimeric complexes and contribute to numerous physiological functions, including nociception, inflammation and hearing. The overall structure of P2XRs is well established, but little is known about the spectrum and prevalence of human genetic variations and the functional implications in specific domains.

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

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

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

Bullet point summary What is already known: O_LIP2X2 receptors are ATP-activated ion channels implicated in hearing and nociceptice pathways
C_LI
 What this study adds: O_LIA structural mapping of missense variants observed in the human population
C_LIO_LIWe identify the intersubunit-interface as a variant hotspot and decipher functional impact of mutations
C_LI
 Clinical significance: O_LIThe development of both inhibitors and activators of P2X2 receptor function may be required
C_LI 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. Differential regulation of degradation and immune pathways underlies adaptation of the ectosymbiotic nematode Laxus oneistus to oxic-anoxic interfaces Eukaryotes may experience oxygen deprivation under both physiological and pathological conditions. Because oxygen shortage leads to a reduction in cellular energy production, all eukaryotes studied so far conserve energy by suppressing their metabolism. However, the molecular physiology of animals that naturally and repeatedly experience anoxia is underexplored. One such animal is the marine nematode Laxus oneistus. It thrives, invariably coated by its sulfur-oxidizing symbiont Candidatus Thiosymbion oneisti, in anoxic sulfidic or hypoxic sand. Here, transcriptomics and proteomics showed that, whether in anoxia or not, L. oneistus mostly expressed genes involved in ubiquitination, energy generation, oxidative stress response, immune response, development, and translation. Importantly, ubiquitination genes were also upregulated when the nematode was subjected to anoxic sulfidic conditions, together with genes involved in autophagy, detoxification and ribosome biogenesis. We hypothesize that these degradation pathways were induced to recycle damaged cellular components (mitochondria) and misfolded proteins into nutrients. Remarkably, when L. oneistus was subjected to anoxic sulfidic conditions, lectin and mucin genes were also upregulated, potentially to promote the attachment of its thiotrophic symbiont. Furthermore, the nematode appeared to survive oxygen deprivation by using an alternative electron carrier (rhodoquinone) and acceptor (fumarate), to rewire the electron transfer chain. On the other hand, under hypoxia, genes involved in costly processes (e.g., amino acid biosynthesis, development, feeding, mating) were upregulated, together with the worms Toll- like innate immunity pathway and several immune effectors (e.g., Bacterial Permeability Increasing proteins, fungicides).

In conclusion, we hypothesize that, in anoxic sulfidic sand, L. oneistus upregulates degradation processes, rewires oxidative phosphorylation and by reinforces its coat of bacterial sulfur-oxidizers. In upper sand layers, instead, it appears to produce broad-range antimicrobials and to exploit oxygen for biosynthesis and development. HTSplotter: an end-to-end data processing, analysis and visualisation tool for chemical and genetic in vitro perturbation screening Biomedical researchers are moving towards high-throughput screening, as this allows for automatization, better reproducibility and more and faster results. High-throughput screening experiments encompass drug, drug combination, genetic perturbagen or a combination of genetic and chemical perturbagen screens. These experiments are conducted in real-time assays over time or in an endpoint assay. The data analysis consists of data cleaning and structuring, as well as further data processing and visualisation, which, due to the amount of data, can easily become laborious, time consuming, and error-prone. Therefore, several tools have been developed to aid researchers in this data analysis, but they focus on specific experimental set-ups and are unable to process data of several time points and genetic-chemical perturbagen screens together. To meet these needs, we developed HTSplotter, available as web tool and Python module, that performs automatic data analysis and visualisation of either endpoint or real-time assays from different high-throughput screening experiments: drug, drug combination, genetic perturbagen and genetic-chemical perturbagen screens. HTSplotter implements an algorithm based on conditional statements in order to identify experiment type and controls. After appropriate data normalization, HTSplotter executes downstream analyses such as dose-response relationship and drug synergism by the Bliss independence method. All results are exported as a text file and plots are saved in a PDF file. The main advantage of HTSplotter over other available tools is the automatic analysis of genetic-chemical perturbagen screens and real-time assays where results are plotted over time. In conclusion, HTSplotter allows for the automatic end-to-end data processing, analysis and visualisation of various high-throughput in vitro cell culture screens, offering major improvements in terms of versatility, convenience and time over existing tools.

Author summaryAcademic researchers are moving towards high-throughput screenings, where the experiments execution follows an automatic approach, such as robotic seeding, liquid dispensing and/or automatics readouts. This grants more flexible and reproducible experimental set ups. The type of high-throughput experiment can vary from drug, drug combination, genetic perturbagen to genetic-chemical perturbagen screens. These can be assessed through endpoint assays, measuring the effect at one time point, or real-time assays, assessing the phenotypic effect over time. High-throughput screening results in large amounts of data, requiring laborious, time consuming and error-prone data handling and analysis. These challenges hamper the biological interpretation of data and the fast progress in drug development programs. Hence, we developed HTSplotter, a web tool and Python module, to allow researchers to conduct a fast and flexible end-to-end data processing, analysis and visualisation of different high-throughput experiments, assessed either by endpoint or real-time assays, relieving the high-throughput analysis bottleneck. In this way, HTSplotter, directly contributes to a faster identification of drugs and drug combinations, with potential practical applications, or novel targets for therapies. 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. Sequence-dependent model of genes with dual σ factor preference Escherichia coli uses {sigma} factors to quickly control large gene cohorts during stress conditions. While most of its genes respond to a single {sigma} factor, approximately 5% of them have dual {sigma} factor preference. The most common are those responsive to both {sigma}70, which controls housekeeping genes, and {sigma}38, which activates genes during stationary growth and stresses. Using RNA-seq and flow-cytometry measurements, we show that  {sigma}70+38 genes are nearly as upregulated in stationary growth as  {sigma}38 genes. Moreover, we find a clear quantitative relationship between their promoter sequence and their response strength to changes in {sigma}38 levels. We then propose and validate a sequence dependent model of {sigma}70+38 genes, with dual sensitivity to {sigma}38 and {sigma}70, that is applicable in the exponential and stationary growth phases, as well in the transient period in between. We further propose a general model, applicable to other stresses and {sigma} factor combinations. Given this, promoters controlling {sigma}70+38 genes (and variants) could become important building blocks of synthetic circuits with predictable, sequence-dependent sensitivity to transitions between the exponential and stationary growth phases. Context coding in the mouse nucleus accumbens modulates motivationally relevant information Neural activity in the nucleus accumbens (NAc) is thought to track fundamentally value-centric quantities such as current or future expected reward, reward prediction errors, the value of work, opportunity cost, and approach vigor. However, the NAc also contributes to flexible behavior in ways that are difficult to explain based on value signals alone, raising the question of if and how non-value signals are encoded in NAc. We recorded NAc neural ensembles while head-fixed mice performed a biconditional discrimination task where context-setting cues modulated the stimulus-outcome association of subsequently presented reward-predictive cues. We extracted single-unit and population-level correlates of task features, and found value-independent coding for the context cues. This context signal occupied a subspace orthogonal to outcome-predictive representations, and was predictive of subsequent behaviorally-relevant coding. Together, these findings support a circuit-level gating model for how the NAc contributes to behavioral flexibility and provide a novel population-level perspective from which to view NAc computations.