Integrative Modelling of Signalling Network Dynamics Identifies Cell Type-selective Therapeutic Strategies for FGFR4-driven Cancers Oncogenic FGFR4 signalling represents a potential therapeutic target in many cancer types, including triple negative breast cancer (TNBC) and hepatocellular carcinoma (HCC). However, resistance to single-agent therapy directed at FGFR4 remains a major challenge, prompting the need to identify more effective combinatorial therapeutic strategies. Here, we integrated computational network modelling and experimental validation to characterise dynamic reprogramming of the FGFR4 signalling network in TNBC following FGFR4 kinase inhibition. We found that AKT, which signals downstream of FGFR4, displayed a rapid and potent reactivation following FGFR4 targeting. Through model-based simulation and systematic prediction of the effect of co-targeting specific network nodes, we predicted, and validated experimentally, strong synergism of co-targeting FGFR4 and particular ErbB kinases or AKT, but not the upstream kinase PI3K. Further, incorporation of protein expression data from hundreds of cancer cell types enabled us to adapt our model to other diverse cellular contexts, leading to the prediction that while AKT rebound occurs frequently, it is not a general phenomenon. Instead, ERK is reactivated in a subset of cell types, including the FGFR4-driven HCC cell line Hep3B. This was subsequently corroborated, and moreover, co-targeting FGFR4 and MEK in Hep3B cells markedly enhanced inhibition of cell proliferation. Overall, these findings provide novel insights into the dynamics of drug-induced network remodelling in cancer cells, highlight the impact of protein expression heterogeneity on network response to targeted therapy and identify candidate cell type-selective combination treatments for FGFR4-driven cancer.
Intron-mediated induction of phenotypic heterogeneity Introns are universally present in the nuclear genomes of eukaryotes1. The budding yeast, an otherwise intron-poor species, preserves two sets of ribosomal protein (RP) genes differing primarily in their introns2-4. Despite recent findings on the role of RP introns under stress and starvation5-7, understanding the contribution of introns to ribosome regulation remains challenging. Here, combining isogrowth profiling8 with single-cell protein measurements9, we found that introns can mediate inducible phenotypic heterogeneity conferring a clear fitness advantage. Osmotic stress leads to bimodal expression of the small ribosomal subunit protein Rps22B, mediated by an intron in the 5 untranslated region of its transcript. The two resulting yeast subpopulations differ in their ability to cope with starvation. Low Rps22B protein levels resulted in prolonged survival under sustained starvation, while high Rps22B levels enabled cells to grow faster after transient starvation. Further, yeast growing at high sugar concentrations - similar to those in ripe grapes - exhibit bimodal Rps22B expression when approaching stationary phase. Differential intron-mediated regulation of RP genes thus provides a way to diversify the population when starvation looms in natural environments. Our findings reveal a new role for introns in inducing phenotypic heterogeneity in changing environments and suggest that duplicated RP genes in yeast contribute to resolving the evolutionary conflict between precise expression control and environmental responsiveness10.
The DC1 domain protein BINUCLEATE POLLEN is required for pollen development in Arabidopsis thaliana Development of the male gametophyte is a tightly regulated process that requires precise control of cell division and gene expression. A relevant aspect to understand the events underlying pollen development regulation constitutes the identification and characterization of the genes required for this process. In this work we showed that the DC1 domain protein BINUCLEATE POLLEN (BNP) is essential for pollen development and germination. Pollen grains carrying the defective BNP allele failed to complete mitosis II and are impaired in pollen germination. By yeast two hybrid screening and bimolecular fluorescence complementation assays, we identified a set of BNP-interacting proteins. Among confirmed interactors we found NAC family transcriptional regulators Vascular plant One-Zinc finger 1 (VOZ1) and VOZ2. VOZ1 localization changes during pollen development, moving to the vegetative nucleus at the tricellular stage. We observed that this relocalization requires BNP, as in their absence in pollen from bnp/BNP plants, VOZ1 nuclear localization resulted impaired. As voz1voz2 double mutants showed the same developmental defect observed in bnp pollen grains, we propose that BNP requirement to complete microgametogenesis could be linked to its interaction with VOZ1/2 proteins. BNP could have a role of scaffold protein, recruiting VOZ1/2 to the endosomal system into assemblies that are required for their further translocation to the nucleus, where they act as transcriptional regulators. 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. Decellularized Normal and Tumor Extracellular Matrix as Scaffold for Cancer Organoid Cultures of Colorectal Peritoneal Metastases Peritoneal metastases (PM) from colorectal cancer (CRC) are associated with poor survival. The extracellular matrix (ECM) plays a fundamental role in modulating the homing of CRC metastases to the peritoneum. The mechanisms underlying the interactions between metastatic cells and the ECM, however, remain poorly understood and the number of in vitro models available for the study of the peritoneal metastatic process is limited. Here, we show that decellularized ECM of the peritoneal cavity allows the growth of organoids obtained from PM, favoring the development of three-dimensional nodules that maintain the characteristics of in vivo PM. Organoids preferentially grow on scaffolds obtained from neoplastic peritoneum, which are characterized by greater stiffness than normal scaffolds. A gene expression analysis of organoids grown on different substrates reflected faithfully the clinical and biological characteristics of the organoids. An impact of the ECM on the response to standard chemotherapy treatment for PM was also observed.

SignificanceEvidence of the value of ex vivo 3D models obtained by combining patient-derived extracellular matrices depleted of cellular components and organoids to mimic the metastatic niche, to be used as a tool to develop new therapeutic strategies in a biologically relevant context, to personalize treatments and increase their efficacy. Co-option of immune effectors by the hormonal signalling system triggering metamorphosis in Drosophila melanogaster Insect metamorphosis is regulated by the production, secretion and degradation of two peripheral hormones: 20-hydroxyecdysone (ecdysone) and juvenile hormone (JH). In addition to their roles in developmental regulation, increasing evidence suggests that these hormones are involved in innate immunity processes, such as phagocytosis and the induction of antimicrobial peptide (AMP) production. AMP regulation includes systemic responses as well as local responses at surface epithelia that contact with the external environment. At pupariation, Drosophila melanogaster increases dramatically the expression of three AMP genes, drosomycin (drs), drosomycin-like 2 (drsl2) and drosomycin-like 5 (drsl5). We show that the systemic action of drs at pupariation is dependent on ecdysone signalling in the fat body and operates via the ecdysone downstream target, Broad-Z4. In parallel, ecdysone also regulates local responses, specifically through the activation of drsl2 expression in the gut. Finally, we confirm the relevance of this ecdysone dependent AMP expression for the control of bacterial persistence by showing that flies lacking drs expression in the fat body have higher bacterial persistence over metamorphosis. Together, our data emphasize the importance of the association between ecdysone signalling and immunity using in vivo studies and establish a new role for ecdysone at pupariation, which impacts developmental success by regulating the immune system in a stage-dependent manner. We speculate that this co-option of immune effectors by the hormonal system may constitute a pre-emptive mechanism to control bacterial numbers in the pupa at the core of metamorphosis evolution. Inhibitors of ERp44, PDIA1, and AGR2 induce disulfide-mediated oligomerization of Death Receptors 4 and 5 and cancer cell death Breast cancer mortality remains unacceptably high, indicating a need for safer and more effective therapeutic agents. Disulfide bond Disrupting Agents (DDAs) were previously identified as a novel class of anticancer compounds that selectively kill cancers that overexpress the Epidermal Growth Factor Receptor (EGFR) or its family member HER2. DDAs kill EGFR+ and HER2+ cancer cells via the parallel downregulation of EGFR, HER2, and HER3 and activation/oligomerization of Death Receptors 4 and 5 (DR4/5). However, the mechanisms by which DDAs mediate these effects are unknown. Affinity purification analyses employing biotinylated-DDAs reveal that the Protein Disulfide Isomerase (PDI) family members AGR2, PDIA1, and ERp44 are DDA target proteins. Further analyses demonstrate that shRNA-mediated knockdown of AGR2 and ERp44, or expression of ERp44 mutants, enhance basal and DDA-induced DR5 oligomerization. DDA treatment of breast cancer cells disrupts PDIA1 and ERp44 mixed disulfide bonds with their client proteins. Together, the results herein reveal DDAs as the first small molecule, active site inhibitors of AGR2 and ERp44, and demonstrate roles for AGR2 and ERp44 in regulating the activity, stability, and localization of DR4 and DR5, and activation of Caspase 8. Spatial regulation of the anaphase-promoting complex/cyclosome by a microtubule adaptor The anaphase-promoting complex/cyclosome (APC/C) coordinates advancement through mitosis via temporally controlled polyubiquitination of effector proteins. Despite the long-appreciated spatial organization of key events in mitosis mediated largely by cytoskeletal networks, the spatial regulation of APC/C, the major mitotic E3 ligase, is poorly understood. Here, we describe a microtubule-resident protein, PLEKHA5, as an interactor of APC/C and spatial regulator of its activity in mitosis. PLEKHA5 knockdown delayed mitotic progression, causing accumulation of APC/C substrates dependent upon the PLEKHA5-APC/C interaction. To assess the spatial importance of this interaction, we developed microtubule-localized proximity biotinylation tools, which revealed that depletion of PLEKHA5 decreased the extent of APC/C association with the microtubule cytoskeleton. This decreased APC/C microtubule-localization in turn prevented efficient loading of APC/C with its co-activator CDC20, leading to defects in E3 ligase catalytic activity. We propose that PLEKHA5 functions as an adaptor of APC/C that promotes its subcellular localization to microtubules and facilitates its activation by CDC20, thus ensuring the timely turnover of key mitotic APC/C substrates and proper progression through mitosis. To knot or not to knot: Multiple conformations of the SARS-CoV-2 frameshifting RNA element The SARS-CoV-2 frameshifting RNA element (FSE) is an excellent target for therapeutic intervention against Covid-19. This small gene element employs a shifting mechanism to pause and backtrack the ribosome during translation between Open Reading Frames 1a and 1b, which code for viral polyproteins. Any interference with this process has profound effect on viral replication and propagation. Pinpointing the structures adapted by the FSE and associated structural transformations involved in frameshifting has been a challenge. Using our graph-theory-based modeling tools for representing RNA secondary structures, "RAG" (RNA-As-Graphs), and chemical structure probing experiments, we show that the 3-stem H-type pseudoknot (3_6 dual graph), long assumed to be the dominant structure has a viable alternative, an HL-type 3-stem pseudoknot (3_3) for longer constructs. In addition, an unknotted 3-way junction RNA (3_5) emerges as a minor conformation. These three conformations share Stems 1 and 3, while the different Stem 2 may be involved in a conformational switch and possibly associations with the ribo-some during translation. For full-length genomes, a stem-loop motif (2_2) may compete with these forms. These structural and mechanistic insights advance our understanding of the SARS-CoV-2 frameshifting process and concomitant virus life cycle, and point to three avenues of therapeutic intervention. Detecting phylodiversity-dependent diviersification with a general phylogenetic inference framework AO_SCPLOWBSTRACTC_SCPLOWDiversity-dependent diversification models have been extensively used to study the effect of ecological limits and feedback of community structure on species diversification processes, such as speciation and extinction. Current diversity-dependent diversification models characterise ecological limits by carrying capacities for species richness. Such ecological limits have been justified by niche filling arguments: as species diversity increases, the number of available niches for diversification decreases.

However, as species diversify they may diverge from one another phenotypically, which may open new niches for new species. Alternatively, this phenotypic divergence may not affect the species diversification process or even inhibit further diversification. Hence, it seems natural to explore the consequences of phylogenetic diversity-dependent (or phylodiversity-dependent) diversification. Current likelihood methods for estimating diversity-dependent diversification parameters cannot be used for this, as phylodiversity is continuously changing as time progresses and species form and become extinct.

Here, we present a new method based on Monte Carlo Expectation-Maximization (MCEM), designed to perform statistical inference on a general class of species diversification models and implemented in the R package emphasis. We use the method to fit phylodiversity-dependent diversification models to 14 phylogenies, and compare the results to the fit of a richness-dependent diversification model. We find that in a number of phylogenies, phylogenetic divergence indeed spurs speciation even though species richness reduces it. Not only do we thus shine a new light on diversity-dependent diversification, we also argue that our inference framework can handle a large class of diversification models for which currently no inference method exists. Dynamic tracing using ultra-bright labelling and multi-photon microscopy identifies endothelial uptake of poloxamer-188 coated poly(lactic-co-glycolic acid)nano-carriers in vivo Poly(lactic-co-glycolic acid) (PLGA)-based drug formulations are approved for the use in humans, however, the potential of PLGA to design nanoparticles (NPs) and target the central nervous system remains to be exploited.

The aim of the current study was design PLGA NPs, loading them with bulky fluorophores thereby increasing single particle fluorescence to a level visible by in vivo microscopy, and investigate their brain biodistribution. We developed, highly fluorescent 70 nm PLGA NPs significantly brighter then quantum dots enabling their visualization by intravital real-time 2-photon microscopy. We found that PLGA NPs coated with pluronic F-68 (PF-68) had a substantially longer plasma half-life than uncoated NPs and were taken up by cerebro-vascular endothelial cells. High resolution confocal microscopy revealed that coated PLGA NPs were present in late endothelial endosomes of cerebral vessels within 1 hour after systemic injection and were more readily taken up by endothelial cells in peripheral organs.

The current data suggest that PF-68 coated PLGA NPs are taken up by mouse cerebral and peripheral endothelial cells in vivo. The combination of ultra-bright NPs and in vivo imaging may thus represent a promising approach to reduce the gap between development and clinical application of nanoparticle-based drug carriers. Alpha variant versus D614G strain in the Syrian hamster model Late 2020, SARS-CoV-2 Alpha variant from lineage B.1.1.7 emerged in United Kingdom and gradually replaced the G614 strains initially involved in the global spread of the pandemic. In this study, we used a Syrian hamster model to compare a clinical strain of Alpha variant with an ancestral G614 strain. The Alpha variant succeeded to infect animals and to induce a pathology that mimics COVID-19. However, both strains replicated to almost the same level and induced a comparable disease and immune response. A slight fitness advantage was noted for the G614 strain during competition and transmission experiments. These data do not corroborate the epidemiological situation observed during the first half of 2021 in humans nor reports that showed a more rapid replication of Alpha variant in human reconstituted bronchial epithelium. RESTORING GLUTAMATE RECEPTOR SIGNALING IN PANCREATIC ALPHA CELLS RESCUES GLUCAGON RESPONSES IN TYPE 1 DIABETES Glucagon secretion from pancreatic alpha cells is crucial to prevent hypoglycemia. For reasons still unknown, people with type 1 diabetes lose this glucoregulatory mechanism and are susceptible to dangerous hypoglycemia. Here we show that alpha cells in living pancreas slices from donors with type 1 diabetes failed to secrete glucagon in response to decreases in glucose concentration, thus mirroring the in vivo unresponsiveness to hypoglycemia. Glucagon content and responses to KCl depolarization were not affected, suggesting that alpha cells retained their secretory potential. By contrast, alpha cells had severely impaired signaling via glutamate receptors of the AMPA/kainate type. Under healthy conditions, activating these receptors was required to elicit full glucagon responses to decreases in glucose levels. In type 1 diabetes, reactivating residual glutamate receptor function with the positive allosteric modulators cyclothiazide and aniracetam restored glucagon secretion in response to hypoglycemia. These positive allosteric modulators are already approved to treat other conditions and could be repurposed to prevent hypoglycemia and improve management of diabetes. Including random effects in statistical models in ecology: fewer than five levels? As generalized linear mixed-effects models (GLMMs) have become a widespread tool in ecology, the need to guide the use of such tools is increasingly important. One common guideline is that one needs at least five levels of a random effect. Having such few levels makes the estimation of the variance of random effects terms (such as ecological sites, individuals, or populations) difficult, but it need not muddy ones ability to estimate fixed effects terms - which are often of primary interest in ecology. Here, I simulate ecological datasets and fit simple models and show that having too few random effects terms does not influence the parameter estimates or uncertainty around those estimates for fixed effects terms. Thus, it should be acceptable to use fewer levels of random effects if one is not interested in making inference about the random effects terms (i.e. they are  nuisance parameters used to group non-independent data). I also use simulations to assess the potential for pseudoreplication in (generalized) linear models (LMs), when random effects are explicitly ignored and find that LMs do not show increased type-I errors compared to their mixed-effects model counterparts. Instead, LM uncertainty (and p values) appears to be more conservative in an analysis with a real ecological dataset presented here. These results challenge the view that it is never appropriate to model random effects terms with fewer than five levels - specifically when inference is not being made for the random effects, but suggest that in simple cases LMs might be robust to ignored random effects terms. Given the widespread accessibility of GLMMs in ecology and evolution, future simulation studies and further assessments of these statistical methods are necessary to understand the consequences of both violating and blindly following simple guidelines. Barbary ground squirrels do not have a sentinel system but instead synchronise vigilance Coordinated behavior, such as hunting in lions and coordinated vigilance as antipredator behavior, are examples of benefits of group-living. Instead of asynchronous vigilance, some social species synchronize their vigilance bouts or take turns acting as sentinels. To increase our knowledge on the evolution of vigilance behavior, we studied whether vigilance is coordinated in Barbary ground squirrels, Atlantoxerus getulus. We show that vigilance was synchronized instead of taking turns. Multiple non-mutually exclusive hypotheses could explain synchronization: Barbary ground squirrels may perch because 1) neighbors are perched (copying effect), 2) perch synchrony may be an emergent property of the ecology as all squirrels may be satiated at the same time (collective behavior), or 3) the benefits are large in terms of evading ambush predators and scanning effectiveness (watch each others back). Particularly, in habitats where the field of view is obstructed by man-made structures and multiple individuals may be necessary to watch for terrestrial predators, synchronized vigilance may have greater fitness benefits than sentinel behavior. We conclude that it is essential to test assumptions of coordination and, thus, to analyze coordination to describe sentinel systems.

Significance StatementVigilance behavior can be vital to an animals survival. Taking turns acting as sentinels or synchronizing vigilance bouts reduces the cost of the trade-off between feeding and predation risk. A sentinel system assumes that sentinels are vigilant from raised positions, warn group members of danger, and alternate vigilance bouts. However, the assumption of alternating vigilance bouts remains poorly tested. We tested this assumption in invasive Barbary ground squirrels. We found that instead of alternating, individuals synchronized their vigilance bouts. Perch synchrony may be 1) a response to perching group members (copying effect), 2) an emergent property of the species ecology, and 3) an adaptation to anthropogenically altered habitats (watch each others back). Increased PHOSPHO1 expression mediates cortical bone mineral density in renal osteodystrophy Patients with advanced chronic kidney disease (CKD) often present with skeletal abnormalities; a condition known as renal osteodystrophy (ROD). While Tissue-nonspecific alkaline phosphatase (TNAP) and PHOSPHO1 are recognized to be critical for bone mineralization, their role in the etiology of ROD is unclear. To address this, ROD was induced in both wild-type and Phospho1 knockout (P1KO) mice using dietary adenine supplementation. The mice presented with hyperphosphatemia, hyperparathyroidism, and elevated levels of FGF23 and bone turnover markers. In particular, we noted that in CKD mice, bone mineral density (BMD) was increased in cortical bone (p < 0.05) but decreased in trabecular bone (p < 0.05). These changes were accompanied by decreased TNAP (p < 0.01) and increased PHOSPHO1 (p < 0.001) expression in wild-type CKD bones. In P1KO CKD mice, the cortical BMD phenotype was rescued, suggesting that the increased cortical BMD of CKD mice was driven by increased PHOSPHO1 expression. Other structural parameters were also improved in P1KO CKD mice. We further investigated the driver of the mineralization defects, by studying the effects of FGF23, PTH, and phosphate administration on PHOSPHO1 and TNAP expression by primary murine osteoblasts. We found both PHOSPHO1 and TNAP expression to be down-regulated in response to phosphate and PTH. While matrix mineralization was increased with phosphate (Pi), it decreased with PTH and FGF23 had no effect. The in vitro data suggest that the TNAP reduction in CKD-MBD is driven by the hyperphosphatemia and/or hyperparathyroidism noted in these mice, while the higher PHOSPHO1 expression may be a compensatory mechanism in an attempt to protect the bone from hypomineralization. We propose that increased PHOSPHO1 expression in ROD may contribute to the disordered skeletal mineralization characteristic of this progressive disorder. Design, Mutate, Screen: High-throughput creation of genetic clocks with different period-amplitude characteristics Directed evolution has become an invaluable tool in protein engineering and has also greatly influenced the construction of synthetic gene circuits. The ability to generate diversity at precise targets for directed evolution approaches has improved vastly, allowing researchers to create large, specific mutant libraries with relative ease. Screening approaches for large mutant libraries have similarly come a long way, especially when the desired behavior can easily be tested for with static, single time-point assays. For more complex gene circuits with dynamic phenotypes that change over time, directed evolution approaches to controlling and tuning circuit behavior have been hindered by the lack of sufficiently high-throughput screening methods to isolate variants with desired characteristics. Here we utilize directed mutagenesis and multiplexed microfluidics to develop a workflow for creating, screening and tuning dynamic gene circuits that operate at the population level. Specifically, we create a mutant library of an existing oscillator, the synchronized lysis circuit, and tune its dynamics while uncovering principles regarding its behavior. Lastly, we utilize this directed evolution workflow to construct a new synchronized genetic oscillator that exhibits robust dynamics over long time scales. Natural Selection as the Sum Over All Histories. As evolution can be connected to the principle of least action, and if it is depicted in evolution-space versus time then it corresponds to the direction of ultimate causation. As an organism evolves and follows a path of proximate causation, if the vector is closely parallel to that of the Ultimate Causation then the changes will confer desirable attributes which will lead to further development. If, however, the variations do not occur in a direction close to that of the ultimate causation vector the evolved organism will quickly die out. Therefore Natural Selection may be viewed as similar to Feynmans "sum over all histories". This approach is compatible with both Neutral Theory and Selection, as it includes both positive and negative mutations and selection. Therefore, the principle of least action gives a direction, but not a purpose, to evolution. Proteogenomics: an analysis of how and when itenhances proteomics BackgroundProteogenomics aims to identify variant or unknown proteins in bottom-up proteomics, by searching transcriptome- or genome-derived custom protein databases. However, empirical observations reveal that these large proteogenomic databases produce lower-sensitivity peptide identifications. Various strategies have been proposed to avoid this, including the generation of reduced transcriptome-informed protein databases (i.e., built from reference protein databases only retaining proteins whose transcripts are detected in the sample-matched transcriptome), which were found to increase peptide identification sensitivity. Here, we present a detailed evaluation of this approach.

ResultsFirst, we established that the increased sensitivity in peptide identification is in fact a statistical artifact, directly resulting from the limited capability of target-decoy competition to accurately model incorrect target matches when using excessively small databases. As anti-conservative FDRs are likely to hamper the robustness of the resulting biological conclusions, we advocate for alternative FDR control methods that are less sensitive to database size. Nevertheless, reduced transcriptome-informed databases are useful, as they reduce the ambiguity of protein identifications, yielding fewer shared peptides. Furthermore, searching the reference database and subsequently filtering proteins whose transcripts are not expressed reduces protein identification ambiguity to a similar extent, but is more transparent and reproducible.

ConclusionIn summary, using transcriptome information is an interesting strategy that has not been promoted for the right reasons. While the increase in peptide identifications from searching reduced transcriptome-informed databases is an artifact caused by the use of an FDR control method unsuitable to excessively small databases, transcriptome information can reduce ambiguity of protein identifications. Neuroinvasion and neurotropism by SARS-CoV-2 variants in the K18-hACE2 mouse Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) not only affects the respiratory tract but also causes neurological symptoms such as loss of smell and taste, headache, fatigue or severe cerebrovascular complications. Using transgenic mice expressing human angiotensin-converting enzyme 2 (hACE2) we investigated the spatiotemporal distribution and pathomorphological features in the CNS following intranasal infection with SARS-CoV-2 variants, also after prior influenza A virus infection. Apart from Omicron, we found all variants to frequently spread to and within the CNS. Infection was restricted to neurons and appeared to spread from the olfactory bulb mainly in basally orientated regions in the brain and into the spinal cord, independent of ACE2 expression and without evidence of neuronal cell death, axonal damage or demyelination. However, microglial activation, microgliosis and a mild macrophage and T cell dominated inflammatory response was consistently observed, accompanied by apoptotic death of endothelial, microglial and immune cells, without their apparent infection. Microgliosis and immune cell apoptosis indicate a potential role of microglia for pathogenesis and viral effect in COVID-19 and possible impairment of neurological functions, especially in long COVID. These data may also be informative for the selection of therapeutic candidates, and broadly support investigation of agents with adequate penetration into relevant regions of the CNS. Optimized Multicolour Immunofluorescence Panel for Cattle T Cell Phenotyping by an 8-Colour, 10-Parameter Panel Purpose and appropriate sample typesThis multiplex staining panel was developed to differentiate cattle T cells into conventional (CD4 and CD8) and unconventional ({gamma}{delta}-TCR) subsets as well as their stage of differentiation and activation. The combination of CD45RO and CD62L allows the identification of naive (TNaive), central memory (TCM), effector memory (TEM) and terminal effector (TTE) T cells. Activated cattle T cells (TAV) can be identified by the cell surface expression of CD25. This panel was developed using cryopreserved cattle peripheral blood mononuclear cells (PBMCs) and tested on fresh as well as stimulated PBMCs. Therefore, this 8-colour, 10-parameter flow cytometry panel simultaneously identifies cattle TNaive, TAV, TCM, TEM, TTE and {gamma}{delta}-TCR cells. This panel will improve our ability to examine T cell response to pathogens and vaccines in cattle including the potential to identify previously undescribed subpopulations. Furthermore, this panel can be readily optimised for other bovid species as many of these reagents are likely to cross react. 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. Analysis of science journalism reveals gender and regional disparities in coverage Science journalism is a critical way in which the public can remain informed and benefit from new scientific findings. Such journalism also shapes the publics view of the current state of scientific findings and legitimizes experts. Those covering science can only cite and quote a limited number of sources. Sources may be identified by the journalists research or by recommendations by other scientists. In both cases, biases may influence who is identified and ultimately included as an expert. To analyze possible biases in science journalism, we analyzed 22,001 non-research articles published by Nature. We chose to analyze Nature non-research articles since its research articles provide a natural comparator. Our analysis considered two possible sources of disparity: gender and name origin. To explore these sources of disparity, we extracted cited authors names as well as extracted names of quoted speakers. While citations and quotations within a piece do not reflect the entire information-gathering process, they can provide insight into the demographics of visible sources. We then used the extracted names to predict gender and name origin of the cited authors and speakers.

In order to appropriately quantify the level of difference, we must identify a suitable reference set for comparison. We chose first and last authors within primary research articles in Nature and a subset of Springer Nature articles in the same time period as our comparator. In our analysis, we found a skew towards male quotation in Nature science journalism-related articles. However, quotation is trending toward equal representation at a faster rate than first and last authorship in academic publishing. Interestingly, we found that the gender disparity in Nature quotes was column-dependent, with the "Career Features" column reaching gender parity. Our name origin analysis found a significant over-representation of names with predicted Celtic/English origin and under-representation of names with a predicted East Asian origin. This finding was observed both in extracted quotes and journal citations, but dampened in citations. Numerical Simulation of Concussive-generated Cortical Spreading Depolarization to Optimize DC-EEG Electrode Spacing for Non-invasive Visual Detection BackgroundCortical Spreading Depolarization (SD) is a propagating depolarization wave of neurons and glial cells in the cerebral gray matter. SD occurs in all forms of severe acute brain injury as documented using invasive detection methods. Based on many experimental studies of mechanical brain deformation and concussion, the occurrence of SDs in human concussion has often been hypothesized. However, this hypothesis cannot be confirmed in humans as SDs can only be detected with invasive detection methods that would require either a craniotomy or a burr hole to be performed on athletes. Typical electroencephalography (EEG) electrodes, placed on the scalp, can detect the possible presence of SD but have not been able to accurately and reliably identify SDs.

MethodsTo explore the possibility of a non-invasive method to resolve this hurdle, we developed a finite element numerical model that simulates scalp voltage changes that are induced by a brain-surface SD. We then compared our simulation results with retrospectively evaluated data in aneurysmal subarachnoid hemorrhage (aSAH) patients from Drenckhahn et al. (Brain 135:853, 2012).

ResultsThe ratio of peak scalp to simulated peak cortical voltage, Vscalp/Vcortex, was 0.0735, whereas the ratio from the retrospectively evaluated data was 0.0316 (0.0221, 0.0527) [median (1st quartile, 3rd quartile), n = 161, p < 0.001, one sample Wilcoxon signed rank test]. These differing values provide validation because their differences can be attributed to differences in shape between concussive- and aSAH-SDs, as well as the inherent limitations in human study voltage measurements. This simulated scalp surface potential was used to design a virtual scalp detection array. Error analysis and visual reconstruction showed that 1 cm is the optimal electrode spacing to visually identify the propagating scalp voltage from a cortical SD. Electrode spacings of 2 cm and above produce distorted images and high errors in the reconstructed image.

ConclusionOur analysis suggests that concussive (and other) SDs can be detected from the scalp, which could confirm SD occurrence in human concussion, provide concussion diagnosis based on an underlying physiological mechanism, and lead to non-invasive SD detection in the setting of severe acute brain injury. Frequency modulated timer regulates mammalian hibernation Mammalian hibernators decrease basal metabolism and body temperature (Tb) to minimize energy expenditure in harsh seasons. During hibernation, Tb drops to low temperature (<10 {degrees}C) and remains constant for days, known as deep torpor in small mammalian hibernators. Spontaneous interbout arousals interrupt torpor bouts, when Tb recovers to euthermic state ~37 {degrees}C. Torpor-interbout arousal event repeats during several months of hibernation. However, little is known about mechanisms governing Tb fluctuation across torpor-interbout arousal cycles during hibernation. Recent improvement in data-logging techniques enables us to monitor Tb for more than hundred days with high precision, opening up new avenues for quantitative analysis to address the principle governing Tb fluctuation. Here, we analyzed Tb fluctuation across torpor-interbout arousal cycle of Syrian hamster, which can hibernate in response to chronic cold and short photoperiod under a laboratory condition, using generalized harmonic analysis and discovered a model with frequency modulation quantitatively reproducing Tb fluctuation. This analysis also identified an unexpectedly longer period of 120-430 days as the period that modulates another period of several days, generating Tb fluctuation for Syrian hamster. We propose that concerted action of two endogenous periods governs torpor-interbout arousal cycles during hibernation. 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.