Correlation between speed and turning naturally arises for sparsely sampled cell movements Mechanisms regulating cell movement are not fully understood. One feature of cell movement that determines how far cells displace from an initial position is persistence, the ability to perform movements in a direction similar to the previous movement direction. Persistence is thus determined by turning angles between two sequential displacements and be characterized by an average turning angle or persistence time. Recent studies found that a cells average speed and turning are negatively correlated, suggesting a fundamental cell-intrinsic program whereby cells with a lower turning ability (i.e., larger persistence time) are intrinsically faster (or faster cells turn less). By simulating correlated or persistent random walks (PRWs) using two different frameworks (one based on von Mises-Fisher (vMF) distribution and another based on Ornstein-Uhlenbeck (OU) process) we show that the negative correlation between speed and turning naturally arises when cell trajectories are sub-sampled, i.e., when the frequency of sampling is lower than frequency at which cells make movements. This effect is strongest when the sampling frequency is on the order of magnitude with the typical cell persistence time and when cells vary in persistence time. Both conditions are observed for datasets of T cell movements in vivo that we have analyzed. In simulations the correlation arises due to randomness of cell movements resulting in highly variable persistence times for individual cells that with sub-sampling leads to large variability of average cell speeds. Interestingly, previously suggested methodology of calculating displacement of cell cohorts with different speeds resulted in similar results whether or not there is a cell-intrinsic correlation between cell speed and persistence. For both vMF- and OU-based simulations of PRWs we could find parameter values (distribution of persistence times, speeds, and sampling frequency) that matched experimentally measured correlations between speed and turning for two datasets of T cell movement in vivo suggesting that such simple correlations are not fully informative on the intrinsic link between speed and persistence. Our results thus suggest that sub-sampling may contribute to (and perhaps fully explains) the observed correlation between speed and turning at least for some cell trajectory data and emphasize the role of sampling frequency in inference of critical cellular parameters of cell motility such as speeds.

Secondary AbstractMeasurement of cell movements often results in a negative correlation between average speed and average turning angle suggesting an existence of a universal, cell-intrinsic movement program. We show that such a negative correlation may arise if cells in the population differ in their ability for persistent movement when the movement data are sub-sampled. We show that sub-sampling of cell trajectories generated using two different frameworks of persistent random walk can match the experimentally observed correlation between speed and turning for T cell movements in vivo.
Lymph node swelling combined with temporary effector T cell retention aids T cell response in a model of adaptive immunity Swelling of lymph nodes is commonly observed during the adaptive immune response, yet the impact on T cell trafficking and subsequent immune response is not well known. To better understand the effect of macro-scale alterations, we developed an agent-based model of the lymph node para-cortex, describing T cell trafficking and response to antigen-presenting dendritic cells alongside swelling-induced changes in T cell recruitment and egress, and regulation of expression of egress-modulating T cell receptor Sphingosine-1-phosphate receptor-1. Analysis of effector T cell response under varying swelling conditions showed that swelling consistently aided T cell activation. However, subsequent effector CD8+ T cell production could be reduced if swelling occurred too early in the T cell proliferative phase or if T cell cognate frequency was low due to increased opportunity for T cell exit. Temporarily extending retention of newly differentiated effector T cells, mediated by Sphingosine-1-phosphate receptor-1 expression, mitigated any negative effects of swelling by allowing facilitation of activation to outweigh increased access to exit areas. These results suggest targeting temporary effector T cell retention and egress associated with swelling may offer new ways to modulate effector TC responses in, for example, immuno-suppressed patients and optimisation of vaccine design.
An Extended Motif in the SARS-CoV-2 Spike Modulates Binding and Release of Host Coatomer in Retrograde Trafficking {beta}-Coronaviruses such as SARS-CoV-2 hijack coatomer protein-I (COPI) for spike protein retrograde trafficking to the progeny assembly site in endoplasmic reticulum-Golgi intermediate compartment (ERGIC). However, limited residue-level details are available into how the spike interacts with COPI. Here we identify an extended COPI binding motif in the spike that encompasses the canonical K-x-H dibasic sequence. This motif demonstrates selectivity for COPI subunit. Guided by an in silico analysis of dibasic motifs in the human proteome, we employ mutagenesis and binding assays to show that the spike motif terminal residues are critical modulators of complex dissociation, which is essential for spike release in ERGIC. COPI residues critical for spike motif binding are elucidated by mutagenesis and crystallography and found to be conserved in the zoonotic reservoirs, bats, pangolins, camels, and in humans. Collectively, our investigation on the spike motif identifies key COPI binding determinants with implications for retrograde trafficking. 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. 3D Reconstructions of Mouse Skeletal Muscle Reveal a Decrease in the MICOS Complex and Altered Mitochondrial Networks Background: Skeletal muscle gradually loses mass, strength, endurance, and oxidative capacity during aging. Studies of bioenergetics and protein turnover show that mitochondria mediate this decline in function. Mitochondria cristae are essential for the production of ATP. While mitochondrial aging is associated with endoplasmic reticulum stress, fragmented mitochondria, and decreased mitochondrial capacity, the genes associated with morphological changes in mitochondria during aging still requires further elucidation Further, it is not completely understood how 3D mitochondrial networks and the specialization of mitochondria are altered during aging. Methods: We measured changes in mitochondrial morphology and mitochondrial connectivity during the aging of the mouse gastrocnemius muscle through serial block facing-scanning electron microscopy and 3D reconstruction. Nanotunnels were also measured through 3D reconstruction. CRISPR/Cas9 KO was performed to examine changes in mitochondria upon loss of the mitochondrial contact site and cristae organizing system (MICOS) complex and optic atrophy 1 (OPA-1). Metabolomics were used to find key metabolite pathways changed upon off of the MICOS complex. Results: We found changes in mitochondrial network configuration, nanotunneling, size, shape, number, contact sites, cristae organizing system (MICOS) dynamics and gene expression in skeletal muscle during aging. We also found an association of OPA-1 and the MICOS complex in the gastrocnemius with mitochondrial aging. Further, the loss of the MICOS complex was linked with decreased oxidative capacity, and altered mitochondrial metabolism. Conclusions: To our knowledge, studies of mitochondrial changes in skeletal muscle throughout aging remains limited. 3D reconstructions of nanotunnels elucidated novel patterns of mitochondrial aging in skeletal muscle. MICOS proteins decreased with age and mitochondrial morphology was similar between aged skeletal muscle and that of young mice with MICOS protein loss. In tandem, our data suggest a relationship between the MICOS complex and aging, which could be potentially linked to disease states with additional 3D reconstruction. Commercial porcine gastric mucin contributes to variation in production of small molecule virulence factors by Pseudomonas aeruginosa when cultured in different formulations of artificial sputum medium Site-specific in vitro culture media are being developed to investigate microbial pathogenicity and ecology in nutrient environments that are more reflective of disease. For microbial ecology research in cystic fibrosis (CF), different artificial sputum media (ASM) formulations have been created to recapitulate the nutrient availability of the CF lung environment. However, these ASM formulations vary in concentration of amino acids, mucin, and other niche-specific compounds. Here, we measured the differential production of small molecule virulence factors by Pseudomonas aeruginosa, the predominant pathogen in CF pulmonary infections, cultured in nine different ASM formulations via liquid chromatography tandem mass spectrometry (LC-MS/MS) and molecular networking. We show that different ASM formulations lead to different phenotypes and metabolic profiles of P. aeruginosa and commercial porcine gastric mucin (PGM) contains a myriad of contaminants, including iron, which affect P. aeruginosa physiology.

IMPORTANCEDifferent media formulations aiming to replicate in vivo infection environments contain different nutrients, which affects interpretation of experimental results. Inclusion of undefined components, such as commercial porcine gastric mucin (PGM), in an otherwise chemically defined medium can alter the nutrient content of the medium in unexpected ways and influence experimental outcomes. A platform of patient-derived microtumors identifies therapeutic vulnerabilities in ovarian cancer. BackgroundIn light of the frequent development of therapeutic resistance in cancer treatment, there is a strong need for personalized model systems representing patient tumor heterogeneity, while enabling parallel drug testing and identification of appropriate treatment responses in individual patients. Using ovarian cancer as a prime example of a heterogeneous tumor disease with complex microenvironment and high recurrence rates, we developed a 3D preclinical ovarian cancer model comprised of patient-derived microtumors (PDM) and autologous tumor-infiltrating lymphocytes (TILs) for identification of treatment vulnerabilities and validation of anti-cancer drug efficacy using immunohistochemistry, immune cell phenotyping, functional assays and protein profiling analyses.

MethodsPDM and TILs were isolated from fresh primary ovarian cancer tissue specimen using mechanical disruption and limited enzymatic digestion and were subsequently cultured in suspension in defined media in the absence of serum. The heterogeneous cellular composition of isolated PDM as well as autologous TILs was analyzed by FFPE immunohistochemistry and multi-color flow cytometry, respectively. For in-depth protein profiling of PDM we established Reverse Phase Protein Array (RPPA) analyses of >110 total and phospho-proteins. Treatment efficacy in response to chemotherapeutics as well as immunotherapeutic compounds was assessed in PDM and PDM-TIL co-cultures using a functional viability assay in microplate format.

ResultsThe enzymatic digestion of primary ovarian cancer tissue and suspension culture in defined serum-free media allowed fast and efficient recovery of patient-derived microtumors (PDM). Immunohistochemical analyses demonstrated histopathological comparability of ovarian cancer PDM with corresponding patient tumor tissue. Reverse Phase Protein Array (RPPA)-based analyses of >110 total and phospho-proteins enabled the identification of patient-specific sensitivities to standard, platinum-based therapy and thereby the prediction of potential treatment-responders. Finally, combining PDM and autologous TILs for individual efficacy testing of immune checkpoint inhibitor treatment demonstrated the potential for patient-specific enhancement of cytotoxic TIL activity by this therapeutic approach.

Conclusion3D patient-derived ovarian cancer microtumors represent a preclinical, ex vivo tumor model that reflects intertumoral heterogeneity and represent the cellular complexity of individual patient tumors. Combining protein pathway analysis and anti-cancer drug efficacy testing of PDM enables drug mode-of-action analyses and therapeutic sensitivity prediction within a clinically relevant time frame after surgery. Follow-up studies in larger cohorts are currently under way to further evaluate the applicability of this platform to support clinical decision-making and personalizing cancer treatment. FACT regulates pluripotency through distal regulation of gene expression in murine embryonic stem cells The FACT complex is a conserved histone chaperone with essential roles in transcription and histone deposition. FACT is essential in pluripotent and cancer cells, but otherwise dispensable for most mammalian cell types. FACT deletion or inhibition can block induction of pluripotent stem cells, yet the mechanism through which FACT regulates cell fate decisions remains unclear. To determine this mechanism, we used inducible depletion of FACT subunit SPT16 in murine embryonic stem cells paired with genomic factor localization, nascent transcription, and chromatin accessibility analyses. Over a timecourse of SPT16 depletion, nucleosomes invade loci bound by master pluripotency factors and gene-distal DNaseI hypersensitive sites. Simultaneously, transcription of Pou5f1 (OCT4), Sox2, Nanog, and enhancer RNAs produced at the genes' associated enhancers are downregulated, suggesting that FACT regulates expression of the pluripotency factors themselves. We find that FACT maintains cellular pluripotency through a precise nucleosome-based regulatory mechanism for appropriate expression of both coding and non-coding transcripts associated with pluripotency. Natural variations in reproductive timing of brood size and X-chromosome nondisjunction in Caenorhabditis elegans Reproductive timing underlies self-sperm production and X-chromosome nondisjunction rate in Caenorhabditis elegans. These traits can be measured as brood size and male-production rate, and brood size decreases and male-production rate increases as worms age. These phenotypic changes occur simultaneously; however, whether these changes share common underlying genetic architectures still remains unclear. To enhance our understanding of reproductive timing in C. elegans, we measured and compared these two reproductive timing-related traits of 96 wild strains during early, late and total reproductive periods. We found that the two traits exhibited natural phenotypic variation with few outliers and a similar reproductive timing pattern as previous reports. The brood size and the male-production rate were not genetically correlated, implying that they might have different genetic architectures and that meiosis quality can be maintained despite increased progeny-production in the late reproductive period caused by more self-sperm. We also identified loci and candidate genetic variants significantly associated with male-production rate in the late and total reproductive periods. Our results provide an insight into life history traits in wild C. elegans strains. Host phenology can drive the evolution of intermediate virulence strategies in some parasites AO_SCPLOWBSTRACTC_SCPLOWThe traditional mechanistic trade-offs resulting in a negative correlation between transmission and virulence are the foundation of nearly all current theory on the evolution of parasite virulence. Several ecological factors have been shown to modulate the optimal virulence strategies predicted from mechanistic trade-off models, but these ecological factors have not yet been shown to be sufficient to explain the intermediate virulence strategies observed in any natural system. The timing of seasonal activity, or phenology, is a common factor that influences the types and impact of many ecological interactions but is difficult to incorporate into virulence evolution studies. We develop a mathematical model of a disease system with seasonal host activity to study the evolutionary consequences of host phenology on the virulence of obligate-killer parasite. Results from this model demonstrated that seasonal host activity is sufficient to drive the evolution of intermediate parasite virulence in some types of natural disease systems, even when a traditional mechanistic trade-off between transmission and virulence is not assumed in the modeling framework. The optimal virulence strategy in these systems can be determined by both the duration of the host activity period as well as the variation in the host emergence timing. Parasites with low virulence strategies are favored in environments with long host activity periods and in environments in which hosts emerge synchronously. The results demonstrate that host phenology can be sufficient to select for intermediate optimal virulence strategies, providing an alternative mechanism to account for virulence evolution in some natural systems. Estimating the fundamental niche: accounting for the uneven availability of existing environmental conditions in the calibration area Studies that question important conceptual and methodological aspects of the field of ecological niche modeling (and species distribution modeling) have cast doubts on whether it is possible to estimate the fundamental niche of a species using presence-only data. The main limitation in niche estimation is that presence data come from the realized niche, which is only a subset of the fundamental niche. Most of the existing methods lack the ability to overcome this limitation and therefore fit niches that resemble the realized niche. To obtain a more accurate estimate of the fundamental niche, we propose using the geographic region that is accessible to a species (based on its dispersal ability) to determine a sampling distribution, in environmental space, from which we can quantify the likelihood of observing a particular environmental combination in a sample of presence points. We incorporate this sampling distribution into a multivariate normal model (i.e., a Mahalanobis distance model) by creating a weight function that modifies the probabilities of observing different environmental combinations in a sample of presences. This modification accounts for the uneven availability of environmental conditions. We show that the parameters of this weighted-normal model can be approximated with a maximum likelihood estimation approach; and then used to draw ellipsoids (confidence regions) that represent the fundamental niche of the species. We illustrate the application of our model with two worked examples. First, we use presence data for an invasive species and an accessible area that includes only its native range to evaluate whether the fitted model predicts confirmed establishments of the species outside its native range. Second, we use presence data for closely related species with known accessible areas to demonstrate how the different dispersal abilities of the species constrain a classic Mahalanobis distance model. Overall, we show that accounting for the distribution of environmental conditions that are accessible to a species indeed affects the estimation of the ellipsoids used to model its fundamental niche. Vesalius: high-resolution in silico anatomization of Spatial Transcriptomic data using Image Analysis Characterization of tissue architecture promises to deliver insights into development, cell communication and disease. In silico spatial domain retrieval methods have been developed for spatial transcriptomics (ST) data assuming transcriptional similarity of neighboring barcodes. However, domain retrieval approaches with this assumption cannot work in complex tissues composed of multiple cell types. This task becomes especially challenging in cellular resolution ST methods. We developed Vesalius to decipher tissue anatomy from ST data by applying image processing technology. Vesalius uniquely detected territories composed of multiple cell types and successfully recovered tissue structures in high-resolution ST data including in mouse brain, embryo, liver, and colon. Utilizing this tissue architecture, Vesalius identified tissue morphology specific gene expression and regional specific gene expression changes for astrocytes, interneuron, oligodendrocytes, and entorhinal cells in the mouse brain. Stepwise progression of beta-selection during T cell development as revealed by histone deacetylation inhibition During T cell development, the first step in creating a unique T Cell Receptor (TCR) is the genetic recombination of the TCR{beta} chain. The quality of the newly recombined TCR{beta} is assessed at the {beta}-selection checkpoint. Most cells fail this checkpoint and die, but the coordination of the complex events that control fate at the {beta}-selection checkpoint is not yet understood. We shed new light on fate determination during {beta}-selection using a selective inhibitor of histone deacetylase 6, ACY1215. ACY1215 disrupted the {beta}-selection checkpoint. Characterising the basis for this disruption revealed a new, pivotal stage in {beta}-selection, bookended by upregulation of the TCR co-receptors, CD28 and CD2 respectively. Within this  DN3bPre stage, CD5 and Lef1 are upregulated to reflect pre-TCR signalling and their expression correlates with proliferation. During this phase, ACY1215-mediated disruption of the organisation of the {beta}-selection immunological synapse is associated with a breakdown in connectivity of expression of pre-TCR, CD5 and Lef1. Subsequent deregulation of pre-TCR-induced proliferation leads to bypass of the {beta}-selection checkpoint and subsequent failure to progress. We propose that the progressive expression of CD28, CD5 and Lef1, then CD2 reports and modulates the pre-TCR signal to orchestrate passage through the {beta}-selection checkpoint. These findings suggest a refined model of {beta}-selection in which a coordinated increase in expression of pre-TCR, CD28, CD5 and Lef1 provides for escalating TCR signalling strength, and culminates in the expression of CD2 to enable exit from the {beta}-selection checkpoint. Analysis of Drosophila cardiac hypertrophy by micro-computerized tomography for genetic dissection of heart growth mechanisms Heart failure is often preceded by pathological cardiac hypertrophy, a thickening of the heart musculature driven by complex gene regulatory and signaling processes. The Drosophila heart has great potential as a genetic model for deciphering the underlying mechanisms of cardiac hypertrophy. However, current methods for evaluating hypertrophy of the Drosophila heart are laborious and difficult to carry out reproducibly. Here we demonstrate that micro-computerized tomography (microCT) is an accessible, highly reproducible method for non-destructive, quantitative analysis of Drosophila heart morphology and size. To validate our microCT approach for analyzing Drosophila cardiac hypertrophy, we show that expression of constitutively active Ras (Ras85DV12), previously shown to cause hypertrophy of the fly heart, results in significant thickening of both adult and larval heart walls when measured from microCT images. We then show using microCT analysis that genetic upregulation of store-operated Ca2+ entry (SOCE) driven by expression of constitutively active Stim (StimCA) or Orai (OraiCA) proteins also results in significant hypertrophy of the Drosophila heart, through a process that specifically depends on Orai Ca2+ influx channels. Intravital imaging of heart contractility revealed significantly reduced end diastolic dimensions in StimCA and OraiCA expressing hearts, consistent with the hypertrophic phenotype. These results demonstrate that increased SOCE activity is an important driver of hypertrophic cardiomyocyte growth, and demonstrate how microCT analysis combined with tractable genetic tools in Drosophila can be used to delineate molecular signaling processes that underlie cardiac hypertrophy and heart failure.

NEW AND NOTEWORTHYGenetic analysis of cardiac hypertrophy in Drosophila holds immense potential for the discovery of new therapeutic targets to prevent and treat heart failure. However, this potential has been hindered by a lack of rapid and effective methods for analysis of heart size in flies. Here we demonstrate that analysis of the Drosophila heart with micro-computerized tomography yields accurate and highly reproducible heart size measurements that can be used to efficiently analyze heart growth and cardiac hypertrophy in Drosophila. Secretory cells are the primary source of pIgR in small airways BackgroundLoss of secretory immunoglobulin A (SIgA) is common in COPD small airways and likely contributes to disease progression. We hypothesized loss of SIgA results from reduced expression of pIgR, a chaperone protein needed for SIgA transcytosis, in the COPD small airway epithelium.

MethodspIgR-expressing cells were defined and quantified at single-cell resolution in human airways using RNA in-situ hybridization, immunostaining, and single-cell RNA sequencing. Complementary studies in mice utilized immunostaining, primary murine tracheal epithelial cell (MTEC) culture, and transgenic mice with secretory or ciliated cell-specific knockout of pIgR. SIgA degradation by human neutrophil elastase or secreted bacterial proteases from non-typeable Haemophilus influenzae (NTHi) was evaluated in vitro.

ResultsWe found that secretory cells are the predominant cell type responsible for pIgR expression in human and murine airways. Loss of SIgA in small airways was not associated with a reduction in secretory cells but rather a reduction in pIgR protein expression despite intact PIGR mRNA expression. Neutrophil elastase and NTHi-secreted proteases are both capable of degrading SIgA in vitro and may also contribute to a deficient SIgA immunobarrier in COPD.

InterpretationLoss of the SIgA immunobarrier in small airways of patients severe COPD is complex and likely results from both pIgR-dependent defects in IgA transcytosis and SIgA degradation.

Key MessagesWhat is the key question? Localized SIgA deficiency in small airways is an established driver of COPD pathogenesis, but the mechanism of loss remains unclear. We hypothesized loss of SIgA is due to reduced numbers of pIgR-expressing cells in SIgA-deficient small airways.

What is the bottom line? pIgR is primarily expressed by secretory cells in human and murine airways. Although numbers of secretory cells are similar between SIgA-deficient and SIgA-replete airways in COPD, there is reduced expression of pIgR protein, but not mRNA, in SIgA-deficient airways. Additionally, host and bacterial proteases degrade SIgA in vitro, suggesting loss of SIgA may relate to both impaired transcytosis and increased degradation.

Why read on? This study highlights the complexity of SIgA immunobarrier maintenance and suggests that strategies aimed at restoring the SIgA immunobarrier will need to account for both impaired transcytosis and degradation by host and/or bacterial proteases. Feline tooth resorption findings from a questionnaire survey on 8115 Finnish cats and a case-control study for a more defined population drawn from the survey Tooth resorption (TR) is one of the most common dental diseases in cats. It is a painful disease characterised by progressive dental destruction, which eventually results in loss of teeth. The aetiology of the disease remains unclear, but associations with old age, breed, other oral and dental diseases, and certain environmental factors have been suspected. We collected health and environmental information of 8115 Finnish cats from 41 breeds through an online survey for cat owners to study TR in Finland with an observational study. The frequency of TR that the owner reported as veterinarian-diagnosed was 3.9% in the entire data (316/8115) and 15% in cats that the owner reported as having veterinarian-diagnosed oral or dental disease (316/2070). We determined the risk factors and breed variation of feline TR by logistic regression with the entire data and a case-control study in a more defined population drawn from the entire data. In the entire data, we defined those cats that the owners had reported as being diagnosed with TR by a veterinarian as TR-cats and compared them to all other cats. In the case-control study, we included only cats with owner-reported veterinarian-diagnosed oral and dental diseases. We defined cats with TR-diagnosis made under sedation as TR-cases. Other veterinarian-diagnosed TR-cats were excluded. Cats with other oral and dental diseases were used as controls. The model with the entire population included 308 TR-cats compared with 7508 cats, and the case-control study included 198 TR-cats and 1683 controls. Although the diagnoses were reportedly made by a veterinarian, the owner answered the survey and might have understood, remembered, or interpreted the diagnoses in a way that could affect the results. The results indicated an increasing risk of TR with age. Dental calculus, gingivitis, and periodontitis were potential risk factors for TR. These findings and the interaction between dental calculus and gingivitis suggest that inflammatory changes caused by dental calculus increase the risk of TR. Keeping food available constantly was a potential protective factor. We found that Cornish Rex, European, and Ragdoll breeds were at higher risk for TR. The observed differences between breeds highlight a genetic contribution. 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. EEG biomarkers of reduced inhibition in human cortical microcircuits in depression Reduced cortical inhibition by somatostatin-expressing (SST) interneurons has been strongly associated with treatment-resistant depression. However, whether the effects of reduced SST interneuron inhibition on microcircuit activity have signatures detectible in electroencephalography (EEG) signals remains unknown. We simulated resting-state activity and EEG using detailed models of human cortical microcircuits with normal (healthy) or reduced SST interneuron inhibition (depression). Healthy microcircuit models showed emergent key features of resting-state EEG, and depression microcircuits exhibited increased theta, alpha and low beta power (4 - 15 Hz). The changes in depression involved a combination of an aperiodic broadband, and periodic theta and low beta components. We then demonstrated the specificity of the EEG signatures of reduced SST interneuron inhibition by showing they were distinct from those corresponding to reduced parvalbumin-expressing (PV) interneuron inhibition. Our study thus links SST interneuron inhibition level to distinct features in EEG simulated from detailed human microcircuits, which can serve to better identify mechanistic subtypes of depression using EEG, and non-invasively monitor modulation of cortical inhibition. Multiple molecular events underlie stochastic switching between two heritable cell states in a eukaryotic system Eukaryotic transcriptional networks are often large and contain several levels of feedback regulation. Many of these networks have the ability to generate and maintain several distinct transcriptional states across multiple cell divisions and to switch between them. In certain instances, switching between cell states is stochastic, occurring in a small subset of cells of an isogenic population in a seemingly homogenous environment. Given the scarcity and unpredictability of switching in these cases, investigating the determining molecular events is challenging. White-opaque switching in the fungal species Candida albicans is an example of stably inherited cell states that are determined by a complex transcriptional network and can serve as an experimentally accessible model system to study characteristics important for stochastic cell fate switching in eukaryotes. In standard lab media, genetically identical cells maintain their cellular identity (either "white" or "opaque") through thousands of cell divisions and switching between the states is rare. By isolating populations of white or opaque cells, previous studies have elucidated the many differences between the two stable cell states and identified a set of transcriptional regulators needed for cell type switching. Yet little is known about the molecular events that determine the rare, stochastic switching events that occur in single cells. We use microfluidics combined with fluorescent reporters to directly observe rare switching events between the white and opaque states. We investigate the stochastic nature of switching by beginning with white cells and monitoring the activation of Wor1, a master regulator and marker for the opaque state, in single cells and throughout cell pedigrees. Our results indicate that switching requires two stochastic steps; first an event occurs that predisposes a lineage of cells to switch. In the second step, some but not all, of those predisposed cells rapidly express high levels of Wor1 and commit to the opaque state. To further understand the rapid rise in Wor1, we used a synthetic inducible system in Saccharomyces cerevisiae into which a controllable C. albicans Wor1 and a reporter for its transcriptional control region have been introduced. We document that Wor1 positive autoregulation is highly cooperative (Hill coefficient > 3), leading to rapid activation and producing an "all or none" rather than a graded response. Taken together, our results suggest that reaching a threshold level of a master regulator is sufficient to drive cell type switching in single cells and that an earlier molecular event increases the probability of reaching that threshold in certain small lineages of cells. Quantitative molecular analysis of the white-opaque circuit can serve as a model for the general understanding of complex circuits. First detection of herpesvirus and prevalence of mycoplasma infection in free-ranging Hermann s tortoises (Testudo hermanni), and in potential pet vectors AO_SCPLOWBSTRACTC_SCPLOWTwo types of pathogens cause highly contagious upper respiratory tract diseases (URTD) in Chelonians: testudinid herpesviruses (TeHV) and a mycoplasma (Mycoplasma agassizii). In captivity, these infections are frequent and can provoke outbreaks. Pet trade generates international flow of tortoises, often without sanitary checking; individuals intentionally or accidentally released in the wild may spread pathogens. A better understanding of the transmission of infectious agents from captivity to wild tortoises is needed. Many exotic individuals have been introduced in populations of the endangered western Hermanns tortoise (Testudo hermanni hermanni), notably spur-thighed tortoises (Testudo graeca). We assessed the presence of TeHV and mycoplasma in native western Hermanns tortoises and in potential pet vectors in south-eastern France. Using a large sample (N=572 tortoises), this study revealed, by PCR, the worrying presence of herpesvirus in 7 free-ranging individuals (3 sub-populations). Additionally, Mycoplasma agassizii was detected, by PCR, in 15 of the 18 populations sampled with a frequency ranging from 3.4% (1 of 29 tortoises) to 25% (3 of 12 tortoises). Exotic spur-thighed tortoises showed high frequency of mycoplasma infection in captivity (18.2%) and in individuals (50%) found in native Hermanns tortoise sub-populations, suggesting that this species could be a significant vector. The paucity of information of TeHV on European tortoise URTD in natural settings, especially in combination with mycoplasma, prompts for further studies. Indeed, sick tortoises remain concealed and may not be easily detected in the field. Our results indicate a good health for most infected tortoise but it should be screened in the field as well as in captivity. Extending Bayesian Elo-rating to quantify the steepness of dominance hierarchies The steepness of dominance hierarchies provides information about the degree of competition within animal social groups and is thus an important concept in socioecology. The currently most widely-used metrics to quantify steepness are based on Davids scores (DS) derived from dominance interaction networks. One serious drawback of these DS-based metrics is that they are biased, i.e., network density systematically decreases steepness values. Here, we provide a novel approach to estimate steepness based on Elo-ratings, implemented in a Bayesian framework (STEER: Steepness estimation with Elo-rating). Our new metric has two key advantages. First, STEER is unbiased, precise and more robust to data density than DS-based steepness. Second, it provides explicit probability distributions of the estimated steepness coefficient, which allows uncertainty assessment. In addition, it relies on the same underlying concept and is on the same scale as the original measure, and thus allows comparison to existing published results. We evaluate and validate performance of STEER by means of experimentation on empirical and artificial data sets and compare its performance to that of several other steepness estimators. Our results suggest that STEER provides a considerable improvement over existing methods to estimate dominance hierarchy steepness. We provide an R package EloSteepness to calculate the new steepness measure, and also show an example of using steepness in a comparative analysis. Intramolecular 13C signals inform about carbon allocation in plants Stable isotopes at natural abundance are key tools to study physiological processes occurring outside the temporal scope of manipulation and monitoring experiments. Whole-molecule carbon isotope ratios (13C/12C) enable assessments of plant carbon uptake yet conceal information about carbon allocation. Here, we identify an intramolecular 13C/12C signal at treering glucose C-5 and C-6 and develop experimentally testable theories on its origin. More specifically, we assess the potential of processes within C3 metabolism for signal introduction based (inter alia) on constraints on signal propagation posed by metabolic networks. We propose that the intramolecular signal reports carbon allocation into major metabolic pathways in actively photosynthesising leaf cells including the anaplerotic, shikimate, and non-mevalonate pathway. We support our theoretical framework by linking it to previously reported whole-molecule 13C/12C increases in cellulose of ozone-treated Betula pendula and a highly significant relationship between the intramolecular signal and tropospheric ozone concentration. Our theory postulates a pronounced preference of leaf-cytosolic triose-phosphate isomerase to catalyse the forward reaction in vivo (dihydroxyacetone phosphate to glyceraldehyde 3-phosphate). In conclusion, intramolecular 13C/12C analysis resolves information about carbon uptake and allocation enabling more comprehensive assessments of carbon metabolism than whole-molecule 13C/12C analysis.

HighlightIntramolecular 13C/12C analysis resolves information about carbon uptake and allocation (and associated environmental controls) enabling more comprehensive assessments of carbon metabolism, plant-environment interactions, and environmental variability than whole-molecule 13C/12C analysis. A novel method for producing functionalized vesicles that efficiently deliver oligonucleotides in vitro in cancer cells and in vivo in mice Nano-based delivery systems have enhanced our ability to administer and target drugs and macromolecules to their targets. Oligonucleotide drugs have great therapeutic potential but often have off-target effects and stability issues. Therefore, they are often encapsulated in vesicles with targeting ligands such as antibodies (Ab) to deliver their cargo. Herein, we describe a novel, scalable and straightforward approach to producing functionalized vesicles called the "Functionalized Lipid Insertion Method," which differs from the older "Detergent-Dialysis Method." The latter method required excess detergent and extensive dialysis over many hours to produce the functionalized vesicles. With our method, only the functionalized lipid is detergent-solubilized during the engineering of the vesicle. The approach reduces the dialysis time, keeps the vesicles intact while orienting the targeting moieties of the functionalized lipid toward the outside of the vesicle. Pilot in vitro and in vivo experiments was performed to show the feasibility of our method. Dynamic light scattering (DLS) experiments suggested that the original vesicular structure was relatively unperturbed, and the functionalized lipid was inserted externally. Our approach efficiently delivered oligonucleotides and affected the function of liver cancer HepG2 cells. Furthermore, functionalized vesicles achieved targeted delivery of oligonucleotides in mice without inducing a significant innate immune response. The industrial and therapeutic significance and implications of functionalized vesicles produced by our method are also discussed. Additional experiments and analyses are recommended to bring out the full potential of this molecular delivery technology. F1F0 ATP Hydrolysis is a Determinant of Metabolic Rate, a Correlate of Lifespan, and a Weakness of Cancer Dramatic rectal temperature decrease, in mice administered with a drug that specifically inhibits F1F0 ATP hydrolysis, suggests that F1F0 ATP hydrolysis is a major determinant of metabolic rate in vivo. Across twelve investigated species, less F1F0 ATP hydrolysis correlates with greater maximal lifespan. Specific drug inhibition of F1F0 ATP hydrolysis exerts potent anti-cancer activity in vitro. Systematic temporal analysis of peripheral blood transcriptomes using TrendCatcher identifies early and persistent neutrophil activation as a hallmark of severe COVID-19 Studying temporal gene expression shifts during disease progression provides important insights into the biological mechanisms that distinguish adaptive and maladaptive responses. Existing tools for the analysis of time course transcriptomic data are not designed to optimally identify distinct temporal patterns when analyzing dynamic differentially expressed genes (DDEGs). Moreover, there is a lack of methods to assess and visualize the temporal progression of biological pathways mapped from time course transcriptomic datasets. In this study, we developed an open-source R package TrendCatcher (https://github.com/jaleesr/TrendCatcher), which applies the smoothing spline ANOVA model and break point searching strategy to identify and visualize distinct dynamic transcriptional gene signatures and biological processes from longitudinal datasets. We used TrendCatcher to perform a systematic temporal analysis of COVID-19 peripheral blood transcriptomes, including bulk RNA-seq and scRNA-seq time course data. TrendCatcher uncovered the early and persistent activation of neutrophils and coagulation pathways as well as impaired type I interferon (IFN-I) signaling in circulating cells as a hallmark of patients who progressed to severe COVID-19, whereas no such patterns were identified in individuals receiving SARS- CoV-2 vaccinations or patients with mild COVID-19. These results underscore the importance of systematic temporal analysis to identify early biomarkers and possible pathogenic therapeutic targets. 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.